Determination of output correction factors in magnetic fields using two methods for two detectors at the central axis

Objective.Commissioning an MR-linac treatment planning system requires output correction factors,kB→,Qclin,Qmsrfclin,fmsr, for detectors to accurately measure the linac's output at various field sizes. In this study,kB→,Qclin,Qmsrfclin,fmsrwas determined at the central axis using two methods: one that combines the corrections for the influence of the magnetic field and the small field in a single factor (kB→,Qclin,Qmsrfclin,fmsr), and a second that isolates the magnetic field's influence, allowing the use of output correction factors without a magnetic field,kQclin,Qmsrfclin,fmsr, from literature for determiningkB→,Qclin,Qmsrfclin,fmsr.Approach.To determinekB→,Qclin,Qmsrfclin,fmsrand examine its behaviour across different photon energies and magnetic flux densitiesBin small fields, measurements with an ionization chamber (0.07 cm3sensitive volume) and a solid-state detector were carried out at an experimental facility for both approaches. Changes in absorbed dose to water with field size were determined via Monte Carlo simulations. To evaluate clinical applicability, additional measurements were conducted on a 1.5 T MR-linac.Main results.Both methods determined comparablekB→,Qclin,Qmsrfclin,fmsrresults. For field sizes >3 × 3 cm2,Branging from -1.5 to 1.5 T and photon energies of 6 and 8 MV, no change ofkQclin,Qmsrfclin,fmsras a function of the magnetic field was observed. Comparison with measurement results from the 1.5 T MR-linac confirm this. For ⩽3 × 3 cm2,kB→,Qclin,Qmsrfclin,fmsrdepends on photon energy andB. For 1.5 T and 6 MV,BreduceskQclin,Qmsrfclin,fmsrup to 3% for the ionization chamber and up to 7% for the solid-state detector.Significance.kB→,Qclin,Qmsrfclin,fmsrwere successfully determined for two detectors, enabling their use at a 1.5 T MR-linac. For field sizes of >3 × 3 cm2,kB→,Qclin,Qmsrfclin,fmsris one for most detectors suitable for small field dosimetry for all available perpendicular MR-linac systems, as confirmed in the literature. For these field sizes and detectors, the correction factor accounting for the dosimeter response change in the reference field due to the magnetic field,kB→,Qmsrfmsr, can be used for cross-calibration. Therefore, future research may only focus on small field sizes.

On the correction factors for small field dosimetry in 1.5T MR-linacs

Objective. Clinical dosimetry in the presence of a 1.5 T magnetic field is challenging, let alone in case small fields are involved. The scope of this study is to determine a set of relevant correction factors for a variety of MR-compatible detectors with emphasis on small fields. Two dosimetry formalisms adopted from the literature are considered.Approach. Six small-cavity ionization chambers (from three manufacturers), four active solid-state detectors and a thermoluminescence dosimeter microcube were modeled in the EGSnrc Monte Carlo code. Phase space files for field sizes down to 1 × 1 cm2of the Unity 1.5 T/7 MV MR-linac (Elekta, UK) were used as source models. Simulations were performed to calculate thekQB,QfB,f(also known askB,Q),kQmsrB,fmsrandkQclin,QmsrB,fclin,fmsrrelevant to two different dosimetry formalisms. Two detector orientations with respect to the magnetic field were considered. Moreover, the effect of the ionization chamber's stem length (a construction parameter) on the correction factor was investigated. Simulations were also carried out to determine whether correction factors obtained in water can be applied in dosimetry procedures involving water-equivalent solid phantoms.Main results. Under thekQB,QfB,f-based formalism, the required corrections to ionization chamber responses did not exceed 1.5% even for the smallest field size considered. A much wider range ofkQB,QfB,fvalues was obtained for the active solid-state detectors included in the simulations. This is the first study to reportkQclin,QmsrB,fclin,fmsrvalues for ionization chambers. The impact of the stem on correction factors is not significant for lengths ⩾0.75 cm. Correction factors determined in water are also valid in dosimetry protocols employing solid phantoms.Significance. This work substantially expands the range of available detectors that can be used in small field dosimetry, enabling more options for commissioning, beam modeling and quality assurance procedures in 1.5 T MR-Linacs. However, more studies are needed to establish a complete and reliable dataset.

Small field measurements using electronic portal imaging device

Purpose/Objective. Small-field measurement poses challenges. Although many high-resolution detectors are commercially available, the EPID for small-field dosimetry remains underexplored. This study aimed to evaluate the performance of EPID for small-field measurements and to derive tailored correction factors for precise small-field dosimetry verification.Material/Methods. Six high-resolution radiation detectors, including W2 and W1 plastic scintillators, Edge-detector, microSilicon, microDiamond and EPID were utilized. The output factors, depth doses and profiles, were measured for various beam energies (6 MV-FF, 6 MV-FFF, 10 MV-FF, and 10 MV-FFF) and field sizes (10 × 10 cm2, 5 × 5 cm2, 4 × 4 cm2, 3 × 3 cm2, 2 × 2 cm2, 1 × 1 cm2, 0.5 × 0.5 cm2) using a Varian Truebeam linear accelerator. During measurements, acrylic plates of appropriate depth were placed on the EPID, while a 3D water tank was used with five-point detectors. EPID measured data were compared with W2 plastic scintillator and measurements from other high-resolution detectors. The analysis included percentage deviations in output factors, differences in percentage for PDD and for the profiles, FWHM, maximum difference in the flat region, penumbra, and 1D gamma were analyzed. The output factor and depth dose ratios were fitted using exponential functions and fractional polynomial fitting in STATA 16.2, with W2 scintillator as reference, and corresponding formulae were obtained. The established correction factors were validated using two Truebeam machines.Results. When comparing EPID and W2-PSD across all field-sizes and energies, the deviation for output factors ranged from 1% to 15%. Depth doses, the percentage difference beyond dmax ranged from 1% to 19%. For profiles, maximum of 4% was observed in the 100%-80% region. The correction factor formulae were validated with two independent EPIDs and closely matched within 3%.Conclusion. EPID can effectively serve as small-field dosimetry verification tool with appropriate correction factors.

Cerenkov free micro-dosimetry in small-field radiation therapy technique

Objective. Optical fiber-based scintillating dosimetry is a recent promising technique owing to the miniature size dosimeter and quality measurement in modern radiation therapy treatment. Despite several advantages, the major issue of using scintillating dosimeters is the Cerenkov effect and predominantly requires extra measurement corrections. Therefore, this work highlighted a novel micro-dosimetry technique to ensure Cerenkov-free measurement in radiation therapy treatment protocol by investigating several dosimetric characteristics.Approach.A micro-dosimetry technique was proposed with the performance evaluation of a novel infrared inorganic scintillator detector (IR-ISD). The detector essentially consists of a micro-scintillating head based on IR-emitting micro-clusters with a sensitive volume of 1.5 × 10-6mm3. The proposed system was evaluated under the 6 MV LINAC beam used in patient treatment. Overall measurements were performed using IBATMwater tank phantoms by following TRS-398 protocol for radiotherapy. Cerenkov measurements were performed for different small fields from 0.5 × 0.5 cm2to 10 × 10 cm2under LINAC. In addition, several dosimetric parameters such as percentage depth dose (PDD), high lateral resolution beam profiling, dose linearity, dose rate linearity, repeatability, reproducibility, and field output factor were investigated to realize the performance of the novel detector.Main results. This study highlighted a complete removal of the Cerenkov effect using a point-like miniature detector, especially for small field radiation therapy treatment. Measurements demonstrated that IR-ISD has acceptable behavior with dose rate variability (maximum standard deviation ∼0.18%) for the dose rate of 20-1000 cGy s-1. An entire linear response (R2= 1) was obtained for the dose delivered within the range of 4-1000 cGy, using a selected field size of 1 × 1 cm2. Perfect repeatability (max 0.06% variation from average) with day-to-day reproducibility (0.10% average variation) was observed. PDD profiles obtained in the water tank present almost identical behavior to the reference dosimeter with a build-up maximum depth dose at 1.5 cm. The small field of 0.5 × 0.5 cm2profiles have been characterized with a high lateral resolution of 100µm.Significance. Unlike recent plastic scintillation detector systems, the proposed micro-dosimetry system in this study requires no Cerenkov corrections and showed efficient performance for several dosimetric parameters. Therefore, it is expected that considering the detector correction factors, the IR-ISD system can be a suitable dose measurement tool, such as in small-field dose measurements, high and low gradient dose verification, and, by extension, in microbeam radiation and FLASH radiation therapy.

Grid/lattice therapy: consideration of small field dosimetry

Small-field dosimetry used in special procedures such as gamma knife, Cyberknife, Tomotherapy, IMRT, and VMAT has been in evolution after several radiation incidences with very significant (70%) errors due to poor understanding of the dosimetry. IAEA-TRS-483 and AAPM-TG-155 have provided comprehensive information on small-fields dosimetry in terms of code of practice and relative dosimetry. Data for various detectors and conditions have been elaborated. It turns out that with a suitable detectors dose measurement accuracy can be reasonably (±3%) achieved for 6 MV beams for fields >1×1 cm2. For grid therapy, even though the treatment is performed with small fields created by either customized blocks, multileaf collimator (MLC), or specialized devices, it is multiple small fields that creates combined treatment. Hence understanding the dosimetry in collection of holes of small field is a separate challenge that needs to be addressed. It is more critical to understand the scattering conditions from multiple holes that form the treatment grid fields. Scattering changes the beam energy (softer) and hence dosimetry protocol needs to be properly examined for having suitable dosimetric parameters. In lieu of beam parameter unavailability in physical grid devices, MLC-based forward and inverse planning is an alternative path for bulky tumours. Selection of detectors in small field measurement is critical and it is more critical in mixed beams created by scattering condition. Ramification of small field concept used in grid therapy along with major consideration of scattering condition is explored. Even though this review article is focussed mainly for dosimetry for low-energy megavoltage photon beam (6 MV) but similar procedures could be adopted for high energy beams. To eliminate small field issues, lattice therapy with the help of MLC is a preferrable choice.

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.

Small field output correction factors at 18 MV

BACKGROUND

The response of various detectors in the radiotherapy energy range has been investigated, especially for 6 and 10 MV energies for small fields, and is summarized in TRS-483. However, data for accelerator energies above 10 MV are sparse or unavailable for many detectors, especially for the energy of 18 MV. Small variations in field output factors for the commissioning of a treatment planning system can have a high impact on calculation of dose distributions.

PURPOSE

Many studies describe an energy dependence of the response for a large number of detectors. We wanted to close the gap for the 18 MV energy regime and determined field output correction factors for different detectors at 18 MV.

METHODS

An ELEKTA Versa HD accelerator at 18 MV was used together with a PTW MP3 water phantom at an SSD of 90 cm. The following detectors were examined: PTW Semiflex 31021, PinPoint 3D 31022, diode 60012, diode 60008 and microDiamond 60019, Sun Nuclear EDGE detector, IBA PFD, SFD, Razor Chamber, Razor Nano Chamber and Razor Diode, Standard Imaging Scintillator Exradin W2 1x3, W2 1x1 and Gafchromic EBT3 film. The dose response was determined at a depth of 10 cm for square fields between 0.5 and 10 cm side length. As reference data a composure of radiochromic film data for small fields (s ≤ 3 $s\le 3$  cm) and data of all compatible chambers for larger fields (s ≥ 3 $s\ge 3$  cm) was used. The effective field sizes of small fields were determined from profiles obtained on radiochromic film. The obtained field output correction factors obey the rules of the TRS-483 protocol.

RESULTS

The W2 1x1 scintillator and the Razor Chamber showed the smallest deviations from the reference curve. The shielded diodes (diode 60008, EDGE detector) showed the highest over-response at small fields, followed by PFD, microDiamond and the unshielded diodes (diode 60012, SFD). The ionization chambers exhibited the well-known volume effect, that is, strong under-response at small fields of up to 9% for the PinPoint 3D, 7% for the Razor Chamber and up to 30% for the Semiflex detector for the smallest studied field size. The small chambers showed a polarity effect in axial orientation, especially the Razor Nano Chamber. Corrections at 18 MV are generally larger than those provided by TRS-483, continuing the trend of increasing corrections between 6 and 10 MV also at a higher accelerator energy. Only the PinPoint 3D Chamber showed a slightly smaller correction.

CONCLUSIONS

Field output correction factors were determined for square field sizes between 0.5 and 10 cm at 18 MV. Most detectors needed a larger correction than at 6 and 10 MV. Thus, the use of correction factors will improve beam data for 18 MV.

Dosimetry in 1.5 T MR-Linacs: Monte Carlo determination of magnetic field correction factors and investigation of the air gap effect

BACKGROUND

In Magnetic Resonance-Linac (MR-Linac) dosimetry formalisms, a new correction factor, k_B,Q , has been introduced to account for corresponding changes to detector readings under the beam quality, Q, and the presence of magnetic field, B.

PURPOSE

This study aims to develop and implement a Monte Carlo (MC)-based framework for the determination of k_B,Q correction factors for a series of ionization chambers utilized for dosimetry protocols and dosimetric quality assurance checks in clinical 1.5 T MR-Linacs. Their dependencies on irradiation setup conditions are also investigated. Moreover, to evaluate the suitability of solid phantoms for dosimetry checks and end-to-end tests, changes to the detector readings due to the presence of small asymmetrical air gaps around the detector's tip are quantified.

METHODS

Phase space files for three irradiation fields of the ELEKTA Unity 1.5 T/7 MV flattening-filter-free MR-Linac were provided by the manufacturer and used as source models throughout this study. Twelve ionization chambers (three farmer-type and nine small-cavity detectors, from three manufacturers) were modeled (including their dead volume) using the EGSnrc MC code package. k_B,Q values were calculated for the 10 × 10 cm² irradiation field and for four cardinal orientations of the detectors' axes with respect to the 1.5 T magnetic field. Potential dependencies of k_B,Q values with respect to field size, depth, and phantom material were investigated by performing additional simulations. Changes to the detectors' readings due to the presence of small asymmetrical air gaps (0.1 up to 1 mm) around the chambers' sensitive volume in an RW3 solid phantom were quantified for three small-cavity chambers and two orientations.

RESULTS

For both parallel (to the magnetic field) orientations, k_B,Q values were found close to unity. The maximum correction needed was 1.1%. For each detector studied, the k_B,Q values calculated for the two parallel orientations agreed within uncertainties. Larger corrections (up to 5%) were calculated when the detectors were oriented perpendicularly to the magnetic field. Results were compared with corresponding ones found in the literature, wherever available. No considerable dependence of k_B,Q with respect to field size (down to 3 × 3 cm² ), depth, or phantom material was noticed, for the detectors investigated. As compared to the perpendicular one, in the parallel to the magnetic field orientation, the air gap effect is minimized but is still considerable even for the smallest air gap considered (0.1 mm).

CONCLUSION

For the 10 × 10 cm² field, magnetic field correction factors for 12 ionization chambers and four orientations were determined. For each detector, the k_B,Q value may be also applied for dosimetry procedures under different irradiation parameters provided that the orientation is taken into account. Moreover, if solid phantoms are used, even the smallest asymmetrical air gap may still bias small-cavity chamber response. This work substantially expands the availability and applicability of k_B,Q correction factors that are detector- and orientation-specific, enabling more options in MR-Linac dosimetry checks, end-to-end tests, and quality assurance protocols.

A novel QA phantom based on scintillating fiber ribbons with implementation of 2D dose tomography for small-field radiotherapy

PURPOSE

To develop a novel instrument for real-time quality assurance (QA) procedures in radiotherapy. The system implements a scintillation-based phantom and associated signal acquisition and processing modules and aims to monitor two-dimensional (2D) dose distributions of small fields.

MATERIALS AND METHODS

For the proposed phantom, we have designed and realized a prototype implementing six high-resolution tissue-equivalent scintillating fiber ribbons stacked with in-plane 30° rotated orientations from each other. Each ribbon output is coupled to a silicon photodiode linear array (with an element pitch of 400 μm) to detect scintillating signal, which represents the projected irradiation profile perpendicular to the ribbon's orientation. For the system providing six acquired projected dose profiles at different orientations, we have developed a two-step signal processing method to perform 2D dose reconstruction. The first step is to determine irradiation field geometry parameters using a tomographic geometry approach, and the second one is to perform specific penumbra estimation. The QA system prototype has been tested on a Novalis TrueBeam STX with a 6-MV photon beam for small elliptic fields defined by 5- and 10-mm cone collimators and for 10 × 10- and 20 × 10-mm² rectangular fields defined by the micro-multileaf collimator. Gamma index analysis using EBT3 films as reference has been carried out with tight 2%-dose-difference (DD)/700-μm-distance-to-agreement (DTA) as well as 1%-DD/1-mm-DTA criteria for evaluating the system performances. The testing also includes an evaluation of the proposed two-step field reconstruction method in comparison with two conventional methods: filtered back projection (FBP) and simultaneous iterative reconstruction technique (SIRT).

RESULTS

The reconstructed 2D dose distributions have gamma index pass rates higher than 95% for all the tested configurations as compared with EBT3 film measurements with both 2%-DD/700-μm-DTA and 1%-DD/1-mm criteria. 2D global gamma analysis shows that the two-step and FBP radiation field reconstruction methods systematically outperform the SIRT approach. Moreover, higher gamma index success rates are obtained with the two-step method than with FBP in the case of the fields defined with the stereotactic cones.

CONCLUSIONS

The proposed small-field QA system makes a use of six water-equivalent scintillating detectors (fiber ribbons) to acquire dose distribution. The developed two-step signal processing method performs tomographic 2D dose reconstruction. A system prototype has been built and tested using hospital facilities with small rectangular and elliptic fields. Testing results show 2D reconstructed dose distributions with high accuracy and resolution. Such a system could potentially be an alternative approach to film dosimetry for small-field QA, which is still widely used as reference in clinical practice.

Mitigating the uncertainty in small field dosimetry by leveraging machine learning strategies

Small field dosimetry is significantly different from the dosimetry of broad beams due to loss of electron side scatter equilibrium, source occlusion, and effects related to the choice of detector. However, use of small fields is increasing with the increase in indications for intensity-modulated radiation therapy and stereotactic body radiation therapy, and thus the need for accurate dosimetry is ever more important. Here we propose to leverage machine learning (ML) strategies to reduce the uncertainties and increase the accuracy in determining small field output factors (OFs). Linac OFs from a Varian TrueBeam STx were calculated either by the treatment planning system (TPS) or measured with a W1 scintillator detector at various multi-leaf collimator (MLC) positions, jaw positions, and with and without contribution from leaf-end transmission. The fields were defined by the MLCs with the jaws at various positions. Field sizes between 5 and 100 mm were evaluated. Separate ML regression models were generated based on the TPS calculated or the measured datasets. Accurate predictions of small field OFs at different field sizes (FSs) were achieved independent of jaw and MLC position. A mean and maximum % relative error of 0.38 ± 0.39% and 3.62%, respectively, for the best-performing models based on the measured datasets were found. The prediction accuracy was independent of contribution from leaf-end transmission. Several ML models for predicting small field OFs were generated, validated, and tested. Incorporating these models into the dose calculation workflow could greatly increase the accuracy and robustness of dose calculations for any radiotherapy delivery technique that relies heavily on small fields.

Beam commissioning of the first clinical biology-guided radiotherapy system

This study reports the beam commissioning results for the first clinical RefleXion Linac.

Methods: The X1 produces a 6 MV photon beam and the maximum clinical field size is 40 × 2 cm² at source‐to‐axis distance of 85 cm. Treatment fields are collimated by a binary multileaf collimator (MLC) system with 64 leaves with width of 0.625 cm and y‐jaw pairs to provide either a 1 or 2 cm opening. The mechanical alignment of the radiation source, the y‐jaw, and MLC were checked with film and ion chambers. The beam parameters were characterized using a diode detector in a compact water tank. In‐air lateral profiles and in‐water percentage depth dose (PDD) were measured for beam modeling of the treatment planning system (TPS). The lateral profiles, PDDs, and output factors were acquired for field sizes from 1.25 × 1 to 40 × 2 cm² field to verify the beam modeling. The rotational output variation and synchronicity were tested to check the gantry angle, couch motion, and gantry rotation.

Results: The source misalignments were 0.049 mm in y‐direction, 0.66% out‐of‐focus in x‐direction. The divergence of the beam axis was 0.36 mm with a y‐jaw twist of 0.03°. Clinical off‐axis treatment fields shared a common center in y‐direction were within 0.03 mm. The MLC misalignment and twist were 0.57 mm and 0.15°. For all measured fields ranging from the size from 1.25 × 1 to 40 × 2 cm², the mean difference between measured and TPS modeled PDD at 10 cm depth was −0.3%. The mean transverse profile difference in the field core was −0.3% ± 1.1%. The full‐width half maximum (FWHM) modeling was within 0.5 mm. The measured output factors agreed with TPS within 0.8%.

Conclusions: This study summarizes our specific experience commissioning the first novel RefleXion linac, which may assist future users of this technology when implementing it into their own clinics.

Dose area product primary standards established by graphite calorimetry at the LNE-LNHB for small radiation fields in radiotherapy

PURPOSE

To present primary standards establishment in terms of Dose Area Product (DAP) for small field sizes.

METHODS

A large section graphite calorimeter and two plane-parallel ionization chambers were designed and built in-house. These chambers were calibrated in a 6MV FFF beam at the maximum dose rate of 1400 UM/min for fields defined by specifically designed circular collimators of 5, 7.5, 10, 13 and 15 mm diameter and jaws of 5, 7, 10, 13 and 15 mm side length on a Varian TrueBeam linac.

RESULTS

The two chambers show the same behaviour regardless of field shape and size. From 5 to 15 mm, calibration coefficients slightly increase with the field size with a magnitude of 1.8% and 1.1% respectively for the two chambers, and are independent of the field shape. This tendency was confirmed by Monte Carlo calculations. The average associated uncertainty of the calibration coefficients is around 0.6% at k=1.

CONCLUSIONS

For the first time, primary standards in terms of DAP were established by graphite calorimetry for an extended range of small field sizes. These promising results open the door for an alternative approach in small fields dosimetry.

A probabilistic approach for determining Monte Carlo beam source parameters: II. Impact of beam modeling uncertainties on dosimetric functions and treatment plans

Objective. The Monte Carlo method is recognized as a valid approach for the evaluation of dosimetric functions for clinical use. This procedure requires the accurate modeling of the considered linear accelerator. In Part I, we propose a new method to extract the probability density function of the beam model physical parameters. The aim of this work is to evaluate the impact of beam modeling uncertainties on Monte Carlo evaluated dosimetric functions and treatment plans in the context of small fields. Approach. Simulations of output factors, output correction factors, dose profiles, percent-depth doses and treatment plans are performed using the CyberKnife M6 model developed in Part I. The optimized pair of electron beam energy and spot size, and eight additional pairs of beam parameters representing a 95% confidence region are used to propagate the uncertainties associated to the source parameters to the dosimetric functions. Main results. For output factors, the impact of beam modeling uncertainties increases with the reduction of the field size and confidence interval half widths reach 1.8% for the 5 mm collimator. The impact on output correction factors cancels in part, leading to a maximum confidence interval half width of 0.44%. The impact is less significant for percent-depth doses in comparison to dose profiles. For these types of measurement, in absolute terms and in comparison to the reference dose, confidence interval half widths less than or equal to 1.4% are observed. For simulated treatment plans, the impact is more significant for the treatment delivered with a smaller field size with confidence interval half widths reaching 2.5% and 1.4% for the 5 and 20 mm collimators, respectively. Significance. Results confirm that AAPM TG-157's tolerances cannot apply to the field sizes studied. This study provides an insight on the reachable dose calculation accuracy in a clinical setup.

Investigation of field output factors using IAEA-AAPM TRS-483 code of practice recommendations and Monte Carlo simulation for 6 MV photon beams

Introduction:

This study aims to experimentally determine field output factors using the methodologies suggested by the IAEA-AAPM TRS-483 for small field dosimetry and compare with the calculation from Monte Carlo (MC) simulation.

Methods:

The IBA-CC01, Sun Nuclear EDGE and IBA-SFD detectors were employed to determine the uncorrected and the corrected field output factors for 6 MV photon beams. Measurements were performed at 100 cm source to axis distance, 10 cm depth in water, and the field sizes ranged from 1 × 1 to 10 × 10 cm2. The use of field output correction factors proposed by the TRS-483 was utilised to determine field output factors. The measured field output factors were compared to that calculated using the egs_chamber user code.

Results:

The decrease in the percentage standard deviation of the measured three detectors was observed after applying the field output correction factors. Measured field output factors using CC01 and EDGE detectors agreed with MC values within 3% for field sizes down to 1 × 1 cm2, except the SFD detector.

Conclusions:

The corrected field output factors agree with the calculation from MC, except the SFD detector. CC01 and EDGE are suitable for determining field output factors, while the SFD may need more implementation of the intermediate field method.

Small-field measurement and Monte Carlo model validation of a novel image-guided radiotherapy system

PURPOSE

The RefleXion™ X1 is a novel radiotherapy system that is designed for image-guided radiotherapy, and eventually, biology-guided radiotherapy (BgRT). BgRT is a treatment paradigm that tracks tumor motion using real-time positron emission signals. This study reports the small-field measurement results and the validation of a Monte Carlo (MC) model of the first clinical RefleXion unit.

METHODS

The RefleXion linear accelerator (linac) produces a 6 MV flattening filter free (FFF) photon beam and consists of a binary multileaf collimator (MLC) system with 64 leaves and two pairs of y-jaws. The maximum clinical field size achievable is 400 × 20 mm² . The y-jaws provide either a 10 or 20 mm opening at source-to-axis distance (SAD) of 850 mm. The width of each MLC leaf at SAD is 6.25 mm. Percentage depth doses (PDDs) and relative beam profiles were acquired using an Edge diode detector in a water tank for field sizes from 12.5 × 10 to 100 × 20 mm² . Beam profiles were also measured using films. Output factors of fields ranging from 6.25 × 10 to 100 × 20 mm² were measured using W2 scintillator detector, Edge detector, and films. Output correction factors k of the Edge detector for RefleXion were calculated. An MC model of the linac including pre-MLC beam sources and detailed structures of MLC and lower y-jaws was validated against the measurements. Simulation codes BEAMnrc and GATE were utilized.

RESULTS

The diode measured PDD at 10 cm depth (PDD10) increases from 53.6% to 56.9% as the field opens from 12.5 × 10 to 100 × 20 mm² . The W2-measured output factor increases from 0.706 to 1 as the field opens from 6.25 × 10 to 100 × 20 mm² (reference field size). The output factors acquired by diode and film differ from the W2 results by 1.65% (std = 1.49%) and 2.09% (std = 1.41%) on average, respectively. The profile penumbra and full-width half-maximum (FWHM) measured by diode agree well with the film results with a deviation of 0.60 mm and 0.73% on average, respectively. The averaged beam profile consistency calculated between the diode- and film-measured profiles among different depths is within 1.72%. By taking the W2 measurements as the ground truth, the output correction factors k for Edge detector ranging from 0.958 to 1 were reported. For the MC model validation, the simulated PDD10 agreed within 0.6% to the diode measurement. The MC-simulated output factor differed from the W2 results by 2.3% on average (std = 3.7%), while the MC simulated beam penumbra differed from the diode results by 0.67 mm on average (std = 0.42 mm). The MC FWHM agreed with the diode results to within 1.40% on average. The averaged beam profile consistency calculated between the diode and MC profiles among different depths is less than 1.29%.

CONCLUSIONS

This study represents the first small-field dosimetry of a clinical RefleXion system. A complete and accurate MC model of the RefleXion linac has been validated.

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.

Determination of field output correction factors of radiophotoluminescence glass dosimeter and CC01 ionization chamber and validation against IAEA-AAPM TRS-483 code of practice

PURPOSE

To determine the field output correction factors of the radiophotoluminescence glass dosimeter (RPLGD) in parallel and perpendicular orientations with reference to CC01, the ionization chamber.

METHODS

The dose to a small water volume and the sensitive volume of the RPLGD and the IBA-CC01 were determined for 6-MV, 100-cm SAD, 10-cm depth using egs_chamber user-code. The RPLGD in perpendicular and parallel orientations to the beam axis were studied. The field output correction factors of each detector for 0.5 × 0.5 to 10 × 10 cm² field sizes were determined. These field output correction factors were validated by comparing field output factors against data determined from IAEA-AAPM TRS-483 code of practice.

RESULTS

The field output correction factors of all detectors were within 5% for field sizes down to 0.8 × 0.8 cm². For 0.5 × 0.5 cm², the field output correction factors of CC01, RPLGD in perpendicular and parallel orientations differed from unity by 14%, 19%, and 5%, respectively. The percentage difference between field output factors determined using RPLGD and CC01 data, corrected using the field output correction factors determined in this work and measurements with CC01 data corrected using TRS-483, was less than 3% for all field sizes, except for the smallest field size of RPLGD in perpendicular orientation and the CC01.

CONCLUSIONS

The field output correction factors of RPLGD and CC01 are reported. The validation proves that RPLGD in parallel orientation combined with the field output correction factors is the most suitable for determining the field output factors for the smallest field used in this study.

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.

Fast Monte Carlo codes for occupational dosimetry in interventional radiology

PURPOSE

Interventional radiology techniques cause radiation exposure both to patient and personnel. The radiation dose to the operator is usually measured with dosimeters located at specific points above or below the lead aprons. The aim of this study is to develop and validate two fast Monte Carlo (MC) codes for radiation transport in order to improve the assessment of individual doses in interventional radiology. The proposed methodology reduces the number of required dosemeters and provides immediate dose results.

METHODS

Two fast MC simulation codes, PENELOPE/penEasyIR and MCGPU-IR, have been developed. Both codes have been validated by comparing fast MC calculations with the multipurpose PENELOPE MC code and with measurements during a realistic interventional procedure.

RESULTS

The new codes were tested with a computation time of about 120 s to estimate operator doses while a standard simulation needs several days to obtain similar uncertainties. When compared with the standard calculation in simple set-ups, MCGPU-IR tends to underestimate doses (up to 5%), while PENELOPE/penEasyIR overestimates them (up to 18%). When comparing both fast MC codes with experimental values in realistic set-ups, differences are within 25%. These differences are within accepted uncertainties in individual monitoring.

CONCLUSION

The study highlights the fact that computational dosimetry based on the use of fast MC codes can provide good estimates of the personal dose equivalent and overcome some of the limitations of occupational monitoring in interventional radiology. Notably, MCGPU-IR calculates both organ doses and effective dose, providing a better estimate of radiation risk.

Dose-area product ratio in external small-beam radiotherapy: beam shape, size and energy dependencies in clinical photon beams

In small-field radiotherapy (RT), a significant challenge is to define the amount of radiation dose absorbed in the patient where the quality of the beam has to be measured with high accuracy. The properties of a proposed new beam quality specifier, namely the dose-area-product ratio at 20 and 10 cm depths in water or DAPR_20,10, were studied to yield more information on its feasibility over the conventional quality specifier tissue-phantom ratio or TPR_20,10. The DAPR_20,10may be measured with a large-area ionization chamber (LAC) instead of small volume chambers or semi-conductors where detector, beam and water phantom positioning and beam perturbations introduce uncertainties. The effects of beam shape, size and energy on the DAPR_20,10were studied and it was shown that the DAPR_20,10increases with increasing beam energy similarly to TPR_20,10but in contrast exhibits a small beam size and shape dependence. The beam profile outside the beam limiting devices has been shown to have a large contribution to the DAPR_20,10. There is potential in large area chambers to be used in DAPR measurement and its use in dosimetry of small-beam RT for beam quality measurements.

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%.

Monte Carlo verification of output correction factors for a TrueBeam STx®

The recent publication of the new code of practice IAEA/AAPM TRS-483 introduces output correction factors to correct detector response changes in relative dosimetry of small photon beams. In TRS-483, average correction factors are reported for several detectors in high-energy photon beams at 6 and 10 MV with and without flattening filter. These correction factors were determined by Monte Carlo simulation or experimental measurements using several linacs of different brands and vendors. The goal of this work was to validate the output correction factors reported in TRS-483 for 6 MV photon beams of a TrueBeam STx® linac. The validation was performed using Monte Carlo simulations of four radiation detectors employed in the dosimetry of small photon beams and whose output correction factors were determined using a different radiation source than TrueBeam STx®. The results show that Monte Carlo calculated output correction factors, and those reported in the code of practice TRS-483 fully agree within ∼1%. The use of generic correction factors for a TrueBeam STx® and the detectors studied in this work is suitable for small field dosimetry static beams within the uncertainties of Monte Carlo calculations and output correction factors reported in TRS-483.

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.

Measurement of the small field output factors for 10 MV photon beam using IAEA TRS-483 dosimetry protocol and implementation in Eclipse TPS commissioning

Dosimetry of small fields (SF) is vital for the success of highly conformal techniques. IAEA along with AAPM recently published a code of practice TRS-483 for SF dosimetry. The scope of this paper is to investigate the performance of three different detectors with 10 MV with-flatting-filter (WFF) beam using TRS-483 for SF dosimetry and subsequent commissioning of the Eclipse treatment planning system (TPS version-13.6) for SF data. SF dosimetry data (beam-quality TPR 20,10(10), cross-calibration, beam-profile, and field-output-factor (F.O.F)) measurements were performed for PTW31006-pinpoint, IBA-CC01 and IBA-EFD-3G diode detectors in nominal field size (F.S) range 0.5 × 0.5cm2 to 10 × 10 cm2 with water and solid water medium using Varian Truebeam linac. However, Eclipse-TPS commissioning data was acquired using IBA-EFD-3G diode, and absolute dose calibration was performed with FC-65G detector. The dosimetric performance of the Eclipse-TPS was validated using TLD-LiF chips, IBA-PFD, and IBA-EFD-3G diodes. Dosimetric performance of the PTW31006-pinpoint, IBA-CC01, and IBA-EFD-3G detectors was successfully tested for SF dosimetry. The F.O.Fs were generated and found in close agreement for all F.S except 0.5 × 0.5cm2. It is also found that TPR20,10(10) value can be derived within 0.5% accuracy from a non-reference field using Palmans equation. Cross-calibration can be performed in F.S 6 × 6 cm2 with a maximum variation of 0.5% with respect to 10 × 10cm2. During profile measurement, the full-width half-maxima (FWHM) of F.S 0.5 × 0.5cm2 was found maximum deviated from the geometric F.S. In addition, Eclipse-TPS was commissioned along with some limitations: F.O.F below F.S 1 × 1cm2 was ignored by TPS, PDD and profiles were dropped from configuration below F.S 2 × 2 cm2, and F.O.F which does not satisfy the condition 0.7 < A/B < 1.4 (A and B are FWHM in cross-line and in-line direction) have higher uncertainty than specified in TRS-483. Validation tests for Eclipse-TPS generated plans were also performed. The measured dose was in close agreement (3%) with TPS calculated dose up to F.S 1.5 × 1.5cm2.

Feasibility study of using Stereotactic Field Diode for field output factors measurement and evaluating three new detectors for small field relative dosimetry of 6 and 10 MV photon beams

This study assesses the feasibility of using stereotactic field diode (SFD) as an alternate to gaf chromic films for field output factor (FF) measurement and further evaluating three new detectors for small field dosimetry. Varian 21EX linear accelerator was used to generate 6 and 10 MV beams of nominal square fields ranging from 0.5 × 0.5 cm² to 10 × 10 cm². One passive (EBT3 films) and five active detectors including IBA RAZOR diode(RD), SFD, RAZOR nanochamber (RNC), pinpoint chamber (PTW31023), and semiflex chamber (PTW31010) were employed. FFs were measured using films and SFD while beam profiles and percentage depth dose (PDD) distribution were acquired with active detectors. Polarity (k_pol) and recombination (k_s) effects of ion chambers were determined and corrected for output ratio measurement. Correction factors (CF) of RD, RNC, and PTW31023 in axial and radial orientation were also measured. Stereotactic field diode measured FFs have shown good agreement with films (with difference of <1%). RD and RNC measured beam profiles were within 3% deviation from the SFD values. Variation in k_pol with field size for RNC and PTW31023 was up to 4% and 0.4% (for fields ≥ 1 × 1 cm²), respectively, while variation in k_s of PTW31023 was <0.2 %. The maximum values of CF have been calculated to be 5.2%, 2.0%, 13.6%, and 25.5% for RD, RNC, PTW31023‐axial, and PTW31023‐radial respectively. This study concludes that SFD with appropriate CFs as given in TRS 483 may be used for measuring FFs as an alternate to EBT3 films. Whereas RD and RNC may be used for beam profile and PDD measurement in small fields. Considering the limit of usability of 2%, RNC may be used without CF for FF measurement in the smallfields investigated in this study.

Small field output correction factors of the microSilicon detector and a deeper understanding of their origin by quantifying perturbation factors

PURPOSE

The aim of this study is the experimental and Monte Carlo-based determination of small field correction factors for the unshielded silicon detector microSilicon for a standard linear accelerator as well as the Cyberknife System. In addition, a detailed Monte Carlo analysis has been performed by modifying the detector models stepwise to study the influences of the detector's components.

METHODS

Small field output correction factors have been determined for the new unshielded silicon diode detector, microSilicon (type 60023, PTW Freiburg, Germany) as well as for the predecessors Diode E (type 60017, PTW Freiburg, Germany) and Diode SRS (type 60018, PTW Freiburg, Germany) for a Varian TrueBeam linear accelerator at 6 MV and a Cyberknife system. For the experimental determination, an Exradin W1 scintillation detector (Standard Imaging, Middleton, USA) has been used as reference. The Monte Carlo simulations have been performed with EGSnrc and phase space files from IAEA as well as detector models according to manufacturer blueprints. To investigate the influence of the detector's components, the detector models have been modified stepwise.

RESULTS

The correction factors for the smallest field size investigated at the TrueBeam linear accelerator (equivalent dosimetric square field side length Sclin  = 6.3 mm) are 0.983 and 0.939 for the microSilicon and Diode E, respectively. At the Cyberknife system, the correction factors of the microSilicon are 0.967 at the smallest 5-mm collimator compared to 0.928 for the Diode SRS. Monte Carlo simulations show comparable results from the measurements and literature.

CONCLUSION

The microSilicon (type 60023) detector requires less correction than its predecessors, Diode E (type 60017) and Diode SRS (type 60018). The detector housing has been demonstrated to cause the largest perturbation, mainly due to the enhanced density of the epoxy encapsulation surrounding the silicon chip. This density has been rendered more water equivalent in case of the microSilicon detector to minimize the associated perturbation. The sensitive volume itself has been shown not to cause observable field size-dependent perturbation except for the volume-averaging effect, where the slightly larger diameter of the sensitive volume of the microSilicon (1.5 mm) is still small at the smallest field size investigated with corrections <2%. The new microSilicon fulfils the 5% correction limit recommended by the TRS 483 for output factor measurements at all conditions investigated in this work.

Monte carlo model and output factors of elekta infinity™ 6 and 10 MV photon beam

Aim

This study aimed to commission the Elekta Infinity™ working in 6 and 10 MV photon beam installed in Concord International Hospital, Singapore, and compare the OFs between MC simulation and measurement using PTW semiflex and microDiamond detector for small field sizes.

Material and Methods

There are two main steps in this study: modelling of Linac 6 and 10 MV photon beam and analysis of the output factors for field size 2 × 2-10 × 10 cm². The EGSnrc/BEAMnrc-DOSXYZnrc code was used to model and characterize the Linac and to calculate the dose distributions in a water phantom. The dose distribution and OFs were compared to the measurement data in the same condition.

Results

The commissioning process was only conducted for a 10 × 10 cm² field size. The PDD obtained from MC simulation showed a good agreement with the measurement. The local dose difference of PDDs was less than 2% for 6 and 10 MV. The initial electron energy was 5.2 and 9.4 MeV for 6 and 10 MV photon beam, respectively. This Linac model can be used for dose calculation in other situations and different field sizes because this Linac has been commissioned and validated using Monte Carlo simulation. The 10 MV Linac produces higher electron contamination than that of 6 MV.

Conclusions

The Linac model in this study was acceptable. The most important result in this work comes from OFs resulted from MC calculation. This value was more significant than the OFs from measurement using semiflex and microDiamond for all beam energy and field sizes because of the CPE phenomenon.

A Novel GPU-based Fast Monte Carlo Photon Dose Calculating Method for Accurate Radiotherapy Treatment Planning

BACKGROUND

An accurate and fast radiation dose calculations method is the main part of treatment planning for successful radiation therapy.

OBJECTIVE

This work aimed to create a novel GPU-based fast Monte Carlo Photon Dose Code (MCPDC) as a fast and accurate tool in dose calculation for radiotherapy treatment planning.

MATERIALS AND METHODS

In this analytical study, MCDPC was written to implement photon MC simulation for energies 0.01 to 20 MeV and run on an NVIDIA GTX970. The code was validated using DOSXYZnrc results and experimental measurements, performed by a Mapcheck dosimeter. Using the innovative definition of photon and electron interactions, mean calculation time for the MCPDC was 5.4 sec for 5e7 source particle history, significantly less than that of DOSXYZnrc which was 400 min.

RESULTS

Considering the simulations in the anthropomorphic phantom with bone and lung inhomogeneity, more than 96.1% of all significant voxels passed the gamma criteria of 3%-3 mm. Compared to the experimental dosimetry results, 97.6% or more of all significant voxels passed the acceptable clinical gamma index of 3%-3 mm.

CONCLUSION

Very fast calculation speed and high accuracy in dose calculation may allow the MCPDC to be used in radiotherapy as a central component of a treatment plan verification system and also as the dose calculation engine for MC-based planning. MCPDC is currently being developed for electron dose calculation module and graphic user interface. In addition, future work on the applicability of the improved version of the MCPDC in transit dosimetry of megavoltage CT is in process.

Characterization of a microSilicon diode detector for small-field photon beam dosimetry

This study characterized a new unshielded diode detector, the microSilicon (model 60023), for small-field photon beam dosimetry by evaluating the photon beams generated by a TrueBeam STx and a CyberKnife. Temperature dependence was evaluated by irradiating photons and increasing the water temperature from 11.5 to 31.3°C. For Diode E, microSilicon, microDiamond and EDGE detectors, dose linearity, dose rate dependence, energy dependence, percent-depth-dose (PDD), beam profiles and detector output factor (OFdet) were evaluated. The OFdet of the microSilicon detector was compared to the field output factors of the other detectors. The microSilicon exhibited small temperature dependence within 0.4%, although the Diode E showed a linear variation with a ratio of 0.26%/°C. The Diode E and EDGE detectors showed positive correlations between the detector reading and dose rate, whereas the microSilicon showed a stable response within 0.11%. The Diode E and microSilicon demonstrated negative correlations with the beam energy. The OFdet of microSilicon was the smallest among all the detectors. The maximum differences between the OFdet of microSilicon and the field output factors of microDiamond were 2.3 and 1.6% for 5 × 5 mm2 TrueBeam and 5 mm φ CyberKnife beams, respectively. The PDD data exhibited small variations in the dose fall-off region. The microSilicon and microDiamond detectors yielded similar penumbra widths, whereas the other detectors showed steeper penumbra profiles. The microSilicon demonstrated favorable characteristics including small temperature and dose rate dependence as well as the small spatial resolution and output factors suitable for small field dosimetry.

The impact of corrected field output factors based on IAEA/AAPM code of practice on small-field dosimetry to the calculated monitor unit in eclipse™ treatment planning system

The objective of this study was to investigate the effect of field output factors (FOFs) according to the current protocol for small-field dosimetry in conjunction to treatment planning system (TPS) commissioning. The calculated monitor unit (MU) for intensity-modulated radiation therapy (IMRT) and volumetric modulated arc therapy (VMAT) plans in Eclipse™ TPS were observed. Micro ion chamber (0.01 CC) (CC01), photon field diode (shielded diode) (PFD), and electron field diode (unshielded diode) (EFD) were used to measure percentage depth doses, beam profiles, and FOFs from 1 × 1 cm2 to 10 × 10 cm2 field sizes of 6 MV photon beams. CC01 illustrated the highest percentage depth doses at 10 cm depth while EFD exhibited the lowest with the difference of 1.6% at 1 × 1 cm2 . CC01 also produced slightly broader penumbra, the difference with other detectors was within 1 mm. For uncorrected FOF of three detectors, the maximum percent standard deviation (%SD) was 5.4% at 1 × 1 cm2 field size. When the correction factors were applied, this value dropped to 2.7%. For the calculated MU in symmetric field sizes, beam commissioning group from uncorrected FOF demonstrated maximum %SD of 6.0% at 1 × 1 cm2 field size. This value decreased to 2.2% when the corrected FOF was integrated. For the calculated MU in IMRT-SRS plans, the impact of corrected FOF reduced the maximum %SD from 6.0% to 2.5% in planning target volume (PTV) less than 0.5 cm3 . Beam commissioning using corrected FOF also decreased %SD for VMAT-SRS plans, although it was less pronounced in comparison to other treatment planning techniques, since the %SD remained less than 2%. The use of FOFs based on IAEA/AAPM TRS 483 has been proven in this research to reduce the discrepancy of calculated MU among three beam commissioning datasets in Eclipse™ TPS. The dose measurement of both symmetric field and clinical cases comparing to the calculation illustrated the dependence of the types of detector commissioning and the algorithm of the treatment planning for small field size.

SciFi detector and associated method for real-time determination of profile and output factor for small fields in stereotactic radiotherapy

PURPOSE

For determining small-field profile and output factor during stereotactic radiotherapy quality assurance (QA) procedures, we propose a novel system based on the scintillating fiber (SciFi) detector with output image acquisition and processing to allow real-time monitoring of profile and output factor.

MATERIALS AND METHODS

The employed detector is a SciFi detector made of tissue-equivalent scintillating plastic fibers arranged in 6-layer fiber ribbons with a fiber pitch of 275 μm in each layer. The scintillating signal at the detector output is acquired by a sCMOS (scientific complementary metal-oxide-semiconductor) camera and represents the projected field profile along the fibers axis. An iterative reconstruction method of the field from its projected profile based on a priori knowledge of some features of the radiation field defined by the stereotactic cones is suggested. The detector with implemented data processing has been tested in clinical conditions, for determining beam profiles and output factors, using cone collimators of different sizes from 4 to 15 mm diameter. The detector under test was placed at 1.4 cm depth and 98.6 cm source to surface distance (SSD) in a water-equivalent phantom and irradiated by a 6 MV photon beam.

RESULTS

The reconstructed field profiles obtained from the detector are coherent with data from EBT3 radiochromic films, with differences within ±0.32 mm for both the FWHM and the penumbra region. For real-time determination of the field output factor, the measured data are also in good agreement with data independently determined by the French Institute for Radiological Protection and Nuclear Safety (IRSN) based on radiochromic films and thermoluminescent 1 × 1 mm2 micro-cubes dosimeters (TLD). The differences are within ±1.6% for all the tested cone sizes.

CONCLUSIONS

We propose and have tested a SciFi plastic scintillating detector with an optimized signal processing method to characterize small fields defined by cone collimators. It allows the determination of key field parameters such as full width at half maximum (FWHM) and field output factors. The results are consistent with those independently measured using TLD and radiochromic films. As the SciFi detector does not require a correction factor, it is in line with the International Atomic Energy Agency (IAEA) and the American Association of Physicists in Medicine (AAPM) TRS-483 recommendations, and can be suitable for online QA of small radiation fields used in photon beam radiotherapy, and is compatible with MRI-LINAC.

Evaluation of dosimetric parameters of small fields of 6 MV flattening filter free photon beam measured using various detectors against Monte Carlo simulation

Purpose:

This study aims to evaluate dosimetric parameters like percentage depth dose, dosimetric field size, depth of maximum dose surface dose, penumbra and output factors measured using IBA CC01 pinpoint chamber, IBA stereotactic field diode (SFD), PTW microDiamond against Monte Carlo (MC) simulation for 6 MV flattening filter-free small fields.

Materials and Methods:

The linear accelerator used in the study was a Varian TrueBeam® STx. All field sizes were defined by jaws. The required shift to effective point of measurement was given for CC01, SFD and microdiamond for depth dose measurements. The output factor of a given field size was taken as the ratio of meter readings normalised to 10 × 10 cm2 reference field size without applying any correction to account for changes in detector response. MC simulation was performed using PRIMO (PENELOPE-based program). The phase space files for MC simulation were adopted from the MyVarian Website.

Results and Discussion:

Variations were seen between the detectors and MC, especially for fields smaller than 2 × 2 cm2 where the lateral charge particle equilibrium was not satisfied. Diamond detector was seen as most suitable for all measurements above 1 × 1 cm2. SFD was seen very close to MC results except for under-response in output factor measurements. CC01 was observed to be suitable for field sizes above 2 × 2 cm2. Volume averaging effect for penumbra measurements in CC01 was observed. No detector was found suitable for surface dose measurement as surface ionisation was different from surface dose due to the effect of perturbation of fluence. Some discrepancies in measurements and MC values were observed which may suggest effects of source occlusion, shift in focal point or mismatch between real accelerator geometry and simulation geometry.

Conclusion:

For output factor measurement, TRS483 suggested correction factor needs to be applied to account for the difference in detector response. CC01 can be used for field sizes above 2 × 2 cm2 and microdiamond detector is suitable for above 1 × 1 cm2. Below these field sizes, perturbation corrections and volume averaging corrections need to be applied.

An analytical formalism for the assessment of dose uncertainties due to positioning uncertainties

PURPOSE

To present an analytical formalism for the in depth assessment of uncertainties of field output factors in small fields related to detector positioning based on dose profile measurements. Additionally, a procedure for the propagation of these uncertainties was developed.

METHODS

Based on the assumption that one dimensional and two dimensional second-order polynomial functions can be fitted to dose profiles of small photon beams, equations for the calculation of the expectation value, the variance, and the standard deviation were developed. The following fitting procedures of the dose profiles were considered: A one-dimensional case (1D), a quasi two-dimensional case (2Dq) based on independently measured line profiles and a full 2D case (2Df) which also considers cross-correlations in a two-dimensional dose distribution. A rectangular and a Gaussian probability density function (PDF) characterizing the probability of possible positions of the detector relative to the maximum dose were used. Uncertainty components such as the finite resolution of the scanning water phantom, the reproducibility of the determination of the position of the maximum dose, and the reproducibility of the collimator system were investigated. This formalism was tested in a 0.5 x 0.5 cm² photon field where dose profiles were measured using a radiochromic film, a synthetic diamond detector, and an unshielded diode detector. Additionally, the dose distribution measured with the radiochromic film was convoluted with a convolution kernel mimicking the active volume of the unshielded diode.

RESULTS

Analytic expressions for the calculation of uncertainties on field output factors were found for the 1D, the 2Dq, and the 2Df case. The uncertainty of the field output factor related to the relative position of the detector to the maximum dose increased quadratically with increasing limits of possible detector positions. Analysis of the radiochromic film showed that the 2Dq case gave a more conservative assessment of the uncertainty compared to the 2Df case with a difference of < 0.1%. The 2Dq case applied to the film measurements agreed well with the same approach as was applied to the unshielded diode. The investigated uncertainty components propagated to an uncertainty of the field output factors of 0.5% and 0.4% for the synthetic diamond and the unshielded diode, respectively. Additionally, the expectation value was lower than the maximum dose. The difference was 0.4% and 0.3% for the synthetic diamond and the unshielded diode, respectively.

CONCLUSIONS

The assessment of uncertainties of field output factors related to detector positioning is feasible using the proposed formalism. The 2Dq case is applicable when using online detectors. Accurate positioning in small fields is essential for accurate dosimetry as its related uncertainty increases quadratically. The observed drop of the expectation value needs to be considered in small field dosimetry.

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 cm2  × 0.5 cm2 to 10 cm2  × 10 cm2 . The present study is the second part of an extensive research conducted by our group.

METHODS

Output correction factors k Q clin , Q ref f clin , f ref 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 cm2  × 0.5 cm2 to 10 cm2  × 10 cm2 , 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 k Q clin , Q ref f clin , f ref 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 < k Q clin , Q ref f clin , f ref < 1.05 , down to the field size 0.8 cm2  × 0.8 cm2 , 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 cm2  × 0.5 cm2 .

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 cm2  × 0.5 cm2 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 cm2  × 1.0 cm2 , 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.

Small-field dosimetry of TrueBeamTM flattened and flattening filter-free beams: A multi-institutional analysis

PURPOSE

Detector-dependent interinstitutional variations of the beam data may lead to uncertainties of the delivered dose to patients. Here we evaluated the inter-unit variability of the flattened and flattening filter-free (FFF) beam data of multiple TrueBeam (Varian Medical Systems) linear accelerators focusing on the small-field dosimetry.

METHODS

The beam data of 6- and 10-MV photon beams with and without flattening filter measured for modeling of an iPLAN treatment planning system (BrainLAB) were collected from 12 institutions - ten HD120 Multileaf Collimator (MLC) and two Millennium120 MLC. Percent-depth dose (PDD), off-center ratio (OCR), and detector output factors (OFdet ) measured with different detectors were evaluated. To investigate the detector-associated effects, we evaluated the inter-unit variations of the OFdet before and after having applied the output correction factors provided by the International Atomic Energy Agency (IAEA) Technical Reports Series no. 483.

RESULTS

PDD measured with a field size of 5 × 5 mm2 showed that the data measured using an ionization chamber had variations exceeding 1% from the median values. The maximum difference from median value was 2.87% for 10 MV photon beam. The maximum variations of the penumbra width for OCR with 10 × 10 mm2 field size were 0.97 mm. The OFdet showed large variations exceeding 15% for a field size of 5 × 5 mm2 . When the output correction factors were applied to the OFdet , the variations were greatly reduced. The relative difference of almost all field output factors were within ± 5% from the median field output factors.

CONCLUSION

In this study, the inter-unit variability of small-field dosimetry was evaluated for TrueBeam linear accelerators. The variations were large at a field size of 5 × 5 mm2 , and most occurred in a detector-dependent manner.

An EGS Monte Carlo model for Varian TrueBEAM treatment units: Commissioning and experimental validation of source parameters

In this work we have created and commissioned a Monte Carlo model of 6FFF Varian TrueBeam linear accelerator using BEAMnrc. For this purpose we have experimentally measured the focal spot size and shape of three Varian TrueBeam treatment units in 6FFF modality with a slit collimator and several depth dose and lateral beam profiles in a water phantom. The Monte Carlo model of a 6FFF TrueBeam machine was implemented with a primary electron source commissioned as a 2D Gaussian with Full Width Half Maximum selected by comparison of simulated and measured narrow beam profiles. The energy of the primary electron beam was optimized through a simultaneous fit to the measured beam depth dose profiles. Special attention was paid to evaluation of uncertainties of the selected Monte Carlo source parameters. These uncertainties were calculated by analysing the sensitivity of the commissioning process to changes in both primary beam size and energy. Both experimental and Monte Carlo commissioned focus size values were compared and found to be in excellent agreement. The commissioned Monte Carlo model reproduces within 1% accuracy the dose distributions of radiation field size from 3 cm × 3 cm to 15 cm × 15 cm.

Radiological tissue equivalence of deformable silicone-based chemical radiation dosimeters (FlexyDos3D)

FlexyDos3D, a silicone‐based chemical radiation dosimeter, has great potential to serve as a three‐dimensional (3D) deformable dosimetric tool to verify complex dose distributions delivered by modern radiotherapy techniques. To facilitate its clinical application, its radiological tissue needs to be clarified. In this study we investigated its tissue‐equivalence in comparison with water and Solid Water (RMI457). We found that its effective and mean atomic numbers were 40% and 20% higher and the total interaction probabilities for kV x‐ray photons were larger than those of water respectively. To assess the influence of its over‐response to kV photons, its HU value was measured by kV computed tomography (CT) and was found higher than all the soft‐tissue substitutes. When applied for dose calculation without correction, this effect led to an 8% overestimation in electron density via HU‐value mapping and 0.65% underestimation in target dose. Furthermore, depth dose curves (PDDs) and off‐axis ratios (profiles) at various beam conditions as well as the dose distribution of a full‐arc VMAT plan in FlexyDos3D and reference materials were simulated by Monte Carlo, where the results showed great agreement. As indicated, FlexyDos3D exhibits excellent radiological water‐equivalence for clinical MV x‐ray dosimetry, while its nonwater‐equivalent effect for low energy x‐ray dosimetry requires necessary correction. The key findings of this study provide pertinent reference for further FlexyDos3D characterization research.

Monte Carlo determination of k Q and k Qmsr values for the exradin A26 ionisation chamber for the Varian TrueBeam

We have calculated conversion factors, k _Q for the A26 micro ionisation chamber along with machine specific reference beam quality factors, k _Qmsr, for a number of field sizes and beam qualities for the Varian TrueBeam accelerator. The A12 ionisation chamber was simulated alongside the A26, so as to validate against known literature values. Both ionisation chambers were modelled from manufacturer data sheets and schematics. The egs_chamber Monte Carlo user code was used to simulate each absorbed dose relevant to the beam quality conversion factors k _Q and k _Qmsr. Tabulated spectra for beam energies of 4 through 25 MV were used in the k _Q calculations for both investigated chambers. Varian TrueBeam phase space files for 6 MV flattened as well as 6 and 10 MV unflattened beams were used in the simulations of the A26 chamber in field sizes from 2  ×  2 cm square to 20  ×  20 cm square in order to determine k _Qmsr values. The PDD(10)x  values of the tabulated spectra were found to be within variation between studies, with an average deviance of 0.4% from one prior study. The simulated A12 k _Q values matched the accepted literature values with an average variation of  <0.1%. The A26 k _Q values match the manufacturer provided values to within 0.5%. For all investigated field sizes the k _Qmsr values are within 0.006 of unity. There is no published data for this chamber for a direct comparison, but there is similarity between these results and results from other chambers regularly used in similar circumstances. Furthermore, the agreement of the simulated k _Q values to knowns, and the agreement of the PDD(10)x  factors would suggest the correctness and accuracy of the study.

Detectors assessment for stereotactic radiosurgery with cones

The purpose of this work is to assess eight detectors performance for output factor (OF), percent depth dose (PDD), and beam profiles in a 6‐MV Clinac stereotactic radiosurgery mode for cone irradiation using Monte Carlo simulation as reference. Cones with diameters comprised between 30 and 4 mm have been studied. The evaluated detectors were ionization chambers: pinpoint and pinpoint 3D, diodes: SRS, P and E, Edge, MicroDiamond and EBT3 radiochromic films. The results showed that pinpoints underestimate OF up to −2.3% for cone diameters ≥10 mm and down to −12% for smaller cones. Both nonshielded (SRS and E) and shielded diodes (P and Edge) overestimate the OF respectively up to 3.3% and 5.2% for cone diameters ≥10 mm and in both cases more than 7% for smaller cones. MicroDiamond slightly overestimates the OF, 3.7% for all the cones and EBT3 film is the closest to Monte Carlo with maximum difference of ±1% whatever the cone size is. For the profiles and the PDD, particularly for the small cones, the size of the detector predominates. All diodes and EBT3 agree with the simulation within ±0.2 mm for beam profiles determination. For PDD curve all the active detectors response agree with simulation up to 1% for all the cones. EBT3 is the more accurate detector for beam profiles and OF determinations of stereotactic cones but it is restrictive to use. Due to respectively inappropriate size of the sensitive volume and composition, pinpoints and diodes do not seem appropriate without OF corrective factors below 10 mm diameter cone. MicroDiamond appears to be the best detector for OF determination regardless all cones. For off‐axis measurements, the size of the detector predominates and for PDD all detectors give promising results.