High-energy neutrino flares are interesting prospective counterparts to photon flares, as their detection would guarantee the presence of accelerated hadrons within a source, provide precious information about cosmic-ray acceleration and interactions, and thus impact the subsequent modeling of non-thermal emissions in explosive transients. In these sources, photomeson production can be efficient, producing a large amount of secondary particles, such as charged pions and muons, that decay and produce high-energy neutrinos. Before their decay, secondary particles can experience energy losses and acceleration, which can impact high-energy neutrino spectra and thus affect their detectability. In this work, we focus on the impact of secondary acceleration. We consider a one zone model, characterized mainly by a variability timescale $t_{\rm var}$, a luminosity $L_{\rm bol}$, a bulk Lorentz factor $\Gamma$. The mean magnetic field $B$ is deduced from these parameters. The photon field is modeled by a broken power-law. This generic model allows to evaluate systematically the maximum energy of high-energy neutrinos in the parameter space of explosive transients, and shows that it could be strongly affected by secondary acceleration for a large number of source categories. In order to determine the impact of secondary acceleration on the high-energy neutrino spectrum and in particular on its peak energy and flux, we complement these estimates by several case studies. We show that secondary acceleration can increase the maximum neutrino flux, and produce a secondary peak at the maximum energy in the case of efficient acceleration. Secondary acceleration could therefore enhance the detectability of very-high-energy neutrinos, that will be the target of next generation neutrino detectors such as KM3NeT, IceCube-Gen2, POEMMA or 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.

We investigate the capability of the Probe Of Extreme Multi-Messenger Astrophysics (POEMMA) in performing Target-of-Opportunity (ToO) neutrino observations. POEMMA will detect tau neutrinos via Cherenkov radiation from their upward-moving extensive air showers. POEMMA will be able to quickly slew ($90^{\circ}$ in 500 s) to the direction of an astrophysical source, which in combination with its orbital speed will provide it with unparalleled capability to follow up transient alerts. We calculate POEMMA's transient sensitivity for two observational modes for its two satellites (stereo and dual) and investigate variations in neutrino sensitivity across the sky arising from POEMMA's orbit. We explore separate scenarios for long ($\sim 10^6$ s) and short ($\sim 10^3$ s) bursts, accounting for intrusion from the Sun and the Moon in long-duration scenarios. For long bursts, POEMMA will improve the average neutrino sensitivity above 300 PeV by up to a factor of 7 with respect to existing experiments (e.g., IceCube, ANTARES, and Pierre Auger), reaching the level of model predictions for neutrino fluences at these energies from several types of long-duration astrophysical transients (e.g., blazar flares, binary neutron star mergers, and tidal disruption events). For short bursts in the optimal case, POEMMA will improve the sensitivity over existing experiments by at least an order of magnitude above 100 PeV. POEMMA's orbit and rapid slewing will provide access to the full celestial sky, including regions not accessible to ground-based experiments. Finally, we discuss the prospects for detecting neutrinos from candidate astrophysical neutrino sources in the nearby universe. Our results demonstrate that with its improved neutrino sensitivity at ultra-high energies and unique full-sky coverage, POEMMA will be an essential component in an expanding multi-messenger network.

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.

The Probe Of Extreme Multi-Messenger Astrophysics (POEMMA) is designed to observe cosmic neutrinos (CNs) above 20 PeV and ultra-high energy cosmic rays (UHECRs) above 20 EeV over the full sky. The POEMMA mission calls for two identical satellites flying in loose formation, each comprised of a 4-meter wide field-of-view (45 degrees) Schmidt photometer. The hybrid focal surface includes a fast (1 ${\mu}$s) ultraviolet camera for fluorescence observations and an ultrafast (10 ns) optical camera for Cherenkov observations. POEMMA will provide new multi-messenger windows onto the most energetic events in the universe, enabling the study of new astrophysics and particle physics at these otherwise inaccessible energies.

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.

The detection of very-high-energy (VHE) neutrinos ($E_\nu \gtrsim 10$ PeV) is a challenge that future generations of experiments are being designed and constructed to address. One promising method relies on using the Earth as a neutrino target for indirect detection of skimming tau neutrinos interacting within the Earth and producing tau leptons that are able to escape and decay in the atmosphere. The tau decay produces upward-moving Extensive Air Showers (EASs). A space-based or suborbital instrument observing the ground near the Earth limb can search for the beamed Cherenkov signal produced by the up-going EAS resulting from the tau-lepton decay. In this paper, we derive the geometrical constraints for such an observation in general and for the specific case of the Probe Of Extreme Multi-Messenger Astrophysics (POEMMA) mission currently under study, focusing on the Cherenkov signal detection. We show that, using reasonable orbital parameters, POEMMA can achieve full-sky coverage to search for potential neutrino sources over the length of its mission. We also show that follow-up of a transient Target-of-Opportunity (ToO), such as a flaring source, can be achieved within an orbit time scale depending on the source location on the celestial sphere and its relative position with respect to the Sun and the Moon.

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.

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