Absolute dose calculations for Monte Carlo simulations of radiotherapy beams

Effect of modulation factor and low dose threshold level on gamma pass rates of single isocenter multi-target SRT treatment plans

PURPOSE

SRS MapCHECK (SMC) is a commercially available patient-specific quality assurance (PSQA) tool for stereotactic radiosurgery (SRS) applications. This study investigates the effects of degree of modulation, location off-axis, and low dose threshold (LDT) selection on gamma pass rates (GPRs) between SMC and treatment planning system, Analytical Anisotropic Algorithm (AAA), or Vancouver Island Monte Carlo (VMC++ algorithm) system calculated dose distributions.

METHODS

Volumetric-modulated arc therapy (VMAT) plans with modulation factors (MFs) ranging from 2.7 to 10.2 MU/cGy were delivered to SMC at isocenter and 6 cm off-axis. SMC measured dose distributions were compared against AAA and VMC++ via gamma analysis (3%/1 mm) with LDT of 10% to 80% using SNC Patient software.

RESULTS

Comparing on-axis SMC dose against AAA and VMC++ with LDT of 10%, all AAA-calculated plans met the acceptance criteria of GPR ≥ 90%, and only one VMC++ calculated plan was marginally outside the acceptance criteria with pass rate of 89.1%. Using LDT of 80% revealed decreasing GPR with increasing MF. For AAA, GPRs reduced from 100% at MF of 2.7 MU/cGy to 57% at MF of 10.2 MU/cGy, and for VMC++ calculated plans, the GPRs reduced from 89% to 60% in the same MF range. Comparison of SMC dose off-axis against AAA and VMC++ showed more pronounced reduction of GPR with increasing MF. For LDT of 10%, AAA GPRs reduced from 100% to 83% in the MF range of 2.7 to 9.8 MU/cGy, and VMC++ GPR reduced from 100% to 91% in the same range. With 80% LDT, GPRs dropped from 100% to 42% for both algorithms.

CONCLUSIONS

MF, dose calculation algorithm, and LDT selections are vital in VMAT-based SRT PSQA. LDT of 80% enhances sensitivity of gamma analysis for detecting dose differences compared to 10% LDT. To achieve better agreement between calculated and SMC dose, it is recommended to limit the MF to 4.6 MU/cGy on-axis and 3.6 MU/cGy off-axis.

The delivery assessment for small targets on Halcyon radiotherapy system - Measured and calculated dose comparison

BACKGROUND

With the ever-increasing requirements of accuracy and personalization of radiotherapy treatments, stereotactic radiotherapy (SRT) with volumetric modulated arc therapy (VMAT) on O-ring Halcyon radiotherapy system could potentially provide a fast, safe, and feasible treatment option.

PURPOSE

The purpose of this study was to assess the delivery of Halcyon VMAT plans for small targets.

METHODS

Well-defined VMAT-SRT plans were created on Halcyon radiotherapy system with the stacked and staggered dual-layer MLC design for the film measurement set-up and the target sizes and shapes designed to emulate the targets of the stereotactic treatments. The planar dose distributions were acquired with film measurements and compared to a current clinical reference dose calculation with AcurosXB (v18.0, Varian Medical Systems) and to Monte Carlo simulations. With the collapsed arc versions of the VMAT-SRT plans, the uncertainty in dose delivery due to the multileaf collimator (MLC) without the gantry rotation could be separated and analyzed.

RESULTS

The target size was mainly limited by the resolution originated from the design of the MLC leaves. The results of the collapsed arc versions of the plans show good consistency among measured, calculated, and simulated dose distributions. With the full VMAT plans, the agreement between calculated and simulated dose distributions was consistent with the collapsed arc versions. The measured dose distribution agreed with the calculated and simulated dose distributions within the target regions, but considerable local differences were observed in the margins of the target. The largest differences located in the steep gradient regions presumably originating from the deviation of the isocenter.

CONCLUSIONS

The potential of the Halcyon radiotherapy system for VMAT-SRT delivery was evaluated and the study revealed valuable insights on the machine characteristics with the delivery.

Assessment of four dose calculation algorithms using IAEA-TECDOC-1583 with medium dependency correction factor (Kmed) application

PURPOSE

This study discusses the measurement of dose in clinical commissioning tests described in IAEA-TECDOC-1583. It explores the application of Monte Carlo (MC) modelled medium dependency correction factors (Kmed) for accurate dose measurement in bone and lung materials using the CIRS phantom.

METHODS

BEAMnrc codes simulate radiation sources and model radiation transport for 6 MV and 15 MV photon beams. CT images of the CIRS phantom are converted to an MC compatible phantom. The PTW 30013 farmer chamber measures doses within modeled CIRS phantom. Kmed are determined by averaging values from four central voxels within the sensitive volume of the farmer chamber. Kmed is calculated for Dm.m and Dw.w algorithm types in bone and lung media for both photon beams.

RESULTS

Average modelled correction factors for Dm.m calculations using the farmer chamber are 0.976 (±0.1 %) for 6 MV and 0.979 (±0.1 %) for 15 MV in bone media. Correspondingly, correction factors for Dw.w calculations are 0.99 (±0.3 %) and 0.992 (±0.4 %), respectively. For lung media, average correction factors for Dm.m calculations are 1.02 (±0.3 %) for 6 MV and 1.022 (±0.4 %) for 15 MV. Correspondingly, correction factors for Dw.w calculations are 1.01 (±0.3 %) and 1.012 (±0.2 %), respectively.

CONCLUSIONS

This study highlights the significant impact of applying Kmed on dose differences between measurement and calculation during the dose audit process.

Dosimetric validation of SmART-RAD Monte Carlo modelling for x-ray cabinet radiobiology irradiators

Objective. Accuracy and reproducibility in the measurement of radiation dose and associated reporting are critically important for the validity of basic and preclinical radiobiological studies performed with kilovolt x-ray radiation cabinets. This is essential to enable results of radiobiological studies to be repeated, as well as enable valid comparisons between laboratories. In addition, the commonly used single point dose value hides the 3D dose heterogeneity across the irradiated sample. This is particularly true for preclinical rodent models, and is generally difficult to measure directly. Radiation transport simulations integrated in an easy to use application could help researchers improve quality of dosimetry and reporting.Approach. This paper describes the use and dosimetric validation of a newly-developed Monte Carlo (MC) tool, SmART-RAD, to simulate the x-ray field in a range of standard commercial x-ray cabinet irradiators used for preclinical irradiations. Comparisons are made between simulated and experimentally determined dose distributions for a range of configurations to assess the potential use of this tool in determining dose distributions through samples, based on more readily available air-kerma calibration point measurements.Main results. Simulations gave very good dosimetric agreement with measured depth dose distributions in phantoms containing both water and bone equivalent materials. Good spatial and dosimetric agreement between simulated and measured dose distributions was obtained when using beam-shaping shielding.Significance. The MC simulations provided by SmART-RAD provide a useful tool to go from a limited number of dosimetry measurements to detailed 3D dose distributions through a non-homogeneous irradiated sample. This is particularly important when trying to determine the dose distribution in more complex geometries. The use of such a tool can improve reproducibility and dosimetry reporting in preclinical radiobiological research.

Organ absorbed doses in the IORT treatment of breast cancer with the INTRABEAM device: a Monte-Carlo study

Purpose: The current prescription and the assessment of the delivered absorbed dose in intraoperative radiation therapy (IORT) with the INTRABEAM system rely mainly on depth-dose measurements in water. The accuracy of this approach is limited because tissue heterogeneity is ignored. It is also difficult to accurately determine the dose delivered to the patient experimentally as the steep dose gradient is highly sensitive to geometric errors. Our goal is to determine the dose to the target volume and the organs at risk of a clinical breast cancer patient from treatment with the system.Methods: A homogeneous water-equivalent CT dataset was derived from the preoperative CT scan of a patient by setting all materials in the patient volume as water-equivalent. This homogeneous CT data represents the current assumption of a homogenous patient, while the original CT data is considered the ground truth. An in-house Monte Carlo algorithm was used to simulate the delivered dose in both setups for a prescribed treatment dose of 20 Gy to the surface of the 3.5 cm diameter spherical applicator.Results: The doses received by 2% (D2%) of the target volume for the homogeneous and heterogeneous geometries are 16.26 Gy and 9.33 Gy, respectively. The D2% for the heart are 0.035 Gy and 0.119 Gy for the homogeneous and heterogeneous geometries, respectively. This trend is also observed for the other organs at risk.Conclusions: The assumption of a homogeneous patient overestimates the dose to the target volume and underestimates the doses to the organs at risk.

Monte Carlo modeling of Halcyon and Ethos radiotherapy beam using CAD geometry: validation and IAEA-compliant phase space

Objective.A Monte Carlo (MC) model of a Halcyon and Ethos (Varian Medical Systems, a Siemens Healthineers Company) radiotherapy beam was validated and field-independent phase space (PHSP) files were recorded above the dual-layer multileaf collimators (MLC).Approach.The treatment head geometry was modeled according to engineering drawings and the dual-layer MLC was imported from CAD (computer-aided design) files. The information for the incident electron beam was achieved from an iterative electromagnetic solver. The validation of the model was performed by comparing the dose delivered by the square MLC fields as well as complex field measurements.Main results.An electron phase space was generated from linac simulations and achieved improved MC results. The output factors for square fields were within 1% and the largest differences of 5% were found in the build-up region of PDDs and the penumbra region of profiles. With the more complicated MLC-shaped field (Fishbone), the largest differences of up to 8% were found in the MLC leaf tip region due to the uncertainty of the MLC positioning and the mechanical leaf gap value. The impact of the collimator rotation on the PHSP solution has been assessed with both small and large fields, confirming negligible effects on in-field and out-of-field dose distributions.Significance.A computational model of the Halcyon and Ethos radiotherapy beam with a high accuracy implementation of the MLC was shown to be able to reproduce the radiation beam characteristics with square fields and more complex MLC-shaped fields. The field-independent PHSP files that were produced can be used as an accurate treatment head model above the MLC, and reduce the time to simulate particle transport through treatment head components.

Absorbed dose calculation for a realistic CT-derived mouse phantom irradiated with a standard Cs-137 cell irradiator using a Monte Carlo method

Computed tomography (CT) derived Monte Carlo (MC) phantoms allow dose determination within small animal models that is not feasible with in-vivo dosimetry. The aim of this study was to develop a CT-derived MC phantom generated from a mouse with a xenograft tumour that could then be used to calculate both the dose heterogeneity in the tumour volume and out of field scattered dose for pre-clinical small animal irradiation experiments. A BEAMnrc Monte-Carlo model has been built of our irradiation system that comprises a lead collimator with a 1 cm diameter aperture fitted to a Cs-137 gamma irradiator. The MC model of the irradiation system was validated by comparing the calculated dose results with dosimetric film measurement in a polymethyl methacrylate (PMMA) phantom using a 1D gamma-index analysis. Dose distributions in the MC mouse phantom were calculated and visualized on the CT-image data. Dose volume histograms (DVHs) were generated for the tumour and organs at risk (OARs). The effect of the xenographic tumour volume on the scattered out of field dose was also investigated. The defined gamma index analysis criteria were met, indicating that our MC simulation is a valid model for MC mouse phantom dose calculations. MC dose calculations showed a maximum out of field dose to the mouse of 7% of Dmax. Absorbed dose to the tumour varies in the range 60%-100% of Dmax. DVH analysis demonstrated that tumour received an inhomogeneous dose of 12 Gy-20 Gy (for 20 Gy prescribed dose) while out of field doses to all OARs were minimized (1.29 Gy-1.38 Gy). Variation of the xenographic tumour volume exhibited no significant effect on the out of field scattered dose to OARs. The CT derived MC mouse model presented here is a useful tool for tumour dose verifications as well as investigating the doses to normal tissue (in out of field) for preclinical radiobiological research.

Dosimetric accuracy of Acuros® XB and AAA algorithms for stereotactic body radiotherapy (SBRT) lung treatments: evaluation with PRIMO Monte Carlo code

Purpose:

The study aimed to compare the dosimetric performance of Acuros® XB (AXB) and anisotropic analytical algorithm (AAA) for lung SBRT plans using Monte Carlo (MC) simulations.

Methods:

We compared the dose calculation algorithms AAA and either of the dose reporting modes of AXB (dose to medium (AXB-Dm) or dose to water (AXB-Dw)) algorithms implemented in Eclipse® (Varian Medical Systems, Palo Alto, CA) Treatment planning system (TPS) with MC. PRIMO code was used for the MC simulations. The TPS-calculated dose profiles obtained with a multi-slab heterogeneity phantom were compared to MC. A lung phantom with a tumour was used to validate TPS algorithms using different beam delivery techniques. 2D gamma values obtained from Gafchromic film measurements in the tumour isocentre plane were compared with TPS algorithms and MC. Ten VMAT SBRT plans generated in TPS with each algorithm were recalculated with a PRIMO MC system for identical beam parameters for the clinical plan validation. A dose–volume histogram (DVH) based plan comparison and a 3D global gamma analysis were performed.

Results:

AXB demonstrated better agreement with MC and film measurements in the lung phantom validation, with good agreement in PDD, profiles and gamma analysis. AAA showed an overestimated PDD, a significant difference in dose profiles and a lower gamma pass rate near the field borders. With AAA, there was a dose overestimation at the periphery of the tumour. For clinical plan validation, AXB demonstrated higher agreement with MC than AAA.

Conclusions:

AXB provided better agreement with MC than AAA in the phantom and clinical plan evaluations.

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 new DOSXYZnrc method for Monte Carlo simulations of 4D dose distributions

The purpose of this study is to present a novel method for generating Monte Carlo 4D dose distributions in a single DOSXYZnrc simulation. During a standard simulation, individual energy deposition events are summed up to generate a 3D dose distribution and their associated temporal information is discarded. This means that in order to determine dose distributions as a function of time, separate simulations would have to be run for each interval of interest. Consequently, it has not been clinically feasible until now to routinely perform Monte Carlo simulations of dose rate, time-resolved dose accumulation, or electronic portal imaging devices (EPID) cine-mode images for volumetric modulated arc therapy (VMAT) plans. To overcome this limitation, we modified DOSXYZnrc and defined new input and output variables that allow a time-like parameter associated with each particle history to be binned in a user-defined manner. Under the new code version, computation times are the same as for a standard simulation, and the time-integrated 4D dose is identical to the standard 3D dose. We present a comparison of scintillator measurements and Monte Carlo simulations for dose rate during a VMAT beam delivery, a study of dose rate in a VMAT total body irradiation plan, and simulations of transit (through-patient) EPID cine-mode images.

Clinical verification of treatment planning dose calculation in lung SBRT with GATE Monte Carlo simulation code

PURPOSE

This study aims to use GATE/Geant4 simulation code to evaluate the performance of dose calculations with Anisotropic Analytical Algorithm (AAA) in the context of lung SBRT for complex treatments considering images of patients.

METHODS

Four cases of non-small cell lung cancer treated with SBRT were selected for this study. Irradiation plans were created with AAA and recalculated end to end using Monte Carlo (MC) method maintaining field configurations identical to the original plans. Each treatment plan was evaluated in terms of PTV and organs at risk (OARs) using dose-volume histograms (DVH). Dosimetric parameters obtained from DVHs were used to compare AAA and MC.

RESULTS

The comparison between the AAA and MC DVH using gamma analysis with the passing criteria of 3%/3% showed an average passing rate of more than 90% for the PTV structure and 97% for the OARs. Tightening the criteria to 2%/2% showed a reduction in the average passing rate of the PTV to 86%. The agreement between the AAA and MC dose calculations for PTV dosimetric parameters (V₁₀₀; V₉₀; Homogeneity index; maximum, minimum and mean dose; CI_Paddick and D_2cm) was within 18.4%. For OARs, the biggest differences were observed in the spinal cord and the great vessels.

CONCLUSIONS

In general, we did not find significant differences between AAA and MC. The results indicate that AAA could be used in complex SBRT cases that involve a larger number of small treatment fields in the presence of tissue heterogeneities.

PENELOPE simulations and experiment for 6 MV clinac iX accelerator for standard and small static fields

The goal of this work was to produce accurate data for use as a 'gold standard' and a valid tool for measurements in reference dosimetry for standard/small static field sizes from 0.5 × 0.5 to 10 × 10 cm². It is based on the accuracy of the phase space files (PSFs) as a key quantity. Because the IAEA general public database provides few PSFs for the Varian iX, we simulated the head through Monte Carlo (MC) simulations and calculated validated PSFs for 12 square field sizes including seven for small static fields. The resulting dosimetric calculations allowed us to reach a good level of agreement in comparison to our relative and absolute dose measurements performed on a Varian iX in water phantom. Measured and MC calculated output factors were investigated for different detectors. Based on the TRS 483 formalism and MC (PENELOPE/penEasy), we calculated output correction factors for the unshielded Diode-E (T60017) and the PinPoint-3D (T31016) micro-chamber according to manufacturers' blueprints. Our MC results were in agreement with the recommended data; they compete with recent measurements and MC simulations and in particular the TRS 483 MC data obtained from similar simulations. Moreover, our MC results provide supplemental data in comparison to TRS 483 data in particular for the PinPoint-3D (T31016). We suggest our MC output correction factors as new datasets for future TRS compilations. The work was substantial, used different robust MC strategies depending on the scoring regions, and led in most cases to uncertainties of less than 1%.

Monte Carlo calculated detector-specific correction factors for Elekta radiosurgery cones

A radiation field is considered small if its dimension is lower than the range of secondary electrons and the collimating devices partially occlude the source. Different detector types, such as unshielded diodes, diamond detectors, and small-volume ion chambers, are used for small-field measurements. Although the active volumes of these detectors are small, their non-water equivalent materials cause response variations. Herein, we aim to calculate the correction factors for our clinical detectors, EDGE detector (Sun Nuclear), 60017 diode (PTW), and CC01 ion chamber (IBA), for stereotactic radiosurgery cones of diameters of 5-15 mm in an Elekta Synergy linear accelerator using a Monte Carlo simulation. An Elekta Synergy linear accelerator treatment head was simulated using BEAMnrc Monte Carlo code as per the manufacturer specification. All three detectors were simulated as per the manufacturer specification. Three EGSnrc user codes were used for the detector simulation based on the detector geometry. The Monte Carlo model of the treatment head was validated against the measured data for a standard field size of 10 × 10 cm². The off-axis profile, percentage depth dose, and tissue phantom ratioTPR1020were verified in the validation procedure. The measured and Monte Carlo calculated relative output factors (ROFs) were not consistent. In a 5 mm field size, EDGE diode overestimated the ROF by 7.06%, and 60017 diode to 4.611%. In a 7.5 mm field size, the variations were 4.295% and 3.691% for EDGE and 60017 diodes, respectively. CC01 ion chamber under-responded up to 10% because of its low-density active volume. The maximum corrections were obtained in the smallest field size, which were 0.939(0.007), 0.962(0.006), and 1.117(0.008) for EDGE, PTW T60017, and CC01 detectors, respectively. After applying the Monte Carlo calculated correction factor to the measured ROF, it became consistent with the Monte Carlo calculated ROF.

Monte Carlo simulation of 6-MV dynamic wave VMAT deliveries by Vero4DRT linear accelerator using EGSnrc moving sources

The commissioning and benchmark of a Monte Carlo (MC) model of the 6‐MV Brainlab‐Mitsubishi Vero4DRT linear accelerator for the purpose of quality assurance of clinical dynamic wave arc (DWA) treatment plans is reported. Open‐source MC applications based on EGSnrc particle transport codes are used to simulate the medical linear accelerator head components. Complex radiotherapy irradiations can be simulated in a single MC run using a shared library format combined with BEAMnrc “source20.” Electron energy tuning is achieved by comparing measured vs simulated percentage depth doses (PDDs) for MLC‐defined field sizes in a water phantom. Electron spot size tuning is achieved by comparing measured and simulated inplane and crossplane beam profiles. DWA treatment plans generated from RayStation (RaySearch) treatment planning system (TPS) are simulated on voxelized (2.5 mm³) patient CT datasets. Planning target volume (PTV) and organs at risk (OAR) dose–volume histograms (DVHs) are compared to TPS‐calculated doses for clinically deliverable dynamic volumetric modulated arc therapy (VMAT) trajectories. MC simulations with an electron beam energy of 5.9 MeV and spot size FWHM of 1.9 mm had the closest agreement with measurement. DWA beam deliveries simulated on patient CT datasets results in DVH agreement with TPS‐calculated doses. PTV coverage agreed within 0.1% and OAR max doses (to 0.035 cc volume) agreed within 1 Gy. This MC model can be used as an independent dose calculation from the TPS and as a quality assurance tool for complex, dynamic radiotherapy treatment deliveries. Full patient CT treatment simulations are performed in a single Monte Carlo run in 23 min. Simulations are run in parallel using the Condor High‐Throughput Computing software¹ on a cluster of eight servers. Each server has two physical processors (Intel Xeon CPU E5‐2650 0 @2.00 GHz), with 8 cores per CPU and two threads per core for 256 calculation nodes.

Monte Carlo calculations of radiotherapy dose in "homogeneous" anatomy

Given the substantial literature on the use of Monte Carlo (MC) simulations to verify treatment planning system (TPS) calculations of radiotherapy dose in heterogeneous regions, such as head and neck and lung, this study investigated the potential value of running MC simulations of radiotherapy treatments of nominally homogeneous pelvic anatomy. A pre-existing in-house MC job submission and analysis system, built around BEAMnrc and DOSXYZnrc, was used to evaluate the dosimetric accuracy of a sample of 12 pelvic volumetric arc therapy (VMAT) treatments, planned using the Varian Eclipse TPS, where dose was calculated with both the Analytical Anisotropic Algorithm (AAA) and the Acuros (AXB) algorithm. In-house TADA (Treatment And Dose Assessor) software was used to evaluate treatment plan complexity, in terms of the small aperture score (SAS), modulation index (MI) and a novel exposed leaf score (ELS/ELA). Results showed that the TPS generally achieved closer agreement with the MC dose distribution when treatments were planned for smaller (single-organ) targets rather than larger targets that included nodes or metastases. Analysis of these MC results with reference to the complexity metrics indicated that while AXB was useful for reducing dosimetric uncertainties associated with density heterogeneity, the residual TPS dose calculation uncertainties resulted from treatment plan complexity and TPS model simplicity. The results of this study demonstrate the value of using MC methods to recalculate and check the dose calculations provided by commercial radiotherapy TPSs, even when the treated anatomy is assumed to be comparatively homogeneous, such as in the pelvic region.

Clinical implementation of a Monte Carlo based independent TPS dose checking system

The increase in complexity of treatment plans over time through modalities such as intensity modulated radiotherapy (IMRT) and volumetric modulated arc therapy (VMAT) has often not been met with an increase in capability of the secondary dose calculation checking systems typically used to verify the treatment planning system. Monte Carlo (MC) codes such as EGSnrc have become easily available and are capable of performing calculations of highly complex radiotherapy treatments. This educational note demonstrates a method for implementing and using a fully automated system for performing and analysing full MC calculations of conformal, IMRT and VMAT radiotherapy plans. Example calculations were based on BEAMnrc/DOSXYZnrc and are performed automatically after either uploading exported plan DICOM data through a Python-based web interface, or exporting DICOM data to a monitored network location. This note demonstrates how completed MC calculations can then be analysed using an automatically generated dose point comparison report, or easily re-imported back into the treatment planning system. Agreement between the TPS and MC calculation was an improvement on agreement between RadCalc and the TPS, with differences ranging from 1.2 to 5.5% between RadCalc and the treatment planning system (TPS), and 0.1-1.7% between MC and TPS. Comparison of the dose-volume histogram (DVH) parameters [Formula: see text], [Formula: see text], [Formula: see text], and [Formula: see text] for the example VMAT plans showed agreement for the mean planning target volume dose within [Formula: see text], [Formula: see text] and [Formula: see text] generally within [Formula: see text] with the exception of a brain case, and [Formula: see text] within [Formula: see text]. Overall, this note provides a demonstration of a system that has been integrated well into existing clinical workflow, and has been shown to be a valuable additional tool in the secondary checking of treatment plan calculations.

External beam patient dose verification based on the integral quality monitor (IQM®) output signals

BACKGROUND

The Integral Quality Monitor (IQM^®) can essentially measure the integral fluence through a segment and provide real-time information about the accuracy of radiation delivery based on comparisons of measured segment signals and pre-calculated reference values. However, the present IQM chamber cannot calculate the dose in the patient.

AIM

This study aims to make use of IQM field output signals to calculate the number of monitor units (MUs) delivered through an arbitrary treatment field in order to convert Monte Carlo (MC)-generated dose distributions in a patient model into absolute dose.

METHODS

XiO and Monaco treatment planning systems (TPSs) were used to define treatment beam portals for cervix and esophagus conformal radiotherapy as well as prostate intensity-modulated radiotherapy for the translation of patient and beam setup information from DICOM to DOSXYZnrc. The planned beams were simulated in a patient model built from actual patient CT images and each simulated integral field/segment was weighted with its MUs before summation to get the total dose in the plan. The segment beam weights (MUs) were calculated as the ratio of the open-field IQM measured signal and the calculated signal per MU extracted from chamber sensitivity maps. These are the actual MUs delivered not just MUs set. The beam weighting method was evaluated by comparing weighted MC doses with original planned doses using profile and isodose comparisons, dose difference maps, γ analysis and dose-volume histogram (DVH) data.

RESULTS

γ pass rates of up to 98% were found, except for the esophagus plan where the γ pass rate was below 45%. DVH comparisons showed good agreement for most organs, with the largest differences observed in low-density lung. However, these discrepancies can result from differences in dose calculation algorithms or differences in MUs used for dose weighting planned by the TPS and MUs calculated using IQM field output signals. To test this, a 4-field box DOSXYZnrc MC simulation weighted with planned (XiO) MUs was compared with the same simulation weighted with IQM-based MUs. Dose differences of up to 5% were found on the isocentre slice. For XiO versus MC, up to 7% dose differences were found, indicating additional error due to limitations of XiO's superposition algorithm. Dose differences between MC Monaco and MC EGSnrc were less than 3%.

CONCLUSIONS

The most valuable comparison was MC versus MC as it eliminated algorithm discrepancies and evaluated dose differences precisely according to beam weighting. For XiO TPS, care must be taken as dose differences may also arise due to limitations in XiO's planning software, not merely due to differences in MUs. Overall, the IQM was successfully used to compute beam dose weights to accurately reconstruct the patient dose using unweighted MC beams. Our technique can be used for pre-treatment QA provided each segment output is known and an accurate linac source model is available.

New capabilities of the Monte Carlo dose engine ARCHER-RT: Clinical validation of the Varian TrueBeam machine for VMAT external beam radiotherapy

PURPOSE

The Monte Carlo radiation transport method is considered the most accurate approach for absorbed dose calculations in external beam radiation therapy. In this study, an efficient and accurate source model of the Varian TrueBeam 6X STx Linac is developed and integrated with a fast Monte Carlo photon-electron transport absorbed dose engine, ARCHER-RT, which is capable of being executed on CPUs, NVIDIA GPUs, and AMD GPUs. This capability of fast yet accurate radiation dose calculation is essential for clinical utility of this new technology. This paper describes the software and algorithmic developments made to the ARCHER-RT absorbed dose engine.

METHODS

AMD's Heterogeneous-Compute Interface for Portability (HIP) was implemented in ARCHER-RT to allow for device independent execution on NVIDIA and AMD GPUs. Architecture-specific atomic-add algorithms have been identified and both more accurate single-precision and double-precision computational absorbed dose calculation methods have been added to ARCHER-RT and validated through a test case to evaluate the accuracy and performance of the algorithms. The validity of the source model and the radiation transport physics were benchmarked against Monte Carlo simulations performed with EGSnrc. Secondary dose-check physics plans, and a clinical prostate treatment plan were calculated to demonstrate the applicability of the platform for clinical use. Absorbed dose difference maps and gamma analyses were conducted to establish the accuracy and consistency between the two Monte Carlo models. Timing studies were conducted on a CPU, an NVIDIA GPU, and an AMD GPU to evaluate the computational speed of ARCHER-RT.

RESULTS

Percent depth doses were computed for different field sizes ranging from 1.5 cm²  × 1.5 cm² to 22 cm²  × 40cm² and the two codes agreed for all points outside high gradient regions within 3%. Axial profiles computed for a 10 cm²  × 10 cm² field for multiple depths agreed for all points outside high gradient regions within 2%. The test case investigating the impact of native single-precision compared to double-precision showed differences in voxels as large as 71.47% and the implementation of KAS single-precision reduced the difference to less than 0.01%. The 3%/3mm gamma pass rates for an MPPG5a multileaf collimator (MLC) test case and a clinical VMAT prostate plan were 94.2% and 98.4% respectively. Timing studies demonstrated the calculation of a VMAT plan was completed in 50.3, 187.9, and 216.8 s on an NVIDIA GPU, AMD GPU, and Intel CPU, respectively.

CONCLUSION

ARCHER-RT is capable of patient-specific VMAT external beam photon absorbed dose calculations and its potential has been demonstrated by benchmarking against a well validated EGSnrc model of a Varian TrueBeam. Additionally, the implementation of AMD's HIP has shown the flexibility of the ARCHER-RT platform for device independent calculations. This work demonstrates the significant addition of functionality added to ARCHER-RT framework which has marked utility for both research and clinical applications and demonstrates further that Monte Carlo-based absorbed dose engines like ARCHER-RT have the potential for widespread clinical implementation.

Monte Carlo simulation using PRIMO code as a tool for checking the credibility of commissioning and quality assurance of 6 MV TrueBeam STx varian LINAC

Aim

To validate and implement Monte Carlo simulation using PRIMO code as a tool for checking the credibility of measurements in LINAC initial commissioning and routine Quality Assurance (QA). Relative and absolute doses of 6 MV photon beam from TrueBeam STx Varian Linear Accelerator (LINAC) were simulated and validated with experimental measurement, Analytical Anisotropic Algorithm (AAA) calculation, and golden beam.

Methods and Materials

Varian phase-space files were imported to the PRIMO code and four blocks of jaws were simulated to determine the field size of the photon beam. Water phantom was modeled in the PRIMO code with water equivalent density. Golden beam data, experimental measurement, and AAA calculation results were imported to PRIMO code for gamma comparison.

Results

PRIMO simulations of Percentage Depth Dose (PDD) and in-plane beam profiles had good agreement with experimental measurements, AAA calculations and golden beam. However, PRIMO simulations of cross-plane beam profiles have a better agreement with AAA calculation and golden beam than the experimental measurement. Furthermore, PRIMO simulations of absolute dose agreed well with experimental results with ±0.8% uncertainty.

Conclusion

The PRIMO code has good accuracy and is appropriate for use as a tool to check the credibility of beam scanning and output measurement in initial commissioning and routine QA.

Technical Note: A dose calculation framework for dynamic electron arc radiotherapy (DEAR) using VirtuaLinac Monte Carlo simulation tool

PURPOSE

Dynamic electron arc radiotherapy (DEAR) is a novel dynamic technique that achieves highly conformal dose through simultaneous couch and gantry motion during delivery. The purpose of this study is to develop a framework integrating a Monte Carlo dose engine (VirtuaLinac) to a treatment planning system (TPS, Eclipse) for DEAR. A quality assurance (QA) procedure is also developed.

METHODS AND MATERIALS

The interfaces include the following: computed tomography image export and conversion for VirtuaLinac; VirtuaLinac computation tasks management through application programming interface (API); and dose matrix processing and evaluation. The framework was validated with both static beam and DEAR plan with a 3 × 3 cm² cutout for both 6 and 9 MeV electrons. Verification plans for DEAR were created on flat phantom and a hybrid dose calculation technique was developed which convolves precalculated small field kernel with the beam trajectory, and the resulting dose was compared with the full VirtuaLinac calculation and film measurement.

RESULTS

Excellent agreement between VirtuaLinac and eMC was observed with three-dimensional γ pass rate of 98% at 1%/1 mm criteria for both 6 and 9 MeV electrons. Film measurement shows two-dimensional (2D) γ passing rate of 99.8 % (6 MeV) and 97.1% (9 MeV) at 2%/2 mm criteria. For DEAR plans the comparison of VirtuaLinac and measurement shows the 2D γ passing rates of 94% at 2%/2 mm for 6 MeV. The dose distributions from hybrid method in phantom are identical to the full VirtuaLinac simulations, but can be done instantly.

CONCLUSIONS

A framework has been developed for DEAR dose calculation using VirtuaLinac Monte Carlo dose engine. The VirtuaLinac calculated dose was validated against measurement. A feasible and practical DEAR QA method has been developed for dose measurement in phantom. The hybrid dose calculation technique is efficient and suitable for DEAR QA purpose.

Electron modulated arc therapy (EMAT) using photon MLC for postmastectomy chest wall treatment I: Monte Carlo-based dosimetric characterizations

PURPOSE

To study the dosimetric properties of electron arc beams delivered by photon-beam multi-leaf collimators (pMLC) in electron modulated arc therapy (EMAT) for postmastectomy chest wall treatments.

METHODS

Using the Monte Carlo method, we simulated a 2100EX Varian linear accelerator and verified the beam models in a water tank. Dosimetric characterizations were performed on cylindrical water phantoms of elliptical bases with various field sizes, arc ranges and source-to-surface distances (SSDs) for 6, 9 and 12 MeV beam energy.

RESULTS

The arc beam has a higher bremsstrahlung dose than the static beam at the isocenter due to crossfire, but choosing a field size greater than 5 cm effectively reduces the bremsstrahlung dose. The depths of the 90% maximum dose located at 1.7, 2.8 and 4.1 cm for 6, 9 and 12 MeV, respectively, are similar to those of the static beams and independent of the field size and arc range.

CONCLUSION

Based on the study, we recommend using the 5 cm field width for electron arc beams considering both bremsstrahlung dose at the isocenter and the arc profile penumbra. To ensure sufficient PTV edge coverage, we recommend a field length extension of at least 4 cm from PTV's edge for all beam energies and an arc extension of around 7°, 5°, and 5° for beam energies 6, 9, and 12 MeV, respectively. These dosimetric characterizations are the basis of pMLC-delivered EMAT treatment planning for postmastectomy chest wall patients.

Dosimetric and Monte Carlo verification of jaws-only IMRT plans calculated by the Collapsed Cone Convolution algorithm for head and neck cancers

AIM

The aim of this study is to verify the Prowess Panther jaws-only intensity modulated radiation therapy (JO-IMRT) treatment planning (TP) by comparing the TP dose distributions for head-and-neck (H&N) cancer with the ones simulated by Monte Carlo (MC).

BACKGROUND

To date, dose distributions planned using JO-IMRT for H&N patients were found superior to the corresponding three-dimensional conformal radiotherapy (3D-CRT) plans. Dosimetry of the JO-IMRT plans were also experimentally verified using an ionization chamber, MapCHECK 2, and Octavius 4D and good agreements were shown.

MATERIALS AND METHODS

Dose distributions of 15 JO-IMRT plans of nasopharyngeal patients were recalculated using the EGSnrc Monte Carlo code. The clinical photon beams were simulated using the BEAMnrc. The absorbed dose to patients treated by fixed-field IMRT was computed using the DOSXYZnrc. The simulated dose distributions were then compared with the ones calculated by the Collapsed Cone Convolution (CCC) algorithm on the TPS, using the relative dose error comparison and the gamma index using global methods implemented in PTW-VeriSoft with 3%/3 mm, 2%/2 mm, 1%/1 mm criteria.

RESULTS

There is a good agreement between the MC and TPS dose. The average gamma passing rates were 93.3 ± 3.1%, 92.8 ± 3.2%, 92.4 ± 3.4% based on the 3%/3 mm, 2%/2 mm, 1%/1 mm criteria, respectively.

CONCLUSIONS

According to the results, it is concluded that the CCC algorithm was adequate for most of the IMRT H&N cases where the target was not immediately adjacent to the critical structures.

Dosimetric verification of lung phantom calculated by collapsed cone convolution: A Monte Carlo and experimental evaluation

OBJECTIVE

To evaluate the dose calculation accuracy in the Prowess Panther treatment planning system (TPS) using the collapsed cone convolution (CCC) algorithm.

METHODS

The BEAMnrc Monte Carlo (MC) package was used to predict the dose distribution of photon beams produced by the Oncor® linear accelerator (linac). The MC model of an 18 MV photon beam was verified by measurement using a p-type diode dosimeter. Percent depth dose (PDD) and dose profiles were used for comparison based on three field sizes: 5×5, 10×10, and 20×20cm2. The accuracy of the CCC dosimetry was also evaluated using a plan composed of a simple parallel-opposed field (11×16cm2) in a lung phantom comprised of four tissue simulating media namely, lung, soft tissue, bone and spinal cord. The CCC dose calculation accuracy was evaluated by MC simulation and measurements according to the dose difference and 3D gamma analysis. Gamma analysis was carried out through comparison of the Monte Carlo simulation and the TPS calculated dose.

RESULTS

Compared to the dosimetric results measured by the Farmer chamber, the CCC algorithm underestimated dose in the planning target volume (PTV), right lung and lung-tissue interface regions by about -0.11%, -1.6 %, and -2.9%, respectively. Moreover, the CCC algorithm underestimated the dose at the PTV, right lung and lung-tissue interface regions in the order of -0.34%, -0.4% and -3.5%, respectively, when compared to the MC simulation. Gamma analysis results showed that the passing rates within the PTV and heterogeneous region were above 59% and 76%. For the right lung and spinal cord, the passing rates were above 80% for all gamma criteria.

CONCLUSIONS

This study demonstrates that the CCC algorithm has potential to calculate dose with sufficient accuracy for 3D conformal radiotherapy within the thorax where a significant amount of tissue heterogeneity exists.

The potential impact of ultrathin filter design on dosimetry and relative biological effectiveness in modern image-guided small animal irradiators

OBJECTIVE:

Modern image-guided small animal irradiators like the Xstrahl Small Animal Radiation Research Platform (SARRP) are designed with ultrathin 0.15 mm Cu filters, which compared with more heavily filtrated traditional cabinet-style biological irradiators, produce X-ray spectra weighted toward lower energies, impacting the dosimetric properties and the relative biological effectiveness (RBE). This study quantifies the effect of ultrathin filter design on relative depth dose profiles, absolute dose output, and RBE using Monte Carlo techniques.

METHODS:

The percent depth-dose and absolute dose output are calculated using kVDoseCalc and EGSnrc, respectively, while a tally based on the induction of double-strand breaks as a function of electron spectra invoked in PENELOPE is used to estimate the RBE.

RESULTS:

The RBE increases by >2.4% in the ultrathin filter design compared to a traditional irradiator. Furthermore, minute variations in filter thickness have notable effects on the dosimetric properties of the X-ray beam, increasing the percent depth dose (at 2 cm in water) by + 0.4%/0.01 mm Cu and decreasing absolute dose (at 2 cm depth in water) by -1.8%/0.01 mm Cu for the SARRP.

CONCLUSIONS:

These results show that modern image-guided irradiators are quite sensitive to small manufacturing variations in filter thickness, and show a small change in RBE compared to traditional X-ray irradiators.

ADVANCES IN KNOWLEDGE:

We quantify the consequences of ultrathin filter design in modern image-guided biological irradiators on relative and absolute dose, and RBE. Our results show these to be small, but not insignificant, suggesting laboratories transitioning between irradiators should carefully design their radiobiological experiments.

A fast jaw-tracking model for VMAT and IMRT Monte Carlo simulations

Modern radiotherapy techniques involve routine use of volumetric arc therapy (VMAT) and intensity modulated radiotherapy (IMRT) with jaw‐tracking – dynamic motion of the secondary collimators (jaws) in tandem with multi‐leaf collimators (MLCs). These modalities require accurate dose calculations for the purposes of treatment planning and dose verification. Monte Carlo (MC) methods for radiotherapy dose calculation are widely accepted as capable of achieving high accuracy. This paper presents an efficiency‐enhancement method for secondary collimator modeling, presented in the context of a tool for MC‐based dose second checks. The model constitutes an accuracy trade‐off in the source model for the sake of efficiency enhancement, but maintains the advantages of MC transport in patient heterogeneities. The secondary collimator model is called Flat‐Absorbing‐Jaw‐Tracking (FAJT). Transmission through and scatter from the secondary collimators is neglected, and jaws are modeled as perfectly absorbing planes. To couple the motion of secondary collimators with MLCs for jaw‐tracking, the FAJT model was built into the VCU‐MLC model. Gamma‐index analysis of the dose distributions from FAJT against the full BEAMnrc MC simulations showed over 99% pass rate for a range of open fields, two clinical IMRT, and one VMAT treatment plan, for 2%/2 mm criteria above 10%. Using FAJT, the simulation speed of the secondary collimators for open fields increased by a factor of 237, 1489, and 1395 for 4 × 4, 10 × 10, and 30 × 30 cm², respectively. In general, clinically oriented simulation times are reduced from “hours” to “minutes” on identical hardware. Results for nine representative clinical cases (seven with jaw‐tracking) are presented. The average 2%/2 mm γ‐test success rate above the 80% isodose was 96.8% when tested against the EPIDose electronic portal image‐based dose reconstruction method and 97.3% against the Eclipse analytical anisotropic algorithm.

Measured and Monte Carlo simulated electron backscatter to the monitor chamber for the Varian TrueBeam Linac

To accurately simulate therapeutic electron beams using Monte Carlo methods, backscatter from jaws into the monitor chamber must be accounted for via the backscatter factor, S _b. Measured and simulated values of S _b for the TrueBeam are investigated. Two approaches for measuring S _b are presented. Both require service mode operation with the dose and pulse forming networking servos turned off in order to assess changes in dose rate with field size. The first approach samples an instantaneous dose rate, while the second approach times the delivery of a fixed number of monitor units to assess dose rate. Dose rates were measured for 6, 12 and 20 MeV electrons for jaw- or MLC-shaped apertures between [Formula: see text] and [Formula: see text] cm². The measurement techniques resulted in values of S _b that agreed within 0.21% for square and asymmetric fields collimated by the jaws. Measured values of S _b were used to calculate the forward dose component in a virtual monitor chamber using BEAMnrc. Based on this forward component, simulated values of S _b were calculated and compared to measurement and Varian's VirtuaLinac simulations. BEAMnrc results for jaw-shaped fields agreed with measurements and with VirtuaLinac simulations within 0.2%. For MLC-shaped fields, the respective measurement techniques differed by as much as 0.41% and BEAMnrc results differed with measurement by as much as 0.4%, however, all measured and simulated values agreed within experimental uncertainty. Measurement sensitivity was not sufficient to capture the small backscatter effect due to the MLC, and Monte Carlo predicted backscatter from the MLC to be no more than 0.3%. Backscatter from the jaws changed the electron dose rate by up to 2.6%. This reinforces the importance of including a backscatter factor in simulations of electron fields shaped with secondary collimating jaws, but presents the option of ignoring it when jaws are retracted and collimation is done with the MLC.

Full Monte Carlo and measurement-based overall performance assessment of improved clinical implementation of eMC algorithm with emphasis on lower energy range

New version 13.6.23 of the electron Monte Carlo (eMC) algorithm in Varian Eclipse™ treatment planning system has a model for 4MeV electron beam and some general improvements for dose calculation. This study provides the first overall accuracy assessment of this algorithm against full Monte Carlo (MC) simulations for electron beams from 4MeV to 16MeV with most emphasis on the lower energy range. Beams in a homogeneous water phantom and clinical treatment plans were investigated including measurements in the water phantom. Two different material sets were used with full MC: (1) the one applied in the eMC algorithm and (2) the one included in the Eclipse™ for other algorithms. The results of clinical treatment plans were also compared to those of the older eMC version 11.0.31. In the water phantom the dose differences against the full MC were mostly less than 3% with distance-to-agreement (DTA) values within 2mm. Larger discrepancies were obtained in build-up regions, at depths near the maximum electron ranges and with small apertures. For the clinical treatment plans the overall dose differences were mostly within 3% or 2mm with the first material set. Larger differences were observed for a large 4MeV beam entering curved patient surface with extended SSD and also in regions of large dose gradients. Still the DTA values were within 3mm. The discrepancies between the eMC and the full MC were generally larger for the second material set. The version 11.0.31 performed always inferiorly, when compared to the 13.6.23.

Monte Carlo simulation of beam characteristics from small fields based on TrueBeam flattening-filter-free mode

Purpose

Through the Monte Carlo (MC) simulation of 6 and 10 MV flattening-filter-free (FFF) beams from Varian TrueBeam accelerator, this study aims to find the best incident electron distribution for further studying the small field characteristics of these beams.

Methods

By incorporating the training materials of Varian on the geometry and material parameters of TrueBeam Linac head, the 6 and 10 MV FFF beams were modelled using the BEAMnrc and DOSXYZnrc codes, where the percentage depth doses (PDDs) and the off-axis ratios (OARs) curves of fields ranging from 4 × 4 to 40 × 40 cm² were simulated for both energies by adjusting the incident beam energy, radial intensity distribution and angular spread, respectively. The beam quality and relative output factor (ROF) were calculated. The simulations and measurements were compared using Gamma analysis method provided by Verisoft program (PTW, Freiburg, Germany), based on which the optimal MC model input parameters were selected and were further used to investigate the beam characteristics of small fields.

Results

The Full Width Half Maximum (FWHM), mono-energetic energy and angular spread of the resultant incident Gaussian radial intensity electron distribution were 0.75 mm, 6.1 MeV and 0.9° for the nominal 6 MV FFF beam, and 0.7 mm, 10.8 MeV and 0.3° for the nominal 10 MV FFF beam respectively. The simulation was mostly comparable to the measurement. Gamma criteria of 1 mm/1 % (local dose) can be met by all PDDs of fields larger than 1 × 1 cm², and by all OARs of no larger than 20 × 20 cm², otherwise criteria of 1 mm/2 % can be fulfilled. Our MC simulated ROFs agreed well with the measured ROFs of various field sizes (the discrepancies were less than 1 %), except for the 1 × 1 cm² field.

Conclusions

The MC simulation agrees well with the measurement and the proposed model parameters can be clinically used for further dosimetric studies of 6 and 10 MV FFF beams.

Quantification of dose differences between two versions of Acuros XB algorithm compared to Monte Carlo simulations--the effect on clinical patient treatment planning

A commercialized implementation of linear Boltzmann transport equation solver, the Acuros XB algorithm (AXB), represents a class of most advanced type 'c' photon radiotherapy dose calculation algorithms. The purpose of the study was to quantify the effects of the modifications implemented in the more recent version 11 of the AXB (AXB11) compared to the first commercial implementation, version 10 of the AXB (AXB10), in various anatomical regions in clinical treatment planning. Both versions of the AXB were part of Varian's Eclipse clinical treatment planning system and treatment plans for 10 patients were created using intensity-modulated radiotherapy (IMRT) and volumetric-modulated arc radiotherapy (VMAT). The plans were first created with the AXB10 and then recalculated with the AXB11 and full Monte Carlo (MC) simulations. Considering the full MC simulations as reference, a DVH analysis for gross tumor and planning target volumes (GTV and PTV) and organs at risk was performed, and also 3D gamma agreement index (GAI) values within a 15% isodose region and for the PTV were determined. Although differences up to 12% in DVH analysis were seen between the MC simulations and the AXB, based on the results of this study no general conclusion can be drawn that the modifications made in the AXB11 compared to the AXB10 would imply that the dose calculation accuracy of the AXB10 would be inferior to the AXB11 in the clinical patient treatment planning. The only clear improvement with the AXB11 over the AXB10 is the dose calculation accuracy in air cavities. In general, no large deviations are present in the DVH analysis results between the two versions of the algorithm, and the results of 3D gamma analysis do not favor one or the other. Thus it may be concluded that the results of the comprehensive studies assessing the accuracy of the AXB10 may be extended to the AXB11.

Jaw position uncertainty and adjacent fields in breast cancer radiotherapy

Locoregional treatment of breast cancer involves adjacent, half blocked fields matched at isocenter. The objective of this work is to study the dosimetric effects of the uncertainties in jaw positioning for such a case, and how a treatment planning protocol including adjacent field overlap of 1 mm affects the dose distribution. A representative treatment plan, involving 6 and 15 photon beams, for a patient treated at our hospital is chosen. Monte Carlo method (EGSnrc/BEAMnrc) is used to simulate the treatment. Uncertainties in jaw positioning of ± 1 mm are addressed, which implies extremes in reality of 2 mm field gap/overlap when planning adjacent fields without overlap and 1 mm gap or 3 mm overlap for a planning protocol with 1 mm overlap. Dosimetric parameters for PTV, lung and body are analyzed. Treatment planning protocol with 1 mm overlap of the adjacent fields does not considerably counteract possible underdosage of the target in the case studied. PTV-V95% is for example reduced from 95% for perfectly aligned fields to 90% and 91% for 2 mm and 1 mm gap, respectively. However, the risk of overdosage in PTV and in healthy soft tissue is increased when following the protocol with 1 mm overlap. A 3 mm overlap compared to 2 mm overlap results in an increase in maximum dose to PTV, PTV-D2%, from 113% to 121%. V120% for 'Body-PTV' is also increased from 5 cm(3) to 14 cm(3). A treatment planning protocol with 1 mm overlap does not considerably improve the coverage of PTV in the case of erroneous jaw positions causing gap between fields, but increases the overdosage in PTV and doses to healthy tissue, in the case of overlapping fields, for the case investigated.

A geometrical model for the Monte Carlo simulation of the TrueBeam linac

Monte Carlo simulation of linear accelerators (linacs) depends on the accurate geometrical description of the linac head. The geometry of the Varian TrueBeam linac is not available to researchers. Instead, the company distributes phase-space files of the flattening-filter-free (FFF) beams tallied at a plane located just upstream of the jaws. Yet, Monte Carlo simulations based on third-party tallied phase spaces are subject to limitations. In this work, an experimentally based geometry developed for the simulation of the FFF beams of the Varian TrueBeam linac is presented. The Monte Carlo geometrical model of the TrueBeam linac uses information provided by Varian that reveals large similarities between the TrueBeam machine and the Clinac 2100 downstream of the jaws. Thus, the upper part of the TrueBeam linac was modeled by introducing modifications to the Varian Clinac 2100 linac geometry. The most important of these modifications is the replacement of the standard flattening filters by ad hoc thin filters. These filters were modeled by comparing dose measurements and simulations. The experimental dose profiles for the 6 MV and 10 MV FFF beams were obtained from the Varian Golden Data Set and from in-house measurements performed with a diode detector for radiation fields ranging from 3  ×  3 to 40  ×  40 cm(2) at depths of maximum dose of 5 and 10 cm. Indicators of agreement between the experimental data and the simulation results obtained with the proposed geometrical model were the dose differences, the root-mean-square error and the gamma index. The same comparisons were performed for dose profiles obtained from Monte Carlo simulations using the phase-space files distributed by Varian for the TrueBeam linac as the sources of particles. Results of comparisons show a good agreement of the dose for the ansatz geometry similar to that obtained for the simulations with the TrueBeam phase-space files for all fields and depths considered, except for the 40  ×  40 cm(2) field where the ansatz geometry was able to reproduce the measured dose more accurately. Our approach overcomes some of the limitations of using the Varian phase-space files. It makes it possible to: (i) adapt the initial beam parameters to match measured dose profiles; (ii) reduce the statistical uncertainty to arbitrarily low values; and (iii) assess systematic uncertainties (type B) by using different Monte Carlo codes. One limitation of using phase-space files that is retained in our model is the impossibility of performing accurate absolute dosimetry simulations because the geometrical description of the TrueBeam ionization chamber remains unknown.

A hybrid phase-space and histogram source model for GPU-based Monte Carlo radiotherapy dose calculation

In GPU-based Monte Carlo simulations for radiotherapy dose calculation, source modelling from a phase-space source can be an efficiency bottleneck. Previously, this has been addressed using phase-space-let (PSL) sources, which provided significant efficiency enhancement. We propose that additional speed-up can be achieved through the use of a hybrid primary photon point source model combined with a secondary PSL source. A novel phase-space derived and histogram-based implementation of this model has been integrated into gDPM v3.0. Additionally, a simple method for approximately deriving target photon source characteristics from a phase-space that does not contain inheritable particle history variables (LATCH) has been demonstrated to succeed in selecting over 99% of the true target photons with only ~0.3% contamination (for a Varian 21EX 18 MV machine). The hybrid source model was tested using an array of open fields for various Varian 21EX and TrueBeam energies, and all cases achieved greater than 97% chi-test agreement (the mean was 99%) above the 2% isodose with 1% / 1 mm criteria. The root mean square deviations (RMSDs) were less than 1%, with a mean of 0.5%, and the source generation time was 4-5 times faster. A seven-field intensity modulated radiation therapy patient treatment achieved 95% chi-test agreement above the 10% isodose with 1% / 1 mm criteria, 99.8% for 2% / 2 mm, a RMSD of 0.8%, and source generation speed-up factor of 2.5.

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.

Automatic commissioning of a GPU-based Monte Carlo radiation dose calculation code for photon radiotherapy

Monte Carlo (MC) simulation is commonly considered as the most accurate method for radiation dose calculations. Commissioning of a beam model in the MC code against a clinical linear accelerator beam is of crucial importance for its clinical implementation. In this paper, we propose an automatic commissioning method for our GPU-based MC dose engine, gDPM. gDPM utilizes a beam model based on a concept of phase-space-let (PSL). A PSL contains a group of particles that are of the same type and close in space and energy. A set of generic PSLs was generated by splitting a reference phase-space file. Each PSL was associated with a weighting factor, and in dose calculations the particle carried a weight corresponding to the PSL where it was from. Dose for each PSL in water was pre-computed, and hence the dose in water for a whole beam under a given set of PSL weighting factors was the weighted sum of the PSL doses. At the commissioning stage, an optimization problem was solved to adjust the PSL weights in order to minimize the difference between the calculated dose and measured one. Symmetry and smoothness regularizations were utilized to uniquely determine the solution. An augmented Lagrangian method was employed to solve the optimization problem. To validate our method, a phase-space file of a Varian TrueBeam 6 MV beam was used to generate the PSLs for 6 MV beams. In a simulation study, we commissioned a Siemens 6 MV beam on which a set of field-dependent phase-space files was available. The dose data of this desired beam for different open fields and a small off-axis open field were obtained by calculating doses using these phase-space files. The 3D γ-index test passing rate within the regions with dose above 10% of dmax dose for those open fields tested was improved averagely from 70.56 to 99.36% for 2%/2 mm criteria and from 32.22 to 89.65% for 1%/1 mm criteria. We also tested our commissioning method on a six-field head-and-neck cancer IMRT plan. The passing rate of the γ-index test within the 10% isodose line of the prescription dose was improved from 92.73 to 99.70% and from 82.16 to 96.73% for 2%/2 mm and 1%/1 mm criteria, respectively. Real clinical data measured from Varian, Siemens, and Elekta linear accelerators were also used to validate our commissioning method and a similar level of accuracy was achieved.

The accuracy of Acuros XB algorithm for radiation beams traversing a metallic hip implant - comparison with measurements and Monte Carlo calculations

In this study, the clinical benefit of the improved accuracy of the Acuros XB (AXB) algorithm, implemented in a commercial radiotherapy treatment planning system (TPS), Varian Eclipse, was demonstrated with beams traversing a high‐Z material. This is also the first study assessing the accuracy of the AXB algorithm applying volumetric modulated arc therapy (VMAT) technique compared to full Monte Carlo (MC) simulations. In the first phase the AXB algorithm was benchmarked against point dosimetry, film dosimetry, and full MC calculation in a water‐filled anthropometric phantom with a unilateral hip implant. Also the validity of the full MC calculation used as reference method was demonstrated. The dose calculations were performed both in original computed tomography (CT) dataset, which included artifacts, and in corrected CT dataset, where constant Hounsfield unit (HU) value assignment for all the materials was made. In the second phase, a clinical treatment plan was prepared for a prostate cancer patient with a unilateral hip implant. The plan applied a hybrid VMAT technique that included partial arcs that avoided passing through the implant and static beams traversing the implant. Ultimately, the AXB‐calculated dose distribution was compared to the recalculation by the full MC simulation to assess the accuracy of the AXB algorithm in clinical setting. A recalculation with the anisotropic analytical algorithm (AAA) was also performed to quantify the benefit of the improved dose calculation accuracy of type ‘c’ algorithm (AXB) over type ‘b’ algorithm (AAA). The agreement between the AXB algorithm and the full MC model was very good inside and in the vicinity of the implant and elsewhere, which verifies the accuracy of the AXB algorithm for patient plans with beams traversing through high‐Z material, whereas the AAA produced larger discrepancies.

PACS numbers: 87.55.‐x, 87.55.D‐, 87.55.K‐, 87.55.kd, 87.55.Qr

Monte Carlo modeling of HD120 multileaf collimator on Varian TrueBeam linear accelerator for verification of 6X and 6X FFF VMAT SABR treatment plans

A Monte Carlo (MC) validation of the vendor-supplied Varian TrueBeam 6 MV flattened (6X) phase-space file and the first implementation of the Siebers-Keall MC MLC model as applied to the HD120 MLC (for 6X flat and 6X flattening filter-free (6X FFF) beams) are described. The MC model is validated in the context of VMAT patient-specific quality assurance. The Monte Carlo commissioning process involves: 1) validating the calculated open-field percentage depth doses (PDDs), profiles, and output factors (OF), 2) adapting the Siebers-Keall MLC model to match the new HD120-MLC geometry and material composition, 3) determining the absolute dose conversion factor for the MC calculation, and 4) validating this entire linac/MLC in the context of dose calculation verification for clinical VMAT plans. MC PDDs for the 6X beams agree with the measured data to within 2.0% for field sizes ranging from 2 × 2 to 40 × 40 cm2. Measured and MC profiles show agreement in the 50% field width and the 80%-20% penumbra region to within 1.3 mm for all square field sizes. MC OFs for the 2 to 40 cm2 square fields agree with measurement to within 1.6%. Verification of VMAT SABR lung, liver, and vertebra plans demonstrate that measured and MC ion chamber doses agree within 0.6% for the 6X beam and within 2.0% for the 6X FFF beam. A 3D gamma factor analysis demonstrates that for the 6X beam, > 99% of voxels meet the pass criteria (3%/3 mm). For the 6X FFF beam, > 94% of voxels meet this criteria. The TrueBeam accelerator delivering 6X and 6X FFF beams with the HD120 MLC can be modeled in Monte Carlo to provide an independent 3D dose calculation for clinical VMAT plans. This quality assurance tool has been used clinically to verify over 140 6X and 16 6X FFF TrueBeam treatment plans.