A dual source photon beam model used in convolution/superposition dose calculations for clinical megavoltage x-ray beams

A phase space model of a Versa HD linear accelerator for application to Monte Carlo dose calculation in a real-time adaptive workflow

PURPOSE

This study aims to develop and validate a simple geometric model of the accelerator head, from which a particle phase space can be calculated for application to fast Monte Carlo dose calculation in real-time adaptive photon radiotherapy. With this objective in view, the study investigates whether the phase space model can facilitate dose calculations which are compatible with those of a commercial treatment planning system, for convenient interoperability.

MATERIALS AND METHODS

A dual-source model of the head of a Versa HD accelerator (Elekta AB, Stockholm, Sweden) was created. The model used parameters chosen to be compatible with those of 6-MV flattened and 6-MV flattening filter-free photon beams in the RayStation treatment planning system (RaySearch Laboratories, Stockholm, Sweden). The phase space model was used to calculate a photon phase space for several treatment plans, and the resulting phase space was applied to the Dose Planning Method (DPM) Monte Carlo dose calculation algorithm. Simple fields and intensity-modulated radiation therapy (IMRT) treatment plans for prostate and lung were calculated for benchmarking purposes and compared with the convolution-superposition dose calculation within RayStation.

RESULTS

For simple square fields in a water phantom, the calculated dose distribution agrees to within ±2% with that from the commercial treatment planning system, except in the buildup region, where the DPM code does not model the electron contamination. For IMRT plans of prostate and lung, agreements of ±2% and ±6%, respectively, are found, with slightly larger differences in the high dose gradients.

CONCLUSIONS

The phase space model presented allows convenient calculation of a phase space for application to Monte Carlo dose calculation, with straightforward translation of beam parameters from the RayStation beam model. This provides a basis on which to develop dose calculation in a real-time adaptive setting.

A polar-coordinate-based pencil beam algorithm for VMAT dose computation with high-resolution gantry angle sampling

PURPOSE

Most commercially available treatment planning systems (TPSs) approximate the continuous delivery of volumetric modulated arc therapy (VMAT) plans with a series of discretized static beams for treatment planning, which can make VMAT dose computation extremely inefficient. In this study, we developed a polar-coordinate-based pencil beam (PB) algorithm for efficient VMAT dose computation with high-resolution gantry angle sampling that can improve the computational efficiency and reduce the dose discrepancy due to the angular under-sampling effect.

METHODS AND MATERIALS

6 MV 1 × 1 mm² pencil beams were simulated on a uniform cylindrical phantom under an EGSnrc Monte Carlo (MC) environment. The MC-generated PB kernels were collected in the polar coordinate system for each bixel on a 40 × 40 cm² fluence map and subsequently fitted via a series of Gaussians. The fluence was calculated using a detectors' eye view with off-axis and MLC transmission factors corrected. Doses of VMAT arc on the phantom were computed by summing the convolution results between the corresponding PB kernels and fluence for each bixel in the polar coordinate system. The convolution was performed using Fast Fourier Transform to expedite the computing speed. The calculated doses were converted to the Cartesian coordinate system and compared with the reference dose computed by a collapsed cone convolution (CCC) algorithm of the TPS. A heterogeneous phantom was created to study the heterogeneity corrections using the proposed algorithm. Ten VMAT arcs were included to evaluate the algorithm performance. Gamma analysis and computation complexity theory were used to measure the dosimetric accuracy and computational efficiency, respectively.

RESULTS

The dosimetric comparisons on the homogeneous phantom between the proposed PB algorithm and the CCC algorithm for ten VMAT arcs demonstrate that the proposed algorithm can achieve a dosimetric accuracy comparable to that of the CCC algorithm with average gamma passing rates of 96% (2%/2mm) and 98% (3%/3mm). In addition, the proposed algorithm can provide better computational efficiency for VMAT dose computation using a PC equipped with a 4-core processor, compared to the CCC algorithm utilizing a dual 10-core server. Moreover, the computation complexity theory reveals that the proposed algorithm has a great advantage with regard to computational efficiency for VMAT dose computation on homogeneous medium, especially when a fine angular sampling rate is applied. This can support a reduction in dose errors from the angular under-sampling effect by using a finer angular sampling rate, while still preserving a practical computing speed. For dose calculation on the heterogeneous phantom, the proposed algorithm with heterogeneity corrections can still offer a reasonable dosimetric accuracy with comparable computational efficiency to that of the CCC algorithm.

CONCLUSIONS

We proposed a novel polar-coordinate-based pencil beam algorithm for VMAT dose computation that enables a better computational efficiency while maintaining clinically acceptable dosimetric accuracy and reducing dose error caused by the angular under-sampling effect. It also provides a flexible VMAT dose computation structure that allows adjustable sampling rates and direct dose computation in regions of interest, which makes the algorithm potentially useful for clinical applications such as independent dose verification for VMAT patient-specific QA. This article is protected by copyright. All rights reserved.

An automated optimization strategy to design collimator geometry for small field radiation therapy systems

PURPOSE

To develop an automated optimization strategy to facilitate collimator design for small-field radiotherapy systems.

METHODS

We developed an objective function that links the dose profile characteristics (FWHM, penumbra, and central dose rate) and the treatment head geometric parameters (collimator thickness/radii, source-to-distal-collimator distance[SDC]) for small-field radiotherapy systems. We performed optimization using a downhill simplex algorithm. We applied this optimization strategy to a linac-based radiosurgery system to determine the optimal geometry of four pencil-beam collimators to produce 5, 10, 15, and 20mm diameter photon beams (from a 6.7MeV, 2.1mmFWHM electron beam). Two different optimizations were performed to prioritize minimum penumbra or maximum central dose rate for each beam size. We compared the optimized geometric parameters and dose distributions to an existing clinical system (CyberKnife).

RESULTS

When minimum penumbra was prioritized, using the same collimator thickness and SDC (40cm) as a CyberKnife system, the optimized collimator upstream and downstream radii agreed with the CyberKnife system within 3-14%, the optimized output factors agreed within 0-8%, and the optimized transverse and percentage depth dose profiles matched those of the CyberKnife with the penumbras agreeing within 2%. However, when maximum dose rate was prioritized, allowing both the collimator thickness and SDC to change, the central dose rate for larger collimator sizes (10, 15, 20mm) could be increased by about 1.5-2 times at the cost of 1.5-2 times larger penumbras. No further improvement in central dose rate for the 5mm beam size could be achieved.

CONCLUSIONS

We developed an automated optimization strategy to design the collimator geometry for small-field radiation therapy systems. Using this strategy, the penumbra-prioritized dose distribution and geometric parameters agree well with the CyberKnife system as an example, suggesting that this system was designed to prioritize sharp penumbra. This represents proof-of-principle that an automated optimization strategy may apply to more complex collimator designs with multiple optimization parameters.

A preliminary study of a photon dose calculation algorithm using a convolutional neural network

The aim of dose calculation algorithm research is to improve the calculation accuracy while maximizing the calculation efficiency. In this study, the three-dimensional distribution of total energy release per unit mass (TERMA) and the electron density (ED) distribution are considered inputs in a method for calculating the three-dimensional dose distribution based on a convolutional neural network (CNN). Attempts are made to improve the efficiency of the collapsed cone convolution/superposition (CCCS) algorithm while providing an approach to improve the efficiency of other traditional dose calculation algorithms. Twelve sets of computed tomography (CT) images were employed for training. Data sets were generated by the CCCS algorithm with a random beam configuration. For each monoenergetic photon model, 7500 samples were generated for the training set, and 1500 samples were generated for the validation set. Training occurred for 0.5 MeV, 1 MeV, 2 MeV, 3 MeV, 4 MeV, 5 MeV, and 6 MeV monoenergetic photon models. To evaluate the usability under linac conditions, a comparison between CCCS and CNN-Dose was performed for the Mohan 6-MV spectrum for 12 additional new sets of CT images with different anatomies. A total of 1512 test samples were generated. For all anatomies, the mean value, 95% lower confidence limit (LCL) and 95% upper confidence limit (UCL) were 99.56%, 99.51% and 99.61%, respectively, at the 3%/2 mm criteria. The mean value, 95% LCL and 95% UCL were 98.57%, 98.46% and 98.67%, respectively, at the 2%/2 mm criteria. The results meet the relevant clinical requirements. In the proposed methods, the dose distribution of clinical energy can be obtained by TERMA, and the electronic density can be obtained with a CNN. This method can also be used for other traditional dose algorithms and displays potential in treatment planning, adaptive radiation therapy, and in vivo verification.

Comparison of gamma index based on dosimetric error and clinically relevant dose-volume index based on three-dimensional dose prediction in breast intensity-modulated radiation therapy

Background

Measurement-guided dose reconstruction has lately attracted significant attention because it can predict the delivered patient dose distribution. Although the treatment planning system (TPS) uses sophisticated algorithm to calculate the dose distribution, the calculation accuracy depends on the particular TPS used. This study aimed to investigate the relationship between the gamma passing rate (GPR) and the clinically relevant dose–volume index based on the predicted 3D patient dose distribution derived from two TPSs (XiO, RayStation).

Methods

Twenty-one breast intensity-modulated radiation therapy plans were inversely optimized using XiO. With the same plans, both TPSs calculated the planned dose distribution. We conducted per-beam measurements on the coronal plane using a 2D array detector and analyzed the difference in 2D GPRs between the measured and planned doses by commercial software. Using in-house software, we calculated the predicted 3D patient dose distribution and derived the predicted 3D GPR, the predicted per-organ 3D GPR, and the predicted clinically relevant dose–volume indices [dose–volume histogram metrics and the value of the tumor-control probability/normal tissue complication probability of the planning target volume and organs at risk]. The results derived from XiO were compared with those from RayStation.

Results

While the mean 2D GPRs derived from both TPSs were 98.1% (XiO) and 100% (RayStation), the mean predicted 3D GPRs of ipsilateral lung (73.3% [XiO] and 85.9% [RayStation]; p < 0.001) had no correlation with 2D GPRs under the 3% global/3 mm criterion. Besides, this significant difference in terms of referenced TPS between XiO and RayStation could be explained by the fact that the error of predicted V_5Gy of ipsilateral lung derived from XiO (29.6%) was significantly larger than that derived from RayStation (− 0.2%; p < 0.001).

Conclusions

GPR is useful as a patient quality assurance to detect dosimetric errors; however, it does not necessarily contain detailed information on errors. Using the predicted clinically relevant dose–volume indices, the clinical interpretation of dosimetric errors can be obtained. We conclude that a clinically relevant dose–volume index based on the predicted 3D patient dose distribution could add to the clinical and biological considerations in the GPR, if we can guarantee the dose calculation accuracy of referenced TPS.

Monte-Carlo simulation of the Siemens Artiste linear accelerator flat 6 MV and flattening-filter-free 7 MV beam line

The aim of our work is to provide the up-to-now missing information on the Siemens Artiste FFF 7 MV beam line using a Monte-Carlo model fit to the realistic dosimetric measurements at the linear accelerator in clinical use at our department. The main Siemens Artiste 6MV and FFF 7MV beams were simulated using the Geant4 toolkit. The simulations were compared with the measurements with an ionization chamber in a water phantom to verify the validation of simulation and tuning the primary electron parameters. Hereafter, other parameters such as surface dose, spectrum, electron contamination, symmetry, flatness/unflatness, slope, and characteristic off-axis changes were discussed for both Flat and FFF mode. The mean electron energy for the FFF beam was 8.8 MeV and 7.5 MeV for Flat 6 MV, the spread energy and spot size of the selected Gaussian distribution source were 0.4 MeV and 1mm, respectively. The dose rate of the FFF beam was 2.8 (2.96) times higher than for the flattened beam for a field size of 10×10 (20×20) cm2. The electron contamination has significant contribution to the surface dose especially for the flattened beam. The penumbra, surface dose and the mean energy of photons decrease by removing the flattening filter. Finally, the results show that off-axis changes have no strong effect on the mean energy of FFF beams, while this effect was more considerable for the flattened beam.

Characterization of the Cerenkov scatter function: a convolution kernel for Cerenkov light dosimetry

Cerenkov light is created in clinical applications involving high-energy radiation such as in radiation therapy. There is considerable interest in using Cerenkov light as a means to perform in vivo dosimetry during radiation therapy; however, a better understanding of the light-to-dose relationship is needed. One such method to solve this relationship is that of a deconvolution formulation, which relies on the Cerenkov scatter function (CSF). The CSF describes the creation of Cerenkov photons by a pencil beam of high-energy radiation, and the subsequent scattering that occurs before emission from the irradiated medium surface. This study investigated the dependence of the CSF on common radiation beam parameters (beam energy and incident angle) and the type of irradiated medium. An analytical equation with fitting coefficients of the CSF was obtained for common beam energies in a stratified skin model and optical phantom. Perturbation analysis was performed to investigate the dependence of the deconvolved Cerenkov images on the full-width at half-maximum and amplitude of the CSF. The irradiated material and beam angle had a large impact on the deconvolution process, whereas the beam energy had little effect.

A photon source model based on particle transport in a parameterized accelerator structure for Monte Carlo dose calculations

PURPOSE

An accurate source model of a medical linear accelerator is essential for Monte Carlo (MC) dose calculations. This study aims to propose an analytical photon source model based on particle transport in parameterized accelerator structures, focusing on a more realistic determination of linac photon spectra compared to existing approaches.

METHODS

We designed the primary and secondary photon sources based on the photons attenuated and scattered by a parameterized flattening filter. The primary photons were derived by attenuating bremsstrahlung photons based on the path length in the filter. Conversely, the secondary photons were derived from the decrement of the primary photons in the attenuation process. This design facilitates these sources to share the free parameters of the filter shape and be related to each other through the photon interaction in the filter. We introduced two other parameters of the primary photon source to describe the particle fluence in penumbral regions. All the parameters are optimized based on calculated dose curves in water using the pencil-beam-based algorithm. To verify the modeling accuracy, we compared the proposed model with the phase space data (PSD) of the Varian TrueBeam 6 and 15 MV accelerators in terms of the beam characteristics and the dose distributions. The EGS5 Monte Carlo code was used to calculate the dose distributions associated with the optimized model and reference PSD in a homogeneous water phantom and a heterogeneous lung phantom. We calculated the percentage of points passing 1D and 2D gamma analysis with 1%/1 mm criteria for the dose curves and lateral dose distributions, respectively.

RESULTS

The optimized model accurately reproduced the spectral curves of the reference PSD both on- and off-axis. The depth dose and lateral dose profiles of the optimized model also showed good agreement with those of the reference PSD. The passing rates of the 1D gamma analysis with 1%/1 mm criteria between the model and PSD were 100% for 4 × 4, 10 × 10, and 20 × 20 cm² fields at multiple depths. For the 2D dose distributions calculated in the heterogeneous lung phantom, the 2D gamma pass rate was 100% for 6 and 15 MV beams. The model optimization time was less than 4 min.

CONCLUSION

The proposed source model optimization process accurately produces photon fluence spectra from a linac using valid physical properties, without detailed knowledge of the geometry of the linac head, and with minimal optimization time.

A mathematical deconvolution formulation for superficial dose distribution measurement by Cerenkov light dosimetry

PURPOSE

Cerenkov photons are created by high-energy radiation beams used for radiation therapy. In this study, we developed a Cerenkov light dosimetry technique to obtain a two-dimensional dose distribution in a superficial region of medium from the images of Cerenkov photons by using a deconvolution method.

METHODS

An integral equation was derived to represent the Cerenkov photon image acquired by a camera for a given incident high-energy photon beam by using convolution kernels. Subsequently, an equation relating the planar dose at a depth to a Cerenkov photon image using the well-known relationship between the incident beam fluence and the dose distribution in a medium was obtained. The final equation contained a convolution kernel called the Cerenkov dose scatter function (CDSF). The CDSF function was obtained by deconvolving the Cerenkov scatter function (CSF) with the dose scatter function (DSF). The GAMOS (Geant4-based Architecture for Medicine-Oriented Simulations) Monte Carlo particle simulation software was used to obtain the CSF and DSF. The dose distribution was calculated from the Cerenkov photon intensity data using an iterative deconvolution method with the CDSF. The theoretical formulation was experimentally evaluated by using an optical phantom irradiated by high-energy photon beams.

RESULTS

The intensity of the deconvolved Cerenkov photon image showed linear dependence on the dose rate and the photon beam energy. The relative intensity showed a field size dependence similar to the beam output factor. Deconvolved Cerenkov images showed improvement in dose profiles compared with the raw image data. In particular, the deconvolution significantly improved the agreement in the high dose gradient region, such as in the penumbra. Deconvolution with a single iteration was found to provide the most accurate solution of the dose. Two-dimensional dose distributions of the deconvolved Cerenkov images agreed well with the reference distributions for both square fields and a multileaf collimator (MLC) defined, irregularly shaped field.

CONCLUSIONS

The proposed technique improved the accuracy of the Cerenkov photon dosimetry in the penumbra region. The results of this study showed initial validation of the deconvolution method for beam profile measurements in a homogeneous media. The new formulation accounted for the physical processes of Cerenkov photon transport in the medium more accurately than previously published methods.

Image guidance in proton therapy for lung cancer

Proton therapy is a promising but challenging treatment modality for the management of lung cancer. The technical challenges are due to respiratory motion, low dose tolerance of adjacent normal tissue and tissue density heterogeneity. Different imaging modalities are applied at various steps of lung proton therapy to provide information on target definition, target motion, proton range, patient setup and treatment outcome assessment. Imaging data is used to guide treatment design, treatment delivery, and treatment adaptation to ensure the treatment goal is achieved. This review article will summarize and compare various imaging techniques that can be used in every step of lung proton therapy to address these challenges.

Application of Geant4 Monte Carlo simulation in dose calculations for small radiosurgical fields

The Geant4 toolkit was used to develop a Monte Carlo (MC)-based engine for accurate dose calculations in small radiation field sizes. The Geant4 toolkit (version 10.1.p02) was used to simulate 6-MV photon beam of a Varian2100C linear accelerator that is being used for stereotactic radiosurgery (SRS) treatment with small radiation fields. Geometric models of 3 in-house designed radiosurgical divergent cones, with the diameters of their projections at the isocenter being 10, 20, and 30 mm, were simulated. The accuracy of the MC simulation technique was examined by reproducing several different simulated dosimetric parameters of the primary beams with the experimental data. The dose distributions are first checked for single beams for each cone, then standard multiple field (SMF) techniques are applied. A sample set of DICOM files from computed tomography (CT) scan imaging of a patient's head was converted to the Geant4 geometry format to implement MC-based engine for a clinical test. To validate the accuracy of the MC-based calculations for SMF arrangements, the isodose lines from MC simulation in water phantom were compared with the measured isodose lines using EBT3 Gafchromic film in Solid Water phantoms. Agreements between measured and simulated depth dose values and beam profiles for SRS cones were generally within 2%/2 mm. For output factors, the largest discrepancy was observed for 10 mm SRS cone, which was 1.7%. For SMF techniques, in SRS cones, the MC simulation and EBT3 Gafchromic film dosimetry were in acceptable agreement (5%/5 mm). Excellent agreement between the results of the MC-based and measured dose values for both single and SMF techniques in SRS cones indicates the ability of the Geant4 toolkit to be applied as the platform for treatment planning of advanced radiotherapy techniques.

Development of a flattening filter free multiple source model for use as an independent, Monte Carlo, dose calculation, quality assurance tool for clinical trials

PURPOSE

The Imaging and Radiation Oncology Core-Houston (IROC-H) Quality Assurance Center (formerly the Radiological Physics Center) has reported varying levels of compliance from their anthropomorphic phantom auditing program. IROC-H studies have suggested that one source of disagreement between institution submitted calculated doses and measurement is the accuracy of the institution's treatment planning system dose calculations and heterogeneity corrections used. In order to audit this step of the radiation therapy treatment process, an independent dose calculation tool is needed.

METHODS

Monte Carlo multiple source models for Varian flattening filter free (FFF) 6 MV and FFF 10 MV therapeutic x-ray beams were commissioned based on central axis depth dose data from a 10 × 10 cm² field size and dose profiles for a 40 × 40 cm² field size. The models were validated against open-field measurements in a water tank for field sizes ranging from 3 × 3 cm² to 40 × 40 cm² . The models were then benchmarked against IROC-H's anthropomorphic head and neck phantom and lung phantom measurements.

RESULTS

Validation results, assessed with a ±2%/2 mm gamma criterion, showed average agreement of 99.9% and 99.0% for central axis depth dose data for FFF 6 MV and FFF 10 MV models, respectively. Dose profile agreement using the same evaluation technique averaged 97.8% and 97.9% for the respective models. Phantom benchmarking comparisons were evaluated with a ±3%/2 mm gamma criterion, and agreement averaged 90.1% and 90.8% for the respective models.

CONCLUSIONS

Multiple source models for Varian FFF 6 MV and FFF 10 MV beams have been developed, validated, and benchmarked for inclusion in an independent dose calculation quality assurance tool for use in clinical trial audits.

Development of a Monte Carlo multiple source model for inclusion in a dose calculation auditing tool

PURPOSE

The Imaging and Radiation Oncology Core Houston (IROC-H) (formerly the Radiological Physics Center) has reported varying levels of agreement in their anthropomorphic phantom audits. There is reason to believe one source of error in this observed disagreement is the accuracy of the dose calculation algorithms and heterogeneity corrections used. To audit this component of the radiotherapy treatment process, an independent dose calculation tool is needed.

METHODS

Monte Carlo multiple source models for Elekta 6 MV and 10 MV therapeutic x-ray beams were commissioned based on measurement of central axis depth dose data for a 10 × 10 cm² field size and dose profiles for a 40 × 40 cm² field size. The models were validated against open field measurements consisting of depth dose data and dose profiles for field sizes ranging from 3 × 3 cm² to 30 × 30 cm² . The models were then benchmarked against measurements in IROC-H's anthropomorphic head and neck and lung phantoms.

RESULTS

Validation results showed 97.9% and 96.8% of depth dose data passed a ±2% Van Dyk criterion for 6 MV and 10 MV models respectively. Dose profile comparisons showed an average agreement using a ±2%/2 mm criterion of 98.0% and 99.0% for 6 MV and 10 MV models respectively. Phantom plan comparisons were evaluated using ±3%/2 mm gamma criterion, and averaged passing rates between Monte Carlo and measurements were 87.4% and 89.9% for 6 MV and 10 MV models respectively.

CONCLUSIONS

Accurate multiple source models for Elekta 6 MV and 10 MV x-ray beams have been developed for inclusion in an independent dose calculation tool for use in clinical trial audits.

Measurement of skin surface dose distributions in radiation therapy using poly(vinyl alcohol) cryogel dosimeters

In external beam radiation therapy (EBRT), skin dose measurement is important to evaluate dose coverage of superficial target volumes. Treatment planning systems (TPSs) are often inaccurate in this region of the patient, so in vivo measurements are necessary for skin surface dose estimation. In this work, superficial dose distributions were measured using radiochromic translucent poly(vinyl alcohol) cryogels. The cryogels simultaneously served as bolus material, providing the necessary buildup to achieve the desired superficial dose. The relationship between dose to the skin surface and dose measured with the bolus was established using a series of oblique irradiations with gantry angles ranging from 0° to 90°. EBT-2 Gafchromic film was placed under the bolus, and the ratio of bolus-film dose was determined ranging from 0.749 ± 0.005 to 0.930 ± 0.002 for 0° and 90° gantry angles, respectively. The average ratio over 0-67.5° (0.800 ± 0.064) was used as the single correction factor to convert dose in bolus to dose to the skin surface. The correction factor was applied to bolus measurements of skin dose from head and neck intensity-modulated radiation therapy (IMRT) treatments delivered to a RANDO phantom. The resulting dose distributions were compared to film measurements using gamma analysis with a 3%/3 mm tolerance and a 10% threshold. The minimum gamma pass rate was 95.2% suggesting that the radiochromic bolus may provide an accurate estimation of skin surface dose using a simple correction factor. This study demonstrates the suitability of radiochromic cryogels for superficial dose measurements in megavoltage photon beams.

Efficient independent planar dose calculation for FFF IMRT QA with a bivariate Gaussian source model

The aim of this study is to perform a direct comparison of the source model for photon beams with and without flattening filter (FF) and to develop an efficient independent algorithm for planar dose calculation for FF‐free (FFF) intensity‐modulated radiotherapy (IMRT) quality assurance (QA). The source model consisted of a point source modeling the primary photons and extrafocal bivariate Gaussian functions modeling the head scatter, monitor chamber backscatter, and collimator exchange effect. The model parameters were obtained by minimizing the difference between the calculated and measured in‐air output factors (S_c). The fluence of IMRT beams was calculated from the source model using a backprojection and integration method. The off‐axis ratio in FFF beams were modeled with a fourth degree polynomial. An analytical kernel consisting of the sum of three Gaussian functions was used to describe the dose deposition process. A convolution‐based method was used to account for the ionization chamber volume averaging effect when commissioning the algorithm. The algorithm was validated by comparing the calculated planar dose distributions of FFF head‐and‐neck IMRT plans with measurements performed with a 2D diode array. Good agreement between the measured and calculated S_c was achieved for both FF beams (<0.25%) and FFF beams (<0.10%). The relative contribution of the head‐scattered photons reduced by 34.7% for 6 MV and 49.3% for 10 MV due to the removal of the FF. Superior agreement between the calculated and measured dose distribution was also achieved for FFF IMRT. In the gamma comparison with a 2%/2 mm criterion, the average passing rate was 96.2 ± 1.9% for 6 MV FFF and 95.5 ± 2.6% for 10 MV FFF. The efficient independent planar dose calculation algorithm is easy to implement and can be valuable in FFF IMRT QA.

Modification and validation of an analytical source model for external beam radiotherapy Monte Carlo dose calculations

PURPOSE

A dose calculation tool, which combines the accuracy of the dose planning method (DPM) Monte Carlo code and the versatility of a practical analytical multisource model, which was previously reported has been improved and validated for the Varian 6 and 10 MV linear accelerators (linacs). The calculation tool can be used to calculate doses in advanced clinical application studies. One shortcoming of current clinical trials that report dose from patient plans is the lack of a standardized dose calculation methodology. Because commercial treatment planning systems (TPSs) have their own dose calculation algorithms and the clinical trial participant who uses these systems is responsible for commissioning the beam model, variation exists in the reported calculated dose distributions. Today's modern linac is manufactured to tight specifications so that variability within a linac model is quite low. The expectation is that a single dose calculation tool for a specific linac model can be used to accurately recalculate dose from patient plans that have been submitted to the clinical trial community from any institution. The calculation tool would provide for a more meaningful outcome analysis.

METHODS

The analytical source model was described by a primary point source, a secondary extra-focal source, and a contaminant electron source. Off-axis energy softening and fluence effects were also included. The additions of hyperbolic functions have been incorporated into the model to correct for the changes in output and in electron contamination with field size. A multileaf collimator (MLC) model is included to facilitate phantom and patient dose calculations. An offset to the MLC leaf positions was used to correct for the rudimentary assumed primary point source.

RESULTS

Dose calculations of the depth dose and profiles for field sizes 4 × 4 to 40 × 40 cm agree with measurement within 2% of the maximum dose or 2 mm distance to agreement (DTA) for 95% of the data points tested. The model was capable of predicting the depth of the maximum dose within 1 mm. Anthropomorphic phantom benchmark testing of modulated and patterned MLCs treatment plans showed agreement to measurement within 3% in target regions using thermoluminescent dosimeters (TLD). Using radiochromic film normalized to TLD, a gamma criteria of 3% of maximum dose and 2 mm DTA was applied with a pass rate of least 85% in the high dose, high gradient, and low dose regions. Finally, recalculations of patient plans using DPM showed good agreement relative to a commercial TPS when comparing dose volume histograms and 2D dose distributions.

CONCLUSIONS

A unique analytical source model coupled to the dose planning method Monte Carlo dose calculation code has been modified and validated using basic beam data and anthropomorphic phantom measurement. While this tool can be applied in general use for a particular linac model, specifically it was developed to provide a singular methodology to independently assess treatment plan dose distributions from those clinical institutions participating in National Cancer Institute trials.

Optical cone beam tomography of Cherenkov-mediated signals for fast 3D dosimetry of x-ray photon beams in water

PURPOSE

To test the use of a three-dimensional (3D) optical cone beam computed tomography reconstruction algorithm, for estimation of the imparted 3D dose distribution from megavoltage photon beams in a water tank for quality assurance, by imaging the induced Cherenkov-excited fluorescence (CEF).

METHODS

An intensified charge-coupled device coupled to a standard nontelecentric camera lens was used to tomographically acquire two-dimensional (2D) projection images of CEF from a complex multileaf collimator (MLC) shaped 6 MV linear accelerator x-ray photon beam operating at a dose rate of 600 MU/min. The resulting projections were used to reconstruct the 3D CEF light distribution, a potential surrogate of imparted dose, using a Feldkamp-Davis-Kress cone beam back reconstruction algorithm. Finally, the reconstructed light distributions were compared to the expected dose values from one-dimensional diode scans, 2D film measurements, and the 3D distribution generated from the clinical Varian ECLIPSE treatment planning system using a gamma index analysis. A Monte Carlo derived correction was applied to the Cherenkov reconstructions to account for beam hardening artifacts.

RESULTS

3D light volumes were successfully reconstructed over a 400 × 400 × 350 mm(3) volume at a resolution of 1 mm. The Cherenkov reconstructions showed agreement with all comparative methods and were also able to recover both inter- and intra-MLC leaf leakage. Based upon a 3%/3 mm criterion, the experimental Cherenkov light measurements showed an 83%-99% pass fraction depending on the chosen threshold dose.

CONCLUSIONS

The results from this study demonstrate the use of optical cone beam computed tomography using CEF for the profiling of the imparted dose distribution from large area megavoltage photon beams in water.

Development and reproducibility evaluation of a Monte Carlo-based standard LINAC model for quality assurance of multi-institutional clinical trials

Technical developments in radiotherapy (RT) have created a need for systematic quality assurance (QA) to ensure that clinical institutions deliver prescribed radiation doses consistent with the requirements of clinical protocols. For QA, an ideal dose verification system should be independent of the treatment-planning system (TPS). This paper describes the development and reproducibility evaluation of a Monte Carlo (MC)-based standard LINAC model as a preliminary requirement for independent verification of dose distributions. The BEAMnrc MC code is used for characterization of the 6-, 10- and 15-MV photon beams for a wide range of field sizes. The modeling of the LINAC head components is based on the specifications provided by the manufacturer. MC dose distributions are tuned to match Varian Golden Beam Data (GBD). For reproducibility evaluation, calculated beam data is compared with beam data measured at individual institutions. For all energies and field sizes, the MC and GBD agreed to within 1.0% for percentage depth doses (PDDs), 1.5% for beam profiles and 1.2% for total scatter factors (Scps.). Reproducibility evaluation showed that the maximum average local differences were 1.3% and 2.5% for PDDs and beam profiles, respectively. MC and institutions' mean Scps agreed to within 2.0%. An MC-based standard LINAC model developed to independently verify dose distributions for QA of multi-institutional clinical trials and routine clinical practice has proven to be highly accurate and reproducible and can thus help ensure that prescribed doses delivered are consistent with the requirements of clinical protocols.

The clinical impact of detector choice for beam scanning

Recently, the developers of Eclipse have recommended the use of ionization chambers for all profile scanning, including for the modeling of VMAT and stereotactic applications. The purpose of this study is to show the clinical impact caused by the choice of detector with respect to its ability to accurately measure dose in the penumbra and tail regions of a scanned profile. Using scan data acquired with several detectors, including an IBA CC13, a PTW 60012, and a Sun Nuclear EDGE Detector, three complete beam models are created, one for each respective detector. Next, using each beam model, dose volumes are retrospectively recalculated from actual anonymous patient plans. These plans include three full‐arc VMAT prostate plans, three left chest wall plans delivered using irregular compensators, two half‐arc VMAT lung plans, three MLC‐collimated static‐field pairs, and two SBRT liver plans. Finally, plans are reweighted to deliver the same number of monitor units, and mean dose‐to‐target volumes and organs at risk are calculated and compared. Penumbra width did not play a role. Dose in the tail region of the profile made the largest difference. By overresponding in the tail region of the profile, the 60012 diode detector scan data affected the beam model in such a way that target doses were reduced by as much as 0.4% (in comparison to CC13 and EDGE data). This overresponse also resulted in an overestimation of dose to peripheral critical structure, whose dose consisted mainly of scatter. This study shows that, for modeling the 6 MV beam of Acuros XB in Eclipse Version 11, the choice to use a CC13 scanning ion chamber or an EDGE Detector was an unimportant choice, providing nearly identical models in the treatment planning system.

PACS number: 87.55.kh

Use of IAEA's phase-space files for the implementation of a clinical accelerator virtual source model

In the present work, phase-space data files (phsp) provided by the International Atomic Energy Agency (IAEA) for different accelerators were used in order to develop a Virtual Source Model (VSM) for clinical photon beams. Spectral energy distributions extracted from supplied phsp files were used to define the radiation pattern of a virtual extended source in a hybrid model which is completed with a virtual diaphragm used to simulate both electron contamination and the shape of the penumbra region. This simple virtual model was used as the radiation source for dosimetry calculations in a water phantom. The proposed model proved easy to build and test, and good agreement with clinical accelerators dosimetry measurements were obtained for different field sizes. Our results suggest this simple method could be useful for treatment planning systems (TPS) verification purposes.

Independent calculation-based verification of IMRT plans using a 3D dose-calculation engine

Independent monitor unit verification of intensity-modulated radiation therapy (IMRT) plans requires detailed 3-dimensional (3D) dose verification. The aim of this study was to investigate using a 3D dose engine in a second commercial treatment planning system (TPS) for this task, facilitated by in-house software. Our department has XiO and Pinnacle TPSs, both with IMRT planning capability and modeled for an Elekta-Synergy 6MV photon beam. These systems allow the transfer of computed tomography (CT) data and RT structures between them but do not allow IMRT plans to be transferred. To provide this connectivity, an in-house computer programme was developed to convert radiation therapy prescription (RTP) files as generated by many planning systems into either XiO or Pinnacle IMRT file formats. Utilization of the technique and software was assessed by transferring 14 IMRT plans from XiO and Pinnacle onto the other system and performing 3D dose verification. The accuracy of the conversion process was checked by comparing the 3D dose matrices and dose volume histograms (DVHs) of structures for the recalculated plan on the same system. The developed software successfully transferred IMRT plans generated by 1 planning system into the other. Comparison of planning target volume (TV) DVHs for the original and recalculated plans showed good agreement; a maximum difference of 2% in mean dose, - 2.5% in D95, and 2.9% in V95 was observed. Similarly, a DVH comparison of organs at risk showed a maximum difference of +7.7% between the original and recalculated plans for structures in both high- and medium-dose regions. However, for structures in low-dose regions (less than 15% of prescription dose) a difference in mean dose up to +21.1% was observed between XiO and Pinnacle calculations. A dose matrix comparison of original and recalculated plans in XiO and Pinnacle TPSs was performed using gamma analysis with 3%/3mm criteria. The mean and standard deviation of pixels passing gamma tolerance for XiO-generated IMRT plans was 96.1 ± 1.3, 96.6 ± 1.2, and 96.0 ± 1.5 in axial, coronal, and sagittal planes respectively. Corresponding results for Pinnacle-generated IMRT plans were 97.1 ± 1.5, 96.4 ± 1.2, and 96.5 ± 1.3 in axial, coronal, and sagittal planes respectively.

Monte Carlo simulation for Neptun 10 PC medical linear accelerator and calculations of output factor for electron beam

AIM

Exact knowledge of dosimetric parameters is an essential pre-requisite of an effective treatment in radiotherapy. In order to fulfill this consideration, different techniques have been used, one of which is Monte Carlo simulation.

MATERIALS AND METHODS

This study used the MCNP-4C to simulate electron beams from Neptun 10 PC medical linear accelerator. Output factors for 6, 8 and 10 MeV electrons applied to eleven different conventional fields were both measured and calculated.

RESULTS

The measurements were carried out by a Wellhofler-Scanditronix dose scanning system. Our findings revealed that output factors acquired by MCNP-4C simulation and the corresponding values obtained by direct measurements are in a very good agreement.

CONCLUSION

In general, very good consistency of simulated and measured results is a good proof that the goal of this work has been accomplished.

GPU-accelerated Monte Carlo convolution/superposition implementation for dose calculation

PURPOSE

Dose calculation is a key component in radiation treatment planning systems. Its performance and accuracy are crucial to the quality of treatment plans as emerging advanced radiation therapy technologies are exerting ever tighter constraints on dose calculation. A common practice is to choose either a deterministic method such as the convolution/superposition (CS) method for speed or a Monte Carlo (MC) method for accuracy. The goal of this work is to boost the performance of a hybrid Monte Carlo convolution/superposition (MCCS) method by devising a graphics processing unit (GPU) implementation so as to make the method practical for day-to-day usage.

METHODS

Although the MCCS algorithm combines the merits of MC fluence generation and CS fluence transport, it is still not fast enough to be used as a day-to-day planning tool. To alleviate the speed issue of MC algorithms, the authors adopted MCCS as their target method and implemented a GPU-based version. In order to fully utilize the GPU computing power, the MCCS algorithm is modified to match the GPU hardware architecture. The performance of the authors' GPU-based implementation on an Nvidia GTX260 card is compared to a multithreaded software implementation on a quad-core system.

RESULTS

A speedup in the range of 6.7-11.4x is observed for the clinical cases used. The less than 2% statistical fluctuation also indicates that the accuracy of the authors' GPU-based implementation is in good agreement with the results from the quad-core CPU implementation.

CONCLUSIONS

This work shows that GPU is a feasible and cost-efficient solution compared to other alternatives such as using cluster machines or field-programmable gate arrays for satisfying the increasing demands on computation speed and accuracy of dose calculation. But there are also inherent limitations of using GPU for accelerating MC-type applications, which are also analyzed in detail in this article.

Evaluation of calculation methods of collimator scatter factors in a linear accelerator equipped with MLC instead of lower collimators

In the monitor unit verification for high-energy radiation therapy, we evaluated methods of calculation of collimator scatter factors (S(c)) in a linear accelerator equipped with MLC instead of lower collimators. Routinely,S(c) is calculated from rectangular fields shaped by upper and lower jaws in the linear accelerator. However, this calculation method should not be used for the linear accelerator equipped with MLC instead of lower collimators. Consequently, we used a backprojected field at the flattening filter plane projected by calculation point's eye view on each MLC. We then attempted to deviseS(c) by using Clarkson's integration for these backprojected irregular fields. This method makes it possible to calculate collimator scatter factors in error of less than +/-0.3% in all of sixteen measured irregular fields.

In vivo verification of superficial dose for head and neck treatments using intensity-modulated techniques

Skin dose is one of the key issues for clinical dosimetry in radiation therapy. Currently planning computer systems are unable to accurately predict dose in the buildup region, leaving ambiguity as to the dose levels actually received by the patient's skin during radiotherapy. This is one of the prime reasons why in vivo measurements are necessary to estimate the dose in the buildup region. A newly developed metal-oxide-semiconductor-field-effect-transistor (MOSFET) detector designed specifically for dose measurements in rapidly changing dose gradients was introduced for accurate in vivo skin dosimetry. The feasibility of this detector for skin dose measurements was verified in comparison with plane parallel ionization chamber and radiochromic films. The accuracy of a commercial treatment planning system (TPS) in skin dose calculations for intensity-modulated radiation therapy treatment of nasopharyngeal carcinoma was evaluated using MOSFET detectors in an anthropomorphic phantom as well as on the patients. Results show that this newly developed MOSFET detector can provide a minimal but highly reproducible intrinsic buildup of 7 mg cm(-2) corresponding to the requirements of personal surface dose equivalent Hp (0.07). The reproducibility of the MOSFET response, in high sensitivity mode, is found to be better than 2% at the phantom surface for the doses normally delivered to the patients. The MOSFET detector agrees well with the Attix chamber and the EBT Gafchromic film in terms of surface and buildup region dose measurements, even for oblique incident beams. While the dose difference between MOSFET measurements and TPS calculations is within measurement uncertainty for the depths equal to or greater than 0.5 cm, an overestimation of up to 8.5% was found for the surface dose calculations in the anthropomorphic phantom study. In vivo skin dose measurements reveal that the dose difference between the MOSFET results and the TPS calculations was on average -7.2%, ranging from -4.3% to -9.2%. The newly designed MOSFET detector encapsulated into a thin water protective film has a minimal reproducible intrinsic buildup recommended for skin dosimetry. This feature makes it very suitable for routine IMRT QA and accurate in vivo skin dosimetry.

VMC++ versus BEAMnrc: a comparison of simulated linear accelerator heads for photon beams

BEAMnrc, a code for simulating medical linear accelerators based on EGSnrc, has been bench-marked and used extensively in the scientific literature and is therefore often considered to be the gold standard for Monte Carlo simulations for radiotherapy applications. However, its long computation times make it too slow for the clinical routine and often even for research purposes without a large investment in computing resources. VMC++ is a much faster code thanks to the intensive use of variance reduction techniques and a much faster implementation of the condensed history technique for charged particle transport. A research version of this code is also capable of simulating the full head of linear accelerators operated in photon mode (excluding multileaf collimators, hard and dynamic wedges). In this work, a validation of the full head simulation at 6 and 18 MV is performed, simulating with VMC++ and BEAMnrc the addition of one head component at a time and comparing the resulting phase space files. For the comparison, photon and electron fluence, photon energy fluence, mean energy, and photon spectra are considered. The largest absolute differences are found in the energy fluences. For all the simulations of the different head components, a very good agreement (differences in energy fluences between VMC++ and BEAMnrc <1%) is obtained. Only a particular case at 6 MV shows a somewhat larger energy fluence difference of 1.4%. Dosimetrically, these phase space differences imply an agreement between both codes at the <1% level, making VMC++ head module suitable for full head simulations with considerable gain in efficiency and without loss of accuracy.

A Monte-Carlo derived dual-source model for helical tomotherapy treatment planning

Full Monte Carlo radiation transport simulations of accelerator heads are impractical for routine treatment planning because of the excessive computational burden and memory requirements. To improve the accuracy and efficiency of treatment plans for helical tomotherapy, we have developed a dual-source model to characterize the radiation emitted from the head of a commercial helical tomotherapy accelerator. Percentage depth dose (PDD) and beam profiles computed using the dual-source model with the EGS/BEAMnrc Monte Carlo package agree within 2% of measurements for a wide range of field sizes, which suggests that the proposed dual-source model provides an adequate representation of the tomotherapy head for dose calculations in routine treatment planning.

Monte carlo evaluation of the AAA treatment planning algorithm in a heterogeneous multilayer phantom and IMRT clinical treatments for an Elekta SL25 linear accelerator

The Anisotropic Analytical Algorithm (AAA) is a new pencil beam convolution/superposition algorithm proposed by Varian for photon dose calculations. The configuration of AAA depends on linear accelerator design and specifications. The purpose of this study was to investigate the accuracy of AAA for an Elekta SL25 linear accelerator for small fields and intensity modulated radiation therapy (IMRT) treatments in inhomogeneous media. The accuracy of AAA was evaluated in two studies. First, AAA was compared both with Monte Carlo (MC) and the measurements in an inhomogeneous phantom simulating lung equivalent tissues and bone ribs. The algorithm was tested under lateral electronic disequilibrium conditions, using small fields (2 x 2 cm(2)). Good agreement was generally achieved for depth dose and profiles, with deviations generally below 3% in lung inhomogeneities and below 5% at interfaces. However, the effects of attenuation and scattering close to the bone ribs were not fully taken into account by AAA, and small inhomogeneities may lead to planning errors. Second, AAA and MC were compared for IMRT plans in clinical conditions, i.e., dose calculations in a computed tomography scan of a patient. One ethmoid tumor, one orophaxynx and two lung tumors are presented in this paper. Small differences were found between the dose volume histograms. For instance, a 1.7% difference for the mean planning target volume dose was obtained for the ethmoid case. Since better agreement was achieved for the same plans but in homogeneous conditions, these differences must be attributed to the handling of inhomogeneities by AAA. Therefore, inherent assumptions of the algorithm, principally the assumption of independent depth and lateral directions in the scaling of the kernels, were slightly influencing AAA's validity in inhomogeneities. However, AAA showed a good accuracy overall and a great ability to handle small fields in inhomogeneous media compared to other pencil beam convolution algorithms.

A virtual photon source model of an Elekta linear accelerator with integrated mini MLC for Monte Carlo based IMRT dose calculation

A dedicated, efficient Monte Carlo (MC) accelerator head model for intensity modulated stereotactic radiosurgery treatment planning is needed to afford a highly accurate simulation of tiny IMRT fields. A virtual source model (VSM) of a mini multi-leaf collimator (MLC) (the Elekta Beam Modulator (EBM)) is presented, allowing efficient generation of particles even for small fields. The VSM of the EBM is based on a previously published virtual photon energy fluence model (VEF) (Fippel et al 2003 Med. Phys. 30 301) commissioned with large field measurements in air and in water. The original commissioning procedure of the VEF, based on large field measurements only, leads to inaccuracies for small fields. In order to improve the VSM, it was necessary to change the VEF model by developing (1) a method to determine the primary photon source diameter, relevant for output factor calculations, (2) a model of the influence of the flattening filter on the secondary photon spectrum and (3) a more realistic primary photon spectrum. The VSM model is used to generate the source phase space data above the mini-MLC. Later the particles are transmitted through the mini-MLC by a passive filter function which significantly speeds up the time of generation of the phase space data after the mini-MLC, used for calculation of the dose distribution in the patient. The improved VSM model was commissioned for 6 and 15 MV beams. The results of MC simulation are in very good agreement with measurements. Less than 2% of local difference between the MC simulation and the diamond detector measurement of the output factors in water was achieved. The X, Y and Z profiles measured in water with an ion chamber (V = 0.125 cm(3)) and a diamond detector were used to validate the models. An overall agreement of 2%/2 mm for high dose regions and 3%/2 mm in low dose regions between measurement and MC simulation for field sizes from 0.8 x 0.8 cm(2) to 16 x 21 cm(2) was achieved. An IMRT plan film verification was performed for two cases: 6 MV head&neck and 15 MV prostate. The simulation is in agreement with film measurements within 2%/2 mm in the high dose regions (> or = 0.1 Gy = 5% D(max)) and 5%/2 mm in low dose regions (<0.1 Gy).

An integrated Monte Carlo dosimetric verification system for radiotherapy treatment planning

An integrated Monte Carlo (MC) dose calculation system, MCRTV (Monte Carlo for radiotherapy treatment plan verification), has been developed for clinical treatment plan verification, especially for routine quality assurance (QA) of intensity-modulated radiotherapy (IMRT) plans. The MCRTV system consists of the EGS4/PRESTA MC codes originally written for particle transport through the accelerator, the multileaf collimator (MLC), and the patient/phantom, which run on a 28-CPU Linux cluster, and the associated software developed for the clinical implementation. MCRTV has an interface with a commercial treatment planning system (TPS) (Eclipse, Varian Medical Systems, Palo Alto, CA, USA) and reads the information needed for MC computation transferred in DICOM-RT format. The key features of MCRTV have been presented in detail in this paper. The phase-space data of our 15 MV photon beam from a Varian Clinac 2300C/D have been developed and several benchmarks have been performed under homogeneous and several inhomogeneous conditions (including water, aluminium, lung and bone media). The MC results agreed with the ionization chamber measurements to within 1% and 2% for homogeneous and inhomogeneous conditions, respectively. The MC calculation for a clinical prostate IMRT treatment plan validated the implementation of the beams and the patient/phantom configuration in MCRTV.