NuRadioMC is a Monte Carlo framework designed to simulate ultra-high energy neutrino detectors that rely on the radio detection method. This method exploits the radio emission generated in the electromagnetic component of a particle shower following a neutrino interaction. NuRadioMC simulates everything from the neutrino interaction in a medium, the subsequent Askaryan radio emission, the propagation of the radio signal to the detector and finally the detector response. NuRadioMC is designed as a modern, modular Python-based framework, combining flexibility in detector design with user-friendliness. It includes a state-of-the-art event generator, an improved modelling of the radio emission, a revisited approach to signal propagation and increased flexibility and precision in the detector simulation. This paper focuses on the implemented physics processes and their implications for detector design. A variety of models and parameterizations for the radio emission of neutrino-induced showers are compared and reviewed. Comprehensive examples are used to discuss the capabilities of the code and different aspects of instrumental design decisions.

We present a search for a neutrino signal from dark matter self-annihilations in the Milky Way using the IceCube Neutrino Observatory (IceCube). In 1005 days of data we found no significant excess of neutrinos over the background of neutrinos produced in atmospheric air showers from cosmic ray interactions. We derive upper limits on the velocity averaged product of the dark matter self-annihilation cross section and the relative velocity of the dark matter particles ⟨ σAv⟩. Upper limits are set for dark matter particle candidate masses ranging from 10 GeV up to 1 TeV while considering annihilation through multiple channels. This work sets the most stringent limit on a neutrino signal from dark matter with mass between 10 and 100 GeV, with a limit of 1.18·10-23cm3s-1 for 100 GeV dark matter particles self-annihilating via τ+τ- to neutrinos (assuming the Navarro–Frenk–White dark matter halo profile).

We present a search for coincidence between IceCube TeV neutrinos and fast radio bursts (FRBs). During the search period from 2010 May 31 to 2016 May 12, a total of 29 FRBs with 13 unique locations have been detected in the whole sky. An unbinned maximum likelihood method was used to search for spatial and temporal coincidence between neutrinos and FRBs in expanding time windows, in both the northern and southern hemispheres. No significant correlation was found in six years of IceCube data. Therefore, we set upper limits on neutrino fluence emitted by FRBs as a function of time window duration. We set the most stringent limit obtained to date on neutrino fluence from FRBs with an E -2 energy spectrum assumed, which is 0.0021 GeV cm-2 per burst for emission timescales up to ∼102 s from the northern hemisphere stacking search.

The IceCube Collaboration has observed a high-energy astrophysical neutrino flux and recently found evidence for neutrino emission from the blazar TXS 0506+ 056. These results open a new window into the high-energy universe. However, the source or sources of most of the observed flux of astrophysical neutrinos remains uncertain. Here, a search for steady point-like neutrino sources is performed using an unbinned likelihood analysis. The method searches for a spatial accumulation of muon-neutrino events using the very high-statistics sample of about 497,000 neutrinos recorded by IceCube between 2009 and 2017. The median angular resolution is ∼ 1 ∘ at 1 TeV and improves to ∼ 0. 3 ∘ for neutrinos with an energy of 1 PeV. Compared to previous analyses, this search is optimized for point-like neutrino emission with the same flux-characteristics as the observed astrophysical muon-neutrino flux and introduces an improved event-reconstruction and parametrization of the background. The result is an improvement in sensitivity to the muon-neutrino flux compared to the previous analysis of ∼ 35 % assuming an E- 2 spectrum. The sensitivity on the muon-neutrino flux is at a level of E2dN/dE=3·10-13TeVcm-2s-1. No new evidence for neutrino sources is found in a full sky scan and in an a priori candidate source list that is motivated by gamma-ray observations. Furthermore, no significant excesses above background are found from populations of sub-threshold sources. The implications of the non-observation for potential source classes are discussed.

The origins of high-energy astrophysical neutrinos remain a mystery despite extensive searches for their sources. We present constraints from seven years of IceCube Neutrino Observatory muon data on the neutrino flux coming from the Galactic plane. This flux is expected from cosmic-ray interactions with the interstellar medium or near localized sources. Two methods were developed to test for a spatially extended flux from the entire plane, both of which are maximum likelihood fits but with different signal and background modeling techniques. We consider three templates for Galactic neutrino emission based primarily on gamma-ray observations and models that cover a wide range of possibilities. Based on these templates and in the benchmark case of an unbroken E-2.5 power-law energy spectrum, we set 90% confidence level upper limits, constraining the possible Galactic contribution to the diffuse neutrino flux to be relatively small, less than 14% of the flux reported in Aartsen et al. above 1 TeV. A stacking method is also used to test catalogs of known high-energy Galactic gamma-ray sources.

Neutrinos interact only very weakly, so they are extremely penetrating. The theoretical neutrino-nucleon interaction cross-section, however, increases with increasing neutrino energy, and neutrinos with energies above 40 teraelectronvolts (TeV) are expected to be absorbed as they pass through the Earth. Experimentally, the cross-section has been determined only at the relatively low energies (below 0.4 TeV) that are available at neutrino beams from accelerators. Here we report a measurement of neutrino absorption by the Earth using a sample of 10,784 energetic upward-going neutrino-induced muons. The flux of high-energy neutrinos transiting long paths through the Earth is attenuated compared to a reference sample that follows shorter trajectories. Using a fit to the two-dimensional distribution of muon energy and zenith angle, we determine the neutrino-nucleon interaction cross-section for neutrino energies 6.3-980 TeV, more than an order of magnitude higher than previous measurements. The measured cross-section is about 1.3 times the prediction of the standard model, consistent with the expectations for charged- and neutral-current interactions. We do not observe a large increase in the cross-section with neutrino energy, in contrast with the predictions of some theoretical models, including those invoking more compact spatial dimensions or the production of leptoquarks. This cross-section measurement can be used to set limits on the existence of some hypothesized beyond-standard-model particles, including leptoquarks.

As atmospheric neutrinos propagate through the Earth, vacuumlike oscillations are modified by Standard Model neutral- and charged-current interactions with electrons. Theories beyond the Standard Model introduce heavy, TeV-scale bosons that can produce nonstandard neutrino interactions. These additional interactions may modify the Standard Model matter effect producing a measurable deviation from the prediction for atmospheric neutrino oscillations. The result described in this paper constrains nonstandard interaction parameters, building upon a previous analysis of atmospheric muon-neutrino disappearance with three years of IceCube DeepCore data. The best fit for the muon to tau flavor changing term is ϵμτ=-0.0005, with a 90% C.L. allowed range of -0.0067<ϵμτ<0.0081. This result is more restrictive than recent limits from other experiments for ϵμτ. Furthermore, our result is complementary to a recent constraint on ϵμτ using another publicly available IceCube high-energy event selection. Together, they constitute the world's best limits on nonstandard interactions in the μ-τ sector.

The IceCube neutrino observatory has established the existence of a flux of high-energy astrophysical neutrinos, which is inconsistent with the expectation from atmospheric backgrounds at a significance greater than 5σ. This flux has been observed in analyses of both track events from muon neutrino interactions and cascade events from interactions of all neutrino flavors. Searches for astrophysical neutrino sources have focused on track events due to the significantly better angular resolution of track reconstructions. To date, no such sources have been confirmed. Here we present the first search for astrophysical neutrino sources using cascades interacting in IceCube with deposited energies as small as 1 TeV. No significant clustering was observed in a selection of 263 cascades collected from 2010 May to 2012 May. We show that compared to the classic approach using tracks, this statistically independent search offers improved sensitivity to sources in the southern sky, especially if the emission is spatially extended or follows a soft energy spectrum. This enhancement is due to the low background from atmospheric neutrinos forming cascade events and the additional veto of atmospheric neutrinos at declinations ≲-30.

We present a measurement of the atmospheric neutrino oscillation parameters using three years of data from the IceCube Neutrino Observatory. The DeepCore infill array in the center of IceCube enables the detection and reconstruction of neutrinos produced by the interaction of cosmic rays in Earth's atmosphere at energies as low as ∼5  GeV. That energy threshold permits measurements of muon neutrino disappearance, over a range of baselines up to the diameter of the Earth, probing the same range of L/E_{ν} as long-baseline experiments but with substantially higher-energy neutrinos. This analysis uses neutrinos from the full sky with reconstructed energies from 5.6 to 56 GeV. We measure Δm_{32}^{2}=2.31_{-0.13}^{+0.11}×10^{-3}  eV^{2} and sin^{2}θ_{23}=0.51_{-0.09}^{+0.07}, assuming normal neutrino mass ordering. These results are consistent with, and of similar precision to, those from accelerator- and reactor-based experiments.

Lorentz symmetry is a fundamental spacetime symmetry underlying both the standard model of particle physics and general relativity. This symmetry guarantees that physical phenomena are observed to be the same by all inertial observers. However, unified theories, such as string theory, allow for violation of this symmetry by inducing new spacetime structure at the quantum gravity scale. Thus, the discovery of Lorentz symmetry violation could be the first hint of these theories in nature. Here we report the results of the most precise test of spacetime symmetry in the neutrino sector to date. We use high-energy atmospheric neutrinos observed at the IceCube Neutrino Observatory to search for anomalous neutrino oscillations as signals of Lorentz violation. We find no evidence for such phenomena. This allows us to constrain the size of the dimension-four operator in the standard-model extension for Lorentz violation to the 1 0 - 28 level and to set limits on higher-dimensional operators in this framework. These are among the most stringent limits on Lorentz violation set by any physical experiment.