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.

Efficiency enhancements of a Monte Carlo beamlet based treatment planning process: implementation and parameter study.

Objective.The computational effort to perform beamlet calculation, plan optimization and final dose calculation of a treatment planning process (TPP) generating intensity modulated treatment plans is enormous, especially if Monte Carlo (MC) simulations are used for dose calculation. The goal of this work is to improve the computational efficiency of a fully MC based TPP for static and dynamic photon, electron and mixed photon-electron treatment techniques by implementing multiple methods and studying the influence of their parameters.Approach.A framework is implemented calculating MC beamlets efficiently in parallel on each available CPU core. The user can specify the desired statistical uncertainty of the beamlets, a fractional sparse dose threshold to save beamlets in a sparse format and minimal distances to the PTV surface from which 2 × 2 × 2 = 8 (medium) or even 4 × 4 × 4 = 64 (large) voxels are merged. The compromise between final plan quality and computational efficiency of beamlet calculation and optimization is studied for several parameter values to find a reasonable trade-off. For this purpose, four clinical and one academic case are considered with different treatment techniques.Main results.Setting the statistical uncertainty to 5% (photon beamlets) and 15% (electron beamlets), the fractional sparse dose threshold relative to the maximal beamlet dose to 0.1% and minimal distances for medium and large voxels to the PTV to 1 cm and 2 cm, respectively, does not lead to substantial degradation in final plan quality compared to using 2.5% (photon beamlets) and 5% (electron beamlets) statistical uncertainty and no sparse format nor voxel merging. Only OAR sparing is slightly degraded. Furthermore, computation times are reduced by about 58% (photon beamlets), 88% (electron beamlets) and 96% (optimization).Significance.Several methods are implemented improving computational efficiency of beamlet calculation and plan optimization of a fully MC based TPP without substantial degradation in final plan quality.

Enabling non-isocentric dynamic trajectory radiotherapy by integration of dynamic table translations.

Objective.The purpose of this study is to develop a treatment planning process (TPP) for non-isocentric dynamic trajectory radiotherapy (DTRT) using dynamic gantry rotation, collimator rotation, table rotation, longitudinal, vertical and lateral table translations and intensity modulation and to validate the dosimetric accuracy.Approach.The TPP consists of two steps. First, a path describing the dynamic gantry rotation, collimator rotation and dynamic table rotation and translations is determined. Second, an optimization of the intensity modulation along the path is performed. We demonstrate the TPP for three use cases. First, a non-isocentric DTRT plan for a brain case is compared to an isocentric DTRT plan in terms of dosimetric plan quality and delivery time. Second, a non-isocentric DTRT plan for a craniospinal irradiation (CSI) case is compared to a multi-isocentric intensity modulated radiotherapy (IMRT) plan. Third, a non-isocentric DTRT plan for a bilateral breast case is compared to a multi-isocentric volumetric modulated arc therapy (VMAT) plan. The non-isocentric DTRT plans are delivered on a TrueBeam in developer mode and their dosimetric accuracy is validated using radiochromic films.Main results.The non-isocentric DTRT plan for the brain case is similar in dosimetric plan quality and delivery time to the isocentric DTRT plan but is expected to reduce the risk of collisions. The DTRT plan for the CSI case shows similar dosimetric plan quality while reducing the delivery time by 45% in comparison with the IMRT plan. The DTRT plan for the breast case showed better treatment plan quality in comparison with the VMAT plan. The gamma passing rates between the measured and calculated dose distributions are higher than 95% for all three plans.Significance.The versatile benefits of non-isocentric DTRT are demonstrated with three use cases, namely reduction of collision risk, reduced setup and delivery time and improved dosimetric plan quality.

TriB-RT: Simultaneous optimization of photon, electron and proton beams.

PURPOSE: To develop a novel treatment planning process (TPP) with simultaneous optimization of modulated photon, electron and proton beams for improved treatment plan quality in radiotherapy. METHODS: A framework for fluence map optimization of Monte Carlo (MC) calculated beamlet dose distributions is developed to generate treatment plans consisting of photon, electron and spot scanning proton fields. Initially, in-house intensity modulated proton therapy (IMPT) plans are compared to proton plans created by a commercial treatment planning system (TPS). A triple beam radiotherapy (TriB-RT) plan is generated for an exemplary academic case and the dose contributions of the three particle types are investigated. To investigate the dosimetric potential, a TriB-RT plan is compared to an in-house IMPT plan for two clinically motivated cases. Benefits of TriB-RT for a fixed proton beam line with a single proton field are investigated. RESULTS: In-house optimized IMPT are of at least equal or better quality than TPS-generated proton plans, and MC-based optimization shows dosimetric advantages for inhomogeneous situations. Concerning TriB-RT, for the academic case, the resulting plan shows substantial contribution of all particle types. For the clinically motivated case, improved sparing of organs at risk close to the target volume is achieved compared to IMPT (e.g. myelon and brainstem [Formula: see text] -37%) at cost of an increased low dose bath (healthy tissue V 10% +22%). In the scenario of a fixed proton beam line, TriB-RT plans are able to compensate the loss in degrees of freedom to substantially improve plan quality compared to a single field proton plan. CONCLUSION: A novel TPP which simultaneously optimizes photon, electron and proton beams was successfully developed. TriB-RT shows the potential for improved treatment plan quality and is especially promising for cost-effective single-room proton solutions with a fixed beamline in combination with a conventional linac delivering photon and electron fields.

Development of an extended Macro Monte Carlo method for efficient and accurate dose calculation in magnetic fields.

MOTIVATION: Progress in the field of magnetic resonance (MR)-guided radiotherapy has triggered the need for fast and accurate dose calculation in presence of magnetic fields. The aim of this work is to satisfy this need by extending the macro Monte Carlo (MMC) method to enable dose calculation for photon, electron, and proton beams in a magnetic field. METHODS: The MMC method is based on the transport of particles in macroscopic steps through an absorber by sampling the relevant physical quantities from a precalculated database containing probability distribution functions. To enable MMC particle transport in a magnetic field, a transformation accounting for the Lorentz force is applied for each macro step by rotating the sampled position and direction around the magnetic field vector. The transformed position and direction distributions on local geometries are validated against full MC for electron and proton pencil beams. To enable photon dose calculation, an in-house MC algorithm is used for photon transport and interaction. Emerging secondary charged particles are passed to MMC for transport and energy deposition. The extended MMC dose calculation accuracy and efficiency is assessed by comparison with EGSnrc (photon and electron beams) and Geant4 (proton beam) calculated dose distributions of different energies and homogeneous magnetic fields for broad beams impinging on water phantoms with bone and lung inhomogeneities. RESULTS: The geometric transformation on the local geometries is able to reproduce the results of full MC for all investigated settings (difference in mean value and standard deviation <1%). Macro Monte Carlo calculated dose distributions in a homogeneous magnetic field are in agreement with EGSnrc and Geant4, respectively, with gamma passing rates >99.6% (global 2%, 2 mm and 10% threshold criteria) for all situations. MMC achieves a substantial efficiency gain of up to a factor of 21 (photon beam), 66 (electron beam), and 356 (proton beam) compared to EGSnrc or Geant4. CONCLUSION: Efficient and accurate dose calculation in magnetic fields was successfully enabled by utilizing the developed extended MMC transport method for photon, electron, and proton beams.

Benchmarking Monte-Carlo dose calculation for MLC CyberKnife treatments.

BACKGROUND: Vendor-independent Monte Carlo (MC) dose calculation (IDC) for patient-specific quality assurance of multi-leaf collimator (MLC) based CyberKnife treatments is used to benchmark and validate the commercial MC dose calculation engine for MLC based treatments built into the CyberKnife treatment planning system (Precision MC). METHODS: The benchmark included dose profiles in water in 15 mm depth and depth dose curves of rectangular MLC shaped fields ranging from 7.6 mm × 7.7 mm to 115.0 mm × 100.1 mm, which were compared between IDC, Precision MC and measurements in terms of dose difference and distance to agreement. Dose distributions of three phantom cases and seven clinical lung cases were calculated using both IDC and Precision MC. The lung PTVs ranged from 14 cm3 to 93 cm3. Quantitative comparison of these dose distributions was performed using dose-volume parameters and 3D gamma analysis with 2% global dose difference and 1 mm distance criteria and a global 10% dose threshold. Time to calculate dose distributions was recorded and efficiency was assessed. RESULTS: Absolute dose profiles in 15 mm depth in water showed agreement between Precision MC and IDC within 3.1% or 1 mm. Depth dose curves agreed within 2.3% / 1 mm. For the phantom and clinical lung cases, mean PTV doses differed from - 1.0 to + 2.3% between IDC and Precision MC and gamma passing rates were > =98.1% for all multiple beam treatment plans. For the lung cases, lung V20 agreed within ±1.5%. Calculation times ranged from 2.2 min (for 39 cm3 PTV at 1.0 × 1.0 × 2.5 mm3 native CT resolution) to 8.1 min (93 cm3 at 1.1 × 1.1 × 1.0 mm3), at 2% uncertainty for Precision MC for the 7 examined lung cases and 4-6 h for IDC, which, however, is not optimized for efficiency but used as a gold standard for accuracy. CONCLUSIONS: Both accuracy and efficiency of Precision MC in the context of MLC based planning for the CyberKnife M6 system were benchmarked against MC based IDC framework. Precision MC is used in clinical practice at our institute.

Investigation on the resolution of a micro cone beam CT scanner scintillating detector using Monte Carlo methods.

The impact of several physical quantities on the spatial resolution of an X-ray scintillating pixel detector for a micro cone beam CT (µCBCT) is investigated and discussed. The XtremeCT from SCANCO Medical AG was simulated using the EGSnrc/EGS++ Monte Carlo (MC) framework and extensively benchmarked in a previous work. The resolution of the detector was determined by simulating a titanium knife-edge to obtain the edge spread function (ESF) and the modulation transfer function (MTF). Propagation of the scintillation light through the scintillator and its coupling into the fiber optics system was taken into account. The contribution of particles scattered in the main scanner components to the detector signal is very low and does not affect the spatial resolution of the detector. The resolution obtained from the energy deposition in the scintillator without any blurring due to the propagation of the scintillation light into the fiber optics array was 31 µm. By assuming isotropic light propagation in the scintillator, the resolution degraded to 360 µm. A simple light propagation model taking into account the impact of the scintillator's columnar microstructures was developed and compared with the MANTIS Monte Carlo simulation package. By reducing the width of the model's light propagation kernel by a factor of 2 compared to the isotropic case, the detector resolution can be improved to 83 µm, which corresponds well to the measured resolution of 86 µm. The resolution of the detector is limited mainly by the propagation of the scintillation light through the scintillator layer. It offers the greatest potential to improve the resolution of the µCBCT imaging system.

Part 2: Dynamic mixed beam radiotherapy (DYMBER): Photon dynamic trajectories combined with modulated electron beams.

PURPOSE: The purpose of this study was to develop a treatment technique for dynamic mixed beam radiotherapy (DYMBER) utilizing increased degrees of freedom (DoF) of a conventional treatment unit including different particle types (photons and electrons), intensity and energy modulation and dynamic gantry, table, and collimator rotations. METHODS: A treatment planning process has been developed to create DYMBER plans combining photon dynamic trajectories (DTs) and step and shoot electron apertures collimated with the photon multileaf collimator (pMLC). A gantry-table path is determined for the photon DTs with minimized overlap of the organs at risk (OARs) with the target. In addition, an associated dynamic collimator rotation is established with minimized area between the pMLC leaves and the target contour. pMLC sequences of photon DTs and electron pMLC apertures are then simultaneously optimized using direct aperture optimization (DAO). Subsequently, the final dose distribution of the electron pMLC apertures is calculated using the Swiss Monte Carlo Plan (SMCP). The pMLC sequences of the photon DTs are then re-optimized with a finer control point resolution and with the final electron dose distribution taken into account. Afterwards, the final photon dose distribution is calculated also using the SMCP and summed together with the one of the electrons. This process is applied for a brain and two head and neck cases. The resulting DYMBER dose distributions are compared to those of dynamic trajectory radiotherapy (DTRT) plans consisting only of photon DTs and clinically applied VMAT plans. Furthermore, the deliverability of the DYMBER plans is verified in terms of dosimetric accuracy, delivery time and collision avoidance. For this purpose, The DYMBER plans are delivered to Gafchromic EBT3 films placed in an anthropomorphic head phantom on a Varian TrueBeam linear accelerator. RESULTS: For each case, the dose homogeneity in the target is similar or better for DYMBER compared to DTRT and VMAT. Averaged over all three cases, the mean dose to the parallel OARs is 16% and 28% lower, D2% to the serial OARs is 17% and 37% lower and V10% to normal tissue is 12% and 4% lower for the DYMBER plans compared to the DTRT and VMAT plans, respectively. The DYMBER plans are delivered without collision and with a 4-5 min longer delivery time than the VMAT plans. The absolute dose measurements are compared to calculation by gamma analysis using 2% (global)/2 mm criteria with passing rates of at least 99%. CONCLUSIONS: A treatment technique for DYMBER has been successfully developed and verified for its deliverability. The dosimetric superiority of DYMBER over DTRT and VMAT indicates utilizing increased DoF to be the key to improve brain and head and neck radiation treatments in future.

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.

Simultaneous optimization of photons and electrons for mixed beam radiotherapy.

The aim of this work is to develop and investigate an inverse treatment planning process (TPP) for mixed beam radiotherapy (MBRT) capable of performing simultaneous optimization of photon and electron apertures. A simulated annealing based direct aperture optimization (DAO) is implemented to perform simultaneous optimization of photon and electron apertures, both shaped with the photon multileaf collimator (pMLC). Validated beam models are used as input for Monte Carlo dose calculations. Consideration of photon pMLC transmission during DAO and a weight re-optimization of the apertures after deliverable dose calculation are utilized to efficiently reduce the differences between optimized and deliverable dose distributions. The TPP for MBRT is evaluated for an academic situation with a superficial and an enlarged PTV in the depth, a left chest wall case including the internal mammary chain and a squamous cell carcinoma case. Deliverable dose distributions of MBRT plans are compared to those of modulated electron radiotherapy (MERT), photon IMRT and if available to those of clinical VMAT plans. The generated MBRT plans dosimetrically outperform the MERT, photon IMRT and VMAT plans for all investigated situations. For the clinical cases of the left chest wall and the squamous cell carcinoma, the MBRT plans cover the PTV similarly or more homogeneously than the VMAT plans, while OARs are spared considerably better with average reductions of the mean dose to parallel OARs and D 2% to serial OARs by 54% and 26%, respectively. Moreover, the low dose bath expressed as V 10% to normal tissue is substantially reduced by up to 45% compared to the VMAT plans. A TPP for MBRT including simultaneous optimization is successfully implemented and the dosimetric superiority of MBRT plans over MERT, photon IMRT and VMAT plans is demonstrated for academic and clinical situations including superficial targets with and without deep-seated part.

Assessing dose rate distributions in VMAT plans.

Dose rate is an essential factor in radiobiology. As modern radiotherapy delivery techniques such as volumetric modulated arc therapy (VMAT) introduce dynamic modulation of the dose rate, it is important to assess the changes in dose rate. Both the rate of monitor units per minute (MU rate) and collimation are varied over the course of a fraction, leading to different dose rates in every voxel of the calculation volume at any point in time during dose delivery. Given the radiotherapy plan and machine specific limitations, a VMAT treatment plan can be split into arc sectors between Digital Imaging and Communications in Medicine control points (CPs) of constant and known MU rate. By calculating dose distributions in each of these arc sectors independently and multiplying them with the MU rate, the dose rate in every single voxel at every time point during the fraction can be calculated. Independently calculated and then summed dose distributions per arc sector were compared to the whole arc dose calculation for validation. Dose measurements and video analysis were performed to validate the calculated datasets. A clinical head and neck, cranial and liver case were analyzed using the tool developed. Measurement validation of synthetic test cases showed linac agreement to precalculated arc sector times within ±0.4 s and doses ±0.1 MU (one standard deviation). Two methods for the visualization of dose rate datasets were developed: the first method plots a two-dimensional (2D) histogram of the number of voxels receiving a given dose rate over the course of the arc treatment delivery. In similarity to treatment planning system display of dose, the second method displays the dose rate as color wash on top of the corresponding computed tomography image, allowing the user to scroll through the variation over time. Examining clinical cases showed dose rates spread over a continuous spectrum, with mean dose rates hardly exceeding 100 cGy min(-1) for conventional fractionation. A tool to analyze dose rate distributions in VMAT plans with sub-second accuracy was successfully developed and validated. Dose rates encountered in clinical VMAT test cases show a continuous spectrum with a mean less than or near 100 cGy min(-1) for conventional fractionation.

Beamlet based direct aperture optimization for MERT using a photon MLC.

PURPOSE: A beamlet based direct aperture optimization (DAO) for modulated electron radiotherapy (MERT) using photon multileaf collimator (pMLC) shaped electron fields is developed and investigated. METHODS: The Swiss Monte Carlo Plan (SMCP) allows the calculation of dose distributions for pMLC shaped electron beams. SMCP is interfaced with the Eclipse TPS (Varian Medical Systems, Palo Alto, CA) which can thus be included into the inverse treatment planning process for MERT. This process starts with the import of a CT-scan into Eclipse, the contouring of the target and the organs at risk (OARs), and the choice of the initial electron beam directions. For each electron beam, the number of apertures, their energy, and initial shape are defined. Furthermore, the DAO requires dose-volume constraints for the structures contoured. In order to carry out the DAO efficiently, the initial electron beams are divided into a grid of beamlets. For each of those, the dose distribution is precalculated using a modified electron beam model, resulting in a dose list for each beamlet and energy. Then the DAO is carried out, leading to a set of optimal apertures and corresponding weights. These optimal apertures are now converted into pMLC shaped segments and the dose calculation for each segment is performed. For these dose distributions, a weight optimization process is launched in order to minimize the differences between the dose distribution using the optimal apertures and the pMLC segments. Finally, a deliverable dose distribution for the MERT plan is obtained and loaded back into Eclipse for evaluation. For an idealized water phantom geometry, a MERT treatment plan is created and compared to the plan obtained using a previously developed forward planning strategy. Further, MERT treatment plans for three clinical situations (breast, chest wall, and parotid metastasis of a squamous cell skin carcinoma) are created using the developed inverse planning strategy. The MERT plans are compared to clinical standard treatment plans using photon beams and the differences between the optimal and the deliverable dose distributions are determined. RESULTS: For the idealized water phantom geometry, the inversely optimized MERT plan is able to obtain the same PTV coverage, but with an improved OAR sparing compared to the forwardly optimized plan. Regarding the right-sided breast case, the MERT plan is able to reduce the lung volume receiving more than 30% of the prescribed dose and the mean lung dose compared to the standard plan. However, the standard plan leads to a better homogeneity within the CTV. The results for the left-sided thorax wall are similar but also the dose to the heart is reduced comparing MERT to the standard treatment plan. For the parotid case, MERT leads to lower doses for almost all OARs but to a less homogeneous dose distribution for the PTV when compared to a standard plan. For all cases, the weight optimization successfully minimized the differences between the optimal and the deliverable dose distribution. CONCLUSIONS: A beamlet based DAO using multiple beam angles is implemented and successfully tested for an idealized water phantom geometry and clinical situations.

Forward treatment planning for modulated electron radiotherapy (MERT) employing Monte Carlo methods.

PURPOSE: This paper describes the development of a forward planning process for modulated electron radiotherapy (MERT). The approach is based on a previously developed electron beam model used to calculate dose distributions of electron beams shaped by a photon multi leaf collimator (pMLC). METHODS: As the electron beam model has already been implemented into the Swiss Monte Carlo Plan environment, the Eclipse treatment planning system (Varian Medical Systems, Palo Alto, CA) can be included in the planning process for MERT. In a first step, CT data are imported into Eclipse and a pMLC shaped electron beam is set up. This initial electron beam is then divided into segments, with the electron energy in each segment chosen according to the distal depth of the planning target volume (PTV) in beam direction. In order to improve the homogeneity of the dose distribution in the PTV, a feathering process (Gaussian edge feathering) is launched, which results in a number of feathered segments. For each of these segments a dose calculation is performed employing the in-house developed electron beam model along with the macro Monte Carlo dose calculation algorithm. Finally, an automated weight optimization of all segments is carried out and the total dose distribution is read back into Eclipse for display and evaluation. One academic and two clinical situations are investigated for possible benefits of MERT treatment compared to standard treatments performed in our clinics and treatment with a bolus electron conformal (BolusECT) method. RESULTS: The MERT treatment plan of the academic case was superior to the standard single segment electron treatment plan in terms of organs at risk (OAR) sparing. Further, a comparison between an unfeathered and a feathered MERT plan showed better PTV coverage and homogeneity for the feathered plan, with V95% increased from 90% to 96% and V107% decreased from 8% to nearly 0%. For a clinical breast boost irradiation, the MERT plan led to a similar homogeneity in the PTV compared to the standard treatment plan while the mean body dose was lower for the MERT plan. Regarding the second clinical case, a whole breast treatment, MERT resulted in a reduction of the lung volume receiving more than 45% of the prescribed dose when compared to the standard plan. On the other hand, the MERT plan leads to a larger low-dose lung volume and a degraded dose homogeneity in the PTV. For the clinical cases evaluated in this work, treatment plans using the BolusECT technique resulted in a more homogenous PTV and CTV coverage but higher doses to the OARs than the MERT plans. CONCLUSIONS: MERT treatments were successfully planned for phantom and clinical cases, applying a newly developed intuitive and efficient forward planning strategy that employs a MC based electron beam model for pMLC shaped electron beams. It is shown that MERT can lead to a dose reduction in OARs compared to other methods. The process of feathering MERT segments results in an improvement of the dose homogeneity in the PTV.

Monte Carlo based beam model using a photon MLC for modulated electron radiotherapy.

PURPOSE: Modulated electron radiotherapy (MERT) promises sparing of organs at risk for certain tumor sites. Any implementation of MERT treatment planning requires an accurate beam model. The aim of this work is the development of a beam model which reconstructs electron fields shaped using the Millennium photon multileaf collimator (MLC) (Varian Medical Systems, Inc., Palo Alto, CA) for a Varian linear accelerator (linac). METHODS: This beam model is divided into an analytical part (two photon and two electron sources) and a Monte Carlo (MC) transport through the MLC. For dose calculation purposes the beam model has been coupled with a macro MC dose calculation algorithm. The commissioning process requires a set of measurements and precalculated MC input. The beam model has been commissioned at a source to surface distance of 70 cm for a Clinac 23EX (Varian Medical Systems, Inc., Palo Alto, CA) and a TrueBeam linac (Varian Medical Systems, Inc., Palo Alto, CA). For validation purposes, measured and calculated depth dose curves and dose profiles are compared for four different MLC shaped electron fields and all available energies. Furthermore, a measured two-dimensional dose distribution for patched segments consisting of three 18 MeV segments, three 12 MeV segments, and a 9 MeV segment is compared with corresponding dose calculations. Finally, measured and calculated two-dimensional dose distributions are compared for a circular segment encompassed with a C-shaped segment. RESULTS: For 15 × 34, 5 × 5, and 2 × 2 cm(2) fields differences between water phantom measurements and calculations using the beam model coupled with the macro MC dose calculation algorithm are generally within 2% of the maximal dose value or 2 mm distance to agreement (DTA) for all electron beam energies. For a more complex MLC pattern, differences between measurements and calculations are generally within 3% of the maximal dose value or 3 mm DTA for all electron beam energies. For the two-dimensional dose comparisons, the differences between calculations and measurements are generally within 2% of the maximal dose value or 2 mm DTA. CONCLUSIONS: The results of the dose comparisons suggest that the developed beam model is suitable to accurately reconstruct photon MLC shaped electron beams for a Clinac 23EX and a TrueBeam linac. Hence, in future work the beam model will be utilized to investigate the possibilities of MERT using the photon MLC to shape electron beams.

SU-E-I-03: Scatter and Beam Hardening Correction for µCBCT Scanners Using Monte Carlo.

PURPOSE: Micro CBCT scanners have a broad spectrum of applications in medicine and material science. However, CBCT suffers from scatter radiation and spectral effects such as beam hardening (BH). In this work an iterative BH- and scatter-correction algorithm was developed using Monte Carlo (MC) methods. METHODS: Two µCBCT scanner models (XtremeCT and µCT100 from SCANCO Medical AG) were simulated using EGSnrc/EGS++. The scanner measures the attenuation of X-rays passing through the phantom and interacting in a scintillating detector. The MC method is used to characterize the influence of the scattering and BH- effects. In addition, an analytical model is developed in order to correct for the scattering effect. For this purpose, by using MC methods, different scatter components were analyzed with respect to the number of scatter interactions within the given geometry. For the BH-effect, after assessing the detector response for an equivalent mono-energetic and scatter-free system using MC methods, an analytical correction model was developed. Both correction methods were implemented as an iterative reconstruction correction algorithm and were tested for various phantoms. RESULTS: MC simulations show that the amounts of correction for the scattering and BH- effects are in the same order of magnitude. The correction term for scattering effects is a function of the scatter to primary ratio only and is mainly due to single scattered particles. The contribution of multiple scattered particles to the total scatter signal is small and can be approximated by a constant. In all cases tested, the reconstructed linear attenuation coefficients converge to the mono-energetic reference values after 2-3 iteration steps with a deviation of about 1%. CONCLUSIONS: By using an iterative correction algorithm using single scatter approximation, BH- and scatter correction can be performed accurately for µCBCT scanners. This work was supported by CTI-10629.1 and SCANCO Medical AG.

MO-F-BRB-02: Macro Monte Carlo for Proton Dose Calculation in Different Materials.

PURPOSE: Although the Monte Carlo (MC) method allows accurate dose calculation its usage is limited due to long computing time. In order to gain efficiency, a new macro MC (MMC) technique for proton dose calculations in homogeneous materials has been developed. METHODS: The macro MC is based on a local to global MC approach. The local simulations using Geant4 consist of mono-energetic proton pencil beams (10 to 250 MeV) impinging perpendicularly on slabs of different thicknesses (1-10 mm) and different materials (water, lung, muscle, adipose, bone). During the local simulation multiple scattering, ionization, elastic and inelastic interactions have been taken into account and the physical characteristics such as lateral displacement and energy loss have been scored for primary and secondary particles. The scored data from appropriate slabs is then used for the stepwise transport of the protons in the MMC simulation while calculating the energy loss along the path between entrance and exit position. Additionally, ions and neutrons are taken into account for the dose calculation. In order to validate the MMC, calculated dose distributions using the MMC transport and Geant4 have been compared for different mono-energetic proton pencil beams impinging on phantoms with different homogeneous materials. RESULTS: The agreement of calculated integral depth dose curves is better than 1% or 1 mm for all pencil beams and materials considered. For the dose profiles the agreement is within 1% or 1 mm for all energies, depths and materials. The efficiency of MMC is about 200 times higher than for Geant4. CONCLUSIONS: The dose comparisons demonstrate that the new MMC results in very accurate and efficient dose calculations for proton beams in homogeneous materials. In future, the MMC method will be extended to inhomogeneous situations in order to allow patient dose calculations for proton beams. This work was supported by Varian Medical Systems.

A virtual source model for kilo-voltage cone beam CT: source characteristics and model validation.

PURPOSE: The purpose of this investigation was to study the source characteristics of a clinical kilo-voltage cone beam CT unit and to develop and validate a virtual source model that could be used for treatment planning purposes. METHODS: We used a previously commissioned full Monte Carlo model and new bespoke software to study the source characteristics of a clinical kilo-voltage cone beam CT (CBCT) unit. We identified the main particle sources, their spatial, energy and angular distribution for all the image acquisition presets currently used in our clinical practice. This includes a combination of two energies (100 and 120 kVp), two filters (neutral and bowtie), and eight different x-ray beam apertures. We subsequently built a virtual source model which we validated against full Monte Carlo calculations. RESULTS: We found that the radiation output of the clinical kilo-voltage cone beam CT unit investigated in this study could be reproduced with a virtual model comprising of two sources (target and filtration cone) or three sources (target, filtration cone and bowtie filter) when additional filtration was used. With this model, we accounted for more than 97% of the photons exiting the unit. Each source in our model was characterised by a origin distribution in both X and Y directions, a fluence map, a single energy spectrum for unfiltered beams and a two dimensional energy spectrum for bowtie filtered beams. The percentage dose difference between full Monte Carlo and virtual source model based dose distributions was well within the statistical uncertainty associated with the calculations ( ± 2%, one standard deviation) in all cases studied. CONCLUSIONS: The virtual source that we developed is accurate in calculating the dose delivered from a commercial kilo-voltage cone beam CT unit operating with routine clinical image acquisition settings. Our data have also shown that target, filtration cone, and bowtie filter sources needed to be all included in the model in order to accurately replicate the dose distribution from the clinical radiation beam.

Comparison of monte carlo collimator transport methods for photon treatment planning in radiotherapy.

PURPOSE: The aim of this work was a Monte Carlo (MC) based investigation of the impact of different radiation transport methods in collimators of a linear accelerator on photon beam characteristics, dose distributions, and efficiency. Thereby it is investigated if it is possible to use different simplifications in the radiation transport for some clinical situations in order to save calculation time. METHODS: Within the Swiss Monte Carlo Plan, a GUI-based framework for photon MC treatment planning, different MC methods are available for the radiation transport through the collimators [secondary jaws and multileaf collimator (MLC)]: EGSnrc (reference), VMC++, and Pin (an in-house developed MC code). Additional nonfull transport methods were implemented in order to provide different complexity levels for the MC simulation: Considering collimator attenuation only, considering Compton scatter only or just the firstCompton process, and considering the collimators as totally absorbing. Furthermore, either a simple or an exact geometry of the collimators can be selected for the absorbing or attenuation method. Phasespaces directly above and dose distributions in a water phantom are analyzed for academic and clinical treatment fields using 6 and 15 MV beams, including intensity modulated radiation therapy with dynamic MLC. RESULTS: For all MC transport methods, differences in the radial mean energy and radial energy fluence are within 1% inside the geometric field. Below the collimators, the energy fluence is underestimated for nonfull MC transport methods ranging from 5% for Compton to 100% for Absorbing. Gamma analysis using EGSnrc calculated doses as reference shows that the percentage of voxels fulfilling a 1% /1 mm criterion is at least 98% when using VMC++, Compton, or firstCompton transport methods. When using the methods Pin, Transmission, Flat-Transmission, Flat-Absorbing or Absorbing, the mean value of points fulfilling this criterion over all tested cases is 97%, 88%, 74%, 68%, or 65%, respectively. However, compared to EGSnrc calculations, the gain in efficiency is a factor of up to 10 for VMC++ and up to 48 for the absorbing method. CONCLUSIONS: The results of this investigation suggest that it is an option to use a simple transport technique in the initial treatment planning process and use more accurate transport methods for the final dose calculation accepting longer calculation times.

The exchange of radiotherapy data as part of an electronic patient-referral system.

PURPOSE: To describe the implementation and use of an electronic patient-referral system as an aid to the efficient referral of patients to a remote and specialized treatment center. METHODS AND MATERIALS: A system for the exchange of radiotherapy data between different commercial planning systems and a specially developed planning system for proton therapy has been developed through the use of the PAPYRUS diagnostic image standard as an intermediate format. To ensure the cooperation of the different TPS manufacturers, the number of data sets defined for transfer has been restricted to the three core data sets of CT, VOIs, and three-dimensional dose distributions. As a complement to the exchange of data, network-wide application-sharing (video-conferencing) technologies have been adopted to provide methods for the interactive discussion and assessment of treatments plans with one or more partner clinics. RESULTS: Through the use of evaluation plans based on the exchanged data, referring clinics can accurately assess the advantages offered by proton therapy on a patient-by-patient basis, while the practicality or otherwise of the proposed treatments can simultaneously be assessed by the proton therapy center. Such a system, along with the interactive capabilities provided by video-conferencing methods, has been found to be an efficient solution to the problem of patient assessment and selection at a specialized treatment center, and is a necessary first step toward the full electronic integration of such centers with their remotely situated referral centers.