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: 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: Measurements of air showers made using the hybrid technique developed with the fluorescence and surface detectors of the Pierre Auger Observatory allow a sensitive search for point sources of EeV photons anywhere in the exposed sky. A multivariate analysis reduces the background of hadronic cosmic rays. The search is sensitive to a declination band from -85 degrees to +20 degrees, in an energy range from 10(17.3) eV to 10(18.5) eV. No photon point source has been detected. An upper limit on the photon flux has been derived for every direction. The mean value of the energy flux limit that results from this, assuming a photon spectral index of -2, is 0.06 eV cm(-2) s(-1), and no celestial direction exceeds 0.25 eV cm(-2) s(-1). These upper limits constrain scenarios in which EeV cosmic ray protons are emitted by non-transient sources in the Galaxy.

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: Data collected by the Pierre Auger Observatory through 31 August 2007 showed evidence for anisotropy in the arrival directions of cosmic rays above the Greisen-Zatsepin-Kuz'min energy threshold, 6 x 10(19) eV. The anisotropy was measured by the fraction of arrival directions that are less than 3.1 degrees from the position of an active galactic nucleus within 75 Mpc (using the Veron-Cetty and Veron 12th catalog). An updated measurement of this fraction is reported here using the arrival directions of cosmic rays recorded above the same energy threshold through 31 December 2009. The number of arrival directions has increased from 27 to 69, allowing a more precise measurement. The correlating fraction is (38(-6)(+7))%, compared with 21% expected for isotropic cosmic rays. This is down from the early estimate of (69-(+11)(13))%. The enlarged set of arrival directions is examined also in relation to other populations of nearby extragalactic objects: galaxies in the 2 Microns All Sky Survey and active galactic nuclei detected in hard X-rays by the Swift Burst Alert Telescope. A celestial region around the position of the radiogalaxy Cen A has the largest excess of arrival directions relative to isotropic expectations. The 2-point autocorrelation function is shown for the enlarged set of arrival directions and compared to the isotropic expectation.