Investigation of systematic uncertainties in Monte Carlo-calculated beam quality correction factors

Dose-response dependencies of OSL dosimeters in conventional linacs and 1.5T MR-linacs: an experimental and Monte Carlo study

Objective. This work focuses on the optically stimulated luminescence dosimetry (OSLD) dose-response characterization, with emphasis on 1.5T MR-Linacs.Approach. Throughout this study, the nanoDots OSLDs (Landauer, USA) were considered. In groups of three, the mean OSLD response was measured in a conventional linac and an MR-Linac under various irradiation conditions to investigate (i) dose-response linearity with and without the 1.5T magnetic field, (ii) signal fading rate and its dependencies, (iii) beam quality, detector orientation and dose rate dependencies in a conventional linac, (iii) potential MR imaging related effects on OSLD response and (iv) detector orientation dependence in an MR-Linac. Monte Carlo calculations were performed to further quantify angular dependence after rotating the detector around its central axis parallel to the magnetic field, and determine the magnetic field correction factors,kB,Q,for all cardinal detector orientations.Main results. OSLD dose-response supralinearity in an MR-Linac setting was found to agree within uncertainties with the corresponding one in a conventional linac, for the axial detector orientation investigated. Signal fading rate does not depend on irradiation conditions for the range of 3-30 d considered. OSLD angular (orientation) dependence is more pronounced under the presence of a magnetic field. OSLDs irradiated with and without real-time T2w MR imaging enabled during irradiation yielded the same response within uncertainties.kB,Qvalues were determined for all three cardinal orientations. Corrections needed reached up to 6.4%. However, if OSLDs are calibrated in the axial orientation and then irradiated in an MR-Linac placed again in the axial orientation (perpendicular to the magnetic field), then simulations suggest thatkB,Qcan be considered unity within uncertainties, irrespective of the incident beam angle.Significance. This work contributes towards OSLD dose-response characterization and relevant correction factors availability. OSLDs are suitable for QA checks in MR-based beam gating applications andin vivodosimetry in MR-Linacs.

Quantitative evaluation of dosimetric uncertainties in electron therapy by measurement and calculation using the electron Monte Carlo dose algorithm in the Eclipse treatment planning system

In the electron beam radiation therapy, customized blocks are mostly used to shape treatment fields to generate conformal doses. The goal of this study is to investigate quantitatively dosimetric uncertainties associated with heterogeneities, detectors used in the measurement of the beam data commissioning, and modeling of the interactions of high energy electrons with tissue. These uncertainties were investigated both by measurements with different detectors and calculations using electron Monte Carlo algorithm (eMC) in the Eclipse treatment planning system. Dose distributions for different field sizes were calculated using eMC and measured with a multiple-diode-array detector (MapCheck2) for cone sizes ranging from 6 to 25 cm. The dose distributions were calculated using the CT images of the MapCheck2 and water-equivalent phantoms. In the umbra region (<20% isodose line), the eMC underestimated dose by a factor of 3 for high energy electron beams due to lack of consideration of bremsstrahlung emitted laterally that was not accounted by eMC in the low dose region outside the field. In the penumbra (20%-80% isodose line), the eMC overestimated dose (40%) for high energy 20 MeV electrons compared to the measured dose with small diodes in the high gradient dose region. This was mainly due to lack of consideration of volume averaging of the ion chamber used in beam data commissioning which was input to the eMC dose calculation algorithm. Large uncertainties in the CT numbers (25%) resulted from the image artifacts in the CT images of the MapCheck2 phantom due to metal artifacts. The eMC algorithm used the electron and material densities extracted from the CT numbers which resulted large dosimetric uncertainties (10%) in the material densities and corresponding stopping power ratios. The dose calculations with eMC are associated with large uncertainties particularly in penumbra and umbra regions and around heterogeneities which affect the low dose level that cover nearby normal tissue or critical structures.

Monte Carlo calculation of detector perturbation and quality correction factors in a 1.5 T magnetic resonance guided radiation therapy small photon beams

Objective.With future advances in magnetic resonance imaging-guided radiation therapy, small photon beams are expected to be included regularly in clinical treatments. This study provides physical insights on detector dose-response to multiple megavoltage photon beam sizes coupled to magnetic fields and determines optimal orientations for measurements.Approach.Monte Carlo simulations determine small-cavity detector (solid-state: PTW60012 and PTW60019, ionization chambers: PTW31010, PTW31021, and PTW31022) dose-responses in water to an Elekta Unity 7 MV FFF photon beam. Investigations are performed for field widths between 0.25 and 10 cm in four detector axis orientations with respect to the 1.5 T magnetic field and the photon beam. The magnetic field effect on the overall perturbation factor (P_MC) accounting for the extracameral components, atomic composition, and density is quantified in each orientation. The density (P_ρ) and volume averaging (P_vol) perturbation factors and quality correction factors (kQB,QfB,f) accounting for the magnetic field are also calculated in each orientation.Main results.Results show thatP_volremains the most significant perturbation both with and without magnetic fields. In most cases, the magnetic field effect onP_volis 1% or less. The magnetic field effect onP_ρis more significant on ionization chambers than on solid-state detectors. This effect increases up to 1.564 ± 0.001 with decreasing field size for chambers. On the contrary, the magnetic field effect on the extracameral perturbation factor is higher on solid-state detectors than on ionization chambers. For chambers, the magnetic field effect onP_MCis only significant for field widths <1 cm, while, for solid-state detectors, this effect exhibits different trends with orientation, indicating that the beam incident angle and geometry play a crucial role.Significance.Solid-state detectors' dose-response is strongly affected by the magnetic field in all orientations. The magnetic field impact on ionization chamber response increases with decreasing field size. In general, ionization chambers yieldkQB,QfB,fcloser to unity, especially in orientations where the chamber axis is parallel to the magnetic field.

Properties of IBA Razor Nano Chamber in small-field radiation therapy using 6 MV FF, 6 MV FFF, and 10 MV FFF photon beams

BACKGROUND

Small megavoltage photon fields are increasingly used in modern radiotherapy techniques such as stereotactic radiotherapy. Therefore, it is important to study the reliability of dosimetry in the small-field conditions. The IBA Razor Nano Chamber (Nano chamber) ionization chamber is particularly intended for small-field measurements. In this work, properties of the Nano chamber were studied with both measurements and Monte Carlo (MC) simulations.

MATERIAL AND METHODS

The measurements and MC simulations were performed with 6 MV, 6 MV FFF and 10 MV FFF photon beams from the Varian TrueBeam linear accelerator. The source-to-surface distance was fixed at 100 cm. The measurements and MC simulations included profiles, percentage depth doses (PDD), and output factors (OF) in square jaw-collimated fields. The MC simulations were performed with the EGSnrc software system in a large water phantom.

RESULTS

The measured profiles and PDDs obtained with the Nano chamber were compared against IBA Razor Diode, PTW microDiamond and the PTW Semiflex ionization chamber. These results indicate that the Nano chamber is a high-resolution detector and thus suitable for small field profile measurements down to field sizes 2 × 2 cm² and appropriate for the PDD measurements. The field output correction factors kQclin, Qmsrfclin, fmsr and field OFs ΩQclin, Qmsrfclin, fmsr were determined according to TRS-483 protocol In the 6 MV FF and FFF beams, the determined correction factors kQclin, Qmsrfclin, fmsr were within 1.2% for the field sizes of 1 × 1 cm²-3 × 3 cm² and the experimental and MC defined field output factors ΩQclin,Qmsrfclin,fmsr showed good agreement.

CONCLUSION

The Nano chamber with its small cavity volume is a potential detector for the small-field dosimetry. In this study, the properties of this detector were characterized with measurements and MC simulations. The determined correction factors kQclin, Qmsrfclin, fmsr are novel results for the NC in the TrueBeam fields.

Monte Carlo-derived ionization chamber correction factors in therapeutic carbon ion beams

The accuracy of electromagnetic transport in the GEANT4 Monte Carlo (MC) code was investigated for carbon ion beams and ionization chamber (IC)-specific beam quality correction factors were calculated. This work implemented a Fano cavity test for carbon ion beams in the 100-450 MeV/u energy range to assess the accuracy of the default electromagnetic physics parameters. TheUrbanand theWentzel-VImultiple Coulomb scattering models were evaluated and the impact ofmaxStep,dRover,andfinal rangeparameters on the accuracy of the transport algorithm was investigated. The optimal production thresholds for an accurate calculation offQvalues, which is the product of the water-to-air stopping power ratio and the IC-specific perturbation correction factor, were also studied. ThefQcorrection factors were calculated for a cylindrical and a parallel-plate IC using carbon ions in the 150-450 MeV/u energy range. Modifying the default electromagnetic physics parameters resulted in a maximum deviation from theory of 0.3%. Therefore, the default EM parameters were used for the remainder of this work. ThefQfactors were found to converge for both ICs with decreasing production threshold distance below 5μm. ThefQvalues obtained in this work agreed with the TRS-398 stopping power ratios and other previously reported results within uncertainty. This study highlights an accurate MC-based technique to calculate the combined stopping power ratio and the perturbation correction factor for any IC in carbon ion beams.

Monte Carlo methods for device simulations in radiation therapy

Monte Carlo (MC) simulations play an important role in radiotherapy, especially as a method to evaluate physical properties that are either impossible or difficult to measure. For example, MC simulations (MCSs) are used to aid in the design of radiotherapy devices or to understand their properties. The aim of this article is to review the MC method for device simulations in radiation therapy. After a brief history of the MC method and popular codes in medical physics, we review applications of the MC method to model treatment heads for neutral and charged particle radiation therapy as well as specific in-room devices for imaging and therapy purposes. We conclude by discussing the impact that MCSs had in this field and the role of MC in future device design.

SURFACE DOSE ESTIMATION BY A KAP METER FOR KILOVOLTAGE X-RAY BEAMS

This study aims to estimate the entrance surface dose (ESD) of a water phantom for kilovoltage x-ray beams using an air kerma area product meter (KAP meter) equipped in an x-ray unit. The KAP meter was calibrated in terms of the ESD determined by a plane-parallel ionization chamber based on a 60Co absorbed dose-to-water calibration coefficient, ${N}_{D,w}^{{}^{60}\mathrm{C}\mathrm{o}}$. The ESD measured using the KAP meter was verified by comparing it with that estimated by the air kerma calibration coefficient, NK, for x-ray beam qualities. The ratio of ESDs based on ${N}_{D,w}^{{}^{60}\mathrm{C}\mathrm{o}}$ and NK was 1.003 on average and independent of the beam quality. The ESD by the KAP meter was an agreement within ±1.5% with that measured using the plane-parallel chamber for 10 × 10-30 × 30 cm2 fields with a source-surface distance of 75-150 cm. It was possible to estimate the ESD directly in a water phantom for x-ray beams without correction factors compared to the existing air kerma calibration, using a KAP meter calibrated based on ${N}_{D,w}^{{}^{60}\mathrm{C}\mathrm{o}}$.

Uncertainty in tissue equivalent proportional counter assessments of microdosimetry and RBE estimates in carbon radiotherapy

Microdosimetry is an important tool for assessing energy deposition distributions from ionizing radiation at cellular and cellular nucleus scales. It has served as an input parameter for multiple common mathematical models, including evaluation of relative biological effectiveness (RBE) of carbon ion therapy. The most common detector used for microdosimetry is the tissue-equivalent proportional counter (TEPC). Although it is widely applied, TEPC has various inherent uncertainties. Therefore, this work quantified the magnitude of TEPC measurement uncertainties and their impact on RBE estimates for therapeutic carbon beams. Microdosimetric spectra and frequency-, dose-, and saturation-corrected dose-mean lineal energy (****) were calculated using the Monte Carlo toolkit Geant4 for five monoenergetic and three spread-out Bragg peak carbon beams in water at every millimeter along the central beam axis. We simulated the following influences on these spectra from eight sources of uncertainty: wall effects, pulse pile-up, electronics, gas pressure, W-value, gain instability, low energy cut-off, and counting statistics. Statistic uncertainty was quantified as the standard deviation of perturbed values for each source. Bias was quantified as the difference between default lineal energy values and the mean of perturbed values for each systematic source. Uncertainties were propagated to RBE using the modified microdosimetric kinetic model (MKM). Variance introduced by statistic sources iny¯Fandy¯Daveraged 3.8% and 3.4%, respectively, and 1.5% iny*across beam depths and energies. Bias averaged 6.2% and 7.3% iny¯Fandy¯D,and 4.8% iny*.These uncertainties corresponded to 1.2 ± 0.9% on average in RBE_MKM. The largest contributors to variance and bias were pulse pile-up and wall effects. This study established an error budget for microdosimetric carbon measurements by quantifying uncertainty inherent to TEPC measurements. It is necessary to understand how robust the measurement of RBE model input parameters are against this uncertainty in order to verify clinical model implementation.

The dose response of PTW microDiamond and microSilicon in transverse magnetic field under small field conditions

The aim of the present work is to investigate the behavior of two diode-type detectors (PTW microDiamond 60019 and PTW microSilicon 60023) in transverse magnetic field under small field conditions. A formalism based on TRS 483 has been proposed serving as the framework for the application of these high-resolution detectors under these conditions. Measurements were performed at the National Metrology Institute of Germany (PTB, Braunschweig) using a research clinical linear accelerator facility. Quadratic fields corresponding to equivalent square field sizesSbetween 0.63 and 4.27 cm at the depth of measurement were used. The magnetic field strength was varied up to 1.4 T. Experimental results have been complemented with Monte Carlo simulations up to 1.5 T. Detailed simulations were performed to quantify the small field perturbation effects and the influence of detector components on the dose response. The does response of both detectors decreases by up to 10% at 1.5 T in the largest field size investigated. InS = 0.63 cm, this reduction at 1.5 T is only about half of that observed in field sizesS > 2 cm for both detectors. The results of the Monte Carlo simulations show agreement better than 1% for all investigated conditions. Due to normalization at the machine specific reference field, the resulting small field output correction factors for both detectors in magnetic fieldkQclin,QmsrBare smaller than those in the magnetic field-free case, where correction up to 6.2% at 1.5 T is required for the microSilicon in the smallest field size investigated. The volume-averaging effect of both detectors was shown to be nearly independent of the magnetic field. The influence of the enhanced-density components within the detectors has been identified as the major contributors to their behaviors in magnetic field. Nevertheless, the effect becomes weaker with decreasing field size that may be partially attributed to the deficiency of low energy secondary electrons originated from distant locations in small fields.

Report of AAPM Task Group 155: Megavoltage photon beam dosimetry in small fields and non-equilibrium conditions

Small-field dosimetry used in advance treatment technologies poses challenges due to loss of lateral charged particle equilibrium (LCPE), occlusion of the primary photon source, and the limited choice of suitable radiation detectors. These challenges greatly influence dosimetric accuracy. Many high-profile radiation incidents have demonstrated a poor understanding of appropriate methodology for small-field dosimetry. These incidents are a cause for concern because the use of small fields in various specialized radiation treatment techniques continues to grow rapidly. Reference and relative dosimetry in small and composite fields are the subject of the International Atomic Energy Agency (IAEA) dosimetry code of practice that has been published as TRS-483 and an AAPM summary publication (IAEA TRS 483; Dosimetry of small static fields used in external beam radiotherapy: An IAEA/AAPM International Code of Practice for reference and relative dose determination, Technical Report Series No. 483; Palmans et al., Med Phys 45(11):e1123, 2018). The charge of AAPM task group 155 (TG-155) is to summarize current knowledge on small-field dosimetry and to provide recommendations of best practices for relative dose determination in small megavoltage photon beams. An overview of the issue of LCPE and the changes in photon beam perturbations with decreasing field size is provided. Recommendations are included on appropriate detector systems and measurement methodologies. Existing published data on dosimetric parameters in small photon fields (e.g., percentage depth dose, tissue phantom ratio/tissue maximum ratio, off-axis ratios, and field output factors) together with the necessary perturbation corrections for various detectors are reviewed. A discussion on errors and an uncertainty analysis in measurements is provided. The design of beam models in treatment planning systems to simulate small fields necessitates special attention on the influence of the primary beam source and collimating devices in the computation of energy fluence and dose. The general requirements for fluence and dose calculation engines suitable for modeling dose in small fields are reviewed. Implementations in commercial treatment planning systems vary widely, and the aims of this report are to provide insight for the medical physicist and guidance to developers of beams models for radiotherapy treatment planning systems.

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.

Experimental determination of magnetic field correction factors for ionization chambers in parallel and perpendicular orientations

Magnetic field correction factors are needed for absolute dosimetry in magnetic resonance (MR)-linacs. Currently experimental data for magnetic field correction factors, especially for small volume ionization chambers, are largely lacking. The purpose of this work is to establish, independent methods for the experimental determination of magnetic field correction factors [Formula: see text] in an orientation in which the ionization chamber is parallel to the magnetic field. The aim is to confirm previous experiments on the determination of Farmer type ionization chamber correction factors and to gather information about the usability of small-volume ionization chambers for absolute dosimetry in MR-linacs. The first approach to determine [Formula: see text] is based on a cross-calibration of measurements using a conventional linac with an electromagnet and an MR-linac. The absolute influence of the magnetic field in perpendicular orientation is quantified with the help of the conventional linac and the electromagnet. The correction factors for the parallel orientation are then derived by combining these measurements with relative measurements in the MR-linac. The second technique utilizes alanine electron paramagnetic resonance dosimetry. The alanine system as well as several ionization chambers were directly calibrated with the German primary standard for absorbed dose to water. Magnetic field correction factors for the ionization chambers were determined by a cross-calibration with the alanine in an MR-linac. Important quantities like [Formula: see text] for Farmer type ionization chambers in parallel orientation and the change of the dose to water due the magnetic field [Formula: see text] have been confirmed. In addition, magnetic field correction factors have been determined for small volume ionization chambers in parallel orientation. The electromagnet-based measurements of [Formula: see text] for [Formula: see text] MR-linacs and parallel ionization chamber orientations resulted in 0.9926(22), 0.9935(31) and 0.9841(27) for the PTW 30013, the PTW 31010 and the PTW 31021, respectively. The measurements based on the second technique resulted in values for [Formula: see text] of 0.9901(72), 0.9955(72), and 0.9885(71). Both methods show excellent accuracy and reproducibility and are therefore suitable for the determination of magnetic field correction factors. Small-volume ionization chambers showed a variation in the resulting values for [Formula: see text] and should be cross-calibrated instead of using tabulated values for correction factors.

Small-cavity chamber dose response in megavoltage photon beams coupled to magnetic fields

In MRgRT, dosimetry measurements are performed in the presence of magnetic fields. For high-resolution measurements, small-cavity ionization chambers are required. While Monte Carlo simulations are essential to determine dosimetry correction factors, models of small-chambers require careful validation with experimental measurements. The aim of this study is to characterize small-cavity chamber response coupled to magnetic fields. Small-cavity chambers (PTW31010, PTW31016, PTW31021 and PTW3022) are irradiated by a 6 MV photon beam for 9 magnetic field strengths between -1.5 T and +1.5 T. The chamber axis is orientated either parallel or perpendicular to the irradiation beam, with the magnetic field always perpendicular to the beam. MC simulations are performed in EGSnrc. The sensitive volume of the chambers is reduced to account for the inefficiency adjacent to the guard electrode (dead volume) based on COMSOL calculations of electric potentials. The magnetic field affects the chamber response by up to 4.1% and 4.5% in the parallel and perpendicular orientations, respectively, compared to no magnetic field. The maximal difference in dose response between experiments and simulations is up to 6.1% and 4.5% for parallel and perpendicular orientation, respectively. When the dead volume is removed, which accounts for the 15%-23% of the nominal volume, the difference, in most cases, is within the stated uncertainties. Nevertheless, for a particular chamber, the reduced nominal volume barely improved the agreement between the experimental and calculated relative response (4.53% to 4.13%). This disagreement may be due to the imperfect chamber geometry model, as was found from microCT images. A detailed uncertainty analysis is presented. The characterization of small-cavity ion chamber response coupled to magnetic fields is complex. Small differences between real and model chamber geometry that normally would be insignificant become an issue in the presence of magnetic fields. Accurate characterization of the nominal volume is essential for small-cavity ion chamber modelling.

Determination of the absorbed dose to water for medium-energy x-ray beams using Fricke dosimetry

PURPOSE

For x-ray beams in the low and medium energy range, reference dosimetry is established in terms of air kerma. Fricke dosimetry has shown great potential in the absolute measurements of the absorbed dose to water for high-energy ranges. Therefore, the main purpose of this work was to compare the absorbed dose to water for medium-energy x-ray beams obtained through Fricke dosimetry with that obtained from the air kerma rate.

METHODS

To determine the absorbed dose to water using Fricke dosimetry, the polyethylene bags methodology was chosen. Fricke solution was irradiated at four different beam qualities. The absorbed dose to water values obtained using Fricke dosimetry were compared to those obtained using the standard protocol, using the Z-score.

RESULTS

Values of the Z-score were <2 for all measurements of absorbed dose to water, which means that the values obtained using Fricke dosimetry are equivalent to those obtained using the reference protocol. The combined standard uncertainty for the absorbed dose to water obtained by Fricke dosimetry was lower than that obtained with the ionization chamber.

CONCLUSIONS

Chemical dosimetry using a standard FeSO₄ solution has been demonstrated to be a potential option as a standard for the quantity absorbed dose to water for medium kV x-ray qualities.

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.

Effect of ICRU report 90 recommendations on Monte Carlo calculated kQ for ionization chambers listed in the Addendum to AAPM's TG-51 protocol

PURPOSE

The ICRU has published new recommendations for ionizing radiation dosimetry. In this work, the effect of recommendations on the water-to-air and graphite-to-air restricted mass electronic stopping power ratios (s_{w, air} and s_{g, air} ) and the individual perturbation correction factors P_i was calculated. The effect on the beam quality conversion factors k_Q for reference dosimetry of high-energy photon beams was estimated for all ionization chambers listed in the Addendum to AAPM's TG-51 protocol.

METHODS

The s_{w, air} , s_{g, air} , individual P_i, and k_Q were calculated using EGSnrc Monte Carlo code system and key data of both ICRU report 37 and ICRU report 90. First, the P_i and k_Q were calculated using precise models of eight ionization chambers: NE2571 (Nuclear Enterprise), 30013, 31010, 31021 (PTW), Exradin A12, A12S, A1SL (Standard imaging), and FC-65P (IBA). In this simulation, the radiation sources were one ⁶⁰ Co beam and ten photon beams with nominal energy between 4 MV and 25 MV. Then, the change in k_Q for ionization chambers listed in the Addendum to AAPM's TG-51 protocol was calculated by changing the specification of the simple-model of ionization chamber. The simple-models were made with only cylindrical component modules. In this simulation, the radiation sources of ⁶⁰ Co beam and 24 MV photon beam were used.

RESULTS

The significant changes (p < 0.05) were observed for s_{w, air} , s_{g, air} , the wall correction factor P_wall , and the waterproofing sleeve correction factor P_sleeve . The decrease in s_{w, air} varied from -0.57% for a ⁶⁰ Co beam to -0.36% for the highest beam quality. The decrease in s_{g, air} varied from -0.72% to -1.12% in the same range. The changes in P_wall and P_sleeve were up to 0.41% and 0.14% and those maximum changes were observed for the ⁶⁰ Co beam. All changes in the central electrode correction factor P_cel , the stem correction factor P_stem , and the replacement correction factor P_repl were from -0.02% to 0.12%. Those changes were statistically insignificant (p = 0.07 or more) and were independent of photon energy. The change in k_Q was mainly characterized by the change in s_{w, air} , P_wall , and P_sleeve . The relationship between the change in k_Q and the beam quality index was linear approximately. The changes in k_Q of the simple-models were agreed with those of the precise-models within 0.08%.

CONCLUSION

The effects of ICRU-90 recommendations on k_Q for the ionization chambers listed in the Addendum to AAPM's TG-51 protocol were from -0.15% to 0.30%. To remove the known systematic effect on the clinical reference dosimetry, the k_Q based on ICRU-37 should be updated to the k_Q based on ICRU-90.

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.

On the impact of ICRU report 90 recommendations on kQ factors for high-energy photon beams

PURPOSE

To assess the impact of the ICRU report 90 recommendations on the beam-quality conversion factor, kQ , used for clinical reference dosimetry of megavoltage linac photon beams.

METHODS

The absorbed dose to water and the absorbed dose to the air in ionization chambers representative of those typically used for linac photon reference dosimetry are calculated at the reference depth in a water phantom using Monte Carlo simulations. Depth-dose calculations in water are also performed to investigate changes in beam quality specifiers. The calculations are performed in a cobalt-60 beam and MV photon beams with nominal energy between 6 MV and 25 MV using the EGSnrc simulation toolkit. Inputs to the calculations use stopping-power data for graphite and water from the original ICRU-37 report and the new proposed values from the recently published ICRU-90 report. Calculated kQ factors are compared using the two different recommendations for key dosimetry data and measured kQ factors.

RESULTS

Less than about 0.1% effects from ICRU-90 recommendations on the beam quality specifiers, the photon component of the percentage depth-dose at 10 cm, %dd(10)x , and the tissue-phantom ratio at 20 cm and 10 cm, TPR1020, are observed. Although using different recommendations for key dosimetric data impact water-to-air stopping-power ratios and ion chamber perturbation corrections by up to 0.54% and 0.40%, respectively, we observe little difference (≤0.14%) in calculated kQ factors. This is contradictory to the predictions in ICRU-90 that suggest differences up to 0.5% in high-energy photon beams. A slightly better agreement with experimental values is obtained when using ICRU-90 recommendations.

CONCLUSION

Users of the addendum to the TG-51 protocol for reference dosimetry of high-energy photon beams, which recommends Monte Carlo calculated kQ factors, can rest assured that the recommendations of ICRU report 90 on basic data have little impact on this central dosimetric parameter.

Monte Carlo and experimental determination of correction factors for gamma knife perfexion small field dosimetry measurements

Detector-, field size- and machine-specific correction factors are required for precise dosimetry measurements in small and non-standard photon fields. In this work, Monte Carlo (MC) simulation techniques were used to calculate the [Formula: see text] and [Formula: see text] correction factors for a series of ionization chambers, a synthetic microDiamond and diode dosimeters, used for reference and/or output factor (OF) measurements in the Gamma Knife Perfexion photon fields. Calculations were performed for the solid water (SW) and ABS plastic phantoms, as well as for a water phantom of the same geometry. MC calculations for the [Formula: see text] correction factors in SW were compared against corresponding experimental results for a subset of ionization chambers and diode detectors. Reference experimental OF data were obtained through the weighted average of corresponding measurements using TLDs, EBT-2 films and alanine pellets. [Formula: see text] values close to unity (within 1%) were calculated for most of ionization chambers in water. Greater corrections of up to 6.0% were observed for chambers with relatively large air-cavity dimensions and steel central electrode. A phantom correction of 1.006 and 1.024 (breaking down to 1.014 from the ABS sphere and 1.010 from the accompanying ABS phantom adapter) were calculated for the SW and ABS phantoms, respectively, adding up to [Formula: see text] corrections in water. Both measurements and MC calculations for the diode and microDiamond detectors resulted in lower than unit [Formula: see text] correction factors, due to their denser sensitive volume and encapsulation materials. In comparison, higher than unit [Formula: see text] results for the ionization chambers suggested field size depended dose underestimations (being significant for the 4 mm field), with magnitude depending on the combination of contradicting phenomena associated with volume averaging and electron fluence perturbations. Finally, the presence of 0.5 mm air-gap between the diodes' frontal surface and their phantom-inserts may considerably influence OF measurements, reaching 4.6% for the Razor diode.

Detector dose response in megavoltage small photon beams. II. Pencil beam perturbation effects

PURPOSE

To quantify detector perturbation effects in megavoltage small photon fields and support the theoretical explanation on the nature of quality correction factors in these conditions.

METHODS

In this second paper, a modern approach to radiation dosimetry is defined for any detector and applied to small photon fields. Fano's theorem is adapted in the form of a cavity theory and applied in the context of nonstandard beams to express four main effects in the form of perturbation factors. The pencil-beam decomposition method is detailed and adapted to the calculation of perturbation factors and quality correction factors. The approach defines a perturbation function which, for a given field size or beam modulation, entirely determines these dosimetric factors. Monte Carlo calculations are performed in different cavity sizes for different detection materials, electron densities, and extracameral components.

RESULTS

Perturbation effects are detailed with calculated perturbation functions, showing the relative magnitude of the effects as well as the geometrical extent to which collimating or modulating the beam impacts the dosimetric factors. The existence of a perturbation zone around the detector cavity is demonstrated and the approach is discussed and linked to previous approaches in the literature to determine critical field sizes.

CONCLUSIONS

Monte Carlo simulations are valuable to describe pencil beam perturbation effects and detail the nature of dosimetric factors in megavoltage small photon fields. In practice, it is shown that dosimetric factors could be avoided if the field size remains larger than the detector perturbation zone. However, given a detector and beam quality, a full account for the detector geometry is necessary to determine critical field sizes.

On the development of a comprehensive MC simulation model for the Gamma Knife Perfexion radiosurgery unit

This work presents a comprehensive Monte Carlo (MC) simulation model for the Gamma Knife Perfexion (PFX) radiosurgery unit. Model-based dosimetry calculations were benchmarked in terms of relative dose profiles (RDPs) and output factors (OFs), against corresponding EBT2 measurements. To reduce the rather prolonged computational time associated with the comprehensive PFX model MC simulations, two approximations were explored and evaluated on the grounds of dosimetric accuracy. The first consists in directional biasing of the (60)Co photon emission while the second refers to the implementation of simplified source geometric models. The effect of the dose scoring volume dimensions in OF calculations accuracy was also explored. RDP calculations for the comprehensive PFX model were found to be in agreement with corresponding EBT2 measurements. Output factors of 0.819  ±  0.004 and 0.8941  ±  0.0013 were calculated for the 4 mm and 8 mm collimator, respectively, which agree, within uncertainties, with corresponding EBT2 measurements and published experimental data. Volume averaging was found to affect OF results by more than 0.3% for scoring volume radii greater than 0.5 mm and 1.4 mm for the 4 mm and 8 mm collimators, respectively. Directional biasing of photon emission resulted in a time efficiency gain factor of up to 210 with respect to the isotropic photon emission. Although no considerable effect on relative dose profiles was detected, directional biasing led to OF overestimations which were more pronounced for the 4 mm collimator and increased with decreasing emission cone half-angle, reaching up to 6% for a 5° angle. Implementation of simplified source models revealed that omitting the sources' stainless steel capsule significantly affects both OF results and relative dose profiles, while the aluminum-based bushing did not exhibit considerable dosimetric effect. In conclusion, the results of this work suggest that any PFX simulation model should be benchmarked in terms of both RDP and OF results.

[Benchmark experiment to verify radiation transport calculations for dosimetry in radiation therapy]

Monte Carlo simulations are regarded as the most accurate method of solving complex problems in the field of dosimetry and radiation transport. In (external) radiation therapy they are increasingly used for the calculation of dose distributions during treatment planning. In comparison to other algorithms for the calculation of dose distributions, Monte Carlo methods have the capability of improving the accuracy of dose calculations - especially under complex circumstances (e.g. consideration of inhomogeneities). However, there is a lack of knowledge of how accurate the results of Monte Carlo calculations are on an absolute basis. A practical verification of the calculations can be performed by direct comparison with the results of a benchmark experiment. This work presents such a benchmark experiment and compares its results (with detailed consideration of measurement uncertainty) with the results of Monte Carlo calculations using the well-established Monte Carlo code EGSnrc. The experiment was designed to have parallels to external beam radiation therapy with respect to the type and energy of the radiation, the materials used and the kind of dose measurement. Because the properties of the beam have to be well known in order to compare the results of the experiment and the simulation on an absolute basis, the benchmark experiment was performed using the research electron accelerator of the Physikalisch-Technische Bundesanstalt (PTB), whose beam was accurately characterized in advance. The benchmark experiment and the corresponding Monte Carlo simulations were carried out for two different types of ionization chambers and the results were compared. Considering the uncertainty, which is about 0.7 % for the experimental values and about 1.0 % for the Monte Carlo simulation, the results of the simulation and the experiment coincide.

Towards a quantitative, measurement-based estimate of the uncertainty in photon mass attenuation coefficients at radiation therapy energies

In this study, a quantitative estimate is derived for the uncertainty in the XCOM photon mass attenuation coefficients in the energy range of interest to external beam radiation therapy-i.e. 100 keV (orthovoltage) to 25 MeV-using direct comparisons of experimental data against Monte Carlo models and theoretical XCOM data. Two independent datasets are used. The first dataset is from our recent transmission measurements and the corresponding EGSnrc calculations (Ali et al 2012 Med. Phys. 39 5990-6003) for 10-30 MV photon beams from the research linac at the National Research Council Canada. The attenuators are graphite and lead, with a total of 140 data points and an experimental uncertainty of ∼0.5% (k = 1). An optimum energy-independent cross section scaling factor that minimizes the discrepancies between measurements and calculations is used to deduce cross section uncertainty. The second dataset is from the aggregate of cross section measurements in the literature for graphite and lead (49 experiments, 288 data points). The dataset is compared to the sum of the XCOM data plus the IAEA photonuclear data. Again, an optimum energy-independent cross section scaling factor is used to deduce the cross section uncertainty. Using the average result from the two datasets, the energy-independent cross section uncertainty estimate is 0.5% (68% confidence) and 0.7% (95% confidence). The potential for energy-dependent errors is discussed. Photon cross section uncertainty is shown to be smaller than the current qualitative 'envelope of uncertainty' of the order of 1-2%, as given by Hubbell (1999 Phys. Med. Biol 44 R1-22).

A feasibility study of Fricke dosimetry as an absorbed dose to water standard for 192Ir HDR sources

High dose rate brachytherapy (HDR) using ¹⁹²Ir sources is well accepted as an important treatment option and thus requires an accurate dosimetry standard. However, a dosimetry standard for the direct measurement of the absolute dose to water for this particular source type is currently not available. An improved standard for the absorbed dose to water based on Fricke dosimetry of HDR ¹⁹²Ir brachytherapy sources is presented in this study. The main goal of this paper is to demonstrate the potential usefulness of the Fricke dosimetry technique for the standardization of the quantity absorbed dose to water for ¹⁹²Ir sources. A molded, double-walled, spherical vessel for water containing the Fricke solution was constructed based on the Fricke system. The authors measured the absorbed dose to water and compared it with the doses calculated using the AAPM TG-43 report. The overall combined uncertainty associated with the measurements using Fricke dosimetry was 1.4% for k = 1, which is better than the uncertainties reported in previous studies. These results are promising; hence, the use of Fricke dosimetry to measure the absorbed dose to water as a standard for HDR ¹⁹²Ir may be possible in the future.

Addendum to the AAPM's TG-51 protocol for clinical reference dosimetry of high-energy photon beams

An addendum to the AAPM's TG-51 protocol for the determination of absorbed dose to water in megavoltage photon beams is presented. This addendum continues the procedure laid out in TG-51 but new kQ data for photon beams, based on Monte Carlo simulations, are presented and recommendations are given to improve the accuracy and consistency of the protocol's implementation. The components of the uncertainty budget in determining absorbed dose to water at the reference point are introduced and the magnitude of each component discussed. Finally, the consistency of experimental determination of ND,w coefficients is discussed. It is expected that the implementation of this addendum will be straightforward, assuming that the user is already familiar with TG-51. The changes introduced by this report are generally minor, although new recommendations could result in procedural changes for individual users. It is expected that the effort on the medical physicist's part to implement this addendum will not be significant and could be done as part of the annual linac calibration.

Correction factors for source strength determination in HDR brachytherapy using the in-phantom method

For the purpose of clinical source strength determination for HDR brachytherapy sources, the German society for Medical Physics (DGMP) recommends in their report 13 the usage of a solid state phantom (Krieger-phantom) with a thimble ionization chamber. In this work, the calibration chain for the determination of the reference air-kerma rate Ka,100 and reference dose rate to waterDw,1 by ionization chamber measurement in the Krieger-phantom was modeled via Monte Carlo simulations. These calculations were used to determine global correction factors k(tot), which allows a user to directly convert the reading of an ionization chamber calibrated in terms of absorbed dose to water, into the desired quantity Ka,100 or Dw,1. The factor k(tot) was determined for four available (192)Ir sources and one (60)Co source with three different thimble ionization chambers. Finally, ionization chamber measurements on three μSelectron V2 HDR sources within the Krieger-phantom were performed and Ka,100 was determined according to three different methods: 1) using a calibration factor in terms of absorbed dose to water with the global correction factor [Formula: see text] according DGMP 13 2) using a global correction factor calculated via Monte Carlo 3) using a direct reference air-kerma rate calibration factor determined by the national metrology institute PTB. The comparison of Monte Carlo based [Formula: see text] with those from DGMP 13 showed that the DGMP data were systematically smaller by about 2-2.5%. The experimentally determined [Formula: see text] , based on the direct Ka,100 calibration were also systematically smaller by about 1.5%. Despite of these systematical deviations, the agreement of the different methods was in almost all cases within the 1σ level of confidence of the interval of their respective uncertainties in a Gaussian distribution. The application of Monte Carlo based [Formula: see text] for the determination of Ka,100 for three μSelectron V2 sources revealed the smallest deviation to the manufacturer's source certificate. With the calculated [Formula: see text] for a (60)Co source, the user is now able to accurately determine Ka,100 of a HDR (60)Co source via in-phantom measurement. Moreover, using the presented global correction factor [Formula: see text] , the user is able to determine the future source specification quantity Dw,1 with the same in-phantom setup.