At early times dark matter has a thermal velocity dispersion of unknown amplitude which, for warm dark matter models, can influence the formation of nonlinear structure on observable scales. We propose a new scheme to simulate cosmologies with a small-scale suppression of perturbations that combines two previous methods in a way that avoids the numerical artefacts which have so far prevented either from producing fully reliable results. At low densities and throughout most of the cosmological volume, we represent the dark matter phase-sheet directly using high-accuracy interpolation, thereby avoiding the artificial fragmentation which afflicts particle-based methods in this regime. Such phase-sheet methods are, however, unable to follow the rapidly increasing complexity of the denser regions of dark matter haloes, so for these we switch to an N-body scheme which uses the geodesic deviation equation to track phase-sheet properties local to each particle. In addition, we present a novel high-resolution force calculation scheme based on an oct-tree of cubic force resolution elements which is well suited to approximate the force-field of our combined sheet+particle distribution. Our hybrid simulation scheme enables the first reliable simulations of the internal structure of low-mass haloes in a warm dark matter cosmology.

We present the BACCO project, a simulation framework specially designed to provide highly-accurate predictions for the distribution of mass, galaxies, and gas as a function of cosmological parameters. In this paper, we describe our main suite of simulations (L $\sim2$ Gpc and $4320^3$ particles) and present various validation tests. Using a cosmology-rescaling technique, we predict the nonlinear mass power spectrum over the redshift range 0 < z < 1.5 and over scales $10^{-2} < k/(h Mpc^{-1} ) < 5$ for 800 points in an 8-dimensional cosmological parameter space. For an efficient interpolation of the results, we build a Gaussian emulator which we test against other widely-used methods. Over the whole range of scales considered, we expect our predictions to be accurate at the 2% level for parameters in the minimal ${\Lambda}$CDM model and to 3% when extended to dynamical dark energy and massive neutrinos. We make our emulator publicly available under http://www.dipc.org/bacco

We present anisotropic ``separate universe'' simulations which modify the N-body code Gadget-4 in order to represent a large-scale tidal field through an anisotropic expansion factor. These simulations are used to measure the linear, quasi-linear and nonlinear response of the matter power spectrum to a spatially uniform trace-free tidal field up to wavenumber $k = 10 h \text{Mpc}^{-1}$. Together with the response to a large-scale overdensity measured in previous work, this completely describes the nonlinear matter bispectrum in the squeezed limit. We find that the response amplitude does not approach zero on small scales in physical coordinates, but rather a constant value, $R_K\approx 0.5$ for $k \geq 3 h\text{Mpc}^{-1}$ at $z=0$. This shows that even the inner regions of haloes are affected by the large-scale tidal field. We also measure directly the alignment of halo shapes with the tidal field, finding a clear signal which increases with halo mass.

The interpretation of cosmological observables requires the use of increasingly sophisticated theoretical models. Since these models are becoming computationally very expensive and display non-trivial uncertainties, the use of standard Bayesian algorithms for cosmological inferences, such as MCMC, might become inadequate. Here, we propose a new approach to parameter estimation based on an iterative Gaussian emulation of the target likelihood function. This requires a minimal number of likelihood evaluations and naturally accommodates for stochasticity in theoretical models. We apply the algorithm to estimate 9 parameters from the monopole and quadrupole of a mock power spectrum in redshift space. We obtain accurate posterior distribution functions with approximately 100 times fewer likelihood evaluations than an affine invariant MCMC, roughly independently from the dimensionality of the problem. We anticipate that our parameter estimation algorithm will accelerate the adoption of more accurate theoretical models in data analysis, enabling more comprehensive exploitation of cosmological observables.

In this paper we carry out anisotropic "separate universe" simulations by including a large-scale tidal field in the N-body code \textsc{gadget}4 using an anisotropic expansion factor $A_{ij}$. We use the code in a pure \textit{particle-mesh} (PM) mode to simulate the evolution of 16 realizations of an initial density field with and without a large-scale tidal field, which are then used to measure the \textit{response function} describing how the tidal field influences structure formation in the linear and non-linear regimes. Together with the previously measured response to a large scale overdensity, this completely describes the nonlinear matter bispectrum in the squeezed limit. We find that, contrary to the density response, the tidal response never significantly exceeds the large-scale perturbation-theory prediction even on nonlinear scales for the redshift range we discuss. We develop a simple halo model that takes into account the effect of the tidal field and compare it with our direct measurement from the anisotropic N-body simulations.

Despite the fact that the mean matter density of the universe has been measured to an accuracy of a few percent within the standard $\Lambda$CDM paradigm, its median density is not known even to order of magnitude. Typical points lie in low-density regions and are not part of a collapsed structure of any scale. Locally, the dark matter distribution is then simply a stretched version of that in the early universe. In this single-stream regime, the distribution of unsmoothed density is sensitive to the initial power spectrum on all scales, in particular on very small scales, and hence to the nature of the dark matter. It cannot be estimated reliably using conventional cosmological simulations because of the enormous dynamic range involved, but a suitable excursion set procedure can be used instead. For the Planck cosmological parameters, a 100 GeV WIMP, corresponding to a free-streaming mass $\sim 10^{-6}M_\odot$, results in a median density of $\sim 4\times 10^{-3}$ in units of the mean density, whereas a 10 $\mu$eV axion with free-streaming mass $\sim 10^{-12}M_\odot$ gives $\sim 3\times 10^{-3}$, and Warm Dark Matter with a (thermal relic) mass of 1 keV gives $\sim 8\times 10^{-2}$. In CDM (but not in WDM) universes, single-stream regions are predicted to be topologically isolated by the excursion set formalism. A test by direct N-Body simulations seems to confirm this prediction, although it is still subject to finite size and resolution effects. Unfortunately, it is unlikely that any of these properties is observable and so suitable for constraining the properties of dark matter.