Neutrino Density Research Articles

Neutrino halo profiles: HR-DEMNUni simulation analysis

Using the high-resolution HR-DEMNUni simulations, we computed neutrino profiles within virialized dark matter haloes. These new high-resolution simulations allowed us to revisit fitting formulas proposed in the literature and provided updated fitting parameters that extend to less massive haloes and lower neutrino masses than previously in the literature, in accordance with new cosmological limits. The trend we observe for low neutrino masses is that, for dark matter halo masses below ∼ 4 × 1014 h -1 M ⊙, the presence of the core becomes weaker and the profile over the whole radius is closer to a simple power law. We also characterised the neutrino density profile dependence on the solid angle within clustered structures: a forward-backward asymmetry larger than 10% was found when comparing the density profiles from neutrinos along the direction of motion of cold dark matter particles within the same halo. In addition, we looked for neutrino wakes around halo centres produced by the peculiar motion of the halo itself. Our results suggest that the wakes effect is present in haloes with masses greater than 3 × 1014 h -1 M ⊙ where a mean displacement of 0.06h -1Mpc was found.

Local clustering of relic neutrinos: comparison of kinetic field theory and the Vlasov equation

Gravitational clustering in our cosmic vicinity is expected to lead to an enhancement of the local density of relic neutrinos. We derive expressions for the neutrino density, using a perturbative approach to kinetic field theory and perturbative solutions of the Vlasov equation up to second order. Our work reveals that both formalisms give exactly the same results and can thus be considered equivalent. Numerical evaluation of the local relic neutrino density at first and second order provides some fundamental insights into the frequently applied approach of linear response to neutrino clustering (also known as the Gilbert equation). Against the naive expectation, including the second-order contribution does not lead to an improvement of the prediction for the local relic neutrino density but to a dramatic overestimation. This is because perturbation theory breaks down in a momentum-dependent fashion and in particular for densities well below unity.

Solar Models and Astrophysical S-factors Constrained by Helioseismic Results and Updated Neutrino Fluxes

The ratio of metal abundance to hydrogen abundance of the solar photosphere, (Z/X) s , has been revised several times. Standard solar models, based on these revised solar abundances, are in disagreement with seismically inferred results. Recently, Magg et al. introduced a new value for (Z/X) s , which is still under debate in the community. The solar abundance problem or solar modeling problem remains a topic of ongoing debate. We constructed rotating solar models in accordance with various abundance scales where the effects of convection overshoot and enhanced diffusion were included. Among these models, those utilizing Magg’s abundance scale exhibit superior sound speed and density profiles compared to models using other abundance scales. Additionally, they reproduce the observed frequency separation ratios r 02 and r 13. These models also match the seismically inferred surface helium abundance and convection zone depth within the 1σ level. Furthermore, the calculated neutrino fluxes from these models agree with detected ones at the level of 1σ. We found that neutrino fluxes and density profile are influenced by nuclear reactions, allowing us to use the combination of detected neutrino fluxes and seismically inferred density for diagnosing astrophysical S-factors. This diagnostic approach shows that S 11 may be underestimated by 2%, while S 33 may be overestimated by about 3% in previous determinations. The S-factors favored by updated neutrino fluxes and helioseismic results can lead to significant improvements in solar models.

Chaos in inhomogeneous neutrino fast flavor instability

In dense neutrino gases, the neutrino-neutrino coherent forward scattering gives rise to a complex flavor oscillation phenomenon not fully incorporated in simulations of neutron star mergers (NSM) and core collapse supernovae (CCSNe). Moreover, it has been proposed to be chaotic, potentially limiting our ability to predict neutrino flavor transformations in simulations. To address this issue, we explore how small flavor perturbations evolve in the nonlinear regime of the neutrino quantum kinetic equation within a narrow centimeter-scale region inside a NSM and a toy neutrino distribution. Our findings reveal that paths in the flavor state space of solutions with similar initial conditions diverge exponentially, exhibiting chaos. This inherent chaos makes the microscopic scales of neutrino flavor transformations unpredictable. However, the domain-averaged neutrino density matrix remains relatively stable, with chaos minimally affecting it. This particular property suggests that domain-averaged quantities remain reliable despite the exponential amplification of errors. Published by the American Physical Society 2024

Once-in-a-lifetime encounter models for neutrino media: From coherent oscillations to flavor equilibration

Collective neutrino oscillations are typically studied using the lowest-order quantum kinetic equation, also known as the mean-field approximation. However, some recent quantum many-body simulations suggest that quantum entanglement among neutrinos may be important and may result in flavor equilibration of the neutrino gas. In this work, we develop new quantum models for neutrino gases in which any pair of neutrinos can interact at most once in their lifetimes. A key parameter of our models is γ=μΔz, where μ is the neutrino coupling strength, which is proportional to the neutrino density, and Δz is the duration over which a pair of neutrinos can interact each time. Our models reduce to the mean-field approach in the limit γ→0 and achieve flavor equilibration in time t≫(γμ)−1. These models demonstrate the emergence of coherent flavor oscillations from the particle perspective and may help elucidate the role of quantum entanglement in collective neutrino oscillations. Published by the American Physical Society 2024

Cosmological constraints from density-split clustering in the BOSS CMASS galaxy sample

ABSTRACT We present a clustering analysis of the BOSS DR12 CMASS galaxy sample, combining measurements of the galaxy two-point correlation function and density-split clustering down to a scale of $1 \, h^{-1}\, \text{Mpc}$. Our theoretical framework is based on emulators trained on high-fidelity mock galaxy catalogues that forward model the cosmological dependence of the clustering statistics within an extended-ΛCDM framework, including redshift-space and Alcock–Paczynski distortions. Our base-ΛCDM analysis finds ωcdm = 0.1201 ± 0.0022, σ8 = 0.792 ± 0.034, and ns = 0.970 ± 0.018, corresponding to fσ8 = 0.462 ± 0.020 at z ≈ 0.525, which is in agreement with Planck 2018 predictions and various clustering studies in the literature. We test single-parameter extensions to base-ΛCDM, varying the running of the spectral index, the dark energy equation of state, and the density of mass-less relic neutrinos, finding no compelling evidence for deviations from the base model. We model the galaxy–halo connection using a halo occupation distribution framework, finding signatures of environment-based assembly bias in the data. We validate our pipeline against mock catalogues that match the clustering and selection properties of CMASS, showing that we can recover unbiased cosmological constraints even with a volume 84 times larger than the one used in this study.

Neutrino at Different Epochs of the Friedmann Universe

At least two relics of the Big Bang have survived: the cosmological microwave background (CMB) and the cosmological neutrino background (CνB). Being the second most abundant particle in the universe, the neutrino has a significant impact on its evolution from the Big Bang to the present day. Neutrinos affect the following cosmological processes: the expansion rate of the universe, its chemical and isotopic composition, the CMB anisotropy and the formation of the large-scale structure of the universe. Another relic neutrino background is theoretically predicted, it consists of non-equilibrium antineutrinos of Primordial Nucleosynthesis arising as a result of the decay of neutrons and tritium nuclei. Such antineutrinos are an indicator of the baryon asymmetry of the universe. In addition to experimentally detectable active neutrinos, the existence of sterile neutrinos is theoretically predicted to generate neutrino masses and explain their oscillations. Sterile neutrinos can also solve such cosmological problems as the baryonic asymmetry of the universe and the nature of dark matter. The recent results of several independent experiments point to the possibility of the existence of a light sterile neutrino. However, the existence of such a neutrino is inconsistent with the predictions of the Standard Cosmological Model. The inclusion of a non-zero lepton asymmetry of the universe and/or increasing the energy density of active neutrinos can eliminate these contradictions and reconcile the possible existence of sterile neutrinos with Primordial Nucleosynthesis, the CMB anisotropy, and also reduce the H0-tension. In this brief review, we discuss the influence of the physical properties of active and sterile neutrinos on the evolution of the universe from the Big Bang to the present day.

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

ABSTRACT Cosmology is poised to measure the neutrino mass sum Mν and has identified several smaller-scale observables sensitive to neutrinos, necessitating accurate predictions of neutrino clustering over a wide range of length scales. The FlowsForTheMasses non-linear perturbation theory for the the massive neutrino power spectrum, $\Delta ^2_\nu (k)$, agrees with its companion N-body simulation at the $10~{{\ \rm per\ cent}}-15~{{\ \rm per\ cent}}$ level for k ≤ 1 h Mpc−1. Building upon the Mira-Titan IV emulator for the cold matter, we use FlowsForTheMasses to construct an emulator for $\Delta ^2_\nu (k)$, Cosmic-Eν, which covers a large range of cosmological parameters and neutrino fractions Ων, 0h2 ≤ 0.01 (Mν ≤ 0.93 eV). Consistent with FlowsForTheMasses at the 3.5 per cent level, it returns a power spectrum in milliseconds. Ranking the neutrinos by initial momenta, we also emulate the power spectra of momentum deciles, providing information about their perturbed distribution function. Comparing a Mν = 0.15 eV model to a wide range of N-body simulation methods, we find agreement to 3 per cent for k ≤ 3kFS = 0.17 h Mpc−1 and to 19 per cent for k ≤ 0.4 h Mpc−1. We find that the enhancement factor, the ratio of $\Delta ^2_\nu (k)$ to its linear-response equivalent, is most strongly correlated with Ων, 0h2, and also with the clustering amplitude σ8. Furthermore, non-linearities enhance the free-streaming-limit scaling $\partial \log (\Delta ^2_\nu /\Delta ^2_{\rm m}) / \partial \log (M_\nu)$ beyond its linear value of 4, increasing the Mν-sensitivity of the small-scale neutrino density.

Quantum algorithm for the Vlasov simulation of the large-scale structure formation with massive neutrinos

Investigating the cosmological implication of the fact that neutrino has finite mass is of importance for fundamental physics. In particular, massive neutrino affects the formation of the large-scale structure (LSS) of the universe, and, conversely, observations of the LSS can give constraints on the neutrino mass. Numerical simulations of the LSS formation including massive neutrino along with conventional cold dark matter is thus an important task. For this, calculating the neutrino distribution in the phase space by solving the Vlasov equation is a suitable approach, but it requires solving the PDE in the (6+1)-dimensional space and is thus computationally demanding: Configuring ngr grid points in each coordinate and nt time grid points leads to O(ngr6) memory space and O(ntngr6) queries to the coefficients in the discretized PDE. We propose a quantum algorithm for this task. Linearizing the Vlasov equation by neglecting the relatively weak self-gravity of the neutrino, we perform the Hamiltonian simulation to produce quantum states that encode the phase-space distribution of neutrino. We also propose a way to extract the power spectrum of the neutrino density perturbations as classical data from the quantum state by quantum amplitude estimation with accuracy ε and query complexity of order Õ[(ngr+nt)/ε]. Our method also reduces the space complexity to O[polylog(ngr/ε)] in terms of the qubit number, while using quantum random access memories with O(ngr3) entries. As far as we know, this is the first quantum algorithm for the LSS simulation that outputs the quantity of practical interest with guaranteed accuracy. Published by the American Physical Society 2024

Light sterile neutrinos in the early universe: effects of altered dispersion relations and a coupling to axion-like dark matter

We investigate the cosmological consequences of light sterile neutrinos with altered dispersion relations (ADRs) and couplings to an ultra-light, axion-like scalar field. In particular we study the impact on the number of additional, light, fermionic degrees of freedom and primordial nucleosynthesis. While the ADR leads to a new potential term in the Hamiltonian, the coupling to the scalar field results in a time dependent, effective mass contribution. We solve the quantum kinetic equations (QKEs) for the neutrino density matrix and find that in certain parameter regions both new physics effects can individually yield a suppressed population of sterile neutrino species and the correct observed amount of helium in nucleosynthesis. Combining both effects opens up new patches of parameter space excluded by experimental bounds applying to models featuring only one of the effects.

Neutrino winds on the sky

We develop a first-principles formalism to compute the distortion to the relic neutrino density field caused by the peculiar motions of large-scale structures. This distortion slows halos down due to dynamical friction, causes a local anisotropy in the neutrino-CDM cross-correlation, and reduces the global cross-correlation between neutrinos and CDM. The local anisotropy in the neutrino-CDM cross-spectrum is imprinted in the three point cross-correlations of matter and galaxies, or the bispectrum in Fourier space, producing a signal peaking at squeezed triangle configurations. This bispectrum signature of neutrino masses is not limited by cosmic variance or potential inaccuracies in the modeling of complicated nonlinear and galaxy formation physics, and it is not degenerate with the optical depth to reionization. We show that future surveys have the potential to detect the distortion bispectrum.

Where shadows lie: reconstruction of anisotropies in the neutrino sky

The Cosmic Neutrino Background (CNB) encodes a wealth of information, but has not yet been observed directly. To determine the prospects of detection and to study its information content, we reconstruct the phase-space distribution of local relic neutrinos from the three-dimensional distribution of matter within 200 h -1 Mpc of the Milky Way. Our analysis relies on constrained realization simulations and forward modelling of the 2M++ galaxy catalogue. We find that the angular distribution of neutrinos is anti-correlated with the projected matter density, due to the capture and deflection of neutrinos by massive structures along the line of sight. Of relevance to tritium capture experiments, we find that the gravitational clustering effect of the large-scale structure on the local number density of neutrinos is more important than that of the Milky Way for neutrino masses less than 0.1 eV. Nevertheless, we predict that the density of relic neutrinos is close to the cosmic average, with a suppression or enhancement over the mean of (-0.3%, +7%, +27%) for masses of (0.01, 0.05, 0.1) eV. This implies no more than a marginal increase in the event rate for tritium capture experiments like PTOLEMY. We also predict that the CNB and CMB rest frames coincide for 0.01 eV neutrinos, but that neutrino velocities are significantly perturbed for masses larger than 0.05 eV. Regardless of mass, we find that the angle between the neutrino dipole and the ecliptic plane is small, implying a near-maximal annual modulation in the bulk velocity. Along with this paper, we publicly release our simulation data, comprising more than 100 simulations for six different neutrino masses.

Exploring ultralight scalar assistance in sterile neutrino dark matter: cold spectrum and unusual X/gamma-ray signatures

We present a scalar-driven sterile neutrino production model where the interaction with the ultralight scalar field modifies the oscillation production of sterile neutrinos in the early universe. The model effectively suppresses the production of sterile neutrinos at low temperatures due to the heavy scalar mass, resulting in a colder matter power spectrum that avoids constraints from small-scale structure observations. In this model, the dominant dark matter relic is from sterile neutrinos, with only a small fraction originating from the ultralight scalar. Furthermore, the model predicts a detectable X/γ-ray flux proportional to the cubic density of local sterile neutrinos for a light scalar mass due to the light scalar coupling to sterile neutrinos. This distinguishes our model from normal decaying dark matter, which has a linear dependence on the density. In addition, the model predicts a potential low-energy monochromatic neutrino signal that can be detectable by future neutrino telescopes.

Toward Powerful Probes of Neutrino Self-Interactions in Supernovae.

Neutrinos remain mysterious. As an example, enhanced self-interactions (νSI), which would have broad implications, are allowed. At the high neutrino densities within core-collapse supernovae, νSI should be important, but robust observables have been lacking. We show that νSI make neutrinos form a tightly coupled fluid that expands under relativistic hydrodynamics. The outflow becomes either a burst or a steady-state wind; which occurs here is uncertain. Though the diffusive environment where neutrinos are produced may make a wind more likely, further work is needed to determine when each case is realized. In the burst-outflow case, νSI increase the duration of the neutrino signal, and even a simple analysis of SN 1987A data has powerful sensitivity. For the wind-outflow case, we outline several promising ideas that may lead to new observables. Combined, these results are important steps toward solving the 35-year-old puzzle of how νSI affect supernovae.

Local clustering of relic neutrinos with kinetic field theory

The density of relic neutrinos is expected to be enhanced due to clustering in our local neighbourhood at Earth. We introduce a novel analytical technique to calculate the neutrino overdensity, based on kinetic field theory. Kinetic field theory is a particle-based theory for cosmic structure formation and in this work we apply it for the first time to massive neutrinos. The gravitational interaction is expanded in a perturbation series and we take into account the first-order contribution to the local density of relic neutrinos. For neutrino masses that are consistent with cosmological neutrino mass bounds, our results are in excellent agreement with state-of-the-art calculations.

Cross-correlation power spectra and cosmic birefringence of the CMB via photon-neutrino interaction

In the context of the standard model of particles, the weak interaction of cosmic microwave background (CMB) and cosmic neutrino background (CνB), can generate non-vanishing TB and EB power spectra in the order of one loop forward scattering, in the presence of scalar perturbation, which is in contrast with the standard scenario cosmology. Comparing our results with the current experimental data may provide, significant information about the nature of CνB, including CMB-CνB forward scattering for TB, TE, and EB power spectra. To this end, different cases were studied, including Majorana CνB and Dirac CνB. On the other hand, it was shown that the mean opacity due to cosmic neutrino background could behave as an anisotropic birefringent medium and change the linear polarization rotation angle. Considering the contributions from neutrino and anti-neutrino forward scattering with CMB photons (in the case of Dirac neutrino), we introduce relative neutrino and anti-neutrino density asymmetry (δν = Δnν /nν = nν -nν̅ /nν ). Then, using the cosmic birefringence angle reported by the Planck data release β = 0.30° ± 0.11° (68%C.L.), some constraints can be put on δν . Also, the value of cosmic birefringence due to Majorana CνB medium is estimated at about β| ν ≃ 0.2 rad. In this respect, since Majorana neutrino and anti-neutrino are exactly the same, both CB contributions will be added together. However, this value is at least two orders larger than the cosmic birefringence angle reported by the Planck data release, β = 0.30° ± 0.11° (68%C.L.). Finally, we shortly discussed this big inconsistency. It is noteworthy that to calculate the contribution of photon-neutrino forward scattering for cosmic birefringence, we just consider the standard model of particles and the standard scenario of cosmology.

Improving initialization and evolution accuracy of cosmological neutrino simulations

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 \U0001d4aa(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ν = 0.15 eV are in very good agreement with the expected level of errors in the emulator and in the direct N-body simulation.

Collision-induced flavor instability in dense neutrino gases with energy-dependent scattering

We investigate the collision-induced flavor instability in homogeneous, isotropic, dense neutrino gases in the two-flavor mixing scenario with energy-dependent scattering. We uncover a simple expression of the growth rate of this instability in terms of the flavor-decohering collision rates and the electron lepton number distribution of the neutrino. This growth rate is common to the neutrinos and antineutrinos of different energies, and is independent of the mass splitting and vacuum-mixing angle of the neutrino, the matter density, and the neutrino density, although the initial amplitude of the unstable oscillation mode can be suppressed by a large matter density. Our results suggest that neutrinos are likely to experience collision-induced flavor conversions deep inside a core-collapse supernova even when both the fast and slow collective flavor oscillations are suppressed.

A seesaw model for large neutrino masses in concordance with cosmology

Cosmological constraints on the sum of the neutrino masses can be relaxed if the number density of active neutrinos is reduced compared to the standard scenario, while at the same time keeping the effective number of neutrino species Neff ≈ 3 by introducing a new component of dark radiation. We discuss a UV complete model to realise this idea, which simultaneously provides neutrino masses via the seesaw mechanism. It is based on a U(1) symmetry in the dark sector, which can be either gauged or global. In addition to heavy seesaw neutrinos, we need to introduce \U0001d4aa(10) generations of massless sterile neutrinos providing the dark radiation. Then we can accommodate active neutrino masses with ∑mν ~ 1 eV, in the sensitivity range of the KATRIN experiment. We discuss the phenomenology of the model and identify the allowed parameter space. We argue that the gauged version of the model is preferred, and in this case the typical energy scale of the model is in the 10 MeV to few GeV range.

Time-Dependent and Quasisteady Features of Fast Neutrino-Flavor Conversion.

Despite the theoretical indication that fast neutrino-flavor conversion (FFC) ubiquitously occurs iin core-collapse supernovae and binary neutron star mergers, the lack of global simulations has been the greatest obstacle to study their astrophysical consequences. In this Letter, we present large-scale (50km) simulations of FFC in spherical symmetry by using a novel approach. We effectively rescale the oscillation scale of FFC by reducing the number of injected neutrinos in the simulation box, and then extrapolate back to the case of the target density of neutrinos with a convergence study. We find that FFC in all models achieves a quasisteady state in the nonlinear regime, and its saturation property of FFC is universal. We also find that temporal- and spatial variations of FFC are smeared out at large radii due to phase cancellation through neutrino self-interactions. Finally, we provide a new diagnostic quantity, electron neutrino lepton number subtracted by heavy one angular crossing, to assess the nonlinear saturation of FFC.