Journal cover Journal topic
Geoscientific Model Development An interactive open-access journal of the European Geosciences Union
Journal topic

Journal metrics

Journal metrics

  • IF value: 4.252 IF 4.252
  • IF 5-year value: 4.890 IF 5-year 4.890
  • CiteScore value: 4.49 CiteScore 4.49
  • SNIP value: 1.539 SNIP 1.539
  • SJR value: 2.404 SJR 2.404
  • IPP value: 4.28 IPP 4.28
  • h5-index value: 40 h5-index 40
  • Scimago H index value: 51 Scimago H index 51
Volume 11, issue 11 | Copyright
Geosci. Model Dev., 11, 4515-4535, 2018
© Author(s) 2018. This work is distributed under
the Creative Commons Attribution 4.0 License.

Model evaluation paper 13 Nov 2018

Model evaluation paper | 13 Nov 2018

Evaluation of Monte Carlo tools for high-energy atmospheric physics II: relativistic runaway electron avalanches

David Sarria1, Casper Rutjes2, Gabriel Diniz2,3, Alejandro Luque4, Kevin M. A. Ihaddadene5, Joseph R. Dwyer5, Nikolai Østgaard1, Alexander B. Skeltved1, Ivan S. Ferreira3, and Ute Ebert2,6 David Sarria et al.
  • 1Birkeland Centre for Space Science, Department of Physics and Technology, University of Bergen, Bergen, Norway
  • 2Centrum Wiskunde & Informatica (CWI), Amsterdam, the Netherlands
  • 3Instituto de Física, Universidade de Brasília, Brasília, Brazil
  • 4Instituto de Astrofísica de Andalucía (IAA-CSIC), P.O. Box 3004, Granada, Spain
  • 5University of New Hampshire Main Campus, Department of Physics, Durham, NH, USA
  • 6Eindhoven University of Technology, Eindhoven, the Netherlands

Abstract. The emerging field of high-energy atmospheric physics studies how high-energy particles are produced in thunderstorms, in the form of terrestrial γ-ray flashes and γ-ray glows (also referred to as thunderstorm ground enhancements). Understanding these phenomena requires appropriate models of the interaction of electrons, positrons and photons with air molecules and electric fields. We investigated the results of three codes used in the community – Geant4, GRanada Relativistic Runaway simulator (GRRR) and Runaway Electron Avalanche Model (REAM) – to simulate relativistic runaway electron avalanches (RREAs). This work continues the study of Rutjes et al. (2016), now also including the effects of uniform electric fields, up to the classical breakdown field, which is about 3.0MVm−1 at standard temperature and pressure.

We first present our theoretical description of the RREA process, which is based on and incremented over previous published works. This analysis confirmed that the avalanche is mainly driven by electric fields and the ionisation and scattering processes determining the minimum energy of electrons that can run away, which was found to be above  ≈ 10keV for any fields up to the classical breakdown field.

To investigate this point further, we then evaluated the probability to produce a RREA as a function of the initial electron energy and of the magnitude of the electric field. We found that the stepping methodology in the particle simulation has to be set up very carefully in Geant4. For example, a too-large step size can lead to an avalanche probability reduced by a factor of 10 or to a 40% overestimation of the average electron energy. When properly set up, both Geant4 models show an overall good agreement (within  ≈ 10%) with REAM and GRRR. Furthermore, the probability that particles below 10keV accelerate and participate in the high-energy radiation is found to be negligible for electric fields below the classical breakdown value. The added value of accurately tracking low-energy particles (< 10keV) is minor and mainly visible for fields above 2MVm−1.

In a second simulation set-up, we compared the physical characteristics of the avalanches produced by the four models: avalanche (time and length) scales, convergence time to a self-similar state and energy spectra of photons and electrons. The two Geant4 models and REAM showed good agreement on all parameters we tested. GRRR was also found to be consistent with the other codes, except for the electron energy spectra. That is probably because GRRR does not include straggling for the radiative and ionisation energy losses; hence, implementing these two processes is of primary importance to produce accurate RREA spectra. Including precise modelling of the interactions of particles below 10keV (e.g. by taking into account molecular binding energy of secondary electrons for impact ionisation) also produced only small differences in the recorded spectra.

Publications Copernicus
Short summary
We evaluate three models (Geant4, REAM, GRRR) used in the field of high-energy atmospheric physics that are able to simulate relativistic runaway electron avalanches. Several models have been used by the community, but there was, up until now, no study evaluating their consistency in this context. We conclude that there are no major differences to report, and we discuss minor ones. We also provide advice on how to properly set up the general purpose code (Geant4) in this context.
We evaluate three models (Geant4, REAM, GRRR) used in the field of high-energy atmospheric...