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

Evaluation of radiation detectors for the determination of field output factors in Leksell Gamma Knife dosimetry using 3D printed phantom inserts

The Leksell Gamma Knife® Perfexion™ and Icon™ have a unique geometry, containing 192 60Co sources with collimation for field sizes of 4 mm, 8 mm, and 16 mm. 4 mm and 8 mm collimated fields lack lateral charged particle equilibrium, so accurate field output factors are essential. This study performs field output factor measurements for the microDiamond, microSilicon, and RAZOR™ Nano detectors. 3D printed inserts for the spherical Solid Water® Phantom were fabricated for microDiamond detector, the microSilicon unshielded diode and the RAZOR™ Nano micro-ionisation chamber. Detectors were moved iteratively to identify the peak detector signal for each collimator, representing the effective point of measurement of the chamber. In addition, field output correction factors were calculated for each detector relative to vendor supplied Monte Carlo simulated field output factors and field output factors measured with a W2 scintillator. All field output factors where within 1.1 % for the 4 mm collimator and within 2.3 % for the 8 mm collimator. The 3D printed phantom inserts were suitable for routine measurements if the user identifies the effective point of measurement, and ensures a reproducible setup by marking the rotational alignment of the cylindrical print. Measurements with the microDiamond and microSilicon can be performed faster compared to the RAZOR™ Nano due to differences in the signal to noise ratio. All detectors are suitable for field output factor measurements for the Leksell Gamma Knife® Perfexion™ and Icon™.

Feasibility of Isodose-shaped scintillation detectors for the measurement of gamma Knife® output factors

PURPOSE

Scintillation detectors were 3D printed based on a gamma knife (GK) dose distribution to calculate the volume averaging effect. The collimator output factors were measured using isodose-shaped scintillators (ISSs) and compared with those of a micro-diamond detector and previous reports.

METHODS

An absorbed dose distribution in a spherical dosimetry phantom with a radius of 8 cm was obtained from GK treatment planning software (Leksell GammaPlan (LGP), Elekta AB, Stockholm, Sweden). Two types of ISSs were fabricated to fit the 97.2% (ISS-1) and 95.6% (ISS-2) isodose surfaces. The volume averaging correction factors were obtained by dividing the absorbed dose to water in the central voxel (CV) by that in the ISS. The correction effect due to the difference between the ISS and water was calculated by Monte Carlo simulations. Ten ISS detectors, five of each type, were used to measure the output factors of the 4 and 8 mm collimators of a GK Icon^TM to assess system consistency. The output factors of seven GKs were measured using two ISS detectors, one of each type, and a PTW T60019 (PTW, Freiburg, Germany) micro-diamond detector.

RESULTS

The detector output ratios (DORs) measured using the five ISSs of each type were consistent, with standard uncertainties less than 0.2%. In the 4 mm field, the volume averaging correction factor ratios were 1.018 and 1.026, and the output factors after all corrections were 0.827 (0.006) and 0.825 (0.006) for ISS-1 and ISS-2, respectively. In the 8 mm field, the volume averaging correction factor ratios were 1.000 for both ISS types, and the output factors were 0.898 (0.003) and 0.900 (0.003) for ISS-1 and ISS-2, respectively. The ISS detectors could measure the output factors of a GK with uncertainties comparable to that of the PTW 60019 detector. The output factors of all detectors decreased with the dose rate.

CONCLUSION

The volume averaging effect of an ISS developed in-house could be calculated using known dose distributions. The collimator output factors of the GK Perfexion/Icon™ models measured using ISS detectors were consistent with those of a commercial synthetic micro-diamond detector and recent studies. This article is protected by copyright. All rights reserved.

Dose kernel decomposition for spot-based radiotherapy treatment planning

PURPOSE

Pre-calculation of accurate dose deposition kernels for treatment planning of spot-based radiotherapies, such as Gamma Knife (GK) and Gamma Pod (GP), can be very time-consuming and may require large data storage with an enormous number of possible spots. We proposed a novel kernel decomposition (KD) model to address accurate and fast (real-time) dose calculation with reduced data storage requirements for spot-based treatment planning. The application of the KD model was demonstrated for clinical GK and GP radiotherapy platforms.

METHODS

The dose deposition kernel at each spot (shot position) is modeled as the product of a shift-invariant kernel based on a reference kernel and spatially variant scale factor. The reference kernel, one for each collimator, is defined at the center of the commissioning phantom for GK and at the center of the treatment target for GP and calculated using the Monte Carlo (MC) method. The spatially variant scale factor is defined as the ratio of the mean tissue maximum ratio (TMR) at the candidate shot position to that at the reference kernel position, and the mean TMR map is calculated within the entire volume through parallel beam ray tracing on the density image followed by averaging over all source directions. The proposed KD dose calculations were compared with the MC method and with the GK and GP treatment planning system (TPS) computations for various shot positions and collimator sizes utilizing a phantom and 14 and 12 clinical plans for GK and GP, respectively.

RESULTS

For the phantom study, the KD Gamma index (3%/1 mm) passing rates were greater than 99% (median 100%) relative to the MC doses, except for the shots close to the boundary. The passing rates dropped below 90% for 8 mm (16 mm) shots positioned within ∼1 cm (∼2 cm) of the boundary. For the clinical GK plans, the KD Gamma passing rates were greater than 99% (median 100%) compared to the MC and greater than 92% (median 99%) compared to the TPS. For the clinical GP plans, the KD Gamma passing rates were greater than 95% (median 98%) compared to the MC and greater than 91% (median 97%) compared to the TPS. The scale factors were calculated in sub-seconds with GPU implementation and only need to be calculated once before treatment plan optimization. The calculation of the dose kernel was also within sub-seconds without requiring beam-by-beam calculation commonly done in the TPS.

CONCLUSION

The proposed model can provide an accurate dose and enables real-time dose and derivative calculations by kernel shifting and scaling without pre-calculating or requiring large data storage for GK and GP dose deposition kernels during treatment planning. This model could be useful for spot-based radiotherapy treatment planning by allowing an efficient global fine search for optimal spots.

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.

Small-field beam data acquisition, detector dependency, and film-based validation for a novel self-shielded stereotactic radiosurgery system

PURPOSE

This study reports a single-institution experience with beam data acquisition and film-based validation for a novel self-shielded sterotactic radiosurgery unit and investigates detector dependency on field output factors (OFs), off-axis ratios (OARs), and percent depth dose (PDD) measurements within the context of small-field dosimetry.

METHODS

The delivery platform for this unit consists of a 2.7-MV S-band linear accelerator mounted on coupled gimbals that rotate around a common isocenter (source-to-axis distance [SAD] = 450 mm), allowing for more than 260 noncoplanar beam angles. Beam collimation is achieved via a tungsten collimator wheel with eight circular apertures ranging from 4 mm to 25 mm in diameter. Three diodes (PTW 60012 Diode E, PTW 60018 SRS Diode, and Sun Nuclear EDGE) and a synthetic diamond detector (PTW 60019 micro Diamond [µD] detector) were used for OAR, PDD, and OF measurements. OFs were also acquired with a PTW 31022 PinPoint ionization chamber. Beam scanning was performed using a 3D water tank at depths of 7, 50, 100, 200, and 250 mm with a source-to-surface distance of 450 mm. OFs were measured at the depth of maximum dose (d_max  = 7 mm) with the SAD at 450 mm. Gafchromic EBT3 film was used to validate OF and profile measurements and as a reference detector for estimating correction factors for active detector OFs. Deviations in field size, penumbra, and PDDs across the different detectors were quantified.

RESULTS

Relative OFs (ROFs) for the diodes were within 1.4% for all collimators except for 5 and 7.5 mm, for which SRS Diode measurements were higher by 1.6% and 2.6% versus Diode E. The µD ROFs were within 1.4% of the diode measurements. PinPoint ROFs were lower by >10% for the 4-mm and 5-mm collimators versus the Diode E and µD. Corrections to OFs using EBT3 film as a reference were within 1.2% for all diodes and the µD detector for collimators 10 mm and greater and within 2.0%, 2.8%, and 1.1% for the 7.5-, 5-, and 4-mm collimators, respectively. The maximum difference in full width at half maximum (FWHM) between the Diode E and the other active detectors was for the 25-mm collimator and was 0.09 mm (µD), 0.16 mm (SRS Diode), and 0.65 mm (EDGE). Differences seen in PDDs beyond the depth of d_max were <1% across the three diodes and the µD. FWHM and penumbra measurements made using EBT3 film were within 1.34% and 3.26%, respectively, of the processed profile data entered into the treatment planning system.

CONCLUSIONS

Minimal differences were seen in OAR and PDD measurements acquired with the diodes and the µD. ROFs measured with the three diodes were within 2.6% and within 1.4% versus the µD. Gafchromic Film measurements provided independent verification of the OAR and OF measurements. Estimated corrections to OFs using film as a reference were <1.6% for the Diode E, EDGE, and µD detector.

Experimental investigation of TRS-483 reference dosimetry correction factors for Leksell Gamma Knife® Icon™ beams

Purpose

Radiosurgery using the Leksell Gamma Knife® (LGK) Icon™ is an established technique used for treating intracranial lesions. The largest beam field size the LGK Icon can produce is a 16 mm diameter sphere. Despite this, reference dosimetry on the LGK Icon is typically performed using ionization chambers calibrated in 10 × 10 cm² fields. Furthermore, plastic phantoms are widely used instead of liquid water phantoms. In an effort to resolve these issues, the International Atomic Energy Agency (IAEA) in collaboration with American Association of Physicists in medicine (AAPM) recently published Technical Report Series No. 483 (TRS‐483) as a Code of Practice for small‐field dosimetry. TRS‐483 includes small‐field correction factors, kQmsr,Q0fmsr,fref, intended to account for the differences between setups when using small‐field modalities such as the LGK Icon, and conventional setups. Since the publication of TRS‐483, at least three new sets of values of kQmsr,Q0fmsr,fref for the LGK Icon have been published. The purpose of this study was to experimentally investigate the published values of kQmsr,Q0fmsr,fref for commonly used phantom and ionization chamber (IC) models for the LGK Icon.

Methods

Dose‐rates from two LGK units were determined using acrylonitrile butadiene styrene (ABS) and Certified Medical Grade Solid Water® (SW) phantoms, and PTW 31010 and PTW 31016 ICs. Correction factors were applied, and the resulting dose‐rates compared. Relative validity of the correction factors was investigated by taking the ratios of dose‐rate correction factor products. Additionally, dose‐rates from the individual sectors were determined in order to calculate the beam attenuation caused by the ABS phantom adapter.

Results and Conclusions

It was seen that the dose‐rate is underestimated by at least 1% when using the ABS phantom, which was attributed to fluence perturbation caused by the IC and phantom adapter. Published correction factors kQmsr,Q0fmsr,fref account for these effects to varying degree and should be used. The SW phantom is unlikely to underestimate the dose‐rate by more than 1%, and applying kQmsr,Q0fmsr,fref could not be shown to be necessary. Out of the two phantom models, the ABS phantom is not recommended for use in LGK reference dosimetry. The use of newly published values of kQmsr,Q0fmsr,fref should be considered.

Dosimetric evaluation of the Leksell GammaPlan™ Convolution dose calculation algorithm

The dosimetric accuracy of the Leksell GammaPlan Convolution calculation algorithm was evaluated through comparison with corresponding Monte Carlo (MC) dosimetric results. MC simulations were based on generated sector phase space files for the 4 mm, 8 mm and 16 mm collimator sizes, using a previous comprehensive Gamma Knife Perfexion^™ source model and validated using film dosimetry. Test cases were designed for the evaluation of the Convolution algorithm involving irradiation of homogeneous and inhomogeneous phantom geometries mimicking clinical cases, with radiation fields created using one sector (single sector), all sectors with the same (single shot) or different (composite shot) collimator sizes. Dose calculations using the Convolution algorithm were found to be in excellent agreement (gamma pass rate greater than 98%, applying 1%/1 mm local dose difference and distance agreement criteria), with corresponding MC calculations, indicating the accuracy of the Convolution algorithm in homogeneous and heterogeneous model geometries. While of minor clinical importance, large deviations were observed for the voxels laying inside air media. The calculated beam on times using the Convolution algorithm were found to increase (up to 7%) relative to the TMR 10 algorithm currently used in clinical practice, especially in a test case mimicking a brain metastasis close to the skull, in excellent agreement with corresponding MC calculations.

Output correction factors for small static fields in megavoltage photon beams for seven ionization chambers in two orientations - perpendicular and parallel

Purpose

The goal of the present work was to provide a large set of detector‐specific output correction factors for seven small volume ionization chambers on two linear accelerators in four megavoltage photon beams utilizing perpendicular and parallel orientation of ionization chambers in the beam for nominal field sizes ranging from 0.5 cm² × 0.5 cm² to 10 cm² × 10 cm². The present study is the second part of an extensive research conducted by our group.

Methods

Output correction factors kQclin,Qreffclin,fref were experimentally determined on two linacs, Elekta Versa HD and Varian TrueBeam for 6 and 10 MV beams with and without flattening filter for nine square fields ranging from 0.5 cm² × 0.5 cm² to 10 cm² × 10 cm², for seven mini and micro ionization chambers, IBA CC04, IBA Razor, PTW 31016 3D PinPoint, PTW 31021 3D Semiflex, PTW 31022 3D PinPoint, PTW 31023 PinPoint, and SI Exradin A16. An Exradin W1 plastic scintillator and EBT3 radiochromic films were used as the reference detectors.

Results

For all ionization chambers, values of output correction factors kQclin,Qreffclin,fref were lower for parallel orientation compared to those obtained in the perpendicular orientation. Five ionization chambers from our study set, IBA Razor, PTW 31016 3D PinPoint, PTW 31022 3D PinPoint, PTW 31023 PinPoint, and SI Exradin A16, fulfill the requirement recommended in the TRS‐483 Code of Practice, that is, 0.95<kQclin,Qreffclin,fref<1.05, down to the field size 0.8 cm² × 0.8 cm², when they are positioned in parallel orientation; two of the ionization chambers, IBA Razor and PTW 31023 PinPoint, satisfy this condition down to the field size of 0.5 cm² × 0.5 cm².

Conclusions

The present paper provides experimental results of detector‐specific output correction factors for seven small volume ionization chambers. Output correction factors were determined in 6 and 10 MV photon beams with and without flattening filter down to the square field size of 0.5 cm² × 0.5 cm² for two orientations of ionization chambers — perpendicular and parallel. Our main finding is that output correction factors are smaller if they are determined in a parallel orientation compared to those obtained in a perpendicular orientation for all ionization chambers regardless of the photon beam energy, filtration, or linear accelerator being used. Based on our findings, we recommend using ionization chambers in parallel orientation, to minimize corrections in the experimental determination of field output factors. Latter holds even for field sizes below 1.0 cm² × 1.0 cm², whenever necessary corrections remain within 5%, which was the case for several ionization chambers from our set.

TRS‐483 recommended perpendicular orientation of ionization chambers for the determination of field output factors. The present study presents results for both perpendicular and parallel orientation of ionization chambers. When validated by other researchers, the present results for parallel orientation can be considered as a complementary dataset to those given in TRS‐483.

Quality assurance for Gamma Knife Perfexion using the Exradin W1 plastic scintillation detector and Lucy phantom

The goal of this work is to validate the use of the Exradin W1 plastic scintillation detector (PSD) to measure profiles and output factors from Gamma Knife Perfexion collimators in a Lucy phantom. The Exradin W1 PSD has a small-volume, near-water-equivalent, energy-independent sensitive element. Output measurements were performed for all 3 collimators (4 mm, 8 mm, and 16 mm) of the Gamma Knife Perfexion system, and these measurements were compared to measurements made with an A16 ion chamber and an EBT3 film and to the nominal values. We showed that a configuration in which the focus or 'shot' moves while the detector remains fixed is essentially equivalent to a configuration in which the focus is fixed while the detector moves. A Lucy phantom containing a PSD was moved in small steps to acquire profiles in all three dimensions. EBT3 film was inserted in the Lucy phantom and exposed to a single shot for each collimator. The relative values for output factors measured with the PSD were 1.000, 0.892, and 0.795, for the 16 mm, 8 mm, and 4 mm collimators, respectively. The values measured with EBT3 film were 1.000, 0.881, and 0.793, and the values measured with the A16 ion chamber were 1.000, 0.883, and 0.727. The nominal output factors for the Gamma Knife Perfexion are 1.000, 0.900, and 0.814, respectively. There was excellent agreement between all profiles measured with the PSD and EBT3 as well as with the treatment planning system data provided by the vendor. In light of our results, the Exradin W1 PSD is well suited for beam quality assurance of a Gamma Knife Perfexion irradiator.

A novel method for the determination of field output factors and output correction factors for small static fields for six diodes and a microdiamond detector in megavoltage photon beams

Purpose

The goal of this work is to provide a large and consistent set of data for detector‐specific output correction factors, kQclin,Qreffclin,fref, for small static fields for seven solid‐state detectors and to determine field output factors, ΩQclin,Qreffclin,fref, using EBT3 radiochromic films and W1 plastic scintillator as reference detectors on two different linear accelerators and four megavoltage photon beams. Consistent measurement conditions and recommendations given in the International Code of Practice TRS‐483 for small‐field dosimetry were followed throughout the study.

Methods

ΩQclin,Qreffclin,fref were determined on two linacs, Elekta Versa HD and Varian TrueBeam, for 6 and 10 MV beams with and without flattening filter and for nine fields ranging from 0.5 × 0.5 cm² to 10 × 10 cm². Signal readings obtained with EBT3 radiochromic films and W1 plastic scintillator were fitted by an analytical function. Volume averaging correction factors, determined from two‐dimensional (2D) dose matrices obtained with EBT3 films and fitted to bivariate Gaussian function, were used to correct measured signals. kQclin,Qreffclin,fref were determined empirically for six diodes, IBA SFD, IBA Razor, PTW 60008 P, PTW 60012 E, PTW 60018 SRS, and SN EDGE, and a PTW 60019 microDiamond detector.

Results

Field output factors and detector‐specific kQclin,Qreffclin,fref are presented in the form of analytical functions as well as in the form of discrete values. It is found that in general, for a given linac, small‐field output factors need to be determined for every combination of beam energy and filtration (WFF or FFF) and field size as the differences between them can be statistically significant (P < 0.05). For different beam energies, the present data for kQclin,Qreffclin,fref are found to differ significantly (P < 0.05) from the corresponding data published in TRS‐483 mostly for the smallest fields (<1.5 cm). For the PTW microDiamond detector, statistically significant differences (P < 0.05) between kQclin,Qreffclin,fref values were found for all investigated beams on an Elekta Versa HD linac for field sizes 0.5 × 0.5 cm² and 0.8 × 0.8 cm². Significant differences in kQclin,Qreffclin,fref between beams of a given energy but with and without flattening filters are found for measurements made in small fields (<1.5 cm) at a given linac. Differences in kQclin,Qreffclin,fref are also found when measurements are made at different linacs using the same beam energy filtration combination; for the PTW microDiamond detector, these differences were found to be around 6% and were considered as significant.

Conclusions

Selection of two reference detectors, EBT3 films and W1 plastic scintillator, and use of an analytical function, is a novel approach for the determination of ΩQclin,Qreffclin,fref for small static fields in megavoltage photon beams. Large set of kQclin,Qreffclin,fref data for seven solid‐state detectors and four beam energies determined on two linacs by a single group of researchers can be considered a valuable supplement to the literature and the TRS‐483 dataset.

Electrometer offset current due to scattered radiation

Relative dose measurements with small ionization chambers in combination with an electrometer placed in the treatment room ("internal electrometer") show a large dependence on the polarity used. While this was observed previously for percent depth dose curves (PDDs), the effect has not been understood or preventable. To investigate the polarity dependence of internal electrometers used in conjunction with a small-volume ionization chamber, we placed an internal electrometer at a distance of 1 m from the isocenter and exposed it to different amounts of scattered radiation by varying the field size. We identified irradiation of the electrometer to cause a current of approximately -1 pA, regardless of the sign of the biasing voltage. For low-sensitivity detectors, such a current noticeably distorts relative dose measurements. To demonstrate how the current systematically changes PDDs, we collected measurements with nine ionization chambers of different volumes. As the chamber volume decreased, signal ratios at 20 and 10 cm depth (M20/M10) became smaller for positive bias voltage and larger for negative bias voltage. At the size of the iba CC04 (40 mm³) the difference of M20/M10 was around 1% and for the smallest studied chamber, the iba CC003 chamber (3 mm³), around 7% for a 10 × 10 cm² field. When the electrometer was moved further from the source or shielded, the additional current decreased. Consequently, PDDs at both polarities were brought into alignment at depth even for the 3 mm³ ionization chamber. The apparent polarity effect on PDDs and lateral beam profiles was reduced considerably by shielding the electrometer. Due to normalization the effect on output values was low. When measurements with a low-sensitivity probe are carried out in conjunction with an internal electrometer, we recommend careful monitoring of the particular setup by testing both polarities, and if deemed necessary, we suggest shielding the electrometer.

Determination of kQmsr,Q0fmsr,fref factors for ion chambers used in the calibration of Leksell Gamma Knife Perfexion model using EGSnrc and PENELOPE Monte Carlo codes

PURPOSE

To calculate the kQmsr,Q0fmsr,fref factors for nine common ionization chamber types following the small fields dosimetry formalism for the calibration of the Leksell Gamma Knife^® (LGK) Perfexion^TM using Monte Carlo simulation. This study also provides the first independent comparison of EGSnrc and PENELOPE for the calculation of kQmsr,Q0fmsr,fref correction factors and proposes a practical method to predict these factors based on chamber type, chamber orientation and phantom electron density.

METHODS

The ionization chambers are modeled using the EGSnrc and PENELOPE Monte Carlo codes based on the blueprints provided by the manufacturers. The chambers are placed in a half-sphere water phantom and five spherical phantoms made of liquid water, solid water, ABS, polystyrene, and PMMA, respectively. Dose averaged over the air cavity of the chambers and a small water volume are calculated using EGSnrc and PENELOPE Monte Carlo codes for both conventional and machine specific reference (msr) setups. Using the calculated dose ratio, the kQmsr,Q0fmsr,fref factor is determined for all phantom materials and two possible orientations of chamber. The calculated kQmsr,Q0fmsr,fref factors are compared to a previous Monte Carlo study. A relationship between the kQmsr,Q0fmsr,fref factor and the electron density of the phantom material is derived to predict the kQmsr,Q0fmsr,fref factor for any phantom material type. Applying the calculated kQmsr,Q0fmsr,fref factors to the measured dose rate of a recent round robin study improves consistency of reference dosimetry of the Leksell Gamma Knife^® (LGK) Perfexion^TM .

RESULTS

Agreement within uncertainty is observed between kQmsr,Q0fmsr,fref values determined in this study and the previous PEGASOS/PENELOPE study in a liquid water phantom. The difference between kQmsr,Q0fmsr,fref values in parallel and perpendicular detector orientations is most significant for the PTW 31010 (1.8%) chamber. The percentage root-mean-square (%RMS) deviation between EGSnrc and PENELOPE calculated kQmsr,Q0fmsr,fref values for Exradin-A1SL, A14 and A14SL chambers studies in this work was found to be 0.4%. The kQmsr,Q0fmsr,fref values increase linearly with electron density of the phantom material for all chamber types mainly due to the linear dependency of photon energy fluence ratios on electron density. The average percentage difference between the calculated and predicted kQmsr,Q0fmsr,fref values using two methods is found to be 0.15% and 0.16%. Previously measured dose rates corrected with the kQmsr,Q0fmsr,fref values determined in this work leads to absorbed dose values consistent to within 0.8%.

CONCLUSIONS

The calculated kQmsr,Q0fmsr,fref values in this work will enable users to apply the appropriate correction for their own specific phantom material only knowing the electron density of the phantom material.