Astrophysical Neutrinos Research Articles

AM3: An Open-source Tool for Time-dependent Lepto-hadronic Modeling of Astrophysical Sources

We present the Astrophysical Multimessenger Modeling (AM 3 ) software. AM 3 is a documented open-source software (source code at https://gitlab.desy.de/am3/am3; user guide and documentation at https://am3.readthedocs.io/en/latest/) that efficiently solves the coupled integro-differential equations describing the temporal evolution of the spectral densities of particles interacting in astrophysical environments, including photons, electrons, positrons, protons, neutrons, pions, muons, and neutrinos. The software has been extensively used to simulate the multiwavelength and neutrino emission from active galactic nuclei (including blazars), gamma-ray bursts, and tidal disruption events. The simulations include all relevant nonthermal processes, namely synchrotron emission, inverse Compton scattering, photon–photon annihilation, proton–proton and proton–photon pion production, and photo-pair production. The software self-consistently calculates the full cascade of primary and secondary particles, including nonlinear feedback processes and predictions in the time domain. It also allows the user to track separately the particle densities produced by means of each distinct interaction process, including the different hadronic channels. With its efficient hybrid solver combining analytical and numerical techniques, AM 3 combines efficiency and accuracy at a user-adjustable level. We describe the technical details of the numerical framework and present three examples of applications to different astrophysical environments.

Association of the IceCube neutrinos with blazars in the CGRaBS sample

The origin of high-energy (HE) astrophysical neutrinos has remained an elusive hot topic in the field of HE astrophysics for the past decade. Apart from a handful of individual associations, the vast majority of HE neutrinos arise from unknown sources. While there are theoretically motivated candidate populations, such as blazars – a subclass of active galactic nuclei with jets pointed toward our line of sight – they have not been convincingly linked to HE neutrino production yet. Here, we perform a spatio-temporal association analysis between a sample of blazars (from the CGRaBS catalog) in the radio and optical bands and the most up-to-date IceCube HE neutrino catalog. We find that if the IceCube error regions are enlarged by 1° in quadrature, to account for unknown systematic errors at a maximal level, a spatio-temporal correlation between the multiwavelength light curves of the CGRaBS blazars and the IceCube HE neutrinos is hinted at, least at a 2.17σ significance level. On the other hand, when the IceCube error regions are taken as their published values, we do not find any significant correlations. A discrepancy in the blazar-neutrino correlation strengths, when using such minimal and enlarged error region scenarios, was also obtained in a recent study by the IceCube collaboration. In our study, this difference arises because several flaring blazars – coinciding with a neutrino arrival time – happen to narrowly miss the published 90%-likelihood error region of the nearest neutrino event. For all of the associations driving our most significant correlations, the flaring blazar is much less than 1° away from the published error regions. Therefore, our results indicate that the question of the blazar-neutrino connection is highly sensitive to the reconstruction of the neutrino error regions, whose reliability is expected to improve with the next generation of neutrino observatories.

Towards multi-messenger observations of core-collapse supernovae harbouring choked jets

Context. Over the past decade, choked jets have attracted particular attention as potential sources of high-energy cosmic neutrinos. It is challenging to test this hypothesis because of the missing gamma-ray counterpart. An identification of other electromagnetic signatures is therefore crucial. Extended H envelopes surrounding collapsing massive stars might choke launched jets. In addition, the same progenitors are expected to produce a shock-breakout signal in the ultraviolet (UV) and optical that lasts several days. Early UV radiation in particular carries important information about the presence and nature of choked jets. Aims. While UV observations of core-collapse supernovae have so far been limited, the full potential of observations in this spectral band will soon be transformed by the ULTRASAT satellite mission with its unprecedented field of view. We investigated the detection prospects of choked jet progenitors by ULTRASAT in relation to their visibility in the optical band by the currently operating telescope ZTF. In addition, as choked jets can produce neutrinos via hadronic and photohadronic interactions in choked jets, we also investigated how neutrino observations by existing Cherenkov high-energy neutrino telescopes (e.g. IceCube and KM3NeT) can be used in association with electromagnetic signals from shock-breakout events. Methods. By considering fiducial parameters of the source population and instrument performances, we estimated the maximum redshift up to which ULTRASAT and ZTF are able to detect ultraviolet and optical signals from these explosions, respectively. Furthermore, we discuss coordinated multi-messenger observations using ULTRASAT, ZTF, and high-energy neutrino telescopes. Results. We find that ULTRASAT will double the volume of the sky that is currently visible by ZTF for the same emitting sources. This will enlarge the sample of observed Type II supernovae by ∼60%. For optimised multi-messenger detections, the delay between neutrinos produced at the shock breakout (during the jet propagation inside the stellar envelope) and ULTRASAT observations should be of ∼4 (5) days, with subsequent follow-up by instruments such as ZTF about one week later. We estimate that fewer than 1% of the core-collapse supernovae from red supergiant stars that are detectable in UV with ULTRASAT might host a choked jet and release TeV neutrinos. Electromagnetic and neutrino detections, if accompanied by additional photometric and spectroscopic follow-up with compelling evidence for a relativistic jet launched by the central engine of the source, would suggest that core-collapse supernovae harbouring choked jets are the main contributors to the diffuse astrophysical high-energy neutrino flux.

Track signals at IceCube from subleading channels

Tracks events at the IceCube Observatory are characterized by an energetic muon crossing several kilometers before decaying. Such muons are dominantly produced in charged current (CC) muon neutrino-hadron interactions. However, muons are also produced through W-boson production and in the decay of tau leptons and heavy mesons created in neutral and charged current interactions induced by all neutrino flavors. In this paper, we investigate the contribution of these subleading channels to events characterized as tracks at the IceCube. Our results indicate that these channels correspond to a non-negligible fraction of the high-energy starting events track events. In addition, we show that its contributions are concentrated in muons that are less energetic than those arising from muonic neutrino CC interactions for the same visible energies of the process. Finally, we investigate the impact of these additional channels on the description of the astrophysical neutrino flux, and we find that the inclusion of these subleading processes are important in determining the parameters of the astrophysical neutrino flux. Published by the American Physical Society 2024

Prospects for measuring the time variation of astrophysical neutrino sources at dark matter detectors

We study the prospects for measuring the time variation of solar and atmospheric neutrino fluxes at future large-scale xenon and argon dark matter detectors. For solar neutrinos, a yearly time variation arises from the eccentricity of Earth’s orbit and, for charged current interactions, from a smaller energy-dependent day-night variation due to flavor regeneration as neutrinos travel through Earth. For a 100-ton xenon detector running for ten years with a xenon-136 fraction of ≲0.1%, in the electron recoil channel a time-variation amplitude of about 0.8% is detectable with a power of 90% and the level of significance of 10%. This is sufficient to detect time variation due to eccentricity, which has amplitude of ∼3%. In the nuclear recoil channel, the detectable amplitude is about 10% under current detector resolution and efficiency conditions, and this generally reduces to about 1% for improved detector resolution and efficiency, the latter of which is sufficient to detect time variation due to eccentricity. Our analysis assumes both known and unknown periods. We provide scalings to determine the sensitivity to an arbitrary time-varying amplitude as a function of detector parameters. Identifying the time variation of the neutrino fluxes will be important for distinguishing neutrinos from dark matter signals and other detector-related backgrounds and extracting properties of neutrinos that can be uniquely studied in dark matter experiments. Published by the American Physical Society 2024

Refractive index in the JUNO liquid scintillator

In the field of rare event physics, it is common to have huge masses of organic liquid scintillator as detection medium. In particular, they are widely used to study neutrino properties or astrophysical neutrinos. Thanks to its safety properties (such as low toxicity and high flash point) and easy scalability, linear alkyl benzene is the most common solvent used to produce liquid scintillators for large mass experiments. The knowledge of the refractive index is a pivotal point to understand the detector response, as this quantity (and its wavelength dependence) affects the Cherenkov radiation and photon propagation in the medium. In this paper, we report the measurement of the refractive index of the JUNO liquid scintillator between 260–1064 nm performed with two different methods to cover wide range of refractive index (an ellipsometer and a refractometer), with a sub percent level precision. In addition, we used an interferometer to measure the group velocity in the JUNO liquid scintillator and verify the expected value derived from the refractive index measurements.

Solar flare observations with the Radio Neutrino Observatory Greenland (RNO-G)

The Radio Neutrino Observatory – Greenland (RNO-G) seeks discovery of ultra-high energy neutrinos from the cosmos through their interactions in ice. The science program extends beyond particle astrophysics to include radioglaciology and, as we show herein, solar observations, as well. Currently seven of 35 planned radio-receiver stations (24 antennas/station) are operational. These stations are sensitive to impulsive radio signals with frequencies between 80 and 700 MHz and feature a neutrino trigger threshold for recording data close to the thermal floor. RNO-G can also trigger on elevated signals from the Sun, resulting in nanosecond resolution time-domain flare data; such temporal resolution is significantly shorter than from most dedicated solar observatories. In addition to possible RNO-G solar flare polarization measurements, the Sun also represents an extremely useful above-surface calibration source.Using RNO-G data recorded during the summers of 2022 and 2023, we find signal excesses during solar flares reported by the solar-observing Callisto network and also in coincidence with ∼2/3 of the brightest excesses recorded by the SWAVES satellite. These observed flares are characterized by significant time-domain impulsivity. Using the known position of the Sun, the flare sample is used to calibrate the RNO-G absolute pointing on the radio signal arrival direction to sub-degree resolution. We thus establish the Sun as a regularly observed astronomical calibration source to provide the accurate absolute pointing required for neutrino astronomy.

New limits on neutrino decay from high-energy astrophysical neutrinos

New limits on neutrino decay from high-energy astrophysical neutrinos

Search for neutrino emission from GRB 221009A using the KM3NeT ARCA and ORCA detectors

Gamma-ray bursts are promising candidate sources of high-energy astrophysical neutrinos. The recent GRB 221009A event, identified as the brightest gamma-ray burst ever detected, provides a unique opportunity to investigate hadronic emissions involving neutrinos. The KM3NeT undersea neutrino detectors participated in the worldwide follow-up effort triggered by the event, searching for neutrino events. In this paper, we summarize subsequent searches, in a wide energy range from MeV up to a few PeVs. No neutrino events are found in any of the searches performed. Upper limits on the neutrino emission associated with GRB 221009A are computed.

Challenges and solutions in the mechanical design of the KM3NeT ARCA Base Module

The goal of KM3NeT ARCA is to search for astrophysical neutrinos with a neutrino telescope operating on the bottom of the Mediterranean Sea, at 3500 m depth, for more than 15 years. A key component of the neutrino telescope is the Base Module (BM) which is built with a titanium cylindrical container hosting optics, electronics and power boards. From a mechanical point of view the BM must be water tight and resistant to an external water pressure up to 350 bar, be able to stand all mechanical stresses along its life, provide excellent corrosion resistance for 15 years at least in sea water and be able to host all internal components safely, by guaranteeing an adequate heat dissipation. The tests performed for the qualification of the mechanical design are presented here. The first BM is successfully operating underwater since 2015.

Potential of water-Cherenkov air shower arrays for detecting transient sources of high-energy astrophysical neutrinos

Potential of water-Cherenkov air shower arrays for detecting transient sources of high-energy astrophysical neutrinos

Constraints on quantum spacetime-induced decoherence from neutrino oscillations

We investigate the implications of decoherence induced by quantum spacetime properties on neutrino oscillation phenomena. We develop a general formalism where the evolution of neutrinos is governed by a Lindblad-type equation and we compute the oscillation damping factor for various models that have been proposed in the literature. Furthermore, we discuss the sensitivity to these effects of different types of neutrino oscillation experiments, encompassing astrophysical, atmospheric, solar, and reactor neutrino experiments. By using neutrino oscillation data from long-baseline reactors and atmospheric neutrino observations, we establish stringent constraints on the energy scale governing the strength of the decoherence induced by stochastic metric fluctuations, amounting to, respectively, EQG≥2.6×1034 GeV and EQG≥2.5×1055 GeV. Published by the American Physical Society 2024

Characterization of the astrophysical diffuse neutrino flux using starting track events in IceCube

Characterization of the astrophysical diffuse neutrino flux using starting track events in IceCube

AT2021lwx: Another Neutrino-coincident Tidal Disruption Event with a Strong Dust Echo?

We discuss the possible association of an astrophysical neutrino (IC220405B) with the recently reported, extremely energetic tidal disruption event (TDE) candidate AT2021lwx (ZTF20abrbeie, aka “Scary Barbie”) at redshift z = 0.995. Although the TDE is about 2.°6 off the direction of the reconstructed neutrino event (outside the 90% confidence level localization region), the TDE candidate shares some important characteristics with so-far-reported neutrino–TDE associations: a strong infrared dust echo, high bolometric luminosity, a neutrino time delay with respect to the peak mass accretion rate of the order of a hundred days, and a high observed X-ray luminosity. We interpret this new association using an isotropic emission model, where neutrinos are produced by the collision of accelerated protons with infrared photons. After accounting for the high redshift of AT2021lwx (by interpreting the data in the supermassive black hole (SMBH) frame), we find that the expected neutrino fluences and neutrino time delays are qualitatively comparable to the other TDEs. Since data are only available up to 300 days postpeak in the SMBH frame, significant uncertainties exist in the dust echo interpretation, and therefore in the predicted number of neutrinos detected, Nν≃3.0×10−3−0.012 . We recommend further follow-up of this object for an extended period and suggest refining the reconstruction of the neutrino arrival direction in this particular case.

Testing secret interaction with astrophysical neutrino point sources

Recently, the IceCube collaboration observed a neutrino excess in the direction of NGC 1068 with high statistical significance.This constitutes the second detection of an astrophysical neutrino point source after the discovery of a variable emission originating from the blazar TXS 0506+056. Neutrinos emitted by these sources traverse huge, well-determined distances on their way to Earth. This makes them a promising tool to test new physics in the neutrino sector. We consider secret interactions with the cosmic neutrino background and discuss their impact on the flux of neutrino point sources. The observation of emission from NGC 1068 and TXS 0506+056 can then be used to putlimits on the strength of the interaction. We find that our ignorance of the absolute neutrino masses has a strong impact and, therefore, we present limits in two benchmark scenarios with the sum of the neutrino masses around their lower and upper limits.

A hybrid approach to event reconstruction for atmospheric Cherenkov Telescopes combining machine learning and likelihood fitting

The imaging atmospheric Cherenkov technique provides potentially the highest angular resolution achievable in astronomy at energies above the X-ray waveband. High-resolution measurements provide the key to progress on many of the major questions in high energy astrophysics, including the sites and mechanisms of particle acceleration to PeV energies. The huge potential of the next-generation CTA observatory in this regard can be realised with the help of improved algorithms for the reconstruction of the air-shower direction and energy.Hybrid methods combining maximum-likelihood-fitting techniques with neural networks represent a particularly promising approach and have recently been successfully applied for the reconstruction of astrophysical neutrinos. Here, we present the FreePACT algorithm, a hybrid reconstruction method for IACTs. In this, making use of the neural ratio estimation technique from the field of likelihood-free inference, the analytical likelihood used in traditional image likelihood fitting is replaced by a neural network that approximates the charge probability density function for each pixel in the camera.The performance of this improved algorithm is demonstrated using simulations of the planned CTA southern array. For this setupFreePACT provides significant performance improvements over analytical likelihood techniques, with improvements in angular and energy resolution of 25% or more over a wide energy range and an angular resolution as low as 40′′ at energies above 50TeV for observations at 20° zenith angle. It also yields more accurate estimations of the uncertainties on the reconstructed parameters and significantly speeds up the reconstruction compared to analytical likelihood techniques while showing the same stability with respect to changes in the observation conditions. Therefore, the FreePACT method is a promising upgrade over the current state-of-the-art likelihood event reconstruction techniques.

Generating Baikal-GVD high energy cascade-like neutrino events with a GEANT4-based simulation toolkit

Abstract Using GEANT4 and Cosmic Ray Monte Carlo (CRMC) software packages, we developed a new simulation toolkit for astrophysical neutrino telescopes. By configuring the Baikal-GVD detector and comparing the vertex position and direction of incident particles, as well as the channel-by-channel signals, to the events detected by Baikal-GVD, we successfully generated 13 high-energy cascades neutrino events with the toolkit. Our analysis revealed a systematic offset between the reconstructed shower position and the true interaction position, with a distance close to the scale of the shower maximum of -0.64±1.32 m. The results of the linear fit between real data and simulation demonstrate reliable outcomes, indicating that the simulation toolkit is capable of studying the performance of astrophysical neutrino telescopes.

High-energy neutrino deep inelastic scattering cross sections

We present a state-of-the-art prediction for cross sections of neutrino deep inelastic scattering (DIS) from nucleon at high neutrino energies, Eν, up to 1000 EeV (1012 GeV). Our calculations are based on the latest CT18 NNLO parton distribution functions (PDFs) and their associated uncertainties. To make predictions for the highest energies, we extrapolate the PDFs to small x according to several procedures and assumptions, thus affecting the uncertainties at ultrahigh Eν; we quantify the uncertainties corresponding to these choices. Similarly, we quantify the uncertainties introduced by the nuclear corrections that are required to evaluate neutrino-nuclear cross sections for the neutrino observatories. These results can be applied to currently running astrophysical neutrino observatories, such as IceCube and KM3NeT, as well as various future experiments that have been proposed. Published by the American Physical Society 2024

A search for AGN sources of the IceCube diffuse neutrino flux

The origin of the diffuse astrophysical neutrino flux measured by the IceCube Observatory remains largely unknown.Although NGC 1068 and TXS 0506+056 have been identified as potential neutrino sources, the diffuse flux of neutrinos must have additional sources that have not yet been identified.Here we investigate potential correlations between IceCube's neutrino events and the Fermi and MOJAVE source catalogs, using the publicly-available IceCube data set. We perform three separate spatially-dependent, energy-dependent, and time-dependent searches, and find no statistically significant sources outside of NGC 1068. We find that, under the most optimistic assumptions of a spectral index of 2.0 and a neutrino flux uncorrelated with the gamma ray flux, no more than 13% of IceCube's neutrino flux originates from blazars over the whole sky. Then, using an energy-dependent likelihood analysis, the limit on neutrinos originating from blazars reduces to 9% in the Northern hemisphere under the same spectral index and flux assumptions. Finally, we set limits on individual sources from the MOJAVE radio catalog after finding no statistically significant time-flaring sources.

Collective Neutrino Oscillations and Heavy-element Nucleosynthesis in Supernovae: Exploring Potential Effects of Many-body Neutrino Correlations

In high-energy astrophysical processes involving compact objects, such as core-collapse supernovae or binary neutron star mergers, neutrinos play an important role in the synthesis of nuclides. Neutrinos in these environments can experience collective flavor oscillations driven by neutrino–neutrino interactions, including coherent forward scattering and incoherent (collisional) effects. Recently, there has been interest in exploring potential novel behaviors in collective oscillations of neutrinos by going beyond the one-particle effective or “mean-field” treatments. Here, we seek to explore implications of collective neutrino oscillations, in the mean-field treatment and beyond, for the nucleosynthesis yields in supernova environments with different astrophysical conditions and neutrino inputs. We find that collective oscillations can impact the operation of the ν p-process and r-process nucleosynthesis in supernovae. The potential impact is particularly strong in high-entropy, proton-rich conditions, where we find that neutrino interactions can nudge an initial ν p-process neutron-rich, resulting in a unique combination of proton-rich low-mass nuclei as well as neutron-rich high-mass nuclei. We describe this neutrino-induced neutron-capture process as the “ν i-process.” In addition, nontrivial quantum correlations among neutrinos, if present significantly, could lead to different nuclide yields compared to the corresponding mean-field oscillation treatments, by virtue of modifying the evolution of the relevant one-body neutrino observables.