Dosimetric validation of Monaco treatment planning system on an Elekta VersaHD linear accelerator

Validation of RayStation Monte Carlo dose calculation algorithm for multiple LINACs

PURPOSE

A photon Monte Carlo (MC) model was commissioned for flattened (FF) and flattening filter free (FFF) 6 MV beam energy. The accuracy of this model, as a single model to be used for three beam matched LINACs, was evaluated.

METHODS

Multiple models were created in RayStation v.10A for three linacs equipped with Elekta "Agility" collimator. A clinically commissioned collapsed cone (CC) algorithm (GoldCC), a MC model automatically created from the CC algorithm without further optimization (CCtoMC) and an optimized MC model (GoldMC) were compared with measurements. The validation of the model was performed by following the recommendations of IAEA TRS 430 and comprised of basic validation in a water tank, validation in a heterogeneous phantom and validation of complex IMRT/VMAT paradigms using gamma analysis of calculated and measured dose maps in a 2D-Array.

RESULTS

Dose calculation with the GoldMC model resulted in a confidence level of 3% for point measurements in water tank and heterogeneous phantom for measurements performed in all three linacs. The same confidence level resulted for GoldCC model. Dose maps presented an agreement for all models on par to each other with γ criteria 2%/2mm.

CONCLUSIONS

The GoldMC model showed a good agreement with measured data and is determined to be accurate for clinical use for all three linacs in this study.

Performance evaluation of Monaco radiotherapy treatment planning system using CIRS Thorax Phantom: Dosimetric assessment of flattened and non-flattened photon beams

Aim

The present study was undertaken to evaluate the performance of different algorithms for flattening filter-free (FFF) and flattened (FF) photon beams in three different in-homogeneities.

Materials and Method

Computed tomography (CT) image sets of the CIRS phantom maintained in the SAD setup by placing the ionization chamber in the lung, bone, and tissue regions, respectively, were acquired. The treatment planning system (TPS) calculated and the ionization chamber measured the doses at the center of the chamber (in the three mediums) were recorded for the flattened and non-flattened photon beams.

Results

The results were reported for photon energies of 6 MV, 10 MV, 15 MV, 6 FFF, and 10 FFF of field sizes 5 × 5 cm2, 10 × 10 cm2, and 15 × 15 cm2. In the bone inhomogeneity, the pencil beam algorithm predicted that the maximum dose variation was 4.88% of measured chamber dose in 10-MV photon energy for the field size 10 × 10 cm2. In water inhomogeneity, both the collapsed cone and Monte Carlo algorithm predicted that the maximum dose variation was ± 3% of measured chamber dose in 10-MV photon energy for the field size 10 × 10 cm2 and in 10-MV FFF photon energy for the field size 5 × 5 cm2, whereas in lung inhomogeneity, the pencil beam algorithm predicted that the highest dose variation was - 6.9% of measured chamber dose in 10-MV FFF photon energy for the field size 5 × 5 cm2.

Conclusion

FF and FFF beams performed differently in lung, water, and bone mediums. The assessment of algorithms was conducted using the anthropomorphic phantom; therefore, these findings may help in the selection of appropriate algorithms for particular clinical settings in radiation delivery.

Comparison of four commercial dose calculation algorithms in different evaluation tests

BACKGROUND

Accurate and fast dose calculation is crucial in modern radiation therapy. Four dose calculation algorithms (AAA, AXB, CCC, and MC) are available in Varian Eclipse and RaySearch Laboratories RayStation Treatment Planning Systems (TPSs).

OBJECTIVES

This study aims to evaluate and compare dosimetric accuracy of the four dose calculation algorithms applying to homogeneous and heterogeneous media, VMAT plans (based on AAPM TG-119 test cases), and the surface and buildup regions.

METHODS

The four algorithms are assessed in homogeneous (IAEA-TECDOCE 1540) and heterogeneous (IAEA-TECDOC 1583) media. Dosimetric evaluation accuracy for VMAT plans is then analyzed, along with the evaluation of the accuracy of algorithms applying to the surface and buildup regions.

RESULTS

Tests conducted in homogeneous media revealed that all algorithms exhibit dose deviations within 5% for various conditions, with pass rates exceeding 95% based on recommended tolerances. Additionally, the tests conducted in heterogeneous media demonstrate high pass rates for all algorithms, with a 100% pass rate observed for 6 MV and mostly 100% pass rate for 15 MV, except for CCC, which achieves a pass rate of 94%. The results of gamma index pass rate (GIPR) for dose calculation algorithms in IMRT fields show that GIPR (3% /3 mm) for all four algorithms in all evaluated tests based on TG119, are greater than 97%. The results of the algorithm testing for the accuracy of superficial dose reveal variations in dose differences, ranging from -11.9% to 7.03% for 15 MV and -9.5% to 3.3% for 6 MV, respectively. It is noteworthy that the AXB and MC algorithms demonstrate relatively lower discrepancies compared to the other algorithms.

CONCLUSIONS

This study shows that generally, two dose calculation algorithms (AXB and MC) that calculate dose in medium have better accuracy than other two dose calculation algorithms (CCC and AAA) that calculate dose to water.

Dosimetric evaluation of a treatment planning system using the AAPM Medical Physics Practice Guideline 5.a (MPPG 5.a) validation tests

Verifying the accuracy of the dose calculation algorithm is considered one of the most critical steps in radiotherapy treatment for delivering an accurate dose to the patient. This work aimed to evaluate the dosimetric performance of the treatment planning system (TPS) algorithms; the AAA (v. 15.6), AXB (v. 15.6) and eMC (v. 15.6) following the AAPM medical physics practice guideline 5.a (MPPG 5.a) validation tests package in a Varian iX Linear Accelerator (Linac). A series of tests were developed based on the MPPG 5.a. on a Varian's Eclipse TPS (v. 15.6) (Varian Medical Systems). First, the basic photon and electron tests were validated by comparing the TPS calculated dose with the measurements. Next, for heterogeneity tests, we verified the Computed Tomography number to electron density (CT-to-ED) curve by comparing it with the baseline values, and TPS calculated point doses beyond heterogeneous media were compared to the measurements. Finally, for IMRT/VMAT dose validation tests, clinical reference plans were re-calculated on ArcCheck's virtual phantom (Sun Nuclear Corporation, Melbourne, FL, USA) and exported to the Linac for delivery using the ArcCheck dosimetry system. All validation tests were evaluated following the MPPG 5.a recommended tolerances. In basic dose validation tests, the TPS calculated depth dose profiles agreed well with the measurements, with a minimum gamma passing rate of 95% at 2%/2 mm criteria. However, disagreements are seen in the build-up and penumbra region. Results for most point doses in homogeneous water phantoms were within the MPPG 5.a tolerance. For the heterogeneity tests, the CT-to-ED curve was established, and calculated point doses were all within 3% of the measurements for heterogeneous media for both photon algorithms at three energies. These results are within the MPPG5.a the recommended tolerance of 3%. Moreover, for electron beams, the differences between the calculated and measured point doses averaged 5% and 7%, but were just within the MPPG 5.a tolerance of 7%. For IMRT and VMAT validation tests using a gamma criteria of a 2%/2 mm, IMRT plans showed maximum and minimum passing rates of 98.2% and 97.4%, respectively. Whereas VMAT plans showed maximum and minimum passing rates of 100% and 94.3%, respectively. We conclude that the dosimetric accuracy of the Eclipse TPS (v15.6) algorithm is adequate for clinical use. The MPPG 5.a tests are valuable for evaluating dose calculation accuracy and are very useful for TPS upgrade checks, commissioning tests, and routine TPS QA.

Planning System for the Optimization of Electric Field Delivery using Implanted Electrodes for Brain Tumor Control

BACKGROUND

The use of non-ionizing electric fields from low intensity voltage sources (<10 V) to control malignant tumor growth is showing increasing potential as a cancer treatment modality. A method of applying these low intensity electric fields using multiple implanted electrodes within or adjacent to tumor volumes has been termed as intratumoral modulation therapy (IMT).

PURPOSE

This study explores advancements in the previously established IMT optimization algorithm, and the development of a custom treatment planning system for patient specific IMT. The practicality of the treatment planning system is demonstrated by implementing the full optimization pipeline on a brain phantom with robotic electrode implantation, post-operative imaging, and treatment stimulation.

METHODS

The integrated planning pipeline in 3D Slicer begins with importing and segmenting patient magnetic resonance images (MRI) or computed tomography (CT) images. The segmentation process is manual, followed by a semi-automatic smoothing step that allows the segmented brain and tumor mesh volumes to be smoothed and simplified by applying selected filters. Electrode trajectories are planned manually on the patient MRI or CT by selecting insertion and tip coordinates for a chosen number of electrodes. The electrode tip positions, and stimulation parameters (phase shift and voltage) can then be optimized with the custom semi-automatic IMT optimization algorithm where users can select the prescription electric field, voltage amplitude limit, tissue electrical properties, nearby organs at risk, optimization parameters (electrode tip location, individual contact phase shift and voltage), desired field coverage percent, and field conformity optimization. Tables of optimization results are displayed, and the resulting electric field is visualized as a field-map superimposed on the MR or CT image, with 3D renderings of the brain, tumor, and electrodes. Optimized electrode coordinates are transferred to robotic electrode implantation software to enable planning and subsequent implantation of the electrodes at the desired trajectories.

RESULTS

An IMT treatment planning system was developed that incorporates patient specific MRI or CT, segmentation, volume smoothing, electrode trajectory planning, electrode tip location and stimulation parameter optimization, and results visualization. All previous manual pipeline steps operating on diverse software platforms were coalesced into a single semi-automated 3D Slicer based user interface. Brain phantom validation of the full system implementation was successful in pre-operative planning, robotic electrode implantation, and post-operative treatment planning to adjust stimulation parameters based on actual implant locations. Voltage measurements were obtained in the brain phantom to determine the electrical parameters of the phantom and validate the simulated electric field distribution.

CONCLUSIONS

A custom treatment planning and implantation system for IMT has been developed in this study, and validated on a phantom brain model, providing an essential step in advancing IMT technology towards future clinical safety and efficacy investigations. This article is protected by copyright. All rights reserved.

Beam modeling and commissioning for Monte Carlo photon beam on an Elekta Versa HD LINAC

PURPOSE

This study aims at analyzing beam data commissioning along with verifying Monte Carlo-based treatment planning system on the basis of the manufacturer guidelines for Elekta Versa HD Linear Accelerator. Moreover, evaluating the beam match process in terms of quality index, photon profile and multi leaf collimator (MLC) offset is aimed as well.

MATERIALS AND METHODS

The process of collecting beam data for Monaco 5.51 Treatment Planning System commissioning was done based on the instructions provided by the manufacturer as well as AAPM TG-106. Monte Carlo analysis was done for output factors in water, percent depth dose and beam profiles. A set of eight static and intensity modulated radiation therapy fields were used to verify the MLC parameters.

RESULTS

The difference between the measured and modeled penetrative quality (D₁₀) was achieved to be 0.54%. The output factors for 6 MV photon energy were measured and the difference between the measured and Monte Carlo output results was smaller than 1% for all the fields. The average percentage of passing the gamma criteria for commissioning test fields was (95+-4)%, however, the minimum agreement was 87.5% belonging to "7SEGA". Additionally, the agreement between both LINAC is 96%, however, the second LINAC reveals a positive offset in the point of approximately -4 cm on the x-axis at the crossplane.

CONCLUSION

Test commissioning was successfully verified using a homogeneous phantom for point dose measurements, post modelling MLC parameters and patient QA plans. All plan parameters pass the gamma criteria. 6 MV photon beam was successfully commissioned for Elekta VersaHD LINAC and is ready for clinical implementation.

Comparison between flattening filter-free (FFF) and flattened photon beam VMAT plans for the whole brain radiotherapy (WBRT) with hippocampus sparing

Purpose

To evaluate and investigate the feasibility of flattening filter‐free (FFF) beam for the whole‐brain radiotherapy (WBRT) with hippocampus sparing.

Methods

Eighteen patients with volumetric‐modulated arc therapy (VMAT) plans in FFF and conventional beam modes were included in this study. The prescribed dose was 30 Gy in 10 fractions. The conformity index (CI), heterogeneity index reported by TPS (HI‐M), and homogeneity index (HI) for planning target volume (PTV) were evaluated. Subsequently, the following parameters for PTV were calculated and compared: D_2%, D_98%; the mean dose, maximum dose, and minimal dose for OARs. Plan modulation index, total MUs, and the gamma index were used to evaluate the plan quality.

Results

HI‐M results were similar for the two techniques (1.239 vs. 1.247, respectively, p = 0.048); FFF beam plans yielded lower D2% compared to FF beam plans (3,416.3 cGy vs. 3,437.2 cGy, p = 0.22), mean dose (3,177.5 cGy vs. 3,195.2 cGy, p = 0.009), and CI (0.884 vs. 0.876, p = 0.001) for PTV. Significant differences were observed between the two beam modes (FF model vs. FFF model) for the maximum dose (1,612.9 cGy vs. 1,470.2 cGy, respectively, p < 0.001), minimum dose (987.6 cGy vs. 898.8 cGy, respectively, p < 0.001), and the mean dose (1144.4 cGy vs. 1047.3 cGy, respectively, p < 0.001) to the hippocampus, and the maximum dose to the eyes (2,792.6 cGy vs. 2,751.3 cGy, respectively, p < 0.001). The average total MUs for FFF‐VMAT plans was significantly greater than FF‐VMAT plans. However, differences for the plan modulation index and the gamma index were negligible.

Conclusion

In comparison with FF beam, the FFF beam mode offers a clear benefit with respect to WBRT with hippocampal sparing.

Accurate dosimetric measurement of large extended SSD fields for comparison to TPS models

PURPOSE

There is little evidence in the literature which quantifies the accuracy of Treatment Planning Systems (TPSs) using large fields at extended SSD (eSSD). This paper introduces the approach taken at Christchurch Hospital, New Zealand to validate the use of the Monaco TPS for Total Body Irradiation (TBI) treatments.

METHODS

A purpose-built device for allowing precise movements of block-like phantoms called a Phantom Mobility Device (PMD) was used for collecting measurements at eSSD. These measurements were used for determining the ability of the Monaco TPS (originally validated for SSDs between 80 and 110 cm) to accurately model dose distributions for TBI treatments at Christchurch Hospital on either treatment machine one (T1) or two (T2) with SSD values of 341 and 432.6 and clinically useful field sizes of 120 and 170 cm, respectively.

RESULTS

We found that within the limits of measurement uncertainty the PMD contributed no determinable scatter to the measurements and proved a reliable approach for eSSD dose measurements. Additionally, by applying depth and off-axis distance constraints of use for TPS information it is possible to use the existing Monaco CCC model at eSSD for block phantom geometries. Dose Difference (DD) analysis showed a clinically acceptable agreement between the CCC model and measured data over a range of depths and off-axis distances.

CONCLUSIONS

The PMD was determined to be a useful tool for accurate measurement of extended SSD treatment fields. Monaco TPS CCC model agreed well for block phantoms so future comparisons to anthropomorphic phantoms or patient data are feasible.

Comparison of the RayStation photon Monte Carlo dose calculation algorithm against measured data under homogeneous and heterogeneous irradiation geometries

PURPOSE

This work compares Monte Carlo dose calculations performed using the RayStation treatment planning system against data measured on a Varian Truebeam linear accelerator with 6 MV and 10 MV FFF photon beams.

METHODS

The dosimetric performance of the RayStation Monte Carlo calculations was evaluated in a variety of irradiation geometries employing homogeneous and heterogeneous phantoms. Profile and depth dose comparisons against measurement were carried out in relative mode using the gamma index as a quantitative measure of similarity within the central high dose regions.

RESULTS

The results demonstrate that the treatment planning system dose calculation engine agrees with measurement to within 2%/1 mm for more than 95% of the data points in the high dose regions for all test cases. A systematic underestimation was observed at the tail of the profile penumbra and out of field, with mean differences generally <0.5 mm or 1% of curve dose maximum respectively. Out of field agreement varied between evaluated beam models.

CONCLUSIONS

The RayStation implementation of photon Monte Carlo dose calculations show good agreement with measured data for the range of scenarios considered in this work and is deemed sufficiently accurate for introduction into clinical use.

Interplay effect modeling in stereotactic body radiotherapy treatment of liver cancer using volumetric modulated arc therapy

To model the interplay effect and minimize it by a selection of optimum parameters value using a predictive model for SBRT of liver cancers. Ten cases of liver tumors treated with the VMAT technique were selected retrospectively. The dosimetric error due to the interplay effect was measured with a micro ionization chamber (0.015cm³) in a Quasar phantom simulating the moving tumor. The interplay effect dependent parameter's viz. patient breaths per minute, the amplitude of respiration, fractional dose (FD), plan complexity due to different energies (Relative degree of modulation), degree of modulation due to a different level of dose optimization constraints, and dose rate (DR) were measured. For the predictive model, mathematical equations were modeled in python from 300 combinations of proposed parameters using multivariate regression analysis. It was observed that the dose variation reduced from -8.44% to -5.16% for change in the BPM values from 7 to 31 and similarly for amplitude, the dose variation reduced from -9.44% to -4.93% for change in amplitude value from 16 mm to 2 mm. The DR and FD have a prominent effect with R² values of 0.990 and 0.880 respectively. The calculated mean square errors of equations excluding amplitude for the predictive model were 0.90 and 0.82 whereas those for equations excluding BPM were 1.31 and 1.41 for 6 MV and 10 MV beams respectively. The values of the parameters can be prospectively optimized by the use of the predictive model according to clinical situations, so dose variation can be minimized.

Implementation of carbon fibre treatment couches in the XiO® and Monaco® Treatment Planning Systems

Purpose: Carbon fibre treatment couches on linear accelerators provide a strong, rigid framework for patient support. Patient safety is a priority, therefore the dosimetric properties of treatment couches need to be accurately incorporated in treatment plans, to minimize differences between planned and delivered dose. This study aims to determine the attenuation effect of treatment couches for 3-D Conformal Radiotherapy (3-D CRT) and to validate the implementation thereof in the XiO and Monaco treatment planning systems (TPS).

Material and methods: Attenuation measurements were performed on the ELEKTA Connexion couches of the ELEKTA Precise and Synergy-Agility linear accelerators. Measurements were made at 10° intervals in RMI-457 Solid water (30 cm x 30 cm x 30 cm) using a PTW Farmer-type ionization chamber (TW30013) positioned at the accelerator’s isocentre. The percentage attenuation was calculated as the ratio of the electrometer readings for parallel-opposed fields. The Computed Tomography (CT) data sets of the set-ups were obtained on a Philips Big Bore 16-slice CT scanner and exported to the TPS. The individual couch structures were delineated and electron density (ED) values were assigned using the commissioned CT-to-ED curve. Test treatment plans were generated with 100MU per field at 10° gantry intervals.

Results: The percentage attenuation was determined to be within 2% and 3% for beams perpendicular to the couch surface for XiO and Monaco, respectively. The maximum attenuation was observed for oblique fields which was significantly higher than the manufacturer specified values. TPS validation showed an agreement to 1% for XiO and Monaco. At extreme oblique angles, both planning systems overestimated this effect up to a maximum of 4%.

Conclusions: Couch attenuation differs significantly with gantry angle and beam energy. As a result, the treatment couch models should be included in all treatment planning calculations.

Derivative-based gamma index: a novel methodology for stringent patient-specific quality assurance in the stereotactic treatment planning of liver cancer

Objective:The development of a stringent derivative-based gamma (DBG) index for patient-specific QA in stereotactic radiotherapy treatment planning (SRTP) to account for the spatial change in dose.Methods:Twenty-five patients of liver SBRT were selected retrospectively for this study. Deliberately, two different kinds of treatment planning approaches were used for each patient. Firstly, the treatment plans were generated using a conventional treatment planning (CTP) approach in which the target was covered with a homogeneous dose along with the nominal dose fall-off around the treatment field. Subsequently, the other treatment plans were generated using an SRTP approach with the intent of heterogeneous dose within the target region along with a steeper dose gradient outside the treatment field as much as possible. For both kinds of treatment plans, two dimensional (2D) conventional gamma (CG) and DBG analysis were performed using the 2D ion chamber array and radiochromic film.Results:Difference in the DBG index was statistically significant whereas, for CG analysis, the difference in CG index was insignificant for both types of treatment plans (CTP and SRTP). A significant positive correlation was observed between the difference in the DBG index and the difference in HI for high gamma criteria.Conclusion:The DBG evaluation is found to be more rigorous, and sensitive to the only SRTP. The proposed method could be opted-in the routine clinical practice in addition to CG.Advances in knowledge:DBG is more sensitive to detect the spatial change of dose, especially in high dose gradient regions.

Modeling Elekta VersaHD using the Varian Eclipse treatment planning system for photon beams: A single-institution experience

The aim of this study was to report a single‐institution experience and commissioning data for Elekta VersaHD linear accelerators (LINACs) for photon beams in the Eclipse treatment planning system (TPS). Two VersaHD LINACs equipped with 160‐leaf collimators were commissioned. For each energy, the percent‐depth‐dose (PDD) curves, beam profiles, output factors, leaf transmission factors and dosimetric leaf gaps (DLGs) were acquired in accordance with the AAPM task group reports No. 45 and No. 106 and the vendor‐supplied documents. The measured data were imported into Eclipse TPS to build a VersaHD beam model. The model was validated by creating treatment plans spanning over the full‐spectrum of treatment sites and techniques used in our clinic. The quality assurance measurements were performed using MatriXX, ionization chamber, and radiochromic film. The DLG values were iteratively adjusted to optimize the agreement between planned and measured doses. Mobius, an independent LINAC logfile‐based quality assurance tool, was also commissioned both for routine intensity‐modulated radiation therapy (IMRT) QA and as a secondary check for the Eclipse VersaHD model. The Eclipse‐generated VersaHD model was in excellent agreement with the measured PDD curves and beam profiles. The measured leaf transmission factors were less than 0.5% for all energies. The model validation study yielded absolute point dose agreement between ionization chamber measurements and Eclipse within ±4% for all cases. The comparison between Mobius and Eclipse, and between Mobius and ionization chamber measurements lead to absolute point dose agreement within ±5%. The corresponding 3D dose distributions evaluated with 3%global/2mm gamma criteria resulted in larger than 90% passing rates for all plans. The Eclipse TPS can model VersaHD LINACs with clinically acceptable accuracy. The model validation study and comparisons with Mobius demonstrated that the modeling of VersaHD in Eclipse necessitates further improvement to provide dosimetric accuracy on par with Varian LINACs.

The commissioning and validation of Monaco treatment planning system on an Elekta VersaHD linear accelerator

Accurate beam modeling is essential to help ensure overall accuracy in the radiotherapy process. This study describes our experience with beam model validation of a Monaco treatment planning system on a Versa HD linear accelerator. Data were collected such that Monaco beam models could be generated using three algorithms: collapsed cone (CC) and photon Monte Carlo (MC) for photon beams, and electron Monte Carlo (eMC) for electron beams. Validations are performed on measured percent depth doses (PDDs) and profiles, for open‐field point‐doses in homogenous and heterogeneous media, and for obliquely incident electron beams. Gamma analysis is used to assess the agreement between calculation and measurement for intensity modulated radiation therapy (IMRT) and volumetric modulated arc therapy (VMAT) plans, including volumetric modulated arc therapy for stereotactic body radiation therapy (VMAT SBRT). For all relevant conditions, gamma index values below 1 are obtained when comparing Monaco calculated PDDs and profiles with measured data. Point‐doses in a water medium are found to be within 2% agreement of commissioning data in 99.5% and 98.6% of the points computed by MC and CC, respectively. All point‐dose calculations for the eMC algorithm in water are within 4% agreement of measurement, and 92% of measurements are within 3%. In heterogeneous media of air and cortical bone, both CC and MC yielded better than 3% agreement with ion chamber measurements. eMC yielded 3% agreement to measurement downstream of air with oblique beams of up to 27°, 5% agreement distal to bone, and within 4% agreement at extended source to surface distance (SSD) for all electron energies except 6 MeV. The 6‐MeV point of measurement is on a steep dose gradient which may impact the magnitude of discrepancy measured. The average gamma passing rate for IMRT/VMAT plans is 96.9% (±2.1%) and 98.0% (±1.9%) for VMAT SBRT when evaluated using 3%/2 mm criteria. Monaco beam models for the Versa HD linac were successfully commissioned for clinical use.