Many Galactic sources of gamma rays, such as supernova remnants, are expected to produce neutrinos with a typical energy cutoff well below 100 TeV. For the IceCube Neutrino Observatory located at the South Pole, the southern sky, containing the inner part of the Galactic plane and the Galactic Center, is a particularly challenging region at these energies, because of the large background of atmospheric muons. In this paper, we present recent advancements in data selection strategies for track-like muon neutrino events with energies below 100 TeV from the southern sky. The strategies utilize the outer detector regions as veto and features of the signal pattern to reduce the background of atmospheric muons to a level which, for the first time, allows IceCube searching for point-like sources of neutrinos in the southern sky at energies between 100 GeV and several TeV in the muon neutrino charged current channel. No significant clustering of neutrinos above background expectation was observed in four years of data recorded with the completed IceCube detector. Upper limits on the neutrino flux for a number of spectral hypotheses are reported for a list of astrophysical objects in the southern hemisphere.
In order to identify the sources of the observed diffuse high-energy neutrino flux, it is crucial to discover their electromagnetic counterparts. To increase the sensitivity of detecting counterparts of transient or variable sources by telescopes with a limited field of view, IceCube began releasing alerts for single high-energy (E-v > 60 TeV) neutrino detections with sky localisation regions of order 1 degrees radius in 2016. We used Pan-STARRS1 to follow-up five of these alerts during 2016-2017 to search for any optical transients that may be related to the neutrinos. Typically 10-20 faint m(ip1) less than or similar to 22.5 mag) extragalactic transients are found within the Pan-STARRS1 footprints and are generally consistent with being unrelated field supernovae (SNe) and AGN. We looked for unusual properties of the detected transients, such as temporal coincidence of explosion epoch with the IceCube timestamp, or other peculiar light curve and physical properties. We found only one transient that had properties worthy of a specific follow-up. In the Pan-STARRS1 imaging for IceCube-160427A (probability to be of astrophysical origin of similar to 50%), we found a SN PS16cgx, located at 10.0' from the nominal IceCube direction. Spectroscopic observations of PS16cgx showed that it was an H-poor SN at redshift z = 0.2895 +/- 0.0001. The spectra and light curve resemble some high-energy Type Ic SNe, raising the possibility of a jet driven SN with an explosion epoch temporally coincident with the neutrino detection. However, distinguishing Type Ia and Type Ic SNe at this redshift is notoriously difficult. Based on all available data we conclude that the transient is more likely to be a Type Ia with relatively weak Sin absorption and a fairly normal rest-frame r-band light curve. If, as predicted, there is no high-energy neutrino emission from Type Ia SNe, then PS16cgx must be a random coincidence, and unrelated to the IceCube-160427A. We find no other plausible optical transient for any of the five IceCube events observed down to a 5 sigma limiting magnitude of mip1 approximate to 22 mag, between 1 day and 25 days after detection.
In the analysis published in 1, constraints on the number of signal events can be interpreted as constraints on the volumetric neutrino to muon conversion rate.
Gamma-ray blazars have been suggested as promising contributors to the diffuse flux of high-energy neutrinos detected by IceCube. After the identification of the gamma-ray blazar TXS 0506+056 as the first compelling IceCube neutrino source candidate, a systematic analysis of all high-energy neutrino events satisfying the IceCube realtime trigger criteria was performed. One additional known gamma-ray source, the blazar GB6 J1040+0617, was found to be in spatial coincidence with a neutrino in this sample. The chance probability of this coincidence is 30$%$ after trials correction. For the first time, we present a systematic study of the gamma-ray flux, spectral and optical variability and multi-wavelength behavior of GB6 J1040+0617 and compare it to TXS 0506+056. We discuss a possible source confusion with a new emerging gamma-ray source, the blazar 4C+06.41. We find that TXS 0506+056 shows strong flux variability in the $textit{Fermi}$-LAT gamma-ray band, being in an active state around the arrival of IceCube-170922A, but in a low state during the archival IceCube neutrino flare in 2014/15. In both cases the spectral shape is statistically compatible ($leq 2sigma$) with the average spectrum. Assuming the reported redshift estimate of z=0.73, GB6 J1040+0617 has an average gamma-ray luminosity comparable with that of TXS 0506+056. It shows a bright optical flare recorded by ASAS-SN accompanied by modest gamma-ray activity at the neutrino arrival time. While the association with the neutrino is consistent with background expectations, GB6 J1040+0617 appears to be a plausible neutrino source candidate based on its energetics and multi-wavelength features. Finding one or two neutrinos originating from gamma-ray blazars in the given sample of high-energy neutrinos is consistent with previously derived limits of neutrino emission from gamma-ray blazars.
We present a measurement of atmospheric tau neutrino appearance from oscillations with three years of data from the DeepCore sub-array of the IceCube Neutrino Observatory. This analysis uses atmospheric neutrinos from the full sky with reconstructed energies between 5.6 GeV and 56 GeV to search for a statistical excess of cascade-like neutrino events which are the signature of nutau interactions. For CC+NC (CC-only) interactions, we measure the tau neutrino normalization to be 0.73 +0.30 -0.24 (0.57 +0.36 -0.30) and exclude the absence of tau neutrino oscillations at a significance of 3.2 sigma (2.0 sigma) These results are consistent with, and of similar precision to, a confirmatory IceCube analysis also presented, as well as measurements performed by other experiments.
We present the first full-sky analysis of the cosmic ray arrival direction distribution with data collected by the High-Altitude Water Cherenkov and IceCube observatories in the northern and southern hemispheres at the same median primary particle energy of 10 TeV. The combined sky map and angular power spectrum largely eliminate biases that result from partial sky coverage and present a key to probe into the propagation properties of TeV cosmic rays through our local interstellar medium and the interaction between the interstellar and heliospheric magnetic fields. From the map, we determine the horizontal dipole components of the anisotropy δ 0 h = 9.16 × 10 −4 and δ 6 h = 7.25 × 10 −4 (±0.04 × 10 −4 ). In addition, we infer the direction (229.°2 ± 3.°5 R.A., 11.°4 ± 3.°0 decl.) of the interstellar magnetic field from the boundary between large-scale excess and deficit regions from which we estimate the missing corresponding vertical dipole component of the large-scale anisotropy to be
High-energy neutrino emission has been predicted for several short-lived astrophysical transients including gamma-ray bursts (GRBs), core-collapse supernovae with choked jets, and neutron star mergers. IceCube's optical and x-ray follow-up program searches for such transient sources by looking for two or more muon neutrino candidates in directional coincidence and arriving within 100 s. The measured rate of neutrino alerts is consistent with the expected rate of chance coincidences of atmospheric background events and no likely electromagnetic counterparts have been identified in Swift follow-up observations. Here, we calculate generic bounds on the neutrino flux of short-lived transient sources. Assuming an E^{-2.5} neutrino spectrum, we find that the neutrino flux of rare sources, like long gamma-ray bursts, is constrained to <5% of the detected astrophysical flux and the energy released in neutrinos (100 GeV to 10 PeV) by a median bright GRB-like source is <10^{52.5} erg. For a harder E^{-2.13} neutrino spectrum up to 30% of the flux could be produced by GRBs and the allowed median source energy is <10^{52} erg. A hypothetical population of transient sources has to be more common than 10^{-5} Mpc^{-3} yr^{-1} (5×10^{-8} Mpc^{-3} yr^{-1} for the E^{-2.13} spectrum) to account for the complete astrophysical neutrino flux.
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 E^2 d N /d E = 3· 10^-13 TeV cm^-2 s^-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.
Inelasticity, the fraction of a neutrino's energy transferred to hadrons, is a quantity of interest in the study of astrophysical and atmospheric neutrino interactions at multi-TeV energies with IceCube. In this work, a sample of contained neutrino interactions in IceCube is obtained from five years of data and classified as 2650 tracks and 965 cascades. Tracks arise predominantly from charged-current nu(mu) interactions, and we demonstrate that we can reconstruct their energy and inelasticity. The inelasticity distribution is found to be consistent with the calculation of Cooper-Sarkar et al. across the energy range from similar to 1 to similar to 100 TeV. Along with cascades from neutrinos of all flavors, we also perform a fit over the energy, zenith angle, and inelasticity distribution to characterize the flux of astrophysical and atmospheric neutrinos. The energy spectrum of diffuse astrophysical neutrinos is described well by a power law in both track and cascade samples, and a best-fit index gamma = 2.62 +/- 0.07 is found in the energy range from 3.5 TeV to 2.6 PeV. Limits are set on the astrophysical flavor composition and are compatible with a ratio of (1/3 : 1/3 : 1/3)(circle plus). Exploiting the distinct inelasticity distribution of nu(mu) and (nu) over bar (mu) interactions, the atmospheric nu(mu) to (nu) over bar (mu) flux ratio in the energy range from 770 GeV to 21 TeV is found to be 0.77(-0.25)(+0.44) times the calculation by Honda et al. Lastly, the inelasticity distribution is also sensitive to neutrino charged-current charm production. The data are consistent with a leading-order calculation, with zero charm production excluded at 91% confidence level. Future analyses of inelasticity distributions may probe new physics that affects neutrino interactions both in and beyond the Standard Model.
We report on the observation of a deficit in the cosmic ray flux from the directions of the Moon and Sun with five years of data taken by the IceCube Neutrino Observatory. Between 2010 May and 2011 May the IceCube detector operated with 79 strings deployed in the glacial ice at the South Pole, and with 86 strings between 2011 May and 2015 May. A binned analysis is used to measure the relative deficit and significance of the cosmic ray shadows. Both the cosmic ray Moon and Sun shadows are detected with high statistical significance (>10σ) for each year. The results for the Moon shadow are consistent with previous analyses and verify the stability of the IceCube detector over time. This work represents the first observation of the Sun shadow with the IceCube detector. We show that the cosmic ray shadow of the Sun varies with time. These results make it possible to study cosmic ray transport near the Sun with future data from IceCube.
We present the first full-sky analysis of the cosmic ray arrival direction distribution with data collected by the High-Altitude Water Cherenkov and IceCube observatories in the northern and southern hemispheres at the same median primary particle energy of 10 TeV. The combined sky map and angular power spectrum largely eliminate biases that result from partial sky coverage and present a key to probe into the propagation properties of TeV cosmic rays through our local interstellar medium and the interaction between the interstellar and heliospheric magnetic fields. From the map we determine the horizontal dipole components of the anisotropy δ0h = 9.16×10−4 and δ6h = 7.25×10−4 (±0.04×10−4). In addition, we infer the direction (229.2±3◦.5 R.A. , 11.4±3◦.0 decl.) of the interstellar magnetic field from the boundary between large scale excess and deficit regions from which we estimate the missing corresponding vertical dipole component of the large scale anisotropy to be δN ∼ −3.97 −2.0 × 10−4.
In the analyses, published in Ref. [1], the exclusion limits are calculated in dependence of the mean free path of the magnetic monopole - nucleon catalysis interaction.
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 similar to 10(2) s from the northern hemisphere stacking search.
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.
The existence of diffuse Galactic neutrino production is expected from cosmic-ray interactions with Galactic gas and radiation fields. Thus, neutrinos are a unique messenger offering the opportunity to test the products of Galactic cosmic-ray interactions up to energies of hundreds of TeV. Here we present a search for this production using ten years of Astronomy with a Neutrino Telescope and Abyss environmental RESearch (ANTARES) track and shower data, as well as seven years of IceCube track data. The data are combined into a joint likelihood test for neutrino emission according to the KRA(gamma) model assuming a 5 PeV per nucleon Galactic cosmic-ray cutoff. No significant excess is found. As a consequence, the limits presented in this Letter start constraining the model parameter space for Galactic cosmic-ray production and transport.
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 epsilon(mu tau) = -0.0005, with a 90% C.L. allowed range of -0.0067 < epsilon(mu tau) < 0.0081. This result is more restrictive than recent limits from other experiments for.mu t. Furthermore, our result is complementary to a recent constraint on epsilon(mu tau) using another publicly available IceCube high-energy event selection. Together, they constitute the world's best limits on nonstandard interactions in the mu - tau sector.
The core mission of the IceCube neutrino observatory is to study the origin and propagation of cosmic rays. IceCube, with its surface component IceTop, observes multiple signatures to accomplish this mission. Most important are the astrophysical neutrinos that are produced in interactions of cosmic rays, close to their sources and in interstellar space. IceCube is the first instrument that measures the properties of this astrophysical neutrino flux and constrains its origin. In addition, the spectrum, composition, and anisotropy of the local cosmic-ray flux are obtained from measurements of atmospheric muons and showers. Here we provide an overview of recent findings from the analysis of IceCube data, and their implications to our understanding of cosmic rays. (C) 2017 COSPAR. Published by Elsevier Ltd. All rights reserved.
With the observation of high-energy astrophysical neutrinos by the IceCube Neutrino Observatory, interest has risen in models of PeV-mass decaying dark matter particles to explain the observed flux. We present two dedicated experimental analyses to test this hypothesis. One analysis uses 6 years of IceCube data focusing on muon neutrino track' events from the Northern Hemisphere, while the second analysis uses 2 years of cascade' events from the full sky. Known background components and the hypothetical flux from unstable dark matter are fitted to the experimental data. Since no significant excess is observed in either analysis, lower limits on the lifetime of dark matter particles are derived: we obtain the strongest constraint to date, excluding lifetimes shorter than 10(28) s at 90% CL for dark matter masses above 10 TeV.
We report the quasi-differential upper limits of extremely high energy (EHE) neutrino flux above 10 PeV based on the analysis of nine years of IceCube data. A complete frequentist approach to calculate the differential limit using the Poisson binned likelihood is developed. It enables the limit to be set in the presence of unknown astrophysical neutrino flux. An event with deposited energy clearly above 1 PeV was detected in addition to two events found in the previous EHE neutrino search. They are consistent with the astrophysical neutrino flux of a power-law-like spectrum but incompatible with predictions of cosmogenic neutrino fluxes with spectrum peaks at energies well above the PeV range. Thus, they are considered as the bulk of background events in setting the limits on EHE neutrino fluxes. The resultant differential upper limit is the most stringent to date in the energy range between $5 \times 10^{6}$ and $5 \times 10^{10}$ GeV. This result indicates that cosmogenic neutrino models that predict a three-flavor neutrino flux of $E_\nu^2\phi_{\nu_e+\nu_\mu+\nu_\tau}\simeq2\times 10^{-8}\ {\rm GeV}/{\rm cm}^2\ \sec\ {\rm sr}$ at $10^9\ {\rm GeV}$ are constrained, bounding a significant parameter space on EHE neutrino models, which assumes a composition of proton-dominated ultra-high-energy cosmic rays.
Individual astrophysical sources previously detected in neutrinos are limited to the Sun and the supernova 1987A, whereas the origins of the diffuse flux of high-energy cosmic neutrinos remain unidentified. On 22 September 2017 we detected a high-energy neutrino, IceCube-170922A, with an energy of approximately 290 TeV. Its arrival direction was consistent with the location of a known gamma-ray blazar TXS 0506+056, observed to be in a flaring state. An extensive multi-wavelength campaign followed, ranging from radio frequencies to gamma-rays. These observations characterize the variability and energetics of the blazar and include the first detection of TXS 0506+056 in very-high-energy gamma-rays. This observation of a neutrino in spatial coincidence with a gamma-ray emitting blazar during an active phase suggests that blazars may be a source of high-energy neutrinos.