The observation of electromagnetic radiation from radio to γ-ray wavelengths has provided a wealth of information about the Universe. However, at PeV (1015 eV) energies and above, most of the Universe is impenetrable to photons. New messengers, namely cosmic neutrinos, are needed to explore the most extreme environments of the Universe where black holes, neutron stars, and stellar explosions transform gravitational energy into non-thermal cosmic rays. These energetic particles have millions of times higher energies than those produced in the most powerful particle accelerators on Earth. As neutrinos can escape from regions otherwise opaque to radiation, they allow an unique view deep into exploding stars and the vicinity of the event horizons of black holes. The discovery of cosmic neutrinos with IceCube has opened this new window on the Universe. IceCube has been successful in finding first evidence for cosmic particle acceleration in the jet of an active galactic nucleus. Yet, ultimately, its sensitivity is too limited to detect even the brightest neutrino sources with high significance, or to detect populations of less luminous sources. In this white paper, we present an overview of a next-generation instrument, IceCube-Gen2, which will sharpen our understanding of the processes and environments that govern the Universe at the highest energies. IceCube-Gen2 is designed to: (a) Resolve the high-energy neutrino sky from TeV to EeV energies (b) Investigate cosmic particle acceleration through multi-messenger observations (c) Reveal the sources and propagation of the highest energy particles in the Universe (d) Probe fundamental physics with high-energy neutrinos IceCube-Gen2 will enhance the existing IceCube detector at the South Pole. It will increase the annual rate of observed cosmic neutrinos by a factor of ten compared to IceCube, and will be able to detect sources five times fainter than its predecessor. Furthermore, through the addition of a radio array, IceCube-Gen2 will extend the energy range by several orders of magnitude compared to IceCube. Construction will take 8 years and cost about $350M. The goal is to have IceCube-Gen2 fully operational by 2033. IceCube-Gen2 will play an essential role in shaping the new era of multi-messenger astronomy, fundamentally advancing our knowledge of the high-energy Universe. This challenging mission can be fully addressed only through the combination of the information from the neutrino, electromagnetic, and gravitational wave emission of high-energy sources, in concert with the new survey instruments across the electromagnetic spectrum and gravitational wave detectors which will be available in the coming years.

LS I +61 303 is one of only a few high-mass X-ray binaries currently detected at high significance in very high energy -rays. The system was observed over several orbital cycles (between 2006 September and 2007 February) with the VERITAS array of imaging air Cerenkov telescopes. A signal of -rays with energies above 300 GeV is found with a statistical significance of 8.4 standard deviations. The detected flux is measured to be strongly variable; the maximum flux is found during most orbital cycles at apastron. The energy spectrum for the period of maximum emission can be characterized by a power law with a photon index of À ¼ 2:40 AE 0:16 stat AE 0:2 sys and a flux above 300 GeV corresponding to 15%-20% of the flux from the Crab Nebula.

We report on TeV γ -ray observations of the blazar Mrk 421 (redshift of 0.031) with the VERITAS observatory and the Whipple 10 m Cherenkov telescope. The excellent sensitivity of VERITAS allowed us to sample the TeV γ -ray fluxes and energy spectra with unprecedented accuracy where Mrk 421 was detected in each of the pointings. A total of 47.3 hr of VERITAS and 96 hr of Whipple 10 m data were acquired between 2006 January and 2008 June. We present the results of a study of the TeV γ -ray energy spectra as a function of time and for different flux levels. On 2008 May 2 and 3, bright TeV γ -ray flares were detected with fluxes reaching the level of 10 Crab. The TeV γ -ray data were complemented with radio, optical, and X-ray observations, with flux variability found in all bands except for the radio wave band. The combination of the Rossi X-ray Timing Explorer and Swift X-ray data reveal spectral hardening with increasing flux levels, often correlated with an increase of the source activity in TeV 1 The Astrophysical Journal, 738:25 (19pp), 2011 September 1 Acciari et al.γ -rays. Contemporaneous spectral energy distributions were generated for 18 nights, each of which are reasonably described by a one-zone synchrotron self-Compton model.

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 stateof-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 a are used to discuss the capabilities of the code and different aspects of instrumental design decisions.Keywords Neutrino astronomy · radio detection · simulation · signal processing · ice propagation · Askaryan · PACS 07.05.Kf · 95.85.Ry · 95.55.Vj · 95.85.Bh · 07.05.Tp

The ANtarctic Impulsive Transient Antenna (ANITA) NASA long-duration balloon payload completed its fourth flight in December 2016, after 28 days of flight time. ANITA is sensitive to impulsive broadband radio emission from interactions of ultra-high-energy neutrinos in polar ice (Askaryan emission). We present the results of two separate blind analyses searching for signals from Askaryan emission in the data from the fourth flight of ANITA. The more sensitive analysis, with a better expected limit, has a background estimate of 0.64 +0.69 −0.45 and an analysis efficiency of 82±2%. The second analysis has a background estimate of 0.34 +0.66 −0.16 and an analysis efficiency of 71±6%. Each analysis found one event in the signal region, consistent with the background estimate for each analysis. The resulting limit further tightens the constraints on the diffuse flux of ultra-high-energy neutrinos at energies above 10 19.5 eV.

The first flight of the Antarctic Impulsive Transient Antenna (ANITA) experiment recorded 16 radio signals that were emitted by cosmic-ray induced air showers. The dominant contribution to the radiation comes from the deflection of positrons and electrons in the geomagnetic field, which is beamed in the direction of motion of the air shower. For 14 of these events, this radiation is reflected from the ice and subsequently detected by the ANITA experiment at a flight altitude of ∼36 km. In this paper, we estimate the energy of the 14 individual events and find that the mean energy of the cosmic-ray sample is 2.9 × 10 18 eV, which is significantly lower than the previous estimate. By simulating the ANITA flight, we calculate its exposure for ultra-high energy cosmic rays. We estimate for the first time the cosmic-ray flux derived only from radio observations and find agreement with measurements performed at other observatories. In addition, we find that the ANITA data set is consistent with Monte-Carlo simulations for the total number of observed events and with the properties of those events.

The Antarctic Impulsive Transient Antenna, a NASA long-duration balloon payload, searches for radio emission from interactions of ultrahigh-energy neutrinos in polar ice. The third flight of the Antarctic Impulsive Transient Antenna was launched in December 2014 and completed a 22-day flight. We present the results of three analyses searching for Askaryan radio emission of neutrino origin. In the most sensitive of the analyses, we find one event in the signal region on an expected background of 0.7 þ0.5 −0.3 . Though consistent with the background estimate, the event remains compatible with a neutrino hypothesis even after additional postunblinding scrutiny.

Indirect dark matter searches with ground-based gamma-ray observatories provide an alternative for identifying the particle nature of dark matter that is complementary to that of direct search or accelerator production experiments. We present the results of observations of the dwarf spheroidal galaxies Draco, Ursa Minor, Boötes 1, and Willman 1 conducted by the Very Energetic Radiation Imaging Telescope Array System (VERITAS). These galaxies are nearby dark matter dominated objects located at a typical distance of several tens of kiloparsecs for which there are good measurements of the dark matter density profile from stellar velocity measurements. Since the conventional astrophysical background of very high energy gamma rays from these objects appears to be negligible, they are good targets to search for the secondary gamma-ray photons produced by interacting or decaying dark matter particles. No significant gamma-ray flux above 200 GeV was detected from these four dwarf galaxies for a typical exposure of ∼20 hr. The 95% confidence upper limits on the integral gamma-ray flux are in the range (0.4-2.2) × 10 −12 photons cm −2 s −1 . We interpret this limiting flux in the context of pair annihilation of weakly interacting massive particles (WIMPs) and derive constraints on the thermally averaged product of the total self-annihilation cross section and the relative velocity of the WIMPs ( σ v 10 −23 cm 3 s −1 for m χ 300 GeV c −2 ). This limit is obtained under conservative assumptions regarding the dark matter distribution in dwarf galaxies and is approximately 3 orders of magnitude above the generic theoretical prediction for WIMPs in the minimal supersymmetric standard model framework. However, significant uncertainty exists in the dark matter distribution as well as the neutralino cross sections which under favorable assumptions could further lower this limit.