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: In the past decade, a new class of bright transient radio sources with millisecond duration has been discovered. The origin of these so-called fast radio bursts (FRBs) is still a mystery, despite the growing observational efforts made by various multiwavelength and multimessenger facilities. To date, many models have been proposed to explain FRBs, but neither the progenitors nor the radiative and the particle acceleration processes at work have been clearly identified. In this paper, we assess whether hadronic processes may occur in the vicinity of the FRB source. If they do, FRBs may contribute to the high-energy cosmic-ray and neutrino fluxes. A search for these hadronic signatures was carried out using the ANTARES neutrino telescope. The analysis consists in looking for high-energy neutrinos, in the TeV-PeV regime, that are spatially and temporally coincident with the detected FRBs. Most of the FRBs discovered in the period 2013-2017 were in the field of view of the ANTARES detector, which is sensitive mostly to events originating from the Southern hemisphere. From this period, 12 FRBs were selected and no coincident neutrino candidate was observed. Upper limits on the per-burst neutrino fluence were derived using a power-law spectrum, dN/DE nu proportional to E-nu(-gamma), for the incoming neutrino flux, assuming spectral indexes gamma = 1.0, 2.0, 2.5. Finally, the neutrino energy was constrained by computing the total energy radiated in neutrinos, assuming different distances for the FRBs. Constraints on the neutrino fluence and on the energy released were derived from the associated null results.

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.