Detection Of High-energy Neutrinos Research Articles

EUSO-Offline: A comprehensive simulation and analysis framework

The complexity of modern cosmic ray observatories and the rich data sets they capture often require a sophisticated software framework to support the simulation of physical processes, detector response, as well as reconstruction and analysis of real and simulated data. Here we present the EUSO-Offline framework. The code base was originally developed by the Pierre Auger Collaboration, and portions of it have been adopted by other collaborations to suit their needs. We have extended this software to fulfill the requirements of Ultra-High Energy Cosmic Ray detectors and very high energy neutrino detectors developed for the Joint Exploratory Missions for an Extreme Universe Observatory (JEM-EUSO). These path-finder instruments constitute a program to chart the path to a future space-based mission like POEMMA. For completeness, we describe the overall structure of the framework developed by the Auger collaboration and continue with a description of the JEM-EUSO simulation and reconstruction capabilities. The framework is written predominantly in modern C++ (compliled against C++17) and incorporates third-party libraries chosen based on functionality and our best judgment regarding support and longevity. Modularity is a central notion in the framework design, a requirement for large collaborations in which many individuals contribute to a common code base and often want to compare different approaches to a given problem. For the same reason, the framework is designed to be highly configurable, which allows us to contend with a variety of JEM-EUSO missions and observation scenarios. We also discuss how we incorporate broad, industry-standard testing coverage which is necessary to ensure quality and maintainability of a relatively large code base, and the tools we employ to support a multitude of computing platforms and enable fast, reliable installation of external packages. Finally, we provide a few examples of simulation and reconstruction applications using EUSO-Offline.

Detecting High-energy Neutrino Minibursts from Local Supernovae with Multiple Neutrino Observatories

Growing evidence from multiwavelength observations of extragalactic supernovae (SNe) has established the presence of dense circumstellar material in Type II SNe. Interaction between the SN ejecta and the circumstellar material should lead to diffusive shock acceleration of cosmic rays and associated high-energy emission. Observation of high-energy neutrinos along with the MeV neutrinos from SNe will provide unprecedented opportunities to understand unanswered questions in cosmic-ray and neutrino physics. We show that current and future neutrino detectors can identify high-energy neutrinos from an extragalactic SN in the neighborhood of the Milky Way. We present the prospects for detecting high-energy neutrino minibursts from SNe in known local galaxies, and demonstrate how the network of multiple high-energy neutrino detectors will extend the horizon for the identification of high-energy SN neutrinos. We also discuss high-energy neutrino emission from SN 2023ixf.

Disentangling the Hadronic Components in NGC 1068

The recent detection of high-energy neutrinos by IceCube in the direction of the nearby Seyfert/starburst galaxy NGC 1068 implies that radio-quiet active galactic nuclei can accelerate cosmic-ray ions. Dedicated multimessenger analyses suggest that the interaction of these high-energy ions with ambient gas or photons happens in a region of the galaxy that is highly opaque for GeV–TeV gamma rays. Otherwise, the GeV–TeV emission would violate existing constraints provided by the Fermi Large Area Telescope (LAT) and the Major Atmospheric Gamma Imaging Cherenkov. The conditions of high optical depth are realized near the central supermassive black hole (SMBH). At the same time, the GeV emission detected by the Fermi LAT is likely related to the galaxy’s sustained star formation activity. In this work, we derive a 20 MeV–1 TeV spectrum of NGC 1068 using 14 yr of Fermi LAT observations. We find that the starburst hadronic component is responsible for NGC 1068's emission above ∼500 MeV. However, below this energy, an additional component is required. In the 20–500 MeV range, the Fermi LAT data are consistent with hadronic emission initiated by non-thermal ions interacting with gas or photons in the vicinity of the central SMBH. This highlights the importance of the MeV band to discover hidden cosmic-ray accelerators.

Hadronic versus Leptonic Origin of Gamma-Ray Emission from Supernova Remnants

GeV and TeV emission from the forward shocks of supernova remnants (SNRs) indicates that they are capable particle accelerators, making them promising sources of Galactic cosmic rays (CRs). However, it remains uncertain whether this γ-ray emission arises primarily from the decay of neutral pions produced by very-high-energy hadrons, or from inverse-Compton and/or bremsstrahlung emission from relativistic leptons. By applying a semi-analytic approach to non-linear diffusive shock acceleration, and calculating the particle and photon spectra produced in different environments, we parameterize the relative strength of hadronic and leptonic emission. We show that even if CR acceleration is likely to occur in all SNRs, the observed photon spectra may primarily reflect the environment surrounding the SNR: the emission is expected to look hadronic unless the ambient density is particularly low (with proton number density ≲0.1 cm−3) or the photon background is enhanced with respect to average Galactic values (with radiation energy density u rad ≳ 10 eV cm−3). We introduce a hadronicity parameter to characterize how hadronic or leptonic we expect a source to look based on its environment, which can be used to guide the interpretation of current γ-ray observations and the detection of high-energy neutrinos from SNRs.

High-energy Neutrino Emission Associated with Gravitational-wave Signals: Effects of Cocoon Photons and Constraints on Late-time Emission

We investigate prospects for the detection of high-energy neutrinos produced in the prolonged jets of short gamma-ray bursts (sGRBs). The X-ray light curves of sGRBs show extended emission components lasting for 100–1000 s, which are considered to be produced by prolonged engine activity. Jets produced by such activity should interact with photons in the cocoon formed by the propagation of the jet inside the ejecta of neutron star mergers. We calculate neutrino emission from jets produced by prolonged engine activity, taking account of the interaction between photons provided from the cocoon and cosmic rays accelerated in the jets. We find that IceCube-Gen2, a future neutrino telescope, with second-generation gravitational-wave detectors will probably be able to observe neutrino signals associated with gravitational waves with around 10 years of operation, regardless of the assumed value of the Lorentz factor of the jets. Neutrino observations may enable us to constrain the dissipation region of the jets. We apply this model to GRB 211211A, a peculiar long GRB whose origin may be a binary neutron star merger. Our model predicts that IceCube is unlikely to detect any associated neutrinos, but a few similar events will be able to put a meaningful constraint on the physical quantities of the prolonged engine activities.

First-principle calculation of birefringence effects for in-ice radio detection of neutrinos

The detection of high-energy neutrinos in the EeV range requires new detection techniques to cope with the small expected flux. The radio detection method, utilizing Askaryan emission, can be used to detect these neutrinos in polar ice. The propagation of the radio pulses has to be modeled carefully to reconstruct the energy, direction, and flavor of the neutrino from the detected radio flashes. Here, we study the effect of birefringence in ice, which splits up the radio pulse into two orthogonal polarization components with slightly different propagation speeds. This provides useful signatures to determine the neutrino energy and is potentially important to determine the neutrino direction to degree precision. We calculated the effect of birefringence from first principles where the only free parameter is the dielectric tensor as a function of position. Our code, for the first time, can propagate full RF waveforms, taking interference due to changing polarization eigenvectors during propagation into account. The model is available open-source through the NuRadioMC framework. We compare our results to in-situ calibration data from the ARA and ARIANNA experiments and find good agreement for the available time delay measurements. This indicates a significant improvement of the prediction power of birefringence effects compared to previous models. Finally, the implications and opportunities for neutrino detection are discussed.

Detection of ultra high energy cosmic rays and neutrinos with lunar orbital radio telescope

Particle cascades induced by ultra-high-energy (UHE) cosmic rays and neutrinos impacting on the lunar regolith usually radiate Cherenkov radio emissions due to the presence of excess negative charge, which is known as Askaryan effect. Several experiments have been carried out to detect the Cherenkov radio emissions in the lunar regolith. To prepare for future lunar Ultra-Long Wavelength (ULW, frequencies below 30 MHz) radio astronomy missions, we study the detection of the Cherenkov radio emissions with the ULW radio telescope that are operating at the lunar orbit. We have carried out instrument modelling and analytic calculations for the analysis of aperture, flux and event rate, and the analyses show the detectability of the Cherenkov radiation. Based on the properties of the Cherenkov radiation, we have demonstrated that the cosmic ray and neutrino events could be reconstructed with the three ULW vector antennas onboard the lunar satellites via measurements of the Askaryan radio pulse intensity, polarizations, etc. The results obtained by this study would be useful for future lunar radio explorer mission, where the detections of UHE cosmic rays and neutrinos could be successfully attempted.

Model Constraints Based on the IceCube Neutrino Nondetection of GRB 221009A

GRB 221009A is a bright gamma-ray burst (GRB) with isotropic energy larger than 1054 erg. Its fairly low redshift makes it a promising candidate for high-energy neutrino detection. However, a neutrino search for this GRB reported by the IceCube collaboration yielded a null result. In this paper, we utilize the upper limit from the IceCube observation to test different GRB prompt emission models. We find that, at least for this specific burst, the dissipative photosphere model could be ruled out in a large parameter space. The internal-shock model can survive only with a large bulk motion Lorentz factor Γ, where the most stringent and conservative constraints are Γ > ∼ 450 and Γ > ∼ 200, respectively. Also, the ratio of the total dissipated energy that goes into the protons and electrons (ϵ p /ϵ e ) can be constrained with a given Γ. For Γ < 400, ϵ p /ϵ e < 10 is required. For the Internal-collision-induced Magnetic Reconnection and Turbulence (ICMART) model, the constraint from GRB 221009A is modest. Under the ICMART model, only for extreme situations when most dissipated energy deposit into protons and all accelerated protons are suitable for producing neutrinos, a slightly large bulk motion (Γ > ∼ 250) is required.

Charmed Galaxies

The quest for finding the origins of cosmic rays has been going on for many decades. Cosmic rays as charged particles react to cosmic magnetic fields and typically travel in diffusive motion through the Universe. Their imprint on Earth therefore holds little information on their origin, and finding the sources of cosmic rays is a major challenge. The question of their origins has been a leading questions in physics and astrophysics in the past decades. To solve this riddle, a multimessenger approach is used, including cosmic-ray interaction products in the searches, specifically gamma-rays and neutrinos produced in the resulting particle showers. In this multimessenger picture, the detection of high-energy neutrinos from the cosmos by IceCube - as a unique tracer of cosmic rays - is an important piece of the puzzle. First evidence for neutrino emission from the active galaxies TXS0506+056 and NGC1068 indicates that a significant fraction comes from such sources. In this paper, the intriguing fact that gamma-rays seem to be absorbed in these sources discussed. The possibility of neutrinos being produced in regions of high photon or gas densities, together with the possibility of in the future revealing neutrinos from the decay of charmed particles, will be investigated, possibly opening a window to Charmed Galaxies.

Toise: a framework to describe the performance of high-energy neutrino detectors

Neutrinos offer a unique window to the distant, high-energy universe. Several next-generation instruments are being designed and proposed to characterize the flux of TeV–EeV neutrinos. The projected physics reach of the detectors is often quantified with simulation studies. However, a complete Monte Carlo estimate of detector performance is costly from a computational perspective, restricting the number of detector configurations considered when designing the instruments. In this paper, we present a new Python-based software framework, toise, which forecasts the performance of a high-energy neutrino detector using parameterizations of the detector performance, such as the effective areas, angular and energy resolutions, etc. The framework can be used to forecast performance of a variety of physics analyses, including sensitivities to diffuse fluxes of neutrinos and sensitivity to both transient and steady state point sources. This parameterized approach reduces the need for extensive simulation studies in order to estimate detector performance, and allows the user to study the influence of single performance metrics, like the angular resolution, in isolation. The framework is designed to allow for multiple detector components, each with different responses and exposure times, and supports paramterization of both optical- and radio-Cherenkov (Askaryan) neutrino telescopes. In the paper, we describe the mathematical concepts behind toise and introduce the reader to the use of the framework.

Propagation of Cosmic Rays in Plasmoids of AGN Jets-Implications for Multimessenger Predictions

After the successful detection of cosmic high-energy neutrinos, the field of multiwavelength photon studies of active galactic nuclei (AGN) is entering an exciting new phase. The first hint of a possible neutrino signal from the blazar TXS 0506+056 leads to the anticipation that AGN could soon be identified as point sources of high-energy neutrino radiation, representing another messenger signature besides the established photon signature. To understand the complex flaring behavior at multiwavelengths, a genuine theoretical understanding needs to be developed. These observations of the electromagnetic spectrum and neutrinos can only be interpreted fully when the charged, relativistic particles responsible for the different emissions are modeled properly. The description of the propagation of cosmic rays in a magnetized plasma is a complex question that can only be answered when analyzing the transport regimes of cosmic rays in a quantitative way. In this paper, therefore, a quantitative analysis of the propagation regimes of cosmic rays is presented in the approach that is most commonly used to model non-thermal emission signatures from blazars, i.e., the existence of a high-energy cosmic-ray population in a relativistic plasmoid traveling along the jet axis. It is shown that in the considered energy range of high-energy photon and neutrino emission, the transition between diffusive and ballistic propagation takes place, significantly influencing not only the spectral energy distribution, but also the lightcurve of blazar flares.

The supernova of the MAGIC gamma-ray burst GRB 190114C

We observed GRB 190114C (redshiftz = 0.4245), the first gamma-ray burst (GRB) ever detected at TeV energies, at optical and near-infrared wavelengths with several ground-based telescopes and theHubbleSpace Telescope, with the primary goal of studying its underlying supernova, SN 2019jrj. The monitoring spanned the time interval between 1.3 and 370 days after the burst, in the observer frame. We find that the afterglow emission can be modelled with a forward shock propagating in a uniform medium modified by time-variable extinction along the line of sight. A jet break could be present after 7 rest-frame days, and accordingly the maximum luminosity of the underlying supernova (SN) ranges between that of stripped-envelope core-collapse SNe of intermediate luminosity and that of the luminous GRB-associated SN 2013dx. The observed spectral absorption lines of SN 2019jrj are not as broad as in classical GRB SNe and are instead more similar to those of less-luminous core-collapse SNe. Taking the broad-lined stripped-envelope core-collapse SN 2004aw as an analogue, we tentatively derive the basic physical properties of SN 2019jrj. We discuss the possibility that a fraction of the TeV emission of this source might have had a hadronic origin and estimate the expected high-energy neutrino detection level with IceCube.

Modeling Particle Transport in Astrophysical Outflows and Simulations of Associated Emissions from Hadronic Microquasar Jets

In this work, after improving the formulation of the model on particle transport within astrophysical plasma outflows and constructing the appropriate algorithms, we test the reliability and effectiveness of our method through numerical simulations on well-studied galactic microquasars as the SS 433 and the Cyg X-1 systems. Then, we concentrate on predictions of the associated emissions, focusing on detectable high-energy neutrinos and γ -rays originated from the extragalactic M33 X-7 system, which is an X-ray binary discovered in 2006, located in the neighboring galaxy Messier 33, and has not yet been modeled in detail. The particle and radiation energy distributions, produced from magnetized hadronic astrophysical jets in the context of our method, are assumed to originate from decay and scattering processes taking place among the secondary particles created when hot (relativistic) protons of the jet scatter on thermal (cold) ones (p-p interaction mechanism inside the jet). These distributions are computed by solving the system of coupled integrodifferential transport equations of multiparticle processes (reactions chain) following the inelastic proton-proton (p-p) collisions. For the detection of such high-energy neutrinos as well as multiwavelength (radio, X-ray, and gamma-ray) emissions, extremely sensitive space telescopes and other γ -ray and neutrino detection instruments are in operation or have been designed like the CTA, IceCube, ANTARES, KM3NeT, and IceCube-Gen-2.

The spectra of IceCube neutrino (SIN) candidate sources – II. Source characterization

ABSTRACT Eight years after the first detection of high-energy astrophysical neutrinos by IceCube, we are still almost clueless as regards to their origin, although the case for blazars being neutrino sources is getting stronger. After the first significant association at the $3\!-\!3.5\, \sigma$ level in time and space with IceCube neutrinos, i.e. the blazar TXS 0506+056 at z = 0.3365, some of us have in fact selected a unique sample of 47 blazars, out of which ∼16 could be associated with individual neutrino track events detected by IceCube. Building upon our recent spectroscopy work on these objects, here we characterize them to determine their real nature and check if they are different from the rest of the blazar population. For the first time we also present a systematic study of the frequency of masquerading BL Lacs, i.e. flat-spectrum radio quasars with their broad lines swamped by non-thermal jet emission, in a γ-ray- and IceCube-selected sample, finding a fraction &amp;gt;24 per cent and possibly as high as 80 per cent. In terms of their broad-band properties, our sources appear to be indistinguishable from the rest of the blazar population. We also discuss two theoretical scenarios for neutrino emission, one in which neutrinos are produced in interactions of protons with jet photons and one in which the target photons are from the broad-line region. Both scenarios can equally account for the neutrino–blazar correlation observed by some of us. Future observations with neutrino telescopes and X-ray satellites will test them out.

Cosmic rays, neutrinos, and GeV-TeV gamma rays from starburst galaxy NGC 4945

The detection of high-energy astrophysical neutrinos and ultra-high-energy cosmic rays (UHECRs) provides a new way to explore sources of cosmic rays. One of the highest energy neutrino events detected by IceCube, tagged as IC35, is close to the UHECR anisotropy region detected by Pierre Auger Observatory. The nearby starburst (SB) galaxy, NGC 4945, is close to this anisotropic region and inside the mean angular error of the IC35 event. Considering the hypernovae contribution located in the SB region of NGC 4945, which can accelerate protons up to $\ensuremath{\sim}{10}^{17}\text{ }\text{ }\mathrm{eV}$ and inject them into the interstellar medium, we investigate the origin of this event around this starburst galaxy. We show that the interaction of these protons with the SB region's gas density could explain Fermi-LAT gamma-ray and radio observations if the magnetic field's strength in the SB region is the order of $\ensuremath{\sim}\mathrm{mG}$. Our estimated PeV neutrino events, in ten years, for this source is approximately 0.01 ($4\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}4}$) if a proton spectral index of 2.4 (2.7) is considered, which would demonstrate that IC35 is not produced in the central region of this SB galaxy. Additionally, we consider the superwind region of NGC 4945 and show that protons can hardly be accelerated in it up to UHECRs.

Soft gamma rays from low accreting supermassive black holes and connection to energetic neutrinos

The Universe is filled with a diffuse background of MeV gamma-rays and PeV neutrinos, whose origins are unknown. Here, we propose a scenario that can account for both backgrounds simultaneously. Low-luminosity active galactic nuclei have hot accretion flows where thermal electrons naturally emit soft gamma rays via Comptonization of their synchrotron photons. Protons there can be accelerated via turbulence or reconnection, producing high-energy neutrinos via hadronic interactions. We demonstrate that our model can reproduce the gamma-ray and neutrino data. Combined with a contribution by hot coronae in luminous active galactic nuclei, these accretion flows can explain the keV – MeV photon and TeV – PeV neutrino backgrounds. This scenario can account for the MeV background without non-thermal electrons, suggesting a higher transition energy from the thermal to nonthermal Universe than expected. Our model is consistent with X-ray data of nearby objects, and testable by future MeV gamma-ray and high-energy neutrino detectors.

Modeling of the tau and muon neutrino-induced optical Cherenkov signals from upward-moving extensive air showers

We present a detailed modeling and computation methodology to determine the optical Cherenkov signals produced by upward-moving extensive air showers (EASs) induced by $\ensuremath{\tau}$-leptons and muons, sourced from the interaction of high-energy astrophysical neutrinos interacting in the Earth. Following and extending the physics modeling and Cherenkov signal simulations performed in M. H. Reno et al. [Phys. Rev. D 100, 063010 (2019)], this scheme encompasses a new, state-of-the-art computation of the muon neutrino propagation inside the Earth and the contribution to the $\ensuremath{\tau}$-lepton muon decay channel. The modeling takes into account all possible $\ensuremath{\tau}$-lepton decay and muon energy loss channels that feed the optical Cherenkov emission, produced by both tau and muon initiated EASs. The EAS modeling uses the electron energy, angular, and lateral distributions in the EAS and their evolution as well as the wavelength dependence of the Cherenkov emission and its atmospheric attenuation. The results presented here are focused on the detection capabilities of suborbital (balloon-borne) and orbital (satellite) based instruments. The latter case was calculated for POEMMA [The Probe Of Extreme MultiMessenger Astrophysics] to compare to that presented in M. H. Reno et al. [Phys. Rev. D 100, 063010 (2019)], specifically including the muon-decay channel of $\ensuremath{\tau}$-leptons and the muonic EAS Cherenkov signal from muon neutrino interactions in the Earth. By detailing all these individual contributions to the optical Cherenkov emission and detection, we show how the ensemble that includes muonic channels provides a large detection capability for space-based, high-energy cosmic neutrino detection. Specifically, we show that for neutrino energies $\ensuremath{\lesssim}10\text{ }\text{ }\mathrm{PeV}$, the upward-EAS sensitivity due to muon neutrino interactions in the Earth begin to dominate over that for tau neutrino interactions, effectively extending the neutrino sensitivity to lower energies.

A Search for Ultra-high-energy Neutrinos from TXS 0506+056 Using the Pierre Auger Observatory

Abstract Results of a search for ultra-high-energy neutrinos with the Pierre Auger Observatory from the direction of the blazar TXS 0506+056 are presented. They were obtained as part of the follow-up that stemmed from the detection of high-energy neutrinos and gamma rays with IceCube, Fermi-LAT, MAGIC, and other detectors of electromagnetic radiation in several bands. The Pierre Auger Observatory is sensitive to neutrinos in the energy range from 100 PeV to 100 EeV and in the zenith-angle range from θ = 60° to θ = 95°, where the zenith angle is measured from the vertical direction. No neutrinos from the direction of TXS 0506+056 have been found. The results were analyzed in three periods: one of 6 months around the detection of IceCube-170922 A, coinciding with a flare period of TXS 0506+056, a second one of 110 days during which the IceCube collaboration found an excess of 13 neutrinos from a direction compatible with TXS 0506+056, and a third one from 2004 January 1 up to 2018 August 31, over which the Pierre Auger Observatory has been taking data. The sensitivity of the Observatory is addressed for different spectral indices by considering the fluxes that would induce a single expected event during the observation period. For indices compatible with those measured by the IceCube collaboration the expected number of neutrinos at the Observatory is well below one. Spectral indices as hard as 1.5 would have to apply in this energy range to expect a single event to have been detected.

Probing cosmic-ray accelerated light dark matter with IceCube

The direct detection of particle dark matter through its scattering with nucleons is of fundamental importance to understand the nature of DM . In this work, we propose that the high-energy neutrino detectors like IceCube can be used to uniquely probe the DM-nucleon cross-section for high-energy DM of ∼ PeV, up-scattered by the high-energy cosmic rays. We derive for the first time strong constraints on the DM-nucleon cross-section down to ∼ 10−32 cm2 at this energy scale for sub-GeV DM candidates. Such independent probe at energy scale far exceeding other existing direct detection experiments can therefore provide useful insights complementary to other searches.

Neutrino Counterparts of Fast Radio Bursts

Abstract The discovery of a luminous radio burst, FRB 200428, with properties similar to those of fast radio bursts (FRBs), in coincidence with an X-ray flare from the Galactic magnetar SGR 1935+2154, supports magnetar models for cosmological FRBs. The burst’s X-ray to radio fluence ratio, as well as the X-ray spectral shape and peak energy, are consistent with FRB 200428 being the result of an ultra-relativistic shock (powered, e.g., by an ejected plasmoid) propagating into a magnetized baryon-rich external medium; the shock simultaneously generates X-ray/gamma-rays via thermal synchrotron emission from electrons heated behind the shock, and coherent radio emission via the synchrotron maser mechanism. Here, we point out that a unique consequence of this baryon-loaded shock scenario is the generation of a coincident burst of high-energy neutrinos, generated by photohadronic interaction of relativistic ions—heated or accelerated at the shock—with thermal synchrotron photons. We estimate the properties of these neutrino burst FRB counterparts and find that a fraction ∼10−8–10−5 of the flare energy (or ∼10−4–10−1 of the radio isotropic energy) is channeled into production of neutrinos with typical energies ∼TeV–PeV. We conclude by discussing prospects for detecting this signal with IceCube and future high-energy neutrino detectors.