2018/06/15 by Jacob Brandbyge, Steen Hannestad, Thomas Tram · 12 citations
Computer Science · Physics and Astronomy · #Computational Physics and Python Applications #Cosmic neutrino background #Dark Matter and Cosmic Phenomena #Distribution (mathematics) #Momentum (technical analysis) #Neutrino #Particle physics theoretical and experimental studies #Position and momentum space #Redshift #Spectral density #Spurious relationship #astro-ph.CO
paper · pdf · doi:10.1088/1475-7516/2019/03/047
published in Journal of Cosmology and Astroparticle Physics 2019(03), 047 (Institute of Physics) · 14 pages, 8 figures
arxiv created 2018/06/15 · openalex created_date 2018/06/21 · openalex publication_date 2019/03/29 · arxiv updated 2019/04/10 · openalex updated_date 2026/08/05
Including massive neutrinos in N -body simulations is a challenging task due to the large thermal velocities of the neutrinos. In particle based codes this leads to problems of shot-noise due to insufficient sampling of the neutrino momentum distribution function. In this paper we investigate the benefits and drawbacks of a scheme first suggested in a paper by Banerjee et al. [1] in which the initial neutrino distribution is symmetrised in momentum space. Using a somewhat different, but qualitatively similar implementation, we confirm that this method reduces shot-noise significantly, but we also find that it generates some spurious power in the neutrino power spectrum at intermediate and small scales. We speculate that this happens because many neutrinos in the simulation sample the same underlying dark matter structures while moving through the simulation. By carefully tuning the number of directions in momentum space of the initial neutrino distribution we show that some improvements can be made over the case where initial neutrino directions are purely random. At redshifts z ≳ 3 the method works very well, but at smaller redshifts significant improvements are not possible due to the spurious power generation.