Cosmic Neutrinos Research Articles

Monte-carlo simulation of the effective lunar aperture for detection of ultra-high energy neutrinos with LOFAR

Ultra-high-energy (UHE) cosmic neutrinos interacting with the Moon’s regolith generate particle showers that emit Askaryan radiation. This radiation can be observed from the Earth using ground-based radio telescopes like LOFAR. We simulate the effective detection aperture for UHE neutrinos hitting the Moon. Under the same assumptions, results from this work are in good agreement with previous analytic parameterizations and Monte Carlo codes. The dependence of the effective detection aperture on the observing parameters, such as observing frequency and minimum detection threshold, and lunar characteristics like surface topography have been studied. Using a Monte Carlo simulation, we find that the detectable neutrino energy threshold is lowered when we include a realistic treatment of the inelasticity, transmission coefficient, and surface roughness. Lunar surface roughness at large scales enhances the total aperture for higher observation frequencies (ν≥1GHz\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\ u \\ge 1~\ extrm{GHz}$$\\end{document}) but has no significant effect on the LOFAR aperture. However, roughness at scales small compared to the wavelength reduces the aperture at all frequencies.

A Configurable 64-Channel ASIC for Cherenkov Radiation Detection from Space

This work presents the development of a 64-channel application-specific integrated circuit (ASIC), implemented to detect the optical Cherenkov light from sub-orbital and orbital altitudes. These kinds of signals are generated by ultra-high energy cosmic rays (UHECRs) and cosmic neutrinos (CNs). The purpose of this front-end electronics is to provide a readout unit for a matrix of silicon photo-multipliers (SiPMs) to identify extensive air showers (EASs). Each event can be stored into a configurable array of 256 cells where the on-board digitization can take place with a programmable 12-bits Wilkinson analog-to-digital converter (ADC). The sampling, the conversion process, and the main digital logic of the ASIC run at 200 MHz, while the readout is managed by dedicated serializers operating at 400 MHz in double data rate (DDR). The chip is designed in a commercial 65 nm CMOS technology, ensuring a high configurability by selecting the partition of the channels, the resolution in the interval 8–12 bits, and the source of its trigger. The production and testing of the ASIC is planned for the forthcoming months.

A Compact Particle Detector for Space-Based Applications: Development of a Low-Energy Module (LEM) for the NUSES Space Mission

NUSES is a planned space mission aiming to test new observational and technological approaches related to the study of relatively low-energy cosmic rays, gamma rays, and high-energy astrophysical neutrinos. Two scientific payloads will be hosted onboard the NUSES space mission: Terzina and Zirè. Terzina will be an optical telescope readout by SiPM arrays, for the detection and study of Cerenkov light emitted by Extensive Air Showers generated by high-energy cosmic rays and neutrinos in the atmosphere. Zirè will focus on the detection of protons and electrons up to a few hundred MeV and to 0.1–10 MeV photons and will include the Low Energy Module (LEM). The LEM will be a particle spectrometer devoted to the observation of fluxes of relatively low-energy electrons in the 0.1–7-MeV range and protons in the 3–50 MeV range along the Low Earth Orbit (LEO) followed by the hosting platform. The detection of Particle Bursts (PBs) in this Physics channel of interest could give new insight into the understanding of complex phenomena such as eventual correlations between seismic events or volcanic activity with the collective motion of particles in the plasma populating van Allen belts. With its compact sizes and limited acceptance, the LEM will allow the exploration of hostile environments such as the South Atlantic Anomaly (SAA) and the inner Van Allen Belt, in which the anticipated electron fluxes are on the order of 106 to 107 electrons per square centimeter per steradian per second. Concerning the vast literature of space-based particle spectrometers, the innovative aspect of the LEM resides in its compactness, within 10 × 10 × 10 cm3, and in its “active collimation” approach dealing with the problem of multiple scattering at these very relatively low energies. In this work, the geometry of the detector, its detection concept, its operation modes, and the hardware adopted will be presented. Some preliminary results from the Monte Carlo simulation (Geant4) will be shown.

Gravitational wave memory produced by cosmic background radiation

It is well known that energy fluxes will produce gravitational wave memory. The gravitational wave memory produced by background including cosmic microwave background (CMB), cosmic neutrino background (CνB), and gravitational wave background is investigated in this work. We construct a theory relating the gravitational wave memory strength to the energy flux of a stochastic background. We find that the resulted gravitational wave memory behaves as a constantly varying metric tensor. Such a varying metric tensor will introduce a quadrupole structure to the universe expansion. The gravitational wave memory due to the CMB is too small to be detected. But the gravitational wave memory due to the CνB and the gravitational wave background is marginally detectable. Interestingly, such detection can be used to estimate the neutrino masses and the properties of the gravitational wave background.

Influence of local structure on relic neutrino abundances and anisotropies

Gravitational potentials of the Milky Way and extragalactic structures can influence the propagation of the cosmic neutrino background (CNB). Of particular interest to future CNB observatories, such as PTOLEMY, is the CNB number density on Earth. In this study, we have developed a simulation framework that maps the trajectories of relic neutrinos as they move through the local gravitational environment. The potentials are based on the dark matter halos found in state-of-the-art cosmological N-body simulations, resulting in a more nuanced and realistic input than the previously employed analytical models. We find that the complex dark matter distributions, along with their dynamic evolution, influence the abundance and anisotropies of the CNB in ways unaccounted for by earlier analytical methods. Importantly, these cosmological simulations contain multiple instances of Milky Way-like halos that we employ to model a variety of gravitational landscapes. Consequently, we notice a variation in the CNB number densities that can be primarily attributed to the differences in the masses of these individual halos. For neutrino masses between 0.01 and 0.3 eV, we note clustering factors within the range of 1 + 𝒪(10-3) to 1 + 𝒪(1). Furthermore, the asymmetric nature of the underlying dark matter distributions within the halos results in not only overdense, but intriguingly, underdense regions within the full-sky anisotropy maps. Gravitational clustering appears to have a significant impact on the angular power spectra of these maps, leading to orders of magnitude more power on smaller scales beyond multipoles of ℓ = 3 when juxtaposed against predictions by primordial fluctuations. We discuss how our results reshape our understanding of relic neutrino clustering and how this might affect observability of future CNB observatories such as PTOLEMY.GitHub: our simulation code will be made visible here.

Simulation of RPC responses in mini-Iron Calorimeter for cosmic ray muon studies

A Prototype of the Iron Calorimeter for the India-based Neutrino Observatory is currently running in Madurai, India. This consists of twenty large area single gap Resistive Plate Chambers (RPCs) of size, ∼ 185 cm × 175 cm sandwiched between 11 layers of iron plates of thickness 5.6 cm. The total size of the mini-Iron Calorimeter(mini-ICAL) is 4 m × 4 m × 1.2 m and having a weight of 85 tons is a (1/600) scaled-down version of the ICAL detector of INO [1]. The detector is magnetized using iron plates and two sets of copper coils, each coil has 18 turns. The central region of the mini-ICAL, where the RPCs are placed, has a nearly uniform magnetic field of 1.4 T, similar to the magnetic field in the final ICAL detector. This prototype serves as the basis for the design of the demonstrators, which closely mimics the characteristics of a future ICAL. This detector is built to test the final electronics in the fringe field of the magnet and also to develop the experience to construct the ICAL detector. The mini-ICAL has been operational since 2018 and collects cosmic muon data with different configurations of RPC, and electronics, which are continuously changing for various R&D efforts. Dedicated efforts are made in parallel to measure the momentum and azimuthal angle of μ+ and μ- independently, which could be an important input to the cosmic neutrino event generators. To improve the precision of those measurements a dedicated simulation effort was made, particularly in the digitisation of the RPC signal for a real detector simulation, where efficiency, noise rate, strip multiplicity, etc. were matched with the real data collected during runs at mini-ICAL.

Neutrinos from GRB 221009A: producing ALPs and explaining LHAASO anomalous γ event

We propose a novel explanation for the 18 TeV gamma ray from GRB 221009A observed by LHAASO. High-energy neutrinos are converted into axion-like particles (ALPs) via their interaction with the cosmic neutrino background. Subsequently, ALPs are converted into high-energy photons in the magnetic field of our galaxy. We compute the fluxes of neutrinos, ALPs, and photons reaching Earth. IceCube's constraints on the neutrino flux from GRB 221009A translate into a severe upper bound on the photon flux. We find a range of parameters where all existing bounds are satisfied and the 18 TeV LHAASO photon can be explained. In the future, the specific correlation between the photon and neutrino flux reaching Earth from powerful neutrino sources with energies larger than 10 TeV such as GRBs or AGNs, can be used as a tool to differentiate our explanation from the alternatives suggested in the literature. We discuss how the interactions of our scenario can be embedded within electroweak gauge-invariant models, avoiding various cosmological and terrestrial bounds. We comment on the possibility of explaining the 251 TeV photon observed by the Carpet-2 detector, taking into account the bounds from the observation of high-energy neutrinos from TXS 0506+056.

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.

Best-case scenarios for neutrino capture experiments

A direct discovery of the cosmic neutrino background would bring to a closure the searches for relic left-over radiation predicted by the Hot Big Bang cosmology. Recently, the KATRIN experiment put a limit on the local relic neutrino overdensity with respect to the cosmological predicted average value at η ≲ 1011 [Phys. Rev. Lett. 129 (2022) 011806]. In this work, we first examine to what extent such values of η are conceivable. We show that even under cavalier assumptions, a cosmic origin of η ≳ 104 seems out of reach (with the caveat of forming bound objects under a new force,) but find that a hypothetical local source of low-energy neutrinos could achieve η ∼ 1011. Second, when such values are considered, we point out that the experimental signature in KATRIN and other neutrino-capture experiments changes, contrary to what has hitherto been assumed.Our results are model-independent and maximally accommodating as they only assume the Pauli exclusion principle.As intermittent physics target in the quest for CνB detection, we identify an experimental sensitivity to η ∼ 104 for which conceivable sources exist; to resolve the effect of a degenerate Fermi gas for such overdensity an energy resolution of 10 meV is required.

Bremsstrahlung from neutrino scattering via magnetic dipole moments

In this paper, we discuss bremsstrahlung induced by neutrino scattering. This process should exist since neutrinos are expected to couple to photons via magnetic dipole and transition moments. These moments are loop-induced and tiny in the Standard Model with neutrino masses but could be significantly enhanced in extended theories. As concrete example, we study the scattering of the two largest neutrino fluxes on earth, solar neutrinos and Cosmic Neutrino Background (CNB). It is tempting to consider this as a potential signature for CNB searches but it turns out that the signal is extremely small and unlikely to be observed.

High-energy Electromagnetic, Neutrino, and Cosmic-Ray Emission by Stellar-mass Black Holes in Disks of Active Galactic Nuclei

Some Seyfert galaxies are detected in high-energy gamma rays, but the mechanism and site of gamma-ray emission are unknown. Also, the origins of the cosmic high-energy neutrino and MeV gamma-ray backgrounds have been veiled in mystery since their discoveries. We propose emission from stellar-mass BHs (sBHs) embedded in disks of active galactic nuclei as their possible sources. These sBHs are predicted to launch jets due to the Blandford–Znajek mechanism, which can produce intense electromagnetic, neutrino, and cosmic-ray emissions. We investigate whether these emissions can be the sources of cosmic high-energy particles. We find that emission from internal shocks in the jets can explain gamma rays from nearby radio-quiet Seyfert galaxies including NGC 1068, if the Lorentz factor of the jets (Γj) is high. On the other hand, for moderate Γj, the emission can significantly contribute to the background gamma-ray and neutrino intensities in the ~MeV and ≲PeV bands, respectively. Furthermore, for moderate Γj with efficient amplification of the magnetic field and cosmic-ray acceleration, the neutrino emission from NGC 1068 and the ultrahigh-energy cosmic rays can be explained. These results suggest that the neutrino flux from NGC 1068 as well as the background intensities of MeV gamma rays, neutrinos, and the ultrahigh-energy cosmic rays can be explained by a unified model. Future MeV gamma-ray satellites will test our scenario for neutrino emission.

Saving CNB assisted EDE model at the expense of quantum corrections?

Abstract As it is suggested in [1,2], one can dynamically introduce the coincidence time-scale for EDE in the framework of a particular mass-varying-neutrino-model as a time at which neutrinos constituting the cosmic neutrino background enter the non-relativistic regime. The model does not predict, however, the right amount of EDE density because of smallness of neutrino masses. One may hope to adjust the parameters in such a way as to ensure that the two-loop contributions are kept small while at the same time the effective mass for scalar field that enters the expression of zero-point-energy (for the field trapped in the minimum of effective potential) is sufficient for explaining the needed amount of EDE. Unfortunately, the answer is not in the affirmative.

Cored dark matter halos in the cosmic neutrino background

We study the impact of the interaction between DM and the cosmic neutrino background on the evolution of galactic dark matter halos. The energy transfer from the cosmic background neutrinos to galactic dark matter can heat the center of the galaxy and make it cored. The core formation is efficient for the small galaxies such as the satellite galaxies of the Milky Way and we can put a conservative constraint on the non-relativistic elastic scattering cross section as σχν ≲ 10-31 cm2 for 0.1 keV dark matter.

The dark Stodolsky effect: constraining effective dark matter operators with spin-dependent interactions

We present a comprehensive discussion of the Stodolsky effect for dark matter (DM), and discuss two techniques to measure the effect and constrain the DM parameter space. The Stodolsky effect is the spin-dependent shift in the energy of a Standard Model (SM) fermion sitting in a bath of neutrinos. This effect, which scales linearly in the effective coupling, manifests as a small torque on the SM fermion spin and has historically been proposed as a method of detecting the cosmic neutrino background. We generalise this effect to DM, and give expressions for the induced energy shifts for DM candidates from spin-0 to spin-3/2, considering all effective operators up to mass dimension-6. In all cases, the effect scales inversely with the DM mass, but requires an asymmetric background. We show that a torsion balance experiment is sensitive to energy shifts of ΔE ≳ 10-28 eV, whilst a more intricate setup using a SQUID magnetometer is sensitive to shifts of ΔE ≳ 10-32 eV. Finally, we compute the energy shifts for a model of scalar DM, and demonstrate that the Stodolsky effect can be used to constrain regions of parameter space that are not presently excluded.

Constraints on Populations of Neutrino Sources from Searches in the Directions of IceCube Neutrino Alerts

Beginning in 2016, the IceCube Neutrino Observatory has sent out alerts in real time containing the information of high-energy (E ≳ 100 TeV) neutrino candidate events with moderate to high (≳30%) probability of astrophysical origin. In this work, we use a recent catalog of such alert events, which, in addition to events announced in real time, includes events that were identified retroactively and covers the time period of 2011–2020. We also search for additional, lower-energy neutrinos from the arrival directions of these IceCube alerts. We show how performing such an analysis can constrain the contribution of rare populations of cosmic neutrino sources to the diffuse astrophysical neutrino flux. After searching for neutrino emission coincident with these alert events on various timescales, we find no significant evidence of either minute-scale or day-scale transient neutrino emission or of steady neutrino emission in the direction of these alert events. This study also shows how numerous a population of neutrino sources has to be to account for the complete astrophysical neutrino flux. Assuming that sources have the same luminosity, an E −2.5 neutrino spectrum, and number densities that follow star formation rates, the population of sources has to be more numerous than 7 × 10−9 Mpc−3. This number changes to 3 × 10−7 Mpc−3 if number densities instead have no cosmic evolution.

Possible contribution of X-ray binary jets to the Galactic cosmic ray and neutrino flux

ABSTRACT For over a century, the identification of high-energy cosmic ray (CR) sources remains an open question. For Galactic CRs with energy up to 1015 eV, supernova remnants (SNRs) have traditionally been thought the main candidate source. However, recent TeV γ-ray observations have questioned the SNR paradigm. Propagating CRs are deflected by the Galactic magnetic field, hence, γ-rays and neutrinos produced via inelastic hadronic interactions are the only means for unveiling the CR sources. In this work, we study the γ-ray and neutrino emission produced by CRs accelerated inside Galactic jets of stellar-mass black holes in X-ray binaries (BHXBs). We calculate the intrinsic neutrino emission of two prototypical BHXBs , Cygnus X–1 and GX 339–4, for which we have high-quality, quasi-simultaneous multiwavelength spectra. Based on these prototypical sources, we discuss the likelihood of the 35 known Galactic BHXBs to be efficient CR accelerators. Moreover, we estimate the potential contribution to the CR spectrum of a viable population of BHXBs that reside in the Galactic plane. When these BHXBs go into outburst, they may accelerate particles up to hundreds of TeV that contribute to the diffuse γ-ray and neutrino spectra while propagating in the Galactic medium. Using HERMES, an open-source code that calculates the hadronic processes along the line of sight, we discuss the contribution of BHXBs to the diffuse γ-ray and neutrino fluxes, and compare these to their intrinsic γ-ray and neutrino emissions. Finally, we discuss the contribution of BHXBs to the observed spectrum of Galactic CRs.

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.

Calibration Facility for Detector Strings for the KM3NeT/ARCA Neutrino Telescope at the Capacity Laboratory

KM3NeT is a network of submarine Cherenkov neutrino telescopes under construction at two different sites in the Mediterranean Sea. ARCA, near Sicily in Italy, is optimised for the detection of cosmic neutrinos while ORCA, near Toulon in France, is optimized for atmospheric neutrinos. ARCA and ORCA are both arrays of thousands of optical sensors (Digital Optical Modules - DOMs), each consisting of 31 small photomultipliers (PMTs) housed inside a glass sphere, which detect the Cherenkov light produced by the secondary particles generated in the neutrino interactions. 18 DOMs are arranged on flexible strings, referred to as vertical Detection Units (DUs), anchored to the sea floor. Once completed, ARCA and ORCA will consist of 230 and 115 DUs, respectively. Each DOM of a string communicates at a dedicated wavelength to the shore station via a network of optical fibers transmitting optical and acoustic signal information as well as orientation information. Before the deployment, each DU is tested and calibrated in a dark room. The test bench is equipped with a full data acquisition system for communication, data processing, and time synchronisation. Several steps are needed to accomplish the DU calibration, including the High-Voltage (HV) tuning of the PMTs, checking of the acoustic receivers and calibration light sources in the DOM, and time calibration using laser signals distributed to all DOMs. Here, we describe the DU test and calibration facility at the CAPACITY laboratory (Campania AstroPArtiCle InfrastrucTure facilitY) in Caserta, focusing on the functional tests and calibrations performed at the end of the DU integration.

Light curves of BSM-induced neutrino echoes in the optically thin limit

High-energy neutrinos from astrophysical transients serve as a probe of neutrino physics beyond the Standard Model. In particular, nonstandard interaction of neutrinos with the cosmic neutrino background or dark matter (DM) may have imprints on not only their spectra but also the arrival and time-delay distributions. Assuming that the interaction occurs at most once during the neutrino propagation, we provide general analytic formulas for light curves of the neutrino echoes induced by BSM. The formulas can be used for constraining neutrino-neutrino scattering, neutrino-DM scattering, and other scattering processes experienced by relativistic particles.

Neutrino decoherence from generalised uncertainty

Quantum gravity models predict a minimal measurable length which gives rise to a modification in the uncertainty principle. One of the simplest manifestations of these generalised uncertainty principles is the linear quadratic generalised uncertainty principle which leads to a modified Heisenberg algebra. This can alter the usual von-Neumann evolution of density matrix to a Lindblad-type equation. We show how this can give rise to a decoherence in neutrino propagation in vacuum. The decoherence effects due to the linear quadratic generalised uncertainty principle are extremely minimal and is unlikely to be detectable in the existing or upcoming experimental facilities for any of the natural sources of neutrinos. We also show that, in principle, there can be other variants of generalised uncertainty principle which predicts verifiable decoherence effects for the cosmic neutrino background.