Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new user$ u$CO$N$CEPT: Cosmological neutrino simulations from the non-linear Boltzmann hierarchy
50
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
In this paper the non-linear effect of massive neutrinos on cosmological structures is studied in a conceptually new way. We have solved the non-linear continuity and Euler equations for the neutrinos on a grid in real space in $N$-body simulations, and closed the Boltzmann hierarchy at the non-linear Euler equation using the stress and pressure perturbations from linear theory. By comparing with state-of-the art cosmological neutrino simulations, we are able to simulate the non-linear neutrino power spectrum very accurately. This translates into a negligible error in the matter power spectrum, and so our CONCEPT code is ideally suited for extracting the neutrino mass from future high precision non-linear observational probes such as Euclid.
rate research
Read More
Calculations of the evolution of cosmological perturbations generally involve solution of a large number of coupled differential equations to describe the evolution of the multipole moments of the distribution of photon intensities and polarization. However, this Boltzmann hierarchy communicates with the rest of the system of equations for the other perturbation variables only through the photon-intensity quadrupole moment. Here I develop an alternative formulation wherein this photon-intensity quadrupole is obtained via solution of two coupled integral equations -- one for the intensity quadrupole and another for the linear-polarization quadrupole -- rather than the full Boltzmann hierarchy. This alternative method of calculation provides some physical insight and a cross-check for the traditional approach. I describe a simple and efficient iterative numerical solution that converges fairly quickly. I surmise that this may allow current state-of-the-art cosmological-perturbation codes to be accelerated.
The impact of dark matter-neutrino interactions on the measurement of the cosmological parameters has been investigated in the past in the context of massless neutrinos exclusively. Here we revisit the role of a neutrino-dark matter coupling in light of ongoing cosmological tensions by implementing the full Boltzmann hierarchy for three massive neutrinos. Our tightest 95% CL upper limit on the strength of the interactions, parameterized via $u_chi =frac{sigma_0}{sigma_{Th}}left(frac{m_chi}{100 text{GeV}}right)^{-1}$, is $u_chileq3.34 cdot 10^{-4}$, arising from a combination of Planck TTTEEE data, Planck lensing data and SDSS BAO data. This upper bound is, as expected, slightly higher than previous results for interacting massless neutrinos, due to the correction factor associated with neutrino masses. We find that these interactions significantly relax the lower bounds on the value of $sigma_8$ that is inferred in the context of $Lambda$CDM from the Planck data, leading to agreement within 1-2$sigma$ with weak lensing estimates of $sigma_8$, as those from KiDS-1000. However, the presence of these interactions barely affects the value of the Hubble constant $H_0$.
The non-zero mass of neutrinos suppresses the growth of cosmic structure on small scales. Since the level of suppression depends on the sum of the masses of the three active neutrino species, the evolution of large-scale structure is a promising tool to constrain the total mass of neutrinos and possibly shed light on the mass hierarchy. In this work, we investigate these effects via a large suite of N-body simulations that include massive neutrinos using an analytic linear-response approximation: the Cosmological Massive Neutrino Simulations (MassiveNuS). The simulations include the effects of radiation on the background expansion, as well as the clustering of neutrinos in response to the nonlinear dark matter evolution. We allow three cosmological parameters to vary: the neutrino mass sum M_nu in the range of 0-0.6 eV, the total matter density Omega_m, and the primordial power spectrum amplitude A_s. The rms density fluctuation in spheres of 8 comoving Mpc/h (sigma_8) is a derived parameter as a result. Our data products include N-body snapshots, halo catalogues, merger trees, ray- traced galaxy lensing convergence maps for four source redshift planes between z_s=1-2.5, and ray-traced cosmic microwave background lensing convergence maps. We describe the simulation procedures and code validation in this paper. The data are publicly available at http://columbialensing.org.
We study in detail how neutrino perturbations can be followed in linear theory by using only terms up to $l=2$ in the Boltzmann hierarchy. We provide a new approximation to the third moment and demonstrate that the neutrino power spectrum can be calculated to a precision of better than $sim$ 5% for masses up to $sim$ 1 eV. The matter and CMB power spectra can be calculated far more precisely and typically at least a factor of a few better than with existing approximations. We then proceed to study how the neutrino power spectrum can be reliably calculated even in the presence of non-linear gravitational clustering by using the full non-linear gravitational potential derived from semi-analytic methods based on $N$-body simulations in the Boltzmann evolution hierarchy. Our results agree extremely well with results derived from $N$-body simulations.
Constraining neutrino mass remains an elusive challenge in modern physics. Precision measurements are expected from several upcoming cosmological probes of large-scale structure. Achieving this goal relies on an equal level of precision from theoretical predictions of neutrino clustering. Numerical simulations of the non-linear evolution of cold dark matter and neutrinos play a pivotal role in this process. We incorporate neutrinos into the cosmological N-body code CUBEP3M and discuss the challenges associated with pushing to the extreme scales demanded by the neutrino problem. We highlight code optimizations made to exploit modern high performance computing architectures and present a novel method of data compression that reduces the phase-space particle footprint from 24 bytes in single precision to roughly 9 bytes. We scale the neutrino problem to the Tianhe-2 supercomputer and provide details of our production run, named TianNu, which uses 86% of the machine (13,824 compute nodes). With a total of 2.97 trillion particles, TianNu is currently the worlds largest cosmological N-body simulation and improves upon previous neutrino simulations by two orders of magnitude in scale. We finish with a discussion of the unanticipated computational challenges that were encountered during the TianNu runtime.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
