Monte Carlo Beam Model Research Articles

First Monte Carlo beam model for ultra-high dose rate radiotherapy with a compact electron LINAC.

FLASH radiotherapy based on ultra-high dose rate (UHDR) is actively being studied by the radiotherapy community. Dedicated UHDR electron devices are currently a mainstay for FLASH studies. To present the first Monte Carlo (MC) electron beam model for the UHDR capable Mobetron (FLASH-IQ) as a dose calculation and treatment planning platform for preclinical research and FLASH-radiotherapy (RT) clinical trials. The initial beamline geometry of the Mobetron was provided by the manufacturer, with the first-principal implementation realized in the Geant4-based GAMOS MC toolkit. The geometry and electron source characteristics, such as energy spectrum and beamline parameters, were tuned to match the central-axis percentage depth dose (PDD) and lateral profiles for the pristine beam measured during machine commissioning. The thickness of the small foil in secondary scatter affected the beam model dominantly and was fine tuned to achieve the best agreement with commissioning data. Validation of the MC beam modeling was performed by comparing the calculated PDDs and profiles with EBT-XD radiochromic film measurements for various combinations of applicators and inserts. The nominal 9 MeV electron FLASH beams were best represented by a Gaussian energy spectrum with mean energy of 9.9 MeV and variance (σ) of 0.2 MeV. Good agreement between the MC beam model and commissioning data were demonstrated with maximal discrepancy<3% for PDDs and profiles. Hundred percent gamma pass rate was achieved for all PDDs and profiles with the criteria of 2mm/3%. With the criteria of 2mm/2%, maximum, minimum and mean gamma pass rates were (100.0%, 93.8%, 98.7%) for PDDs and (100.0%, 96.7%, 99.4%) for profiles, respectively. A validated MC beam model for the UHDR capable Mobetron is presented for the first time. The MC model can be utilized for direct dose calculation or to generate beam modeling input required for treatment planning systems for FLASH-RT planning. The beam model presented in this work should facilitate translational and clinical FLASH-RT for trials conducted on the Mobetron FLASH-IQ platform.

Evaluation of monte carlo to support commissioning of the treatment planning system of new pencil beam scanning proton therapy facilities

ABSTRACT Objective. To demonstrate the potential of Monte Carlo (MC) to support the resource-intensive measurements that comprise the commissioning of the treatment planning system (TPS) of new proton therapy facilities. Approach. Beam models of a pencil beam scanning system (Varian ProBeam) were developed in GATE (v8.2), Eclipse proton convolution superposition algorithm (v16.1, Varian Medical Systems) and RayStation MC (v12.0.100.0, RaySearch Laboratories), using the beam commissioning data. All models were first benchmarked against the same commissioning data and validated on seven spread-out Bragg peak (SOBP) plans. Then, we explored the use of MC to optimise dose calculation parameters, fully understand the performance and limitations of TPS in homogeneous fields and support the development of patient-specific quality assurance (PSQA) processes. We compared the dose calculations of the TPSs against measurements (DDTPSvs.Meas.) or GATE (DDTPSvs.GATE) for an extensive set of plans of varying complexity. This included homogeneous plans with varying field-size, range, width, and range-shifters (RSs) (n = 46) and PSQA plans for different anatomical sites (n = 11). Main results. The three beam models showed good agreement against the commissioning data, and dose differences of 3.5% and 5% were found for SOBP plans without and with RSs, respectively. DDTPSvs.Meas. and DDTPSvs.GATE were correlated in most scenarios. In homogeneous fields the Pearson’s correlation coefficient was 0.92 and 0.68 for Eclipse and RayStation, respectively. The standard deviation of the differences between GATE and measurements (±0.5% for homogeneous and ±0.8% for PSQA plans) was applied as tolerance when comparing TPSs with GATE. 72% and 60% of the plans were within the GATE predicted dose difference for both TPSs, for homogeneous and PSQA cases, respectively. Significance. Developing and validating a MC beam model early on into the commissioning of new proton therapy facilities can support the validation of the TPS and facilitate comprehensive investigation of its capabilities and limitations.

Commissioning and Initial Validation of Commercial Treatment Planning System for the Electron FLASH Research Extension

The aim of this study was to investigate the feasibility of commissioning the 16 MeV electron FLASH beam in a commercial treatment planning system (TPS) for pre-clinical research purposes. The delivery system consisted of a new commercial solution for which a linear accelerator was modified into a FLASH Research Extension platform. Additionally, preliminary radiation biology results were highlighted to showcase the future use of this system. To commission a commercial electron Monte Carlo (MC) for dose calculation of a 16 MeV FLASH beam in the TPS, radiochromic film was used to measure the vendor-required beam data, e.g., profiles and percent depth dose (PDD) curves for cone sizes of 6 × 6 cm2, 10 × 10 cm2, and 15 × 15 cm2 as well as an in-air profile for a 40 × 40 cm2 open field (no cone). Once the electron MC beam model was generated, additional measurements were collected for validation and compared against the calculated dose from the TPS. A treatment planning comparison between the newly commissioned FLASH beam and the conventional electron beam was conducted. Specifically, the dose-volume histograms (DVHs) for target volumes and organs at risk were investigated for skin cancer cases previously treated with conventional electron beams. Lastly, the FLASH dose distribution predicted by the electron MC for an in vitro cell study setup was validated with radiochromic film measurements, and initial radiobiology tests were conducted using FLASH and conventional dose-rate electron beams. The electron MC calculated dose for the 16 MeV electron FLASH beam agreed with measured PDDs within 1% for all field sizes. The beam profile characteristics, such as penumbra, shape, and full width at half maximum, demonstrated good agreement with less than 0.5 mm difference between the TPS and measurements. There were noticeable differences in the profiles of large fields between the FLASH and conventional dose-rate beam models due to the more forward-peaked FLASH beam. For treatment planning, Regarding DVH, the FLASH dose-rate plan provided comparable plan quality to the conventional dose-rate plan, achieving adequate coverage for the target volumes and sparing the healthy organs and tissues. The electron MC dose prediction for the FLASH beam was also found to be in good agreement with the film measurements of the in vitro cell study setup. Furthermore, the FLASH beam was observed to be more effective with a 20 % increase in killing pancreatic cancer cells compared to the conventional dose rate. The study successfully incorporated the 16 MeV electron FLASH Research Extension into the commercial TPS using electron Monte Carlo for dose calculation. This will be valuable for pre-clinical cell and animal studies. This research also enables FLASH treatment planning studies, a key component for the future implementation of FLASH into clinical care. Further research is necessary to incorporate the radiation biology effect of FLASH into the treatment planning system.

Auto-commissioning of a Monte Carlo electron beam model with application to photon MLC shaped electron fields

Objective. Presently electron beam treatments are delivered using dedicated applicators. An alternative is the usage of the already installed photon multileaf collimator (pMLC) enabling efficient electron treatments. Currently, the commissioning of beam models is a manual and time-consuming process. In this work an auto-commissioning procedure for the Monte Carlo (MC) beam model part representing the beam above the pMLC is developed for TrueBeam systems with electron energies from 6 to 22 MeV. Approach. The analytical part of the electron beam model includes a main source representing the primary beam and a jaw source representing the head scatter contribution each consisting of an electron and a photon component, while MC radiation transport is performed for the pMLC. The auto-commissioning of this analytical part relies on information pre-determined from MC simulations, in-air dose profiles and absolute dose measurements in water for different field sizes and source to surface distances (SSDs). For validation calculated and measured dose distributions in water were compared for different field sizes, SSDs and beam energies for eight TrueBeam systems. Furthermore, a sternum case in an anthropomorphic phantom was considered and calculated and measured dose distributions were compared at different SSDs. Main results. Instead of the manual commissioning taking up to several days of calculation time and several hours of user time, the auto-commissioning is carried out in a few minutes. Measured and calculated dose distributions agree generally within 3% of maximum dose or 2 mm. The gamma passing rates for the sternum case ranged from 96% to 99% (3% (global)/2 mm criteria, 10% threshold). Significance. The auto-commissioning procedure was successfully implemented and applied to eight TrueBeam systems. The newly developed user-friendly auto-commissioning procedure allows an efficient commissioning of an MC electron beam model and eases the usage of advanced electron radiotherapy utilizing the pMLC for beam shaping.

Sensitivity Analysis of a 6 MeV Photon Beam Monte Carlo Model.

This study aimed to optimize efficiency in Monte Carlo (MC) simulation using sensitivity analysis of a beam model. The BEAMnrc-based model of 6 MV beam of a Siemens Primus linac was developed. For sensitivity analysis, the effect of the electron source, treatment head, and virtual phantom specifications on calculated percent depth dose (PDD) and lateral dose profiles was evaluated. The optimum mean energy (E) and the full width at half maximum (FWHM) of the intensity distribution of the electron beam were calculated as 6.7 MeV and 3 mm, respectively. Increasing E from 6.1 to 6.7 MeV, increased the PDD in the fall-off region by 4.70% and decreased the lateral profile by 8.76%. Changing the FWHM had a significant effect on the buildup region of PDD and the horns and out-of-field regions of the lateral profile. Increasing the collimators opening by 0.5 mm, PDD increased by 2.13% and the central and penumbra regions of profiles decreased by 1.98% and 11.40% respectively. Collimator properties such as thickness and density were effective in changing the penumbra (11.32% for 0.25 cm increment) and the out-of-field (22.82% for 3 g/cm3) regions of the lateral profiles. Analysis of a 6 MV model showed that PDD profiles were more sensitive to changes in energy than to FWHM of the electron source. The lateral profiles were sensitive to E, FWHM, and collimator opening. The density of the collimator affected only the out-of-field region of lateral profiles. The findings of this study may be used to make benchmarking of an MC beam model more efficient.

Characterization of the PTB ultra-high pulse dose rate reference electron beam

Purpose. This investigation aims to present the characterisation and optimisation of an ultra-high pulse dose rate (UHPDR) electron beam at the PTB facility in Germany. A Monte Carlo beam model has been developed for dosimetry study for future investigation in FLASH radiotherapy and will be presented. Material and methods. The 20 MeV electron beams generated by the research linear accelerator has been characterised both in-beamline with profile monitors and magnet spectrometer, and in-water with a diamond detector prototype. The Monte Carlo model has been used to investigate six different setups to enable different dose per pulse (DPP) ranges and beam sizes in water. The properties of the electron radiation field in water have also been characterised in terms of beam size, quality specifier R 50 and flatness. The beam stability has also been studied. Results. The difference between the Monte-Carlo simulated and measured R 50 was smaller than 0.5 mm. The simulated beam sizes agreed with the measured ones within 2 mm. Two suitable setups have been identified for delivering reference UHPDR electron beams. The first one is characterised by a SSD of 70 cm, while in the second one an SSD of 90 cm is used in combination with a 2 mm aluminium scattering plates. The two set-ups are quick and simple to install and enable an expected overall DPP range from 0.13 Gy up to 6.7 Gy per pulse. Conclusion. The electron beams generated by the PTB research accelerator have shown to be stable throughout the four-months length of this investigation. The Monte Carlo models have shown to be in good agreement for beam size and depth dose and within 1% for the beam flatness. The diamond detector prototype has shown to be a promising tool to be used for relative measurements in UHPDR electron beams.

The Evolution of Lateral Dose Distributions of Helium Ion Beams in Air: From Measurement and Modeling to Their Impact on Treatment Planning

To start clinical trials with the first clinical treatment planning system supporting raster-scanned helium ion therapy, a comprehensive database of beam characteristics and parameters was required for treatment room-specific beam physics modeling at the Heidelberg Ion-Beam Therapy Center (HIT). At six different positions in the air gap along the beam axis, lateral beam profiles were systematically measured for 14 initial beam energies covering the full range of available energies at HIT. The 2D-array of liquid-filled ionization chambers OCTAVIUS from PTW was irradiated by a pencil beam focused at the central axis. With a full geometric representation of HIT’s monitoring chambers and beamline elements in FLUKA, our Monte Carlo beam model matches the measured lateral beam profiles. A second set of measurements with the detector placed in a water tank was used to validate the adjustments of the initial beam parameters assumed in the FLUKA simulation. With a deviation between simulated and measured profiles below ±0.8 mm for all investigated beam energies, the simulated profiles build part of the database for the first clinical treatment planning system for helium ions. The evolution of beamwidth was also compared to similar simulations of the clinically available proton and carbon beam. This allows a choice of treatment modality based on quantitative estimates of the physical beam properties. Finally, we investigated the influence of beamwidth variation on patient treatment plans in order to estimate the relevance and necessary precision limits for lateral beam width models.

Treatment Planning System for Electron FLASH Radiation Therapy: Open-Source for Clinical Implementation

To present a Monte Carlo (MC) beam model and its implementation in a clinical treatment planning system (TPS, Varian Eclipse) for a modified ultrahigh dose-rate electron FLASH radiation therapy (eFLASH-RT) linear accelerator (LINAC) using clinical accessories and geometry. The gantry head without scattering foils or targets, representative of the LINAC modifications, was modeled in the Geant4-based GAMOS MC toolkit. The energy spectrum (σE) and beam source emittance cone angle (θcone) were varied to match the calculated open-field central-axis percent depth dose (PDD) and lateral profiles with Gafchromic film measurements. The beam model and its Eclipse configuration were validated with measured profiles of the open field and nominal fields for clinical applicators. An MC forward dose calculation was conducted for a mouse whole-brain treatment, and an eFLASH-RT plan was compared with a conventional (Conv-) RT electron plan in Eclipse for a human patient with metastatic renal cell carcinoma. The eFLASH beam model agreed best with measurements at σE=0.5 MeV and θcone=3.9° ± 0.2°. The model and its Eclipse configuration were validated to clinically acceptable accuracy (the absolute average error was within 1.5% for in-water lateral, 3% for in-air lateral, and 2% for PDDs). The forward calculation showed adequate dose delivery to the entire mouse brain while sparing the organ at risk (lung). The human patient case demonstrated the planning capability with routine accessories to achieve an acceptable plan (90% of the tumor volume receiving 95% and 90% of the prescribed dose for eFLASH and Conv-RT, respectively). To our knowledge, this is the first functional beam model commissioned in a clinical TPS for eFLASH-RT enabling planning and evaluation with minimal deviation from the Conv-RT workflow. It facilitates the clinical translation because eFLASH-RT and Conv-RT plan quality were comparable for a human patient involving complex geometries and tissue heterogeneity. The methods can be expanded to model other eFLASH irradiators with different beam characteristics.

Treatment Planning System for Clinical Translation of Electron FLASH Radiotherapy

Electron ultrahigh dose rate FLASH radiotherapy (UHDR/eFLASH-RT) is currently being extensively tested and modeled, and was studied in several rodent and large animal experiments and one clinical cancer patient. These studies suggest prospective dosimetry assessments are necessary to ensure safe delivery and widespread clinical translation of eFLASH-RT would clearly benefit from prior dose prediction and a treatment planning process that accounts for anatomical heterogeneity and complex geometries. In this study, a GAMOS Monte Carlo beam model was developed and implemented in a clinical treatment planning system for a modified UHDR electron LINAC to demonstrate eFLASH-RT planning of patients.The gantry head was modelled without scattering foils or targets, representative of the LINAC modifications. Centered at the nominal energy of 10 MeV, the energy spectrum (σE) was varied to match the central axis percent depth dose (PDD) profiles of an open field measured with radiographic film. The beam's emittance cone angle (θcone) was determined by matching calculated and film measured lateral profiles, at several depths. The beam model and its configuration in the treatment planning system were validated with film measured lateral profiles and PDD's of the open field and nominal fields for clinical applicators. A GAMOS forward dose calculation was conducted for a mouse whole brain treatment and an eFLASH-RT plan was compared to a conventional electron plan for a human patient with a rib metastasis.Our eFLASH beam model agrees best with measurements at σE = 0.5MeV and θcone = 3.9o ± 0.2o. The model and its treatment planning system configuration were validated to clinically acceptable accuracy (the absolute average error was within 1.5% for in water lateral, 3% for in air lateral, and 2% for PDD's). The forward dose calculation showed dose was delivered to the entire mouse brain with adequate conformality (90% of the tumor volume receiving 90% of the prescribed dose). The comparison plans for the patient undergoing treatment for a rib metastasis demonstrated the capability to achieve an acceptable eFLASH-RT plan with routine accessories (90% of the tumor receiving 95% and 90% of the prescribed dose for eFLASH and conventional, respectively) in relatively complex geometry (right posterior oblique treatment with patient specific applicator cutout).To the best of our knowledge, this study presents the first functional beam model commissioned in a clinical TPS for eFLASH-RT, enabling treatment planning and evaluation with minimal deviation from conventional radiotherapy workflow. By producing plans with the quality comparable to conventional plans, specifically for human patients involving complex geometries and tissue heterogeneity, the clinical translation of eFLASH-RT can be expedited. The open-source methods (described in https://github.com/mr3536/eFLASHBeamModeltoTPS) can model other eFLASH irradiators with different beam characteristics.

Computer-assisted beam modeling for particle therapy.

PurposeTo develop a computer‐driven and thus less user‐dependent method, allowing for a simple and straightforward generation of a Monte Carlo (MC) beam model of a scanned proton and carbon ion beam delivery system.MethodsIn a first step, experimental measurements were performed for proton and carbon ion energies in the available energy ranges. Data included depth dose profiles measured in water and spot sizes in air at various isocenter distances. Using an automated regularization‐based optimization process (AUTO‐BEAM), GATE/Geant4 beam models of the respective beam lines were generated. These were obtained sequentially by using least square weighting functions with and without regularization, to iteratively tune the beam parameters energy, energy spread, beam sigma, divergence, and emittance until a user‐defined agreement was reached. Based on the parameter tuning for a set of energies, a beam model was semi‐automatically generated. The resulting beam models were validated for all centers comparing to independent measurements of laterally integrated depth dose curves and spot sizes in air. For one representative center, three‐dimensional dose cubes were measured and compared to simulations. The method was applied on one research as well as four different clinical beam lines for proton and carbon ions of three different particle therapy centers using synchrotron or cyclotron accelerator systems: (a) MedAustron ion therapy center, (b) University Proton Therapy Dresden, and (c) Center Antoine Lacassagne Nice.ResultsParticle beam ranges in the MC beam models agreed on average within 0.2 mm compared to measurements for all energies and beam lines. Spot sizes in air (full‐width at half maximum) at all positions differed by less than 0.4% from the measurements. Dose calculation with the beam model for the clinical beam line at MedAustron agreed better than 1.7% in absolute dose for a representative clinical case treated with protons. For protons, beam model generation, including geometry creation, data conversion, and validation, was possible within three working days. The number of iterations required for the optimization process to converge, was found to be similar for all beam line geometries and particle types.ConclusionThe presented method was demonstrated to work independently of the beam optics behavior of the different beam lines, particle types, and geometries. Furthermore, it is suitable for non‐expert users and requires only limited user interaction. Beam model validation for different beam lines based on different beam delivery systems, showed good agreement.

Commissioning and validation of TOPAS beam model for IBA Proteus-ONE at UFHPTI

Proteus-ONE (P-ONE) is a novel compact single-room system developed by Ion Beam Application (IBA, Belgium) for proton therapy. The objective of this research was to commission and validate Monte Carlo (MC) beam model for P-ONE using the TOol of PArticle Simulation (TOPAS) Monte Carlo software package. Thirty-three energies were determined within P-ONE energy range (70–227 MeV) to measure integral depth dose (IDD), absolute dose and in-air lateral profiles. The beam model was established by matching TOPAS simulation with measurement for each energy. To validate the accuracy of obtained TOPAS beam model, we compared TOPAS simulation with measurement under various scenarios, including different spread-out Bragg peaks, lateral dose curves and gamma passing rates (GPRs) for 3 × 3x3 cm3, 5 × 5x5 cm3, and 7 × 7x7 cm3 uniform target in a water phantom. Good agreement was observed and GPRs reached 100% among all measurement planes with 3%/3 mm criteria and above 95% with 2%/2 mm criteria. Three patient cases including head and neck, breast and prostate were used for gamma analysis and at least GPRs of 99.8% and 95.4% were reached among all measurement planes with the criteria of 3%/3 mm and 2%/2 m, respectively. The results indicated the accuracy and reliability of the obtained TOPAS beam model which has the potential to serve clinical or research purposes.

Development of a Monte Carlo beam model for raster scanning proton beams and dosimetric comparison

Purpose To develop a Monte Carlo (MC) beam model for raster scanning proton beams for dose verification purposes. Methods and materials MC program FLUKA was used in the model. The nominal energy, momentum spread and beam angular distribution in the model were determined by matching the simulation profiles with the measured integral depth dose (IDD) and in air spot size. Dosimetric comparison was done by comparing the measured and simulated dose distributions. The 1 D dose profile of cubic Spread Out Bragg Peak (SOBP) plans, and the 2 D dose distribution of previously treated breast cancer patients’ clinical plans were measured by using Pinpoint chambers and 2 D array ionization chambers, respectively. Corresponding DICOM plan information was utilized for MC simulation. Results The MC results showed good agreement with measurements for the SOBP plans. The absolute comparison of the absorbed dose difference between the MC and the measurement was 0.93%±0.88%. For the patient plans, the overall passing rate of the gamma index analysis (γ-PR) between the MC simulation and measurement with the 2%-2 mm criteria was 97.78%, and only 1 case had a γ-PR less than 90%. With the 3%-3 mm criteria, γ-PR was never below 99% for all cases with and without the range shifter. Conclusions This work described a method for adapting a MC simulation model for a raster scanning proton beam. The good concordance between the simulations and measurements shows that the MC model is an accurate and reliable method. It has the potential to be used for patient specific quality assurance (PSQA) to reduce the beam time for the measurements in water.

Automated Monte-Carlo re-calculation of proton therapy plans using Geant4/Gate: implementation and comparison to plan-specific quality assurance measurements.

Software re-calculation of proton pencil beam scanning plans provides a method of verifying treatment planning system (TPS) dose calculations prior to patient treatment. This study describes the implementation of AutoMC, a Geant4 v10.3.3/Gate v8.1 (Gate-RTion v1.0)-based Monte-Carlo (MC) system for automated plan re-calculation, and presents verification results for 153 patients (730 fields) planned within year one of the proton service at The Christie NHS Foundation Trust. A MC beam model for a Varian ProBeam delivery system with four range-shifter options (none, 2 cm, 3 cm, 5 cm) was derived from beam commissioning data and implemented in AutoMC. MC and TPS (Varian Eclipse v13.7) calculations of 730 fields in solid-water were compared to physical plan-specific quality assurance (PSQA) measurements acquired using a PTW Octavius 1500XDR array and PTW 31021 Semiflex 3D ion chamber. TPS and MC showed good agreement with array measurements, evaluated using γ analyses at 3%, 3 mm with a 10% lower dose threshold:>94% of fields calculated by the TPS and >99% of fields calculated by MC had γ ≤ 1 for>95% of measurement points within the plane. TPS and MC also showed good agreement with chamber measurements of absolute dose, with systematic differences of <1.5% for all range-shifter options. Reliable independent verification of the TPS dose calculation is a valuable complement to physical PSQA and may facilitate reduction of the physical PSQA workload alongside a thorough delivery system quality assurance programme. A Gate/Geant4-based MC system is thoroughly validated against an extensive physical PSQA dataset for 730 clinical fields, showing that clinical implementation of MC for PSQA is feasible.

A GATE/Geant4 beam model for the MedAustron non-isocentric proton treatment plans quality assurance

PurposeTo present a reference Monte Carlo (MC) beam model developed in GATE/Geant4 for the MedAustron fixed beam line. The proposed model includes an absolute dose calibration in Dose-Area-Product (DAP) and it has been validated within clinical tolerances for non-isocentric treatments as routinely performed at MedAustron. Material and MethodsThe proton beam model was parametrized at the nozzle entrance considering optic and energy properties of the pencil beam. The calibration in terms of absorbed dose to water was performed exploiting the relationship between number of particles and DAP by mean of a recent formalism. Typical longitudinal dose distribution parameters (range, distal penumbra and modulation) and transverse dose distribution parameters (spot sizes, field sizes and lateral penumbra) were evaluated. The model was validated in water, considering regular-shaped dose distribution as well as clinical plans delivered in non-isocentric conditions. ResultsSimulated parameters agree with measurements within the clinical requirements at different air gaps. The agreement of distal and longitudinal dose distribution parameters is mostly better than 1 mm. The dose difference in reference conditions and for 3D dose delivery in water is within 0.5% and 1.2%, respectively. Clinical plans were reproduced within 3%. ConclusionA full nozzle beam model for active scanning proton pencil beam is described using GATE/Geant4. Absolute dose calibration based on DAP formalism was implemented. The beam model is fully validated in water over a wide range of clinical scenarios and will be inserted as a reference tool for research and for independent dose calculation in the clinical routine.

Benchmarking a GATE/Geant4 Monte Carlo model for proton beams in magnetic fields.

PurposeMagnetic resonance guidance in proton therapy (MRPT) is expected to improve its current performance. The combination of magnetic fields with clinical proton beam lines poses several challenges for dosimetry, treatment planning and dose delivery. Proton beams are deflected by magnetic fields causing considerable changes in beam trajectories and also a retraction of the Bragg peak positions. A proper prediction and compensation of these effects is essential to ensure accurate dose calculations. This work aims to develop and benchmark a Monte Carlo (MC) beam model for dose calculation of MRPT for static magnetic fields up to 1 T.MethodsProton beam interactions with magnetic fields were simulated using the GATE/Geant4 toolkit. The transport of charged particle in custom 3D magnetic field maps was implemented for the first time in GATE. Validation experiments were done using a horizontal proton pencil beam scanning system with energies between 62.4 and 252.7 MeV and a large gap dipole magnet (B = 0–1 T), positioned at the isocenter and creating magnetic fields transverse to the beam direction. Dose was measured with Gafchromic EBT3 films within a homogeneous PMMA phantom without and with bone and tissue equivalent material slab inserts. Linear energy transfer (LET) quenching of EBT3 films was corrected using a linear model on dose‐averaged LET method to ensure a realistic dosimetric comparison between simulations and experiments. Planar dose distributions were measured with the films in two different configurations: parallel and transverse to the beam direction using single energy fields and spread‐out Bragg peaks. The MC model was benchmarked against lateral deflections and spot sizes in air of single beams measured with a Lynx PT detector, as well as dose distributions using EBT3 films. Experimental and calculated dose distributions were compared to test the accuracy of the model.ResultsMeasured proton beam deflections in air at distances of 465, 665, and 1155 mm behind the isocenter after passing the magnetic field region agreed with MC‐predicted values within 4 mm. Differences between calculated and measured beam full width at half maximum (FWHM) were lower than 2 mm. For the homogeneous phantom, measured and simulated in‐depth dose profiles showed range and average dose differences below 0.2 mm and 1.2%, respectively. Simulated central beam positions and widths differed <1 mm to the measurements with films. For both heterogenous phantoms, differences within 1 mm between measured and simulated central beam positions and widths were obtained, confirming a good agreement of the MC model.ConclusionsA GATE/Geant4 beam model for protons interacting with magnetic fields up to 1 T was developed and benchmarked to experimental data. For the first time, the GATE/Geant4 model was successfully validated not only for single energy beams, but for SOBP, in homogeneous and heterogeneous phantoms. EBT3 film dosimetry demonstrated to be a powerful dosimetric tool, once the film response function is LET corrected, for measurements in‐line and transverse to the beam direction in magnetic fields. The proposed MC beam model is foreseen to support treatment planning and quality assurance (QA) activities toward MRPT.

Electron beam collimation with a photon MLC for standard electron treatments

Standard electron treatments are currently still performed using standard or molded patient-specific cut-outs placed in the electron applicator. Replacing cut-outs and electron applicators with a photon multileaf collimator (pMLC) for electron beam collimation would make standard electron treatments more efficient and would facilitate advanced treatment techniques like modulated electron radiotherapy (MERT) and mixed beam radiotherapy (MBRT). In this work, a multiple source Monte Carlo beam model for pMLC shaped electron beams commissioned at a source-to-surface distance (SSD) of 70 cm is extended for SSDs of up to 100 cm and validated for several Varian treatment units with field sizes typically used for standard electron treatments. Measurements and dose calculations agree generally within 3% of the maximal dose or 2 mm distance to agreement. To evaluate the dosimetric consequences of using pMLC collimated electron beams for standard electron treatments, pMLC-based and cut-out-based treatment plans are created for a left and a right breast boost, a sternum, a testis and a parotid gland case. The treatment plans consist of a single electron field, either alone (1E) or in combination with two 3D conformal tangential photon fields (1E2X). For each case, a pMLC plan with similar treatment plan quality in terms of dose homogeneity to the target and absolute mean dose values to the organs at risk (OARs) compared to a cut-out plan is found. The absolute mean dose to an OAR is slightly increased for pMLC-based compared to cut-out-based 1E plans if the OAR is located laterally close to the target with respect to beam direction, or if a 6 MeV electron beam is used at an extended SSD. In conclusion, treatment plans using cut-out collimation can be replaced by plans of similar treatment plan quality using pMLC collimation with accurately calculated dose distributions.

Independent Monte-Carlo dose calculation for MLC based CyberKnife radiotherapy

This work aims to develop, implement and validate a Monte Carlo (MC)-based independent dose calculation (IDC) framework to perform patient-specific quality assurance (QA) for multi-leaf collimator (MLC)-based CyberKnife® (Accuray Inc., Sunnyvale, CA) treatment plans. The IDC framework uses an XML-format treatment plan as exported from the treatment planning system (TPS) and DICOM format patient CT data, an MC beam model using phase spaces, CyberKnife MLC beam modifier transport using the EGS++ class library, a beam sampling and coordinate transformation engine and dose scoring using DOSXYZnrc.The framework is validated against dose profiles and depth dose curves of single beams with varying field sizes in a water tank in units of cGy/Monitor Unit and against a 2D dose distribution of a full prostate treatment plan measured with Gafchromic EBT3 (Ashland Advanced Materials, Bridgewater, NJ) film in a homogeneous water-equivalent slab phantom. The film measurement is compared to IDC results by gamma analysis using 2% (global)/2 mm criteria. Further, the dose distribution of the clinical treatment plan in the patient CT is compared to TPS calculation by gamma analysis using the same criteria.Dose profiles from IDC calculation in a homogeneous water phantom agree within 2.3% of the global max dose or 1 mm distance to agreement to measurements for all except the smallest field size. Comparing the film measurement to calculated dose, 99.9% of all voxels pass gamma analysis, comparing dose calculated by the IDC framework to TPS calculated dose for the clinical prostate plan shows 99.0% passing rate. IDC calculated dose is found to be up to 5.6% lower than dose calculated by the TPS in this case near metal fiducial markers.An MC-based modular IDC framework was successfully developed, implemented and validated against measurements and is now available to perform patient-specific QA by IDC.

Experimental investigation on the accuracy of plastic scintillators and of the spectrum discrimination method in small photon fields.

To investigate the accuracy of output factor measurements using a commercial (Exradin W1, SI) and a prototype, "in-house" developed, plastic scintillation dosimeter (PSD) in small photon fields. Repetitive detector-specific output factor OFdet measurements were performed in water (parallel to the CAX) using two W1 PSDs (SI), a PTW microLion, a PTW microDiamond and an unshielded diode D1V (SI) to which Monte Carlo calculated corrections factors were applied. Four sets of repetitive measurements were performed with the W1 PSD positioned parallel and perpendicular to the CAX, each set on a different day, and with analytically calculated volume averaging corrections applied. The W1 OFdet measurements were compared to measurements using an "in-house" developed PSD in water (CHUQ) and both were validated against a previously commissioned Monte Carlo beam model in small photon fields. The performance of the spectrum discrimination calibration procedure was evaluated under different fiber orientations and wavelength threshold choices and the impact on the respective OFdet was reported. For all detectors in the study an excellent agreement was observed down to a field size of 1×1 cm2 . For the smallest field size of 0.5×0.5 cm2 , the W1 PSDs presented OFdet readings higher by 3.8 to 5.0% relative to the mean corrected OFdet of the rest of the detectors and by 5.8 to 6.1% relative to the CHUQ PSD. The repetitive W1 OFdet measurements in water (parallel CAX) were higher by 3.9% relative to the OFdet measurements in Solid WaterTM (perpendicular CAX) even after volume averaging corrections were applied, indicating a potential fiber orientation dependency in small fields. Uncertainties in jaw and detector repositioning as well as source variations with time were estimated to be less than 0.9% (1 σ) for the W1 under both orientations. The CHUQ PSD agreed with the MC dose calculations in water, for the smallest field size, within 1.1-1.7% before any corrections and within 0.3-0.8% after volume averaging corrections. The spectrum discrimination method provided reproducible Cherenkov spectra under the different calibration set-ups with noisier spectra extracted if the calibration is performed in water and parallel to the CAX. The impact of fiber orientation and wavelength threshold during calibration on OFdet was in general minimal. Clinically relevant differences were observed between similar scintillator dosimeters in photon fields smaller than 1× 1 cm2 . Further research on PSDs is needed that can explain the origin of these differences especially related to the Cherenkov spectrum dependencies on the optical fiber technical characteristics.

SU-F-T-619: Dose Evaluation of Specific Patient Plans Based On Monte Carlo Algorithm for a CyberKnife Stereotactic Radiosurgery System

Purpose: This study will use Monte Carlo to simulate the Cyberknife system, and intend to develop the third-party tool to evaluate the dose verification of specific patient plans in TPS. Methods: By simulating the treatment head using the BEAMnrc and DOSXYZnrc software, the comparison between the calculated and measured data will be done to determine the beam parameters. The dose distribution calculated in the Raytracing, Monte Carlo algorithms of TPS (Multiplan Ver4.0.2) and in-house Monte Carlo simulation method for 30 patient plans, which included 10 head, lung and liver cases in each, were analyzed. The γ analysis with the combined 3mm/3% criteria would be introduced to quantitatively evaluate the difference of the accuracy between three algorithms. Results: More than 90% of the global error points were less than 2% for the comparison of the PDD and OAR curves after determining the mean energy and FWHM.The relative ideal Monte Carlo beam model had been established. Based on the quantitative evaluation of dose accuracy for three algorithms, the results of γ analysis shows that the passing rates (84.88±9.67% for head,98.83±1.05% for liver,98.26±1.87% for lung) of PTV in 30 plans between Monte Carlo simulation and TPS Monte Carlo algorithms were good. And the passing rates (95.93±3.12%,99.84±0.33% in each) of PTV in head and liver plans between Monte Carlo simulation and TPS Ray-tracing algorithms were also good. But the difference of DVHs in lung plans between Monte Carlo simulation and Ray-tracing algorithms was obvious, and the passing rate (51.263±38.964%) of γ criteria was not good. It is feasible that Monte Carlo simulation was used for verifying the dose distribution of patient plans. Conclusion: Monte Carlo simulation algorithm developed in the CyberKnife system of this study can be used as a reference tool for the third-party tool, which plays an important role in dose verification of patient plans. This work was supported in part by the grant from Chinese Natural Science Foundation (Grant No. 11275105). Thanks for the support from Accuray Corp.

SU‐D‐BRC‐01: An Automatic Beam Model Commissioning Method for Monte Carlo Simulations in Pencil‐Beam Scanning Proton Therapy

Purpose:Monte Carlo (MC) simulation is typically regarded as the most accurate dose calculation method for proton therapy. Yet for real clinical cases, the overall accuracy also depends on that of the MC beam model. Commissioning a beam model to faithfully represent a real beam requires finely tuning a set of model parameters, which could be tedious given the large number of pencil beams to commmission. This abstract reports an automatic beam‐model commissioning method for pencil‐beam scanning proton therapy via an optimization approach.Methods:We modeled a real pencil beam with energy and spatial spread following Gaussian distributions. Mean energy, and energy and spatial spread are model parameters. To commission against a real beam, we first performed MC simulations to calculate dose distributions of a set of ideal (monoenergetic, zero‐size) pencil beams. Dose distribution for a real pencil beam is hence linear superposition of doses for those ideal pencil beams with weights in the Gaussian form. We formulated the commissioning task as an optimization problem, such that the calculated central axis depth dose and lateral profiles at several depths match corresponding measurements. An iterative algorithm combining conjugate gradient method and parameter fitting was employed to solve the optimization problem. We validated our method in simulation studies.Results:We calculated dose distributions for three real pencil beams with nominal energies 83, 147 and 199 MeV using realistic beam parameters. These data were regarded as measurements and used for commission. After commissioning, average difference in energy and beam spread between determined values and ground truth were 4.6% and 0.2%. With the commissioned model, we recomputed dose. Mean dose differences from measurements were 0.64%, 0.20% and 0.25%.Conclusion:The developed automatic MC beam‐model commissioning method for pencil‐beam scanning proton therapy can determine beam model parameters with satisfactory accuracy.