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.

Abstract: Research in many areas of modern physics such as, e. g., indirect searches for dark matter and particle acceleration in supernova remnant shocks rely heavily on studies of cosmic rays (CRs) and associated diffuse emissions (radio, microwave, X-rays, gamma-rays). While very detailed numerical models of CR propagation exist, a quantitative statistical analysis of such models has been so far hampered by the large computational effort that those models require. Although statistical analyses have been carried out before using semi-analytical models (where the computation is much faster), the evaluation of the results obtained from such models is difficult, as they necessarily suffer from many simplifying assumptions. The main objective of this paper is to present a working method for a full Bayesian parameter estimation for a numerical CR propagation model. For this study, we use the GALPROP code, the most advanced of its kind, which uses astrophysical information, and nuclear and particle data as inputs to self-consistently predict CRs, gamma-rays, synchrotron, and other observables. We demonstrate that a full Bayesian analysis is possible using nested sampling and Markov Chain Monte Carlo methods (implemented in the SuperBayeS code) despite the heavy computational demands of a numerical propagation code. The best-fit values of parameters found in this analysis are in agreement with previous, significantly simpler, studies also based on GALPROP.

Abstract: The surface detector array of the Pierre Auger Observatory can detect neutrinos with energy E-nu between 10(17) eV and 10(20) eV from point-like sources across the sky south of +55 degrees and north of -65 degrees declinations. A search has been performed for highly inclined extensive air showers produced by the interaction of neutrinos of all flavors in the atmosphere (downward-going neutrinos), and by the decay of tau leptons originating from tau neutrino interactions in Earth's crust (Earth-skimming neutrinos). No candidate neutrinos have been found in data up to 2010 May 31. This corresponds to an equivalent exposure of similar to 3.5 years of a full surface detector array for the Earth-skimming channel and similar to 2 years for the downward-going channel. An improved upper limit on the diffuse flux of tau neutrinos has been derived. Upper limits on the neutrino flux from point-like sources have been derived as a function of the source declination. Assuming a differential neutrino flux k(PS) . E-nu(-2). from a point-like source, 90% confidence level upper limits for k(PS) at the level of approximate to 5x10(-7) and 2.5x10(-6) GeV cm(-2) s(-1) have been obtained over a broad range of declinations from the searches for Earth-skimming and downward-going neutrinos, respectively.

Abstract: In this paper, a time-integrated search for point sources of cosmic neutrinos is presented using the data collected from 2007 to 2010 by the ANTARES neutrino telescope. No statistically significant signal has been found and upper limits on the neutrino flux have been obtained. Assuming an E-nu(-2). spectrum, these flux limits are at 1-10x10(-8) GeV cm(-2) s(-1) for declinations ranging from -90 degrees to 40 degrees. Limits for specific models of RX J1713.7-3946 and Vela X, which include information on the source morphology and spectrum, are also given.

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.