BullFrog: Multi-step perturbation theory as a time integrator for cosmological simulations

Abstract

Modelling the cosmic large-scale structure can be done through numerical N-body simulations or by using perturbation theory. Here, we present an N-body approach that effectively implements a multi-step forward model based on Lagrangian Perturbation Theory (LPT) in a $\Lambda$CDM Universe. This is achieved by introducing the second-order accurate BullFrog integrator, which automatically performs 2LPT time steps to second order without requiring the explicit computation of 2LPT displacements. Importantly, we show that BullFrog trajectories rapidly converge to the exact solution as the number of time steps increases, at any moment in time, even though 2LPT becomes invalid after shell-crossing. As a validation test, we compare BullFrog against other N-body integrators and high-order LPT, both for a realistic $\Lambda$CDM cosmology and for simulations with a sharp UV cutoff in the initial conditions. The latter scenario enables controlled experiments against LPT and, in practice, is particularly relevant for modelling coarse-grained fluids arising in the context of effective field theory. We demonstrate that BullFrog significantly improves upon other LPT-inspired integrators, such as FastPM and COLA, without incurring any computational overhead compared to standard N-body integrators. Implementing BullFrog in any existing N-body code is straightforward, particularly if FastPM is already integrated.

31 pages plus appendices and references, 14 figures, v2: added clarifications and results for cross-correlation, JCAP accepted

References in corpus (100)

• The cosmological simulation code GADGET-2
• Nine-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Cosmological Parameter Results
• LSST: from Science Drivers to Reference Design and Anticipated Data Products
• The 6dF Galaxy Survey: Baryon Acoustic Oscillations and the Local Hubble Constant
• Cosmological Hydrodynamics with Adaptive Mesh Refinement: a new high resolution code called RAMSES
• Large-Scale Structure of the Universe and Cosmological Perturbation Theory
• The WiggleZ Dark Energy Survey: mapping the distance-redshift relation with baryon acoustic oscillations
• Large-Scale Galaxy Bias
• Cosmic Microwave Background Anisotropies
• Transients from Initial Conditions in Cosmological Simulations
• Cosmological Non-Linearities as an Effective Fluid
• The Effective Field Theory of Cosmological Large Scale Structures
• Resumming Cosmological Perturbations via the Lagrangian Picture: One-loop Results in Real Space and in Redshift Space
• Simulating cosmic structure formation with the GADGET-4 code
• Recovery of the Power Spectrum of Mass Fluctuations from Observations of the Lyman-alpha Forest
• Fundamental Physics with the Square Kilometre Array
• Solving Large Scale Structure in Ten Easy Steps with COLA
• The Quijote simulations
• The CAMELS project: Cosmology and Astrophysics with MachinE Learning Simulations
• The Uchuu Simulations: Data Release 1 and Dark Matter Halo Concentrations
• HACC: Simulating Sky Surveys on State-of-the-Art Supercomputing Architectures
• Convolution Lagrangian perturbation theory for biased tracers
• AbacusSummit: A Massive Set of High-Accuracy, High-Resolution $N$-Body Simulations
• Learning to Predict the Cosmological Structure Formation
• FastPM: a new scheme for fast simulations of dark matter and halos
• Flowing with Time: a New Approach to Nonlinear Cosmological Perturbations
• PTHalos: A fast method for generating mock galaxy distributions
• Large-scale dark matter simulations
• Reconstruction of the early Universe as a convex optimization problem
• Lagrangian theory of gravitational instability of Friedman-Lemaitre cosmologies - generic third-order model for non-linear clustering
• Newtonian Cosmology in Lagrangian Formulation: Foundations and Perturbation Theory
• Cosmology with the Wide-Field Infrared Survey Telescope -- Multi-Probe Strategies
• Cosmological Perturbation Theory Using the FFTLog: Formalism and Connection to QFT Loop Integrals
• SWIFT: A modern highly-parallel gravity and smoothed particle hydrodynamics solver for astrophysical and cosmological applications
• ICE-COLA: Towards fast and accurate synthetic galaxy catalogues optimizing a quasi $N$-body method
• Accurate initial conditions for cosmological N-body simulations: Minimizing truncation and discreteness errors
• The Euclid mission design
• Simulation techniques for cosmological simulations
• The cosmology dependence of galaxy clustering and lensing from a hybrid $N$-body-perturbation theory model
• Coarse-Grained Cosmological Perturbation Theory
• An adaptively refined phase-space element method for cosmological simulations and collisionless dynamics
• Lagrangian perturbations and the matter bispectrum I: fourth-order model for non-linear clustering
• Lagrangian perturbation theory at one loop order: successes, failures, and improvements
• On the reach of perturbative descriptions for dark matter displacement fields
• Fast and accurate mock catalogue generation for low-mass galaxies
• Time-analyticity of Lagrangian particle trajectories in ideal fluid flow
• A rigorous EFT-based forward model for large-scale structure
• Improving Initial Conditions for Cosmological $N$-Body Simulations
• Cosmological Parameter Estimation with Large Scale Structure Observations
• On the non-linear scale of cosmological perturbation theory
• Behaviour of Lagrangian Approximations in Spherical Voids
• Linear perturbative theory of the discrete cosmological N-body problem
• Cosmological inference from Bayesian forward modelling of deep galaxy redshift surveys
• How much information can be extracted from galaxy clustering at the field level?
• Multiscale Inference of Matter Fields and Baryon Acoustic Oscillations from the Ly-alpha Forest
• An $n$-th order Lagrangian Forward Model for Large-Scale Structure
• CARPool: fast, accurate computation of large-scale structure statistics by pairing costly and cheap cosmological simulations
• GridSPT: Grid-based calculation for perturbation theory of large-scale structure
• FlowPM: Distributed TensorFlow Implementation of the FastPM Cosmological N-body Solver
• Perturbation theory approach for the power spectrum: from dark matter in real space to massive haloes in redshift space
• Cosmology Inference from Biased Tracers using the EFT-based Likelihood
• Constraining the Baryon-Dark Matter Relative Velocity with the Large-Scale 3-Point Correlation Function of the SDSS BOSS DR12 CMASS Galaxies
• A fast particle-mesh simulation of non-linear cosmological structure formation with massive neutrinos
• Contact variational integrators
• Post-collapse perturbation theory in 1D cosmology -- beyond shell-crossing
• Lagrangian perturbations and the matter bispectrum II: the resummed one-loop correction to the matter bispectrum
• Accelerating Large-Scale-Structure data analyses by emulating Boltzmann solvers and Lagrangian Perturbation Theory
• Priors on red galaxy stochasticity from hybrid effective field theory
• Consistency tests of field level inference with the EFT likelihood
• Higher-order initial conditions with massive neutrinos
• Sigma-Eight at the Percent Level: The EFT Likelihood in Real Space
• Science with the X-ray Imaging and Spectroscopy Mission (XRISM)
• Vlasov-Poisson in 1D for initially cold systems: post-collapse Lagrangian perturbation theory
• Shell-crossing in quasi-one-dimensional flow
• Simulating the Complexity of the Dark Matter Sheet I: Numerical Algorithms
• Lagrangian cosmological perturbation theory at shell-crossing
• Structure formation beyond shell-crossing: nonperturbative expansions and late-time attractors
• On Second-order Perturbation Theories of Gravitational Instability in Friedmann-Lemaitre Models
• Modeling features in the redshift-space halo power spectrum with perturbation theory
• Shell-crossing in a $\Lambda$CDM Universe
• EFTofLSS meets simulation-based inference: ${\sigma }_8$ from biased tracers
• The Cauchy-Lagrangian method for numerical analysis of Euler flow
• One line to run them all: SuperEasy massive neutrino linear response in $N$-body simulations
• Unveiling the singular dynamics in the cosmic large-scale structure
• Fast simulations of cosmic large-scale structure with massive neutrinos
• Perturbation theory with dispersion and higher cumulants: non-linear regime
• Cosmological Vlasov-Poisson equations for dark matter: Recent developments and connections to selected plasma problems
• Differentiable Cosmological Simulation with Adjoint Method
• Perturbation theory of LSS in the $\Lambda$CDM Universe: exact time evolution and the two-loop power spectrum
• Grid-based calculations of redshift-space matter fluctuations from perturbation theory: UV sensitivity and convergence at the field level
• DISCO-DJ I: a differentiable Einstein-Boltzmann solver for cosmology
• Eye of the Tyger: early-time resonances and singularities in the inviscid Burgers equation
• Analytical growth functions for cosmic structures in a $\Lambda$CDM Universe
• A Minimal Model for Massive Neutrinos in Newtonian N-body Simulations
• Renormalization group and UV completion of cosmological perturbations: Gravitational collapse as a critical phenomenon
• COmoving Computer Acceleration (COCA): $N$-body simulations in an emulated frame of reference
• Perturbation-theory informed integrators for cosmological simulations
• Starting Cosmological Simulations from the Big Bang
• The gravitational force field of proto-pancakes
• Fast and accurate collapse-time predictions for collisionless matter

Cited by in corpus (6)

• PySCo: A fast Particle-Mesh $N$-body code for modified gravity simulations in Python
• Towards an optimal extraction of cosmological parameters from galaxy cluster surveys using convolutional neural networks
• Simulation budgeting for hybrid effective field theories
• Post-collapse Lagrangian perturbation theory in three dimensions
• Benchmarking field-level cosmological inference from galaxy redshift surveys
• Cosmological Analysis with Calibrated Neural Quantile Estimation and Approximate Simulators