Total Neutrino Research Articles

Core-collapse supernova neutrino emission and detection informed by state-of-the-art three-dimensional numerical models

ABSTRACT Based on our recent three-dimensional core-collapse supernova (CCSN) simulations including both exploding and non-exploding models, we study the detailed neutrino signals in representative terrestrial neutrino observatories, namely Super-Kamiokande (Hyper-Kamiokande), DUNE, JUNO, and IceCube. We find that the physical origin of difference in the neutrino signals between 1D and 3D is mainly proto-neutron-star convection. We study the temporal and angular variations of the neutrino signals and discuss the detectability of the time variations driven by the spiral standing accretion shock instability (spiral SASI) when it emerges for non-exploding models. In addition, we determine that there can be a large angular asymmetry in the event rate (${\gtrsim} 50 {{\ \rm per\ cent}}$), but the time-integrated signal has a relatively modest asymmetry (${\lesssim} 20 {{\ \rm per\ cent}}$). Both features are associated with the lepton-number emission self-sustained asymmetry and the spiral SASI. Moreover, our analysis suggests that there is an interesting correlation between the total neutrino energy (TONE) and the cumulative number of neutrino events in each detector, a correlation that can facilitate data analyses of real observations. We demonstrate the retrieval of neutrino energy spectra for all flavours of neutrino by applying a novel spectrum reconstruction technique to the data from multiple detectors. We find that this new method is capable of estimating the TONE within the error of ∼20 per cent if the distance to the CCSN is ≲6 kpc.

Update constraints on neutrino mass and mass hierarchy in light of dark energy models

Combining cosmic microwave (CMB) background data from Planck satellite data, Baryon Acoustic Oscillations (BAO) measurements and Type Ia supernovae (SNe Ia) data, we obtain the bounds on total neutrino masses [Formula: see text] with the approximation of degenerate neutrino masses and for three dark energy models: the cosmological constant ([Formula: see text]CDM) model, a phenomenological emergent dark energy (PEDE) model and a model-independent quintessential parametrization (HBK). The bounds on the sum of neutrino masses [Formula: see text] depend on the dark energy (DE) models. In the HBK model, we confirm the conclusion from some previous work that the quintessence prior of DE tends to tighten the cosmological constraint on [Formula: see text]. On the other hand, the PEDE model leads to larger [Formula: see text] and a nonzero lower bound. Besides, we also explore the correlation between three different neutrino hierarchies and DE models.

What will it take to measure individual neutrino mass states using cosmology?

We study the impact of assumptions made about the neutrino mass ordering on cosmological parameter estimation with the purpose of understanding whether in the future it will be possible to infer the specific neutrino mass distribution from cosmological data. We find that although the commonly used assumption of a degenerate neutrino hierarchy is manifestly wrong and leads to changes in cosmological observables such as the cosmic microwave background and large scale structure compared to the correct (normal or inverted) neutrino hierarchy, the induced changes are so small that even with extremely optimistic assumptions about future data they will remain undetectable. We are thus able to conclude that while cosmology can probe the neutrino contribution to the cosmic energy density extremely precisely (and hence provide a detection of a non-zero total neutrino mass at high significance), it will not be possible to directly measure the individual neutrino masses.

A precision calculation of relic neutrino decoupling

We study the distortions of equilibrium spectra of relic neutrinos due to the interactions with electrons, positrons, and neutrinos in the early Universe. We solve the integro-differential kinetic equations for the neutrino density matrix, including three-flavor oscillations and finite temperature corrections from QED up to the next-to-leading order O(e3) for the first time. In addition, the equivalent kinetic equations in the mass basis of neutrinos are directly solved, and we numerically evaluate the distortions of the neutrino spectra in the mass basis as well, which can be easily extrapolated into those for non-relativistic neutrinos in the current Universe. In both bases, we find the same value of the effective number of neutrinos, Neff=3.044, which parameterizes the total neutrino energy density. The estimated error for the value of Neff due to the numerical calculations and the choice of neutrino mixing parameters would be at most 0.0005.

IceCube Search for High-energy Neutrino Emission from TeV Pulsar Wind Nebulae

Abstract Pulsar wind nebulae (PWNe) are the main gamma-ray emitters in the Galactic plane. They are diffuse nebulae that emit nonthermal radiation. Pulsar winds, relativistic magnetized outflows from the central star, shocked in the ambient medium produce a multiwavelength emission from the radio through gamma-rays. Although the leptonic scenario is able to explain most PWNe emission, a hadronic contribution cannot be excluded. A possible hadronic contribution to the high-energy gamma-ray emission inevitably leads to the production of neutrinos. Using 9.5 yr of all-sky IceCube data, we report results from a stacking analysis to search for neutrino emission from 35 PWNe that are high-energy gamma-ray emitters. In the absence of any significant correlation, we set upper limits on the total neutrino emission from those PWNe and constraints on hadronic spectral components.

Dynamical dark sectors and neutrino masses and abundances

We investigate generalized interacting dark matter-dark energy scenarios with a time-dependent coupling parameter, allowing also for freedom in the neutrino sector. The models are tested in the phantom and quintessence regimes, characterized by an equation of state $w_x<-1$ and $w_x>-1$, respectively. Our analyses show that for some of the scenarios the existing tensions on the Hubble constant $H_0$ and on the clustering parameter $S_8$ can be significantly alleviated. The relief is either due to \textit{(a)} a dark energy component which lies within the phantom region; or \textit{(b)} the presence of a dynamical coupling in quintessence scenarios. The inclusion of massive neutrinos into the interaction schemes does not affect neither the constraints on the cosmological parameters nor the bounds on the total number or relativistic degrees of freedom $N_{\rm eff}$, which are found to be extremely robust and, in general, strongly consistent with the canonical prediction $N_{\rm eff}=3.045$. The most stringent bound on the total neutrino mass $M_{\nu}$ is $M_{\nu}<0.116$ eV and it is obtained within a quintessence scenario in which the matter mass-energy density is only mildly affected by the presence of a dynamical dark sector coupling.

Including massive neutrinos in thermal Sunyaev Zeldovich power spectrum and cluster counts analyses

ABSTRACT We consistently include the effect of massive neutrinos in the thermal Sunyaev Zeldovich (SZ) power spectrum and cluster counts analyses, highlighting subtle dependencies on the total neutrino mass and data combination. In particular, we find that using the transfer functions for cold dark matter (CDM) + baryons in the computation of the halo mass function, instead of the transfer functions including neutrino perturbations, as prescribed in recent work, yields an ≈0.25 per cent downward shift of the σ8 constraint from tSZ power spectrum data, with a fiducial neutrino mass Σmν = 0.06 eV. In ΛCDM, with an X-ray mass bias corresponding to the expected hydrostatic mass bias, i.e. (1 − b) ≃ 0.8, our constraints from Planck SZ data are consistent with the latest results from SPT, DES-Y1, and KiDS+VIKING-450. In νΛCDM, our joint analyses of Planck SZ with Planck 2015 primary CMB yield a small improvement on the total neutrino mass bound compared to the Planck 2015 primary CMB constraint, as well as (1 − b) = 0.64 ± 0.04 (68 per cent CL). For forecasts, we find that competitive neutrino mass measurements using cosmic variance limited SZ power spectrum require masking the heaviest clusters and probing the small-scale SZ power spectrum up to ℓmax ≈ 104. Although this is challenging, we find that SZ power spectrum can realistically be used to tightly constrain intracluster medium properties: we forecast a 2 per cent determination of the X-ray mass bias by combining CMB-S4 and our mock SZ power spectrum with ℓmax = 103.

Constraints on the spacetime dynamics of an early dark energy component

We consider an Early Dark Energy (EDE) cosmological model, and perform an analysis which takes into account both background and perturbation effects via the parameters c2eff and c2vis, representing effective sound speed and viscosity, respectively. By using the latest available data we derive constraints on the amount of dark energy at early times and the present value of the equation of state. Our focus is on the effect that early dark energy has on the Cosmic Microwave Background (CMB) data, including polarization and lensing, in a generalized parameter space including a varying total neutrino mass, and tensor to scalar ratio, besides the 6 standard parameters of the minimal cosmological model. We find that the inclusion of Baryonic Acoustic Oscillations (BAO) data and CMB lensing significantly improves the constraints on the EDE parameters, while other high redshift data like the Quasar Hubble diagram and the Lyman-α forest BAO have instead a negligible impact. We find ΩeDE < 0.0039 and w0 < −0.95 at the 95 % C.L. for EDE accounting for its clustering through the inclusion of perturbation dynamics. This limit becomes weaker ΩeDE < 0.0034 if perturbations are neglected. The constraints on the EDE parameters are remarkably stable even when Σ mν, and r parameters are varied, with weak degeneracies between ΩeDE and r or Σ mν. In general we expect smaller values for the upper limits on the total amount of EDE with an increasing neutrino mass, while with a decreasing value of the tensor to scalar ratio we expect the 2σ upper limits on EDE to increase. We compare this EDE model with a simple wCDM with zero dark energy at early times and we find ∼1–2% different upper limits on total neutrino mass and ∼0.1–0.2% difference on the equation of state at the present time. Perturbation parameters are not constrained with current data sets, and tensions between the CMB derived H0 and σ8 values and those measured with local probes are not eased. This work demonstrates the capability of CMB probes to constrain the total amount of EDE well below the percent level.

The Effect of Massive Neutrinos on the Halo Spin Flip Phenomenon

Abstract The halo spin flip refers to the phenomenon that the spin axes of dark matter halos with masses above a certain threshold tend to be preferentially aligned perpendicular to the hosting large-scale filaments, while low-mass halos tend to have their spin axes aligned parallel to such structures. Extensive work has so far been conducted to understand this phenomenon under the assumption of cold dark matter and suggested that its origin should be closely related to the nonlinear evolution of the halo angular momentum in the anisotropic cosmic web. We present, for the first time, a numerical examination of this phenomenon assuming the presence of massive neutrinos, finding a clear and robust dependence of the threshold mass for the spin flip on the total neutrino mass. Our physical explanation is that the presence of more massive neutrinos retard the nonlinear evolution of the cosmic web, which in turn allows the halo spin vectors to better retain their memories of the initial tidal interactions in the nonlinear regime. Our finding implies that the statistical alignment of halo spins with the large-scale structures can be in principle used as a probe of the total neutrino mass.

Fisher for complements: extracting cosmology and neutrino mass from the counts-in-cells PDF

ABSTRACT We comprehensively analyse the cosmology dependence of counts-in-cells statistics. We focus on the shape of the one-point probability distribution function (PDF) of the matter density field at mildly non-linear scales. Based on large-deviation statistics, we parametrize the cosmology dependence of the matter PDF in terms of the linear power spectrum, the growth factor, the spherical collapse dynamics, and the non-linear variance. We extend our formalism to include massive neutrinos, finding that the total matter PDF is highly sensitive to the total neutrino mass Mν and can disentangle it from the clustering amplitude σ8. Using more than a million PDFs extracted from the Quijote simulations, we determine the response of the matter PDF to changing parameters in the νΛCDM model and successfully cross-validate the theoretical model and the simulation measurements. We present the first νΛCDM Fisher forecast for the matter PDF at multiple scales and redshifts, and its combination with the matter power spectrum. We establish that the matter PDF and the matter power spectrum are highly complementary at mildly non-linear scales. The matter PDF is particularly powerful for constraining the matter density Ωm, clustering amplitude σ8 and the total neutrino mass Mν. Adding the mildly non-linear matter PDF to the mildly non-linear matter power spectrum improves constraints on Ωm by a factor of 5 and σ8 by a factor of 2 when considering the three lowest redshifts. In our joint analysis of the matter PDF and matter power spectrum at three redshifts, the total neutrino mass is constrained to better than 0.01 eV with a total volume of 6 (Gpc h−1)3. We discuss how density-split statistics can be used to translate those encouraging results for the matter PDF into realistic observables in galaxy surveys.

Constraining the Neutrino Mass with the Drifting Coefficient of the Field Cluster Mass Function

Abstract A new diagnostics to break the degeneracy between the total neutrino mass (M ν ) and the primordial power spectrum amplitude (σ 8) using the drifting coefficient of the field cluster mass function is presented. Analyzing the data from the Cosmological Massive Neutrino Simulations, we first determine the numerical mass functions of the field clusters at various redshifts. Then, we compare the numerical results with the analytical model characterized by a single parameter called the drifting coefficient, which measures the drifts of the collapse density threshold, δ c, from the Einstein–de Sitter spherical value, δ sc, at a given mass scale. It is found that the analytic model for the field cluster mass function is found to work excellently even in the presence of massive neutrinos, and that its drifting coefficient evolves differently in the cosmologies with different values of M ν . At low redshifts (z ≲ 0.3) the more massive neutrinos drift δ c further from δ sc, while the opposite trend is found at higher redshifts (z ≳ 0.3). Speculating that this distinct redshift-dependent effect of massive neutrinos on the drifting coefficient of the field cluster mass function might help break the σ 8–M ν degeneracy, we also show that the sensitivity of this new diagnostics to M ν is high enough to discriminate the case of M ν = 0.1 eV from that of massless neutrinos.

Constraining sterile neutrino interpretations of the LSND and MiniBooNE anomalies with coherent neutrino scattering experiments

Results from the LSND and MiniBooNE experiments have been interpreted as evidence for a sterile neutrino with a mass near the electronvolt scale. Here we propose to test such a scenario by measuring the coherent elastic scattering rate of neutrinos from a pulsed spallation source. Coherent scattering is universal across all active neutrino flavors, and thus can provide a measurement of the total Standard Model neutrino flux. By performing measurements over different baselines and making use of timing information, it is possible to significantly reduce the systematic uncertainties and to independently measure the fluxes of neutrinos that originate as $\nu_{\mu}$ or as either $\nu_e$ or $\bar{\nu}_{\mu}$. We find that a 100 kg CsI detector would be sensitive to the large fraction of the sterile neutrino parameter space that could potentially account for the LSND and MiniBooNE anomalies.

All-inclusive interacting dark sector cosmologies

In this paper we explore possible extensions of Interacting Dark Energy cosmologies, where Dark Energy and Dark Matter interact non-gravitationally with one another. In particular, we focus on the neutrino sector, analyzing the effect of both neutrino masses and the effective number of neutrino species. We consider the Planck 2018 legacy release data combined with several other cosmological probes, finding no evidence for new physics in the dark radiation sector. The current neutrino constraints from cosmology should be therefore regarded as robust, as they are not strongly dependent on the dark sector physics, once all the available observations are combined. Namely, we find a total neutrino mass $M_\nu<0.15$ eV and a number of effective relativistic degrees of freedom of $N_{\rm eff}=3.03^{+0.33}_{-0.33}$, both at 95\% CL, which are close to those obtained within the $\Lambda$CDM cosmology, $M_\nu<0.12$ eV and $N_{\rm eff}=3.00^{+0.36}_{-0.35}$ for the same data combination.

On Stellar Evolution in a Neutrino Hertzsprung–Russell Diagram

Abstract We explore the evolution of a select grid of solar metallicity stellar models from their pre-main-sequence phase to near their final fates in a neutrino Hertzsprung–Russell diagram, where the neutrino luminosity replaces the traditional photon luminosity. Using a calibrated MESA solar model for the solar neutrino luminosity ( · × 1031 erg s−1) as a normalization, we identify ≃0.3 MeV electron neutrino emission from helium burning during the helium flash (peak , flux (10 pc/d)2 cm−2 s−1 for a star located at a distance of d parsec, timescale ≃3 days) and the thermal pulse (peak , flux × 107 (10 pc/d)2 cm−2 s−1, timescale ≃0.1 yr) phases of evolution in low-mass stars as potential probes for stellar neutrino astronomy. We also delineate the contribution of neutrinos from nuclear reactions and thermal processes to the total neutrino loss along the stellar tracks in a neutrino Hertzsprung–Russell diagram. We find, broadly but with exceptions, that neutrinos from nuclear reactions dominate whenever hydrogen and helium burn, and that neutrinos from thermal processes dominate otherwise.

Hunting the Glashow resonance with PeV neutrino telescopes

The Glashow resonant scattering, i.e. ν̄e + e− → W− → anything, offers us a possibility of disentangling ν̄e from the total astrophysical neutrino fluxes. Meanwhile, a great number of high-energy neutrino telescopes, with various detection mechanisms, are advancing towards a better understanding of one of the most energetic frontiers of the Universe. In this work, we investigate a connection between through-going muons at IceCube and the Glashow resonance signal through the channel W− → μ. We find that for IceCube, muons from ν̄e can induce a ∼20% excess of PeV events around the horizontal direction. However, the current statistic of IceCube is not enough to observe such an excess. We also address the novel possibility of ν̄e detection via W− → τ at telescopes aiming to detect Earth-skimming and mountain-penetrating neutrinos. The subsequent hadronic decay of a tau will induce an extensive air shower which can be detected by telescopes with Cherenkov or fluorescence techniques. Similar to IceCube, it is challenging to observe the Glashow resonance excess from the Earth-skimming neutrinos. Nevertheless, we find it is promising to observe Glashow resonance events with a mountain as the target.

IceCube Search for Galactic Neutrino Sources based on Very High Energy γ-ray Observations

Galactic cosmic rays reach energies of at least several PeV, and their interactions should generate γ-rays and neutrinos from decay of secondary pions. Therefore, Galactic sources have a guaranteed contribution to the total high-energy cosmic neutrino flux observed by IceCube. Assuming that the highest energy γ-rays are pionic, promising neutrino source candidates have been identified based on their spectra, and observing them is likely over the lifetime of the IceCube experiment. Here, we present the search for Galactic sources of high-energy cosmic neutrinos by focusing on sources identified by HAWC’s very high energy γ-ray survey.

Viable secret neutrino interactions with ultralight dark matter

Several anomalies in neutrino oscillation experiments point to the existence of a ∼1 eV sterile neutrino νs mixing with νe at the level of Ue4≅0.1, but such a neutrino is strongly disfavored by constraints on additional light degrees of freedom (δNeff) and total neutrino mass (∑νmν) from cosmology. “Secret neutrino interactions” that have been invoked to suppress the cosmological production of νs typically falter, but recently it was pointed out that νs could get a large mass in the early universe by coupling to ultralight dark matter ϕ, which can robustly suppress its production. The model has essentially two free parameters: mϕ, and ms,0, the mass of the sterile neutrino at early times, enhanced by its coupling to ϕ. I determine the parameter regions allowed by limits on δNeff and ∑νmν from the cosmic microwave background and big bang nucleosynthesis, using a simplified yet accurate treatment of neutrino oscillations in the early universe. This mechanism could have an important impact on laboratory experiments that suggest oscillations with sterile neutrinos.

Studying low−x structure function models with astrophysical tau neutrinos

Despite not have been yet identified by the IceCube detector, events generated from $\nu_{\tau}$ deep inelastic neutrino scattering in ice with varied topologies, such as double cascades (often called \textit{double bangs}), \textit{lollipops} and \textit{sugardaddies}, constitute a potential laboratory for low-x parton studies. Here we investigate these events, analyzing the effect of next-to-next-to-leading order (NNLO) Parton Distribution Function (PDFs) in the total neutrino--nucleon cross section, as compared with the color dipole formalism, where saturation effects play a major role. Energy deposit profiles in the `bangs' are also analysed in terms of virtual $W$-boson and tauon energy distributions and are found to be crucial in establishing a clear signal for gluon distribution determination at very small $x$. By taking the average (all flavor) neutrino flux ($\Phi_{\nu}\sim E_{\nu}^{-2.3}$) into differential cross sections as a function of $\tau$ and $W$ energies, we find significant deviations from pure DGLAP parton interactions for neutrino energies already at a few PeV. With these findings one aims at providing not only possible observables to be measured in large volume neutrino detectors in the near future, but also theoretical ways of unraveling QCD dynamics using unintegrated neutrino-nucleon cross sections in the ultrahigh-energy frontier.

Long-term temperature evolution of neutron stars undergoing episodic accretion outbursts

Context. Transient neutron star low-mass X-ray binaries undergo episodes of accretion, alternated with quiescent periods. During an accretion outburst, the neutron star heats up due to exothermic accretion-induced processes taking place in the crust. Besides the long-known deep crustal heating of nuclear origin, a likely non-nuclear source of heat, dubbed “shallow heating”, is present at lower densities. Most of the accretion-induced heat slowly diffuses into the core on a timescale of years. Over many outburst cycles, a state of equilibrium is reached when the core temperature is high enough that the heating and cooling (photon and neutrino emission) processes are in balance.Aims. We investigate how stellar characteristics and outburst properties affect the long-term temperature evolution of a transiently accreting neutron star. For the first time the effects of crustal properties are considered, particularly that of shallow heating.Methods. Using our codeNSCool, we tracked the thermal evolution of a neutron star undergoing outbursts over a period of 105yr. The outburst sequence is based on the regular outbursts observed from the neutron star transient Aql X-1. For each model we calculated the timescale over which equilibrium was reached and we present these timescales along with the temperature and luminosity parameters of the equilibrium state.Results. We performed several simulations with scaled outburst accretion rates, to vary the amount of heating over the outburst cycles. The results of these models show that the equilibrium core temperature follows a logarithmic decay function with the equilibrium timescale. Secondly, we find that shallow heating significantly contributes to the equilibrium state. Increasing its strength raises the equilibrium core temperature. We find that if deep crustal heating is replaced by shallow heating alone, the core would still heat up, reaching only a 2% lower equilibrium core temperature. Deep crustal heating may therefore not be vital to the heating of the core. Additionally, shallow heating can increase the quiescent luminosity to values higher than previously expected. The thermal conductivity in the envelope and crust, including the potentially low-conductivity pasta layer at the bottom of the crust, is unable to significantly alter the long-term internal temperature evolution. Stellar compactness and nucleon pairing in the core change the specific heat and the total neutrino emission rate as a function of temperature, with the consequences for the properties of the equilibrium state depending on the exact details of the assumed pairing models. The presence of direct Urca emission leads to the lowest equilibrium core temperature and the shortest equilibrium timescale.

First Measurement of Inclusive Muon Neutrino Charged Current Differential Cross Sections on Argon at E_{ν}∼0.8 GeV with the MicroBooNE Detector.

We report the first measurement of the double-differential and total muon neutrino charged current inclusive cross sections on argon at a mean neutrino energy of 0.8GeV. Data were collected using the MicroBooNE liquid argon time projection chamber located in the Fermilab Booster neutrino beam and correspond to 1.6×10^{20} protons on target of exposure. The measured differential cross sections are presented as a function of muon momentum, using multiple Coulomb scattering as a momentum measurement technique, and the muon angle with respect to the beam direction. We compare the measured cross sections to multiple neutrino event generators and find better agreement with those containing more complete treatment of quasielastic scattering processes at low Q^{2}. The total flux integrated cross section is measured to be 0.693±0.010(stat)±0.165(syst)×10^{-38} cm^{2}.