We present the detection of VHE gamma-ray emission above 100 TeV from HAWC J2227+610 with the HAWC observatory. Combining our observations with previously published results by VERITAS, we interpret the gamma-ray emission from HAWC J2227+610 as emission from protons with a lower limit in their cutoff energy of 800 TeV. The most likely source of the protons is the associated supernova remnant G106.3+2.7, making it a good candidate for a Galactic PeVatron. However, a purely leptonic origin of the observed emission cannot be excluded at this time.

Following the coalescence of binary neutron stars, debris from the merger which remains marginally bound to the central compact remnant will fallback at late times, feeding a sustained accretion flow. Unbound winds or a wide-angle jet from this radiatively-inefficient disk may collide with the comparatively slow dense kilonova ejecta released from an earlier phase. Under the assumption that such interaction accelerate cosmic rays to ultra-high energies, we numerically simulate their propagation and interactions through the dynamical ejecta. The hadronuclear and photo-hadronic processes experienced by particles produce isotropic high-energy neutrino fluxes, peaking at times $10^{3-4}\,$s, which we calculate for two sets of parameters. A first set is inspired by the observations of GW170817. In the second scenario, which we call optimistic, parameters are chosen so as to optimize the neutrino flux, within the range allowed by observation and theory. We find that single sources can only be detected with IceCube-Gen2 for optimistic scenarios and if located within $\sim 4\,$Mpc. The cumulative flux could contribute to $\sim 0.5-10\%$ of the diffuse flux observed by the IceCube Observatory, depending on the fall-back power and the cosmic ray composition. The neutrino emission powered by fallback is nearly isotropic, and can be used for future correlation studies with gravitational wave signals.

The hardening and softening features in the DAMPE proton spectrum are very likely to be originated from a nearby supernova remnant (SNR). The proton spectrum from the nearby SNR is required to be very hard below $\approx10$ TeV. To reproduce this feature, we illustrate that anomalously slow-diffusion zone for cosmic rays (CRs) must be existed in the local interstellar medium (ISM) after also taking the dipole anisotropy of CRs into account. Assuming that the diffusion coefficient is homogeneous in the nearby ISM, we show that the diffusion coefficient is constrained to the magnitude of $10^{26}$ cm$^2$ s$^{-1}$ when normalized to 1 GeV, which is about 100 times smaller than the average value in the Galaxy. We further discuss the spatial distribution of the slow diffusion and find two distinct possibilities. In one case, the SNR is several hundred of parsecs away from the solar system, meanwhile both the SNR and the solar system are required to be included in a large slow-diffusion zone. The homogeneous diffusion belongs to this case. In the other case, the SNR is very close with a distance of $\sim50$ pc and the slow-diffusion zone is only limited around the SNR. The required diffusion coefficient is further smaller in the latter case. This work provides a new way of studying the CR diffusion in the local ISM.

The origin of the bulk of the astrophysical neutrinos detected by the IceCube Observatory remains a mystery. Previous source-finding analyses compare the directions of IceCube events and individual sources in astrophysical catalogs. The source association method is technically challenging when the number of source candidates is much larger than the number of the observed astrophysical neutrinos. We show that in this large source number regime, a cross-correlation analysis of neutrino data and source catalog can instead be used to constrain potential source populations for the high-energy astrophysical neutrinos, and provide spatial evidence for the existence of astrophysical neutrinos. We present an analysis of the cross-correlation of the IceCube 2010-2012 point-source data and a WISE-2MASS galaxy sample. While we find no significant detection of cross-correlation with the publicly available neutrino dataset, we show that, when applied to the full IceCube data, which has a longer observation time and higher astrophysical neutrino purity, our method has sufficient statistical power to detect a cross-correlation signal if the neutrino sources trace the Large Scale Structure of the Universe.

The extended jets of the microquasar SS 433 have been observed in optical, radio, X-ray, and recently very-high-energy (VHE) $\gamma$-rays by HAWC. The detection of HAWC $\gamma$-rays with energies as great as 25 TeV motivates searches for high-energy $\gamma$-ray counterparts in the Fermi-LAT data in the 100 MeV--300 GeV band. In this paper, we report on the first-ever joint analysis of Fermi-LAT and HAWC observations to study the spectrum and location of $\gamma$-ray emission from SS~433. Our analysis finds common emission sites of GeV-to-TeV $\gamma$-rays inside the eastern and western lobes of SS 433. The total flux above 1 GeV is $\sim 1\times10^{-10}\,\rm cm^{-2}\,s^{-1}$ in both lobes. The $\gamma$-ray spectrum in the eastern lobe is consistent with inverse-Compton emission by an electron population that is accelerated by jets. To explain both the GeV and TeV flux, the electrons need to have a soft intrinsic energy spectrum, or undergo a quick cooling process due to synchrotron radiation in a magnetized environment.

We investigate prolonged engine activities of short gamma-ray bursts (SGRBs), such as extended and/or plateau emissions, as high-energy gamma-ray counterparts to gravitational waves (GWs). Binary neutron-star mergers lead to relativistic jets and merger ejecta with $r$-process nucleosynthesis, which are observed as SGRBs and kilonovae/macronovae, respectively. Long-term relativistic jets may be launched by the merger remnant as hinted in X-ray light curves of some SGRBs. The prolonged jets may dissipate their kinetic energy within the radius of the cocoon formed by the jet-ejecta interaction. Then the cocoon supplies seed photons to non-thermal electrons accelerated at the dissipation region, causing high-energy gamma-ray production through the inverse Compton scattering process. We numerically calculate high-energy gamma-ray spectra in such a system using a one-zone and steady-state approximation, and show that GeV--TeV gamma-rays are produced with a duration of $10^2-10^5$ seconds. They can be detected by {\it Fermi}/LAT or CTA as gamma-ray counterparts to GWs.

Tidal disruptions are extremely powerful phenomena that have been designated as candidate sources of ultra-high-energy cosmic rays. The disruption of a star by a black hole can naturally provide protons and heavier nuclei, which can be injected and accelerated to ultra-high energies within a jet. Inside the jet, accelerated nuclei are likely to interact with a dense photon field, leading to a significant production of neutrinos and secondary particles. We model numerically the propagation and interactions of high-energy nuclei in jetted tidal disruption events in order to evaluate consistently their signatures in cosmic rays and neutrinos. We propose a simple model of the light curve of tidal disruption events, consisting of two stages: a high state with bright luminosity and short duration and a medium state, less bright and longer lasting. These two states have different impacts on the production of cosmic rays and neutrinos. In order to calculate the diffuse fluxes of cosmic rays and neutrinos, we model the luminosity function and redshift evolution of jetted tidal disruption events. We find that we can fit the latest ultra-high-energy cosmic-ray spectrum and composition results of the Auger experiment for a range of reasonable parameters. The diffuse neutrino flux associated with this scenario is found to be subdominant, but nearby events can be detected by IceCube or next-generation detectors such as IceCube-Gen2.

Geminga pulsar is surrounded by a multi-TeV $\gamma$-ray halo radiated by the high energy electrons and positrons accelerated by the central pulsar wind nebula (PWN). The angular profile of the $\gamma$-ray emission reported by HAWC indicates an anomalously slow diffusion for the cosmic-ray electrons and positrons in the halo region around Geminga. In the paper we study the possible mechanism for the origin of the slow diffusion. At first, we consider the self-generated Alfv\'en waves due to the streaming instability of the electrons and positrons released by Geminga. However, even considering a very optimistic scenario for the wave growth, we find this mechanism DOES NOT work to account for the extremely slow diffusion at the present day if taking the proper motion of Geminga pulsar into account. The reason is straightforward as the PWN is too weak to generate enough high energy electrons and positrons to stimulate strong turbulence at the late time. We then propose an assumption that the strong turbulence is generated by the shock wave of the parent supernova remnant (SNR) of Geminga. Geminga may still be inside the SNR, and we find that the SNR can provide enough energy to generate the slow-diffusion circumstance. The TeV halos around PSR B0656+14, Vela X, and PSR J1826-1334 may also be explained under this assumption.

The Giant Radio Array for Neutrino Detection (GRAND) is a planned large-scale observatory of ultra-high-energy (UHE) cosmic particles, with energies exceeding 10^8 GeV. Its goal is to solve the long-standing mystery of the origin of UHE cosmic rays. To do this, GRAND will detect an unprecedented number of UHE cosmic rays and search for the undiscovered UHE neutrinos and gamma rays associated to them with unmatched sensitivity. GRAND will use large arrays of antennas to detect the radio emission coming from extensive air showers initiated by UHE particles in the atmosphere. Its design is modular: 20 separate, independent sub-arrays, each of 10 000 radio antennas deployed over 10 000 km^2. A staged construction plan will validate key detection techniques while achieving important science goals early. Here we present the science goals, detection strategy, preliminary design, performance goals, and construction plans for GRAND.

As the TeV halos around Geminga and PSR B0656+14 have been confirmed by HAWC, slow diffusion of cosmic rays could be general around pulsars, and the cosmic positron spectrum from pulsars could be significantly changed. As a consequence, the most likely pulsar source of the positron excess, Geminga, is no more a viable candidate under the additional constraint from Fermi-LAT. Moreover, the latest measurement by AMS-02 shows a clear cutoff in the positron spectrum, which sets a strict constraint on the age of the pulsar source. Considering these new developments we reanalyze the scenario in this work. By checking all the observed pulsars under the two-zone diffusion scenario, we propose for the first time that PSR B1055-52 is a very promising source of the positron excess. B1055-52 can well reproduce both the intensity and the high-energy cut of the AMS-02 positron spectrum, and may also explain the H.E.S.S $e^-+e^+$ spectrum around 10 TeV. Moreover, if the slow diffusion is universal in the local interstellar medium, B1055-52 will be the unique reasonable source of the AMS-02 positron spectrum among the observed pulsars.

Newly-born, rapidly-spinning magnetars have been invoked as the power sources of super-luminous transients, including the class of "fast-luminous optical transients" (FBOTs). The extensive multi-wavelength analysis of AT2018cow, the first FBOT discovered in real time, is consistent with the magnetar scenario and offers an unprecedented opportunity to comprehend the nature of these sources and assess their broader implications. Using AT2018cow as a prototype, we investigate high-energy neutrino and cosmic ray production from FBOTs and the more general class of superluminous supernovae (SLSNe). By calculating the interaction of cosmic rays and the time-evolving radiation field and baryon background, we find that particles accelerated in the magnetar wind may escape the ejecta at ultrahigh energies (UHE). The predicted high-energy neutrino fluence from AT2018cow is below the sensitivity of the IceCube Observatory, and estimates of the cosmically-integrated neutrino flux from FBOTs are consistent with the extreme-high-energy upper limits posed by IceCube. High-energy $\gamma$ rays exceeding GeV energies are obscured for the first months to years by thermal photons in the magnetar nebula, but are potentially observable at later times. Given also their potentially higher volumetric rate compared to other engine-powered transients (e.g. SLSNe and gamma-ray bursts), we conclude that FBOTs are favorable targets for current and next-generation multi-messenger observatories.

We report on the highest energy photons from the Crab Nebula observed by the Tibet air shower array with the underground water-Cherenkov-type muon detector array. Based on the criterion of muon number measured in an air shower, we successfully suppress 99.92% of the cosmic-ray background events with energies $E>100$ TeV. As a result, we observed 24 photon-like events with $E>100$ TeV against 5.5 background events, which corresponds to 5.6$\sigma$ statistical significance. This is the first detection of photons with $E>100$ TeV from an astrophysical source.

It has long been suggested that nearby pulsars within $\sim 1 \,{\rm kpc}$ are the leading candidate of the 10-500 GeV cosmic-ray positron excess measured by PAMELA and other experiments. The recent measurement of surface brightness profile of TeV nebulae surrounding Geminga and PSR~B0656+14 by the High-Altitude Water Cherenkov Observatory (HAWC) suggests inefficient diffusion of particles from the sources, giving rise to a debate on the pulsar interpretation of the cosmic-ray positron excess. Here we argue that GeV observations provide more direct constraints on the positron density in the TeV nebulae in the energy range of 10-500 GeV and hence on the origin of the observed positron excess. Motivated by this, we search for GeV emission from the TeV nebulae with the \textsl{Fermi} Large Area Telescope (LAT). No spatially-extended GeV emission is detected from these two TeV nebulae in the framework of two-zone diffusion spatial templates, suggesting a relatively low density of GeV electrons/positrons in the TeV nebulae. A joint modelling of the data from HAWC and \textsl{Fermi}-LAT disfavors Geminga and PSR~B0656+14 as the dominant source of the positron excess at $\sim 50-500$ GeV for the usual Kolmogorov-type diffusion, while for an energy-independent diffusion, a dominant part of the positron excess contributed by them cannot be ruled out by the current data.

We study the time-dependent modulation effect and derive the local interstellar spectra (LIS) for the cosmic ray (CR) proton, helium, boron and carbon. A two-dimensional modulation model including the variation of the interplanetary environment with time is adopted to describe modulation process. The propagation equation of CRs in the heliosphere is numerically solved by the package Solarprop. We derive the LIS by fitting the latest results of several experiments, including Voyager 1, PAMELA, BESS-POLARII and ACE, during low solar activity periods. We further study the modulation in the polarity reversal periods with the PAMELA proton data. We find that the rigidity dependence of the diffusion coefficient is critical to explain the modulation effect during reversal periods. Our results also indicate a power law relation between the diffusion coefficient and the magnitude of the heliospheric magnetic field (HMF) at the Earth.

We review open questions and prospects for progress in ultrahigh-energy cosmic ray (UHECR) research, based on a series of discussions that took place during the `The High-Energy Universe: Gamma-Ray, Neutrino, and Cosmic-ray Astronomy' MIAPP workshop in 2018. Specifically, we overview open questions on the origin of the bulk of UHECRs, the UHECR mass composition, the origin of the end of the cosmic-ray spectrum, the transition from Galactic to extragalactic cosmic rays, the effect of magnetic fields on the trajectories of UHECRs, anisotropy expectations for specific astrophysical scenarios, hadronic interactions, and prospects for discovering neutral particles as well as new physics at ultrahigh energies. We also briefly present upcoming and proposed UHECR experiments and discuss their projected science reach.

The possible association of the blazar TXS 0506+056 with a high-energy neutrino detected by IceCube holds the tantalizing potential to answer three astrophysical questions: 1. Where do high-energy neutrinos originate? 2. Where are cosmic rays produced and accelerated? 3. What radiation mechanisms produce the high-energy {\gamma}-rays in blazars? The MeV gamma-ray band holds the key to these questions, because it is an excellent proxy for photo-hadronic processes in blazar jets, which also produce neutrino counterparts. Variability in MeV gamma-rays sheds light on the physical conditions and mechanisms that take place in the particle acceleration sites in blazar jets. In addition, hadronic blazar models also predict a high level of polarization fraction in the MeV band, which can unambiguously distinguish the radiation mechanism. Future MeV missions with a large field of view, high sensitivity, and polarization capabilities will play a central role in multi-messenger astronomy, since pointed, high-resolution telescopes will follow neutrino alerts only when triggered by an all-sky instrument.

This Astro2020 white paper advocates for a multi-messenger approach that combines high-energy neutrino and broad multi-wavelength electromagnetic observations to study AGN during the coming decade. The unique capabilities of these joint observations promise to solve several long-standing issues in our understanding of AGN as powerful cosmic accelerators.

The characteristic two-component blazar spectral energy distribution (SED) can be of either leptonic and/or hadronic origins. The potential association of the high-energy neutrino event IceCube-170922A with the flaring blazar TXS~0506+056 indicates that hadronic processes may operate in a blazar jet. Despite multi-wavelength follow-ups of the event and extensive theoretical modelings, the radiation mechanisms and the underlying magnetic field strength and configuration remain poorly understood. In this paper, we consider generic leptonic and hadronic blazar spectral models with distinct magnetic field strengths and radiation mechanisms. We analytically reproduce the SEDs and the neutrino flux of hadronic models, and predict their X-ray to $\gamma$-ray polarization degrees. Furthermore, by performing relativistic magnetohydrodynamic (RMHD) simulations taking into account the polarization-dependent radiation transfer, we study the time-dependent multi-wavelength polarization variability of the proton synchrotron model under a shock scenario. Our results suggest that the high-energy polarization degree and the neutrino flux can be jointly used to pinpoint the leptonic and/or hadronic blazar radiation mechanisms in the X-ray and $\gamma$-ray bands, and to infer the magnetic field strength in the emission region. Additionally, the temporal multi-wavelength polarization signatures in the proton synchrotron model shed light on the jet energy composition and the dynamical importance of magnetic fields in the blazar emission region. Future multi-wavelength polarimetry facilities such as {\it IXPE} and {\it AMEGO} together with neutrino telescopes such as {\it IceCube} can provide unprecedented observational constraints to probe the blazar radiation mechanisms and jet dynamics.

The origin of ultrahigh-energy cosmic rays (UHECRs) is a half-century old enigma (Linsley 1963). The mystery has been deepened by an intriguing coincidence: over ten orders of magnitude in energy, the energy generation rates of UHECRs, PeV neutrinos, and isotropic sub-TeV gamma rays are comparable, which hints at a grand-unified picture (Murase and Waxman 2016). Here we report that powerful black hole jets in aggregates of galaxies can supply the common origin of all of these phenomena. Once accelerated by a jet, low-energy cosmic rays confined in the radio lobe are adiabatically cooled; higher-energy cosmic rays leaving the source interact with the magnetized cluster environment and produce neutrinos and gamma rays; the highest-energy particles escape from the host cluster and contribute to the observed cosmic rays above 100 PeV. The model is consistent with the spectrum, composition, and isotropy of the observed UHECRs, and also explains the IceCube neutrinos and the non-blazar component of the Fermi gamma-ray background, assuming a reasonable energy output from black hole jets in clusters.

SS 433 is a binary system containing a supergiant star that is overflowing its Roche lobe with matter accreting onto a compact object (either a black hole or neutron star). Two jets of ionized matter with a bulk velocity of $\sim0.26c$ extend from the binary, perpendicular to the line of sight, and terminate inside W50, a supernova remnant that is being distorted by the jets. SS 433 differs from other microquasars in that the accretion is believed to be super-Eddington, and the luminosity of the system is $\sim10^{40}$ erg s$^{-1}$. The lobes of W50 in which the jets terminate, about 40 pc from the central source, are expected to accelerate charged particles, and indeed radio and X-ray emission consistent with electron synchrotron emission in a magnetic field have been observed. At higher energies (>100 GeV), the particle fluxes of $\gamma$ rays from X-ray hotspots around SS 433 have been reported as flux upper limits. In this energy regime, it has been unclear whether the emission is dominated by electrons that are interacting with photons from the cosmic microwave background through inverse-Compton scattering or by protons interacting with the ambient gas. Here we report TeV $\gamma$-ray observations of the SS 433/W50 system where the lobes are spatially resolved. The TeV emission is localized to structures in the lobes, far from the center of the system where the jets are formed. We have measured photon energies of at least 25 TeV, and these are certainly not Doppler boosted, because of the viewing geometry. We conclude that the emission from radio to TeV energies is consistent with a single population of electrons with energies extending to at least hundreds of TeV in a magnetic field of $\sim16$~micro-Gauss.

Vela X is a nearby pulsar wind nebula (PWN) powered by a $\sim 10^4$ year old pulsar. Modeling of the spectral energy distribution of the Vela X PWN has shown that accelerated electrons have largely escaped from the confinement, which is likely due to the disruption of the initially confined PWN by the SNR reverse shock. The escaped electrons propagate to the earth and contribute to the measured local cosmic-ray (CR) electron spectrum. We find that the escaped CR electrons from Vela X would hugely exceed the measured flux by HESS at $\sim 10$ TeV if the standard diffusion coefficient for the interstellar medium is used. We propose that the diffusion may be highly inefficient around Vela X and find that a spatially-dependent diffusion can lead to CR flux consistent with the HESS measurement. Using a two-zone model for the diffusion around Vela X, we find that the diffusion coefficient in the inner region of a few tens of parsecs should be $<10^{28}{\rm cm^2 s^{-1}}$ for $\sim10$ TeV CR electrons, which is about two orders of magnitude lower than the standard value for ISM. Such inefficient diffusion around PWN resembles the case of the Geminga and Monogem PWNe, suggesting that inefficient diffusion may be common in the vicinity of PWNe spanning a wide range of ages.

The association of the high-energy neutrino event IceCube-170922A with the flaring blazar TXS~0506+056 indicates that hadronic processes may operate in a blazar jet. We perform semi-analytical spectral fitting of the multi-wavelength emission to obtain estimates of the jet physical parameters, and find that the multi-wavelength emission can be explained by either a proton synchrotron scenario or an electron inverse Compton scattering scenario. In the proton synchrotron scenario, a strong magnetic field of $10-100$~G is required, implying that the particle acceleration is likely driven by magnetic energy dissipation such as magnetic reconnection events. The inverse Compton scenario implies a magnetic field of $0.1-1$~G. Thus the particle acceleration is likely driven by the kinetic energy dissipation such as shocks. We also discuss the neutrino production in the context of single-zone and multi-zone models based on the above two scenarios. We demonstrate that the variability and optical polarization signatures can be used to distinguish the two scenarios due to their drastically different magnetic field. Specifically, the proton synchrotron scenario may show orphan fast variability in the low-energy spectral component on top of the active state, with an optical polarization degree $\lesssim 10\%$ throughout the active state. The inverse Compton scattering scenario instead predicts co-variability of the low- and high-energy components on both short and long time scales, as well as a strongly variable optical polarization degree that can reach $\gtrsim 20\%$. Our results suggest that optical polarization measurements and well-resolved multi-wavelength light curves can be used to understand the electromagnetic and high-energy neutrino emissions by TXS~0506+056 and similar events in the future.

The recent HAWC observations of very-high-energy $\gamma$-ray halo around Geminga and Monogem indicate a very slow diffusion of cosmic rays, which results in tiny contribution of positrons from these two pulsars to the local flux. This makes the cosmic positron excess anomaly observed by PAMELA and AMS-02 even more puzzling. However, from the Boron-to-Carbon ratio data one can infer that the average diffusion coefficient in the Galaxy should be much larger. In this work we propose a two-zone diffusion model that the diffusion is slow only in a small region around the source, outside of which the propagation is as fast as usual. We find that such a scenario can naturally explain the positron excess data with parameters even more reasonable than that in the conventional one-zone diffusion model. The reason is that during the life time of Geminga ($\sim 300$ kyr) the electrons/positrons have propagated too far away with a fast diffusion and lead to a low local flux. The slow diffusion region in the two-zone model helps to confine the electrons/positrons for a long time and lead to an enhancement of the local flux. So under the constraint of the HAWC observations, pulsars are still the probable origin of the cosmic-ray positron excess.

The Fermi-LAT detects no significant anisotropy of the cosmic-ray (CR) electrons and positrons ($e^-+e^+$) with seven years of data, which provides the strongest restriction to the $e^-+e^+$ anisotropy up to now. As next generation CR observatory, HERD is expected to have a better capability of anisotropy detection than Fermi-LAT. In this paper, we discuss several models aimed to explain the AMS-02 data by the present and future anisotropy measurements. We find that the upper limits of Fermi-LAT disfavor Vela SNR as the dominant source in sub-TeV, while other cases that remain safe under the constraint of Fermi-LAT are expected to be distinguished from each other by HERD. We then discuss the possibilities of remarkable TeV spectral features, and test the corresponding anisotropies. We find the conditions under which the TeV model can have a prominent spectral feature and avoid the constraint of Fermi-LAT at the same time. Furthermore, the expected performance of HERD is sensitive enough to detect the anisotropies of all these TeV models, and even for the case of a featureless TeV spectrum. Thus HERD may play a crucial part in the study of the origin of cosmic electrons and positrons.

A tentative excess in the electron spectrum at 1.4 TeV was recently reported by the DArk Matter Particle Explorer (DAMPE). A non-astrophysical scenario in which dark matter particles annihilate or decay in a local clump has been invoked to explain the excess. If $e^\pm$ annihilation channels in the final states are mediated by left-handed leptons as a component in the $SU(2)_L$ doublet, neutrinos with similar energies should have been simultaneously produced. We demonstrate that generic dark matter models can be decisively tested by the existing IceCube data. In case of a non-detection, such models would be excluded at the $5\sigma$ level by the five-year data for a point-like source and by the ten-year data for an extended source of dark matter particles with left-handed leptons.

The recent detection of gravitational waves and electromagnetic counterparts from the double neutron star merger event GW+EM170817, supports the standard paradigm of short gamma-ray bursts (SGRBs) and kilonovae/macronovae. It is important to reveal the nature of the compact remnant left after the merger, either a black hole or neutron star, and their physical link to the origin of the long-lasting emission observed in SGRBs. The diversity of the merger remnants may also lead to different kinds of transients that can be detected in future. Here we study the high-energy emission from the long-lasting central engine left after the coalescence, under certain assumptions. In particular, we consider the X-ray emission from a remnant disk and the non-thermal nebular emission from disk-driven outflows or pulsar winds. We demonstrate that late-time X-ray and high-frequency radio emission can provide useful constraints on properties of the hidden compact remnants and their connections to long-lasting SGRB emission, and we discuss the detectability of nearby merger events through late-time observations at $\sim30-100$ d after the coalescence. We also investigate the GeV-TeV gamma-ray emission that occurs in the presence of long-lasting central engines, and show the importance of external inverse-Compton radiation due to up-scattering of X-ray photons by relativistic electrons in the jet. We also search for high-energy gamma-rays from GW170817 in the Fermi-LAT data, and report upper limits on such long-lasting emission. Finally, we consider the implications of GW+EM170817 and discuss the constraints placed by X-ray and high-frequency radio observations.

The DArk Matter Particle Explorer (DAMPE), a space-based high precision cosmic ray detector, has just reported the new measurement of the total electron plus positron energy spectrum up to 4.6 TeV. A notable feature in the spectrum is the spectral break at $\sim 0.9$ TeV, with the spectral index softening from $-3.1$ to $-3.9$. Such a feature is very similar to the knee at the cosmic nuclei energy spectrum. In this work we propose that the knee-like feature can be explained naturally by assuming that the electrons are accelerated at the supernova remnants (SNRs) and released when the SNRs die out with lifetime around $10^5$ years. The cut-off energy of those electrons have already decreased to several TeV due to radiative cooling, which may induce the observed TeV spectral break. Another possibility is that the break is induced by a single nearby old SNR. Such a scenario may bring a large electron flux anisotropy that is observable by the future detectors. We also show that a minor part of electrons escaping during the acceleration in young and nearby SNRs are able to contribute to several TeV or higher energy region of the spectrum.

The DArk Matter Particle Explorer (DAMPE), a high energy cosmic ray and $\gamma$-ray detector in space, has recently reported the new measurement of the total electron plus positron flux between 25 GeV and 4.6 TeV. A spectral softening at $\sim0.9$ TeV and a tentative peak at $\sim1.4$ TeV have been reported. We study the physical implications of the DAMPE data in this work. The presence of the spectral break significantly tightens the constraints on the model parameters to explain the electron/positron excesses. The spectral softening can either be explained by the maximum acceleration limits of electrons by astrophysical sources, or a breakdown of the common assumption of continuous distribution of electron sources at TeV energies in space and time. The tentive peak at $\sim1.4$ TeV implies local sources of electrons/positrons with quasi-monochromatic injection spectrum. We find that the cold, ultra-relativistic $e^+e^-$ winds from pulsars may give rise to such a structure. The pulsar is requird to be middle-aged, relatively slowly-rotated, mildly magnetized, and isolated in a density cavity. The annihilation of DM particles ($m_{\chi}\sim1.5$ TeV) into $e^+e^-$ pairs in a nearby clump or an over-density region may also explain the data. In the DM scenario, the inferred clump mass (or density enhancement) is about $10^7-10^8$ M$_\odot$ (or $17-35$ times of the canonical local density) assuming a thermal production cross section, which is relatively extreme compared with the expectation from numerical simulations. A moderate enhancement of the annihilation cross section via, e.g., the Sommerfeld mechanism or non-thermal production, is thus needed.

The nature of the $\gamma$-ray emission from the \emph{Fermi} bubbles is unknown. Both hadronic and leptonic models have been formulated to explain the peculiar $\gamma$-ray signal observed by the Fermi-LAT between 0.1-500~GeV. If this emission continues above $\sim$30~TeV, hadronic models of the \emph{Fermi} bubbles would provide a significant contribution to the high-energy neutrino flux detected by the IceCube observatory. Even in models where leptonic $\gamma$-rays produce the \emph{Fermi} bubbles flux at GeV energies, a hadronic component may be observable at very high energies. The combination of IceCube and HAWC measurements have the ability to distinguish these scenarios through a comparison of the neutrino and $\gamma$-ray fluxes at a similar energy scale. We examine the most recent four-year dataset produced by the IceCube collaboration and find no evidence for neutrino emission originating from the \emph{Fermi} bubbles. In particular, we find that previously suggested excesses are consistent with the diffuse astrophysical background with a p-value of 0.22 (0.05 in an extreme scenario that all the IceCube events that overlap with the bubbles come from them). Moreover, we show that existing and upcoming HAWC observations provide independent constraints on any neutrino emission from the \emph{Fermi} bubbles, due to the close correlation between the $\gamma$-ray and neutrino fluxes in hadronic interactions. The combination of these results disfavors a significant contribution from the \emph{Fermi} bubbles to the IceCube neutrino flux.

AMIGA (Auger Muons and Infill for the Ground Array) is an upgrade of the Pierre Auger Observatory designed to extend its energy range of detection and to directly measure the muon content of the cosmic ray primary particle showers. The array will be formed by an infill of surface water-Cherenkov detectors associated with buried scintillation counters employed for muon counting. Each counter is composed of three scintillation modules, with a 10 m$^2$ detection area per module. In this paper, a new generation of detectors, replacing the current multi-pixel photomultiplier tube (PMT) with silicon photo sensors (aka. SiPMs), is proposed. The selection of the new device and its front-end electronics is explained. A method to calibrate the counting system that ensures the performance of the detector is detailed. This method has the advantage of being able to be carried out in a remote place such as the one where the detectors are deployed. High efficiency results, i.e. 98 % efficiency for the highest tested overvoltage, combined with a low probability of accidental counting ($\sim$2 %), show a promising performance for this new system.

The merger of a neutron star (NS) binary may result in the formation of a long-lived, or indefinitely stable, millisecond magnetar remnant surrounded by a low-mass ejecta shell. A portion of the magnetar's prodigious rotational energy is deposited behind the ejecta in a pulsar wind nebula, powering luminous optical/X-ray emission for hours to days following the merger. Ions in the pulsar wind may also be accelerated to ultra-high energies, providing a coincident source of high energy cosmic rays and neutrinos. At early times, the cosmic rays experience strong synchrotron losses; however, after a day or so, pion production through photomeson interaction with thermal photons in the nebula comes to dominate, leading to efficient production of high-energy neutrinos. After roughly a week, the density of background photons decreases sufficiently for cosmic rays to escape the source without secondary production. These competing effects result in a neutrino light curve that peaks on a few day timescale near an energy of $\sim10^{18}$ eV. This signal may be detectable for individual mergers out to $\sim$ 10 (100) Mpc by current (next-generation) neutrino telescopes, providing clear evidence for a long-lived NS remnant, the presence of which may otherwise be challenging to identify from the gravitational waves alone. Under the optimistic assumption that a sizable fraction of NS mergers produce long-lived magnetars, the cumulative cosmological neutrino background is estimated to be $\sim 10^{-9}-10^{-8}\,\rm GeV\,cm^{-2}\,s^{-1}\,sr^{-1}$ for a NS merger rate of $10^{-7}\,\rm Mpc^{-3}\,yr^{-1}$, overlapping with IceCube's current sensitivity and within the reach of next-generation neutrino telescopes.

The Giant Radio Array for Neutrino Detection (GRAND) aims at detecting ultra-high energy extraterrestrial neutrinos via the extensive air showers induced by the decay of tau leptons created in the interaction of neutrinos under the Earth's surface. Consisting of an array of $\sim10^5$ radio antennas deployed over $\sim 2\times10^5\,\rm {km}^2$, GRAND plans to reach, for the first time, an all-flavor sensitivity of $\sim1.5\times10^{-10} \,\rm GeV\, cm^{-2} \,s^{-1}\, sr^{-1}$ above $5\times10^{17}$ eV and a sub-degree angular resolution, beyond the reach of other planned detectors. We describe here preliminary designs and simulation results, plans for the ongoing, staged approach to the construction of GRAND, and the rich research program made possible by GRAND's design sensitivity and angular resolution.

In this work we use the newly reported Boron-to-Carbon ratio (B/C) from AMS-02 and the time-dependent proton fluxes from PAMELA and AMS-02 to constrain the source and propagation parameters of cosmic rays in the Milky Way. A linear correlation of the solar modulation parameter with solar activities is assumed to account for the time-varying cosmic ray fluxes. A comprehensive set of propagation models, with/without reacceleration or convection, have been discussed and compared. We find that only the models with reacceleration can self-consistently fit both the proton and B/C data. The rigidity dependence slope of the diffusion coefficient, $\delta$, is found to be about $0.38-0.50$ for the diffusion-reacceleration models. The plain diffusion and diffusion-convection models fit the data poorly. We compare different model predictions of the positron and antiproton fluxes with the data. We find that the diffusion-reacceleration models over-produce low energy positrons, while non-reacceleration models give better fit to the data. As for antiprotons, reacceleration models tend to under-predict low energy antiproton fluxes, unless a phenomenological modification of the velocity-dependence of the diffusion coefficient is applied. Our results suggest that there could be important differences of the propagation for nuclei and leptons, in either the Milky Way or the solar heliosphere.

Powerful jets and outflows generated in tidal disruption events (TDEs) around supermassive black holes have been suggested to be possible sites producing high-energy neutrinos, but it is unclear whether such environment can provide the bulk of the neutrinos detected by the IceCube Observatory. In this work, by considering realistic limits on the non-thermal emission power of a TDE jet and the birth rate of the TDEs with jets pointing towards us, we show that it is hard to use the jetted TDE population to explain the flux and the isotropic arrival directions of the observed TeV-PeV neutrinos. Therefore, TDEs cannot be the dominant sources, unless those without aligned jets can produce wide-angle emission of neutrino particles. Supposing that is the case, we list a few recent jetted and non-jetted TDEs that have the best chance to be detected by IceCube, based on their energetics, distances, and directions. A spatial and temporal association of these predicted events with the IceCube data should provide a decisive test on the full TDE population as origins of the IceCube neutrinos.

Observations by HAWC and Milagro have detected bright and spatially extended TeV gamma-ray sources surrounding the Geminga and Monogem pulsars. We argue that these observations, along with a substantial population of other extended TeV sources coincident with pulsar wind nebulae, constitute a new morphological class of spatially extended TeV halos. We show that HAWCs wide field-of-view unlocks an expansive parameter space of TeV halos not observable by atmospheric Cherenkov telescopes. Under the assumption that Geminga and Monogem are typical middle-aged pulsars, we show that ten-year HAWC observations should eventually observe 37$^{+17}_{-13}$ middle-aged TeV halos that correspond to pulsars whose radio emission is not beamed towards Earth. Depending on the extrapolation of the TeV halo efficiency to young pulsars, HAWC could detect more than 100 TeV halos from mis-aligned pulsars. These pulsars have historically been difficult to detect with existing multiwavelength observations. TeV halos will constitute a significant fraction of all HAWC sources, allowing follow-up observations to efficiently find pulsar wind nebulae and thermal pulsar emission. The observation and subsequent multi-wavelength follow-up of TeV halos will have significant implications for our understanding of pulsar beam geometries, the evolution of PWN, the diffusion of cosmic-rays near energetic pulsars, and the contribution of pulsars to the cosmic-ray positron excess.

Recent measurements of the Geminga and B0656+14 pulsars by the gamma-ray telescope HAWC (along with earlier measurements by Milagro) indicate that these objects generate significant fluxes of very high-energy electrons. In this paper, we use the very high-energy gamma-ray intensity and spectrum of these pulsars to calculate and constrain their expected contributions to the local cosmic-ray positron spectrum. Among models that are capable of reproducing the observed characteristics of the gamma-ray emission, we find that pulsars invariably produce a flux of high-energy positrons that is similar in spectrum and magnitude to the positron fraction measured by PAMELA and AMS-02. In light of this result, we conclude that it is very likely that pulsars provide the dominant contribution to the long perplexing cosmic-ray positron excess.

The Giant Radio Array for Neutrino Detection (GRAND) is a planned array of ~200 000 radio antennas deployed over ~200 000 km2 in a mountainous site. It aims primarly at detecting high-energy neutrinos via the observation of extensive air showers induced by the decay in the atmosphere of taus produced by the interaction of cosmic neutrinos under the Earth surface. GRAND aims at reaching a neutrino sensitivity of 5.10$^{11}$ E$^{-2}$ GeV$^{-1}$cm$^{-2}$s$^{-1}$sr$^{-1}$ above 3.10$^{16}$ eV. This ensures the detection of cosmogenic neutrinos in the most pessimistic source models, and ~50 events per year are expected for the standard models. The instrument will also detect UHECRs and possibly FRBs. Here we show how our preliminary design should enable us to reach our sensitivity goals, and discuss the steps to be taken to achieve GRAND.

On September 14, 2015 the Advanced LIGO detectors observed their first gravitational-wave (GW) transient GW150914. This was followed by a second GW event observed on December 26, 2015. Both events were inferred to have arisen from the merger of black holes in binary systems. Such a system may emit neutrinos if there are magnetic fields and disk debris remaining from the formation of the two black holes. With the surface detector array of the Pierre Auger Observatory we can search for neutrinos with energy above 100 PeV from point-like sources across the sky with equatorial declination from about -65 deg. to +60 deg., and in particular from a fraction of the 90% confidence-level (CL) inferred positions in the sky of GW150914 and GW151226. A targeted search for highly-inclined extensive air showers, produced either by interactions of downward-going neutrinos of all flavors in the atmosphere or by the decays of tau leptons originating from tau-neutrino interactions in the Earth's crust (Earth-skimming neutrinos), yielded no candidates in the Auger data collected within $\pm 500$ s around or 1 day after the coordinated universal time (UTC) of GW150914 and GW151226, as well as in the same search periods relative to the UTC time of the GW candidate event LVT151012. From the non-observation we constrain the amount of energy radiated in ultrahigh-energy neutrinos from such remarkable events.

The AMS-02 has measured the cosmic ray electron (plus positron) spectrum up to ~TeV with an unprecedent precision. The spectrum can be well described by a power law without any obvious features above 10 GeV. The satellite instrument Dark Matter Particle Explorer (DAMPE), which was launched a year ago, will measure the electron spectrum up to 10 TeV with a high energy resolution. The cosmic electrons beyond TeV may be attributed to few local cosmic ray sources, such as supernova remnants. Therefore, spectral features, such as cutoff and bumps, can be expected at high energies. In this work we give a careful study on the perspective of the electron spectrum beyond TeV. We first examine our astrophysical source models on the latest leptonic data of AMS-02 to give a self-consistent picture. Then we focus on the discussion about the candidate sources which could be electron contributors above TeV. Depending on the properties of the local sources (especially on the nature of Vela), DAMPE may detect interesting features in the electron spectrum above TeV in the future.

We report a first measurement for ultra-high energy cosmic rays of the correlation between the depth of shower maximum and the signal in the water Cherenkov stations of air-showers registered simultaneously by the fluorescence and the surface detectors of the Pierre Auger Observatory. Such a correlation measurement is a unique feature of a hybrid air-shower observatory with sensitivity to both the electromagnetic and muonic components. It allows an accurate determination of the spread of primary masses in the cosmic-ray flux. Up till now, constraints on the spread of primary masses have been dominated by systematic uncertainties. The present correlation measurement is not affected by systematics in the measurement of the depth of shower maximum or the signal in the water Cherenkov stations. The analysis relies on general characteristics of air showers and is thus robust also with respect to uncertainties in hadronic event generators. The observed correlation in the energy range around the `ankle' at $\lg(E/{\rm eV})=18.5-19.0$ differs significantly from expectations for pure primary cosmic-ray compositions. A light composition made up of proton and helium only is equally inconsistent with observations. The data are explained well by a mixed composition including nuclei with mass $A > 4$. Scenarios such as the proton dip model, with almost pure compositions, are thus disfavoured as the sole explanation of the ultrahigh-energy cosmic-ray flux at Earth.

The detection of ultrahigh-energy (UHE) neutrino sources would contribute significantly to solving the decades-old mystery of the origin of the highest-energy cosmic rays. We investigate the ability of a future UHE neutrino detector to identify the brightest neutrino point sources, by exploring the parameter space of the total number of observed events and the angular resolution of the detector. The favored parameter region can be translated to requirements for the effective area, sky coverage and angular resolution of future detectors, for a given source number density and evolution history. Moreover, by studying the typical distance to sources that are expected to emit more than one event for a given diffuse neutrino flux, we find that a significant fraction of the identifiable UHE neutrino sources may be located in the nearby Universe if the source number density is above $\sim10^{-6}\,\rm Mpc^{-3}$. If sources are powerful and rare enough, as predicted in blazar scenarios, they can first be detected at distant locations. Our result also suggests that if UHE cosmic-ray accelerators are neither beamed nor transients, it will be possible to associate the detected UHE neutrino sources with nearby UHE cosmic-ray and gamma-ray sources, and that they may also be observed using other messengers, including ones with limited horizons such as TeV gamma rays, UHE gamma rays and cosmic rays. We find that for a $\gtrsim5\sigma$ detection of UHE neutrino sources with a uniform density, $n_s\sim{10}^{-7}-{10}^{-5}~{\rm Mpc}^{-3}$, at least $\sim100-1000$ events and sub-degree angular resolution are needed, and the results depend on the source evolution model.

Ultrahigh energy cosmic ray air showers probe particle physics at energies beyond the reach of accelerators. Here we introduce a new method to test hadronic interaction models without relying on the absolute energy calibration, and apply it to events with primary energy 6-16 EeV (E_CM = 110-170 TeV), whose longitudinal development and lateral distribution were simultaneously measured by the Pierre Auger Observatory. The average hadronic shower is 1.33 +- 0.16 (1.61 +- 0.21) times larger than predicted using the leading LHC-tuned models EPOS-LHC (QGSJetII-04), with a corresponding excess of muons.

The bulk of observed ultrahigh energy cosmic rays could be light or heavier elements, and originate from an either steady or transient population of sources. This leaves us with four general categories of sources. Energetic requirements set a lower limit on single source luminosities, while the distribution of particle arrival directions in the sky sets a lower limit on the source number density. The latter constraint depends on the angular smearing in the skymap due to the magnetic deflections of the charged particles during their propagation from the source to the Earth. We contrast these limits with the luminosity functions from surveys of existing luminous steady objects in the nearby universe, and strongly constrain one of the four categories of source models, namely, steady proton sources. The possibility that cosmic rays with energy $> 8\times 10^{19}\,$eV are dominantly pure protons coming from steady sources is excluded at 95\% confidence level, under the safe assumption that protons experience less than $30^\circ$ magnetic deflection on flight.

Radio observations at multiple frequencies have detected a significant isotropic emission component between 22~MHz and 10~GHz, commonly termed the ARCADE-2 Excess. The origin of this radio emission is unknown, as the intensity, spectrum and isotropy of the signal are difficult to model with either traditional astrophysical mechanisms or novel physics such as dark matter annihilation. We posit a new model capable of explaining the key components of the excess radio emission. Specifically, we show that the re-acceleration of non-thermal electrons via turbulence in merging galaxy clusters are capable of explaining the intensity, spectrum, and isotropy of the ARCADE-2 data. We examine the parameter spaces of cluster re-acceleration, magnetic field, and merger rate, finding that the radio excess can be reproduced assuming reasonable assumptions for each. Finally, we point out that future observations will definitively confirm or rule-out the contribution of cluster mergers to the isotropic radio background.

We present a search for ultra-relativistic magnetic monopoles with the Pierre Auger Observatory. Such particles, possibly a relic of phase transitions in the early universe, would deposit a large amount of energy along their path through the atmosphere, comparable to that of ultrahigh-energy cosmic rays (UHECRs). The air shower profile of a magnetic monopole can be effectively distinguished by the fluorescence detector from that of standard UHECRs. No candidate was found in the data collected between 2004 and 2012, with an expected background of less than 0.1 event from UHECRs. The corresponding 90% confidence level (C.L.) upper limits on the flux of ultra-relativistic magnetic monopoles range from $10^{-19}$ (cm$^{2}$ sr s)$^{-1}$ for a Lorentz factor $\gamma=10^9$ to $2.5 \times10^{-21}$ (cm$^{2}$ sr s)$^{-1}$ for $\gamma=10^{12}$. These results - the first obtained with a UHECR detector - improve previously published limits by up to an order of magnitude.

The IceCube collaboration has reported the first detection of high-energy astrophysical neutrinos including $\sim 50$ high-energy starting events, but no individual sources have been identified. It is therefore important to develop the most sensitive and efficient possible algorithms to identify point sources of these neutrinos. The most popular current method works by exploring a dense grid of possible directions to individual sources, and identifying the single direction with the maximum probability of having produced multiple detected neutrinos. This method has numerous strengths, but it is computationally intensive and, because it focuses on the single best location for a point source, additional point sources are not included in the evidence. We propose a new maximum likelihood method that uses the angular separations between all pairs of neutrinos in the data. Unlike existing autocorrelation methods for this type of analysis, which also use angular separations between neutrino pairs, our method incorporates information about the point spread function and can identify individual point sources. We find that if the angular resolution is a few degrees or better, then this approach reduces both false positive and false negative errors compared to the current method, and is also more computationally efficient up to, potentially, hundreds of thousands of detected neutrinos.

High-energy cosmic rays can be accelerated in clusters of galaxies, by mega-parsec scale shocks induced by accretion of gas during the formation of large-scale structure, or by powerful sources harbored in clusters. Once accelerated, the highest energy particles leave the cluster via almost rectilinear trajectories, while lower energy ones can be confined by the cluster magnetic field up to cosmological time and interact with the intracluster gas. Using a realistic model of the baryon distribution and the turbulent magnetic field in clusters, we studied the propagation and hadronic interaction of high-energy protons in the intracluster medium. We report the cumulative cosmic ray and neutrino spectra generated by galaxy clusters including embedded sources, and demonstrate that clusters can contribute a significant fraction of the observed IceCube neutrinos above 30 TeV while remaining undetected in high-energy cosmic rays and $\gamma$ rays for reasonable choices of parameters and source scenarios.

We measure the energy emitted by extensive air showers in the form of radio emission in the frequency range from 30 to 80 MHz. Exploiting the accurate energy scale of the Pierre Auger Observatory, we obtain a radiation energy of 15.8 \pm 0.7 (stat) \pm 6.7 (sys) MeV for cosmic rays with an energy of 1 EeV arriving perpendicularly to a geomagnetic field of 0.24 G, scaling quadratically with the cosmic-ray energy. A comparison with predictions from state-of-the-art first-principle calculations shows agreement with our measurement. The radiation energy provides direct access to the calorimetric energy in the electromagnetic cascade of extensive air showers. Comparison with our result thus allows the direct calibration of any cosmic-ray radio detector against the well-established energy scale of the Pierre Auger Observatory.

The Auger Engineering Radio Array (AERA) is part of the Pierre Auger Observatory and is used to detect the radio emission of cosmic-ray air showers. These observations are compared to the data of the surface detector stations of the Observatory, which provide well-calibrated information on the cosmic-ray energies and arrival directions. The response of the radio stations in the 30 to 80 MHz regime has been thoroughly calibrated to enable the reconstruction of the incoming electric field. For the latter, the energy deposit per area is determined from the radio pulses at each observer position and is interpolated using a two-dimensional function that takes into account signal asymmetries due to interference between the geomagnetic and charge-excess emission components. The spatial integral over the signal distribution gives a direct measurement of the energy transferred from the primary cosmic ray into radio emission in the AERA frequency range. We measure 15.8 MeV of radiation energy for a 1 EeV air shower arriving perpendicularly to the geomagnetic field. This radiation energy -- corrected for geometrical effects -- is used as a cosmic-ray energy estimator. Performing an absolute energy calibration against the surface-detector information, we observe that this radio-energy estimator scales quadratically with the cosmic-ray energy as expected for coherent emission. We find an energy resolution of the radio reconstruction of 22% for the data set and 17% for a high-quality subset containing only events with at least five radio stations with signal.

The azimuthal asymmetry in the risetime of signals in Auger surface detector stations is a source of information on shower development. The azimuthal asymmetry is due to a combination of the longitudinal evolution of the shower and geometrical effects related to the angles of incidence of the particles into the detectors. The magnitude of the effect depends upon the zenith angle and state of development of the shower and thus provides a novel observable, $(\sec \theta)_\mathrm{max}$, sensitive to the mass composition of cosmic rays above $3 \times 10^{18}$ eV. By comparing measurements with predictions from shower simulations, we find for both of our adopted models of hadronic physics (QGSJETII-04 and EPOS-LHC) an indication that the mean cosmic-ray mass increases slowly with energy, as has been inferred from other studies. However, the mass estimates are dependent on the shower model and on the range of distance from the shower core selected. Thus the method has uncovered further deficiencies in our understanding of shower modelling that must be resolved before the mass composition can be inferred from $(\sec \theta)_\mathrm{max}$.

The Pierre Auger Observatory has begun a major Upgrade of its already impressive capabilities, with an emphasis on improved mass composition determination using the surface detectors of the Observatory. Known as AugerPrime, the upgrade will include new 4 m$^2$ plastic scintillator detectors on top of all 1660 water-Cherenkov detectors, updated and more flexible surface detector electronics, a large array of buried muon detectors, and an extended duty cycle for operations of the fluorescence detectors. This Preliminary Design Report was produced by the Collaboration in April 2015 as an internal document and information for funding agencies. It outlines the scientific and technical case for AugerPrime. We now release it to the public via the arXiv server. We invite you to review the large number of fundamental results already achieved by the Observatory and our plans for the future.

Relativistic winds of fast-spinning pulsars have been proposed as a potential site for cosmic-ray acceleration from very high energies (VHE) to ultrahigh energies (UHE). We re-examine conditions for high-energy neutrino production, considering the interaction of accelerated particles with baryons of the expanding supernova ejecta and the radiation fields in the wind nebula. We make use of the current IceCube sensitivity in diffusive high-energy neutrino background, in order to constrain the parameter space of the most extreme neutron stars as sources of VHE and UHE cosmic rays. We demonstrate that the current non-observation of $10^{18}$ eV neutrinos put stringent constraints on the pulsar scenario. For a given model, birthrates, ejecta mass and acceleration efficiency of the magnetar sources can be constrained. When we assume a proton cosmic ray composition and spherical supernovae ejecta, we find that the IceCube limits almost exclude their significant contribution to the observed UHE cosmic-ray flux. Furthermore, we consider scenarios where a fraction of cosmic rays can escape from jet-like structures piercing the ejecta, without significant interactions. Such scenarios would enable the production of UHE cosmic rays and help remove the tension between their EeV neutrino production and the observational data.

Fast-spinning newborn pulsars are intriguing candidate sources of ultrahigh energy cosmic rays (UHECRs). The acceleration of particles with a given composition in a fraction of the extragalactic pulsar population can give a consistent explanation for the measurements of the Auger Observatory. We calculate the associated diffuse neutrino flux produced while particles cross the supernova ejecta surrounding the stars. We show that in the minimal pulsar scenario that are compatible with the UHECR data, the effective optical depth to hadronuclear interactions is larger than unity at ultrahigh energies. Thus, even in the most pessimistic case, one expects energy fluxes of $\sim 0.1 - 1 \,\rm EeV$ neutrinos that should be detectable with IceCube or Askaryan Radio Array within a decade.

This paper presents the results of different searches for correlations between very high-energy neutrino candidates detected by IceCube and the highest-energy cosmic rays measured by the Pierre Auger Observatory and the Telescope Array. We first consider samples of cascade neutrino events and of high-energy neutrino-induced muon tracks, which provided evidence for a neutrino flux of astrophysical origin, and study their cross-correlation with the ultrahigh-energy cosmic ray (UHECR) samples as a function of angular separation. We also study their possible directional correlations using a likelihood method stacking the neutrino arrival directions and adopting different assumptions on the size of the UHECR magnetic deflections. Finally, we perform another likelihood analysis stacking the UHECR directions and using a sample of through-going muon tracks optimized for neutrino point-source searches with sub-degree angular resolution. No indications of correlations at discovery level are obtained for any of the searches performed. The smallest of the p-values comes from the search for correlation between UHECRs with IceCube high-energy cascades, a result that should continue to be monitored.

A measurement of the cosmic-ray spectrum for energies exceeding $4{\times}10^{18}$ eV is presented, which is based on the analysis of showers with zenith angles greater than $60^{\circ}$ detected with the Pierre Auger Observatory between 1 January 2004 and 31 December 2013. The measured spectrum confirms a flux suppression at the highest energies. Above $5.3{\times}10^{18}$ eV, the "ankle", the flux can be described by a power law $E^{-\gamma}$ with index $\gamma=2.70 \pm 0.02 \,\text{(stat)} \pm 0.1\,\text{(sys)}$ followed by a smooth suppression region. For the energy ($E_\text{s}$) at which the spectral flux has fallen to one-half of its extrapolated value in the absence of suppression, we find $E_\text{s}=(5.12\pm0.25\,\text{(stat)}^{+1.0}_{-1.2}\,\text{(sys)}){\times}10^{19}$ eV.

Neutrinos in the cosmic ray flux with energies near 1 EeV and above are detectable with the Surface Detector array of the Pierre Auger Observatory. We report here on searches through Auger data from 1 January 2004 until 20 June 2013. No neutrino candidates were found, yielding a limit to the diffuse flux of ultra-high energy neutrinos that challenges the Waxman-Bahcall bound predictions. Neutrino identification is attempted using the broad time-structure of the signals expected in the SD stations, and is efficiently done for neutrinos of all flavors interacting in the atmosphere at large zenith angles, as well as for "Earth-skimming" neutrino interactions in the case of tau neutrinos. In this paper the searches for downward-going neutrinos in the zenith angle bins $60^\circ-75^\circ$ and $75^\circ-90^\circ$ as well as for upward-going neutrinos, are combined to give a single limit. The $90\%$ C.L. single-flavor limit to the diffuse flux of ultra-high energy neutrinos with an $E^{-2}$ spectrum in the energy range $1.0 \times 10^{17}$ eV - $2.5 \times 10^{19}$ eV is $E_\nu^2 dN_\nu/dE_\nu < 6.4 \times 10^{-9}~ {\rm GeV~ cm^{-2}~ s^{-1}~ sr^{-1}}$.

Joint contributions of the Pierre Auger Collaboration and the Telescope Array Collaboration to the 34th International Cosmic Ray Conference, 30 July - 6 August 2015, The Hague, The Netherlands.