Monte Carlo-based diode design for correction-less small field dosimetry

The PTW microSilicon diode: Performance in small 6 and 15 MV photon fields and utility of density compensation

PURPOSE

We have experimentally and computationally characterized the PTW microSilicon 60023-type diode's performance in 6 and 15 MV photon fields ≥5 × 5 mm² projected to isocenter. We tested the detector on- and off-axis at 5 and 15 cm depths in water, and investigated whether its response could be improved by including within it a thin airgap.

METHODS

Experimentally, detector readings were taken in fields generated by a Varian TrueBeam linac and compared with doses-to-water measured using Gafchromic film and ionization chambers. An unmodified 60023-type diode was tested along with detectors modified to include 0.6, 0.8, and 1.0 mm thick airgaps. Computationally, doses absorbed by water and detectors' sensitive volumes were calculated using the EGSnrc/BEAMnrc Monte Carlo radiation transport code. Detector response was characterized using k Q c l i n , 4 cm f c l i n , 4 cm , a factor that corrects for differences in the ratio of dose-to-water to detector reading between small fields and the reference condition, in this study 5 cm deep on-axis in a 4 × 4 cm² field.

RESULTS

The greatest errors in measurements of small field doses made using uncorrected readings from the unmodified 60023-type detector were over-responses of 2.6% ± 0.5% and 5.3% ± 2.0% determined computationally and experimentally, relative to the reading-per-dose in the reference field. Corresponding largest errors for the earlier 60017-type detector were 11.9% ± 0.6% and 11.7% ± 1.4% over-responses. Adding even the thinnest, 0.6 mm, airgap to the 60023-type detector over-corrected it, leading to under-responses of up to 4.8% ± 0.6% and 5.0% ± 1.8% determined computationally and experimentally. Further, Monte Carlo calculations indicate that a detector with a 0.3 mm airgap would read correctly to within 1.3% on-axis. The ratio of doses at 15 and 5 cm depths in water in a 6 MV 4 × 4 cm² field was measured more accurately using the unmodified 60023-type detector than using the 60017-type detector, and was within 0.3% of the ratio measured using an ion chamber. The 60023-type diode's sensitivity also varied negligibly as dose-rate was reduced from 13 to 4 Gy min^-1 by decreasing the linac pulse repetition frequency, whereas the sensitivity of the 60017-type detector fell by 1.5%.

CONCLUSIONS

The 60023-type detector performed well in small fields across a wide range of beam energies, field sizes, depths, and off-axis positions. Its response can potentially be further improved by adding a thin, 0.3 mm, airgap.

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.

Technical Note: Small field dose correction factors for radiochromic film in lung phantoms

PURPOSE

Radiochromic film has been established as a detector that can be used without the need for perturbation correction factors for small field dosimetry in water. However, perturbation factors in low density media such as lung have yet to be published. This study calculated the factors required to account for the perturbation of radiochromic film when used for small field dosimetry in lung equivalent material.

METHOD

Monte Carlo simulations were used to calculate dose to Gafchromic EBT3 film when placed inside a lung phantom. The beam simulated had a nominal energy of 6 MV and the field sizes simulated ranged from 10 × 10 mm² to 30 × 30 mm² . The lung density simulated was varied between 0.2 and 0.3 g/cm³ . Each simulation was repeated with the film replaced by lung material (the same as the surrounding medium), and the required correction factors for film dosimetry in lung ( D M e d , Q D D e t , Q ) were calculated by dividing the dose in lung by the dose in film.

RESULTS

For field sizes 30 × 30 mm² and larger, no correction factors were required. At a 20 × 20 mm² field size, small corrections were required, but were within the approximate accuracy of film dosimetry (~2%). For a 10 × 10 mm² field size, significant correction factors need to be applied (0.935 for lung density of 0.20 g/cm³ to 0.963 for lung density of 0.30 g/cm³ ). The values lower than one mean that the film is over-responding. At the "upstream" lung-water interface the correction factors were close to unity; while at the downstream interface the corrections required were marginally smaller to those at the center of lung. One centimeter or more away from the interfaces, the correction factor did not vary as a function distance from the interface (in the beam direction). Away from the central axis (perpendicular to the beam direction), the correction factors increased slightly (away from unity) as a function of off-axis distance, before abruptly changing direction at the penumbra, with the film actually under-responding by ~10% outside the field edges.

CONCLUSION

Accurate dosimetry of very small fields (15 × 15 mm² or smaller) using radiochromic film requires correction factors for the perturbation of the film on the surrounding lung material. This correction factor was as high as 6.5% for a 10 × 10 mm² field size and a density of 0.2 g/cm³ . This will increase if either the density or the field size decrease further. This correction factor does not vary as a function of depth in lung once charged particle equilibrium is established.

Patient-specific quality assurance on a Varian Halcyon linear accelerator using a PTW Octavius 4D device

This study investigates and validates the use of the Octavius 4D system for patient specific quality assurance on Halcyon, which is capable of rotating at 4 revolutions per minute (RPM). A commercially available PTW Octavius 4D system was used for this study which had a maximum rotation speed of 3 RPM. Initial validation included testing the accuracy of the inclinometer, percent depth doses (PDD), output factors, and dose profiles for selected static square fields. The same static fields were also subject to a gamma comparison with the TPS. This was followed by an evaluation of twelve clinical treatment plans and seven non-clinical plans with varying gantry rotation speeds. All testing was completed using detector array measurement times of 200 ms and 100 ms. Inclinometer accuracy was within 0.3° of actual gantry angle. Output factors varied less than 0.6%, PDD differences were no greater than 1.4%, and dose profile differences were less than 2.2%. Gamma pass rates for the static fields were 96.7% (2%/2mm) and 99.7% (3%/3mm). A prototype control unit, which had a maximum rotation speed of 4 RPM was also used to test the clinical and non-clinical plans. For the clinical plans, the mean gamma pass rates (2%/2mm) were 86.1% and 88.1% for the commercial unit and prototype unit respectively. Results using a measurement time of 200 ms were superior to those using 100 ms. For Halcyon deliveries greater than 3 RPM, worst case gamma results for the commercial unit were 28.6% compared to 98.5% using the prototype unit. Accurate patient specific quality assurance results can be obtained using the Octavius 4D system with a Halcyon linac, provided that the system measurement time is kept at 200 ms and the rotation speed of Halycon does not exceed 3 RPM. For higher RPM deliveries, an Octavius 4D unit with 4 RPM rotation capability is recommended.

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.

Density compensated diodes for small field dosimetry: comprehensive testing and implications for design

PURPOSE

In small megavoltage photon fields, the accuracies of an unmodified PTW 60017-type diode dosimeter and six diodes modified by adding airgaps of thickness 0.6-1.6 mm and diameter 3.6 mm have been comprehensively characterized experimentally and computationally. The optimally thick airgap for density compensation was determined, and detectors were micro-CT imaged to investigate differences between experimentally measured radiation responses and those predicted computationally.

METHODS

Detectors were tested on- and off-axis, at 5 and 15 cm depths in 6 and 15 MV fields ≥ 0.5 × 0.5 cm². Computational studies were carried out using the EGSnrc/BEAMnrc Monte Carlo radiation transport code. Experimentally, radiation was delivered using a Varian TrueBeam linac and doses absorbed by water were measured using Gafchromic EBT3 film and ionization chambers, and compared with diode readings. Detector response was characterized via the [Formula: see text] formalism, choosing a 4 × 4 cm² reference field.

RESULTS

For the unmodified 60017 diode, the maximum error in small field doses obtained from diode readings uncorrected by [Formula: see text] factors was determined as 11.9% computationally at +0.25 mm off-axis and 5 cm depth in a 15 MV 0.5 × 0.5 cm² field, and 11.7% experimentally at -0.30 mm off-axis and 5 cm depth in the same field. A detector modified to include a 1.6 mm thick airgap performed best, with maximum computationally and experimentally determined errors of 2.2% and 4.1%. The 1.6 mm airgap deepened the modified dosimeter's effective point of measurement by 0.5 mm. For some detectors significant differences existed between responses in small fields determined computationally and experimentally, micro-CT imaging indicating that these differences were due to within-tolerance variations in the thickness of an epoxy resin layer.

CONCLUSIONS

The dosimetric performance of a 60017 diode detector was comprehensively improved throughout 6 and 15 MV small photon fields via density compensation. For this approach to work well with good detector-to-detector reproducibility, tolerances on dense component dimensions should be reduced to limit associated variations of response in small fields, or these components should be modified to have more water-like densities.

Assessing the fit of 3D printed bolus from CT, optical scanner and photogrammetry methods

Bolus plays an important role in the radiation therapy of superficial lesions and the application of 3D printing to its design can improve fit and dosimetry. This study quantitatively compares the fits of boluses designed from different imaging modalities. A head phantom was imaged using three systems: a CT simulator, a 3D optical scanner, and an interchangeable lens camera. Nose boluses were designed and 3D printed from each modality. A 3D printed phantom with air gaps of known thicknesses was used to calibrate mean HU to measure air gaps of unknown thickness and assess the fit of each bolus on the head phantom. The bolus created from the optical scanner data resulted in the best fit, with a mean air gap of 0.16 mm. Smoothing of the CT bolus resulted in a more clinically suitable model, comparable to that from the optical scanner method. The bolus produced from the photogrammetry method resulted in air gaps larger than 1 mm in thickness. The use of optical scanner and photogrammetry models have many advantages over the conventional bolus-from-CT method, however workflow should be refined to ensure accuracy if implemented clinically.

2D monolithic silicon-diode array detectors in megavoltage photon beams: does the fabrication technology matter? A medical physicist's perspective

A family of prototype 2D monolithic silicon-diode array detectors (MP512, Duo, Octa) has been proposed by the Centre for Medical Radiation Physics, University of Wollongong (Australia) for relative dosimetry in small megavoltage photon beams. These detectors, which differ in the topology of their 512 sensitive volumes, were originally fabricated on bulk p-type substrates. More recently, they have also been fabricated on epitaxial p-type substrates. In the literature, their performance has been individually characterized for quality assurance (QA) applications. The present study directly assessed and compared that of a MP512-bulk and that of a MP512-epitaxial in terms of radiation hardness, long-term stability, response linearity with dose, dose per pulse and angular dependence. Their measurements of output factors, off-axis ratios and percentage depth doses in square radiation fields collimated by the jaws and produced by 6 MV and 10 MV flattened photon beams were then benchmarked against those by commercially available detectors. The present investigation was aimed at establishing, from a medical physicist's perspective, how the bulk and epitaxial fabrication technologies would affect the implementation of the MP512s into a QA protocol. Based on results, the MP512-epitaxial would offer superior radiation hardness, long-term stability and achievable uniformity and reproducibility of the response across the 2D active area.

CyberKnife® fixed cone and Iris™ defined small radiation fields: Assessment with a high-resolution solid-state detector array

Purpose

The challenges of accurate dosimetry for stereotactic radiotherapy (SRT) with small unflattened radiation fields have been widely reported in the literature. In this case, suitable dosimeters would have to offer a submillimeter spatial resolution. The CyberKnife^® (Accuray Inc., Sunnyvale, CA, USA) is an SRT‐dedicated linear accelerator (linac), which can deliver treatments with submillimeter positional accuracy using circular fields. Beams are delivered with the desired field size using fixed cones, the InCise™ multileaf collimator or a dynamic variable‐aperture Iris™ collimator. The latter, allowing for field sizes to be varied during treatment delivery, has the potential to decrease treatment time, but its reproducibility in terms of output factors (OFs) and dose profiles (DPs) needs to be verified.

Methods

A 2D monolithic silicon array detector, the “Octa”, was evaluated for dosimetric quality assurance (QA) for a CyberKnife system. OFs, DPs, percentage depth‐dose (PDD) and tissue maximum ratio (TMR) were investigated, and results were benchmarked against the PTW SRS diode. Cross‐plane, in‐plane and 2 diagonal dose profiles were measured simultaneously with high spatial resolution (0.3 mm). Monte Carlo (MC) simulations with a GEANT4 (GEometry ANd Tracking 4) tool‐kit were added to the study to support the experimental characterization of the detector response.

Results

For fixed cones and the Iris, for all field sizes investigated in the range between 5 and 60 mm diameter, OFs, PDDs, TMRs, and DPs in terms of FWHM measured by the Octa were accurate within 3% when benchmarked against the SRS diode and MC calculations.

Conclusions

The Octa was shown to be an accurate dosimeter for measurements with a 6 MV FFF beam delivered with a CyberKnife system. The detector enabled real‐time dosimetric verification for the variable aperture Iris collimator, yielding OFs and DPs consistent with those obtained with alternative methods.

A novel high-resolution 2D silicon array detector for small field dosimetry with FFF photon beams

PURPOSE

Flattening filter free (FFF) beams are increasingly being considered for stereotactic radiotherapy (SRT). For the first time, the performance of a monolithic silicon array detector under 6 and 10 MV FFF beams was evaluated. The dosimeter, named "Octa" and designed by the Centre for Medical Radiation Physics (CMRP), was tested also under flattened beams for comparison.

METHODS

Output factors (OFs), percentage depth-dose (PDD), dose profiles (DPs) and dose per pulse (DPP) dependence were investigated. Results were benchmarked against commercially available detectors for small field dosimetry.

RESULTS

The dosimeter was shown to be a 'correction-free' silicon array detector for OFs and PDD measurements for all the beam qualities investigated. Measured OFs were accurate within 3% and PDD values within 2% compared against the benchmarks. Cross-plane, in-plane and diagonal DPs were measured simultaneously with high spatial resolution (0.3 mm) and real time read-out. A DPP dependence (24% at 0.021 mGy/pulse relative to 0.278 mGy/pulse) was found and could be easily corrected for in the case of machine specific quality assurance applications.

CONCLUSIONS

Results were consistent with those for monolithic silicon array detectors designed by the CMRP and previously characterized under flattened beams only, supporting the robustness of this technology for relative dosimetry for a wide range of beam qualities and dose per pulses. In contrast to its predecessors, the design of the Octa offers an exhaustive high-resolution 2D dose map characterization, making it a unique real-time radiation detector for small field dosimetry for field sizes up to 3 cm side.

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

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

Determination of small-field correction factors for cylindrical ionization chambers using a semiempirical method

A semiempirical method based on the averaging effect of the sensitive volumes of different air-filled ionization chambers (ICs) was employed to approximate the correction factors for beam quality produced from the difference in the sizes of the reference field and small fields.We measured the output factors using several cylindrical ICs and calculated the correction factors using a mathematical method similar to deconvolution; in the method, we modeled the variable and inhomogeneous energy fluence function within the chamber cavity. The parameters of the modeled function and the correction factors were determined by solving a developed system of equations as well as on the basis of the measurement data and the geometry of the chambers. Further, Monte Carlo (MC) computations were performed using the Monaco® treatment planning system to validate the proposed method.The determined correction factors () were comparable to the values derived from the MC computations performed using Monaco®. For example, for a 6 MV photon beam and a field size of 1  ×  1 cm2, was calculated to be for a PTW 31010 chamber and for a PTW 31016 chamber. On the other hand, the values determined from the MC computations were 1.121 and 1.031, respectively; the difference between the proposed method and the MC computation is less than 2%. In addition, we determined the values for PTW 30013, PTW 31010, PTW 31016, IBA FC23-C, and IBA CC13 chambers as well.We devised a method for determining from both the measurement of the output factors and model-based mathematical computation. The proposed method can be useful in case the MC simulation would not be applicable for the clinical settings.

Small field detector correction factors: effects of the flattening filter for Elekta and Varian linear accelerators

Flattening filter‐free (FFF) beams are becoming the preferred beam type for stereotactic radiosurgery (SRS) and stereotactic ablative radiation therapy (SABR), as they enable an increase in dose rate and a decrease in treatment time. This work assesses the effects of the flattening filter on small field output factors for 6 MV beams generated by both Elekta and Varian linear accelerators, and determines differences between detector response in flattened (FF) and FFF beams. Relative output factors were measured with a range of detectors (diodes, ionization chambers, radiochromic film, and microDiamond) and referenced to the relative output factors measured with an air core fiber optic dosimeter (FOD), a scintillation dosimeter developed at Chris O'Brien Lifehouse, Sydney. Small field correction factors were generated for both FF and FFF beams. Diode measured detector response was compared with a recently published mathematical relation to predict diode response corrections in small fields. The effect of flattening filter removal on detector response was quantified using a ratio of relative detector responses in FFF and FF fields for the same field size. The removal of the flattening filter was found to have a small but measurable effect on ionization chamber response with maximum deviations of less than ±0.9% across all field sizes measured. Solid‐state detectors showed an increased dependence on the flattening filter of up to ±1.6%. Measured diode response was within ±1.1% of the published mathematical relation for all fields up to 30 mm, independent of linac type and presence or absence of a flattening filter. For 6 MV beams, detector correction factors between FFF and FF beams are interchangeable for a linac between FF and FFF modes, providing that an additional uncertainty of up to ±1.6% is accepted.

PACS number(s): 87.55.km, 87.56.bd, 87.56.Da

Dose rate dependence of the PTW 60019 microDiamond detector in high dose-per-pulse pulsed beams

Recombination effects can affect the detectors used for the dosimetry of radiotherapy fields. They are important when using ionization chambers, especially in liquid-filled ionization chambers, and should be corrected for. The introduction of flattening-filter-free accelerators increases the typical dose-per-pulse used in radiotherapy beams, which leads to more important recombination effects. Diamond detectors provide a good solution for the dosimetry and quality assurance of small radiotherapy fields, due to their low energy dependence and small volume. The group of Università di Roma Tor Vergata has developed a synthetic diamond detector, which is commercialized by PTW as microDiamond detector type 60019. In this work we present an experimental characterization of the collection efficiency of the microDiamond detector, focusing on high dose-per-pulse FFF beams. The collection efficiency decreases with dose-per-pulse, down to 0.978 at 2.2 mGy/pulse, following a Fowler-Attix-like curve. On the other hand, we have found no significant dependence of the collection efficiency on the pulse repetition frequency (or pulse period).

A methodological approach to reporting corrected small field relative outputs

PURPOSE

The goal of this work was to set out a methodology for measuring and reporting small field relative output and to assess the application of published correction factors across a population of linear accelerators.

METHODS AND MATERIALS

Measurements were made at 6 MV on five Varian iX accelerators using two PTW T60017 unshielded diodes. Relative output readings and profile measurements were made for nominal square field sizes of side 0.5 to 1.0 cm. The actual in-plane (A) and cross-plane (B) field widths were taken to be the FWHM at the 50% isodose level. An effective field size, defined as √FS eff=A · B, was calculated and is presented as a field size metric. FSeff was used to linearly interpolate between published Monte Carlo (MC) calculated [Formula in text] values to correct for the diode over-response in small fields.

RESULTS

The relative output data reported as a function of the nominal field size were different across the accelerator population by up to nearly 10%. However, using the effective field size for reporting showed that the actual output ratios were consistent across the accelerator population to within the experimental uncertainty of ± 1.0%. Correcting the measured relative output using [Formula in text] at both the nominal and effective field sizes produce output factors that were not identical but differ by much less than the reported experimental and/or MC statistical uncertainties.

CONCLUSIONS

In general, the proposed methodology removes much of the ambiguity in reporting and interpreting small field dosimetric quantities and facilitates a clear dosimetric comparison across a population of linacs.