# Publications

Papers from the IceCube Collaboration

# With Major Contribution - Published and Upcoming

## Study of mass composition of cosmic rays with IceTop and IceCube

The IceCube Neutrino Observatory is a multi-component detector at the South Pole which detects high-energy particles emerging from astrophysical events. These particles provide us with insights into the fundamental properties and behaviour of their sources. Besides its principal usage and merits in neutrino astronomy, using IceCube in conjunction with its surface array, IceTop, also makes it a unique three-dimensional cosmic-ray detector. This distinctive feature helps facilitate detailed cosmic-ray analysis in the transition region from galactic to extragalactic sources. We will present the progress made on multiple fronts to establish a framework for mass-estimation of primary cosmic rays. The first technique relies on a likelihood-based analysis of the surface signal distribution and improves upon the standard reconstruction technique. The second uses advanced methods in graph neural networks to use the full in-ice shower footprint, in addition to global shower-footprint features from IceTop. A comparison between the two methods for composition analysis as well as a possible extension of the analysis techniques for sub-PeV cosmic-ray air-showers will also be discussed.

## Online Masterclass built on the KASCADE Cosmic ray Data Centre

During the ongoing Covid-19 pandemic, people all over the world were forced to think about new ways of interacting with each other and this has especially challenged academics in their outreach activities with pupils. New online formats needed to be developed, and we used this opportunity to design and implement an (not only) online Masterclass using data from the KASCADE experiment. The masterclass is built on the KASCADE Cosmic Ray Data Centre and uses Jupyterhub and Notebooks for data analysis. We gained first practical experience during the International Cosmic Day with students at the age of 14-19 years. The Masterclass includes lectures on cosmic ray physics and data analysis, which are then consolidated in a hands-on part. By performing a cosmic-ray composition analysis on KASCADE data, the participants gain experience in using the KCDC open data web platform, working in the Jupyter environment, preprocessing data from a real astroparticle physics experiment, programming Python and performing exploratory data analysis.

## The Global Cosmic Ray Observatory - GCOS

Nature is providing particles with energies exceeding 100 EeV. Their existence imposes immediate questions: Are they ordinary particles, accelerated in extreme astrophysical environments, or are they annihilation or decay products of super-heavy dark matter or other exotic objects? If the particles are accelerated in extreme astrophysical environments, are their sources related to those of high-energy neutrinos, gamma rays, and/or gravitational waves, such as the recently observed mergers of compact objects? The particles can also be used to study physics processes at extreme energies; is Lorentz invariance still valid? Are the particles interacting according to the Standard Model or are there new physics processes? The particles can be used to study hadronic interactions (QCD) in the kinematic forward direction; what is the cross section of protons at center-of-mass energies s√>100~TeV? These questions are addressed at present by installations like the Telescope Array and the Pierre Auger Observatory. After the year 2030, a next-generation observatory will be needed to study the physics and properties of the highest-energy particles in Nature, building on the knowledge harvested from the existing observatories. It should have an aperture at least an order of magnitude bigger than the existing observatories.Recently, more than 200 scientists from around the world came together to discuss the future of the field of multi-messenger astroparticle physics beyond the year 2030.Ideas have been discussed towards the physics case and possible scenarios for detection concepts of the Global Cosmic Ray Observatory - GCOS. A synopsis of the key results discussed during the brainstorming workshop will be presented.

## Graph Neural Networks and application for Cosmic-Ray Analysis

Deep Learning has emerged as one of the most promising areas of computational research for pattern learning, inference drawing, and decision-making, with wide-ranging applications across various scientific disciplines. This has also made it possible for faster and more precise analysis in astroparticle physics, enabling new insights from massive volumes of input data. Graph Neural Networks have naturally developed as a critical implementation method among the numerous deep-learning architectures over the last few years because of the unique ability to represent complex input data from a wide range of problems in its most natural form. Described using nodes and edges, graphs allow us to efficiently represent relational data and learn hidden representations of input data to obtain better model accuracy. At IceCube Neutrino Observatory, a complex multi-component detector, traditional likelihood-based analysis on a per-event basis to reconstruct cosmic-ray air shower parameters is time-consuming and computationally costly. Using advanced and flexible models based on Graph Neural Networks naturally emerges as a possible solution, reducing the time and computing cost for performing such analysis while boosting sensitivity. This paper will outline Graph Neural Networks and discuss a possible application of using such methods at the IceCube Observatory.

Upcoming

# Published

## A Muon-Track Reconstruction exploiting Stochastic Losses for Large-Scale Cherenkov Detectors

IceCube is a cubic-kilometer Cherenkov telescope operating at the South Pole. The main goal of IceCube is the detection of astrophysical neutrinos and the identification of their sources. High-energy muon neutrinos are observed via the secondary muons produced in charge current interactions with nuclei in the ice. Currently, the best performing muon track directional reconstruction is based on a maximum likelihood method using the arrival time distribution of Cherenkov photons registered by the experiment's photomultipliers. A known systematic shortcoming of the prevailing method is to assume a continuous energy loss along the muon track. However at energies >1 TeV the light yield from muons is dominated by stochastic showers. This paper discusses a generalized ansatz where the expected arrival time distribution is parametrized by a stochastic muon energy loss pattern. This more realistic parametrization of the loss profile leads to an improvement of the muon angular resolution of up to 20% for through-going tracks and up to a factor 2 for starting tracks over existing algorithms. Additionally, the procedure to estimate the directional reconstruction uncertainty has been improved to be more robust against numerical errors.

## Search for GeV Neutrino Emission During Intense Gamma-Ray Solar Flares with the IceCube Neutrino Observatory

Solar flares convert magnetic energy into thermal and non-thermal plasma energy, the latter implying particle acceleration of charged particles such as protons. Protons are injected out of the coronal acceleration region and can interact with dense plasma in the lower solar atmosphere, producing mesons that subsequently decay into gamma rays and neutrinos at O(MeV-GeV) energies. We present the results of the first search for GeV neutrinos emitted during solar flares carried out with the IceCube Neutrino Observatory. While the experiment was originally designed to detect neutrinos with energies between 10 GeV and a few PeV, a new approach allowing for a O(GeV) energy threshold will be presented. The resulting limits allow us to constrain some of the theoretical estimates of the expected neutrino flux.

## Detection of a particle shower at the Glashow resonance with IceCube

The Glashow resonance describes the resonant formation of a W− boson during the interaction of a high-energy electron antineutrino with an electron1, peaking at an antineutrino energy of 6.3 petaelectronvolts (PeV) in the rest frame of the electron. Whereas this energy scale is out of reach for currently operating and future planned particle accelerators, natural astrophysical phenomena are expected to produce antineutrinos with energies beyond the PeV scale. Here we report the detection by the IceCube neutrino observatory of a cascade of high-energy particles (a particle shower) consistent with being created at the Glashow resonance. A shower with an energy of 6.05 ± 0.72 PeV (determined from Cherenkov radiation in the Antarctic Ice Sheet) was measured. Features consistent with the production of secondary muons in the particle shower indicate the hadronic decay of a resonant W− boson, confirm that the source is astrophysical and provide improved directional localization. The evidence of the Glashow resonance suggests the presence of electron antineutrinos in the astrophysical flux, while also providing further validation of the standard model of particle physics. Its unique signature indicates a method of distinguishing neutrinos from antineutrinos, thus providing a way to identify astronomical accelerators that produce neutrinos via hadronuclear or photohadronic interactions, with or without strong magnetic fields. As such, knowledge of both the flavour (that is, electron, muon or tau neutrinos) and charge (neutrino or antineutrino) will facilitate the advancement of neutrino astronomy.

## A search for time-dependent astrophysical neutrino emission with IceCube data from 2012 to 2017

High-energy neutrinos are unique messengers of the high-energy universe, tracing the processes of cosmic-ray acceleration. This paper presents analyses focusing on time-dependent neutrino point-source searches. A scan of the whole sky, making no prior assumption about source candidates, is performed, looking for a space and time clustering of high-energy neutrinos in data collected by the IceCube Neutrino Observatory between 2012 and 2017. No statistically significant evidence for a time-dependent neutrino signal is found with this search during this period since all results are consistent with the background expectation. Within this study period, the blazar 3C 279, showed strong variability, inducing a very prominent gamma-ray flare observed in 2015 June. This event motivated a dedicated study of the blazar, which consists of searching for a time-dependent neutrino signal correlated with the gamma-ray emission. No evidence for a time-dependent signal is found. Hence, an upper limit on the neutrino fluence is derived, allowing us to constrain a hadronic emission model.

## Measurement of the high-energy all-flavor neutrino-nucleon cross section with IceCube

The flux of high-energy neutrinos passing through the Earth is attenuated due to their interactions with matter. The interaction rate is modulated by the neutrino interaction cross section and affects the flux arriving at the IceCube Neutrino Observatory, a cubic-kilometer neutrino detector embedded in the Antarctic ice sheet. We present a measurement of the neutrino cross section between 60 TeV and 10 PeV using the high-energy starting events (HESE) sample from IceCube with 7.5 years of data. The result is binned in neutrino energy and obtained using both Bayesian and frequentist statistics. We find it compatible with predictions from the Standard Model. Flavor information is explicitly included through updated morphology classifiers, proxies for the the three neutrino flavors. This is the first such measurement to use the three morphologies as observables and the first to account for neutrinos from tau decay.

## Follow-up of astrophysical transients in real time with the IceCube Neutrino Observatory

In multi-messenger astronomy, rapid investigation of interesting transients is imperative. As an observatory with a 4*pi steradian field of view and approx. 99% uptime, the IceCube Neutrino Observatory is a unique facility to follow up transients, and to provide valuable insight for other observatories and inform their observing decisions. Since 2016, IceCube has been using low-latency data to rapidly respond to interesting astrophysical events reported by the multi-messenger observational community. Here, we describe the pipeline used to perform these follow up analyses and provide a summary of the 58 analyses performed as of July 2020. We find no significant signal in the first 58 analyses performed. The pipeline has helped inform various electromagnetic observing strategies, and has constrained neutrino emission from potential hadronic cosmic accelerators.

## Multimessenger Gamma-Ray and Neutrino Coincidence Alerts using HAWC and IceCube sub-threshold Data

The High Altitude Water Cherenkov (HAWC) and IceCube observatories, through the Astrophysical Multimessenger Observatory Network (AMON) framework, have developed a multimessenger joint search for extragalactic astrophysical sources. This analysis looks for sources that emit both cosmic neutrinos and gamma rays that are produced in photo-hadronic or hadronic interactions. The AMON system is running continuously, receiving sub-threshold data (i.e. data that is not suited on its own to do astrophysical searches) from HAWC and IceCube, and combining them in real-time. We present here the analysis algorithm, as well as results from archival data collected between June 2015 and August 2018, with a total live-time of 3.0 years. During this period we found two coincident events that have a false alarm rate (FAR) of < 1 coincidence per year, consistent with the background expectations. The real-time implementation of the analysis in the AMON system began on November 20th, 2019, and issues alerts to the community through the Gamma-ray Coordinates Network with a FAR threshold of < 4 coincidences per year.

## Multimessenger Gamma-Ray and Neutrino Coincidence Alerts Using HAWC and IceCube Sub-threshold Data

The High Altitude Water Cherenkov (HAWC) and IceCube observatories, through the Astrophysical Multimessenger Observatory Network (AMON) framework, have developed a multimessenger joint search for extragalactic astrophysical sources. This analysis looks for sources that emit both cosmic neutrinos and gamma rays that are produced in photo-hadronic or hadronic interactions. The AMON system is running continuously, receiving sub-threshold data (i.e. data that is not suited on its own to do astrophysical searches) from HAWC and IceCube, and combining them in real-time. We present here the analysis algorithm, as well as results from archival data collected between June 2015 and August 2018, with a total live-time of 3.0 years. During this period we found two coincident events that have a false alarm rate (FAR) of less than 1 coincidence per year, consistent with the background expectations. The real-time implementation of the analysis in the AMON system began on November 20th, 2019, and issues alerts to the community through the Gamma-ray Coordinates Network with a FAR threshold of less than 4 coincidences per year.

## Measurements of the Time-Dependent Cosmic-Ray Sun Shadow with Seven Years of IceCube Data -- Comparison with the Solar Cycle and Magnetic Field Models

Observations of the time-dependent cosmic-ray Sun shadow have been proven as a valuable diagnostic for the assessment of solar magnetic field models. In this paper, seven years of IceCube data are compared to solar activity and solar magnetic field models. A quantitative comparison of solar magnetic field models with IceCube data on the event rate level is performed for the first time. Additionally, a first energy-dependent analysis is presented and compared to recent predictions. We use seven years of IceCube data for the Moon and the Sun and compare them to simulations on data rate level. The simulations are performed for the geometrical shadow hypothesis for the Moon and the Sun and for a cosmic-ray propagation model governed by the solar magnetic field for the case of the Sun. We find that a linearly decreasing relationship between Sun shadow strength and solar activity is preferred over a constant relationship at the 6.4sigma level. We test two commonly used models of the coronal magnetic field, both combined with a Parker spiral, by modeling cosmic-ray propagation in the solar magnetic field. Both models predict a weakening of the shadow in times of high solar activity as it is also visible in the data. We find tensions with the data on the order of 3σ for both models, assuming only statistical uncertainties. The magnetic field model CSSS fits the data slightly better than the PFSS model. This is generally consistent with what is found previously by the Tibet AS-gamma Experiment, a deviation of the data from the two models is, however, not significant at this point. Regarding the energy dependence of the Sun shadow, we find indications that the shadowing effect increases with energy during times of high solar activity, in agreement with theoretical predictions.

## Cosmic Ray Spectrum from 250 TeV to 10 PeV using IceTop

We report here an extension of the measurement of the all-particle cosmic-ray spectrum with IceTop to lower energy. The new measurement gives full coverage of the knee region of the spectrum and reduces the gap in energy between previous IceTop and direct measurements. With a new trigger that selects events in closely spaced detectors in the center of the array, the IceTop energy threshold is lowered by almost an order of magnitude below its previous threshold of 2 PeV. In this paper, we explain how the new trigger is implemented, and we describe the new machine-learning method developed to deal with events with very few detectors hit. We compare the results with previous measurements by IceTop and others that overlap at higher energy and with HAWC and Tibet in the 100 TeV range.

## Searching for eV-scale sterile neutrinos with eight years of atmospheric neutrinos at the IceCube neutrino telescope

We report in detail on searches for eV-scale sterile neutrinos, in the context of a 3+1 model, using eight years of data from the IceCube neutrino telescope. By analyzing the reconstructed energies and zenith angles of 305,735 atmospheric νμ and ν¯μ events we construct confidence intervals in two analysis spaces: sin2(2θ24) vs. Δm241 under the conservative assumption θ34=0; and sin2(2θ24) vs. sin2(2θ34) given sufficiently large Δm241 that fast oscillation features are unresolvable. Detailed discussions of the event selection, systematic uncertainties, and fitting procedures are presented. No strong evidence for sterile neutrinos is found, and the best-fit likelihood is consistent with the no sterile neutrino hypothesis with a p-value of 8% in the first analysis space and 19\% in the second.

## An eV-scale sterile neutrino search using eight years of atmospheric muon neutrino data from the IceCube Neutrino Observatory

The results of a 3+1 sterile neutrino search using eight years of data from the IceCube Neutrino Observatory are presented. A total of 305,735 muon neutrino events are analyzed in reconstructed energy-zenith space to test for signatures of a matter-enhanced oscillation that would occur given a sterile neutrino state with a mass-squared differences between 0.01\,eV2 and 100\,eV2. The best-fit point is found to be at sin2(2θ24)=0.10 and Δm241=4.5eV2, which is consistent with the no sterile neutrino hypothesis with a p-value of 8.0%.

## IceCube Search for Neutrinos Coincident with Compact Binary Mergers from LIGO-Virgo's First Gravitational-Wave Transient Catalog

Using the IceCube Neutrino Observatory, we search for high-energy neutrino emission coincident with compact binary mergers observed by the LIGO and Virgo gravitational wave (GW) detectors during their first and second observing runs. We present results from two searches targeting emission coincident with the sky localization of each gravitational wave event within a 1000 second time window centered around the reported merger time. One search uses a model-independent unbinned maximum likelihood analysis, which uses neutrino data from IceCube to search for point-like neutrino sources consistent with the sky localization of GW events. The other uses the Low-Latency Algorithm for Multi-messenger Astrophysics, which incorporates astrophysical priors through a Bayesian framework and includes LIGO-Virgo detector characteristics to determine the association between the GW source and the neutrinos. No significant neutrino coincidence is seen by either search during the first two observing runs of the LIGO-Virgo detectors. We set upper limits on the time-integrated neutrino emission within the 1000 second window for each of the 11 GW events. These limits range from 0.02-0.7 GeV cm−2. We also set limits on the total isotropic equivalent energy, Eiso, emitted in high-energy neutrinos by each GW event. These limits range from 1.7 × 10^51 - 1.8 × 10^55 erg. We conclude with an outlook for LIGO-Virgo observing run O3, during which both analyses are running in real time.

# Submitted / Under Review

## All-flavor constraints on nonstandard neutrino interactions and generalized matter potential with three years of IceCube DeepCore data

We report constraints on nonstandard neutrino interactions (NSI) from the observation of atmospheric neutrinos with IceCube, limiting all individual coupling strengths from a single dataset. Furthermore, IceCube is the first experiment to constrain flavor-violating and nonuniversal couplings simultaneously. Hypothetical NSI are generically expected to arise due to the exchange of a new heavy mediator particle. Neutrinos propagating in matter scatter off fermions in the forward direction with negligible momentum transfer. Hence the study of the matter effect on neutrinos propagating in the Earth is sensitive to NSI independently of the energy scale of new physics. We present constraints on NSI obtained with an all-flavor event sample of atmospheric neutrinos based on three years of IceCube DeepCore data. The analysis uses neutrinos arriving from all directions, with reconstructed energies between 5.6 GeV and 100 GeV. We report constraints on the individual NSI coupling strengths considered singly, allowing for complex phases in the case of flavor-violating couplings. This demonstrates that IceCube is sensitive to the full NSI flavor structure at a level competitive with limits from the global analysis of all other experiments. In addition, we investigate a generalized matter potential, whose overall scale and flavor structure are also constrained.

## Search for High-Energy Neutrinos from Ultra-Luminous Infrared Galaxies with IceCube

Ultra-luminous infrared galaxies (ULIRGs) have infrared luminosities LIR≥1012L⊙, 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 years 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 Convolutional Neural Network based Cascade Reconstruction for the IceCube Neutrino Observatory

Continued improvements on existing reconstruction methods are vital to the success of high-energy physics experiments, such as the IceCube Neutrino Observatory. In IceCube, further challenges arise as the detector is situated at the geographic South Pole where computational resources are limited. However, to perform real-time analyses and to issue alerts to telescopes around the world, powerful and fast reconstruction methods are desired. Deep neural networks can be extremely powerful, and their usage is computationally inexpensive once the networks are trained. These characteristics make a deep learning-based approach an excellent candidate for the application in IceCube. A reconstruction method based on convolutional architectures and hexagonally shaped kernels is presented. The presented method is robust towards systematic uncertainties in the simulation and has been tested on experimental data. In comparison to standard reconstruction methods in IceCube, it can improve upon the reconstruction accuracy, while reducing the time necessary to run the reconstruction by two to three orders of magnitude.

## Search for GeV Neutrino Emission During Intense Gamma-Ray Solar Flares with the IceCube Neutrino Observatory

Solar flares convert magnetic energy into thermal and non-thermal plasma energy, the latter implying particle acceleration of charged particles such as protons. Protons are injected out of the coronal acceleration region and can interact with dense plasma in the lower solar atmosphere, producing mesons that subsequently decay into gamma rays and neutrinos at O(MeV-GeV) energies. We present the results of the first search for GeV neutrinos emitted during solar flares carried out with the IceCube Neutrino Observatory. While the experiment was originally designed to detect neutrinos with energies between 10 GeV and a few PeV, a new approach allowing for a O(GeV) energy threshold will be presented. The resulting limits allow us to constrain some of the theoretical estimates of the expected neutrino flux.

## LeptonInjector and LeptonWeighter: A neutrino event generator and weighter for neutrino observatories

We present a high-energy neutrino event generator, called LeptonInjector, alongside an event weighter, called LeptonWeighter. Both are designed for large-volume Cherenkov neutrino telescopes such as IceCube. The neutrino event generator allows for quick and flexible simulation of neutrino events within and around the detector volume, and implements the leading Standard Model neutrino interaction processes relevant for neutrino observatories: neutrino-nucleon deep-inelastic scattering and neutrino-electron annihilation. In this paper, we discuss the event generation algorithm, the weighting algorithm, and the main functions of the publicly available code, with examples.

## A search for time-dependent astrophysical neutrino emission with IceCube data from 2012 to 2017

High-energy neutrinos are unique messengers of the high-energy universe, tracing the processes of cosmic-ray acceleration. This paper presents analyses focusing on time-dependent neutrino point-source searches. A scan of the whole sky, making no prior assumption about source candidates, is performed, looking for a space and time clustering of high-energy neutrinos in data collected by the IceCube Neutrino Observatory between 2012 and 2017. No statistically significant evidence for a time-dependent neutrino signal is found with this search during this period since all results are consistent with the background expectation. Within this study period, the blazar 3C 279, showed strong variability, inducing a very prominent gamma-ray flare observed in 2015 June. This event motivated a dedicated study of the blazar, which consists of searching for a time-dependent neutrino signal correlated with the gamma-ray emission. No evidence for a time-dependent signal is found. Hence, an upper limit on the neutrino fluence is derived, allowing us to constrain a hadronic emission model.

## Measurement of Astrophysical Tau Neutrinos in IceCube's High-Energy Starting Events

We present the results of a search for astrophysical tau neutrinos in 7.5 years of IceCube's high-energy starting event data. At high energies, two energy depositions stemming from the tau neutrino charged-current interaction and subsequent tau lepton decay may be resolved. We report the first detection of two such events, with probabilities of ∼76% and ∼98% of being produced by astrophysical tau neutrinos. The resultant astrophysical neutrino flavor measurement is consistent with expectations, disfavoring a no-astrophysical tau neutrino flux scenario with 2.8σ significance.

## Search for sub-TeV neutrino emission from transient sources with three years of IceCube data

Since the discovery of a flux of high-energy astrophysical neutrinos, searches for their origins have focused primarily at TeV-PeV energies. Compared to sub-TeV searches, high-energy searches benefit from an increase in the neutrino cross section, improved angular resolution on the neutrino direction, and a reduced background from atmospheric neutrinos and muons. However, the focus on high energy does not preclude the existence of subTeV neutrino emission where IceCube retains sensitivity. Here we present the first all-flavor search from IceCube for transient emission of low-energy neutrinos, between 1-100 GeV using three years of data obtained with the IceCube-DeepCore detector. We find no evidence of transient neutrino emission in the data, thus leading to a constraint on the volumetric rate of astrophysical transient sources in the range of approximately 705 − 2301 Gpc−3 yr−1 for sources following a subphotospheric energy spectrum with a mean energy of 100 GeV and a bolometric energy of 10^52 erg