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: 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: Results are presented of a search for cosmic sources of high-energy neutrinos with the ANTARES neutrino telescope. The data were collected during 2007 and 2008 using detector configurations containing between 5 and 12 detection lines. The integrated live time of the analyzed data is 304 days. Muon tracks are reconstructed using a likelihood-based algorithm. Studies of the detector timing indicate a median angular resolution of 0.5 +/- 0.1 deg. The neutrino flux sensitivity is 7.5 x 10(-8)(E(v)/GeV)(-2) GeV(-1) s(-1) cm(-2) for the part of the sky that is always visible (delta < -48 deg), which is better than limits obtained by previous experiments. No cosmic neutrino sources have been observed.

Abstract: A multimessenger analysis optimized for a correlation of arrival directions of ultra-high energy cosmic rays (UHECRs) and neutrinos is presented and applied to 2190 neutrino candidate events detected in 2007-2008 by the ANTARES telescope and 69 UHECRs observed by the Pierre Auger Observatory between 2004 January 1 and 2009 December 31. No significant correlation is observed. Assuming an equal neutrino flux (E-2 energy spectrum) from all UHECR directions, a 90% CL upper limit on the neutrino flux of 5.0 x 10(-8) GeV cm(-2) s(-1) per source is derived.

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.