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.

Output factor comparison of Monte Carlo and measurement for Varian TrueBeam 6 MV and 10 MV flattening filter-free stereotactic radiosurgery system

The dose measurements of the small field sizes, such as conical collimators used in stereotactic radiosurgery (SRS), are a significant challenge due to many factors including source occlusion, detector size limitation, and lack of lateral electronic equilibrium. One useful tool in dealing with the small field effect is Monte Carlo (MC) simulation. In this study, we report a comparison of Monte Carlo simulations and measurements of output factors for the Varian SRS system with conical collimators for energies of 6 MV flattening filter‐free (6 MV) and 10 MV flattening filter‐free (10 MV) on the TrueBeam accelerator. Monte Carlo simulations of Varian's SRS system for 6 MV and 10 MV photon energies with cones sizes of 17.5 mm, 15.0 mm, 12.5 mm, 10.0 mm, 7.5 mm, 5.0 mm, and 4.0 mm were performed using EGSnrc (release V4 2.4.0) codes. Varian's version‐2 phase‐space files for 6 MV and 10 MV of TrueBeam accelerator were utilized in the Monte Carlo simulations. Two small diode detectors Edge (Sun Nuclear) and Small Field Detector (SFD) (IBA Dosimetry) were applied to measure the output factors. Significant errors may result if detector correction factors are not applied to small field dosimetric measurements. Although it lacked the machine‐specific kQclin,Qmsrfclin,fmsr correction factors for diode detectors in this study, correction factors were applied utilizing published studies conducted under similar conditions. For cone diameters greater than or equal to 12.5 mm, the differences between output factors for the Edge detector, SFD detector, and MC simulations are within 3.0% for both energies. For cone diameters below 12.5 mm, output factors differences exhibit greater variations.

PACS number(s): 87.55.k, 87.55.Qr

Monte Carlo validation of the TrueBeam 10XFFF phase-space files for applications in lung SABR

PURPOSE

To establish the clinical acceptability of universal Monte Carlo phase-space data for the 10XFFF (flattening filter free) photon beam on the Varian TrueBeam Linac, including previously unreported data for small fields, output factors, and inhomogeneous media. The study was particularly aimed at confirming the suitability for use in simulations of lung stereotactic ablative radiotherapy treatment plans.

METHODS

Monte Carlo calculated percent depth doses (PDDs), transverse profiles, and output factors for the TrueBeam 10 MV FFF beam using generic phase-space data that have been released by the Varian MC research team were compared with in-house measurements and published data from multiple institutions (ten Linacs from eight different institutions). BEAMnrc was used to create field size specific phase-spaces located underneath the jaws. Doses were calculated with DOSXYZnrc in a water phantom for fields ranging from 1 × 1 to 40 × 40 cm(2). Particular attention was paid to small fields (down to 1 × 1 cm(2)) and dose per pulse effects on dosimeter response for high dose rate 10XFFF beams. Ion chamber measurements were corrected for changes in ion collection efficiency (P(ion)) with increasing dose per pulse. MC and ECLIPSE ANISOTROPIC ANALYTICAL ALGORITHM (AAA) calculated PDDs were compared to Gafchromic film measurement in inhomogeneous media (water, bone, lung).

RESULTS

Measured data from all machines agreed with Monte Carlo simulations within 1.0% and 1.5% for PDDs and in-field transverse profiles, respectively, for field sizes >1 × 1 cm(2) in a homogeneous water phantom. Agreements in the 80%-20% penumbra widths were better than 2 mm for all the fields that were compared. For all the field sizes considered, the agreement between their measured and calculated output factors was within 1.1%. Monte Carlo results for dose to water at water/bone, bone/lung, and lung/water interfaces as well as within lung agree with film measurements to within 2.8% for 10 × 10 and 3 × 3 cm(2) field sizes. This represents a significant improvement over the performance of the ECLIPSE AAA.

CONCLUSIONS

The 10XFFF phase-space data offered by the Varian Monte Carlo research team have been validated for clinical use using measured, interinstitutional beam data in water and with film dosimetry in inhomogeneous media.

A surface energy spectral study on the bone heterogeneity and beam obliquity using the flattened and unflattened photon beams

AIM

Using flattened and unflattened photon beams, this study investigated the spectral variations of surface photon energy and energy fluence in the bone heterogeneity and beam obliquity.

BACKGROUND

Surface dose enhancement is a dosimetric concern when using unflattened photon beam in radiotherapy. It is because the unflattened photon beam contains more low-energy photons which are removed by the flattening filter of the flattened photon beam.

MATERIALS AND METHODS

We used a water and bone heterogeneity phantom to study the distributions of energy, energy fluence and mean energy of the 6 MV flattened and unflattened photon beams (field size = 10 cm × 10 cm) produced by a Varian TrueBEAM linear accelerator. These elements were calculated at the phantom surfaces using Monte Carlo simulations. The photon energy and energy fluence calculations were repeated with the beam angle turned from 0° to 15°, 30° and 45° in the water and bone phantom.

RESULTS

Spectral results at the phantom surfaces showed that the unflattened photon beams contained more photons concentrated mainly in the low-energy range (0-2 MeV) than the flattened beams associated with a flattening filter. With a bone layer of 1 cm under the phantom surface and within the build-up region of the 6 MV photon beam, it is found that both the flattened and unflattened beams had slightly less photons in the energy range <0.4 MeV compared to the water phantom. This shows that the presence of the bone decreased the low-energy photon backscatters to the phantom surface. When both the flattened and unflattened photon beams were rotated from 0° to 45°, the number of photon and mean photon energy increased. This indicates that both photon beams became more hardened or penetrate when the beam angle increased. In the presence of bone, the mean energies of both photon beams increased. This is due to the absorption of low-energy photons by the bone, resulting in more beam hardening.

CONCLUSIONS

This study explores the spectral relationships of surface photon energy and energy fluence with bone heterogeneity and beam obliquity for the flattened and unflattened photon beams. The photon spectral information is important in studies on the patient's surface dose enhancement using unflattened photon beams in radiotherapy.

A Monte Carlo approach to validation of FFF VMAT treatment plans for the TrueBeam linac

PURPOSE

To commission and benchmark a vendor-supplied (Varian Medical Systems) Monte Carlo phase-space data for the 6 MV flattening filter free (FFF) energy mode on a TrueBeam linear accelerator for the purpose of quality assurance of clinical volumetric modulated arc therapy (VMAT) treatment plans. A method for rendering the phase-space data compatible with BEAMnrc/DOSXYZnrc simulation software package is presented.

METHODS

Monte Carlo (MC) simulations were performed to benchmark the TrueBeam 6 MV FFF phase space data that have been released by the Varian MC Research team. The simulations to benchmark the phase space data were done in three steps. First, the original phase space which was created on a cylindrical surface was converted into a format that was compatible with BEAMnrc. Second, BEAMnrc was used to create field size specific phase spaces located underneath the jaws. Third, doses were calculated with DOSXYZnrc in a water phantom for fields ranging from 1 × 1 to 40 × 40 cm(2). Calculated percent depth doses (PDD), transverse profiles, and output factors were compared with measurements for all the fields simulated. After completing the benchmarking study, three stereotactic body radiotherapy (SBRT) VMAT plans created with the Eclipse treatment planning system (TPS) were calculated with Monte Carlo. Ion chamber and film measurements were also performed on these plans. 3D gamma analysis was used to compare Monte Carlo calculation with TPS calculations and with film measurement.

RESULTS

For the benchmarking study, MC calculated and measured values agreed within 1% and 1.5% for PDDs and in-field transverse profiles, respectively, for field sizes >1 × 1 cm(2). Agreements in the 80%-20% penumbra widths were better than 2 mm for all the fields that were compared. With the exception of the 1 × 1 cm(2) field, the agreement between measured and calculated output factors was within 1%. It is of note that excellent agreement in output factors for all field sizes including highly asymmetric fields was achieved without accounting for backscatter into the beam monitor chamber. For the SBRT VMAT plans, the agreement between Monte Carlo and ion chamber point dose measurements was within 1%. Excellent agreement between Monte Carlo, treatment planning system and Gafchromic film dose distribution was observed with over 99% of the points in the high dose volume passing the 3%, 3 mm gamma test.

CONCLUSIONS

The authors have presented a method for making the Varian IAEA compliant 6 MV FFF phase space file of the TrueBeam linac compatible with BEAMnrc/DOSXYZnrc. After benchmarking the modified phase space against measurement, they have demonstrated its potential for use in MC based quality assurance of complex delivery techniques.