Volume 654, October 2021
|Number of page(s)||11|
|Section||Numerical methods and codes|
|Published online||15 October 2021|
Modelling magnetised medium particle transport in the guiding centre limit with GEANT4⋆
Universidade Presbiteriana Mackenzie, Centro de Rádio Astronomia e Astrofísica Mackenzie, São Paulo, Brazil
2 Instituto de Astronomía y Física del Espacio, CONICET-UBA, Buenos Aires, Argentina
3 University of Glasgow, School of Physics and Astronomy, Glasgow, UK
Accepted: 2 August 2021
Monte Carlo codes are a standard tool for studying energetic particle propagation, secondary production, and radiation in astrophysical settings. In magnetised plasmas such as those found in solar active regions, the enormous disparity between particle gyroradii and system scales proves to be a major computational obstacle. To address this problem we have written a new module in Geant4 using the guiding centre (GC) approach in which the particle motion is averaged over a gyrofrequency. We describe the formulation and implementation of this method in particular dealing with the uncertainty in gyrophase so that particle velocities are well-defined for input to the modules handling reactions. As far as feasible, we compare the propagation and slowing down of primary protons, secondary particle production, and run times in the GC limit with the Newton–Lorentz approach, finding very good agreement between the two methods and orders of magnitude improvement in run times in the GC case. Finally, we present an illustrative solar physics application involving two interacting dipoles, which is only achievable using the GC approach.
Key words: astroparticle physics / methods: numerical / Sun: flares / Sun: magnetic fields / Sun: particle emission
The Geant4 codes used in this work, along with running instructions, are available from a public GitHub repository (https://github.com/guigue/Geant4_Codes_for_Solar_Simulations/tree/master).
© ESO 2021
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.