Fluence calculation methods in Monte Carlo dosimetry simulations

Modeling of scavenging systems in water radiolysis with Geant4-DNA

PURPOSE

This paper presents the capabilities of the Geant4-DNA Monte Carlo toolkit to simulate water radiolysis with scavengers using the step-by-step (SBS) or the independent reaction times (IRT) methods. It features two examples of application areas: (1) computing the escape yield of H₂O₂ following a ⁶⁰Co γ-irradiation and (2) computing the oxygen depletion in water irradiated with 1 MeV electrons.

METHODS

To ease the implementation of the chemical stage in Geant4-DNA, we developed a user interface that helps define the chemical reactions and set the concentration of scavengers. The first application area example required two computational steps to perform water radiolysis using NO₂^- and NO₃^- as scavengers and a ⁶⁰Co irradiation. The oxygen depletion computation technique for the second application area example consisted of simulating track segments of 1 MeV electrons and determining the radio-induced loss and gain of oxygen molecules.

RESULTS

The production of H₂O₂ under variable scavenging levels is consistent with the literature; the mean relative difference between the SBS and IRT methods is 7.2 % ± 0.5 %. For the oxygen depletion 1 µs post-irradiation, the mean relative difference between both methods is equal to 9.8 % ± 0.3 %. The results in the microsecond scale depend on the initial partial pressure of oxygen in water. In addition, the computed oxygen depletions agree well with the literature.

CONCLUSIONS

The Geant4-DNA toolkit makes it possible to simulate water radiolysis in the presence of scavengers. This feature offers perspectives in radiobiology, with the possibility of simulating cell-relevant scavenging mechanisms.

Enhancement of the EGSnrc code egs_chamber for fast fluence calculations of charged particles

PURPOSE

Simulation of absorbed dose deposition in a detector is one of the key tasks of Monte Carlo (MC) dosimetry methodology. Recent publications (Hartmann and Zink, 2018; Hartmann and Zink, 2019; Hartmann et al., 2021) have shown that knowledge of the charged particle fluence differential in energy contributing to absorbed dose is useful to provide enhanced insight on how response depends on detector properties. While some EGSnrc MC codes provide output of charged particle spectra, they are often restricted in setup options or limited in calculation efficiency. For detector simulations, a promising approach is to upgrade the EGSnrc code egs_chamber which so far does not offer charged particle calculations.

METHODS

Since the user code cavity offers charged particle fluence calculation, the underlying algorithm was embedded in egs_chamber. The modified code was tested against two EGSnrc applications and DOSXYZnrc which was modified accordingly by one of the authors. Furthermore, the gain in efficiency achieved by photon cross section enhancement was determined quantitatively.

RESULTS

Electron and positron fluence spectra and restricted cema calculated by egs_chamber agreed well with the compared applications thus demonstrating the feasibility of the new code. Additionally, variance reduction techniques are now applicable also for fluence calculations. Depending on the simulation setup, considerable gains in efficiency were obtained by photon cross section enhancement.

CONCLUSION

The enhanced egs_chamber code represents a valuable tool to investigate the response of detectors with respect to absorbed dose and fluence distribution and the perturbation caused by the detector in a reasonable computation time. By using intermediate phase space scoring, egs_chamber offers parallel calculation of charged particle fluence spectra for different detector configurations in one single run.

Cema-based formalism for the determination of absorbed dose for high-energy photon beams

PURPOSE

Determination of absorbed dose is well established in many dosimetry protocols and considered to be highly reliable using ionization chambers under reference conditions. If dosimetry is performed under other conditions or using other detectors, however, open questions still remain. Such questions frequently refer to appropriate correction factors. A converted energy per mass (cema)-based approach to formulate such correction factors offers a good understanding of the specific response of a detector for dosimetry under various measuring conditions and thus an estimate of pros and cons of its application.

METHODS

Determination of absorbed dose requires the knowledge of the beam quality correction factor k_Q,Qo , where Q denotes the quality of a user beam and Qo is the quality of the radiation used for calibration. In modern Monte Carlo (MC)-based methods, k_Q,Qo is directly derived from the MC-calculated dose conversion factor, which is the ratio between the absorbed dose at a point of interest in water and the mean absorbed dose in the sensitive volume of an ion chamber. In this work, absorbed dose is approximated by the fundamental quantity cema. This approximation allows the dose conversion factor to be substituted by the cema conversion factor. Subsequently, this factor is decomposed into a product of cema ratios. They are identified as the stopping power ratio water to the material in the sensitive detector volume, and as the correction factor for the fluence perturbation of the secondary charged particles in the detector cavity caused by the presence of the detector. This correction factor is further decomposed with respect to the perturbation caused by the detector cavity and that caused by external detector properties. The cema-based formalism was subsequently tested by MC calculations of the spectral fluence of the secondary charged particles (electrons and positrons) under various conditions.

RESULTS

MC calculations demonstrate that considerable fluence perturbation may occur particularly under non-reference conditions. Cema-based correction factors to be applied in a 6-MV beam were obtained for a number of ionization chambers and for three solid-state detectors. Feasibility was shown at field sizes of 4 × 4 and 0.5 cm × 0.5 cm. Values of the cema ratios resulting from the decomposition of the dose conversion factor can be well correlated with detector response. Under the small field conditions, the internal fluence correction factor of ionization chambers is considerably dependent on volume averaging and thus on the shape and size of the cavity volume.

CONCLUSIONS

The cema approach is particularly useful at non-reference conditions including when solid-state detectors are used. Perturbation correction factors can be expressed and evaluated by cema ratios in a comprehensive manner. The cema approach can serve to understand the specific response of a detector for dosimetry to be dependent on (a) radiation quality, (b) detector properties, and (c) electron fluence changes caused by the detector. This understanding may also help to decide which detector is best suited for a specific measurement situation.

Towards the standardization of the absorbed dose report mode in high energy photon beams

The benefits of using an algorithm that reports absorbed dose-to-medium have been jeopardized by the clinical experience and the experimental protocols that have mainly relied on absorbed dose-to-water. The aim of the present work was to investigate the physical aspects that govern the dosimetry in heterogeneous media using Monte Carlo method and to introduce a formalism for the experimental validation of absorbed dose-to-medium reporting algorithms. Particle fluence spectra computed within the sensitive volume of two simulated detectors (T31016 Pinpoint 3D ionization chamber and EBT3 radiochromic film) placed in different media (water, RW3, lung and bone) were compared to those in the undisturbed media for 6 MV photon beams. A heterogeneity correction factor that takes into account the difference between the detector perturbation in medium and under reference conditions as well as the stopping-power ratios was then derived for all media using cema calculations. Furthermore, the different conversion approaches and Eclipse treatment planning system algorithms were compared against the Monte Carlo absorbed dose reports. The detectors electron fluence perturbation in RW3 and lung media were close to that in water (≤1.5%). However, the perturbation was greater in bone (∼4%) and impacted the spectral shape. It was emphasized that detectors readings should be corrected by the heterogeneity correction factor that ranged from 0.932 in bone to 0.985 in lung. Significant discrepancies were observed between all the absorbed dose reports and conversions, especially in bone (exceeding 10%) and to a lesser extent in RW3. Given the ongoing advances in dose calculation algorithms, it is essential to standardize the absorbed dose report mode with absorbed dose-to-medium as a favoured choice. It was concluded that a retrospective conversion should be avoided and switching from absorbed dose-to-water to absorbed dose-to-medium reporting algorithm should be carried out by a direct comparison of both algorithms.