|
Johannesson, G., Ruiz de Austri, R., Vincent, A. C., Moskalenko, I. V., Orlando, E., Porter, T. A., et al. (2016). Bayesian analysis of cosmic-ray propagation: evidence against homogeneous diffusion. Astrophys. J., 824(1), 16–19pp.
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.
|
|
|
Trotta, R., Johannesson, G., Moskalenko, I. V., Porter, T. A., Ruiz de Austri, R., & Strong, A. W. (2011). Constraints on Cosmic-Ray Propagation Models from a Global Bayesian Analysis. Astrophys. J., 729(2), 106–16pp.
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.
|
|
|
de los Rios, M., Petac, M., Zaldivar, B., Bonaventura, N. R., Calore, F., & Iocco, F. (2023). Determining the dark matter distribution in simulated galaxies with deep learning. Mon. Not. Roy. Astron. Soc., 525(4), 6015–6035.
Abstract: We present a novel method of inferring the dark matter (DM) content and spatial distribution within galaxies, using convolutional neural networks (CNNs) trained within state-of-the-art hydrodynamical simulations (Illustris-TNG100). Within the controlled environment of the simulation, the framework we have developed is capable of inferring the DM mass distribution within galaxies of mass similar to 10(11)-10(13)M(circle dot) from the gravitationally baryon-dominated internal regions to the DM-rich, baryon-depleted outskirts of the galaxies, with a mean absolute error always below approximate to 0.25 when using photometrical and spectroscopic information. With respect to traditional methods, the one presented here also possesses the advantages of not relying on a pre-assigned shape for the DM distribution, to be applicable to galaxies not necessarily in isolation, and to perform very well even in the absence of spectroscopic observations.
|
|
|
Schiavone, T., Montani, G., & Bombacigno, F. (2023). f(R) gravity in the Jordan frame as a paradigm for the Hubble tension. Mon. Not. Roy. Astron. Soc., 522(1), L72–L77.
Abstract: We analyse the f(R) gravity in the so-called Jordan frame, as implemented to the isotropic Universe dynamics. The goal of the present study is to show that according to recent data analyses of the supernovae Ia Pantheon sample, it is possible to account for an effective redshift dependence of the Hubble constant. This is achieved via the dynamics of a non-minimally coupled scalar field, as it emerges in the f(R) gravity. We face the question both from an analytical and purely numerical point of view, following the same technical paradigm. We arrive to establish that the expected decay of the Hubble constant with the redshift z is ensured by a form of the scalar field potential, which remains essentially constant for z less than or similar to 0.3, independently if this request is made a priori, as in the analytical approach, or obtained a posteriori, when the numerical procedure is addressed. Thus, we demonstrate that an f(R) dark energy model is able to account for an apparent variation of the Hubble constant due to the rescaling of the Einstein constant by the f(R) scalar mode.
|
|
|
De La Torre Luque, P., Gaggero, D., Grasso, D., Fornieri, O., Egberts, K., Steppa, C., et al. (2023). Galactic diffuse gamma rays meet the PeV frontier. Astron. Astrophys., 672, A58–11pp.
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.
|
|