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\sigma$ 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.

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.

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 $\mathrm{GeV~cm^{-2}}$. We also set limits on the total isotropic equivalent energy, $E_{\mathrm{iso}}$, emitted in high-energy neutrinos by each GW event. These limits range from 1.7 $\times$ 10$^{51}$ - 1.8 $\times$ 10$^{55}$ erg. We conclude with an outlook for LIGO-Virgo observing run O3, during which both analyses are running in real time.

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.

Pulsar wind nebulae (PWNe) are main gamma-ray emitters in the Galactic plane. They are diffuse nebulae that emit non-thermal radiation. Pulsar winds, relativistic magnetized outflows from the central star, shocked in the ambient medium produce a multiwavelength emission from the radio through gamma rays. Although the leptonic scenario is able to explain most PWNe emission, a hadronic contribution cannot be excluded. A possible hadronic contribution to the high-energy gamma-ray emission inevitably leads to the production of neutrinos. Using 9.5 years of all-sky IceCube data, we report results from a stacking analysis to search for neutrino emission from 35 PWNe that are high-energy gamma-ray emitters. In the absence of any significant correlation, we set upper limits on the total neutrino emission from those PWNe and constraints on hadronic spectral components.

We present two searches for IceCube neutrino events coincident with 28 fast radio bursts (FRBs) and one repeating FRB. The first improves upon a previous IceCube analysis -- searching for spatial and temporal correlation of events with FRBs at energies greater than roughly 50 GeV -- by increasing the effective area by an order of magnitude. The second is a search for temporal correlation of MeV neutrino events with FRBs. No significant correlation is found in either search, therefore, we set upper limits on the time-integrated neutrino flux emitted by FRBs for a range of emission timescales less than one day. These are the first limits on FRB neutrino emission at the MeV scale, and the limits set at higher energies are an order-of-magnitude improvement over those set by any neutrino telescope.

We present the results of the first combined dark matter search targeting the Galactic Centre using the ANTARES and IceCube neutrino telescopes. For dark matter particles with masses from 50 to 1000 GeV, the sensitivities on the self-annihilation cross section set by ANTARES and IceCube are comparable, making this mass range particularly interesting for a joint analysis. Dark matter self-annihilation through the $\tau^+\tau^-$, $\mu^+\mu^-$, $b\bar{b}$ and $W^+W^-$ channels is considered for both the Navarro-Frenk-White and Burkert halo profiles. In the combination of 2,101.6 days of ANTARES data and 1,007 days of IceCube data, no excess over the expected background is observed. Limits on the thermally-averaged dark matter annihilation cross section $\langle\sigma_A\upsilon\rangle$ are set. These limits present an improvement of up to a factor of two in the studied dark matter mass range with respect to the individual limits published by both collaborations. When considering dark matter particles with a mass of 200 GeV annihilating through the $\tau^+\tau^-$ channel, the value obtained for the limit is $7.44 \times 10^{-24} \text{cm}^{3}\text{s}^{-1}$ for the Navarro-Frenk-White halo profile. For the purpose of this joint analysis, the model parameters and the likelihood are unified, providing a benchmark for forthcoming dark matter searches performed by neutrino telescopes.

We present the variability study of the lowest-luminosity Seyfert 1 galaxy NGC 4395 based on the photometric monitoring campaigns in 2017 and 2018. Using 22 ground-based and space telescopes, we monitored NGC 4395 with a $\sim$5 minute cadence during a period of 10 days and obtained light curves in the UV, V, J, H, and K/Ks bands as well as the H$\alpha$ narrow-band. The RMS variability is $\sim$0.13 mag on \emph{Swift}-UVM2 and V filter light curves, decreasing down to $\sim$0.01 mag on K filter. After correcting for continuum contribution to the H$\alpha$ narrow-band, we measured the time lag of the H$\alpha$ emission line with respect to the V-band continuum as ${55}^{+27}_{-31}$ to ${122}^{+33}_{-67}$ min. in 2017 and ${49}^{+15}_{-14}$ to ${83}^{+13}_{-14}$ min. in 2018, depending on the assumption on the continuum variability amplitude in the H$\alpha$ narrow-band. We obtained no reliable measurements for the continuum-to-continuum lag between UV and V bands and among near-IR bands, due to the large flux uncertainty of UV observations and the limited time baseline. We determined the AGN monochromatic luminosity at 5100\AA\ $\lambda L_\lambda = \left(5.75\pm0.40\right)\times 10^{39}\,\mathrm{erg\,s^{-1}}$, after subtracting the contribution of the nuclear star cluster. While the optical luminosity of NGC 4395 is two orders of magnitude lower than that of other reverberation-mapped AGNs, NGC 4395 follows the size-luminosity relation, albeit with an offset of 0.48 dex ($\geq$2.5$\sigma$) from the previous best-fit relation of Bentz et al. (2013).

The measurement of diffuse PeV gamma-ray emission from the Galactic plane would provide information about the energy spectrum and propagation of Galactic cosmic rays, and the detection of a point-like source of PeV gamma rays would be strong evidence for a Galactic source capable of accelerating cosmic rays up to at least a few PeV. This paper presents several un-binned maximum likelihood searches for PeV gamma rays in the Southern Hemisphere using 5 years of data from the IceTop air shower surface detector and the in-ice array of the IceCube Observatory. The combination of both detectors takes advantage of the low muon content and deep shower maximum of gamma-ray air showers, and provides excellent sensitivity to gamma rays between $\sim$0.6 PeV and 100 PeV. Our measurements of point-like and diffuse Galactic emission of PeV gamma rays are consistent with background, so we constrain the angle-integrated diffuse gamma-ray flux from the Galactic Plane at 2 PeV to $2.61 \times 10^{-19}$ cm$^{-2}$ s$^{-1}$ TeV$^{-1}$ at 90% confidence, assuming an E$^{-3}$ spectrum, and we estimate 90% upper limits on point-like emission at 2 PeV between 10$^{-21}$ - 10$^{-20}$ cm$^{-2}$ s$^{-1}$ TeV$^{-1}$ for an E$^{-2}$ spectrum, depending on declination. Furthermore, we exclude unbroken power-law emission up to 2 PeV for several TeV gamma-ray sources observed by H.E.S.S., and calculate upper limits on the energy cutoffs of these sources at 90% confidence. We also find no PeV gamma rays correlated with neutrinos from IceCube's high-energy starting event sample. These are currently the strongest constraints on PeV gamma-ray emission.

Twenty three years ago, cosmological N-body simulations revealed quasi-universal NFW dark matter halos, whose physical origin is still unclear. This work tries to solve this issue by Landau-Ginzburg (LG) theory in equilibrium statistical mechanics. We replace the order parameter in LG theory by density and consider the dark matter halos as fluctuations from the equilibrium state of the background, which can be the main halos for the subhalos, or the homogeneous and isotropic Universe for the galaxy cluster halos. With the assumption that the fluctuations of Helmholtz free energy depend mainly on the density and its gradient, which is shown to be consistent with the behavior of the power spectrum near the cluster scale, we can easily obtain the universal $r^{-1}$ cusp of the equilibrium dark matte halo without other requirements. This work strongly suggests that more studies should be done for the self-gravitating systems by statistical mechanics.

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$).

A search for point-like and extended sources of cosmic neutrinos using data collected by the ANTARES and IceCube neutrino telescopes is presented. The data set consists of all the track-like and shower-like events pointing in the direction of the Southern Sky included in the nine-year ANTARES point-source analysis, combined with the through-going track-like events used in the seven-year IceCube point-source search. The advantageous field of view of ANTARES and the large size of IceCube are exploited to improve the sensitivity in the Southern Sky by a factor $\sim$2 compared to both individual analyses. In this work, the Southern Sky is scanned for possible excesses of spatial clustering, and the positions of preselected candidate sources are investigated. In addition, special focus is given to the region around the Galactic Centre, whereby a dedicated search at the location of SgrA* is performed, and to the location of the supernova remnant RXJ 1713.7-3946. No significant evidence for cosmic neutrino sources is found and upper limits on the flux from the various searches are presented.

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.

Over the past 20 years, KASCADE and its extension KASCADE-Grande were dedicated to measure high-energy cosmic rays with primary energies of 100 TeV to 1 EeV. The data accumulation was fully completed and all experimental components were dismantled, though the analysis of the high-quality data is still continued. E.g., we investigated the validity of the hadronic interaction model of the new SIBYLL version 2.3c. We also published a new result of a search for large-scale anisotropies performed with the KASCADE-Grande data. Investigation of the attenuation length of the muon in the atmosphere is also updated with the predictions of the SIBYLL 2.3 interaction model. We investigated, in addition, the muon content of high-energy air showers and compared them to all post-LHC interaction models. In this contribution, the new and updated results from KASCADE-Grande will be presented. An update of the web-based data center KCDC offering the original scientific data from KASCADE-Grande to the public will be briefly discussed as well.

In this paper we describe the first results of a compact imaging air-Cherenkov telescope, IceAct, operating in coincidence with the IceCube Neutrino Observatory (IceCube) at the geographic South Pole. An array of IceAct telescopes (referred to as the IceAct project) is under consideration as part of the IceCube-Gen2 extension to IceCube. Surface detectors in general will be a powerful tool in IceCube-Gen2 for distinguishing astrophysical neutrinos from the dominant backgrounds of cosmic-ray induced atmospheric muons and neutrinos: the IceTop array is already in place as part of IceCube, but has a high energy threshold. Although the duty cycle will be lower for the IceAct telescopes than the present IceTop tanks, the IceAct telescopes may prove to be more effective at lowering the detection threshold for air showers. Additionally, small imaging air-Cherenkov telescopes in combination with IceTop, the deep IceCube detector or other future detector systems might improve measurements of the composition of the cosmic ray energy spectrum. In this paper we present measurements of a first 7-pixel imaging air Cherenkov telescope demonstrator, proving the capability of this technology to measure air showers at the South Pole in coincidence with IceTop and the deep IceCube detector.

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.

The past decade has welcomed the emergence of cosmic neutrinos as a new messenger to explore the most extreme environments of the universe. The discovery measurement of cosmic neutrinos, announced by IceCube in 2013, has opened a new window of observation that has already resulted in new fundamental information that holds the potential to answer key questions associated with the high-energy universe, including: what are the sources in the PeV sky and how do they drive particle acceleration; where are cosmic rays of extreme energies produced, and on which paths do they propagate through the universe; and are there signatures of new physics at TeV-PeV energies and above? The planned advancements in neutrino telescope arrays in the next decade, in conjunction with continued progress in broad multimessenger astrophysics, promise to elevate the cosmic neutrino field from the discovery to the precision era and to a survey of the sources in the neutrino sky. The planned detector upgrades to the IceCube Neutrino Observatory, culminating in IceCube-Gen2 (an envisaged $400M facility with anticipated operation in the next decade, described in this white paper) are the cornerstone that will drive the evolution of neutrino astrophysics measurements.

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.

Active galactic nuclei (AGNs) show a correlation between the size of the broad line region (BLR) and the monochromatic continuum luminosity at 5100 \AA, allowing black hole mass estimation based on single-epoch spectra. However, the validity of the correlation is yet to be clearly tested for high-luminosity AGNs. We present the first reverberation-mapping results of the Seoul National University AGN monitoring program (SAMP), which is designed to focus on luminous AGNs for probing the high end of the size-luminosity relation. We report time lag measurements of two AGNs, namely, 2MASS J10261389+5237510 and SDSS J161911.24+501109.2, using the light curves obtained over a $\sim$1000 day period with an average cadence of $\sim$10 and $\sim$20 days, respectively for photometry and spectroscopy monitoring. Based on a cross-correlation analysis and H$\beta$ line width measurements, we determine the H$\beta$ lag as $41.8^{+4.9}_{-6.0}$ and $52.6^{+17.6}_{-14.7}$ days in the observed-frame, and black hole mass as $3.65^{+0.49}_{-0.57} \times 10^7 M_{\odot}$ and $23.02^{+7.81}_{-6.56} \times 10^7 M_{\odot}$, respectively for 2MASS J1026 and SDSS J1619.

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.

Efficient treatment of systematic uncertainties that depend on a large number of nuisance parameters is a persistent difficulty in particle physics experiments. Where low-level effects are not amenable to simple parameterization or re-weighting, analyses often rely on discrete simulation sets to quantify the effects of nuisance parameters on key analysis observables. Such methods may become computationally untenable for analyses requiring high statistics Monte Carlo with a large number of nuisance degrees of freedom, especially in cases where these degrees of freedom parameterize the shape of a continuous distribution. In this paper we present a method for treating systematic uncertainties in a computationally efficient and comprehensive manner using a single simulation set with multiple and continuously varied nuisance parameters. This method is demonstrated for the case of the depth-dependent effective dust distribution within the IceCube Neutrino Telescope.

A data life cycle (DLC) is a high-level data processing pipeline that involves data acquisition, event reconstruction, data analysis, publication, archiving, and sharing. For astroparticle physics a DLC is particularly important due to the geographical and content diversity of the research field. A dedicated and experiment spanning analysis and data centre would ensure that multi-messenger analyses can be carried out using state-of-the-art methods. The German-Russian Astroparticle Data Life Cycle Initiative (GRADLCI) is a joint project of the KASCADE-Grande and TAIGA collaborations, aimed at developing a concept and creating a DLC prototype that takes into account the data processing features specific for the research field. An open science system based on the KASCADE Cosmic Ray Data Centre (KCDC), which is a web-based platform to provide the astroparticle physics data for the general public, must also include effective methods for distributed data storage algorithms and techniques to allow the community to perform simulations and analyses with sophisticated machine learning methods. The aim is to achieve more efficient analyses of the data collected in different, globally dispersed observatories, as well as a modern education to Big Data Scientist in the synergy between basic research and the information society. The contribution covers the status and future plans of the initiative.

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