A computationally efficient moment-preserving Monte Carlo electron transport method with implementation in Geant4

Dixon, David A.; Prinja, Anil K.; Franke, Brian C.

doi:10.1016/j.nimb.2015.07.009

Title: A computationally efficient moment-preserving Monte Carlo electron transport method with implementation in Geant4

Abstract

This paper presents the theoretical development and numerical demonstration of a moment-preserving Monte Carlo electron transport method. Foremost, a full implementation of the moment-preserving (MP) method within the Geant4 particle simulation toolkit is demonstrated. Beyond implementation details, it is shown that the MP method is a viable alternative to the condensed history (CH) method for inclusion in current and future generation transport codes through demonstration of the key features of the method including: systematically controllable accuracy, computational efficiency, mathematical robustness, and versatility. A wide variety of results common to electron transport are presented illustrating the key features of the MP method. In particular, it is possible to achieve accuracy that is statistically indistinguishable from analog Monte Carlo, while remaining up to three orders of magnitude more efficient than analog Monte Carlo simulations. Lastly, it is shown that the MP method can be generalized to any applicable analog scattering DCS model by extending previous work on the MP method beyond analytical DCSs to the partial-wave (PW) elastic tabulated DCS data.

Authors:

Dixon, David A. ^[1]; Prinja, Anil K. ^[2]; Franke, Brian C. ^[3]

Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
Univ. of New Mexico, Albuquerque, NM (United States)
Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)

Publication Date:: Fri Jul 24 00:00:00 EDT 2015

Research Org.:: Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)

Sponsoring Org.:: USDOE National Nuclear Security Administration (NNSA)

OSTI Identifier:: 1512911

Alternate Identifier(s):: OSTI ID: 1422674

Report Number(s):: SAND-015-7698J
Journal ID: ISSN 0168-583X; 665169

Grant/Contract Number:: AC04-94AL85000

Resource Type:: Accepted Manuscript

Journal Name:: Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms

Additional Journal Information:: Journal Volume: 359; Journal Issue: C; Journal ID: ISSN 0168-583X

Publisher:: Elsevier

Country of Publication:: United States

Language:: English

Subject:: 46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND TECHNOLOGY; Electron transport; Moment preservation; Partial-wave differential cross sections; Monte Carlo simulation

Citation Formats


                    Dixon, David A., Prinja, Anil K., and Franke, Brian C. A computationally efficient moment-preserving Monte Carlo electron transport method with implementation in Geant4.  United States: N. p., 2015. 
Web.  doi:10.1016/j.nimb.2015.07.009.

Copy to clipboard


                    Dixon, David A., Prinja, Anil K., & Franke, Brian C. A computationally efficient moment-preserving Monte Carlo electron transport method with implementation in Geant4.  United States.  https://doi.org/10.1016/j.nimb.2015.07.009

Copy to clipboard


                    Dixon, David A., Prinja, Anil K., and Franke, Brian C. Fri .  
"A computationally efficient moment-preserving Monte Carlo electron transport method with implementation in Geant4".  United States.  https://doi.org/10.1016/j.nimb.2015.07.009.  https://www.osti.gov/servlets/purl/1512911.

Copy to clipboard


                    
@article{osti_1512911,

  title        = {A computationally efficient moment-preserving Monte Carlo electron transport method with implementation in Geant4},

  author       = {Dixon, David A. and Prinja, Anil K. and Franke, Brian C.},

  abstractNote = {This paper presents the theoretical development and numerical demonstration of a moment-preserving Monte Carlo electron transport method. Foremost, a full implementation of the moment-preserving (MP) method within the Geant4 particle simulation toolkit is demonstrated. Beyond implementation details, it is shown that the MP method is a viable alternative to the condensed history (CH) method for inclusion in current and future generation transport codes through demonstration of the key features of the method including: systematically controllable accuracy, computational efficiency, mathematical robustness, and versatility. A wide variety of results common to electron transport are presented illustrating the key features of the MP method. In particular, it is possible to achieve accuracy that is statistically indistinguishable from analog Monte Carlo, while remaining up to three orders of magnitude more efficient than analog Monte Carlo simulations. Lastly, it is shown that the MP method can be generalized to any applicable analog scattering DCS model by extending previous work on the MP method beyond analytical DCSs to the partial-wave (PW) elastic tabulated DCS data.},

  doi          = {10.1016/j.nimb.2015.07.009},

  journal      = {Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms},

  number       = C,

  volume       = 359,

  place        = {United States},

  year         = {Fri Jul 24 00:00:00 EDT 2015},

  month        = {Fri Jul 24 00:00:00 EDT 2015}

}

Copy to clipboard

Journal Article:

Free Publicly Available Full Text

Accepted Manuscript (Publisher)

Accepted Manuscript (DOE)

Publisher's Version of Record

https://doi.org/10.1016/j.nimb.2015.07.009

Other availability

Search WorldCat to find libraries that may hold this journal

Citation Metrics:

Cited by: 2 works

Citation information provided by
Web of Science

Figures / Tables:

Figure 1: Partial-wave DCSs for elastic scattering of electrons with energies between 1-keV and 20000-keV by (a) aluminum and (b) gold nuclei.

All figures and tables (21 total)

Save / Share:

Export Metadata

Save to My Library

Figures / Tables found in this record:

Figures/Tables have been extracted from DOE-funded journal article accepted manuscripts.

Similar Records in DOE PAGES and OSTI.GOV collections:

Application of the GBFP Method to Electron Partial-Wave Expansion Elastic Scattering Differential Cross Sections Within the Geant4 Toolkit - Paper 69

Conference Dixon, David A. ; Prinja, Anil K. ; Franke, Brian C.

Analog Monte Carlo simulation of electrons in complex material geometries is generally impractical because the long range Coulomb interactions result in highly peaked differential cross sections (DCS) for scattering and energy-loss and extremely small collision mean free paths (mfp). The condensed history (CH) Monte Carlo method attempts to alleviate this limitation by advancing electrons in fixed steps and sampling angular deflections and energy losses from precomputed distributions obtained from approximate solutions to the underlying transport equation that are valid only for thin layers. Though CH is efficient and is widely employed in production codes (ITS, MCNP, EGS, Geant), step sizemore »« less
A moment-preserving nonanalog method for charged particle transport.

Conference Prinja, Anil K ; Franke, Brian Claude ; Fichtl, Erin ; ...

Extremely short collision mean free paths and near-singular elastic and inelastic differential cross sections (DCS) make analog Monte Carlo simulation an impractical tool for charged particle transport. The widely used alternative, the condensed history method, while efficient, also suffers from several limitations arising from the use of precomputed smooth distributions for sampling. There is much interest in developing computationally efficient algorithms that implement the correct transport mechanics. Here we present a nonanalog transport-based method that incorporates the correct transport mechanics and is computationally efficient for implementation in single event Monte Carlo codes. Our method systematically preserves important physics and ismore »« less
The GEANT4 toolkit for microdosimetry calculations: Application to microbeam radiation therapy (MRT)

Journal Article Spiga, J ; Siegbahn, E A ; Braeuer-Krisch, E ; ... - Medical Physics

Theoretical dose distributions for microbeam radiation therapy (MRT) are computed in this paper using the GEANT4 Monte Carlo (MC) simulation toolkit. MRT is an innovative experimental radiotherapy technique carried out using an array of parallel microbeams of synchrotron-wiggler-generated x rays. Although the biological mechanisms underlying the effects of microbeams are still largely unknown, the effectiveness of MRT can be traced back to the natural ability of normal tissues to rapidly repair small damages to the vasculature, and on the lack of a similar healing process in tumoral tissues. Contrary to conventional therapy, in which each beam is at least severalmore »« less
https://doi.org/10.1118/1.2794170
Region-oriented CT image representation for reducing computing time of Monte Carlo simulations

Journal Article Sarrut, David ; Guigues, Laurent ; CREATIS Laboratory - Medical Physics

Purpose. We propose a new method for efficient particle transportation in voxelized geometry for Monte Carlo simulations. We describe its use for calculating dose distribution in CT images for radiation therapy. Material and methods. The proposed approach, based on an implicit volume representation named segmented volume, coupled with an adapted segmentation procedure and a distance map, allows us to minimize the number of boundary crossings, which slows down simulation. The method was implemented with the GEANT4 toolkit and compared to four other methods: One box per voxel, parameterized volumes, octree-based volumes, and nested parameterized volumes. For each representation, we comparedmore »« less
https://doi.org/10.1118/1.2884854
Using Geant4 in the BaBar Simulation

Technical Report Wright, Dennis H

BABAR was the first large experiment to incorporate Geant4 into its detector simulation. Since July 2001, 1.5 x 10{sup 9} BABAR events have been produced using this simulation. In a typical e{sup +}e{sup -} {yields} {Upsilon} (4s) {yields} B{sup 0}{bar B}{sup 0} event, between 30 and 60 tracks are produced in the generator and propagated through the detector, using decay, electromagnetic and hadronic processes provided by the Geant4 toolkit. The material model of the detector is very detailed and a special particle transportation module was developed so that minute features (on the few micron scale) would be sampled in themore »« less
https://doi.org/10.2172/813141

Full Text Available

Similar Records