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.

Adopting the Standard Halo Model (SHM) of an isotropic Maxwellian velocity distribution for dark matter (DM) particles in the Galaxy, the most stringent current constraints on their spin-dependent scattering cross-section with nucleons come from the IceCube neutrino observatory and the PICO-60 C$_3$F$_8$ superheated bubble chamber experiments. The former is sensitive to high energy neutrinos from the self-annihilation of DM particles captured in the Sun, while the latter looks for nuclear recoil events from DM scattering off nucleons. Although slower DM particles are more likely to be captured by the Sun, the faster ones are more likely to be detected by PICO. Recent N-body simulations suggest significant deviations from the SHM for the smooth halo component of the DM, while observations hint at a dominant fraction of the local DM being in substructures. We use the method of Ferrer et al. (2015) to exploit the complementarity between the two approaches and derive conservative constraints on DM-nucleon scattering. Our results constrain $\sigma_{\mathrm{SD}} \lesssim 3 \times 10^{-39} \mathrm{cm}^2$ (6 $ \times 10^{-38} \mathrm{cm}^2$) at $\gtrsim 90\%$ C.L. for a DM particle of mass 1~TeV annihilating into $\tau^+ \tau^-$ ($b\bar{b}$) with a local density of $\rho_{\mathrm{DM}} = 0.3~\mathrm{ GeV/cm}^3$. The constraints scale inversely with $\rho_{\mathrm{DM}}$ and are independent of the DM velocity distribution.

During the first three flights of the Antarctic Impulsive Transient Antenna (ANITA) experiment, the collaboration detected several neutrino candidates. Two of these candidate events were consistent with an ultra-high-energy up-going air shower and compatible with a tau neutrino interpretation. A third neutrino candidate event was detected in a search for Askaryan radiation in the Antarctic ice, although it is also consistent with the background expectation. The inferred emergence angle of the first two events is in tension with IceCube and ANITA limits on isotropic cosmogenic neutrino fluxes. Here, we test the hypothesis that these events are astrophysical in origin, possibly caused by a point source in the reconstructed direction. Given that any ultra-high-energy tau neutrino flux traversing the Earth should be accompanied by a secondary flux in the TeV-PeV range, we search for these secondary counterparts in seven years of IceCube data using three complementary approaches. In the absence of any significant detection, we set upper limits on the neutrino flux from potential point sources. We compare these limits to ANITA's sensitivity in the same direction and show that an astrophysical explanation of these anomalous events under standard model assumptions is severely constrained regardless of source spectrum.

We report on the first measurement of the astrophysical neutrino flux using particle showers (cascades) in IceCube data from 2010 -- 2015. Assuming standard oscillations, the astrophysical neutrinos in this dedicated cascade sample are dominated ($\sim 90 \%$) by electron and tau flavors. The flux, observed in the energy range from $16\,\mathrm{TeV} $ to $2.6\,\mathrm{PeV}$, is consistent with a single power-law as expected from Fermi-type acceleration of high energy particles at astrophysical sources. We find the flux spectral index to be $\gamma=2.53\pm0.07$ and a flux normalization for each neutrino flavor of $\phi_{astro} = 1.66^{+0.25}_{-0.27}$ at $E_{0} = 100\, \mathrm{TeV}$. This flux of electron and tau neutrinos is in agreement with IceCube muon neutrino results and with all-neutrino flavor results. Results from fits assuming more complex neutrino flux models suggest a flux softening at high energies and a flux hardening at low energies (p-value $\ge 0.06$).

Cosmic-ray interactions with the solar atmosphere are expected to produce particle showers which in turn produce neutrinos from weak decays of mesons. These solar atmospheric neutrinos (SA$\nu$s) have never been observed experimentally. A detection would be an important step in understanding cosmic-ray propagation in the inner solar system and the dynamics of solar magnetic fields. SA$\nu$s also represent an irreducible background to solar dark matter searches and a detection would allow precise characterization of this background. Here, we present the first experimental search based on seven years of data collected from May 2010 to May 2017 in the austral winter with the IceCube Neutrino Observatory. An unbinned likelihood analysis is performed for events reconstructed within 5 degrees of the center of the Sun. No evidence for a SA$\nu$ flux is observed. After inclusion of systematic uncertainties, we set a 90\% upper limit of $1.02^{+0.20}_{-0.18}\cdot10^{-13}$~$\mathrm{GeV^{-1}cm^{-2}s^{-1}}$ at 1 TeV.

The current and upcoming generation of Very Large Volume Neutrino Telescopes---collecting unprecedented quantities of neutrino events---can be used to explore subtle effects in oscillation physics, such as (but not restricted to) the neutrino mass ordering. The sensitivity of an experiment to these effects can be estimated from Monte Carlo simulations. With the high number of events that will be collected, there is a trade-off between the computational expense of running such simulations and the inherent statistical uncertainty in the determined values. In such a scenario, it becomes impractical to produce and use adequately-sized sets of simulated events with traditional methods, such as Monte Carlo weighting. In this work we present a staged approach to the generation of binned event distributions in order to overcome these challenges. By combining multiple integration and smoothing techniques which address limited statistics from simulation it arrives at reliable analysis results using modest computational resources.

The distribution of galaxies within the local universe is characterized by anisotropic features. Observatories searching for the production sites of astrophysical neutrinos can take advantage of these features to establish directional correlations between a neutrino dataset and overdensities in the galaxy distribution in the sky. The results of two correlation searches between a seven-year time-integrated neutrino dataset from the IceCube Neutrino Observatory, and the 2MASS Redshift Survey (2MRS) catalog are presented here. The first analysis searches for neutrinos produced via interactions between diffuse intergalactic Ultra-High Energy Cosmic Rays (UHECRs) and the matter contained within galaxies. The second analysis searches for low-luminosity sources within the local universe, which would produce subthreshold multiplets in the IceCube dataset that directionally correlate with galaxy distribution. No significant correlations were observed in either analyses. Constraints are presented on the flux of neutrinos originating within the local universe through diffuse intergalactic UHECR interactions, as well as on the density of standard candle sources of neutrinos at low luminosities.

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.

We report on measurements of the all-particle cosmic ray energy spectrum and composition in the PeV to EeV energy range using three years of data from the IceCube Neutrino Observatory. The IceTop detector measures cosmic ray induced air showers on the surface of the ice, from which the energy spectrum of cosmic rays is determined by making additional assumptions about the mass composition. A separate measurement is performed when IceTop data are analyzed in coincidence with the high-energy muon energy loss information from the deep in-ice IceCube detector. In this measurement, both the spectrum and the mass composition of the primary cosmic rays are simultaneously reconstructed using a neural network trained on observables from both detectors. The performance and relative advantages of these two distinct analyses are discussed, including the systematic uncertainties and the dependence on the hadronic interaction models, and both all-particle spectra as well as individual spectra for elemental groups are presented.

This paper presents the results from point-like neutrino source searches using ten years of IceCube data collected between Apr.~6, 2008 and Jul.~10, 2018. We evaluate the significance of an astrophysical signal from a point-like source looking for an excess of clustered neutrino events with energies typically above $\sim1\,$TeV among the background of atmospheric muons and neutrinos. We perform a full-sky scan, a search within a selected source catalog, a catalog population study, and three stacked Galactic catalog searches. The most significant point in the Northern hemisphere from scanning the sky is coincident with the Seyfert II galaxy NGC 1068, which was included in the source catalog search. The excess at the coordinates of NGC 1068 is inconsistent with background expectations at the level of $2.9\,\sigma$ after accounting for statistical trials. The combination of this result along with excesses observed at the coordinates of three other sources, including TXS 0506+056, suggests that, collectively, correlations with sources in the Northern catalog are inconsistent with background at 3.3$\,\sigma$ significance. These results, all based on searches for a cumulative neutrino signal integrated over the ten years of available data, motivate further study of these and similar sources, including time-dependent analyses, multimessenger correlations, and the possibility of stronger evidence with coming upgrades to the detector.

The presence of a population of point sources in a dataset modifies the underlying neutrino-count statistics from the Poisson distribution. This deviation can be exactly quantified using the non-Poissonian template fitting technique, and in this work we present the first application this approach to the IceCube high-energy neutrino dataset. Using this method, we search in 7 years of IceCube data for point-source populations correlated with the disk of the Milky Way, the Fermi bubbles, the Schlegel, Finkbeiner, and Davis dust map, or with the isotropic extragalactic sky. No evidence for such a population is found in the data using this technique, and in the absence of a signal we establish constraints on population models with source count distribution functions that can be described by a power-law with a single break. The derived limits can be interpreted in the context of many possible source classes. In order to enhance the flexibility of the results, we publish the full posterior from our analysis, which can be used to establish limits on specific population models that would contribute to the observed IceCube neutrino flux.

Contributions from the IceCube Collaboration presented at the 36th International Cosmic Ray Conference, 24 July - 1 August 2019, Madison, Wisconsin, USA.

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 5 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 $\sim$ 1 TeV to $\sim$ 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 well-described by a power-law in both track and cascade samples, and a best-fit index $\gamma=2.62\pm0.07$ is found in the energy range from 3.5 TeV to 2.6 PeV. Limits are set on the astrophysical flavor composition that are compatible with a ratio of $\left(\frac{1}{3}:\frac{1}{3}:\frac{1}{3}\right)_{\oplus}$. Exploiting the distinct inelasticity distribution of $\nu_{\mu}$ and $\bar{\nu}_{\mu}$ interactions, the atmospheric $\nu_{\mu}$ to $\bar{\nu}_{\mu}$ flux ratio in the energy range from 770 GeV to 21 TeV is found to be $0.77^{+0.44}_{-0.25}$ 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 May 2010 and May 2011 the IceCube detector operated with 79 strings deployed in the glacial ice at the South Pole, and with 86 strings between May 2011 and May 2015. 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\sigma$) 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 open the possibility to study cosmic ray transport near the Sun with future data from IceCube.

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 $\sim1^\circ$ at 1 TeV and improves to $\sim0.3^\circ$ 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 $\sim35\%$ assuming an $E^{-2}$ spectrum. The sensitivity on the muon-neutrino flux is at a level of $E^2 \mathrm{d} N /\mathrm{d} E = 3\cdot 10^{-13}\,\mathrm{TeV}\,\mathrm{cm}^{-2}\,\mathrm{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.

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 six years of IceCube data focusing on muon neutrino 'track' events from the Northern Hemisphere, while the second analysis uses two 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 a quasi-differential upper limit on the extremely-high-energy (EHE) neutrino flux above $5\times 10^{6}$ GeV based on an analysis of nine years of IceCube data. The astrophysical neutrino flux measured by IceCube extends to PeV energies, and it is a background flux when searching for an independent signal flux at higher energies, such as the cosmogenic neutrino signal. We have developed a new method to place robust limits on the EHE neutrino flux in the presence of an astrophysical background, whose spectrum has yet to be understood with high precision at PeV energies. A distinct event with a deposited energy above $10^{6}$ GeV was found in the new two-year sample, in addition to the one event previously found in the seven-year EHE neutrino search. These two events represent a neutrino flux that is incompatible with predictions for a cosmogenic neutrino flux and are considered to be an astrophysical background in the current study. The obtained limit is the most stringent to date in the energy range between $5 \times 10^{6}$ and $5 \times 10^{10}$ GeV. This result constrains neutrino models predicting 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}$. A significant part of the parameter-space for EHE neutrino production scenarios assuming a proton-dominated composition of ultra-high-energy cosmic rays is excluded.

A high-energy neutrino event detected by IceCube on 22 September 2017 was coincident in direction and time with a gamma-ray flare from the blazar TXS 0506+056. Prompted by this association, we investigated 9.5 years of IceCube neutrino observations to search for excess emission at the position of the blazar. We found an excess of high-energy neutrino events with respect to atmospheric backgrounds at that position between September 2014 and March 2015. Allowing for time-variable flux, this constitutes 3.5{\sigma} evidence for neutrino emission from the direction of TXS 0506+056, independent of and prior to the 2017 flaring episode. This suggests that blazars are the first identifiable sources of the high-energy astrophysical neutrino flux.

We report constraints on the sources of ultra-high-energy cosmic ray (UHECR) above $10^{9}$ GeV, based on an analysis of seven years of IceCube data. This analysis efficiently selects very high energy neutrino-induced events which have deposited energies from $\sim 10^6$ GeV to above $10^{11}$ GeV. Two neutrino-induced events with an estimated deposited energy of $(2.6 \pm 0.3) \times 10^6$ GeV, the highest neutrino energies observed so far, and $(7.7 \pm 2.0) \times 10^5$ GeV were detected. The atmospheric background-only hypothesis of detecting these events is rejected at 3.6$\sigma$. The hypothesis that the observed events are of cosmogenic origin is also rejected at $>$99% CL because of the limited deposited energy and the non-observation of events at higher energy, while their observation is consistent with an astrophysical origin. Our limits on cosmogenic neutrino fluxes disfavor the UHECR sources having cosmological evolution stronger than the star formation rate, e.g., active galactic nuclei and $\gamma$-ray bursts, assuming proton-dominated UHECRs. Constraints on UHECR sources including mixed and heavy UHECR compositions are obtained for models of neutrino production within UHECR sources. Our limit disfavors a significant part of parameter space for active galactic nuclei and new-born pulsar models.

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 \textasciitilde10$^2$ seconds from the northern hemisphere stacking search.

Neutrinos interact only very weakly, so they are extremely penetrating. However, the theoretical neutrino-nucleon interaction cross section rises with energy such that, at energies above 40 TeV, neutrinos are expected to be absorbed as they pass through the Earth. Experimentally, the cross section has been measured only at the relatively low energies (below 400 GeV) available at neutrino beams from accelerators \cite{Agashe:2014kda, Formaggio:2013kya}. Here we report the first measurement of neutrino absorption in the Earth, using a sample of 10,784 energetic upward-going neutrino-induced muons observed with the IceCube Neutrino Observatory. The flux of high-energy neutrinos transiting long paths through the Earth is attenuated compared to a reference sample that follows shorter trajectories through the Earth. Using a fit to the two-dimensional distribution of muon energy and zenith angle, we determine the cross section for neutrino energies between 6.3 TeV and 980 TeV, more than an order of magnitude higher in energy than previous measurements. The measured cross section is $1.30^{+0.21}_{-0.19}$ (stat.) $^{+0.39}_{-0.43}$ (syst.) times the prediction of the Standard Model \cite{CooperSarkar:2011pa}, consistent with the expectation for charged and neutral current interactions. We do not observe a dramatic increase in the cross section, expected in some speculative models, including those invoking new compact dimensions \cite{AlvarezMuniz:2002ga} or the production of leptoquarks \cite{Romero:2009vu}.

Papers on the properties of the atmospheric and astrophysical neutrino flux submitted to the 35th International Cosmic Ray Conference (ICRC 2017, Busan, South Korea) by the IceCube Collaboration

Papers on solar flares, supernovae, event reconstruction and education & outreach, submitted to the 35th International Cosmic Ray Conference (ICRC 2017, Busan, South Korea) by the IceCube Collaboration

Papers on searches for beyond the standard model physics, submitted to the 35th International Cosmic Ray Conference (ICRC 2017, Busan, South Korea) by the IceCube Collaboration

Papers on cosmic-ray measurements submitted to the 35th International Cosmic Ray Conference (ICRC 2017, Busan, South Korea) by the IceCube Collaboration