Abstract:Neutrino mass constraints are a primary focus of current and future large-scale structure (LSS) surveys. Non-linear LSS models rely heavily on cosmological simulations -- the impact of massive neutrinos should therefore be included in these simulations in a realistic, computationally tractable, and controlled manner. A recent proposal to reduce the related computational cost employs a symmetric neutrino momentum sampling strategy in the initial conditions. We implement a modified version of this strategy into the Hardware/Hybrid Accelerated Cosmology Code (HACC) and perform convergence tests on its internal parameters. We illustrate that this method can impart $\mathcal{O}(1\%)$ numerical artifacts on the total matter field on small scales, similar to previous findings, and present a method to remove these artifacts using Fourier-space filtering of the neutrino density field. Moreover, we show that the converged neutrino power spectrum does not follow linear theory predictions on relatively large scales at early times at the $15\%$ level, prompting a more careful study of systematics in particle-based neutrino simulations. We also present an improved method for backscaling linear transfer functions for initial conditions in massive neutrino cosmologies that is based on achieving the same relative neutrino growth as computed with Boltzmann solvers. Our self-consistent backscaling method yields sub-percent accuracy in the total matter growth function. Comparisons for the non-linear power spectrum with the Mira-Titan emulator at a neutrino mass of $m_{\nu}=0.15~\mathrm{eV}$ are in very good agreement with the expected level of errors in the emulator and in the direct N-body simulation.

What problem does this paper attempt to address?

### Problems the paper attempts to solve This paper mainly focuses on how to more accurately initialize and evolve systems containing massive neutrinos in cosmological simulations. Specifically, the main challenges faced by the researchers include: 1. **Reducing numerical noise**: The high - thermal motion of massive neutrinos in the early universe leads to a significant "shot noise" effect, which appears as noise in the power spectrum. To reduce this noise, the researchers introduced a symmetric momentum sampling strategy, the "tiling" method, which reduces noise by symmetrically sampling the neutrino momentum under initial conditions. 2. **Eliminating numerical artifacts**: Although the "tiling" method can effectively reduce noise, it may also introduce new numerical artifacts, especially on small scales. The paper proposes a method based on Fourier - space filtering to eliminate these artifacts. 3. **Improving the background scaling method**: When constructing initial conditions, it is necessary to trace the linear theory transfer function from the initial redshift back to the final redshift to match the observed large - scale linear matter power spectrum. Traditional background scaling methods have deficiencies when dealing with massive neutrinos because they assume that neutrinos are a fluid. The paper proposes an iterative background scaling method that more accurately sets the initial conditions by directly using the output of the Boltzmann solver to calculate the growth of neutrinos relative to cold matter. ### Main contributions - **Symmetric momentum sampling strategy**: By symmetrically sampling the neutrino momentum under initial conditions, the "shot noise" is significantly reduced, allowing simulations to be carried out with fewer neutrino particle trackers. - **Fourier - space filtering**: A method is proposed to eliminate the numerical artifacts introduced by the "tiling" method by filtering the neutrino density field in Fourier - space. - **Iterative background scaling**: An iterative background scaling method has been developed that can more accurately set the initial conditions, especially when dealing with massive neutrinos, avoiding the errors of traditional methods. ### Conclusion Through these improvements, the researchers are able to achieve higher precision and reliability in cosmological simulations containing massive neutrinos. These results are of great significance for understanding the neutrino mass and the evolution of large - scale structures, and are helpful for future cosmological observations and experiments.

Improving initialization and evolution accuracy of cosmological neutrino simulations

An efficient implementation of massive neutrinos in non-linear structure formation simulations

Cosmic-Eν: An emulator for the non-linear neutrino power spectrum

Hybrid multi-fluid-particle simulations of the cosmic neutrino background

An accurate fluid approximation for massive neutrinos in cosmology

Learning neutrino effects in Cosmology with Convolutional Neural Networks

How new physics affects primordial neutrinos decoupling: Direct Simulation Monte Carlo approach

Longitudinal kinematic imbalances in neutrino and antineutrino interactions for improved measurements of neutrino energy and the axial vector form factor

Rapid Methods for Modeling Overdensities of Massive Neutrinos and Other Non-Cold Relics

Beware of Fake $ν$s: The Effect of Massive Neutrinos on the Nonlinear Evolution of Cosmic Structure

Revisiting coupled CDM-massive neutrino perturbations in diverse cosmological backgrounds

Leveraging neutrino flavor physics for supernova model differentiation

The Neutrino Puzzle: Anomalies, Interactions, and Cosmological Tensions

Euclid: Modelling massive neutrinos in cosmology -- a code comparison

Effects of primordial fluctuations on relic neutrino simulations

Aemulus $ν$: Precise Predictions for Matter and Biased Tracer Power Spectra in the Presence of Neutrinos

From ray to spray: augmenting amplitudes and taming fast oscillations in fully numerical neutrino codes

Fast Exact Algorithm for Neutrino Oscillation in Constant Matter Density

Modeling neutrino-induced scale-dependent galaxy clustering for photometric galaxy surveys

Measuring Oscillations with A Million Atmospheric Neutrinos

syren-new: Precise formulae for the linear and nonlinear matter power spectra with massive neutrinos and dynamical dark energy