Big Data Cosmology meets AI
An Invitation to Backpropagate Through the Origins of the Universe
IAIFI Fellow
Carol Cuesta-Lazaro
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/11147457/pasted-from-clipboard.png)
Video Credit: N-body simulation Francisco Villaescusa-Navarro
![](https://galaxies-cosmology-2015.wdfiles.com/local--files/baryon-acoustic-oscillations/baryon_acoustic_peak.png)
![](https://d3i71xaburhd42.cloudfront.net/7be5f2f28a985c181d7bdcd3a01bb8dce1d86297/3-Figure1-1.png)
![](https://www.darkenergysurvey.org/wp-content/uploads/2021/07/y1bao_powerspectrum_figure2.png)
![](https://astro.ucla.edu/~wright/DL-vs-z-26Mar2015.gif)
The era of Big Data Astrophysics
![](https://upload.wikimedia.org/wikipedia/commons/3/3c/Ilc_9yr_moll4096.png)
![](https://hubblesite.org/files/live/sites/hubble/files/home/resource-gallery/articles/_images/hs-article-0720a-2400x1840.jpg?t=tn2400)
![](https://news.fnal.gov/wp-content/uploads/2024/01/1-Deep-image-SN-only.jpg)
1-Dimensional
![](https://www.ligo.caltech.edu/system/pages/images/24/page/Gravity_Waves_StillImage.jpg?1699659823)
![](https://upload.wikimedia.org/wikipedia/commons/3/3c/Ilc_9yr_moll4096.png)
![](https://news.fnal.gov/wp-content/uploads/2024/01/1-Deep-image-SN-only.jpg)
![](https://www.ligo.caltech.edu/system/pages/images/24/page/Gravity_Waves_StillImage.jpg?1699659823)
![](https://hubblesite.org/files/live/sites/hubble/files/home/resource-gallery/articles/_images/hs-article-0720a-2400x1840.jpg?t=tn2400)
Machine Learning
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/11144498/pasted-from-clipboard.png)
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/11144498/pasted-from-clipboard.png)
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/11147279/pasted-from-clipboard.png)
Cosmic Cartography
Galaxy Clustering
Galaxy Imaging
Lensing
Cosmic Microwave Background
Gravitational Waves
Time domain
Early Universe Inflation
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/11144498/pasted-from-clipboard.png)
Late Universe
![](https://upload.wikimedia.org/wikipedia/commons/3/3c/Ilc_9yr_moll4096.png)
What's the Universe made of?
Evolution
Dark matter
Dark energy
Non-Gaussianity
...
Multifield Inflation
Initial Conditions
The Universe's forward model
Observables
Why Astrophysics is hard 101
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/11144498/pasted-from-clipboard.png)
Dataset Size = 1
Can't poke it in the lab
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/11335545/pasted-from-clipboard.png)
Simulations
Bayesian statistics
How well can we simulate the Universe?
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/11335519/pasted-from-clipboard.png)
Very interested on ideas in the area of model mispecification!
How do we learn what is the robust information?
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/11335519/pasted-from-clipboard.png)
Simulating dark matter is easy!
"Atoms" are hard" :(
Hybrid ML - Physics Simulators
Unsupervised searches
Cosmological (field level) Inference for Galaxy Surveys
DESI
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/11144498/pasted-from-clipboard.png)
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/11298827/pasted-from-clipboard.png)
High dimensional data
Unknown
Simple summary statistic
estimated with Perturbation Theory
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/11144498/pasted-from-clipboard.png)
![](https://galaxies-cosmology-2015.wdfiles.com/local--files/baryon-acoustic-oscillations/baryon_acoustic_peak.png)
Probability pair of galaxy
Pair separation
Forward Model
Parameters
Observable
Likelihood
Simulator
+ MCMC hammer
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/11130639/pasted-from-clipboard.png)
Dark matter
Dark energy
Inflation
Perturbation Theory
Pen and paper
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/11144507/pasted-from-clipboard.png)
![](https://galaxies-cosmology-2015.wdfiles.com/local--files/baryon-acoustic-oscillations/baryon_acoustic_peak.png)
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/11144498/pasted-from-clipboard.png)
+ Density Estimation
+ Sampler
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/11040096/diffusion_fig-1.png)
"A point cloud approach to generative modeling for galaxy surveys at the field level"
Cuesta-Lazaro and Mishra-Sharma
arXiv:2311.17141
Base Distribution
Target Distribution
- Sample
- Evaluate
![](https://www.smsharma.io/assets/img/siddharth_headshot.png)
Siddharth Mishra-Sharma
Long range correlations
Huge pointclouds (20M)
Homogeneity and isotropy
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/11139914/darkcosmology__1_.gif)
Fixed Initial Conditions
Varying Cosmology
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/10597125/pasted-from-clipboard.png)
Trained on only 5000 positions!
Real observations 20 Million points :(
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/11300112/ll_prof-1.png)
Learning in 5000 dimensions with only 2000 simulations
Symmetries?
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/11335642/pasted-from-clipboard.png)
Julia Balla
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/11335645/pasted-from-clipboard.png)
Loss
Step
Pair counting
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/11335664/pasted-from-clipboard.png)
MP GNN
Hierarchical
Symmetries
"GalaxyBench: A Long- and Short-Range Benchmark for Symmetry-Preserving Data Processing" Balla et al (in prep.)
1 to Many:
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/11298868/pasted-from-clipboard.png)
Can we run the Universe backwards?
Today
Initial Conditions
"Probabilistic Forecasting with Stochastic Interpolants and Follmer Processes" Chen et al
arXiv:2403.13724
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/11338245/pasted-from-clipboard.png)
Sampling SDE
Interpolant
Drift
Regression loss
"Probabilistic Forecasting with Stochastic Interpolants and Follmer Processes" Chen et al
arXiv:2403.13724
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/11338031/pasted-from-clipboard.png)
Current model is not very good when ran forwards!
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/11338031/pasted-from-clipboard.png)
3D U-Nets are annoying :(
True
Initial
Final
Predicted
Can we run larger simulations? (Observable volumes)
At high resolution?
Faster?
All this works depends on simulations, but...
Thousands of them?
![](https://media1.tenor.com/m/VqoQZNqYdkQAAAAC/more-i-want-more.gif)
Hybrid Physical / ML simulators
Gravitational evolution ODE
Particle-mesh
"Nbodyify: Adaptive mesh corrections for PM simulations" Cuesta-Lazaro, Modi in preps
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/10812245/densities-1.png)
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/10812245/densities-1.png)
Particle-mesh
Full Nbody
Hybrid Simulator - on the fly
Gravitational evolution ODE
Trained to match particle velocities and positions: DIFFERENTIABLE
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/10826038/pasted-from-clipboard.png)
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/10826039/pasted-from-clipboard.png)
Density
Gravitational Potential
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/10826217/pasted-from-clipboard.png)
1. CNN
2. Read features at position using attention
3. Compute force correction
4. Run corrected simulation
Learn features
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/10812245/densities-1.png)
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/10812245/densities-1.png)
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/10812245/densities-1.png)
Particle-mesh
Full Nbody
Hybrid ML-Simulator
"Nbodyify: Adaptive mesh corrections for PM simulations" Cuesta-Lazaro, Modi in preps
Video credit: Francisco Villaescusa-Navarro
Gas density
Gas temperature
Finding missing physics with differentiable simulators?
What is the space of plausible solutions and how do we search it?
Differentiable Galaxies ODEs
Humans best bet
Neural Network correction
Are there problems in cosmology that bypass a forward model?
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/11298796/pasted-from-clipboard.png)
Parity violation cannot be originated by gravity
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/11298817/pasted-from-clipboard.png)
"Measurements of parity-odd modes in the large-scale 4-point function of SDSS..." Hou, Slepian, Chan arXiv:2206.03625
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/11298827/pasted-from-clipboard.png)
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/11298828/tomirror.jpg)
"Could sample variance be responsible for the parity-violating signal seen in the BOSS galaxy survey?" Philcox, Ereza arXiv:2401.09523
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/11298827/pasted-from-clipboard.png)
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/11298828/tomirror.jpg)
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/11144498/pasted-from-clipboard.png)
Train
Test
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/11300865/pasted-from-clipboard.png)
Me: I can't wait to work with observations
Me working with observations:
Finding interesting objects:
Very small galaxies (dwarf galaxies)
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/11337211/pasted-from-clipboard.png)
Interesting in Astrophysics: How we define an anomaly and how do we find it?
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/11338104/pasted-from-clipboard.png)
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/11338104/pasted-from-clipboard.png)
Background
Region of Interest
Conclusions
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/10597125/pasted-from-clipboard.png)
1. There is a lot of information in galaxy surveys that ML methods can access
2. We can tackle high dimensional inference problems so far unatainable
3. Our ability to simulate will limit the amount of information we can extract
Hybrid simulators, forward models, robustness
Unsupervised problems: parity violation
Finding anomalies for new physics?
Finding the Initial Conditions of the Universe, let's get creative!
Field level inference
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/10824461/fiducial_imshows-1.png)
![](https://s3.amazonaws.com/media-p.slid.es/uploads/993552/images/11040606/pasted-from-clipboard.png)
MIT
By carol cuesta
MIT
- 117