On the forward-backward-in-time approach for Monte Carlo solution of Parker's transport equation
P. Bobik ; M. J. Boschini ; S. Della Torre ; et al. ; - ASI Sponsor
Jun - 2016
DOI: 10.1002/2015JA022237

journal : Journal of Geophysical Research: Space Physics

Volume : 121 ; Issue : 5
type: Article Journal

The cosmic rays propagation inside the heliosphere is well described by a transport equation introduced by Parker in 1965. To solve this equation, several approaches were followed in the past. Recently, a Monte Carlo approach became widely used in force of its advantages with respect to other numerical methods. In this approach the transport equation is associated to a fully equivalent set of stochastic differential equations (SDE). This set is used to describe the stochastic path of quasi-particle from a source, e.g., the interstellar space, to a specific target, e.g., a detector at Earth. We present a comparison of forward-in-time and backward-in-time methods to solve the cosmic rays transport equation in the heliosphere. The Parker equation and the related set of SDE in the several formulations are treated in this paper. For the sake of clarity, this work is focused on the one-dimensional solutions. Results were compared with an alternative numerical solution, namely, Crank-Nicolson method, specifically developed for the case under study. The methods presented are fully consistent each others for energy greater than 400 MeV. The comparison between stochastic integrations and Crank-Nicolson allows us to estimate the systematic uncertainties of Monte Carlo methods. The forward-in-time stochastic integrations method showed a systematic uncertainty <5%, while backward-in-time stochastic integrations method showed a systematic uncertainty <1% in the studied energy range.

keywords : Monte Carlo;cosmic rays propagation;heliosphere

Notes : Research funded by Agenzia Spaziale Italiana ASI-INFN I/002/13/0)