Higher-order initial conditions for mixed baryon-CDM simulations

Abstract

We present a novel approach to generate higher-order initial conditions (ICs) for cosmological simulations that take into account the distinct evolution of baryons and dark matter. We focus on the numerical implementation and the validation of its performance, based on both collisionless N-body simulations and full hydrodynamic Eulerian and Lagrangian simulations. We improve in various ways over previous approaches that were limited to first-order Lagrangian perturbation theory (LPT). Specifically, we (1) generalize nth-order LPT to multi-fluid systems, allowing 2LPT or 3LPT ICs for two-fluid simulations, (2) employ a novel propagator perturbation theory to set up ICs for Eulerian codes that are fully consistent with 1LPT or 2LPT, (3) demonstrate that our ICs resolve previous problems of two-fluid simulations by using variations in particle masses that eliminate spurious deviations from expected perturbative results, (4) show that the improvements achieved by going to higher-order PT are comparable to those seen for single-fluid ICs, and (5) demonstrate the excellent (i.e., few per cent level) agreement between Eulerian and Lagrangian simulations, once high-quality initial conditions are used. The rigorous development of the underlying perturbation theory is presented in a companion paper. All presented algorithms are implemented in the Monofonic Music-2 package that we make publicly available.

arXiv comment20 pages, 14 figures, MNRAS in press, revised version discusses how relative velocity mode can be included to leading order, code available from https://bitbucket.org/ohahn/monofonic

References in corpus (74)

• The cosmological simulation code GADGET-2
• Efficient Computation of CMB anisotropies in closed FRW models
• The EAGLE project: Simulating the evolution and assembly of galaxies and their environments
• E pur si muove: Galiliean-invariant cosmological hydrodynamical simulations on a moving mesh
• The Cosmic Linear Anisotropy Solving System (CLASS) II: Approximation schemes
• Cosmological Hydrodynamics with Adaptive Mesh Refinement: a new high resolution code called RAMSES
• Large-Scale Structure of the Universe and Cosmological Perturbation Theory
• First results from the IllustrisTNG simulations: matter and galaxy clustering
• GIZMO: A New Class of Accurate, Mesh-Free Hydrodynamic Simulation Methods
• Dancing in the dark: galactic properties trace spin swings along the cosmic web
• Enzo: An Adaptive Mesh Refinement Code for Astrophysics
• Multi-scale initial conditions for cosmological simulations
• Cosmic Microwave Background Anisotropies
• Adaptive Refinement Tree - a new high-resolution N-body code for cosmological simulations
• Transients from Initial Conditions in Cosmological Simulations
• Relative velocity of dark matter and baryonic fluids and the formation of the first structures
• The Arepo public code release
• Fundamental Physics with the Square Kilometre Array
• Transients from Initial Conditions: A Perturbative Analysis
• Nyx: A Massively Parallel AMR Code for Computational Cosmology
• Do Dark Matter Axions Form a Condensate with Long-Range Correlation?
• Gasoline2: A Modern SPH Code
• Accurate Estimators of Correlation Functions in Fourier Space
• A new method to quantify the effects of baryons on the matter power spectrum
• Cosmological $N$-body simulations with suppressed variance
• Quantifying baryon effects on the matter power spectrum and the weak lensing shear correlation
• The Bispectrum: From Theory to Observations
• Impact of the Relative Motion between the Dark Matter and Baryons on the First Stars
• Modelling baryonic feedback for survey cosmology
• The Delay of Population III Star Formation by Supersonic Streaming Velocities
• Large-scale BAO signatures of the smallest galaxies
• Modeling baryonic physics in future weak lensing surveys
• Growth of Linear Perturbations before the Era of the First Galaxies
• Accurate initial conditions for cosmological N-body simulations: Minimizing truncation and discreteness errors
• Accurate estimators of power spectra in N-body simulations
• Modelling the large scale structure of the Universe as a function of cosmology and baryonic physics
• Initial Conditions for Accurate N-Body Simulations of Massive Neutrino Cosmologies
• The formation of the first cosmic structures and the physics of the z~20 Universe
• Recursive Solutions of Lagrangian Perturbation Theory
• Non-Linear Approximations to Gravitational Instability: A Comparison in the Quasi-Linear Regime
• The effect of relative velocity and density perturbations between baryons and dark matter on the clustering of galaxies
• Time-analyticity of Lagrangian particle trajectories in ideal fluid flow
• Resummed propagators in multi-component cosmic fluids with the eikonal approximation
• The recursion relation in Lagrangian perturbation theory
• On the signature of the baryon-dark matter relative velocity in the two and three-point galaxy correlation functions
• Supersonic Relative Velocity Effect on the Baryonic Acoustic Oscillation Measurements
• Linear perturbative theory of the discrete cosmological N-body problem
• How closely do baryons follow dark matter on large scales?
• Cosmological perturbation theory for baryons and dark matter I: one-loop corrections in the RPT framework
• Streaming velocities and the baryon-acoustic oscillation scale
• GridSPT: Grid-based calculation for perturbation theory of large-scale structure
• Baryonic effects on the matter bispectrum
• SWIFT: Using task-based parallelism, fully asynchronous communication, and graph partition-based domain decomposition for strong scaling on more than 100,000 cores
• Numerical convergence of hydrodynamical simulations of galaxy formation: the abundance and internal structure of galaxies and their cold dark matter haloes
• Quantification of discreteness effects in cosmological N-body simulations: I. Initial Conditions
• The Evolution of Baryon Density Fluctuations in Multi-Component Cosmological Simulations
• Constraining the Baryon-Dark Matter Relative Velocity with the Large-Scale 3-Point Correlation Function of the SDSS BOSS DR12 CMASS Galaxies
• How smooth are particle trajectories in a $\Lambda$CDM Universe?
• Relativistic initial conditions for N-body simulations
• Biased Tracers of Two Fluids in the Lagrangian Picture
• Gravitational evolution of a perturbed lattice and its fluid limit
• More accurate simulations with separate initial conditions for baryons and dark matter
• Cosmological $N$-body simulations including radiation perturbations
• Cosmological perturbations for two cold fluids in $\Lambda$CDM
• Semiclassical path to cosmic large-scale structure
• Large-Scale Velocity Dispersion and the Cosmic Web
• How Density Environment Changes the Influence of the Dark Matter-Baryon Streaming Velocity on the Cosmological Structure Formation
• Shell-crossing in a $\Lambda$CDM Universe
• Accurate initial conditions in mixed Dark Matter--Baryon simulations
• A hierarchical field-level inference approach to reconstruction from sparse Lyman-$\alpha$ forest data
• Toward Accurate Modeling of the Nonlinear Matter Bispectrum: Standard Perturbation Theory and Transients from Initial Conditions
• The nonlinear evolution of baryonic overdensities in the early universe: Initial conditions of numerical simulations
• Efficient parallel algorithm for estimating higher-order polyspectra
• The Borg Cube Simulation: Cosmological Hydrodynamics with CRK-SPH

Cited by in corpus (30)

• The FLAMINGO project: cosmological hydrodynamical simulations for large-scale structure and galaxy cluster surveys
• Large-scale dark matter simulations
• Hydrodynamical simulations of the galaxy population: enduring successes and outstanding challenges
• SWIFT: A modern highly-parallel gravity and smoothed particle hydrodynamics solver for astrophysical and cosmological applications
• Accurate initial conditions for cosmological N-body simulations: Minimizing truncation and discreteness errors
• Simulating intergalactic gas for DESI-like small scale Lymanα forest observations
• Higher-order initial conditions with massive neutrinos
• Cosmological perturbations for two cold fluids in $\Lambda$CDM
• SPICE: the connection between cosmic reionisation and stellar feedback in the first galaxies
• Shell-crossing in a $\Lambda$CDM Universe
• Cosmological constraints from the cross-correlation of DESI Luminous Red Galaxies with CMB lensing from Planck PR4 and ACT DR6
• The FLAMINGO project: Baryon effects on the matter power spectrum
• Simulating Hydrodynamics in Cosmology with CRK-HACC
• Rhapsody-C simulations -- Anisotropic thermal conduction, black hole physics, and the robustness of massive galaxy cluster scaling relations
• Cosmological feedback from a halo assembly perspective
• Analytical growth functions for cosmic structures in a $\Lambda$CDM Universe
• Particle initialization effects on Lyman-$\alpha$ forest statistics in cosmological SPH simulations
• The effect of local universe constraints on halo abundance and clustering
• Improving initialization and evolution accuracy of cosmological neutrino simulations
• Radiative transfer of Lyman-$\alpha$ photons at cosmic dawn with realistic gas physics
• COLIBRE: calibrating subgrid feedback in cosmological simulations that include a cold gas phase
• Numerical Discreteness Errors in Multi-Species Cosmological N-body Simulations
• Effect of Separate Initial Conditions on the Lyman-$\alpha$ Forest in Simulations
• Non-linear structure formation with elastic interactions in the dark sector
• The FLAMINGO simulations data release
• Baryonic effects in the Effective Field Theory of Large-Scale Structure and an analytic recipe for lensing in CMB-S4
• Parity Violation in Galaxy Shapes: Primordial Non-Gaussianity
• Cosmic voids and BAO with relative baryon-CDM perturbations
• Particle loads for cosmological simulations with equal-mass dark matter and baryonic particles
• Kinematic dynamos and resolution limits for Smoothed Particle Magnetohydrodynamics