Diffuse Flux Of Astrophysical Neutrinos Research Articles

High-Energy and Ultra-High-Energy Neutrino Astrophysics

The origin of high-energy cosmic rays, and their behavior in astrophysical sources, remains an open question. Recently, new ways to address this question have been made possible by the observation of a new astrophysical messenger, namely neutrinos. The IceCube telescope has detected a diffuse flux of astrophysical neutrinos in the TeV-PeV energy range, likely produced in astrophysical sources accelerating cosmic rays, and more recently it has reported on a few candidate individual neutrino sources. Future experiments will be able to improve on these measurements quantitatively, by the detection of more events, and qualitatively, by extending the measurement into the EeV energy range. In this paper, we review the main features of the neutrino emission and sources observed by IceCube, as well as the main candidate sources that could contribute to the diffuse neutrino flux. As a parallel question, we review the status of high-energy neutrinos as a probe of Beyond the Standard Model physics coupling to the neutrino sector.

Neutrino production in starburst galaxies

ABSTRACT Understanding the origin of the diffuse flux of high-energy astrophysical neutrinos detected by IceCube has become a challenging issue within present High Energy Astrophysics. In this work, we present a model to explore the potential neutrino emission of starburst galaxies (SBG) by considering three different neutrino production zones that can be associated to a typical single SBG. The first zone is the starburst nucleus, where due to the high rate of supernova explosions, a significant amount of protons can be accelerated to high energies and undergo pp interactions with cold protons of the interstellar medium. The second zone we consider is the corresponding to the starburst wind, which is formed by the hot gas that emerges from the nucleus and interacts with the intergalactic medium generating shocks. Protons accelerated there can undergo pp interactions with the ambient matter. The third neutrino production zone we consider, is an external one, where we account for the possibility that protons escaping from the whole system interact with the cosmic microwave background. Finally, adding the neutrino contributions of the three zones, we calculate the diffuse neutrino flux and the diffuse photon flux by integration on the redshift range appropriate for SBG. We find that the model behaves well applied to nearby galaxies such as M82 and NGC 253. The contributions made to the diffuse neutrino flux are able to explain part of the data provided by IceCube if typical parameters are considered.

Leptohadronic multi-messenger modeling of 324 gamma-ray blazars

Context. The origin of the diffuse astrophysical neutrino flux observed by the IceCube experiment is still under debate. Multiple associations have been reported between high-energy neutrino events and individual bla/ars, such as the source TXS 0506+056, which are active galaxies with relativistic jets pointing toward Earth. From a theoretical perspective, the properties of these sources as neutrino emitters are not yet well understood. Aims. By systematically modeling the effect of cosmic-ray protons on the multiwavelength data from the largest sample of bright gamma-ray bla/ars to date, we expect to learn about the multi-messenger nature of the active galaxy population as a whole, as well as the relationship between neutrino production and the multiwavelength spectrum of these sources. Methods. We predict the emitted multiwavelength and neutrino spectrum using a self-consistent numerical radiation model applied individually to each source in the sample. We then study the properties of the full population and identify empirical relations. We focus on public multiwavelength data from the radio to the gamma-ray bands from a sample of 324 bla/ars detected by the Fermi Large Area Telescope (LAT), most of which are flat-spectrum radio quasars (FSRQs). This amounts to 34% of all FSRQs in the latest Fermi catalog. Results. We demonstrate that the optical and gigaelectronvolt gamma-ray broadband features are generally well described by electron emission, which helps for the location of the emission region relative to the central black hole to be constrained. For 33% of the bla/ars in our sample, a description of the observed X-ray spectrum benefits from an additional component from proton interactions, in agreement with recent studies of individual IceCube candidate bla/ars. We show that, on average, bla/ars that are brighter in gigaelectronvolt gamma rays have a higher neutrino production efficiency but a lower best-fit baryonic loading. The predicted neutrino luminosity shows a positive correlation both with the observed flux of gigaelectronvolt gamma rays and with the predicted flux of megaelectronvolt gamma rays. We also estimate the diffuse neutrino flux from gamma-ray bla/ars by extrapolating the result to the Fermi population, and we show that it may be at the level of ~20% of the diffuse neutrino flux observed by IceCube, in agreement with current limits from stacking analyses. We discuss the implications of our results for future neutrino searches and suggest promising sources for potential detections.

A multi-cubic-kilometre neutrino telescope in the western Pacific Ocean

Next-generation neutrino telescopes with substantially improved sensitivity are required to pinpoint the sources of the diffuse astrophysical neutrino flux detected by IceCube and uncover the century-old puzzle of cosmic-ray origins. A detector near the Equator will provide a unique viewpoint of the neutrino sky, complementing IceCube and other neutrino telescopes in the Northern Hemisphere. Here we present results from an expedition to the northeastern region of the South China Sea, in the western Pacific Ocean. A favourable neutrino telescope site was found on an abyssal plain at a depth of ~3.5 km. At depths below 3 km, the sea current speed, water absorption and scattering lengths for Cherenkov light were measured to be vc < 10 cm s−1, λabs ≈ 27 m and λsca ≈ 63 m, respectively. Accounting for these measurements, we present the design and expected performance of a next-generation neutrino telescope, Tropical Deep-sea Neutrino Telescope (TRIDENT). With its advanced photon-detection technology and large dimensions, TRIDENT expects to observe the IceCube steady source candidate NGC 1068 with 5σ significance within 1 year of operation. This level of sensitivity will open a new arena for diagnosing the origin of cosmic rays and probing fundamental physics over astronomical baselines.

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.

Search for directional associations between baikal gigaton volume detector neutrino-induced cascades and high-energy astrophysical sources

ABSTRACT Baikal-GVD has recently published its first measurement of the diffuse astrophysical neutrino flux, performed using high-energy cascade-like events. We further explore the Baikal-GVD cascade data set collected in 2018–2022, with the aim to identify possible associations between the Baikal-GVD neutrinos and known astrophysical sources. We leverage the relatively high angular resolution of the Baikal-GVD neutrino telescope (2–3 deg.), made possible by the use of liquid water as the detection medium, enabling the study of astrophysical point sources even with cascade events. We estimate the telescope’s sensitivity in the cascade channel for high-energy astrophysical sources and refine our analysis prescriptions using Monte-Carlo simulations. We primarily focus on cascades with energies exceeding 100 TeV, which we employ to search for correlation with radio-bright blazars. Although the currently limited neutrino sample size provides no statistically significant effects, our analysis suggests a number of possible associations with both extragalactic and Galactic sources. Specifically, we present an analysis of an observed triplet of neutrino candidate events in the Galactic plane, focusing on its potential connection with certain Galactic sources, and discuss the coincidence of cascades with several bright and flaring blazars.

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.

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.

Energy-dependent flavor ratios, cascade/track spectrum tension and high-energy neutrinos from magnetospheres of supermassive black holes

The IceCube neutrino observatory measures the diffuse flux of high-energy astrophysical neutrinos by means of various techniques, and there exists a mild tension between spectra obtained in different analyses. The spectrum derived from reconstruction of muon tracks is harder than that from cascades, dominated by electron and tau neutrinos. If confirmed, this tension may provide a clue to the origin of these neutrinos, which remains uncertain. Here we investigate the possibility that this tension may be caused by the change of the flavor content of astrophysical neutrinos with energy. We assume that at higher energies, the flux contains more muon neutrinos than expected in the usually assumed flavor equipartition. This may happen if the neutrinos are produced in regions of the magnetic field so strong that muons, born in pi-meson decays, cool by synchrotron radiation faster than decay. The magnetic field of ∼104 G is required for this mechanism to be relevant for the IceCube results. We note that these field values are reachable in the immediate vicinity of supermassive black holes in active galactic nuclei and present a working toy model of the population of these potential neutrino sources. While this model predicts the required flavor ratios and describes the high-energy spectrum, it needs an additional component to explain the observed neutrino flux at lower energies.

Neutrino predictions from choked Gamma-Ray Bursts and comparison with the observed cosmic diffuse neutrino flux

The strong constraints from the Fermi-LAT data on the isotropic gamma-ray background suggest that the neutrinos observed by IceCube might possibly come from sources that are hidden to gamma-ray observations. A possibility emerged in recent years is that neutrinos may come from jets of collapsing massive stars which fail to break out of the stellar envelope, and for this reason they are known as choked jets, or choked Gamma-Ray Bursts (GRBs). We here show our predictions of neutrino flux and spectrum expected from these sources, focusing on Type II SNe, through detailed calculations of pγ interactions and accounting for all the neutrino production channels and scattering angles. We provide predictions of expected event rates for ANTARES, IceCube, and the next generation neutrino telescope KM3NeT.We also compute the contribution of the choked GRB population to the diffuse astrophysical neutrino flux, thus providing constraints on the local rate of this source population as to reproduce the observed neutrino flux.

Extragalactic neutrino-emission induced by supermassive and stellar mass black hole mergers

ABSTRACT The recent detections of binary stellar mass black hole mergers by the LIGO and Virgo Collaborations suggest that such mergers are common occurrences. Galaxy mergers further indicate that supermassive black holes in centres of galaxies also merge and are typically expected to have had at least one merger in their lifetime, possibly many. In the presence of a jet, these mergers are almost always accompanied by a change of the jet direction and a connected jet precession motion, leading to interactions of the jet with ambient matter and producing high-energy particles, and consequently high-energy gamma-rays and neutrinos. In this work, we investigate the possibility under which conditions such mergers could be the sources of the diffuse astrophysical neutrino flux measured by the IceCube Neutrino Observatory. The main free parameters in the calculation concern the frequency of the mergers and the fraction of energy that is transferred from the gravitationally released energy to neutrinos. We show that the merger rate for SMBBHs must lie between ∼10−7 and 10−5 Gpc−3 yr−1. The ratio of energy going to neutrinos during such mergers lies then between ∼10−6 and 3 × 10−4. For stellar mass BBH mergers, the rate needs to be ∼10–100 Gpc−3 yr−1 and the expected ratio of neutrino to GW energy lies in a comparable range as for SMBBHs, ∼2 × 10−5–10−3. These values lie in a reasonable parameter range, so that the production of neutrinos at the level of the detected neutrino flux is a realistic possibility.

Searches for Neutrinos from Gamma-Ray Bursts Using the IceCube Neutrino Observatory

Gamma-ray bursts (GRBs) are considered as promising sources of ultra-high-energy cosmic rays (UHECRs) due to their large power output. Observing a neutrino flux from GRBs would offer evidence that GRBs are hadronic accelerators of UHECRs. Previous IceCube analyses, which primarily focused on neutrinos arriving in temporal coincidence with the prompt gamma-rays, found no significant neutrino excess. The four analyses presented in this paper extend the region of interest to 14 days before and after the prompt phase, including generic extended time windows and targeted precursor searches. GRBs were selected between 2011 May and 2018 October to align with the data set of candidate muon-neutrino events observed by IceCube. No evidence of correlation between neutrino events and GRBs was found in these analyses. Limits are set to constrain the contribution of the cosmic GRB population to the diffuse astrophysical neutrino flux observed by IceCube. Prompt neutrino emission from GRBs is limited to ≲1% of the observed diffuse neutrino flux, and emission on timescales up to 104 s is constrained to 24% of the total diffuse flux.

ASAS-SN follow-up of IceCube high-energy neutrino alerts

ABSTRACT We report on the search for optical counterparts to IceCube neutrino alerts released between 2016 April and 2021 August with the All-Sky Automated Survey for SuperNovae (ASAS-SN). Despite the discovery of a diffuse astrophysical high-energy neutrino flux in 2013, the source of those neutrinos remains largely unknown. Since 2016, IceCube has published likely astrophysical neutrinos as public real-time alerts. Through a combination of normal survey and triggered target-of-opportunity observations, ASAS-SN obtained images within 1 h of the neutrino detection for 20 per cent (11) of all observable IceCube alerts and within one day for another 57 per cent (32). For all observable alerts, we obtained images within at least two weeks from the neutrino alert. ASAS-SN provides the only optical follow-up for about 17 per cent of IceCube’s neutrino alerts. We recover the two previously claimed counterparts to neutrino alerts, the flaring-blazar TXS 0506 + 056 and the tidal disruption event AT2019dsg. We investigate the light curves of previously detected transients in the alert footprints, but do not identify any further candidate neutrino sources. We also analysed the optical light curves of Fermi 4FGL sources coincident with high-energy neutrino alerts, but do not identify any contemporaneous flaring activity. Finally, we derive constraints on the luminosity functions of neutrino sources for a range of assumed evolution models.

Search for High-energy Neutrinos from Ultraluminous Infrared Galaxies with IceCube

Abstract Ultraluminous infrared galaxies (ULIRGs) have infrared luminosities L IR ≥ 1012 L ⊙, making them the most luminous objects in the infrared sky. These dusty objects are generally powered by starbursts with star formation rates that exceed 100 M ⊙ yr−1, possibly combined with a contribution from an active galactic nucleus. Such environments make ULIRGs plausible sources of astrophysical high-energy neutrinos, which can be observed by the IceCube Neutrino Observatory at the South Pole. We present a stacking search for high-energy neutrinos from a representative sample of 75 ULIRGs with redshift z ≤ 0.13 using 7.5 yr of IceCube data. The results are consistent with a background-only observation, yielding upper limits on the neutrino flux from these 75 ULIRGs. For an unbroken E −2.5 power-law spectrum, we report an upper limit on the stacked flux Φ ν μ + ν ¯ μ 90 % = 3.24 × 10 − 14 TeV − 1 cm − 2 s − 1 ( E / 10 TeV ) − 2.5 at 90% confidence level. In addition, we constrain the contribution of the ULIRG source population to the observed diffuse astrophysical neutrino flux as well as model predictions.

A new and improved IceCube point source analysis

The IceCube Neutrino Observatory, a cubic kilometer scale Cherenkov detector deployed in the deep ice at the geographic South Pole, investigates extreme astrophysical phenomena by studying the corresponding high-energy neutrino signal. Its discovery of a diffuse flux of astrophysical neutrinos with energies up to the PeV scale in 2013 has triggered a vast effort to identify the mostly unknown sources of these high energy neutrinos. Here, we present a new IceCube point-source search that improves the accuracy of the statistical analysis, especially at energies of a few TeV and below. The new approach is based on multidimensional kernel density estimation for the probability density functions and new estimators for the observables, namely the reconstructed energy and the estimated angular uncertainty on the reconstructed arrival direction. The more accurate analysis provides an improvement in discovery potential up to ∼30% over previous works for hard spectrum sources.

Constraining the contribution of Gamma-Ray Bursts to the high-energy diffuse neutrino flux with 10 years of ANTARES data

Gamma-Ray Bursts (GRBs) are considered very interesting astrophysical sources to be studied in the context of neutrino astronomy. Indeed, their jet composition continues to be an open issue. Within the framework of the fireball model, mesons can be produced during photo-hadronic interactions occurring in the internal shocks between shells emitted by the central engine. From their decays, high-energy gamma rays and neutrinos are expected to be generated. Neutrino telescopes are particularly interested in GRBs: being the most powerful explosions observable in the Universe, they are potentially able to achieve the energetics required to contribute to the diffuse astrophysical neutrino flux measured few years ago by IceCube. The analysis here presented relies on the search of time and space coincidences between neutrinos and the GRB emission. The undersea neutrino telescope ANTARES, operational in the Northern hemisphere since 2008 in its full configuration, plays an important role in the cosmic neutrino searches. In this contribution, the results of a stacked search for astrophysical muon neutrinos are presented for 784 GRBs in the period 2007-2017 using the ANTARES data. Given the absence of coincidences between a neutrino and a GRB in the considered sample, this analysis has allowed to constrain the contribution of the detected GRB population to the neutrino diffuse flux to be less than 10% around 100 TeV, within the context of the internal shock model. In addition, for the first time in this kind of analysis, the uncertainties on the not well characterised GRB parameters, needed to compute the neutrino flux expectations, are taken into account for each individual burst and then propagated to the diffuse flux from GRBs and to the estimated stacked upper limit.

Search for high-energy neutrino emission from radio-bright AGN

We investigate the possibility that radio-bright active galactic nuclei (AGN) are responsible for the TeV--PeV neutrinos detected by IceCube. We use an unbinned maximum-likelihood-ratio method, 10 yr of IceCube muon-track data, and 3388 radio-bright AGN selected from the Radio Fundamental Catalog. None of the AGN in the catalog have a large global significance. The two most significant sources have global significance of $\ensuremath{\simeq}1.5\ensuremath{\sigma}$ and $0.8\ensuremath{\sigma}$, though $4.1\ensuremath{\sigma}$ and $3.8\ensuremath{\sigma}$ local significance. Our stacking analyses show no significant correlation between the whole catalog and IceCube neutrinos. We infer from the null search that this catalog can account for at most 30% (95% C.L.) of the diffuse astrophysical neutrino flux measured by IceCube. Moreover, our results disagree with recent work that claimed a $4.1\ensuremath{\sigma}$ detection of neutrinos from the sources in this catalog, and we discuss the reasons of the difference.

Constraining the contribution of Gamma-Ray Bursts to the high-energy diffuse neutrino flux with 10 yr of ANTARES data

ABSTRACT Addressing the origin of the astrophysical neutrino flux observed by IceCube is of paramount importance. Gamma-Ray Bursts (GRBs) are among the few astrophysical sources capable of achieving the required energy to contribute to such neutrino flux through pγ interactions. In this work, ANTARES data have been used to search for upward going muon neutrinos in spatial and temporal coincidence with 784 GRBs occurred from 2007 to 2017. For each GRB, the expected neutrino flux has been calculated in the framework of the internal shock model and the impact of the lack of knowledge on the majority of source redshifts and on other intrinsic parameters of the emission mechanism has been quantified. It is found that the model parameters that set the radial distance where shock collisions occur have the largest impact on neutrino flux expectations. In particular, the bulk Lorentz factor of the source ejecta and the minimum variability time-scale are found to contribute significantly to the GRB-neutrino flux uncertainty. For the selected sources, ANTARES data have been analysed by maximizing the discovery probability of the stacking sample through an extended maximum-likelihood strategy. Since no neutrino event passed the quality cuts set by the optimization procedure, 90 per cent confidence level upper limits (with their uncertainty) on the total expected diffuse neutrino flux have been derived, according to the model. The GRB contribution to the observed diffuse astrophysical neutrino flux around 100 TeV is constrained to be less than 10 per cent.

Bounds on secret neutrino interactions from high-energy astrophysical neutrinos

Neutrinos offer a window to physics beyond the Standard Model. In particular, high-energy astrophysical neutrinos, with TeV-PeV energies, may provide evidence of new, "secret" neutrino-neutrino interactions that are stronger than ordinary weak interactions. During their propagation over cosmological distances, high-energy neutrinos could interact with the cosmic neutrino background via secret interactions, developing characteristic energy-dependent features in their observed energy distribution. For the first time, we look for signatures of secret neutrino interactions in the diffuse flux of high-energy astrophysical neutrinos, using 6 years of publicly available IceCube High Energy Starting Events (HESE). We find no significant evidence for secret neutrino interactions, but place competitive upper limits on the coupling strength of the new mediator through which they occur, in the mediator mass range of 1-100 MeV.

Interpretation of the Diffuse Astrophysical Neutrino Flux in Terms of the Blazar Sequence

Abstract We study whether the diffuse astrophysical neutrino flux can come from blazar jets—a subclass of active galactic nuclei—while it, at the same time, respects the blazar stacking limit based on source catalogs and is consistent with the observation from TXS 0506+056. We compute the neutrino flux from resolved and unresolved sources using an averaged, empirical relationship between electromagnetic spectrum and luminosity, known as the blazar sequence, for two populations of blazars (BL Lacs and FSRQs). Using a source model with realistic neutrino flux computations, we demonstrate that blazars can indeed power the diffuse neutrino flux at the highest energies and obey the stacking limit at the same time, and we derive the conditions for the baryonic loading (proton versus γ-ray luminosity) evolving over the blazar sequence. Under the hypothesis that low-luminosity blazars power the diffuse astrophysical neutrino flux, we find that the dominant contribution of the diffuse flux up to PeV energies must come from unresolved BL Lacs with baryonic loadings larger than about 105—while only a very small contribution may come from resolved high-luminosity BL Lacs or FSRQs, which can be directly tested by the stacking limit. We find that the blazar TXS 0506+056 is on the verge of these populations in our baseline scenario, at a relatively high luminosity and redshift; as a consequence, we predict about 0.3 γ-ray-neutrino associations per year from the whole population, dominated by BL Lacs with and z ∼ 0.1.