Ionization chamber dosimetry of small photon fields: a Monte Carlo study on stopping-power ratios for radiosurgery and IMRT beams

Efficient dose-rate correction of silicon diode relative dose measurements

PURPOSE

Silicon diodes are often the detector of choice for relative dose measurements, particularly in the context of radiotherapy involving small photon beams. However, a major drawback lies in their dose-rate dependency. Although ionization chambers are often too large for small field output factor measurements, they are valuable instruments to provide reliable percent-depth dose curves in reference beams. The aim of this work is to propose a practical and accurate method for the characterization of silicon diode dose-rate dependence correction factors using ionization chamber measurements as a reference.

METHODS

The robustness of ionization chambers for percent-depth dose measurements is used to quantify the dose-rate dependency of a diode detector. A mathematical formalism, which exploits the error induced in percent-depth ionization curves for diodes by their dose-rate dependency, is developed to derive a dose-rate correction factor applicable to diode relative measurements. The method is based on the definition of the recombination correction factor given in the addendum to TG 51 and is applied to experimental measurements performed on a CyberKnife M6 radiotherapy unit using a PTW 60012 diode detector. A measurement-based validation is provided by comparing corrected percent-depth ionization curves to measurements performed with a PTW 60019 diamond detector which does not exhibit dose-rate dependence.

RESULTS

Results of dose-rate correction factors for percent-depth ionization curves, off-axis ratios, tissue-phantom ratios and small field output factors are coherent with the expected behavior of silicon diode detectors. For all considered setups and field sizes, the maximum correction and the maximum impact of the uncertainties induced by the correction are obtained for off-axis ratios for the 60 mm collimator, with a correction of 2.5% and an uncertainty of 0.34%. For output factors, corrections range from 0.33% to 0.82% for all field sizes considered, and increase with the reduction of the field size. Comparison of percent-depth ionization curves corrected for dose-rate and for in-depth beam quality variations illustrate excellent agreement with measurements performed using the diamond detector.

CONCLUSIONS

The proposed method allows the efficient and precise correction of the dose-rate dependence of silicon diode detectors in the context of clinical relative dosimetry. This article is protected by copyright. All rights reserved.

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.

Relationship between dose and stopping power values for electrons in skin and muscle tissues

Absorbed dose and stopping power of the target material are two important parameters for determining the radiation effects. In this work, the relationship between these two parameters has been investigated for electron beams incident on skin and muscle tissue. Absorbed dose was obtained by using the EGSnrc code and the stopping power values were calculated by considering the velocity-depended effective charge and mean excitation values. To obtain the relationship between absorbed dose and stopping power values, these parameters were graphed together and simple fitting functions have been obtained. The results obtained show that these parameters are linearly correlated with each other.

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.

Quality assurance of Cyberknife robotic stereotactic radiosurgery using an angularly independent silicon detector

Purpose

The aim of this work was to evaluate the use of an angularly independent silicon detector (edgeless diodes) developed for dosimetry in megavoltage radiotherapy for Cyberknife in a phantom and for patient quality assurance (QA).

Method

The characterization of the edgeless diodes has been performed on Cyberknife with fixed and IRIS collimators. The edgeless diode probes were tested in terms of basic QA parameters such as measurements of tissue‐phantom ratio (TPR), output factor and off‐axis ratio. The measurements were performed in both water and water‐equivalent phantoms. In addition, three patient‐specific plans have been delivered to a lung phantom with and without motion and dose measurements have been performed to verify the ability of the diodes to work as patient‐specific QA devices. The data obtained by the edgeless diodes have been compared to PTW 60016, SN edge, PinPoint ionization chamber, Gafchromic EBT3 film, and treatment planning system (TPS).

Results

The TPR measurement performed by the edgeless diodes show agreement within 2.2% with data obtained with PTW 60016 diode for all the field sizes. Output factor agrees within 2.6% with that measured by SN EDGE diodes corrected for their field size dependence. The beam profiles’ measurements of edgeless diodes match SN EDGE diodes with a measured full width half maximum (FWHM) within 2.3% and penumbra widths within 0.148 mm. Patient‐specific QA measurements demonstrate an agreement within 4.72% in comparison with TPS.

Conclusion

The edgeless diodes have been proved to be an excellent candidate for machine and patient QA for Cyberknife reproducing commercial dosimetry device measurements without need of angular dependence corrections. However, further investigation is required to evaluate the effect of their dose rate dependence on complex brain cancer dose verification.

Simple index for validity of the evaluation point for dosimetric verification results of intensity-modulated radiation therapy using a Farmer-type ionization chamber

PURPOSE

Based on a retrospective analysis, this study aims to develop a simple index for validity of the evaluation point for the dosimetric verification of intensity-modulated radiation therapy (IMRT).

METHODS

The results for the dosimetric verifications of a total of 69 IMRT plans were analyzed in this study. A Farmer-type ion chamber was used as a dose detector, and a solid water-equivalent phantom was used. Index values were obtained by dividing the difference between the maximum and minimum dosages by the mean dosage of the 69 plans, and the values were classified into five groups with index value <4, 4-8, 8-12, 12-16, and >16. A t-test was used to assess the statistical significance of the mean differences of the absolute values of the relative errors among these groups.

RESULTS

We found that there was no significant difference between the groups with index value <4 and 4-8 (p = 0.152); however, there were significant differences between the other groups (p < 0.01). In addition, when the index values were smaller than 8, the pass ratio of 3% tolerance was 96.2% and the pass ratio of 5% tolerance was 99.9%. We observed that the smaller the index value, the smaller the uncertainty of the dose measurement.

CONCLUSIONS

The results obtained in this study may prove to be useful for accurate dosimetric verifications of IMRTs when ion chambers are used.

Monte Carlo systems used for treatment planning and dose verification

General-purpose radiation transport Monte Carlo codes have been used for estimation of the absorbed dose distribution in external photon and electron beam radiotherapy patients since several decades. Results obtained with these codes are usually more accurate than those provided by treatment planning systems based on non-stochastic methods. Traditionally, absorbed dose computations based on general-purpose Monte Carlo codes have been used only for research, owing to the difficulties associated with setting up a simulation and the long computation time required. To take advantage of radiation transport Monte Carlo codes applied to routine clinical practice, researchers and private companies have developed treatment planning and dose verification systems that are partly or fully based on fast Monte Carlo algorithms. This review presents a comprehensive list of the currently existing Monte Carlo systems that can be used to calculate or verify an external photon and electron beam radiotherapy treatment plan. Particular attention is given to those systems that are distributed, either freely or commercially, and that do not require programming tasks from the end user. These systems are compared in terms of features and the simulation time required to compute a set of benchmark calculations.

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.

Analysis of small field percent depth dose and profiles: Comparison of measurements with various detectors and effects of detector orientation with different jaw settings

The advent of modern technologies in radiotherapy poses an increased challenge in the determination of dosimetric parameters of small fields that exhibit a high degree of uncertainty. Percent depth dose and beam profiles were acquired using different detectors in two different orientations. The parameters such as relative surface dose (D_S), depth of dose maximum (D_max), percentage dose at 10 cm (D₁₀), penumbral width, flatness, and symmetry were evaluated with different detectors. The dosimetric data were acquired for fields defined by jaws alone, multileaf collimator (MLC) alone, and by MLC while the jaws were positioned at 0, 0.25, 0.5, and 1.0 cm away from MLC leaf-end using a Varian linear accelerator with 6 MV photon beam. The accuracy in the measurement of dosimetric parameters with various detectors for three different field definitions was evaluated. The relative D_S(38.1%) with photon field diode in parallel orientation was higher than electron field diode (EFD) (27.9%) values for 1 cm ×1 cm field. An overestimation of 5.7% and 8.6% in D₁₀ depth were observed for 1 cm ×1 cm field with RK ion chamber in parallel and perpendicular orientation, respectively, for the fields defined by MLC while jaw positioned at the edge of the field when compared to EFD values in parallel orientation. For this field definition, the in-plane penumbral widths obtained with ion chamber in parallel and perpendicular orientation were 3.9 mm, 5.6 mm for 1 cm ×1 cm field, respectively. Among all detectors used in the study, the unshielded diodes were found to be an appropriate choice of detector for the measurement of beam parameters in small fields.

Small fields measurements with radiochromic films

The small fields in radiotherapy are widely used due to the development of techniques such as intensity-modulated radiotherapy and stereotactic radio surgery. The measurement of the dose distributions for small fields is a challenge. A perfect dosimeter should be independent of the radiation energy and the dose rate and should have a negligible volume effect. The radiochromic (RC) film characteristics fit well to these requirements. However, the response of RC films and their digitizing processes present a significant spatial inhomogeneity problem. The present work uses a method for two-dimensional (2D) measurement with RC films based on the reduction of the spatial inhomogeneity of both the film and the film digitizing process. By means of registering and averaging several measurements of the same field, the inhomogeneities are mostly canceled. Measurements of output factors (OFs), dose profiles (in-plane and cross-plane), and 2D dose distributions are presented. The field sizes investigated are 0.5 × 0.5 cm², 0.7 × 0.7 cm², 1 × 1 cm², 2 × 2 cm², 3 × 3 cm², 6 × 6 cm², and 10 × 10 cm² for 6 and 15 MV photon beams. The OFs measured with the RC film are compared with the measurements carried out with a PinPoint ionization chamber (IC) and a Semiflex IC, while the measured transversal dose profiles were compared with Monte Carlo simulations. The results obtained for the OFs measurements show a good agreement with the values obtained from RC films and the PinPoint and Semiflex chambers when the field size is greater or equal than 2 × 2 cm². These agreements give confidence on the accuracy of the method as well as on the results obtained for smaller fields. Also, good agreement was found between the measured profiles and the Monte Carlo calculated profiles for the field size of 1 × 1 cm². We expect, therefore, that the presented method can be used to perform accurate measurements of small fields.

Ionizing radiation detectors based on Ge-doped optical fibers inserted in resonant cavities

The measurement of ionizing radiation (IR) is a crucial issue in different areas of interest, from environmental safety and industrial monitoring to aerospace and medicine. Optical fiber sensors have recently proven good candidates as radiation dosimeters. Here we investigate the effect of IR on germanosilicate optical fibers. A piece of Ge-doped fiber enclosed between two fiber Bragg gratings (FBGs) is irradiated with gamma radiation generated by a 6 MV medical linear accelerator. With respect to other FBG-based IR dosimeters, here the sensor is only the bare fiber without any special internal structure. A near infrared laser is frequency locked to the cavity modes for high resolution measurement of radiation induced effects on the fiber optical parameters. In particular, we observe a variation of the fiber thermo-optic response with the radiation dose delivered, as expected from the interaction with Ge defect centers, and demonstrate a detection limit of 360 mGy. This method can have an impact in those contexts where low radiation doses have to be measured both in small volumes or over large areas, such as radiation therapy and radiation protection, while bare optical fibers are cheap and disposable.

The influence of field size on stopping-power ratios in- and out-of-field: quantitative data for the BrainLAB m3 micro-multileaf collimator

The objective of this work is to quantify the systematic errors introduced by the common assumption of invariant secondary electron spectra with changing field sizes, as relevant to stereotactic radiotherapy and other treatment modes incorporating small beam segments delivered with a linac‐based stereotactic unit. The EGSnrc/BEAMnrc Monte Carlo radiation transport code was used to construct a dosimetrically‐matched model of a Varian 600C linear accelerator with mounted BrainLAB micro‐multileaf collimator. Stopping‐power ratios were calculated for field sizes ranging from 6×6 mm2 up to the maximum (98×98 mm2), and differences between these and the reference field were computed. Quantitative stopping power data for the BrainLAB micro‐multileaf collimator has been compiled. Field size dependent differences to reference conditions increase with decreasing field size and increasing depth, but remain a fraction of a percent for all field sizes studied. However, for dosimetry outside the primary field, errors induced by the assumption of invariant electron spectra can be greater than 1%, increasing with field size. It is also shown that simplification of the Spencer‐Attix formulation by ignoring secondary electrons below the cutoff kinetic energy applied to the integration results in underestimation of stopping‐power ratios of about 0.3% (and is independent of field size and depth). This work is the first to quantify stopping powers from a BrainLAB micro‐multileaf collimator. Many earlier studies model simplified beams, ignoring collimator scatter, which is shown to significantly influence the spectrum. Importantly, we have confirmed that the assumption of unchanging electron spectra with varying field sizes is justifiable when performing (typical) in‐field dosimetry of stereotactic fields. Clinicians and physicists undertaking precise out‐of‐field measurements for the purposes of risk estimation, ought to be aware that the more pronounced spectral variation results in stopping powers (and hence doses) that differ more than for in‐field dosimetry.

PACS number: 87.10.RT; 87.53.Ly; 87.56.jf; 87.56.jk

A patient-specific quality assurance study on absolute dose verification using ionization chambers of different volumes in RapidArc treatments

The recalculation of 1 fraction from a patient treatment plan on a phantom and subsequent measurements have become the norms for measurement-based verification, which combines the quality assurance recommendations that deal with the treatment planning system and the beam delivery system. This type of evaluation has prompted attention to measurement equipment and techniques. Ionization chambers are considered the gold standard because of their precision, availability, and relative ease of use. This study evaluates and compares 5 different ionization chambers: phantom combinations for verification in routine patient-specific quality assurance of RapidArc treatments. Fifteen different RapidArc plans conforming to the clinical standards were selected for the study. Verification plans were then created for each treatment plan with different chamber-phantom combinations scanned by computed tomography. This includes Medtec intensity modulated radiation therapy (IMRT) phantom with micro-ionization chamber (0.007 cm(3)) and pinpoint chamber (0.015 cm(3)), PTW-Octavius phantom with semiflex chamber (0.125 cm(3)) and 2D array (0.125 cm(3)), and indigenously made Circular wax phantom with 0.6 cm(3) chamber. The measured isocenter absolute dose was compared with the treatment planning system (TPS) plan. The micro-ionization chamber shows more deviations when compared with semiflex and 0.6 cm(3) with a maximum variation of -4.76%, -1.49%, and 2.23% for micro-ionization, semiflex, and farmer chambers, respectively. The positive variations indicate that the chamber with larger volume overestimates. Farmer chamber shows higher deviation when compared with 0.125 cm(3). In general the deviation was found to be <1% with the semiflex and farmer chambers. A maximum variation of 2% was observed for the 0.007 cm(3) ionization chamber, except in a few cases. Pinpoint chamber underestimates the calculated isocenter dose by a maximum of 4.8%. Absolute dose measurements using the semiflex ionization chamber with intermediate volume (0.125 cm(3)) shows good agreement with the TPS calculated among the detectors used in this study. Positioning is very important when using smaller volume chambers because they are more sensitive to geometrical errors within the treatment fields. It is also suggested to average the dose over the sensitive volume for larger-volume chambers. The ionization chamber-phantom combinations used in this study can be used interchangeably for routine RapidArc patient-specific quality assurance with a satisfactory accuracy for clinical practice.

On the quantification of the dosimetric accuracy of collapsed cone convolution superposition (CCCS) algorithm for small lung volumes using IMRT

Specialized techniques that make use of small field dosimetry are common practice in today's clinics. These new techniques represent a big challenge to the treatment planning systems due to the lack of lateral electronic equilibrium. Because of this, the necessity of planning systems to overcome such difficulties and provide an accurate representation of the true value is of significant importance. Pinnacle3 is one such planning system. During the IMRT optimization process, Pinnacle3 treatment planning system allows the user to specify a minimum segment size which results in multiple beams composed of several subsets of different widths. In this study, the accuracy of the engine dose calculation, collapsed cone convolution superposition algorithm (CCCS) used by Pinnacle3, was quantified by Monte Carlo simulations, ionization chamber, and Kodak extended dose range film (EDR2) measurements for 11 SBRT lung patients. Lesions were < 3.0 cm in maximal diameter and <27.0cm3 in volume. The Monte Carlo EGSnrc\BEAMnrc and EGS4\MCSIM were used in the comparison. The minimum segment size allowable during optimization had a direct impact on the number of monitor units calculated for each beam. Plans with the smallest minimum segment size (0.1 cm2 to 2.0 cm2) had the largest number of MUs. Although PTV coverage remained unaffected, the segment size did have an effect on the dose to the organs at risk. Pinnacle3-calculated PTV mean doses were in agreement with Monte Carlo-calculated mean doses to within 5.6% for all plans. On average, the mean dose difference between Monte Carlo and Pinnacle3 for all 88 plans was 1.38%. The largest discrepancy in maximum dose was 5.8%, and was noted for one of the plans using a minimum segment size of 1.0 cm2. For minimum dose to the PTV, a maximum discrepancy between Monte Carlo and Pinnacle3 was noted of 12.5% for a plan using a 6.0 cm2 minimum segment size. Agreement between point dose measurements and Pinnacle3-calculated doses were on average within 0.7% in both phantoms. The profiles show a good agreement between Pinnacle3, Monte Carlo, and EDR2 film. The gamma index and the isodose lines support the result.

Current status and future perspective of flattening filter free photon beams

PURPOSE

Flattening filters (FFs) have been considered as an integral part of the treatment head of a medical accelerator for more than 50 years. The reasons for the longstanding use are, however, historical ones. Advanced treatment techniques, such as stereotactic radiotherapy or intensity modulated radiotherapy have stimulated the interest in operating linear accelerators in a flattening filter free (FFF) mode. The current manuscript reviews treatment head physics of FFF beams, describes their characteristics and the resulting potential advantages in their medical use, and closes with an outlook.

METHODS

A number of dosimetric benefits have been determined for FFF beams, which range from increased dose rate and dose per pulse to favorable output ratio in-air variation with field size, reduced energy variation across the beam, and reduced leakage and out-of-field dose, respectively. Finally, the softer photon spectrum of unflattened beams has implications on imaging strategies and radiation protection.

RESULTS

The dosimetric characteristics of FFF beams have an effect on treatment delivery, patient comfort, dose calculation accuracy, beam matching, absorbed dose determination, treatment planning, machine specific quality assurance, imaging, and radiation protection. When considering conventional C-arm linacs in a FFF mode, more studies are needed to specify and quantify the clinical advantages, especially with respect to treatment plan quality and quality assurance.

CONCLUSIONS

New treatment units are already on the market that operate without a FF or can be operated in a dedicated clinical FFF mode. Due to the convincing arguments of removing the FF, it is expected that more vendors will offer dedicated treatment units for advanced photon beam therapy in the near future. Several aspects related to standardization, dosimetry, treatment planning, and optimization need to be addressed in more detail in order to facilitate the clinical implementation of unflattened beams.

A contemporary review of stereotactic radiotherapy: inherent dosimetric complexities and the potential for detriment

OBJECTIVE

The advantages of highly localised, conformal treatments achievable with stereotactic radiotherapy (SRT) are increasingly being extended to extracranial sites as stereotactic body radiotherapy with advancements in imaging and beam collimation. One of the challenges in stereotactic treatment lies in the significant complexities associated with small field dosimetry and dose calculation. This review provides a comprehensive overview of the complexities associated with stereotactic radiotherapy and the potential for detriment.

METHODS

This study is based on a comprehensive review of literature accessible via PubMed and other sources, covering stereotactic radiotherapy, small-field dosimetry and dose calculation.

FINDINGS

Several key issues were identified in the literature. They pertain to dose prescription, dose measurement and dose calculation within and beyond the treatment field. Field-edge regions and penumbrae occupy a significant portion of the total field size. Spectral and dosimetric characteristics are difficult to determine and are compounded by effects of tissue inhomogeneity. Measurement of small-fields is made difficult by detector volume averaging and energy response. Available dosimeters are compared, and emphasis is given to gel dosimetry which offers the greatest potential for three-dimensional small-field dosimetry. The limitations of treatment planning system algorithms as applied to small-fields (particularly in the presence of heterogeneities) is explained, and a review of Monte Carlo dose calculation is provided, including simplified treatment planning implementations. Not incorporated into treatment planning, there is evidence that far from the primary field, doses to patients (and corresponding risks of radiocarcinogenesis) from leakage/scatter in SRT are similar to large fields.

CONCLUSIONS

Improved knowledge of dosimetric issues is essential to the accurate measurement and calculation of dose as well as the interpretation and assessment of planned and delivered treatments. This review highlights such issues and the potential benefit that may be gained from Monte Carlo dose calculation and verification via three-dimensional dosimetric methods (such as gel dosimetry) being introduced into routine clinical practice.

Statistical analysis of IMRT dosimetry quality assurance measurements for local delivery guideline

Purpose

To establish our institutional guideline for IMRT delivery, we statistically evaluated the results of dosimetry quality assurance (DQA) measurements and derived local confidence limits using the concept confidence limit of |mean|+1.96σ.

Materials and methods

From June 2006 to March 2009, 206 patients with head and neck cancer, prostate cancer, liver cancer, or brain tumor were treated using LINAC-based IMRT technique. In order to determine site specific DQA tolerances at a later stage, a hybrid plan with the same fluence maps as in the treatment plan was generated on CT images of a cylindrical phantom of acryl. Points of measurement using a 0.125 cm³ ion-chamber were typically located in the region of high and uniform doses. The planar dose distributions perpendicular to the central axis were measured by using a diode array in solid water with all fields delivered, and assessed using gamma criteria of 3%/3 mm. The mean values and standard deviations were used to develop the local confidence and tolerance limits. The dose differences and gamma pass rates for the different treatment sites were also evaluated in terms of total monitor uints (MU), MU/cGy, and the number of PTV's pieces.

Results

The mean values and standard deviations of ion-chamber dosimetry differences between calculated and measured doses were -1.6 ± 1.2% for H&N cancer, -0.4 ± 1.2% for prostate and abdominal cancer, and -0.6 ± 1.5% for brain tumor. Most of measured doses (92.2%) agreed with the calculated doses within a tolerance limit of ±3% recommended in the literature. However, we found some systematic under-dosage for all treatment sites. The percentage of points passing the gamma criteria, averaged over all treatment sites was 97.3 ± 3.7%. The gamma pass rate and the agreement of ion-chamber dosimetry generally decreased with increasing the number of PTV's pieces, the degree of modulation (MU/cGy), and the total MU beyond 700. Our local confidence limits were comparable to those of AAPM TG 119 and ESTRO guidelines that were provided as a practical baseline for center-to-center commissioning comparison. Thus, our institutional confidence and action limits for IMRT delivery were set into the same levels of those guidelines.

Discussion and Conclusions

The systematic under-dosage were corrected by tuning up the MLC-related factors (dosimetric gap and transmission) in treatment planning system (TPS) and further by incorporating the tongue-and groove effect into TPS. Institutions that have performed IMRT DQA measurements over a certain period of time need to analyze their accrued DQA data. We confirmed the overall integrity of our IMRT system and established the IMRT delivery guideline during this procedure. Dosimetric corrections for the treatment plans outside of the action level can be suggested only with such rigorous DQA and statistical analysis.

From the limits of the classical model of sensitometric curves to a realistic model based on the percolation theory for GafChromic EBT films

PURPOSE

Modern radiotherapy uses complex treatments that necessitate more complex quality assurance procedures. As a continuous medium, GafChromic EBT films offer suitable features for such verification. However, its sensitometric curve is not fully understood in terms of classical theoretical models. In fact, measured optical densities and those predicted by the classical models differ significantly. This difference increases systematically with wider dose ranges. Thus, achieving the accuracy required for intensity-modulated radiotherapy (IMRT) by classical methods is not possible, plecluding their use. As a result, experimental parametrizations, such as polynomial fits, are replacing phenomenological expressions in modern investigations. This article focuses on identifying new theoretical ways to describe sensitometric curves and on evaluating the quality of fit for experimental data based on four proposed models.

METHODS

A whole mathematical formalism starting with a geometrical version of the classical theory is used to develop new expressions for the sensitometric curves. General results from the percolation theory are also used. A flat-bed-scanner-based method was chosen for the film analysis. Different tests were performed, such as consistency of the numeric results for the proposed model and double examination using data from independent researchers.

RESULTS

Results show that the percolation-theory-based model provides the best theoretical explanation for the sensitometric behavior of GafChromic films. The different sizes of active centers or monomer crystals of the film are the basis of this model, allowing acquisition of information about the internal structure of the films. Values for the mean size of the active centers were obtained in accordance with technical specifications. In this model, the dynamics of the interaction between the active centers of GafChromic film and radiation is also characterized by means of its interaction cross-section value.

CONCLUSIONS

The percolation model fulfills the accuracy requirements for quality-control procedures when large ranges of doses are used and offers a physical explanation for the film response.

The change of response of ionization chambers in the penumbra and transmission regions: impact for IMRT verification

Significant deviations from the expected dose have been reported in the absolute dosimetry validation of an intensity modulated radiation therapy treatment when individual segments are analyzed. However, when full treatment is considered and all segment doses are added together, these discrepancies fade out, leading to overall dose deviations below a 5% action level. This contradictory behavior may be caused by a partial compensation between detector over-responding and under-responding for measurement conditions far from radiation equilibrium. We consider three treatment verification scenarios that may lead to ionization chamber miss-responding, namely: narrow beam irradiation, field penumbra location and multi-leaf collimator transmission contribution. In this work we have analyzed the response of three different ionization chambers with different active volume under these conditions by means of Monte Carlo (MC) simulation methods. Correction factors needed to convert the detector readout into actual dose to water were calculated by inserting the specific detector geometry (carefully modeled) into the simulations. This procedure required extensive use of parallel computing resources in order to achieve the desired level of uncertainty in the final results. The analysis of the simulations shows the relative contribution of each of the three previously mentioned miss-responding scenarios. Additionally, we provide some evidence on dose deviation compensation in multi-segment radiotherapy treatment verification.

Determination of zero-field size percent depth doses and tissue maximum ratios for stereotactic radiosurgery and IMRT dosimetry: Comparison between experimental measurements and Monte Carlo simulation

In this study, zero-field percent depth dose (PDD) and tissue maximum ratio (TMR) for 6 MV x rays have been determined by extrapolation from dosimetric measurements over the field size range 1 x 1-10 x 10 cm2. The key to small field dosimetry is the selection of a proper dosimeter for the measurements, as well as the alignment of the detector with the central axis (CAX) of beam. The measured PDD results are compared with those obtained from Monte Carlo (MC) simulation to examine the consistency and integrity of the measured data from which the zero-field PDD is extrapolated. Of the six most commonly used dosimeters in the clinic, the stereotactic diode field detector (SFD), the PTW Pinpoint, and the Exradin A14 are the most consistent and produce results within 2% of each other over the entire field size range 1 x 1-40 x 40 cm2. Although the diamond detector has the smallest sensitive volume, it is the least stable and tends to disagree with all other dosimeters by more than 10%. The zero-field PDD data extrapolated from larger field measurements obtained with the SFD are in good agreement with the MC results. The extrapolated and MC data agree within 2.5% over the clinical depth range (dmax-30 cm), when the MC data for the zero field are derived from a 1 X 1 cm2 field simulation using a miniphantom (1 x 1 x 48 cm3). The agreement between the measured PDD and the MC data based on a full phantom (48 x 48 x 48 cm3) simulation is fairly good within 1% at shallow depths to approximately 5% at 30 cm. Our results seem to indicate that zero-field TMR can be accurately calculated from PDD measurements with a proper choice of detector and a careful alignment of detector axis with the CAX.

Uncertainty estimation in intensity-modulated radiotherapy absolute dosimetry verification

PURPOSE

Intensity-modulated radiotherapy (IMRT) represents an important method for improving RT. The IMRT relative dosimetry checks are well established; however, open questions remain in reference dosimetry with ionization chambers (ICs). The main problem is the departure of the measurement conditions from the reference ones; thus, additional uncertainty is introduced into the dose determination. The goal of this study was to assess this effect systematically.

METHODS AND MATERIALS

Monte Carlo calculations and dosimetric measurements with five different detectors were performed for a number of representative IMRT cases, covering both step-and-shoot and dynamic delivery.

RESULTS

Using ICs with volumes of about 0.125 cm(3) or less, good agreement was observed among the detectors in most of the situations studied. These results also agreed well with the Monte Carlo-calculated nonreference correction factors (c factors). Additionally, we found a general correlation between the IC position relative to a segment and the derived correction factor c, which can be used to estimate the expected overall uncertainty of the treatment.

CONCLUSION

The increase of the reference dose relative standard uncertainty measured with ICs introduced by nonreference conditions when verifying an entire IMRT plan is about 1-1.5%, provided that appropriate small-volume chambers are used. The overall standard uncertainty of the measured IMRT dose amounts to about 2.3%, including the 0.5% of reproducibility and 1.5% of uncertainty associated with the beam calibration factor. Solid state detectors and large-volume chambers are not well suited to IMRT verification dosimetry because of the greater uncertainties. An action level of 5% is appropriate for IMRT verification. Greater discrepancies should lead to a review of the dosimetric procedure, including visual inspection of treatment segments and energy fluence.

Monte Carlo modeling of the ModuLeaf miniature MLC for small field dosimetry and quality assurance of the clinical treatment planning system

The purpose of this investigation was the verification of both the measured data and quality of the implementation of the add-on ModuLeaf miniature multileaf collimator (ML mMLC) into the clinical treatment planning system for conformal stereotactic radiosurgery treatment. To this end the treatment head with ML mMLC was modeled in the BEAMnrc Monte Carlo (MC) code. The 6 MV photon beams used in the setup were first benchmarked with a set of measurements. A total ML mMLC transmission of 1.13% of the 10 x 10 cm2 open field dose was measured and reproduced with the BEAMnrc/DOSXYZnrc code. Correspondence between calculated and measured output factors (OFs) was within 2%. Correspondence between MC and measured profiles was within 2% dose and 2 mm distance, only for the smallest 0.5 x 0.5 cm2 field the results were within 3% dose. In the next step, the MC model was compared with Gafchromic film measurements and Pinnacle(3) 7.4 f (convolution superposition algorithm) calculated dose distributions, using a gamma evaluation comparison, for a multi-beam patient setup delivered to a Lucytrade mark phantom. The gamma evaluation of the MC versus Gafchromic film resulted in 3.4% of points not fulfilling gamma <or= 1 for a 2%/2 mm criterion, the Pinnacle(3) 7.4 f versus Gafchromic results 3.8% and Pinnacle versus MC less than 1%. For specific patients with lesions of 8 cc and 0.2 cc, Monte Carlo and Pinnacle simulations of the plans were performed and compared using DVH evaluation. DVHs corresponded within 2% dose and 2% volume.

Measurement of output factors for small photon beams

A variety of detectors and procedures for the measurement of small field output factors are discussed in the current literature. Different detectors with or without corrections are recommended. Correction factors are often derived by Monte Carlo methods, where the bias due to approximations in the model is difficult to judge. Over that, results appear to be contradictory in some cases. In this work, output factors were measured for field sizes from 4 mm up to 180 mm side length with different detectors. A simple linear correction for the energy response of solid state detectors is proposed. This led to identical values down to 8 mm field size, as long as the size of the detector is small against the field size. The correction was of the order of a few percent. For a shielded silicon diode it was well below 1%. A physically meaningful function is proposed in order to calculate output factors for arbitrary field sizes with high accuracy.