Higher order Hamiltonian Monte Carlo sampling for cosmological large-scale structure analysis

Hernández-Sánchez, Mónica; Kitaura, Francisco-Shu; Ata, Metin; Dalla Vecchia, Claudio
Bibliographical reference

Monthly Notices of the Royal Astronomical Society

Advertised on:
Number of authors
IAC number of authors
Refereed citations
We investigate higher order symplectic integration strategies within Bayesian cosmic density field reconstruction methods. In particular, we study the fourth-order discretization of Hamiltonian equations of motion (EoM). This is achieved by recursively applying the basic second-order leap-frog scheme (considering the single evaluation of the EoM) in a combination of even numbers of forward time integration steps with a single intermediate backward step. This largely reduces the number of evaluations and random gradient computations, as required in the usual second-order case for high-dimensional cases. We restrict this study to the lognormal-Poisson model, applied to a full volume halo catalogue in real space on a cubical mesh of 1250 h-1 Mpc side and 2563 cells. Hence, we neglect selection effects, redshift space distortions, and displacements. We note that those observational and cosmic evolution effects can be accounted for in subsequent Gibbs-sampling steps within the COSMIC BIRTH algorithm. We find that going from the usual second to fourth order in the leap-frog scheme shortens the burn-in phase by a factor of at least ∼30. This implies that 75-90 independent samples are obtained while the fastest second-order method converges. After convergence, the correlation lengths indicate an improvement factor of about 3.0 fewer gradient computations for meshes of 2563 cells. In the considered cosmological scenario, the traditional leap-frog scheme turns out to outperform higher order integration schemes only when considering lower dimensional problems, e.g. meshes with 643 cells. This gain in computational efficiency can help to go towards a full Bayesian analysis of the cosmological large-scale structure for upcoming galaxy surveys.
Related projects
Project Image
Numerical Astrophysics: Galaxy Formation and Evolution
How galaxies formed and evolved through cosmic time is one of the key questions of modern astronomy and astrophysics. Cosmological time- and length-scales are so large that the evolution of individual galaxies cannot be directly observed. Only through numerical simulations can one follow the emergence of cosmic structures within the current
Dalla Vecchia
 The Invisible Scaffolding of Space
Cosmology with Large Scale Structure Probes
The Cosmic Microwave Background (CMB) contains the statistical information about the early seeds of the structure formation in our Universe. Its natural counterpart in the local universe is the distribution of galaxies that arises as a result of gravitational growth of those primordial and small density fluctuations. The characterization of the