A Markov Chain Monte Carlo for Galactic Cosmic Ray physics: I. Method and results for the Leaky Box Model

Abstract

Propagation of charged cosmic-rays in the Galaxy depends on the transport parameters, whose number can be large depending on the propagation model under scrutiny. A standard approach for determining these parameters is a manual scan, leading to an inefficient and incomplete coverage of the parameter space. We implement a Markov Chain Monte Carlo (MCMC), which is well suited to multi-parameter determination. Its specificities (burn-in length, acceptance, and correlation length) are discussed in the phenomenologically well-understood Leaky-Box Model. From a technical point of view, a trial function based on binary-space partitioning is found to be extremely efficient, allowing a simultaneous determination of up to nine parameters, including transport and source parameters, such as slope and abundances. Our best-fit model includes both a low energy cut-off and reacceleration, whose values are consistent with those found in diffusion models. A Kolmogorov spectrum for the diffusion slope (delta=1/3) is excluded. The marginalised probability-density function for delta and alpha (the slope of the source spectra) are delta~0.55-0.60 and alpha~2.14-2.17, depending on the dataset used and the number of free parameters in the fit. All source-spectrum parameters (slope and abundances) are positively correlated among themselves and with the reacceleration strength, but are negatively correlated with the other propagation parameters. A forthcoming study will extend our analysis to more physical diffusion models.

18 pages, 15 figures, one brief section 2 added on the method (matches the A&A accepted version)

References in corpus (29)

• Cosmological parameters from CMB and other data: a Monte-Carlo approach
• Cosmic-ray propagation and interactions in the Galaxy
• Propagation of cosmic-ray nucleons in the Galaxy
• PAMELA - A Payload for Antimatter Matter Exploration and Light-nuclei Astrophysics
• Antiprotons in cosmic rays from neutralino annihilation
• Cosmic Rays below Z=30 in a diffusion model: new constraints on propagation parameters
• Cosmic-Ray Nuclei, Antiprotons and Gamma-rays in the Galaxy: a New Diffusion Model
• Precise Measurement of Cosmic-Ray Proton and Helium Spectra with the BESS Spectrometer
• Precision Measurement of Cosmic-Ray Antiproton Spectrum
• Fast and reliable MCMC for cosmological parameter estimation
• Galactic Cosmic Rays from Supernova Remnants: II Shock Acceleration of Gas and Dust
• Positrons from dark matter annihilation in the galactic halo: theoretical uncertainties
• Antiprotons from spallations of cosmic rays on interstellar matter
• Galactic Cosmic Rays from Supernova Remnants (I) - a Cosmic Ray Composition controlled by Volatility and Mass-to-Charge Ratio
• Protons in near earth orbit
• Measurements of 0.2 to 20 GeV/n cosmic-ray proton and helium spectra from 1997 through 2002 with the BESS spectrometer
• The Cosmic-Ray antiproton flux between 3 and 49 GeV
• Measurements of Cosmic-ray Low-energy Antiproton and Proton Spectra in a Transient Period of the Solar Field Reversal
• The Modified Weighted Slab Technique: Models and Results
• II: Bayesian Methods for Cosmological Parameter Estimation from Cosmic Microwave Background Measurements
• Successive Measurements of Cosmic-Ray Antiproton Spectrum in a Positive Phase of the Solar Cycle
• Composition of Primary Cosmic-Ray Nuclei at High Energies
• Measurement of cosmic-ray low-energy antiproton spectrum with the first BESS-Polar Antarctic flight
• Flux of light antimatter nuclei near earth, induced by Cosmic Rays in the Galaxy and in the atmosphere
• New results on source and diffusion spectral features of Galactic cosmic rays: I- B/C ratio
• Cosmic ray production in supernova remnants including reacceleration: the secondary to primary ratio
• Launch of the Space experiment PAMELA
• Observations of the Li, Be, and B isotopes and constraints on cosmic-ray propagation
• Uncertainties of Cosmic Ray Spectra and Detectability of Antiproton mSUGRA Contributions With PAMELA

Cited by in corpus (36)

• Gamma rays from molecular clouds illuminated by cosmic rays escaping from interacting supernova remnants
• A Markov Chain Monte Carlo technique to sample transport and source parameters of Galactic cosmic rays: II. Results for the diffusion model combining B/C and radioactive nuclei
• Escape-limited Model of Cosmic-ray Acceleration Revisited
• Propagation of cosmic rays in the AMS-02 era
• Dark matter profiles and annihilation in dwarf spheroidal galaxies: prospectives for present and future gamma-ray observatories - I. The classical dSphs
• Spherical Jeans analysis for dark matter indirect detection in dwarf spheroidal galaxies - Impact of physical parameters and triaxiality
• Systematic uncertainties on the cosmic-ray transport parameters: Is it possible to reconcile B/C data with delta = 1/3 or delta = 1/2?
• Systematic Investigation of Solar Modulation of Galactic Protons for Solar Cycle 23 using a Monte Carlo Approach with Particle Drift Effects and Latitudinal Dependence
• Propagation of Cosmic Rays in Heliosphere: the HelMod Model
• Dark Matter Sees The Light
• Implications of the AMS-02 positron fraction in cosmic rays
• Galactic cosmic-ray propagation in the light of AMS-02: I. Analysis of protons, helium, and antiprotons
• Constraining Galactic cosmic-ray parameters with Z<=2 nuclei
• Markov Chain Monte Carlo analysis to constrain Dark Matter properties with directional detection
• AMS-02 antiprotons' consistency with a secondary astrophysical origin
• USINE: semi-analytical models for Galactic cosmic-ray propagation
• On the e${}^{+}$e${}^{-}$ excesses and the knee of the cosmic ray spectra -- hints of cosmic rays acceleration at young supernova remnants
• CosRayMC: a global fitting method in studying the properties of the new sources of cosmic e${}^{\pm }$ excesses
• p, He, and C to Fe cosmic-ray primary fluxes in diffusion models: Source and transport signatures on fluxes and ratios
• Cosmic-ray acceleration and escape from post-adiabatic Supernova remnants
• Deriving the velocity distribution of Galactic Dark Matter particles from rotation curve data
• New minimal, median, and maximal propagation models for dark matter searches with Galactic cosmic rays
• Studies on cosmic ray nuclei with Voyager, ACE and AMS-02: I. local interstellar spectra and solar modulation
• Where are Milky Way's Hadronic PeVatrons?
• The impact of propagation uncertainties on the potential Dark Matter contribution to the Fermi LAT mid-latitude gamma-ray data
• Geminga SNR: Possible candidate of local cosmic-ray factory
• Synchrotron X-ray diagnostics of cutoff shape of nonthermal electron spectrum at young supernova remnants
• The initial conditions of observed star clusters - I. Method description and validation
• A study on the sharp knee and fine structures of cosmic ray spectra
• On the PeV knee of cosmic rays spectrum and TeV cutoff of electron spectrum
• Markov Chain Monte Carlo technics applied to Parton Distribution Functions determination: proof of concept
• Constraining primordial black holes in dark matter with kinematics of dwarf galaxies
• The HelMod Model in the Works for Inner and Outer Heliosphere: from AMS to Voyager Probes Observations
• Neutrino emission from dark matter annihilation/decay in light of cosmic $e^{\pm }$ and $\bar{p}$ data
• Determination of proton PDF uncertainties with Markov chain Monte Carlo
• The Discrepancy in Modelling of Galactic Cosmic Ray Primaries and Secondaries