Search

Several sources of total cross sections for e+e− pair production by photon interaction with matter, used by major Monte Carlo codes for particle transport, are quantitatively evaluated with respect to experimental measurements collected from the literature. They include data libraries deriving from theoretical calculations and analytical formulae representing empirical fits to tabulations. Although most of these sources of cross sections are based on one theoretical reference, the statistical data analysis documented in this paper highlights differences in their compatibility with experiment. The cross sections collected in the 1997 version of the Evaluated Photon Data Library (EPDL) exhibit the lowest incompatibility with experiment; the paper details the validation results for all the examined cross section sources. The scarcity of experimental data does not allow the validation process to discern pair production in the field of the nucleus and in the field of atomic electrons; nevertheless, this paper documents relevant clues regarding the contribution of the latter to reproduce experimental measurements. [ABSTRACT FROM AUTHOR]

Theoretical and semi-empirical methods to calculate electron impact ionization cross sections for atomic shells are subject to validation tests with respect to a wide collection of experimental measurements to identify the state of the art for Monte Carlo particle transport. The validation process applies rigorous statistical analysis methods. Cross sections based on the EEDL Evaluated Electron Data Library, widely used by Monte Carlo codes, and on calculations by Bote and Salvat, used in the Penelope code, are generally equivalent in compatibility with experiment. Results are also reported for various formulations of the Binary-Encounter-Bethe and Deutsch-Märk models. [ABSTRACT FROM AUTHOR]

Evolutions of Geant4 code have affected the simulation of electron backscattering with respect to previously published results. Their effects are quantified by analyzing the compatibility of the simulated electron backscattering fraction with a large collection of experimental data for a wide set of physics configuration options available in Geant4. Special emphasis is placed on two electron scattering implementations first released in Geant4 version 10.2: the Goudsmit-Saunderson multiple scattering model and a single Coulomb scattering model based on Mott cross section calculation. The new Goudsmit-Saunderson multiple scattering model appears to perform equally or less accurately than the model implemented in previous Geant4 versions, depending on the electron energy. The new Coulomb scattering model was flawed from a physics point of view, but computationally fast in Geant4 version 10.2; the physics correction released in Geant4 version 10.2p01 severely degrades its computational performance. Problems observed in electron backscattering simulation in previous publications have been addressed by evolutions in the Geant4 geometry domain. [ABSTRACT FROM AUTHOR]

We propose a new type of computational phantom, the “4D voxel phantom,” for realistic modeling of continuous respiratory motion in Monte Carlo dose calculation. In this phantom, continuous respiratory motion is realized by linear interpolation of the deformation vector fields (DVFs) between the neighboring original phases in the 4D CT data of a patient and by subsequent application of the DVFs to the phase images or to the reference image to produce multiple inter-phase images between the neighboring original phase images. A 4D voxel phantom is a combination of high-temporal-resolution voxel phantoms and on-the-fly dose registration to the reference phase image. In the course of particle transport simulation, the dose or deposited energy is directly registered to the reference phase image on-the-fly (i.e., after each event) using a DVF for dose registration. In the present study, we investigated two methods — DRP (DIR [deformable image registration] with respect to Reference Phase image) and DNP (DIR with respect to Neighboring original Phase image) — for production of multiple inter-phase images or high-temporal-resolution voxel phantoms. Utilizing these two methods, two 4D voxel phantoms each with 100 phases were produced from the original 10-phase images of the 4D CT data of a real patient in order to compare the two methods and to test the feasibility of the 4D voxel phantom methodology in general. We found that it is possible to produce a 4D voxel phantom very rapidly (i.e., <40 min on a 4-core personal computer for a 100-phase phantom) in a fully automated process. The dose calculation results showed that the constructed 100-phase 4D voxel phantoms provide cumulative-dose distributions very similar to those of the conventional 10-phase approach for stationary proton-beam irradiation. The passing rates of the dose distributions of the 4D voxel phantoms were higher than 99.9% according to the 3% and 3 mm gamma criteria, which results validate the 4D voxel phantom methodology. The point-and dose-tracking analysis data showed that the DRP method, which uses the minimal number of DIR operations but uses inverse DVFs, provides significantly better results than those of the DNP method, which uses only DIR to generate the DVFs for inter-phase image generation and dose registration. The present study also showed that the computation time does not significantly increase when the number of phases in the 4D voxel phantom is increased for more realistic representation of continuous respiratory motion; the only significant increase is in the memory occupancy, which grows almost linearly with the number of phases. [ABSTRACT FROM AUTHOR]

Several total and partial photoionization cross section calculations, based on both theoretical and empirical approaches, are quantitatively evaluated with statistical analyses using a large collection of experimental data retrieved from the literature to identify the state of the art for modeling the photoelectric effect in Monte Carlo particle transport. Some of the examined cross section models are available in general purpose Monte Carlo systems, while others have been implemented and subjected to validation tests for the first time to estimate whether they could improve the accuracy of particle transport codes. The validation process identifies Scofield’s 1973 non-relativistic calculations, tabulated in the Evaluated Photon Data Library (EPDL), as the one best reproducing experimental measurements of total cross sections. Specialized total cross section models, some of which derive from more recent calculations, do not provide significant improvements. Scofield’s non-relativistic calculations are not surpassed regarding the compatibility with experiment of K and L shell photoionization cross sections either, although in a few test cases Ebel’s parameterization produces more accurate results close to absorption edges. Modifications to Biggs and Lighthill’s parameterization implemented in Geant4 significantly reduce the accuracy of total cross sections at low energies with respect to its original formulation. The scarcity of suitable experimental data hinders a similar extensive analysis for the simulation of the photoelectron angular distribution, which is limited to a qualitative appraisal. [ABSTRACT FROM AUTHOR]

Backscattering is a sensitive probe of the accuracy of electron scattering algorithms implemented in Monte Carlo codes. The capability of the Geant4 toolkit to describe realistically the fraction of electrons backscattered from a target volume is extensively and quantitatively evaluated in comparison with experimental data retrieved from the literature. The validation test covers the energy range between approximately 100 eV and 20 MeV, and concerns a wide set of target elements. Multiple and single electron scattering models implemented in Geant4, as well as preassembled selections of physics models distributed within Geant4, are analyzed with statistical methods. The evaluations concern Geant4 versions from 9.1 to 10.1. Significant evolutions are observed over the range of Geant4 versions, not always in the direction of better compatibility with experiment. Goodness-of-fit tests complemented by categorical analysis tests identify a configuration based on Geant4 Urban multiple scattering model in Geant4 version 9.1 and a configuration based on single Coulomb scattering in Geant4 10.0 as the physics options best reproducing experimental data above a few tens of keV. At lower energies only single scattering demonstrates some capability to reproduce data down to a few keV. Recommended preassembled physics configurations appear incapable of describing electron backscattering compatible with experiment. With the support of statistical methods, a correlation is established between the validation of Geant4-based simulation of backscattering and of energy deposition. [ABSTRACT FROM PUBLISHER]