Abstract: A flux of neutrons from an astrophysical source in the Galaxy can be detected in the Pierre Auger Observatory as an excess of cosmic-ray air showers arriving from the direction of the source. To avoid the statistical penalty for making many trials, classes of objects are tested in combinations as nine “target sets,” in addition to the search for a neutron flux from the Galactic center or from the Galactic plane. Within a target set, each candidate source is weighted in proportion to its electromagnetic flux, its exposure to the Auger Observatory, and its flux attenuation factor due to neutron decay. These searches do not find evidence for a neutron flux from any class of candidate sources. Tabulated results give the combined p-value for each class, with and without the weights, and also the flux upper limit for the most significant candidate source within each class. These limits on fluxes of neutrons significantly constrain models of EeV proton emission from non-transient discrete sources in the Galaxy.

Abstract: We present the results of the most complete scan of the parameter space for cosmic ray (CR) injection and propagation. We perform a Bayesian search of the main GALPROP parameters, using the MultiNest nested sampling algorithm, augmented by the BAMBI neural network machine-learning package. This is the first study to separate out low-mass isotopes (p, (p) over bar and He) from the usual light elements (Be, B, C, N, and O). We find that the propagation parameters that best-fit p, (p) over bar, and He data are significantly different from those that fit light elements, including the B/C and Be-10/Be-9 secondary-to-primary ratios normally used to calibrate propagation parameters. This suggests that each set of species is probing a very different interstellar medium, and that the standard approach of calibrating propagation parameters using B/C can lead to incorrect results. We present posterior distributions and best-fit parameters for propagation of both sets of nuclei, as well as for the injection abundances of elements from H to Si. The input GALDEF files with these new parameters will be included in an upcoming public GALPROP update.

Abstract: We analyze 7.3 yr of ANTARES high-energy neutrino and Fermi Large Area Telescope (LAT) gamma-ray data in search of cosmic neutrino + gamma-ray (nu + gamma) transient sources or source populations. Our analysis has the potential to detect either individual nu + gamma transient sources (durations delta t less than or similar to 1000 s), if they exhibit sufficient gamma-ray or neutrino multiplicity, or a statistical excess of nu + gamma transients of individually lower multiplicities. Individual high gamma-ray multiplicity events could be produced, for example, by a single ANTARES neutrino in coincidence with a LAT-detected gamma-ray burst. Treating ANTARES track and cascade event types separately, we establish detection thresholds by Monte Carlo scrambling of the neutrino data, and determine our analysis sensitivity by signal injection against these scrambled data sets. We find our analysis is sensitive to nu + gamma transient populations responsible for >5% of the observed gamma-coincident neutrinos in the track data at 90% confidence. Applying our analysis to the unscrambled data reveals no individual nu + gamma events of high significance; two ANTARES track + Fermi gamma-ray events are identified that exceed a once per decade false alarm rate threshold (p = 17%). No evidence for subthreshold nu + gamma source populations is found among the track (p = 39%) or cascade (p = 60%) events. Exploring a possible correlation of high-energy neutrino directions with Fermi gamma-ray sky brightness identified in previous work yields no added support for this correlation. While TXS.0506+056, a blazar and variable (nontransient) Fermi gamma-ray source, has recently been identified as the first source of high-energy neutrinos, the challenges in reconciling observations of the Fermi gamma-ray sky, the IceCube high-energy cosmic neutrinos, and ultrahigh-energy cosmic rays using only blazars suggest a significant contribution by other source populations. Searches for transient sources of high-energy neutrinos thus remain interesting, with the potential for either neutrino clustering or multimessenger coincidence searches to lead to discovery of the first nu + gamma transients.

Abstract: The exploration of the universe has recently entered a new era thanks to the multi-messenger paradigm, characterized by a continuous increase in the quantity and quality of experimental data that is obtained by the detection of the various cosmic messengers (photons, neutrinos, cosmic rays and gravitational waves) from numerous origins. They give us information about their sources in the universe and the properties of the intergalactic medium. Moreover, multi-messenger astronomy opens up the possibility to search for phenomenological signatures of quantum gravity. On the one hand, the most energetic events allow us to test our physical theories at energy regimes which are not directly accessible in accelerators; on the other hand, tiny effects in the propagation of very high energy particles could be amplified by cosmological distances. After decades of merely theoretical investigations, the possibility of obtaining phenomenological indications of Planck-scale effects is a revolutionary step in the quest for a quantum theory of gravity, but it requires cooperation between different communities of physicists (both theoretical and experimental). This review, prepared within the COST Action CA18108 “Quantum gravity phenomenology in the multi-messenger approach”, is aimed at promoting this cooperation by giving a state-of-the art account of the interdisciplinary expertise that is needed in the effective search of quantum gravity footprints in the production, propagation and detection of cosmic messengers.

Abstract: The Tibet AS gamma and LHAASO collaborations recently reported the observation of a gamma-ray diffuse emission with energy up to the PeV level from the Galactic plane.Aims. We discuss the relevance of non-uniform cosmic-ray transport scenarios and the implications of these results for cosmic-ray physics.Methods. We used the DRAGON and HERMES codes to build high-resolution maps and spectral distributions of that emission for several representative models under the condition that they reproduce a wide set of local cosmic-ray data up to 100 PeV.Results. We show that the energy spectra measured by Tibet AS gamma, LHAASO, ARGO-YBJ, and Fermi-LAT in several regions of interest in the sky can all be reasonably described in terms of the emission arising by the Galactic cosmic-ray “sea”. We also show that all our models are compatible with IceTop gamma-ray upper limits.Conclusions. We compare the predictions of conventional and space-dependent transport models with those data sets. Although the Fermi-LAT, ARGO-YBJ, and LHAASO preliminary data slightly favor this scenario, due to the still large experimental errors, the poorly known source spectral shape at the highest energies, the potential role of spatial fluctuations in the leptonic component, and a possible larger-than-expected contamination due to unresolved sources, a solid confirmation requires further investigations. We discuss which measurements will be most relevant in order to resolve the remaining degeneracy.