Neutrinos of astrophysical origin could be detected through the electromagnetic radiation of the particle showers induced in the atmosphere by their interaction in the Earth. This applies in particular for tau neutrinos of energies $E>10^{16}$eV following Earth-skimming trajectories. The $\sim1^{\circ}$ beaming of the radio emission in the forward direction however implies that the radio signal will likely fly above a detector deployed over a flat site and would therefore not be detected. We study here how a non-flat detector topography can improve the detection probability of these neutrino-induced air showers. We do this by computing with three distinct tools the neutrino detection efficiency for a radio array deployed over a toy-model mountainous terrain, also taking into account experimental and topographic constraints. We show in particular that ground topographies inclined by few degrees only induce detection efficiencies typically three times larger than those obtained for flat areas. We conclude that the topography of the area where the detector is deployed will be a key factor for an experiment like GRAND.

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.

Pulsars have been identified as good candidates for the acceleration of cosmic rays, up to ultra-high energies. However, a precise description of the acceleration processes at play is still to be established. Using 2D particle-in-cell simulations, we study proton acceleration in axisymmetric pulsar magnetospheres. Protons and electrons are extracted from the neutron star surface by the strong electric field induced by the rotation of the star, and electrons and positrons are produced in the magnetosphere through pair production process. As pair production has a crucial impact on electromagnetic fields, on gaps and thus on particle acceleration, we study its influence on the maximum energy and luminosity of protons escaping the magnetosphere. Protons are accelerated and escape in all our simulations. However, the acceleration sites are different for the protons and the pairs. As shown in previous studies, pairs are accelerated to their highest energies at the Y-point and in the equatorial current sheet, where magnetic reconnection plays and important role. In contrast, protons gain most of their kinetic energy below the light-cylinder radius within the separatrix current layers, but they are not confined within the equatorial current sheet. Their maximum Lorentz factors can reach $15\%$ to $75\%$ of the maximum Lorentz factor obtained by acceleration through the full vacuum potential drop from pole to equator, and increase with decreasing pair production. Their luminosity can reach $0.2\%$ to $4\%$ of the theoretical spin down luminosity of an aligned pulsar, and the minimum luminosity is obtained at the transition between the force-free and electrosphere regimes. These estimates support that millisecond pulsars could accelerate cosmic rays up to PeV energies and that new born millisecond pulsars could accelerate cosmic rays up to ultra-high energies.

Binary systems composed of a recycled millisecond pulsar and a stellar companion in close orbit could be excellent sites to diagnose pulsar winds. In such systems, the pulsar outflow irradiates and heats up the companion atmosphere, which can lead to the observation of strong day/night modulations in temperature. We demonstrate with particle shower simulations that the particle energy of the wind affects the heating depth in the atmosphere: the wind heat can be deposited above or below the photosphere, leading to different signatures in the observed spectra. We apply our method to four specific systems: We find that systems with cool night side companions showing strong temperature variations can give interesting lower limits on the particle energy in the winds. In particular, if the companion night side of PSR B1957+20 were to be suddenly irradiated, deep heating would only take place if particles with energy > 100 TeV were present. Observational evidence of deep heating in this system thus suggests that i) such particles exist in the pulsar wind and/or ii) binary evolution non-trivially takes the companion to the observed temperature asymmetry. Besides, the observed temperature difference can be maintained only with particle energies of the order of 100 MeV.

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.

Over the last years, radio detection has matured to become a competitive method for the detection of air showers. Arrays of thousands of antennas are now envisioned for the detection of cosmic rays of ultra high energy or neutrinos of astrophysical origin. The data exploitation of such detectors requires to run massive air-shower simulations to evaluate the radio signal at each antenna position. In order to reduce the associated computational cost, we have developed a semi-analytical method for the computation of the emitted radio signal called Radio Morphing. The method consists in computing the radio signal of any air-shower at any location from the simulation of one single reference shower at given positions by i) a scaling of the electric-field amplitude of this reference shower, ii) an isometry on the simulated positions and iii) an interpolation of the radio pulse at the desired position. This technique enables one to compute electric field time traces with characteristics very similar to those obtained with standard computation methods, but with computation times reduced by several orders of magnitude. In this paper, we present this novel tool, explain its methodology, and discuss its limitations. Furthermore, we validate the method on a typical event set for the future GRAND experiment showing that the calculated peak amplitudes are consistent with the results from ZHAireS simulations with a mean offset of +8.5% and a standard deviation of 27.2% in this specific case. This overestimation of the signal strength by Radio Morphing arises mainly from the choice of the underlying reference shower.

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.

High-energy cosmic neutrinos carry unique information about the most energetic non-thermal sources in the Universe. This white paper describes the outstanding astrophysics questions that neutrino astronomy can address in the coming decade. A companion white paper discusses how the observation of cosmic neutrinos can address open questions in fundamental physics. Detailed measurements of the diffuse neutrino flux, measurements of neutrinos from point sources, and multi-messenger observations with neutrinos will enable the discovery and characterization of the most energetic sources in the Universe.

High-energy cosmic neutrinos can reveal new fundamental particles and interactions, probing energy and distance scales far exceeding those accessible in the laboratory. This white paper describes the outstanding particle physics questions that high-energy cosmic neutrinos can address in the coming decade. A companion white paper discusses how the observation of cosmic neutrinos can address open questions in astrophysics. Tests of fundamental physics using high-energy cosmic neutrinos will be enabled by detailed measurements of their energy spectrum, arrival directions, flavor composition, and timing.

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.

The interaction of ultra-high-energy cosmic rays (UHECRs) with pervasive photon fields generates associated cosmogenic fluxes of neutrinos and photons due to photohadronic and photonuclear processes taking place in the intergalactic medium. We perform a fit of the UHECR spectrum and composition measured by the Pierre Auger Observatory for four source emissivity scenarios: power-law redshift dependence with one free parameter, active galactic nuclei, gamma-ray bursts, and star formation history. We show that negative source emissivity evolution is favoured if we treat the source evolution as a free parameter. In all cases, the best fit is obtained for relatively hard spectral indices and low maximal rigidities, for compositions at injection dominated by intermediate nuclei (nitrogen and silicon groups). In light of these results, we calculate the associated fluxes of neutrinos and photons. Finally, we discuss the prospects for the future generation of high-energy neutrino and gamma-ray observatories to constrain the sources of UHECRs.

Over the last decades, radio detection of air showers has been established as a promising detection technique for ultrahigh-energy cosmic rays and neutrinos. Very large or dense antenna arrays are necessary to be proficient at collecting information about these particles and understanding their properties accurately. The exploitation of such arrays requires to run massive air-shower simulations to evaluate the radio signal at each antenna position, taking into account features such as the ground topology. In order to reduce computational costs, we have developed a fast computation of the emitted radio signal on the basis of generic shower simulations, called Radio Morphing. The method consists in the calculation of the radio signal of any air-shower by i) a scaling of the electric-field amplitude of a reference air shower to the target shower, ii) an isometry on the simulated positions and iii) an interpolation of the radio pulse at the desired position. This technique enables one to gain many orders of magnitude in CPU time compared to a standard computation. In this contribution, we present this novel tool and explain its methodology. In particular, Radio Morphing will be a key element for the simulation chain of the Giant Radio Array for Neutrino Detection (GRAND) project, that aims at detecting ultra-high-energy neutrinos with an array of 200 000 radio antennas in mountainous regions.os with an array of 200 000 radio antennas in mountainous regions.

Over the last decades, radio detection of air showers has been established as a detection technique for ultra-high-energy cosmic-rays impinging on the Earth's atmosphere with energies far beyond LHC energies. Today's second-generation of digital radio-detection experiments, as e.g. AERA or LOFAR, are becoming competitive in comparison to already standard techniques e.g. fluorescence light detection. Thanks to a detailed understanding of the physics of the radio emission in extensive air showers, simulations of the radio signal are already successfully tested and applied in the reconstruction of cosmic rays. However, the limits of the computational power resources are easily reached when it comes to computing electric fields at the numerous positions requested by large or dense antenna arrays. In the case of mountainous areas as e.g. for the GRAND array, where 3D shower simulations are necessary, the problem arises with even stronger acuity. Therefore we developed a full parametrisation of the emitted radio signal on the basis of generic shower simulations which will reduce the simulation time by orders of magnitudes. In this talk we will present this concept after a short introduction to the concept of the radio detection of air-shower induced by cosmic rays.

In this white paper we introduce the IMAGINE Consortium and its scientific background, goals and structure. Our purpose is to coordinate and facilitate the efforts of a diverse group of researchers in the broad areas of the interstellar medium, Galactic magnetic fields and cosmic rays, and our goal is to develop more comprehensive insights into the structures and roles of interstellar magnetic fields and their interactions with cosmic rays. To achieve a higher level of self-consistency, depth and rigour can only be achieved by the coordinated efforts of experts in diverse areas of astrophysics involved in observational, theoretical and numerical work. We present our view of the present status of this topic, identify its key unsolved problems and suggest a strategy that will underpin our work. The backbone of the consortium is the Interstellar MAGnetic field INference Engine, a publicly available Bayesian platform that employs robust statistical methods to explore the multi-dimensional likelihood space using any number of modular inputs. It provides an interpretation and modelling framework that has the power and flexibility to include a variety of observational, theoretical and numerical lines of evidence into a self-consistent and comprehensive picture of the thermal and non-thermal interstellar media. An important innovation is that a consistent understanding of the phenomena that are directly or indirectly influenced by the Galactic magnetic field, such as the deflection of ultra-high energy cosmic rays or extragalactic backgrounds, is made an integral part of the modelling. The IMAGINE Consortium, which is informal by nature and open to new participants, hereby presents a methodological framework for the modelling and understanding of Galactic magnetic fields that is available to all communities whose research relies on a state-of-the-art solution to this problem. (Abridged.)

Diffuse TeV emission has been observed by H.E.S.S. in the Galactic Center region, in addition to the GeV gamma rays observed by Fermi. We propose that a population of unresolved millisecond pulsars located around the Galactic Center, suggested as possible candidates for the diffuse Galactic Center excess observed by Fermi, accelerate cosmic rays up to very high energies, and are thus also responsible for the TeV excess. We model analytically the diffusion of these accelerated protons and their interaction with the molecular clouds, producing gamma rays. The spatial and spectral dependences of the gamma rays produced can reproduce the H.E.S.S. observations, for a population of $\sim 10^4-10^5$ millisecond pulsars above the cosmic-ray luminosity $10^{34}\,{\rm erg\,s}^{-1}$, with moderate acceleration efficiency. More precise measurements at the highest energies would allow us to constrain the properties of the pulsar population, such as the magnetic field or initial spin distributions.

We present the first individual and stacking systematic search for $\gamma$-ray emission in the GeV band with the {\it Fermi} Large Area Telescope (LAT), in the directions of 45 super-luminous supernovae (SLSNe). No excess of \gam rays from the SLSN positions was found. We report $\gamma$-ray luminosity upper limits and discuss the implication of these results on the origin of SLSNe in particular on the scenario of central compact object-aided SNe. From the stacking search, we derive an upper limit at 95\% confidence level (CL) to the $\gamma$-ray luminosity (above 600 MeV) $L_{\gamma}<9.1\times10^{41}$\,erg\,s$^{-1}$ for an assumed $E^{-2}$ photon spectrum, for our full SLSN sample. We conclude that the rate of the neutron stars born with millisecond rotation periods {$P\lesssim 2\,$ms and $B\sim10^{12-13}\,$G} must be lower than the rate of the observed SLSNe. The luminosity limits obtained on individual sources are also constraining: in particular, SN2013fc, CSS140222, SN2010kd, and PTF12dam can be born with millisecond periods only if $B\lesssim 10^{13}\,$G.

The Probe Of Extreme Multi-Messenger Astrophysics (POEMMA) mission is being designed to establish charged-particle astronomy with ultra-high energy cosmic rays (UHECRs) and to observe cosmogenic tau neutrinos (CTNs). The study of UHECRs and CTNs from space will yield orders-of-magnitude increase in statistics of observed UHECRs at the highest energies, and the observation of the cosmogenic flux of neutrinos for a range of UHECR models. These observations should solve the long-standing puzzle of the origin of the highest energy particles ever observed, providing a new window onto the most energetic environments and events in the Universe, while studying particle interactions well beyond accelerator energies. The discovery of CTNs will help solve the puzzle of the origin of UHECRs and begin a new field of Astroparticle Physics with the study of neutrino properties at ultra-high energies.

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.

The new generation of powerful instruments is reaching sensitivities and temporal resolutions that will allow multi-messenger astronomy of explosive transient phenomena, with high-energy neutrinos as a central figure. We derive general criteria for the detectability of neutrinos from powerful transient sources, for given instrument sensitivities. In practice, we provide the minimal photon flux necessary for neutrino detection based on two main observables, the bolometric luminosity and the time variability of the emission. This limit can be compared to the observations in specified wavelengths in order to target the most promising sources for follow-ups. Our criteria can also help discriminating false associations of neutrino events with a flaring source. We find that relativistic transient sources such as High- and Low-Luminosity GRBs, Blazar flares, Tidal Disruption Events and magnetar flares could be observed with IceCube, as they have a good chance to occur within a detectable distance. Among non-relativistic transient sources, only luminous supernovae appear as promising candidates. We caution that our criterion should not be directly applied to low-luminosity GRBs and Type Ibc supernovae, as these objects could have hosted a choked GRB, leading to neutrino emission without a relevant counterpart radiation. We treat a set of concrete examples and show that several transients, some of which are being monitored by IceCube, are far from meeting the criterion for detectability (e.g., Crab flares or Swift J1644+57).

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.

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.

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.

The recent detection of the gravitational wave source GW150914 by the LIGO collaboration motivates a speculative source for the origin of ultrahigh energy cosmic rays as a possible byproduct of the immense energies achieved in black hole mergers, provided that the black holes have spin as seems inevitable and there are relic magnetic fields and disk debris remaining from the formation of the black holes or from their accretion history. We argue that given the modest efficiency $< 0.01$ required per event per unit of gravitational wave energy release, merging black holes potentially provide an environment for accelerating cosmic rays to ultrahigh energies.

We investigate the evaporation of close-by pulsar companions, such as planets, asteroids, and white dwarfs, by induction heating. Assuming that the outflow energy is dominated by a Poynting flux (or pulsar wave) at the location of the companions, we calculate their evaporation timescales, by applying the Mie theory. Depending on the size of the companion compared to the incident electromagnetic wavelength, the heating regime varies and can lead to a total evaporation of the companion. In particular, we find that inductive heating is mostly inefficient for small pulsar companions, although it is generally considered the dominant process. Small objects like asteroids can survive induction heating for $10^4\,$years at distances as small as $1\,R_\odot$ from the neutron star. For degenerate companions, induction heating cannot lead to evaporation and another source of heating (likely by kinetic energy of the pulsar wind) has to be considered. It was recently proposed that bodies orbiting pulsars are the cause of fast radio bursts; the present results explain how those bodies can survive in the pulsar's highly energetic environment.

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.

Pulsar wind nebulae (PWNe) are outstanding accelerators in Nature, in the sense that they accelerate electrons up to the radiation reaction limit. Motivated by this observation, this paper examines the possibility that young pulsar wind nebulae can accelerate ions to ultra-high energies at the termination shock of the pulsar wind. We consider here powerful PWNe, fed by pulsars born with $\sim\,$millisecond periods. Assuming that such pulsars exist, at least during a few years after the birth of the neutron star, and that they inject ions into the wind, we find that protons could be accelerated up to energies of the order of the Greisen-Zatsepin-Kuzmin cut-off, for a fiducial rotation period $P\,\sim\,1\,$msec and a pulsar magnetic field $B_{\star}\,\sim\,10^{13}\,$G, implying a fiducial wind luminosity $L_{\rm p}\,\sim\,10^{45}\,$erg/s and a spin-down time $t_{\rm sd}\,\sim 3\times 10^7\,$s. The main limiting factor is set by synchrotron losses in the nebula and by the size of the termination shock; ions with $Z\geq 1$ may therefore be accelerated to even higher energies. We derive an associated neutrino flux produced by interactions in the source region. For a proton-dominated composition, our maximum flux lies slightly below the 5-year sensitivity of IceCube-86 and above the 3-year sensitivity of the projected Askaryan Radio Array. It might thus become detectable in the next decade, depending on the exact level of contribution of these millisecond pulsar wind nebulae to the ultra-high energy cosmic ray flux.

High-energy neutrino astronomy will probe the working of the most violent phenomena in the Universe. The Giant Radio Array for Neutrino Detection (GRAND) project consists of an array of $\sim10^5$ radio antennas deployed over $\sim$200000km$^2$ in a mountainous site. It aims at detecting high-energy neutrinos via the measurement of air showers induced by the decay in the atmosphere of $\tau$ leptons produced by the interaction of the cosmic neutrinos under the Earth surface. Our objective with GRAND is to reach a neutrino sensitivity of $3\times10^{-11}E^{-2}$GeV$^{-1}$cm$^{-2}$s$^{-1}$sr$^{-1}$ above $3 \times10^{16}$eV. This sensitivity ensures the detection of cosmogenic neutrinos in the most pessimistic source models, and about 100 events per year are expected for the standard models. GRAND would also probe the neutrino signals produced at the potential sources of UHECRs. We show how our preliminary design should enable us to reach our sensitivity goals, and present the experimental characteristics. We assess the possibility to adapt GRAND to other astrophysical radio measurements. We discuss in this token the technological options for the detector and the steps to be taken to achieve the GRAND project.

We examine whether charged particles injected by self-annihilating Dark Matter into regions undergoing Diffuse Shock Acceleration (DSA) can be accelerated to high energies. We consider three astrophysical sites where shock acceleration is supposed to occur, namely the Galactic Centre, galaxy clusters and Active Galactic Nuclei (AGN). For the Milky Way, we find that the acceleration of cosmic rays injected by dark matter could lead to a bump in the cosmic ray spectrum provided that the product of the efficiency of the acceleration mechanism and the concentration of DM particles is high enough. Among the various acceleration sources that we consider (namely supernova remnants (SNRs), Fermi bubbles and AGN jets), we find that the Fermi bubbles are a potentially more efficient accelerator than SNRs. However both could in principle accelerate electrons and protons injected by dark matter to very high energies. At the extragalactic level, the acceleration of dark matter annihilation products could be responsible for enhanced radio emission from colliding clusters and prediction of an increase of the anti-deuteron flux generated near AGNs.

Young, fast-rotating neutron stars are promising candidate sources for the production of ultrahigh energy cosmic rays (UHECRs). The interest in this model has recently been boosted by the latest chemical composition measurements of cosmic rays, that seem to show the presence of a heavy nuclear component at the highest energies. Neutrons stars, with their metal-rich surfaces, are potentially interesting sources of such nuclei, but some open issues remain: 1) is it possible to extract these nuclei from the star's surface? 2) Do the nuclei survive the severe conditions present in the magnetosphere of the neutron star? 3) What happens to the surviving nuclei once they enter the wind that is launched outside the light cylinder? In this paper we address these issues in a quantitative way, proving that for the most reasonable range of neutron star surface temperatures ($T<10^7\,$K), a large fraction of heavy nuclei survive photo-disintegration losses. These processes, together with curvature losses and acceleration in the star's electric potential, lead to injection of nuclei with a chemical composition that is mixed, even if only iron is extracted from the surface. We show that under certain conditions the chemical composition injected into the wind region is compatible with that required in previous work based on purely phenomenological arguments (typically $\sim 50\%$ protons, $\sim 30\%$ CNO and $\sim 20\%$ Fe), and provides a reasonable explanation of the mass abundance inferred from ultra high energy data.

Recent results from the AMS-02 data have confirmed that the cosmic ray positron fraction increases with energy between 10 and 200GeV. This quantity should not exceed 50%, and it is hence expected that it will either converge towards 50% or fall. We study the possibility that future data may show the positron fraction dropping down abruptly to the level expected with only secondary production, and forecast the implications of such a feature in term of possible injection mechanisms that include both Dark Matter and pulsars.

We explore the potential of a future, ultra-high energy cosmic ray (UHECR) experiment, that is able to overcome the limitation of low statistics, to detect anisotropy in the arrival directions of UHECRs. We concentrate on the lower energy range of future instruments (E > 50 EeV), where, if the UHECR source number density is not too low, the sources should be numerous enough to imprint a clustering pattern in the sky, and thus possibly in the UHECR arrival directions. Under these limits, the anisotropy signal should be dominated by the clustering of astrophysical sources per-se in the large-scale structures, and not the clustering of events around individual sources. We study the potential for a statistical discrimination between different astrophysical models which we parametrise by the number density of UHECR sources, the possible bias of the UHECR accelerators with respect to the galaxy distribution, and the unknown fraction of UHECRs that have been deflected by large angles. We demonstrate that an order-of-magnitude increase in statistics would allow to discriminate between a variety of astrophysical models, provided that a sub-sample of light elements can be extracted, and that it represents at least ~70% of the overall flux, sensitive to the UHECR source number density. Even without knowledge of the composition, an anisotropy at the 99.7% level should be detectable when the number of detected events exceeds 2000 beyond 50 EeV, as long as the composition is proton dominated, and the number density of UHECR sources is relatively high. If the UHECR sources are strongly biased relative to the galaxy distribution, as are for example galaxy clusters, an anisotropy at the 99.7% should be detectable once the number of detected events exceeds 1000, if the fraction of protons at the highest energies is ~60% or higher.

Recent results from the AMS-02 data have confirmed that the cosmic ray positron fraction increases with energy between 10 and 200GeV. This quantity should not exceed 50%, and it is hence expected that it will either converge towards 50% or fall. We study the possibility that future data may show the positron fraction dropping down abruptly to the level expected with only secondary production, and forecast the implications of such a feature in term of possible injection mechanisms that include both dark matter and pulsars. {Were a sharp steepening to be found, rather surprisingly, we conclude that pulsar models would do at least as well as dark matter scenarios in terms of accounting for any spectral cut-off.

High frequency peaked high redshift blazars, are extreme in the sense that their spectrum is particularly hard and peaks at TeV energies. Standard leptonic scenarios require peculiar source parameters and/or a special setup in order to account for these observations. Electromagnetic cascades seeded by ultra-high energy cosmic rays (UHECR) in the intergalactic medium have also been invoked, assuming a very low intergalactic magnetic field (IGMF). Here we study the synchrotron emission of UHECR secondaries produced in blazars located in magnetised environments, and show that it can provide an alternative explanation to these challenged channels, for sources embedded in structured regions with magnetic field strengths of the order of $10^{-7}$ G. To demonstrate this, we focus on three extreme blazars: 1ES 0229+200, RGB J0710+591, and 1ES 1218+304. We model the expected gamma-ray signal from these sources through a combination of numerical Monte Carlo simulations and solving the kinetic equations of the particles in our simulations, and explore the UHECR source and intergalactic medium parameter space to test the robustness of the emission. We show that the generated synchrotron pair halo/echo flux at the peak energy is not sensitive to variations in the overall IGMF strength. This signal is unavoidable in contrast to the inverse Compton pair halo/echo intensity, which is appealing in view of the large uncertainties on the IGMF in voids of large scale structure. It is also shown that the variability of blazar gamma-ray emission can be accommodated by the synchrotron emission of secondary products of UHE neutral beams if these are emitted by UHECR accelerators inside magnetised regions.

We discuss the possibility that some recently measured anisotropic cosmic ray components in the TeV-PeV energy range may be an indication of the ejection of a peculiar type of matter. We present a model where a neutron star internal transition with nuclear deconfinement of the quark content takes place. This catastrophic event may cause a mass ejection process seeding the insterstelar medium with droplets of quark matter, so called nuclearites. Neutralization of these droplets in molecular clouds may drive the anisotropy since quasi-rectilinear trajectories are allowed. Complementary information from current experimental settings on earth or magnetic spectrometers on the ISS may shed light on this exotic form of matter.

We explore the effect of pulsars, in particular those born with millisecond periods, on their surrounding supernova ejectas. While they spin down, fast-spinning pulsars release their tremendous rotational energy in the form of a relativistic magnetized wind that can affect the dynamics and luminosity of the supernova. We estimate the thermal and non thermal radiations expected from these specific objects, concentrating at times a few years after the onset of the explosion. We find that the bolometric light curves present a high luminosity plateau (that can reach $10^{43-44}\,$erg/s) over a few years. An equally bright TeV gamma-ray emission, and a milder X-ray peak (of order $10^{40-42}\,$erg/s) could also appear a few months to a few years after the explosion, as the pulsar wind nebula emerges, depending on the injection parameters. The observations of these signatures by following the emission of a large number of supernovae could have important implications for the understanding of core-collapse supernovae and reveal the nature of the remnant compact object.

Several experiments (e.g., Milagro and IceCube) have reported the presence in the sky of regions with significant excess in the arrival direction distributions of Galactic cosmic rays in the TeV to PeV energy range. Here we study the possibility that these hotspots are a manifestation of the peculiar nature of these cosmic rays, and of the presence of molecular clouds near the sources. We propose that stable quark matter lumps or so-called "strangelets" can be emitted in the course of the transition of a neutron star to a more compact astrophysical object. A fraction of these massive particles would lose their charge by spallation or electron capture in molecular clouds located in the immediate neighborhood of their source, and propagate rectilinearly without decaying further, hence inducing anisotropies of the order of the cloud size. With reasonable astrophysical assumptions regarding the neutron star transition rate, strangelet injection and neutralization rates, we can reproduce successfully the observed hotspot characteristics and their distribution in the sky.

The acceleration of ultrahigh energy nuclei in fast spinning newborn pulsars can explain the observed spectrum of ultrahigh energy cosmic rays and the trend towards heavier nuclei for energies above $10^{19}\,$eV as reported by the Auger Observatory. Pulsar acceleration implies a hard injection spectrum ($\sim E^{-1}$) due to pulsar spin down and a maximum energy $E_{\rm max} \sim Z \, 10^{19}$ eV due to the limit on the spin rate of neutron stars. We have previously shown that the escape through the young supernova remnant softens the spectrum, decreases slightly the maximum energy, and generates secondary nuclei. Here we show that the distribution of pulsar birth periods and the effect of propagation in the interstellar and intergalactic media modifies the combined spectrum of all pulsars. By assuming a normal distribution of pulsar birth periods centered at 300 ms, we show that the contribution of extragalactic pulsar births to the ultrahigh energy cosmic ray spectrum naturally gives rise to a contribution to very high energy cosmic rays (VHECRs, between $10^{16}$ and $10^{18}$ eV) by Galactic pulsar births. The required injected composition to fit the observed spectrum depends on the absolute energy scale, which is uncertain, differing between Auger Observatory and Telescope Array. The contribution of Galactic pulsar births can also bridge the gap between predictions for cosmic ray acceleration in supernova remnants and the observed spectrum just below the ankle, depending on the composition of the cosmic rays that escape the supernova remnant and the diffusion behavior of VHECRs in the Galaxy.

The Pierre Auger Observatory is exploring the potential of the radio detection technique to study extensive air showers induced by ultra-high energy cosmic rays. The Auger Engineering Radio Array (AERA) addresses both technological and scientific aspects of the radio technique. A first phase of AERA has been operating since September 2010 with detector stations observing radio signals at frequencies between 30 and 80 MHz. In this paper we present comparative studies to identify and optimize the antenna design for the final configuration of AERA consisting of 160 individual radio detector stations. The transient nature of the air shower signal requires a detailed description of the antenna sensor. As the ultra-wideband reception of pulses is not widely discussed in antenna literature, we review the relevant antenna characteristics and enhance theoretical considerations towards the impulse response of antennas including polarization effects and multiple signal reflections. On the basis of the vector effective length we study the transient response characteristics of three candidate antennas in the time domain. Observing the variation of the continuous galactic background intensity we rank the antennas with respect to the noise level added to the galactic signal.

The Pierre Auger Observatory is a facility built to detect air showers produced by cosmic rays above 10^17 eV. During clear nights with a low illuminated moon fraction, the UV fluorescence light produced by air showers is recorded by optical telescopes at the Observatory. To correct the observations for variations in atmospheric conditions, atmospheric monitoring is performed at regular intervals ranging from several minutes (for cloud identification) to several hours (for aerosol conditions) to several days (for vertical profiles of temperature, pressure, and humidity). In 2009, the monitoring program was upgraded to allow for additional targeted measurements of atmospheric conditions shortly after the detection of air showers of special interest, e.g., showers produced by very high-energy cosmic rays or showers with atypical longitudinal profiles. The former events are of particular importance for the determination of the energy scale of the Observatory, and the latter are characteristic of unusual air shower physics or exotic primary particle types. The purpose of targeted (or "rapid") monitoring is to improve the resolution of the atmospheric measurements for such events. In this paper, we report on the implementation of the rapid monitoring program and its current status. The rapid monitoring data have been analyzed and applied to the reconstruction of air showers of high interest, and indicate that the air fluorescence measurements affected by clouds and aerosols are effectively corrected using measurements from the regular atmospheric monitoring program. We find that the rapid monitoring program has potential for supporting dedicated physics analyses beyond the standard event reconstruction.

The Surface Detector of the Pierre Auger Observatory is sensitive to neutrinos of all flavours above 0.1 EeV. These interact through charged and neutral currents in the atmosphere giving rise to extensive air showers. When interacting deeply in the atmosphere at nearly horizontal incidence, neutrinos can be distinguished from regular hadronic cosmic rays by the broad time structure of their shower signals in the water-Cherenkov detectors. In this paper we present for the first time an analysis based on down-going neutrinos. We describe the search procedure, the possible sources of background, the method to compute the exposure and the associated systematic uncertainties. No candidate neutrinos have been found in data collected from 1 January 2004 to 31 May 2010. Assuming an E^-2 differential energy spectrum the limit on the single flavour neutrino is (E^2 * dN/dE) < 1.74x10^-7 GeV cm^-2 s^-1 sr^-1 at 90% C.L. in the energy range 1x10^17 eV < E < 1x10^20 eV.

Newly-born pulsars offer favorable sites for the injection of heavy nuclei, and for their further acceleration to ultrahigh energies. Once accelerated in the pulsar wind, nuclei have to escape from the surrounding supernova envelope. We examine this escape analytically and numerically, and discuss the pulsar source scenario in light of the latest ultrahigh energy cosmic ray (UHECR) data. Our calculations show that, at early times, when protons can be accelerated to energies E>10^20 eV, the young supernova shell tends to prevent their escape. In contrast, because of their higher charge, iron-peaked nuclei are still accelerated to the highest observed energies at later times, when the envelope has become thin enough to allow their escape. Ultrahigh energy iron nuclei escape newly-born pulsars with millisecond periods and dipole magnetic fields of ~10^(12-13) G, embedded in core-collapse supernovae. Due to the production of secondary nucleons, the envelope crossing leads to a transition of composition from light to heavy elements at a few EeV, as observed by the Auger Observatory. The escape also results in a softer spectral slope than that initially injected via unipolar induction, which allows for a good fit to the observed UHECR spectrum. We conclude that the acceleration of iron-peaked elements in a reasonably small fraction (< 0.01%) of extragalactic rotation-powered young pulsars would reproduce satisfactorily the current UHECR data. Possible signatures of this scenario are also discussed.

Atmospheric conditions at the site of a cosmic ray observatory must be known for reconstructing observed extensive air showers. The Global Data Assimilation System (GDAS) is a global atmospheric model predicated on meteorological measurements and numerical weather predictions. GDAS provides altitude-dependent profiles of the main state variables of the atmosphere like temperature, pressure, and humidity. The original data and their application to the air shower reconstruction of the Pierre Auger Observatory are described. By comparisons with radiosonde and weather station measurements obtained on-site in Malarg\"ue and averaged monthly models, the utility of the GDAS data is shown.

We present a comprehensive study of the influence of the geomagnetic field on the energy estimation of extensive air showers with a zenith angle smaller than $60^\circ$, detected at the Pierre Auger Observatory. The geomagnetic field induces an azimuthal modulation of the estimated energy of cosmic rays up to the ~2% level at large zenith angles. We present a method to account for this modulation of the reconstructed energy. We analyse the effect of the modulation on large scale anisotropy searches in the arrival direction distributions of cosmic rays. At a given energy, the geomagnetic effect is shown to induce a pseudo-dipolar pattern at the percent level in the declination distribution that needs to be accounted for.

In this paper we introduce the concept of Lateral Trigger Probability (LTP) function, i.e., the probability for an extensive air shower (EAS) to trigger an individual detector of a ground based array as a function of distance to the shower axis, taking into account energy, mass and direction of the primary cosmic ray. We apply this concept to the surface array of the Pierre Auger Observatory consisting of a 1.5 km spaced grid of about 1600 water Cherenkov stations. Using Monte Carlo simulations of ultra-high energy showers the LTP functions are derived for energies in the range between 10^{17} and 10^{19} eV and zenith angles up to 65 degs. A parametrization combining a step function with an exponential is found to reproduce them very well in the considered range of energies and zenith angles. The LTP functions can also be obtained from data using events simultaneously observed by the fluorescence and the surface detector of the Pierre Auger Observatory (hybrid events). We validate the Monte-Carlo results showing how LTP functions from data are in good agreement with simulations.

We present the results of an analysis of data recorded at the Pierre Auger Observatory in which we search for groups of directionally-aligned events (or `multiplets') which exhibit a correlation between arrival direction and the inverse of the energy. These signatures are expected from sets of events coming from the same source after having been deflected by intervening coherent magnetic fields. The observation of several events from the same source would open the possibility to accurately reconstruct the position of the source and also measure the integral of the component of the magnetic field orthogonal to the trajectory of the cosmic rays. We describe the largest multiplets found and compute the probability that they appeared by chance from an isotropic distribution. We find no statistically significant evidence for the presence of multiplets arising from magnetic deflections in the present data.

Astrophysical observations of ultra-high-energy cosmic rays with the Pierre Auger Observatory

Studies of the cosmic ray energy spectrum at the highest energies with the Pierre Auger Observatory

Technical reports on operations and monitoring of the Pierre Auger Observatory

Studies of the composition of the highest energy cosmic rays with the Pierre Auger Observatory, including examination of hadronic physics effects on the structure of extensive air showers.

Ongoing and planned enhancements of the Pierre Auger Observatory

The Pierre Auger Collaboration has reported evidence for anisotropy in the distribution of arrival directions of the cosmic rays with energies $E>E_{th}=5.5\times 10^{19}$ eV. These show a correlation with the distribution of nearby extragalactic objects, including an apparent excess around the direction of Centaurus A. If the particles responsible for these excesses at $E>E_{th}$ are heavy nuclei with charge $Z$, the proton component of the sources should lead to excesses in the same regions at energies $E/Z$. We here report the lack of anisotropies in these directions at energies above $E_{th}/Z$ (for illustrative values of $Z=6,\ 13,\ 26$). If the anisotropies above $E_{th}$ are due to nuclei with charge $Z$, and under reasonable assumptions about the acceleration process, these observations imply stringent constraints on the allowed proton fraction at the lower energies.

Newly born magnetars are good candidate sources of ultrahigh energy cosmic rays. These objects can in principle easily accelerate particles to the highest energies required to satisfy the ultrahigh energy cosmic ray scenario (E~10^{20-21} eV), thanks to their important rotational and magnetic energy reservoirs. Their acceleration mechanism, based on unipolar induction, predicts however a hard particle injection that does not fit the observed ultrahigh energy cosmic ray spectrum. Here we show that an adequate distribution of initial voltages among magnetar winds can be found to soften the spectrum. We discuss the effect of these distributions for the stochastic gravitational wave background signature produced by magnetars. The magnetar population characteristics needed to fit the ultrahigh energy cosmic ray spectrum could lead in most optimistic cases to gravitational wave background signals enhanced of up to four orders of magnitudes in the range of frequency 1-100 Hz, compared to the standard predictions. These signals could reach the sensitivities of future detectors such as DECIGO or BBO.