High-energy Neutrinos Research Articles

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.

IceCube Search for Neutrinos Coincident with Gravitational Wave Events from LIGO/Virgo Run O3

Using data from the IceCube Neutrino Observatory, we searched for high-energy neutrino emission from the gravitational-wave events detected by the advanced LIGO and Virgo detectors during their third observing run. We did a low-latency follow-up on the public candidate events released during the detectors’ third observing run and an archival search on the 80 confident events reported in the GWTC-2.1 and GWTC-3 catalogs. An extended search was also conducted for neutrino emission on longer timescales from neutron star containing mergers. Follow-up searches on the candidate optical counterpart of GW190521 were also conducted. We used two methods; an unbinned maximum likelihood analysis and a Bayesian analysis using astrophysical priors, both of which were previously used to search for high-energy neutrino emission from gravitational-wave events. No significant neutrino emission was observed by any analysis, and upper limits were placed on the time-integrated neutrino flux as well as the total isotropic equivalent energy emitted in high-energy neutrinos.

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.

The NUSES space mission

NUSES is a new space mission aiming to test innovative observational and technological approaches related to the study of low energy cosmic and gamma rays, high energy astrophysical neutrinos, Sun-Earth environment, Space weather and magnetosphere-ionosphere-lithosphere coupling (MILC). The satellite will host two payloads, named TERZINA and ZIRÉ. ZIRÉ will perform measurements of electrons, protons and light nuclei from few up to hundreds of MeV, also testing new tools for the detection of cosmic MeV photons, and monitoring of MILC signals. TERZINA will observe the Cherenkov light produced by EAS generated by cosmic ray primaries at very high energies and will monitor the light emissions from the Earth limb in the near UV and visible ranges at the ns time scale, thus testing the observational concept of detecting Earth skimming astrophysical high energy neutrinos. The scientific objectives and development status of the mission will be presented.

A new multimessenger study of Starburst galaxies: a closer look to expecting neutrinos

Considering the high star formation rate (up to 100 M⊙/year) of starburst galaxies (SBGs), they are well posed between the guaranteed “factories” of high energy neutrinos, since they can contain accelerated cosmic rays in the central region where the high-density gas is present. A more accurate description of their possible hadronic emission could help to better explain the diffuse astrophysical flux measured by IceCube as well as the observed point-like excess, like the case of NGC1068. With this in mind we report here a multi-messenger study, looking at diffuse and resolved gamma-ray and neutrino measurements, that explain the very-high-energy (VHE) emission through a calorimetric scenario. For the analysis of the diffuse component we perform a blending of the available spectral indexes and produce a multi-component study of extragalactic background light (EGB), high energy starting events (HESE) and high-energy cascade IceCube data. In contrast to common prototype scenarios, the spectral index blending leads to a non negligible diffuse neutrino component from SBGs, accounting up to 40% of the HESE events, at 95.4% CL. This scenario privileges also a maximal energy within tens of PeV for the accelerated charge particles inside these galaxies. For the point-like study we report the proposed calorimetric description for the gamma-ray resolved SBGs within 100 Mpc, taking into account the star formation rate derived from their infrared emission. These neutrino expectations are then compared with the sensitivity of IceCube, IceCube/Gen2 and the incoming KM3NeT.

Search for Gamma-Ray and Neutrino Coincidences Using HAWC and ANTARES Data

In the quest for high-energy neutrino sources, the Astrophysical Multimessenger Observatory Network has implemented a new search by combining data from the High Altitude Water Cherenkov (HAWC) Observatory and the Astronomy with a Neutrino Telescope and Abyss environmental RESearch (ANTARES) neutrino telescope. Using the same analysis strategy as in a previous detector combination of HAWC and IceCube data, we perform a search for coincidences in HAWC and ANTARES events that are below the threshold for sending public alerts in each individual detector. Data were collected between 2015 July and 2020 February with a live time of 4.39 yr. Over this time period, three coincident events with an estimated false-alarm rate of <1 coincidence per year were found. This number is consistent with background expectations.

Gamma-Ray Emission from Galaxies Hosting Molecular Outflows

Many star-forming galaxies and those hosting active galactic nuclei show evidence of massive outflows of material in a variety of phases including ionized, neutral atomic, and molecular. Molecular outflows in particular have been the focus of recent interest as they may be responsible for removing gas from the galaxy, thereby suppressing star formation. As material is ejected from the cores of galaxies, interactions of the outflowing material with the interstellar medium can accelerate cosmic rays and produce high-energy gamma rays. In this work, we search for gamma-ray emission from a sample of local galaxies known to host molecular outflows using data collected by the Fermi Large Area Telescope. We employ a stacking technique in order to search for and characterize the average gamma-ray emission properties of the sample. Gamma-ray emission is detected from the galaxies in our sample at the 4.4σ level with a power-law photon index of Γ ≈ 2 in the 1–800 GeV energy range. The emission is found to correlate with tracers of star formation activity, namely the 8–1000 μm infrared luminosity. We also find that the observed signal can be predominantly attributed to H ii galaxies hosting energy-driven outflows. While we do not find evidence suggesting that the outflows are accelerating charged particles directly, galaxies with molecular outflows may produce more gamma rays than galaxies without outflows. In particular, the set consisting of gamma-ray-detected galaxies with molecular outflows are nearly perfect calorimeters and may be future targets for searches of high-energy neutrinos.

Multi-messenger implications of the Pierre Auger Observatory measurements

The Pierre Auger Observatory is not only the world’s largest cosmic-ray observatory, but also an efficient detector of cosmic particles other than nuclei. Using its enormous exposure, the Observatory is a competitive player in searches for ultra-high energy neutrinos, photons, and even particles that could be predicted by theories beyond the standard model of particle physics. We present our recent results of such searches, targeting both the diffuse fluxes of particles and specific transient events. Special attention is given to the potential of the Observatory to support or constrain the observations of the ANITA experiment interpreted as upward-going air showers. The search for correlations of arrival directions of ultra-high energy cosmic rays measured by the Pierre Auger Observatory and the Telescope Array with high-energy neutrinos detected by IceCube and ANTARES experiments is also discussed.

High-energy neutrino astronomy — the neutrino connections to the cosmic-ray origin: present and future

High-energy neutrino astronomy has been blooming. In addition to the possible identification of the blazar and Seyfert II galaxies as neutrino emitters, the present neutrino data has indicated some hints to characterize or constrain the origin of cosmic rays. Being motivated by the observational fact that the astrophysical neutrino background energy flux is comparable to that of ultrahigh-energy (UHE) cosmic rays, we derive the generic requirements that a major fraction of UHE cosmic ray sources must meet, if they are also responsible for the ∼ 100 TeV-energy cosmic neutrino background radiation. The source parameters characterizing the cosmic ray – neutrino unified scheme, such as the photon radiation luminosity and the cosmic ray luminosity density, suggest that the yet-unidentified cosmic ray and neutrino origins can be transient objects visible in the optical/NIR wavelength band. We propose a viable scheme of multimessenger observations to identify the sources using the neutrino multiplet detections.

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.

Simulating neutrino echoes induced by secret neutrino interactions

New neutrino interactions beyond the Standard Model (BSM) have been of much interest in not only particle physics but also cosmology and astroparticle physics. We numerically investigate the time delay distribution of astrophysical neutrinos that interact with the cosmic neutrino background. Using the Monte Carlo method, we develop a framework that enables us to simulate the time-dependent energy spectra of high-energy neutrinos that experience even multiple scatterings en route and to handle the sharp increase in the cross section at the resonance energy. As an example, we focus on the case of secret neutrino interactions with a scalar mediator. While we find the excellent agreement between analytical and simulation results for small optical depths, our simulations enable us to study optically thick cases that are not described by the simplest analytic estimates. Our simulations are used to understand effects of cosmological redshifts, neutrino spectra and flavors. The developments will be useful for probing BSM neutrino interactions with not only current neutrino detectors such as IceCube and Super-Kamiokande but also future neutrino detectors such as IceCube-Gen2 and Hyper-Kamiokande.

Upper Limits on the Isotropic Diffuse Flux of Cosmic PeV Photons from Carpet-2 Observations

Isotropic diffuse gamma-ray flux in the PeV energy band is an important tool for multimessenger tests of models of the origin of high-energy astrophysical neutrinos and for new-physics searches. So far, this flux has not yet been observed. Carpet-2 is an air-shower experiment capable of detecting astrophysical gamma rays with energies above 0.1 PeV. Here we report the upper limits on the isotropic gamma-ray flux from Carpet-2 data obtained in 1999-2011 and 2018-2022. These results, obtained with the new statistical method based on the shape of the muon-number distribution, summarize Carpet-2 observations as the upgraded installation, Carpet-3, starts its operation.

Multiwavelength Analysis of the IceCube Neutrino Source Candidate Blazar PKS 1424+240

The true nature of sources of cosmic neutrinos recorded by the Antarctic IceCube Neutrino Detector is still an enigma of high-energy astrophysics. Time-integrated neutrino source searches with the 10 years of IceCube data unfolded neutrino hot-spots of the sky; among them, one is associated with the blazar PKS 1424+240, which is the third most significant neutrino source candidate in the Northern sky. In this paper, we analyze VLBI radio data of PKS 1424+240 taken with the Very Large Baseline Array at 15 GHz as part of the MOJAVE Survey. We generate the adaptively binned gamma-ray light curve of the source, employing Fermi-LAT data between 100 MeV and 300 GeV. We find that the VLBI jet components maintain quasi-stationary core separations at 15 GHz. We find a quiescence and a perturbed phase of the VLBI core of PKS 1424+240, based on that its Doppler factor increased tenfold after 2016 compared to the quiescence phase. We do not find elevated gamma-ray activity after 2016, while archive Swift-XRT measurements show a highly increased 0.3–10 keV X-ray flux in the beginning of 2017. Substantial increase of the activity of the radio core might help us to identify episodes of particle acceleration in lepto-hadronic blazar jets that eventually lead to the emission of high-energy neutrinos.

Strong lensing of high-energy neutrinos

We consider the effects of strong gravitational lensing by galaxy-scale deflectors on the observations of high-energy ($\mathrm{E}\ensuremath{\gg}\mathrm{GeV}$) neutrinos (HENs). For HENs at cosmological distances, the optical depth for multiple imaging is $\ensuremath{\sim}{10}^{\ensuremath{-}3}$, implying that while we do not expect any multiply imaged HENs with present samples, next-generation experiments should be able to detect the first such event. We then present the distribution of expected time delays to aid in the identification of such events, in combination with directional and energy information. In order to assist in the evaluation of HEN production mechanisms, we illustrate how lensing affects the observed number counts for a variety of intrinsic luminosity functions of the source population. Finally, we see that the lensing effects on the cosmic neutrino background flux calculation would be negligible by taking kiloparsec-scale jets as an example.

Multi-Wavelength and Multi-Messenger Studies Using the Next-Generation Event Horizon Telescope

The next-generation Event Horizon Telescope (ngEHT) will provide us with the best opportunity to investigate supermassive black holes (SMBHs) at the highest possible resolution and sensitivity. With respect to the existing Event Horizon Telescope (EHT) array, the ngEHT will provide increased sensitivity and uv-coverage (with the addition of new stations), wider frequency coverage (from 86 GHz to 345 GHz and higher), finer resolution (&lt;15 micro-arcseconds), and better monitoring capabilities. The ngEHT will offer a unique opportunity to deeply investigate the physics around SMBHs, such as the disk-jet connection, the mechanisms responsible for high-energy photon and neutrino events, and the role of magnetic fields in shaping relativistic jets, as well as the nature of binary SMBH systems. In this white paper we describe some ngEHT science cases in the context of multi-wavelength studies and synergies.

Testing hadronic and photo-hadronic interactions as responsible for UHECR and neutrino fluxes from Starburst Galaxies

We test the hypothesis that starburst galaxies are the sources of ultra-high energy cosmic rays and high-energy neutrinos. The computation of interactions of ultra-high energy cosmic rays in the starburst environment as well as in the propagation to the Earth is made using a modified version of the Monte Carlo code SimProp, where hadronic processes are implemented for the first time. Taking into account a star-formationrate distribution of sources, the fluxes of ultra-high energy cosmic rays and high-energy neutrinos are computed and compared with observations, and the explored parameter space for the source characteristics is discussed. We find that, depending on the density of the gas in the source environment, spallation reactions could hide the outcome in neutrinos from photo-hadronic interactions in the source environment and in extra-galactic space. We confirm that source-propagation models constitute a promising way to improve the discrimination power of models considering only ultra-high energy cosmic rays, on the way to unveiling the source class responsible for ultra-high energy cosmic rays and high-energy neutrinos.

High-energy Neutrino Emission from Espresso-reaccelerated Ions in Jets of Active Galactic Nuclei

We present a bottom-up calculation of the flux of ultrahigh-energy cosmic rays (UHECRs) and high-energy neutrinos produced by powerful jets of active galactic nuclei (AGNs). By propagating test particles in 3D relativistic magnetohydrodynamic jet simulations, including a Monte Carlo treatment of sub-grid pitch-angle scattering and attenuation losses due to realistic photon fields, we study the spectrum and composition of the accelerated UHECRs and estimate the amount of neutrinos produced in such sources. We find that UHECRs may not be significantly affected by photodisintegration in AGN jets where the espresso mechanism efficiently accelerates particles, consistent with Auger’s results that favor a heavy composition at the highest energies. Moreover, we present estimates and upper bounds for the flux of high-energy neutrinos expected from AGN jets. In particular, we find that (i) source neutrinos may account for a sizable fraction, or even dominate, the expected flux of cosmogenic neutrinos; (ii) neutrinos from the β-decay of secondary neutrons produced in nucleus photodisintegration end up in the teraelectronvolt to petaelectronvolt band observed by IceCube, but can hardly account for the observed flux; (iii) UHECRs accelerated via the espresso mechanism lead to nearly isotropic neutrino emission, which suggests that nearby radio galaxies may be more promising as potential sources. We discuss our results in light of multimessenger astronomy and current/future neutrino experiments.

Searching for high-energy neutrinos from the most luminous supernovae with the IceCube Neutrino Observatory

The sources of the astrophysical neutrino flux discovered by IceCube remain for the most part unresolved. Extragalactic core-collapse supernovae (CCSNe) have been suggested as potentially able to produce high-energy neutrinos. In recent years, the Zwicky Transient Facility has discovered a population of exceptionally luminous supernovae, whose powering mechanisms have not yet been fully established. A fraction of these objects fall in the broader category of type IIn CCSNe, showing signs of interaction with a dense circumstellar medium. Theoretical models connect the supernova photometric properties to the dynamics of a shock-powered emission, predicting particle acceleration. In this contribution, we outline the plan for a search of high-energy neutrinos targeting the population of superluminous and type IIn supernovae with the IceCube Neutrino Observatory.

Machine learning application for dark matter - background classification in JUNO experiment

Jiangmen underground neutrino observatory (JUNO) has a potential to indirectly detect dark matter (DM), observing neutrino events from annihilations of DM trapped by the gravitational force in the solar core. Weakly interacting massive particle (WIMP) DM candidate with mass > 4 GeV has significant solar capture rate. In this work, we simulate JUNO neutrino events from the most dominated WIMP annihilation channel, . Given the high-energy neutrinos from massive WIMPs, we extend the neutrino-nucleon interactions in the detector to include quasi-elastic (QE) and deep inelastic scatterings (DIS). The most challenging background events in the energy range above 100 MeV is the atmospheric neutrinos (atmo). The pulse shape discrimination (PSD) method is usually applied to distinguish between DM and atmo. In this work, we apply machine learning (ML) techniques to distinguish between the background (BG) atmo and DM neutrino events. We found that the random forest algorithm gives the best results. Using ML, the accuracy of DM-atmo events classification is 93.8% while similar f1-score could be achieved. On the other hand, only 78.1 % accuracy is obtained using conventional PSD with linear regression (LR). The upper limit for spin-dependent DM-proton cross-section with the applied ML is 5.84 × 10−40 cm2 for a 10 GeV WIMP DM which is an improvement over the sensitivity of the LR method of by a factor of ∼ 2.8x.

EUSO-SPB2: A balloon experiment for UHECR and VHE neutrino observation

The Extreme Universe Space Observatory on a Super Pressure Balloon 2 (EUSO-SPB2) experiment will make new measurements from suborbital space as a precursor for future space missions that will address the challenge of the extremely low fluxes of ultra-high energy cosmic rays (UHECR) and very high energy (VHE) neutrinos. The EUSO-SPB2 detector is comprised of two 1m diameter aperture telescopes. The Fluorescence Telescope (FT) will point in nadir and will record fluorescence light from cosmic ray EAS with energies above 1 EeV in its field of view of 36 by 12 degrees. The Cherenkov Telescope (CT) features a silicon photomultiplier focal surface with a field of view of 12 by 6 degrees. The CT will switch between two observation modes: one which points the CT above the limb to measure the Cherenkov emission of cosmic ray EAS with energies above 1 PeV and one which points the CT below the limb to record the Cherenkov emission produced by PeV scale EAS initiated by neutrino-sourced tau decay. As it is the first time such an instrument has been flown, one of the priorities of the CT will be the study of the optical background for observing neutrinos in this way. EUSO-SPB2 is undergoing the final integration steps for launch on a NASA super pressure balloon payload in the spring of 2023 from Wanaka NZ. The CT was field-tested in March 2022 and the FT in August - September 2022. The data collected during the EUSO-SPB2 mission will be essential to advancing the development of a space-based multi-messenger observatory such as the Probe of Extreme Multi-Messenger Astrophysics (POEMMA). In this contribution, we discuss the EUSO-SPB2 science goals, the instruments, the expected performance and the current status.