Baryonification
Inflation
Symmetry-preserving ML
Galaxy Formation
Simulation Based Inference
Epidemiological simulations
Medical Imaging
AI4Science
Exoplanets
LSS
Compute
Simulations
Data
ML
Statistics
Physics
What is dark matter made of?
What is driving the accelerated expansion?
How did the Universe begin?
Late Universe
Early Universe
Tension
Early vs Late
Parametric Extensions
[Image Credit: Prof. Wendy Freedman]
Systematics?
-> Shrink error bars
-> Build methods for attribution
[Image Credit: Claire Lamman (CfA/Harvard) / DESI Collaboration]
Forward Model
Observable
Predict
Infer
Theory Parameters
Inverse mapping
+ MCMC hammer
Dark matter
Dark energy
Inflation
Initial conditions
Carol's optimistic forecast
Reconstructing ALL latent variables:
Dark Matter distribution
Entire formation history
Peculiar velocities
Predictive Cross Validation:
Cross-Correlation with other probes without Cosmic Variance
[Image Credit: Yuuki Omori]
Constraining Inflation:
Inferring primordial non-gaussianity
Data-driven Subgrid models / Data-driven Systematics
1) Likelihood not necessarily Gaussian
2) Forward model no need differentiable
3) Amortized
Marginalizing over ICs
Fixing ICs
HMC: Marginalizing over ICs
True
Reconstructed
Black Hole powered jets regulate star formation
But jets interact with the turbulent interstellar medium!
Hypothesis
Simulate World
Invisible particles
Extra dimensions
Multi Species ...
Simulate World
Invisible particles
Extra dimensions
Multi Species ...
Propose Experiment
Simulate
Text: Conceptual Understanding
Trajectories (.csv)
Science Agent (LLM)
Science Agent (LLM)
Outputs
Python Code: Trajectory MSE
["DiscoverPhysics: Benchmarking LLMs for Out-of-the-Box Scientific Thinking" Wiemann, Smith et al (including CCL)]
World Generator
World Solver
def simulate(
pos1,
pos2,
duration,
**params,
):
"Simulate Universe"
return trajectories
Convergence,
Re-implementation tests....
def discovered_law(pos1, pos2, p1, p2, velocity2, duration, **params):
"""Particle 2 is accelerated toward particle 1 by a radial force
per unit along r^ that combines a static source term G*p1,
a radial-velocity term (analogous to an advective/retardation
coupling), and a centripetal-like tangential kinetic-energy
term, all scaled by 1/r and independent of p2."""
...
return trajectories
def fit_parameters():
return {
"G": {"init": 1.0, "bounds": [0.1, 5.0]},
...
}
Running Experiment...
Reward
Predictiveness
Conceptual Understanding
(MSE)
(Evaluation Score)
"This world consists of ..."
World Definition
Simulation Code
Easy
Medium
Hard
The challenge: Long roll outs / Attribution