Modeling of flattening filter free photon beams with analytical and Monte Carlo TPS

Commissioning of Aktina SRS cones and dosimetric validation of the RayStation photon Monte Carlo dose calculation algorithm

Clinical implementation of SRS cones demands particular experimental care and dosimetric considerations in order to deliver precise and safe radiotherapy to patients. The purpose of this work was to present the commissioning data of recent Aktina cones combined with a 6MV flattened beam produced by an Elekta VersaHD linear accelerator. Additionally, the modelling process, and an assessment of dosimetric accuracy of the RayStation Monte Carlo dose calculation algorithm for cone based SRS was performed. There are currently no studies presenting beam data for this equipment and none that outlines the modelling parameters and validation of dose calculation using RayStation's photon Monte Carlo dose engine with cones. Beam data was measured using an SFD and a microDiamond and benchmarked against EBT3 film for cones of diameter 5-39 mm. Modelling was completed and validated within homogeneous and heterogeneous phantoms. End-to-end image-guided validation was performed using a StereoPHAN™ housing, an SRS MapCHECK and EBT3 film, and calculation time was investigated as a function of statistical uncertainty and field diameter. The TPS calculations agreed with measured data within their estimated uncertainties and clinical treatment plans could be calculated in under a minute. The data presented serves as a reference for others commissioning Aktina stereotactic cones and the modelling parameters serve similarly, while providing a starting point for those commissioning the same TPS algorithm for use with cones. It has been shown in this work that RayStation's Monte Carlo photon dose algorithm performs satisfactorily in the presence of SRS cones.

Acceptance, commissioning and quality assurance of the MRIdian®: Site experience and three years follow-up

PURPOSE

This study presents the methodology and results of the acceptance and periodical quality controls on the MRIdian®.

MATERIALS AND METHODS

The impact of the magnetic field on other machines was investigated by controlling nearby linacs dose profiles. The image quality of the 0.345T MR scanner was evaluated, also assessing the integrated linear accelerator influence. The photon beams lateral and depth dose profiles were measured in motorized water tanks, along dose rate and output factors, and compared to Monte Carlo (MC) calculations. The isocenter position, gantry angles and multi-leaf collimator (MLC) position were controlled using film dosimetry. Gating latency and dosimetric accuracy were controlled with a dynamic phantom.

RESULTS

The magnetic field had no significant impact on other nearby linacs. Image quality was within tolerances and did not vary over time. Dose profiles measured showed good agreement with MC data, with maximum differences of 1.3% in-field. Output factors were within 0.8% of calculated values. Imaging and radiative isocenter matched within 0.9±0.4mm over all monthly controls. Gantry rotation was precise within -0.1±0.2°, with an isocenter variation of 1.4±0.3mm diameter. The average MLC position was within 0.4±0.1mm of theoretical value. Finally, the gating latency was 0.14±0.07sec and the gated dose within 0.3% of base value.

CONCLUSION

All results are within the tolerances fixed by ViewRay® and show low variations over 2 years, comforting the use of small margins and gating for high-dose adaptive treatments.

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.

Validation of a High-Throughput Dicentric Chromosome Assay Using Complex Radiation Exposures

Validation of biodosimetry assays is routinely performed using primarily orthovoltage irradiators at a conventional dose rate of approximately 1 Gy/min. However, incidental/ accidental exposures caused by nuclear weapons can be more complex. The aim of this work was to simulate the DNA damage effects mimicking those caused by the detonation of a several kilotons improvised nuclear device (IND). For this, we modeled complex exposures to: 1. a mixed (photons + IND-neutrons) field and 2. different dose rates that may come from the blast, nuclear fallout, or ground deposition of radionuclides (ground shine). Additionally, we assessed whether myeloid cytokines affect the precision of radiation dose estimation by modulating the frequency of dicentric chromosomes. To mimic different exposure scenarios, several irradiation systems were used. In a mixed field study, human blood samples were exposed to a photon field enriched with neutrons (ranging from 10% to 37%) from a source that mimics Hiroshima's A-bomb's energy spectrum (0.2-9 MeV). Using statistical analysis, we assessed whether photons and neutrons act in an additive or synergistic way to form dicentrics. For the dose rates study, human blood was exposed to photons or electrons at dose rates ranging from low (where the dose was spread over 32 h) to extremely high (where the dose was delivered in a fraction of a microsecond). Potential effects of cytokine treatment on biodosimetry dose predictions were analyzed in irradiated blood subjected to Neupogen or Neulasta for 24 or 48 h at the concentration recommended to forestall manifestation of an acute radiation syndrome in bomb survivors. All measurements were performed using a robotic station, the Rapid Automated Biodosimetry Tool II, programmed to culture lymphocytes and score dicentrics in multiwell plates (the RABiT-II DCA). In agreement with classical concepts of radiation biology, the RABiT-II DCA calibration curves suggested that the frequency of dicentrics depends on the type of radiation and is modulated by changes in the dose rate. The resulting dose-response curves suggested an intermediate dicentric yields and additive effects of photons and IND-neutrons in the mixed field. At ultra-high dose rate (600 Gy/s), affected lymphocytes exhibited significantly fewer dicentrics (P < 0.004, t test). In contrast, we did not find the dose-response modification effects of radiomitigators on the yields of dicentrics (Bonferroni corrected P > 0.006, ANOVA test). This result suggests no bias in the dose predictions should be expected after emergency cytokine treatment initiated up to 48 h prior to blood collection for dicentric analysis.

Evaluation of candidate template beam models for a matched TrueBeam treatment delivery system

PURPOSE

To explore candidate RayStation beam models to serve as a class-specific template for a TrueBeam treatment delivery system.

METHODS

Established validation techniques were used to evaluate three photon beam models: a clinically optimized model from the authors' institution, the built-in RayStation template, and a hybrid consisting of the RayStation template except substituting average MLC parameter values from a recent IROC survey. Comparisons were made for output factors, dose profiles from open fields, as well as representative VMAT test plans.

RESULTS

For jaw-defined output factors, each beam model was within 1.6% of expected published values. Similarly, the majority (57-66%) of jaw-defined dose curves from each model had a gamma pass rate >95% (2% / 3 mm, 20% threshold) when compared to TrueBeam representative beam data. For dose curves from MPPG 5.a MLC-defined fields, average gamma pass rates (1% / 1 mm, 20% threshold) were 92.9%, 85.1%, and 86.0% for the clinical, template, and hybrid models, respectively. For VMAT test plans measured with a diode array detector, median dose differences were 0.6%, 1.3%, and 1.1% for the clinical, template, and hybrid models, respectively. For in-phantom ionization chamber measurements with the same VMAT test plans, the average percent difference was -0.3%, -1.4%, and -1.0% for the clinical, template, and hybrid models, respectively.

CONCLUSION

Beam model templates taken from the vendor and aggregate results within the community were both reasonable starting points, but neither approach was as optimal as a clinically tuned model, the latter producing better agreement with all validation measurements. Given these results, the clinically optimized model represents a better candidate as a consensus template. This can benefit the community by reducing commissioning time and improving dose calculation accuracy for matched TrueBeam treatment delivery systems.

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.

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

Background

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

Methods

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

Results

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

Conclusions

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

Optimizing beam models for dosimetric accuracy over a wide range of treatments

This work presents a systematic approach for testing a dose calculation algorithm over a variety of conditions designed to span the possible range of clinical treatment plans. Using this method, a TrueBeam STx machine with high definition multi-leaf collimators (MLCs) was commissioned in the RayStation treatment planning system (TPS). The initial model parameters values were determined by comparing TPS calculations with standard measured depth dose and profile curves. The MLC leaf offset calibration was determined by comparing measured and calculated field edges utilizing a wide range of MLC retracted and over-travel positions. The radial fluence was adjusted using profiles through both the center and corners of the largest field size, and through measurements of small fields that were located at highly off-axis positions. The flattening filter source was adjusted to improve the TPS agreement for the output of MLC-defined fields with much larger jaw openings. The MLC leaf transmission and leaf end parameters were adjusted to optimize the TPS agreement for highly modulated intensity-modulated radiotherapy (IMRT) plans. The final model was validated for simple open fields, multiple field configurations, the TG 119 C-shape target test, and a battery of clinical IMRT and volumetric-modulated arc therapy (VMAT) plans. The commissioning process detected potential dosimetric errors of over 10% and resulted in a final model that provided in general 3% dosimetric accuracy. This study demonstrates the importance of using a variety of conditions to adjust a beam model and provides an effective framework for achieving high dosimetric accuracy.

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.