The abacus cosmological N-body code Garrison, Lehman H; Eisenstein, Daniel J; Ferrer, Douglas ...
Monthly notices of the Royal Astronomical Society,
11/2021, Volume:
508, Issue:
1
Journal Article
Peer reviewed
Open access
ABSTRACT
We present abacus, a fast and accurate cosmological N-body code based on a new method for calculating the gravitational potential from a static multipole mesh. The method analytically ...separates the near- and far-field forces, reducing the former to direct 1/r2 summation and the latter to a discrete convolution over multipoles. The method achieves 70 million particle updates per second per node of the Summit supercomputer, while maintaining a median fractional force error of 10−5. We express the simulation time-step as an event-driven ‘pipeline’, incorporating asynchronous events such as completion of co-processor work, input/output, and network communication. abacus has been used to produce the largest suite of N-body simulations to date, the abacussummit suite of 60 trillion particles, incorporating on-the-fly halo finding. abacus enables the production of mock catalogues of the volume and resolution required by the coming generation of cosmological surveys.
Abstract
We present a high-fidelity realization of the cosmological N-body simulation from the Schneider et al. code comparison project. The simulation was performed with our AbacusN-body code, which ...offers high-force accuracy, high performance, and minimal particle integration errors. The simulation consists of 20483 particles in a $500\ h^{-1}\, \mathrm{Mpc}$ box for a particle mass of $1.2\times 10^9\ h^{-1}\, \mathrm{M}_\odot$ with $10\ h^{-1}\, \mathrm{kpc}$ spline softening. Abacus executed 1052 global time-steps to z = 0 in 107 h on one dual-Xeon, dual-GPU node, for a mean rate of 23 million particles per second per step. We find Abacus is in good agreement with Ramses and Pkdgrav3 and less so with Gadget3. We validate our choice of time-step by halving the step size and find sub-percent differences in the power spectrum and 2PCF at nearly all measured scales, with ${\lt }0.3{{\ \rm per\ cent}}$ errors at $k\lt 10\ \mathrm{Mpc}^{-1}\, h$. On large scales, Abacus reproduces linear theory better than 0.01 per cent. Simulation snapshots are available at http://nbody.rc.fas.harvard.edu/public/S2016.
We present a public data release of halo catalogs from a suite of 125 cosmological N-body simulations from the Abacus project. The simulations span 40 wCDM cosmologies centered on the Planck 2015 ...cosmology at two mass resolutions, 4 × 1010 h−1 M and 1 × 1010 h−1 M , in 1.1 h−1 Gpc and 720 h−1 Mpc boxes, respectively. The boxes are phase-matched to suppress sample variance and isolate cosmology dependence. Additional volume is available via 16 boxes of fixed cosmology and varied phase; a few boxes of single-parameter excursions from Planck 2015 are also provided. Catalogs spanning z = 1.5 to 0.1 are available for friends-of-friends and Rockstar halo finders and include particle subsamples. All data products are available at https://lgarrison.github.io/AbacusCosmos.
In cosmological N-body simulations, the representation of dark matter as discrete ‘macroparticles’ suppresses the growth of structure, such that simulations no longer reproduce linear theory on small ...scales near k
Nyquist. Marcos et al. demonstrate that this is due to sparse sampling of modes near k
Nyquist and that the often-assumed continuum growing modes are not proper growing modes of the particle system. We develop initial conditions (ICs) that respect the particle linear theory growing modes and then rescale the mode amplitudes to account for growth suppression. These ICs also allow us to take advantage of our very accurate N-body code abacus to implement second-order Lagrangian perturbation theory (2LPT) in configuration space. The combination of 2LPT and rescaling improves the accuracy of the late-time power spectra, halo mass functions, and halo clustering. In particular, we achieve 1 per cent accuracy in the power spectrum down to k
Nyquist, versus k
Nyquist/4 without rescaling or k
Nyquist/13 without 2LPT, relative to an oversampled reference simulation. We anticipate that our 2LPT will be useful for large simulations where fast Fourier transforms are expensive and that rescaling will be useful for suites of medium-resolution simulations used in cosmic emulators and galaxy survey mock catalogues. Code to generate ICs is available at https://github.com/lgarrison/zeldovich-PLT.
We examine the nature of the opacity and radiation transport in Type Ia supernovae. The dominant opacity arises from line transitions. We discuss the nature of line opacities and diffusion in ...expanding media and the appropriateness of various mean and expansion opacities used in light-curve calculations. Fluorescence is shown to be the dominant physical process governing the rate at which energy escapes the supernova. We present a sample light curve that was obtained using a time-dependent solution of the radiative transport equation with a spectral resolution of 80 km s-1 and employing an LTE equation of state. The result compares favorably with light curves and spectra of typical supernovae and is used to illustrate the physics controlling the evolution of the light curve and especially the secondary maxima seen in infrared photometry. (c) 2000 The American Astronomical Society.
We develop an analytic solution of the radiation transport problem for Type Ia supernovae (SNe Ia) and show that it reproduces bolometric light curves produced by more detailed calculations under the ...assumption of a constant-extinction coefficient. This model is used to derive the thermal conditions in the interior of SNe Ia and to study the sensitivity of light curves to various properties of the underlying supernova explosions. Although the model is limited by simplifying assumptions, it is adequate for demonstrating that the relationship between SNe Ia maximum-light luminosity and rate of decline is most easily explained if SNe Ia span a range in mass. The analytic model is also used to examine the size of various terms in the transport equation under conditions appropriate to maximum light. For instance, the Eulerian and advective time derivatives are each shown to be of the same order of magnitude as other order v/c terms in the transport equation. We conclude that a fully time-dependent solution to the transport problem is needed in order to compute SNe Ia light curves and spectra accurate enough to distinguish subtle differences of various explosion models. (c) 2000 The American Astronomical Society.
Using high-resolution cosmological simulations, we study hydrogen and helium gravitational cooling radiation from gas accretion by young galaxies. We focus on the He II cooling lines, which arise ...from gas with a different temperature history (T sub(max) 6 10 super(5) K) than H I line-emitting gas. We examine whether three major atomic cooling lines, H I l1216, He II l1640, and He II l304, are observable, finding that Ly a and He II l1640 cooling emission at z = 2-3 are potentially detectable with deep narrowband (R > 100) imaging and/or spectroscopy from the ground. While the expected strength of H I l1216 cooling emission depends strongly on the treatment of the self-shielded phase of the IGM in the simulations, our predictions for the He II ll640 line are more robust, because the He II emissivity is negligible below T 6 10 super(4.5) K and less sensitive to the UV background. Although He II ll640 cooling emission is fainter than Lya by at least a factor of 10 and, unlike Lya, might not be resolved spatially with current observational facilities, it is more suitable to study gas accretion in the galaxy formation process because it is optically thin and less contaminated by the recombination lines from star-forming galaxies. The He II ll640 line can be used to distinguish among mechanisms for powering the so-called Lya blobs - including gravitational cooling radiation, photoionization by stellar populations, and starburst-driven superwinds -because (1) He II ll640 emission is limited to very low metallicity log(Z/Z sub(z))< -5.3 and Population III stars and (2) the blob's kinematics are probed unambiguously through the He II line width, which for cooling radiation is narrower (s < 400 km s super(-1)) than typical wind speeds.