Astrophysics > Cosmology and Nongalactic Astrophysics
[Submitted on 6 Nov 2019 (v1), last revised 20 Jan 2021 (this version, v5)]
Title:Higher Order Hamiltonian Monte Carlo Sampling for Cosmological Large-Scale Structure Analysis
View PDFAbstract:We investigate higher order symplectic integration strategies within Bayesian cosmic density field reconstruction methods. In particular, we study the fourth-order discretisation 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 256$^3$ 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 $\sim30$. 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 256$^3$ cells. In the considered cosmological scenario, the traditional leap-frog scheme turns out to outperform higher order integration schemes only at lower dimensional problems, e.g. meshes with 64$^3$ 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.
Submission history
From: Francisco-Shu Kitaura [view email][v1] Wed, 6 Nov 2019 23:02:40 UTC (9,372 KB)
[v2] Sat, 23 May 2020 06:40:41 UTC (5,170 KB)
[v3] Sat, 17 Oct 2020 19:40:12 UTC (6,067 KB)
[v4] Tue, 15 Dec 2020 17:21:10 UTC (6,522 KB)
[v5] Wed, 20 Jan 2021 15:14:39 UTC (6,501 KB)
Current browse context:
astro-ph.CO
References & Citations
Bibliographic and Citation Tools
Bibliographic Explorer (What is the Explorer?)
Litmaps (What is Litmaps?)
scite Smart Citations (What are Smart Citations?)
Code, Data and Media Associated with this Article
CatalyzeX Code Finder for Papers (What is CatalyzeX?)
DagsHub (What is DagsHub?)
Gotit.pub (What is GotitPub?)
Papers with Code (What is Papers with Code?)
ScienceCast (What is ScienceCast?)
Demos
Recommenders and Search Tools
Influence Flower (What are Influence Flowers?)
Connected Papers (What is Connected Papers?)
CORE Recommender (What is CORE?)
IArxiv Recommender
(What is IArxiv?)
arXivLabs: experimental projects with community collaborators
arXivLabs is a framework that allows collaborators to develop and share new arXiv features directly on our website.
Both individuals and organizations that work with arXivLabs have embraced and accepted our values of openness, community, excellence, and user data privacy. arXiv is committed to these values and only works with partners that adhere to them.
Have an idea for a project that will add value for arXiv's community? Learn more about arXivLabs.