CORSIKA 8 is a new framework for simulations of particle cascades in air and dense media implemented in modern C++17, based on past experience with existing codes, in particular CORSIKA 7. The flexible and modular structure of the project allows the development of independent modules that can produce a fully customizable particle shower simulation. The radio module in particular is designed to treat the electric field calculation and its propagation through complex media to each observer location in an autonomous and flexible way. It already allows for the simultaneous simulation of the radio emission calculated with two independent time-domain formalisms, the "Endpoint formalism" as previously implemented in CoREAS and the "ZHS" algorithm as ported from ZHAireS. The design acts as the baseline interface for current and future development for the simulation of radio emission from particle showers in standard and complex scenarios, such as cross-media showers penetrating from air into ice. In this work, we present the design and implementation of the radio module in CORSIKA 8, along with validation studies and a direct comparison of the radio emission from air showers simulated with CORSIKA 8, CORSIKA 7 and ZHAireS. We also present the impact of simulation details such as the step size of simulated particle tracks on radio-emission simulations and perform a direct comparison of the "Endpoints" and "ZHS" formalisms for the same underlying air showers. Finally, we present an in-depth comparison of CORSIKA 8 and CORSIKA 7 for optimum simulation settings and discuss the relevance of observed differences in light of reconstruction efforts for the energy and mass of cosmic rays.
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.
Galactic PeV cosmic-ray accelerators (PeVatrons) are Galactic sources theorized to accelerate cosmic rays up to PeV in energy. The accelerated cosmic rays are expected to interact hadronically with nearby ambient gas or the interstellar medium, resulting in γ -rays and neutrinos. Recently, the Large High Altitude Air Shower Observatory (LHAASO) identified 12 γ -ray sources with emissions above 100 TeV, making them candidates for PeVatrons. While at these high energies the Klein–Nishina effect exponentially suppresses leptonic emission from Galactic sources, evidence for neutrino emission would unequivocally confirm hadronic acceleration. Here, we present the results of a search for neutrinos from these γ -ray sources and stacking searches testing for excess neutrino emission from all 12 sources as well as their subcatalogs of supernova remnants and pulsar wind nebulae with 11 yr of track events from the IceCube Neutrino Observatory. No significant emissions were found. Based on the resulting limits, we place constraints on the fraction of γ -ray flux originating from the hadronic processes in the Crab Nebula and LHAASO J2226+6057.
For the in-ice component of the next generation neutrino observatory at the South Pole, IceCube-Gen2, a new sensor module is being developed, which is an evolution of the D-Egg and mDOM sensors developed for the IceCube Upgrade. The sensor design features up to 18 4-inch PMTs distributed homogeneously in a borosilicate glass pressure vessel. Challenges arise for the mechanical design from the tight constraints on the bore hole diameter (which will be 2 inches smaller than for IceCube Upgrade) and from the close packing of the PMTs. The electronics design must meet the space constraints posed by the mechanical design as well as the power consumption and cost considerations from over 10,000 optical modules being deployed. This contribution presents forward-looking solutions to these design considerations. Prototype modules will be installed and integrated in the IceCube Upgrade.
The IceCube neutrino observatory--installed in the Antarctic ice--is the largest neutrino telescope to date. It consists of 5,160 photomultiplier-tubes spread among 86 vertical strings making a total detector volume of more than a cubic kilometer. IceCube detects neutrinos via Cherenkov light emitted by charged relativistic particles produced when a neutrino interacts in or near the detector. The detector is particularly sensitive to high-energy neutrinos of due to its size and photosensor spacing. In this analysis we search for dark matter that annihilates into a metastable mediator that subsequently decays into Standard Model particles. These models yield an enhanced high-energy neutrino flux from dark matter annihilation inside the Sun compared to models without a mediator. Neutrino signals that are produced directly inside the Sun are strongly attenuated at higher energies due to interactions with the solar plasma. In the models considered here, the mediator can escape the Sun before producing any neutrinos, thereby avoiding attenuation. We present the results of an analysis of six years of IceCube data looking for dark matter in the Sun. We consider mediator lifetimes between 1 ms to 10 s and dark matter masses between 200 GeV and 75 TeV.
The IceCube Neutrino Observatory at the South Pole detects Cherenkov light emitted by charged secondary particles created by primary neutrino interactions. Double pulse waveforms can arise from charged current interactions of astrophysical tau neutrinos with nucleons in the ice and the subsequent decay of tau leptons. The previous 8-year tau double pulse analysis found three tau neutrino candidate events. Among them, the most promising one observed in 2014 is located very near the dust layer in the middle of the detector. A posterior analysis on this event will be presented in this paper, using a new ice model treatment with continuously varying nuisance parameters to do the targeted Monte Carlo re-simulation for tau and other background neutrino ensembles. The impact of different ice models on the expected signal and background statistics will also be discussed.
We present performance studies of a segmented optical module for the IceCube-Gen2 detector. Based on the experience gained in sensor development for the IceCube Upgrade, the new sensor will consist of up to eighteen 4 inch PMTs housed in a transparent pressure vessel, providing homogeneous angular coverage. The use of custom molded optical `gel pads' around the PMTs enhances the photon capture rate via total internal reflection at the gel-air interface. This contribution presents simulation studies of various sensor, PMT, and gel pad geometries aimed at optimizing the sensitivity of the optical module in the face of confined space and harsh environmental conditions at the South Pole.
The IceCube Neutrino Observatory is a cubic-kilometer scale neutrino detector embedded in the Antarctic ice of the South Pole. In the near future, the detector will be augmented by extensions, such as the IceCube Upgrade and the planned Gen2 detector. The sparseness of observed light in the detector for low-energy events, and the irregular detector geometry, have always been a challenge to the reconstruction of the detected neutrinos' parameters of interest. This challenge remains with the IceCube Upgrade, currently under construction, which introduces seven new detector strings with novel detector modules. The Upgrade modules will increase the detection rate of low-energy events and allow us to further constrain neutrino oscillation physics. However, the geometry of these modules render existing traditional reconstruction algorithms more difficult to use. We introduce a new reconstruction algorithm based on Graph Neural Networks, which we use to reconstruct neutrino events at much faster processing times than the traditional algorithms, while providing comparable resolution. We show that our algorithm is applicable not only to reconstructing data of the current IceCube detector, but also simulated events for next-generation extensions, such as the IceCube Upgrade.
The field of deep learning has become increasingly important for particle physics experiments, yielding a multitude of advances, predominantly in event classification and reconstruction tasks. Many of these applications have been adopted from other domains. However, data in the field of physics are unique in the context of machine learning, insofar as their generation process and the laws and symmetries they abide by are usually well understood. Most commonly used deep learning architectures fail at utilizing this available information. In contrast, more traditional likelihood-based methods are capable of exploiting domain knowledge, but they are often limited by computational complexity. In this contribution, a hybrid approach is presented that utilizes generative neural networks to approximate the likelihood, which may then be used in a traditional maximum-likelihood setting. Domain knowledge, such as invariances and detector characteristics, can easily be incorporated in this approach. The hybrid approach is illustrated by the example of event reconstruction in IceCube.
Sources of astrophysical neutrinos can potentially be discovered through the detection of neutrinos in coincidence with electromagnetic or gravitational waves. Real-time alerts generated by IceCube play an important role in this search, acting as triggers for follow-up observations with instruments sensitive to other wavelengths. Once a high-energy event is detected by the IceCube real-time program, a complex and time consuming direction reconstruction method is run in order to calculate an accurate localisation. To investigate the effect of systematic uncertainties on the uncertainty estimate of the location, we simulate a set of high-energy events with a wide range of directions for different ice model realisations, the dominant systematic error in our localization uncertainty. This makes use of a novel simulation tool, which allows the treatment of systematic uncertainties with multiple continuously varied nuisance parameters. These events will be reconstructed using various reconstruction methods. This study will enable us to include systematic uncertainties in a robust manner in the real-time direction and error estimates.
The IceCube Neutrino Observatory instruments about 1 km^3 of deep, glacial ice at the geographic South Pole using 5160 photomultipliers to detect Cherenkov light of charged relativistic particles. Most of IceCube's science goals rely heavily on an ever more precise understanding of the optical properties of the instrumented ice. A curious light propagation effect observed by the experiment is an anisotropic attenuation, which is aligned with the local flow of the ice. Having recently identified curved photon trajectories resulting from asymmetric light diffusion in the birefringent polycrystalline microstructure of the ice as the most likely underlying cause of this effect, work is now ongoing to optimize the model parameters (effectively deducing the average crystal size and shape in the detector). We present the parametrization of the birefringence effect in our photon propagation simulation, the fitting procedures and results as well as the impact of the new ice model on data-MC agreement.
Combining observational data from multiple instruments for multi-messenger astronomy can be challenging due to the complexity of the instrument response functions and likelihood calculation. We introduce a python-based unbinned-likelihood analysis package called i3mla (IceCube Maximum Likelihood Analysis). i3mla is designed to be compatible with the Multi-Mission Maximum Likelihood (3ML) framework, which enables multi-messenger astronomy analyses by combining the likelihood across different instruments. By making it possible to use IceCube data in the 3ML framework, we aim to facilitate the use of neutrino data in multi-messenger astronomy. In this work we illustrate how to use the i3mla package with 3ML and present preliminary sensitivities using the i3mla package and 3ML through a joint-fit with HAWC Public dataset.
We summarize initial results for high-energy neutrino counterpart searches coinciding with gravitational-wave events in LIGO/Virgo's GWTC-2 catalog using IceCube's neutrino triggers. We did not find any statistically significant high-energy neutrino counterpart and derived upper limits on the time-integrated neutrino emission on Earth as well as the isotropic equivalent energy emitted in high-energy neutrinos for each event.
The Sun provides an excellent target for studying spin-dependent dark matter-proton scattering due to its high matter density and abundant hydrogen content. Dark matter particles from the Galactic halo can elastically interact with Solar nuclei, resulting in their capture and thermalization in the Sun. The captured dark matter can annihilate into Standard Model particles including an observable flux of neutrinos. We present the results of a search for low-energy ($<$ 500 GeV) neutrinos correlated with the direction of the Sun using 7 years of IceCube data. This work utilizes, for the first time, new optimized cuts to extend IceCube's sensitivity to dark matter mass down to 5 GeV. We find no significant detection of neutrinos from the Sun. Our observations exclude capture by spin-dependent dark matter-proton scattering with cross-section down to a few times $10^{-41}$ cm$^2$, assuming there is equilibrium with annihilation into neutrinos/anti-neutrinos for dark matter masses between 5 GeV and 100 GeV. These are the strongest constraints at GeV energies for dark matter annihilation directly to neutrinos.
Galaxy clusters have the potential to accelerate cosmic rays (CRs) to ultra-high energies via accretion shocks or embedded CR acceleration sites. CRs with energies below the Hillas condition will be confined within the cluster and will eventually interact with the intracluster medium (ICM) gas to produce secondary neutrinos and $\gamma$ rays. Using 9.5 years of muon-neutrino track events from the IceCube Neutrino Observatory, we report the results of a stacking analysis of 1094 galaxy clusters, with masses $\gtrsim 10^{14}$ \(\textup{M}_\odot\) and redshifts between 0.01 and $\sim$1, detected by the {\it Planck} mission via the Sunyaev-Zeldovich (SZ) effect. We find no evidence for significant neutrino emission and report upper limits on the cumulative unresolved neutrino flux from massive galaxy clusters after accounting for the completeness of the catalog up to a redshift of 2, assuming three different weighting scenarios for the stacking and three different power-law spectra. Weighting the sources according to mass and distance, we set upper limits at $90\%$ confidence level that constrain the flux of neutrinos from massive galaxy clusters ($\gtrsim 10^{14}$ \(\textup{M}_\odot\)) to be no more than $4.6\%$ of the diffuse IceCube observations at 100~TeV, assuming an unbroken $E^{-2.5}$ power-law spectrum.
DeepCore, as a densely instrumented sub-detector of IceCube, extends IceCube's energy reach down to about 10 GeV, enabling the search for astrophysical transient sources, e.g., choked gamma-ray bursts. While many other past and on-going studies focus on triggered time-dependent analyses, we aim to utilize a newly developed event selection and dataset for an untriggered all-sky timedependent search for transients. In this work, all-flavor neutrinos are used, where neutrino types are determined based on the topology of the events. We extend the previous DeepCore transient half-sky search to an all-sky search and focus only on short timescale sources (with a duration of 10(2) similar to 10(5) seconds). All-sky sensitivities to transients in an energy range from 10 GeV to 300 GeV will be presented in this poster. We show that DeepCore can be reliably used for all-sky searches for short-lived astrophysical sources.
We present a measurement of the high-energy astrophysical muon – neutrino fl ux with the IceCube Neutrino Observatory. The measurement uses a high-purity selection of 650k neutrino-induced muon tracks from the northern celestial hemisphere, corresponding to 9.5 yr of experimental data. With respect to previous publications, the measurement is improved by the increased size of the event sample and the extended model testing beyond simple power-law hypotheses. An updated treatment of systematic uncertainties and atmospheric background fl uxes has been implemented based on recent models. The best- fi t single power-law parameterization for the astrophysical energy spectrum results in a normalization of in the energy range from 15 TeV to PeV. The model tests include a single power law with a spectral cutoff at high energies, a log-parabola model, several source-class-speci ux predictions from the literature, and a model-independent spectral unfolding. The data are consistent with a single power-law hypothesis, however, spectra with softening above one PeV are statistically more favorable at a two-sigma level.
The IceCube Neutrino Observatory first observed a diffuse flux of high energy astrophysical neutrinos in 2013. Since then, this observation has been confirmed in multiple detection channels such as high energy starting events, cascades, and through-going muon tracks. Combining these event selections into a high statistics global fit of 10 years of IceCube's neutrino data could strongly improve the understanding of the diffuse astrophysical neutrino flux: challenging or confirming the simple unbroken power-law flux model as well as the astrophysical neutrino flux composition. One key component of such a combined analysis is the consistent modelling of systematic uncertainties of different event selections. This can be achieved using the novel SnowStorm Monte Carlo method which allows constraints to be placed on multiple systematic parameters from a single simulation set. We will report on the status of a new combined analysis of through-going muon tracks and cascades. It is based on a consistent all flavor neutrino signal and background simulation using, for the first time, the SnowStorm method to analyze IceCube's high-energy neutrino data. Estimated sensitivities for the energy spectrum of the diffuse astrophysical neutrino flux will be shown.
The first full angular analysis and an updated measurement of the decay-rate CP asymmetry of the D0 ??? x-+x-?????+????? and D0 ??? K+K?????+????? decays are reported. The analysis uses proton-proton collision data collected with the LHCb detector at center-of-mass energies of 7, 8, and 13 TeV. The dataset corresponds to an integrated luminosity of 9 fb???1. The full set of CP -averaged angular observables and their CP asymmetries are measured as a function of the dimuon invariant mass. The results are consistent with expectations from the standard model and with CP symmetry.
The IceCube DeepCore is a dense infill array of the IceCube Neutrino Observatory at the South Pole. While IceCube is best suited for detecting neutrinos with energies of several 100 GeV and above, DeepCore allows to probe neutrinos with lower energies. We focus on a sample of neutrinos with energies above approximately 10 GeV, which was originally optimised for oscillation experiments. Recently, it has been adapted to enable searches for transient sources of astrophysical neutrinos in the sky. In particular, this low-energy dataset can be used to conduct follow-up searches of gravitational wave transients detected by the LIGO-Virgo instruments. A study of this, which complements IceCube's follow-up of gravitational wave events using high-energy neutrino samples, will be discussed here.