High-energy Neutrinos Research Articles

The spectra of IceCube Neutrino (SIN) candidate sources – IV. Spectral energy distributions and multiwavelength variability

ABSTRACT We present hybrid spectral energy distributions, combining photon, and neutrino fluxes, for a sample of blazars, which are candidate IceCube neutrino sources. We furthermore check for differences in our sources’ variability in the near-infrared, optical, X-ray, and γ-ray bands compared to a sample of non-neutrino source candidate blazars, and investigate the state of each blazar at the arrival time of high-energy neutrinos. We find no significant differences when comparing our sample with control sources, also in terms of their spectral energy distributions, and no correlation between flaring states and neutrino arrival times. Looking for signatures of hadronic production, we check for similar strengths of the γ-ray and neutrino fluxes and find a $2.2\, \sigma$ signal for our source candidates. The hybrid spectral energy distributions assembled here will form the basis of the next step of our project, namely lepto-hadronic modelling of these blazars to assess the physical likelihood of a neutrino connection.

Failed supernova simulations beyond black hole formation

ABSTRACT We present an axisymmetric failed supernova simulation beyond black hole formation, for the first time with numerical relativity and two-moment multi-energy neutrino transport. To ensure stable numerical evolution, we use an excision method for neutrino radiation hydrodynamics within the inner part of black hole domain. We demonstrate that our excision method is capable of stably evolving the radiation hydrodynamics in dynamical black hole space–time. As a remarkable signature of the final moment of proto-neutron star (PNS), we find the emergence of high-energy neutrinos. Those high-energy neutrinos are associated with the PNS shock surface being swallowed by the central black hole and could be a possible observable of failed supernovae.

Observation of seasonal variations of the flux of high-energy atmospheric neutrinos with IceCube

Atmospheric muon neutrinos are produced by meson decays in cosmic-ray-induced air showers. The flux depends on meteorological quantities such as the air temperature, which affects the density of air. Competition between decay and re-interaction of those mesons in the first particle production generations gives rise to a higher neutrino flux when the air density in the stratosphere is lower, corresponding to a higher temperature. A measurement of a temperature dependence of the atmospheric nu _{mu } flux provides a novel method for constraining hadronic interaction models of air showers. It is particularly sensitive to the production of kaons. Studying this temperature dependence for the first time requires a large sample of high-energy neutrinos as well as a detailed understanding of atmospheric properties. We report the significant (> 10 ; sigma ) observation of a correlation between the rate of more than 260,000 neutrinos, detected by IceCube between 2012 and 2018, and atmospheric temperatures of the stratosphere, measured by the Atmospheric Infrared Sounder (AIRS) instrument aboard NASA’s AQUA satellite. For the observed 10% seasonal change of effective atmospheric temperature we measure a 3.5(3)% change in the muon neutrino flux. This observed correlation deviates by about 2-3 standard deviations from the expected correlation of 4.3% as obtained from theoretical predictions under the assumption of various hadronic interaction models.

Low- and High-energy Neutrinos from SN 2023ixf in M101

Supernova (SN) 2023ixf in M101 is the closest SN explosion observed in the last decade. Therefore, it is a suitable test bed to study the role of jets in powering the SN ejecta. With this aim, we explored the idea that high-energy neutrinos could be produced during the interaction between the jets and the intense radiation field produced in the SN explosion and eventually be observed by the IceCube neutrino telescope. The lack of detection of such neutrinos has significantly constrained both the fraction of stellar collapses that produce jets and/or the theoretical models for neutrino production. Finally, we investigated the possibility of detecting low-energy neutrinos from SN 2023ixf with the Super- and Hyper-Kamiokande experiments, obtaining, in both cases, subthreshold estimates.

Flare Duty Cycle of Gamma-Ray Blazars and Implications for High-energy Neutrino Emission

Gamma-ray flares of blazars may be accompanied by high-energy neutrinos due to interactions of high-energy cosmic rays in the jet with photons, as suggested by the detection of the high-energy neutrino IceCube-170922A during a major gamma-ray flare from blazar TXS 0506+056 at the ∼3σ significance level. In this work, we present a statistical study of gamma-ray emission from blazars to constrain the contribution of gamma-ray flares to their neutrino output. We construct weekly binned light curves for 145 gamma-ray bright blazars in the Fermi Large Area Telescope Monitored Source List adding TXS 0506+056. We derive the fraction of time spent in the flaring state (flare duty cycle) and the fraction of energy released during each flare from the light curves with a Bayesian blocks algorithm. We find that blazars with lower flare duty cycles and energy fractions are more numerous among our sample. We identify a significant difference in flare duty cycles between blazar subclasses at a significance level of 5%. Then using a general scaling relation for the neutrino and gamma-ray luminosities, with a weighting exponent of γ = 1.0–2.0, normalized to the quiescent gamma-ray or X-ray flux of each blazar, we evaluate the neutrino energy flux of each gamma-ray flare. The gamma-ray flare distribution indicates that blazar neutrino emission may be dominated by flares for γ ≳ 1.5. The neutrino energy fluxes for 1 week and 10 yr bins are compared with the decl.-dependent IceCube sensitivity to constrain the standard neutrino emission models for gamma-ray flares. Finally, we present the upper-limit contribution of blazar gamma-ray flares to the isotropic diffuse neutrino flux.

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.

Flavor composition of neutrinos from choked gamma-ray bursts

The nature of the astrophysical sources responsible for producing the observed high energy neutrinos have yet to be unveiled. Choked gamma-ray bursts (CGRBs) are sources that have been proposed as being capable of generating the flux detected by IceCube, since no accompanying gamma-ray signal is expected from them, as required by observations. We focus on obtaining the neutrino flux and flavor composition corresponding to CGRBs under different assumptions for the target photon density and the magnetic field of the emission region. We consider the injection of both electrons and protons into the internal shocks of CGRBs, and using a steady-state transport equation, we account for all the relevant cooling processes. In particular, we include as a target for pγ interactions the usually adopted background of soft photons, which is a fraction of the thermalized emission originated at the shocked jet head. Additionally, we consider the synchrotron photons emitted by the electrons co-accelerated with the protons at the internal shocks in the jet. We also obtain the distribution of charged pions, kaons, and muons using the transport equation to account for the cooling effects due not only to synchrotron emission but also interactions with the soft photons in the ambient. We integrate the total diffuse flux of neutrinos of different flavors and compute the flavor ratios to be observed on Earth. As a consequence of the losses suffered mainly by pions and muons, we find these ratios to be dependent on the energy: for energies above ∼(105 − 106) GeV (depending on the magnetic field, proton-to-electron ratio, and jet power), we find that the electron flavor ratio decreases and the muon flavor ratio increases, while the tau flavor ratio increases only moderately. Our results are sensitive to the mentioned key physical parameters of the emitting region of CGRBs. Hence, the obtained flavor ratios are to be contrasted with cumulative data from ongoing and future neutrino instruments in order to assess the contribution of these sources to the diffuse flux of astrophysical neutrinos.

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.

Search for Correlations of High-energy Neutrinos Detected in IceCube with Radio-bright AGN and Gamma-Ray Emission from Blazars

Abstract The IceCube Neutrino Observatory sends realtime neutrino alerts with a high probability of being astrophysical in origin. We present a new method to correlate these events and possible candidate sources using 2089 blazars from the Fermi-LAT 4LAC-DR2 catalog and with 3413 active galactic nuclei (AGNs) from the Radio Fundamental Catalog. No statistically significant neutrino emission was found in any of the catalog searches. The result suggests that a small fraction, <1%, of the studied AGNs emit neutrinos that pass the alert criteria, and is compatible with prior evidence for neutrino emission presented by IceCube and other authors from sources such as TXS 0506 + 056 and PKS 1502 + 106. We also present cross-checks to other analyses that claim a significant correlation using similar data samples.

Probing neutrino production in high-energy astrophysical neutrino sources with the Glashow resonance

Probing neutrino production in high-energy astrophysical neutrino sources with the Glashow resonance

Present and future constraints on flavor-dependent long-range interactions of high-energy astrophysical neutrinos

The discovery of new, flavor-dependent neutrino interactions would provide compelling evidence of physics beyond the Standard Model. We focus on interactions generated by the anomaly-free, gauged, abelian lepton-number symmetries, specifically Le – Lμ, Le – Lτ , and Lμ – Lτ , that introduce a new matter potential sourced by electrons and neutrons, potentially impacting neutrino flavor oscillations. We revisit, revamp, and improve the constraints on these interactions that can be placed via the flavor composition of the diffuse flux of high-energy astrophysical neutrinos, with TeV–PeV energies, i.e., the proportion of νe, νμ, and ντ in the flux. Because we consider mediators of these new interactions to be ultra-light, lighter than 10−10 eV, the interaction range is ultra-long, from km to Gpc, allowing vast numbers of electrons and neutrons in celestial bodies and the cosmological matter distribution to contribute to this new potential. We leverage the present-day and future sensitivity of high-energy neutrino telescopes and of oscillation experiments to estimate the constraints that could be placed on the coupling strength of these interactions. We find that, already today, the IceCube neutrino telescope demonstrates potential to constrain flavor-dependent long-range interactions significantly better than existing constraints, motivating further analysis. We also estimate the improvement in the sensitivity due to the next-generation neutrino telescopes such as IceCube-Gen2, Baikal-GVD, KM3NeT, P-ONE, and TAMBO.

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.

Revise Thermal Winds of Remnant Neutron Stars in Gamma-Ray Bursts

It seems that the wealth of information revealed by the multi-messenger observations of the binary neutron star (NS) merger event, GW170817/GRB 170817A/kilonova AT2017gfo, places irreconcilable constraints to models of the prompt emission of this gamma-ray burst (GRB). The observed time delay between the merger of the two NSs and the trigger of the GRB and the thermal tail of the prompt emission can hardly be reproduced by these models simultaneously. We argue that the merger remnant should be an NS (last for, at least, a large fraction of 1 s), and that the difficulty can be alleviated by the delayed formation of the accretion disk due to the absorption of high-energy neutrinos emitted by the NS and the delayed emergence of effective viscosity in the disk. Further, we extend the consideration of the effect of the energy deposition of neutrinos emitted from the NS. If the NS is the central object of a GRB with a distance and duration similar to that of GRB 170817A, thermal emission of the thermal bubble inflated by the NS after the termination of accretion may be detectable. If our scenario is verified, it would be of interest to investigate the cooling of nascent NSs.

Two Candidate Obscured Tidal Disruption Events Coincident with High-energy Neutrinos

Recently, three optical tidal disruption event (TDE) candidates discovered by the Zwicky Transient Facility (ZTF) have been suggested to be coincident with high-energy neutrinos. They all exhibit unusually strong dust infrared echoes, with their peak times matching the neutrino arrival time even better than the optical peaks. We hereby report on two new TDE candidates that are spatially and temporally coincident with neutrinos by matching our sample of mid-infrared outbursts in nearby galaxies (MIRONG) with Gold alerts of IceCube high-energy neutrino events up to 2022 June. The two candidates show negligible optical variability according to their ZTF light curves and can therefore be classified as part of the growing population of obscured TDE candidates. The chance probability of finding two such candidates is about ∼3% by redistributing the MIRONG sources randomly in the Sloan Digital Sky Survey footprint, which will be as low as ∼0.1% (or ∼0.2%) if we limit to sources with increased fluxes (or variability amplitudes) comparable with the two matched sources. Our findings further support the potential connection between high-energy neutrinos and TDEs in dusty environments by increasing the total number of neutrino-associated TDE and TDE candidates to five, although the underlying physics remains poorly understood.

Review of the online analyses of multi-messenger alerts and electromagnetic transient events with the ANTARES neutrino telescope

By constantly monitoring a very large portion of the sky, neutrino telescopes are well-designed to detect neutrinos emitted by transient astrophysical events. Real-time searches with the ANTARES telescope have been performed to look for neutrino candidates coincident with gamma-ray bursts detected by the Swift and Fermi satellites, high-energy neutrino events registered by IceCube, transient events from blazars monitored by HAWC, photon-neutrino coincidences by AMON notices and gravitational wave candidates observed by LIGO/Virgo. By requiring temporal coincidence, this approach increases the sensitivity and the significance of a potential discovery. This paper summarises the results of the follow-up performed of the ANTARES telescope between January 2014 and February 2022, which corresponds to the end of the data-taking period.

High-Energy to Ultrahigh-Energy Neutrino Interactions

The cross sections for neutrino interactions with nucleons have been measured directly in accelerator experiments and through the zenith-angle and energy dependence of neutrino events at the IceCube Neutrino Observatory. Fluxes of high-energy neutrinos are produced at the Large Hadron Collider and by cosmic rays in the atmosphere. High-energy neutrinos also come from astrophysical and cosmic sources. The theory of neutrino interactions is reviewed. Current cross-section measurements and prospects for the future are discussed. The focus here is on neutrino interactions for energies larger than 1 TeV.

Optically informed searches of high-energy neutrinos from interaction-powered supernovae

ABSTRACT The interaction between the ejecta of supernovae (SNe) of Type IIn and a dense circumstellar medium can efficiently generate thermal ultraviolet/optical radiation and lead to the emission of neutrinos in the 1–103 TeV range. We investigate the connection between the neutrino signal detectable at the IceCube Neutrino Observatory and the electromagnetic signal observable by optical wide-field, high-cadence surveys to outline the best strategy for upcoming follow-up searches. We outline a semi-analytical model that connects the optical light-curve properties to the SN parameters and find that a large peak luminosity (${L_{\rm {peak}}\gtrsim 10^{43}{-}10^{44}\, \mathrm{erg \, s^{-1}}}$) and an average rise time (trise ≳ 10−40 d) are necessary for copious neutrino emission. Nevertheless, the most promising Lpeak and trise can be obtained for SN configurations that are not optimal for neutrino emission. Such ambiguous correspondence between the optical light-curve properties and the number of IceCube neutrino events implies that relying on optical observations only, a range of expected neutrino events should be considered (e.g. the expected number of neutrino events can vary up to two orders of magnitude for some among the brightest SNe IIn observed by the Zwicky Transient Facility up to now, SN 2020usa and SN 2020in). In addition, the peak in the high-energy neutrino curve should be expected a few trise after the peak in the optical light curve. Our findings highlight that it is crucial to infer the SN properties from multiwavelength observations rather than focusing on the optical band only to enhance upcoming neutrino searches.

IceCube and the origin of ANITA-IV events

Recently, the ANITA collaboration announced the detection of new, unsettling upgoing Ultra-High-Energy (UHE) events. Understanding their origin is pressing to ensure success of the incoming UHE neutrino program. In this work, we study their internal consistency and the implications of the lack of similar events in IceCube. We introduce a generic, simple parametrization to study the compatibility between these two observatories in Standard Model-like and Beyond Standard Model scenarios: an incoming flux of particles that interact with Earth nucleons with cross section σ, producing particle showers along with long-lived particles that decay with lifetime τ and generate a shower that explains ANITA observations. We find that the ANITA angular distribution imposes significant constraints, and when including null observations from IceCube only τ ~ 10−3–10−2 s and σ ~ 10−33– 10−32 cm2 can explain the data. This hypothesis is testable with future IceCube data. Finally, we discuss a specific model that can realize this scenario. Our analysis highlights the importance of simultaneous observations by high-energy optical neutrino telescopes and new UHE radio detectors to uncover cosmogenic neutrinos or discover new physics.

A Cross-correlation Study between IceCube Neutrino Events and the FERMI Unresolved Gamma-Ray Sky

With the coincident detections of electromagnetic radiation together with gravitational waves (GW170817) or neutrinos (TXS 0506+056), the new era of multimessenger astrophysics has begun. Of particular interest are the searches for correlation between the high-energy astrophysical neutrinos detected by the IceCube Observatory and gamma-ray photons detected by the Fermi Large Area Telescope (LAT). So far, only sources detected by the LAT have been considered in correlation with IceCube neutrinos, neglecting any emission from sources too faint to be resolved individually. Here we present the first cross-correlation analysis considering the unresolved gamma-ray background (UGRB) and IceCube events. We perform a thorough sensitivity study, and, given the lack of identified correlation, we place upper limits on the fraction of the observed neutrinos that would be produced in proton–proton or proton–γ interactions from the population of sources contributing to the UGRB emission and dominating its spatial anisotropy (aka blazars). Our analysis suggests that, under the assumption that there is no intrinsic cutoff and/or hardening of the spectrum above Fermi-LAT energies, and that all gamma rays from the unresolved blazars dominating the UGRB fluctuation field are produced by neutral pions from p–p (p–γ) interactions, up to 60% (30%) of such a population may contribute to the total neutrino events observed by IceCube. This translates into an O (1%) maximum contribution to the astrophysical high-energy neutrino flux observed by IceCube at 100 TeV.

High-energy Neutrino Production from AGN Disk Transients Impacted by the Circum-disk Medium

Various supernovae, compact object coalescences, and tidal disruption events are widely believed to occur embedded in active galactic nucleus (AGN) accretion disks and generate detectable electromagnetic signals. We collectively refer to them as AGN disk transients. The inelastic hadronuclear (pp) interactions between shock-accelerated cosmic rays and AGN disk materials shortly after the ejecta shock breaks out of the disk can produce high-energy neutrinos. However, the expected efficiency of neutrino production would decay rapidly by adopting a pure Gaussian density atmosphere profile applicable for stable gas-dominated disks. On the other hand, AGN outflows and disk winds are commonly found around AGN accretion disks. In this paper, we show that the circum-disk medium would further consume the shock kinetic energy to more efficiently produce high-energy neutrinos, especially for ∼ TeV−PeV neutrinos that IceCube detects. Thanks to the existence of the circum-disk medium, we find that the neutrino production will be enhanced significantly and make a much higher contribution to the diffuse neutrino background. Optimistically, ∼20% of the diffuse neutrino background can be contributed by AGN disk transients.