Electron beam quality correction factors for plane-parallel ionization chambers: Monte Carlo calculations using the PENELOPE system

Investigations on the beam quality correction factor for ionization chambers in high-energy brachytherapy dosimetry

Objective. To enhance the investigations on MC calculated beam quality correction factors of thimble ionization chambers from high-energy brachytherapy sources and to develop reliable reference conditions in source and detector setups in water.Approach. The response of five different ionization chambers from PTW-Freiburg and Standard Imaging was investigated for irradiation by a high dose rate Ir-192 Flexisource in water. For a setup in a Beamscan water phantom, Monte Carlo simulations were performed to calculate correction factors for the chamber readings. After exact positioning of source and detector the absorbed dose rate at the TG-43 reference point at one centimeter nominal distance from the source was measured using these factors and compared to the specification of the calibration certificate. The Monte Carlo calculations were performed using the restricted cema formalism to gain further insight into the chamber response. Calculations were performed for the sensitive volume of the chambers, determined by the methods currently used in investigations of dosimetry in magnetic fields.Main results. Measured dose rates and values from the calibration certificate agreed within the combined uncertainty (k= 2) for all chambers except for one case in which the full air cavity was simulated. The chambers showed a distinct directional dependence. With the restricted cema formalism calculations it was possible to examine volume averaging and energy dependence of the perturbation factors contributing to the beam quality correction factor also differential in energy.Significance. This work determined beam quality correction factors to measure the absorbed dose rate from a brachytherapy source in terms of absorbed dose to water for a variety of ionization chambers. For the accurate dosimetry of brachytherapy sources with ionization chambers it is advisable to use correction factors based on the sensitive volume of the chambers and to take account for the directional dependence of chamber response.

AAPM WGTG51 Report 385: Addendum to the AAPM's TG-51 protocol for clinical reference dosimetry of high-energy electron beams

An Addendum to the AAPM's TG-51 protocol for the determination of absorbed dose to water is presented for electron beams with energies between 4 MeV and 22 MeV (1.70 cm ≤ R 50 ≤ 8.70 cm $1.70\nobreakspace {\rm cm} \le R_{\text{50}} \le 8.70\nobreakspace {\rm cm}$ ). This updated formalism allows simplified calibration procedures, including the use of calibrated cylindrical ionization chambers in all electron beams without the use of a gradient correction. Newk Q $k_{Q}$ data are provided for electron beams based on Monte Carlo simulations. Implementation guidance is provided. Components of the uncertainty budget in determining absorbed dose to water at the reference depth are discussed. Specifications for a reference-class chamber in electron beams include chamber stability, settling, ion recombination behavior, and polarity dependence. Progress in electron beam reference dosimetry is reviewed. Although this report introduces some major changes (e.g., gradient corrections are implicitly included in the electron beam quality conversion factors), they serve to simplify the calibration procedure. Results for absorbed dose per linac monitor unit are expected to be up to approximately 2 % higher using this Addendum compared to using the original TG-51 protocol.

Electron interaction with DNA constituents in aqueous phase

Electron driven chemistry of biomolecules in aqueous phase presents the realistic picture to study molecular processes. In this study we have investigated the interactions of electrons with the DNA constituents in their aqueous phase in order to obtain the quantities useful for DNA damage assessment. We have computed the inelastic mean free path (IMFP), mass stopping power (MSP) and absorbed dose (D) for the DNA constituents (Adenine, Cytosine, Guanine, Thymine and Uracil) in the aqueous medium from ionisation threshold to 5000 eV. We have modified complex optical potential formalism to include band gap of the systems to calculate inelastic cross sections which are used to estimate these entities. This is the maiden attempt to report these important quantities for the aqueous DNA constituents. We have compared our results with available data in gas and other phase and have observed explicable accord for IMFP and MSP. Since these are the first results of absorbed dose (D) for these compounds, we have explored present results vis-a-vis dose absorption in water.

Investigation of ionization chamber perturbation factors using proton beam and Fano cavity test for the Monte Carlo simulation code PHITS

The reference dose for clinical proton beam therapy is based on ionization chamber dosimetry. However, data on uncertainties in proton dosimetry are lacking, and multifaceted studies are required. Monte Carlo simulations are useful tools for calculating ionization chamber dosimetry in radiation fields and are sensitive to the transport algorithm parameters when particles are transported in a heterogeneous region. We aimed to evaluate the proton transport algorithm of the Particle and Heavy Ion Transport Code System (PHITS) using the Fano test. The response of the ionization chamber f Q and beam quality correction factors k Q were calculated using the same parameters as those in the Fano test and compared with those of other Monte Carlo codes for verification. The geometry of the Fano test consisted of a cylindrical gas-filled cavity sandwiched between two cylindrical walls. f Q was calculated as the ratio of the absorbed dose in water to the dose in the cavity in the chamber. We compared the f Q calculated using PHITS with that of a previous study, which was calculated using other Monte Carlo codes (Geant4, FULKA, and PENH) under similar conditions. The flight mesh, a parameter for charged particle transport, passed the Fano test within 0.15%. This was shown to be sufficiently accurate compared with that observed in previous studies. The f Q calculated using PHITS were 1.116 ± 0.002 and 1.124 ± 0.003 for NACP-02 and PTW-30013, respectively, and the k Q were 0.981 ± 0.008 and 1.027 ± 0.008, respectively, at 150 MeV. Our results indicate that PHITS can calculate the f Q and k Q with high precision.

A multicenter study of modified electron beam output calibration

PURPOSE

This study aims to assess the accuracy of a modified electron beam calibration based on the IAEA TRS-398 and AAPM-TG-51 in multicenter radiotherapy.

METHODS

This study was performed using the Elekta and Varian Linear Accelerator electron beams with energies of 4-22 MeV under reference conditions using cylindrical (PTW 30013, IBA FC65-G, and IBA FC65-P) and parallel-plate (PTW 34045, PTW 34001, and IBA PPC-40) chambers. The modified calibration used a cylindrical chamber and an updatedk ' Q $k{^{\prime}}_Q$ based on Monte Carlo calculations, whereas TRS-398 and TG-51 used cylindrical and parallel-plate chambers for reference dosimetry. The dose ratio of the modified calibration procedure, TRS-398 and TG-51 were obtained by comparing the dose at the maximum depth of the modified calibration to TRS-398 and TG-51.

RESULTS

The study found that all cylindrical chambers' beam quality conversion factors determined with the modified calibration( k ' Q ) $( {{{k^{\prime}}}_Q} )$ to the TRS-398 and TG-51 vary from 0.994 to 1.003 and 1.000 to 1.010, respectively. The dose ratio of modified/TRS-398cyl and modified/TRS-398parallel-plate, the variation ranges were 0.980-1.014 and 0.981-1.019, while for the counterpart modified/TG-51cyl was found varying between 0.991 and 1.017 and the ratio of modified/TG-51parallel-plate varied in the range of 0.981-1.019.

CONCLUSION

This multi-institutional study analyzed a modified calibration procedure utilizing new data for electron beam calibrations at multiple institutions and evaluated existing calibration protocols. Based on observed variations, the current calibration protocols should be updated with detailed metrics on the stability of linac components.

The influence of different versions of FLUKA and GEANT4 on the calculation of response functions of ionization chambers in clinical proton beams

Objective.To investigate the influence of different versions of the Monte Carlo codesgeant4 andflukaon the calculation of overall response functionsfQof air-filled ionization chambers in clinical proton beams.Approach. fQfactors were calculated for six plane-parallel and four cylindrical ionization chambers withgeant4 andfluka. These factors were compared to already published values that were derived using older versions of these codes.Main results.Differences infQfactors calculated with different versions of the same Monte Carlo code can be up to ∼1%. Especially forgeant4, the updated version leads to a more pronounced dependence offQon proton energy and to smallerfQfactors for high energies.Significance.Different versions of the same Monte Carlo code can lead to differences in the calculation offQfactors of up to ∼1% without changing the simulation setup, transport parameters, ionization chamber geometry modeling, or employed physics lists. These findings support the statement that the dominant contributor to the overall uncertainty of Monte Carlo calculatedfQfactors are type-B uncertainties.

Monte Carlo-calculated beam quality and perturbation correction factors validated against experiments for Farmer and Markus type ionization chambers in therapeutic carbon-ion beams

Objective. In current dosimetry protocols, the estimated uncertainty of the measured absorbed dose to waterDwin carbon-ion beams is approximately 3%. This large uncertainty is mainly contributed by the standard uncertainty of the beam quality correction factorkQ. In this study, thekQvalues in four cylindrical chambers and two plane-parallel chambers were calculated using Monte Carlo (MC) simulations in the plateau region. The chamber-specific perturbation correction factorPof each chamber was also determined through MC simulations.Approach.kQfor each chamber was calculated using MC code Geant4. The simulatedkQratios in subjected chambers and reference chambers were validated through comparisons against our measured values. In the measurements in Heavy-Ion Medical Accelerator in Chiba,kQratios were obtained fromDwvalues of60Co, 290- and 400 MeV u-1carbon-ion beams that were measured with the subjected ionization chamber and the reference chamber. In the simulations,fQ(the product of the water-to-air stopping power ratio andP) was acquired fromDwand the absorbed dose to air calculated in the sensitive volume of each chamber.kQvalues were then calculated from the simulatedfQand the literature-extractedWairand compared with previous publications.Main results. The calculatedkQratios in the subjected chambers to the reference chamber agreed well with the measuredkQratios. ThekQuncertainty was reduced from the current recommendation of approximately 3% to 1.7%. ThePvalues were close to unity in the cylindrical chambers and nearly 1% above unity in the plane-parallel chambers.Significance. ThekQvalues of carbon-ion beams were accurately calculated in MC simulations and thekQratios were validated through ionization chamber measurements. The results indicate a need for updating the current recommendations, which assume a constantPof unity in carbon-ion beams, to recommendations that consider chamber-induced differences.

Monte Carlo calculated ionization chamber correction factors in clinical proton beams - deriving uncertainties from published data

For the update of the IAEA TRS-398 Code of Practice (CoP), global ionization chamber factors (fQ) and beam quality correction factors (kQ) for air-filled ionization chambers in clinical proton beams have been calculated with different Monte Carlo codes. In this study, average Monte Carlo calculated fQ and kQ factors are provided and the uncertainty of these factors is estimated. Average fQ factors in monoenergetic proton beams with energies between 60 MeV and 250 MeV were derived from Monte Carlo calculated fQ factors published in the literature. Altogether, 195 fQ factors for six plane-parallel and three cylindrical ionization chambers calculated with penh, fluka and geant4 were incorporated. Additionally, a weighted standard deviation of fQ factors was calculated, where the same weight was assigned to each Monte Carlo code. From average fQ factors, kQ factors were derived and compared to the values from the IAEA TRS-398 CoP published in 2000 as well as to the values of the upcoming version. Average Monte Carlo calculated fQ factors are constant within 0.6% over the energy range investigated. In general, the different Monte Carlo codes agree within 1% for low energies and show larger differences up to 2% for high energies. As a result, the standard deviation of fQ factors increases with energy and is ∼0.3% for low energies and ∼0.8% for high energies. kQ factors derived from average Monte Carlo calculated fQ factors differ from the values presented in the IAEA TRS-398 CoP by up to 2.4%. The overall estimated uncertainty of Monte Carlo calculated kQ factors is ∼0.5%-1% smaller than the uncertainties estimated in IAEA TRS-398 CoP since the individual ionization chamber characteristics (e.g. fluence perturbations) are considered in detail in Monte Carlo calculations. The agreement between Monte Carlo calculated kQ factors and the values of the upcoming version of IAEA TRS-398 CoP is better with deviations smaller than 1%.

Current best estimates of beam quality correction factors for reference dosimetry of clinical proton beams

Objective. To review the currently available data on beam quality correction factors, k Q , for ionization chambers in clinical proton beams and derive their current best estimates for the updated recommendations of the IAEA TRS-398 Code of Practice. Approach. The reviewed data come from 20 publications from which k Q values can be derived either directly from calorimeter measurements, indirectly from comparison with other chambers or from Monte Carlo calculated overall chamber factors, f Q . For cylindrical ionization chambers, a distinction is made between data obtained in the centre of a spread-out Bragg peak and those obtained in the plateau region of single-energy fields. For the latter, the effect of depth dose gradients has to be considered. To this end an empirical model for previously published displacement correction factors for single-layer scanned beams was established, while for unmodulated scattered beams experimental data were used. From all the data, chamber factors, f Q , and chamber perturbation correction factors, p Q , were then derived and analysed. Main results. The analysis showed that except for the beam quality dependence of the water-to-air mass stopping power ratio and, for cylindrical ionization chambers in unmodulated beams, of the displacement correction factor, there is no remaining beam quality dependence of the chamber perturbation correction factors p Q . Based on this approach, average values of the beam quality independent part of the perturbation factors were derived to calculate k Q values consistent with the data in the literature. Significance. The resulting data from this analysis are current best estimates of k Q values for modulated scattered beams and single-layer scanned beams used in proton therapy. Based on this, a single set of harmonized values is derived to be recommended in the update of IAEA TRS-398.

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.

Monte Carlo calculation of beam quality correction factors in proton beams using FLUKA

Purpose.To provide Monte Carlo calculated beam quality correction factors (k_Q) for monoenergetic proton beams using the Monte Carlo codefluka.Materials and methods.The Monte Carlo codeflukawas used to calculate the dose absorbed in a water-filled reference volume and the air-filled cavities of six plane-parallel and four cylindrical ionization chambers. The chambers were positioned at the entrance region of monoenergetic proton beams with energies between 60 and 250 MeV. Based on these dose values,f_Qas well ask_Qfactors were calculated whilefQ0factors were taken from Andreoet al(2020Phys. Med. Biol.65095011).Results. k_Qfactors calculated in this work were found to agree with experimentally determinedk_Qfactors on the 1%-level, with only two exceptions with deviations of 1.4% and 1.9%. The comparison off_Qfactors calculated usingflukawithf_Qfactors calculated using the Monte Carlo codesgeant4 andpenhshowed a general good agreement for low energies, while differences for higher energies were pronounced. For high energies, in most cases the Monte Carlo codesflukaandgeant4 lead to comparable results while thef_Qfactors calculated withpenhare larger.Conclusion.flukacan be used to calculatek_Qfactors in clinical proton beams. The divergence of Monte Carlo calculatedk_Qfactors for high energies suggests that the role of nuclear interaction models implemented in the different Monte Carlo codes needs to be investigated in more detail.

Stopping-power ratios for electron beams used in total skin electron therapy

PURPOSE

The electron beams for total skin electron therapy (TSET) are often degraded by a scatter plate in addition to extended distances. For electron dosimetry, both the AAPM TG-51 and IAEA TRS-398 recommend the use of two formulas developed by Burns et al [Med. Phys. 23, 489-501 (1996)] to estimate the water-to-air stopping-power ratios (SPRs). Both formulas are based on a fit to SPRs calculated for standard electron beams. This study aims to find: (1) if the formulas are applicable to beams used in TSET and (2) the impact of the ICRU report 90 recommendations on the SPRs for these beams.

METHODS

The EGSnrc Monte Carlo code system is used to generate 6 MeV high dose rate total skin electron (HDTSe) beams used in TSET. The simulated beams are used to calculate dose distributions and SPRs as a function of depth in a water phantom. The fitted SPRs using the empirical formulas are compared with MC-calculated SPRs.

RESULTS

The electron beam quality specifier, the depth in water at which the absorbed dose falls to 50% of its maximum value, R₅₀ , decreases approximately 1 mm for each additional 100-cm extended distance ranging from 2.24 cm at SSD = 100 to 1.72 cm at SSD = 700 cm. For beams passing through a scatter plate, R₅₀ is 1.76 cm (1.14) at SSD = 300 and 1.48 cm (0.85 cm) at SSD = 600 cm with an Acrylic plate thickness of 3 mm (9 mm), respectively. The discrepancy between fitted and MC-calculated SPRs at d_ref as a function of R₅₀ is <0.8%, and in many cases <0.4%. The difference between fitted and MC-calculated SPRs as a function of depth and R₅₀ is within 1% at depths <0.8R₅₀ for beams with R₅₀  ≥ 1.14 cm. The ICRU-90 recommendations decrease SPRs by 0.3%-0.4% compared to the use of data recommended in ICRU-37.

CONCLUSION

The formulas used by the major protocols are accurate enough for clinical beams used in TSET and the error caused using the formulas is <1% to estimate SPRs as a function of depth and R₅₀ for depths <0.8R₅₀ for beams used in TSET with R₅₀  ≥ 1.14 cm. The impact of the ICRU-90 recommendations shows a decrease of SPRs by a fraction of a percent for beams used in TSET.

PENELOPE simulations and experiment for 6 MV clinac iX accelerator for standard and small static fields

The goal of this work was to produce accurate data for use as a 'gold standard' and a valid tool for measurements in reference dosimetry for standard/small static field sizes from 0.5 × 0.5 to 10 × 10 cm². It is based on the accuracy of the phase space files (PSFs) as a key quantity. Because the IAEA general public database provides few PSFs for the Varian iX, we simulated the head through Monte Carlo (MC) simulations and calculated validated PSFs for 12 square field sizes including seven for small static fields. The resulting dosimetric calculations allowed us to reach a good level of agreement in comparison to our relative and absolute dose measurements performed on a Varian iX in water phantom. Measured and MC calculated output factors were investigated for different detectors. Based on the TRS 483 formalism and MC (PENELOPE/penEasy), we calculated output correction factors for the unshielded Diode-E (T60017) and the PinPoint-3D (T31016) micro-chamber according to manufacturers' blueprints. Our MC results were in agreement with the recommended data; they compete with recent measurements and MC simulations and in particular the TRS 483 MC data obtained from similar simulations. Moreover, our MC results provide supplemental data in comparison to TRS 483 data in particular for the PinPoint-3D (T31016). We suggest our MC output correction factors as new datasets for future TRS compilations. The work was substantial, used different robust MC strategies depending on the scoring regions, and led in most cases to uncertainties of less than 1%.

Reference dosimetry of modulated and dynamic photon beams

In the late 1980s, a new technique was proposed that would revolutionize radiotherapy. Now referred to as intensity-modulated radiotherapy, it is at the core of state-of-the-art photon beam delivery techniques, such as helical tomotherapy and volumetric modulated arc therapy. Despite over two decades of clinical application, there are still no established guidelines on the calibration of dynamic modulated photon beams. In 2008, the IAEA-AAPM work group on nonstandard photon beam dosimetry published a formalism to support the development of a new generation of protocols applicable to nonstandard beam reference dosimetry (Alfonso et al 2008 Med. Phys. 35 5179-86). The recent IAEA Code of Practice TRS-483 was published as a result of this initiative and addresses exclusively small static beams. But the plan-class specific reference calibration route proposed by Alfonso et al (2008 Med. Phys. 35 5179-86) is a change of paradigm that is yet to be implemented in radiotherapy clinics. The main goals of this paper are to provide a literature review on the dosimetry of nonstandard photon beams, including dynamic deliveries, and to discuss anticipated benefits and challenges in a future implementation of the IAEA-AAPM formalism on dynamic photon beams.

Calculated beam quality correction factors for ionization chambers in MV photon beams

The beam quality correction factor, [Formula: see text], which corrects for the difference in the ionization chamber response between the reference and clinical beam quality, is an integral part of radiation therapy dosimetry. The uncertainty of [Formula: see text] is one of the most significant sources of uncertainty in the dose determination. To improve the accuracy of available [Formula: see text] data, four partners calculated [Formula: see text] factors for 10 ionization chamber models in linear accelerator beams with accelerator voltages ranging from 6 MV to 25 MV, including flattening-filter-free (FFF) beams. The software used in the calculations were EGSnrc and PENELOPE, and the ICRU report 90 cross section data for water and graphite were included in the simulations. Volume averaging correction factors were calculated to correct for the dose averaging in the chamber cavities. A comparison calculation between partners showed a good agreement, as did comparison with literature. The [Formula: see text] values from TRS-398 were higher than our values for each chamber where data was available. The [Formula: see text] values for the FFF beams did not follow the same [Formula: see text], [Formula: see text] relation as beams with flattening filter (values for 10 MV FFF beams were below fits made to other data on average by 0.3%), although our FFF sources were only for Varian linacs.

A modified formalism for electron beam reference dosimetry to improve the accuracy of linac output calibration

PURPOSE

To present and demonstrate the accuracy of a modified formalism for electron beam reference dosimetry using updated Monte Carlo calculated beam quality conversion factors.

METHODS

The proposed, simplified formalism allows the use of cylindrical ionization chambers in all electron beams (even those with low beam energies) and does not require a measured gradient correction factor. Data from a previous publication are used for beam quality conversion factors. The formalism is tested and compared to the present formalism in the AAPM TG-51 protocol with measurements made in Elekta Precise electron beams with energies between 4 MeV and 22 MeV and with fields shaped with a 10 × 10 cm² clinical applicator as well as a 20 × 20 cm² clinical applicator for the 18 MeV and 22 MeV beams. A set of six ionization chambers are used for measurements (two cylindical reference-class chambers, two scanning-type chambers and two parallel-plate chambers). Dose per monitor unit is derived using the data and formalism provided in the TG-51 protocol and with the proposed formalism and data and compared to that obtained using ionization chambers calibrated directly against primary standards for absorbed dose in electron beams.

RESULTS

The standard deviation of results using different chambers when TG-51 is followed strictly is on the order of 0.4% when parallel-plate chambers are cross-calibrated against cylindrical chambers. However, if parallel-plate chambers are directly calibrated in a cobalt-60 beam, the difference between results for these chambers is up to 2.2%. Using the proposed formalism and either directly calibrated or cross-calibrated parallel-plate chambers gives a standard deviation using different chambers of 0.4%. The difference between results that use TG-51 and the primary standard measurements are on the order of 0.6% with a maximum difference in the 4 MeV beam of 2.8%. Comparing the results obtained with the proposed formalism and the primary standard measurements are on the order of 0.4% with a maximum difference of 1.0% in the 4 MeV beam.

CONCLUSIONS

The proposed formalism and the use of updated data for beam quality conversion factors improves the consistency of results obtained with different chamber types and improves the accuracy of reference dosimetry measurements. Moreover, it is simpler than the present formalism and will be straightforward to implement clinically.

Comparison of penh, fluka, and Geant4/topas for absorbed dose calculations in air cavities representing ionization chambers in high-energy photon and proton beams

PURPOSE

The purpose of this work is to analyze whether the Monte Carlo codes penh, fluka, and geant4/topas are suitable to calculate absorbed doses andf Q / f Q 0 ratios in therapeutic high-energy photon and proton beams.

METHODS

We used penh, fluka, geant4/topas, and egsnrc to calculate the absorbed dose to water in a reference water cavity and the absorbed dose to air in two air cavities representative of a plane-parallel and a cylindrical ionization chamber in a 1.25 MeV photon beam and a 150 MeV proton beam - egsnrc was only used for the photon beam calculations. The physics and transport settings in each code were adjusted to simulate the particle transport as detailed as reasonably possible. From these absorbed doses, f Q 0 factors, f Q factors, andf Q / f Q 0 ratios (which are the basis of Monte Carlo calculated beam quality correction factors k Q , Q 0 ) were calculated and compared between the codes. Additionally, we calculated the spectra of primary particles and secondary electrons in the reference water cavity, as well as the integrated depth-dose curve of 150 MeV protons in water.

RESULTS

The absorbed doses agreed within 1.4% or better between the individual codes for both the photon and proton simulations. The f Q 0 and f Q factors agreed within 0.5% or better for the individual codes for both beam qualities. The resultingf Q / f Q 0 ratios for 150 MeV protons agreed within 0.7% or better. For the 1.25 MeV photon beam, the spectra of photons and secondary electrons agreed almost perfectly. For the 150 MeV proton simulation, we observed differences in the spectra of secondary protons whereas the spectra of primary protons and low-energy delta electrons also agreed almost perfectly. The first 2 mm of the entrance channel of the 150 MeV proton Bragg curve agreed almost perfectly while for greater depths, the differences in the integrated dose were up to 1.5%.

CONCLUSION

penh, fluka, and geant4/topas are capable of calculating beam quality correction factors in proton beams. The differences in the f Q 0 and f Q factors between the codes are 0.5% at maximum. The differences in thef Q / f Q 0 ratios are 0.7% at maximum.

Impact of new ICRU Report 90 recommendations on calculated correction factors for reference dosimetry

In 2016 the ICRU published a new report dealing with key data for ionizing radiation dosimetry (ICRU Report 90). New recommendations have been made for the mean excitation energies I for air, graphite and liquid water as well as for the graphite density to use when evaluating the density effect. In addition, the ICRU Report 90 discusses renormalized photoelectric cross sections, but refuses to give a recommendation on the use of renormalization factors. However, the Consultative Committee for Ionizing Radiation recommends to use renormalized photoeffect cross sections. Goal of the present work is to evaluate the impact of these new recommendations on clinical reference dosimetry for high energy photon and electron beams. The beam quality correction factor k _Q was calculated by Monte Carlo simulations for compact and parallel plate ionization chambers. In case of photons seven phase space files from clinical accelerators and twelve spectra taken from literature from 4 MV to 24 MV and additionally a ⁶⁰Co source were applied. As electron source thirteen electron spectra available in literature were used in the range of 4 MeV-21 MeV. The new ICRU recommendations have a small impact on Monte Carlo calculated k _Q values for the chosen ionization chambers in the range of 0.1%-0.35% only-the difference increases for higher photon energies. The impact of the ICRU Report 90 recommendations on Monte Carlo calculated stopping power ratios s _w,a , perturbation factors p and beam quality correction factors k _Q was investigated and confirmed a decrese of s _w,a by a fraction of a percent for photon and electron beams. This study indicates that the impact of the new ICRU recommendation is within 0.35%. The determined deviations should be taken into account, when widely published Monte Carlo calculated values are examined.

Monte Carlo simulations in radiotherapy dosimetry

BACKGROUND

The use of the Monte Carlo (MC) method in radiotherapy dosimetry has increased almost exponentially in the last decades. Its widespread use in the field has converted this computer simulation technique in a common tool for reference and treatment planning dosimetry calculations.

METHODS

This work reviews the different MC calculations made on dosimetric quantities, like stopping-power ratios and perturbation correction factors required for reference ionization chamber dosimetry, as well as the fully realistic MC simulations currently available on clinical accelerators, detectors and patient treatment planning.

CONCLUSIONS

Issues are raised that include the necessity for consistency in the data throughout the entire dosimetry chain in reference dosimetry, and how Bragg-Gray theory breaks down for small photon fields. Both aspects are less critical for MC treatment planning applications, but there are important constraints like tissue characterization and its patient-to-patient variability, which together with the conversion between dose-to-water and dose-to-tissue, are analysed in detail. Although these constraints are common to all methods and algorithms used in different types of treatment planning systems, they make uncertainties involved in MC treatment planning to still remain "uncertain".

Technical Note: On the use of cylindrical ionization chambers for electron beam reference dosimetry

PURPOSE

To investigate the use of cylindrical chambers for electron beam dosimetry independent of energy by studying the variability of relative ion chamber perturbation corrections, one of the main concerns for electron beam dosimetry with cylindrical chambers.

METHODS

Measurements are made with sets of cylindrical and plane-parallel reference-class chambers as a function of depth in water in 8 MeV and 18 MeV electron beams. The ratio of chamber readings for similar chambers is normalized in a high-energy electron beam and can be thought of as relative perturbation corrections. Data are plotted as a function of mean electron energy at depth for a range of depths close to the phantom surface to R₈₀ , the depth at which the ionization falls to 80% of its maximum value. Additional, similar measurements are made in a Virtual Water® phantom with cylindrical chambers at the reference depth in a 4 MeV electron beam.

RESULTS

The variability of relative ion chamber perturbation corrections for nominally identical cylindrical Farmer-type chambers is found to be less than 0.4%, no worse than plane-parallel chambers with similar specifications.

CONCLUSIONS

This work discusses several issues related to the use of plane-parallel ion chambers and suggests that reference-class cylindrical chambers may be appropriate for reference dosimetry of all electron beams. This would simplify the reference dosimetry procedure and improve accuracy of beam calibration.

A GPU-accelerated Monte Carlo dose calculation platform and its application toward validating an MRI-guided radiation therapy beam model

PURPOSE

The clinical commissioning of IMRT subject to a magnetic field is challenging. The purpose of this work is to develop a GPU-accelerated Monte Carlo dose calculation platform based on penelope and then use the platform to validate a vendor-provided MRIdian head model toward quality assurance of clinical IMRT treatment plans subject to a 0.35 T magnetic field.

METHODS

penelope was first translated from fortran to c++ and the result was confirmed to produce equivalent results to the original code. The c++ code was then adapted to cuda in a workflow optimized for GPU architecture. The original code was expanded to include voxelized transport with Woodcock tracking, faster electron/positron propagation in a magnetic field, and several features that make gpenelope highly user-friendly. Moreover, the vendor-provided MRIdian head model was incorporated into the code in an effort to apply gpenelope as both an accurate and rapid dose validation system. A set of experimental measurements were performed on the MRIdian system to examine the accuracy of both the head model and gpenelope. Ultimately, gpenelope was applied toward independent validation of patient doses calculated by MRIdian's kmc.

RESULTS

An acceleration factor of 152 was achieved in comparison to the original single-thread fortran implementation with the original accuracy being preserved. For 16 treatment plans including stomach (4), lung (2), liver (3), adrenal gland (2), pancreas (2), spleen(1), mediastinum (1), and breast (1), the MRIdian dose calculation engine agrees with gpenelope with a mean gamma passing rate of 99.1% ± 0.6% (2%/2 mm).

CONCLUSIONS

A Monte Carlo simulation platform was developed based on a GPU- accelerated version of penelope. This platform was used to validate that both the vendor-provided head model and fast Monte Carlo engine used by the MRIdian system are accurate in modeling radiation transport in a patient using 2%/2 mm gamma criteria. Future applications of this platform will include dose validation and accumulation, IMRT optimization, and dosimetry system modeling for next generation MR-IGRT systems.

Accuracy of the electron transport in mcnp5 and its suitability for ionization chamber response simulations: A comparison with the egsnrc and penelope codes

PURPOSE

In this work, accuracy of the mcnp5 code in the electron transport calculations and its suitability for ionization chamber (IC) response simulations in photon beams are studied in comparison to egsnrc and penelope codes.

METHODS

The electron transport is studied by comparing the depth dose distributions in a water phantom subdivided into thin layers using incident energies (0.05, 0.1, 1, and 10 MeV) for the broad parallel electron beams. The IC response simulations are studied in water phantom in three dosimetric gas materials (air, argon, and methane based tissue equivalent gas) for photon beams ((60)Co source, 6 MV linear medical accelerator, and mono-energetic 2 MeV photon source). Two optional electron transport models of mcnp5 are evaluated: the ITS-based electron energy indexing (mcnp5(ITS)) and the new detailed electron energy-loss straggling logic (mcnp5(new)). The electron substep length (ESTEP parameter) dependency in mcnp5 is investigated as well.

RESULTS

For the electron beam studies, large discrepancies (>3%) are observed between the MCNP5 dose distributions and the reference codes at 1 MeV and lower energies. The discrepancy is especially notable for 0.1 and 0.05 MeV electron beams. The boundary crossing artifacts, which are well known for the mcnp5(ITS), are observed for the mcnp5(new) only at 0.1 and 0.05 MeV beam energies. If the excessive boundary crossing is eliminated by using single scoring cells, the mcnp5(ITS) provides dose distributions that agree better with the reference codes than mcnp5(new). The mcnp5 dose estimates for the gas cavity agree within 1% with the reference codes, if the mcnp5(ITS) is applied or electron substep length is set adequately for the gas in the cavity using the mcnp5(new). The mcnp5(new) results are found highly dependent on the chosen electron substep length and might lead up to 15% underestimation of the absorbed dose.

CONCLUSIONS

Since the mcnp5 electron transport calculations are not accurate at all energies and in every medium by general clinical standards, caution is needed, if mcnp5 is used with the current electron transport models for dosimetric applications.