Cosmological inference at the field-level

Hugo SIMON-ONFROY, 
PhD student supervised by
Arnaud DE MATTIA and François LANUSSE

DPhP, 2025/01/16

Current estimation

A toy example

\(\Omega := \{ \Omega_c, \Omega_b, \Omega_\Lambda, H_0, \sigma_8, n_s,...\}\)

Linear matter spectrum

Structure growth

Cosmological modeling and inference

\(\Omega\)

\(\delta_L\)

\(\delta_g\)

inference

$$\boldsymbol{p}(\Omega \mid \delta_g) \propto \int \boldsymbol{p}(\Omega, \delta_L, \delta_g) \;\mathrm d \delta_L$$

\(\Omega := \{ \Omega_c, \Omega_b, \Omega_\Lambda, H_0, \sigma_8, n_s,...\}\)

Linear matter spectrum

Structure growth

Can't integrate? Then compress

\(\Omega\)

\(\delta_L\)

\(\delta_g\)

\(P\)

Power spectrum \(P(\delta_g)\) allows to analytically compute $$\boldsymbol{p}(\Omega \mid P) \propto \int \boldsymbol{p}(\Omega, \delta_L, P) \;\mathrm d \delta_L$$

Compression

 If \(\delta_g\) is Gaussian, \(P\) is a lossless compression and then $$\boldsymbol{p}(\Omega \mid P) = \boldsymbol{p}(\Omega \mid \delta_g)$$

But...

  • At large scales, matter field almost Gaussian so power spectrum is almost lossless compression
  • At smaller scales however, matter field is non-Gaussian

Gaussianity and beyond

2 fields, 1 power spectrum: Gaussian or N-body?

Field-level inference

Idea: sample simultaneously cosmology and initial field, yielding posterior on full universe history

  • High dimensional sampling: \(1024^3\) params! 
  • Use gradient-based MCMC such as Hamiltonian Monte Carlo, thanks to auto-diff

2024Mines

By hsimonfroy

2024Mines

  • 87