Sensitivity analysis of Monte Carlo model of a gantry-mounted passively scattered proton system

Purpose

This study aimed to present guidance on the correlation between treatment nozzle and proton source parameters, and dose distribution of a passive double scattering compact proton therapy unit, known as Mevion S250.

Methods

All 24 beam options were modeled using the MCNPX MC code. The calculated physical dose for pristine peak, profiles, and spread out Bragg peak (SOBP) were benchmarked with the measured data. Track‐averaged LET (LET_t) and dose‐averaged LET (LET_d) distributions were also calculated. For the sensitivity investigations, proton beam line parameters including Average Energy (AE), Energy Spread (ES), Spot Size (SS), Beam Angle (BA), Beam Offset (OA), and Second scatter Offset (SO) from central Axis, and also First Scatter (FS) thickness were simulated in different stages to obtain the uncertainty of the derived results on the physical dose and LET distribution in a water phantom.

Results

For the physical dose distribution, the MCNPX MC model matched measurements data for all the options to within 2 mm and 2% criterion. The Mevion S250 was found to have a LET_t between 0.46 and 8.76 keV.μm^–1 and a corresponding LET_d between 0.84 and 15.91 keV.μm^–1. For all the options, the AE and ES had the greatest effect on the resulting depth of pristine peak and peak‐to‐plateau ratio respectively. BA, OA, and SO significantly decreased the flatness and symmetry of the profiles. The LETs were found to be sensitive to the AE, ES, and SS, especially in the peak region.

Conclusions

This study revealed the importance of considering detailed beam parameters, and identifying those that resulted in large effects on the physical dose distribution and LETs for a compact proton therapy machine.

TOPAS Simulation of the Mevion S250 compact proton therapy unit

As proton therapy becomes increasingly popular, so does the need for Monte Carlo simulation studies involving accurate beam line modeling of proton treatment units. In this study, the 24 beam configurations of the Mevion S250 proton therapy system installed recently at our institution were modeled using the TOolkit for PArticle Simulation (TOPAS) code. Pristine Bragg peak, spread out Bragg peak (SOBP), and lateral beam profile dose distributions were simulated and matched to the measurements taken during commissioning of the unit. Differences in the range for all Percent Depth Dose (PDD) curves between measured and simulated data agreed to within 0.1 cm. For SOBP scans, the SOBP widths all agreed to within 0.3 cm. With regards to lateral beam profile comparisons between the measured and simulated data, the penumbras differed by less than 1 mm and the flatness differed by less than 1% in nearly all cases. This study shows that Monte Carlo simulation studies involving the Mevion S250 proton therapy unit can be a viable tool in commissioning and verification of the proton treatment planning system.

The physics of proton therapy

The physics of proton therapy has advanced considerably since it was proposed in 1946. Today analytical equations and numerical simulation methods are available to predict and characterize many aspects of proton therapy. This article reviews the basic aspects of the physics of proton therapy, including proton interaction mechanisms, proton transport calculations, the determination of dose from therapeutic and stray radiations, and shielding design. The article discusses underlying processes as well as selected practical experimental and theoretical methods. We conclude by briefly speculating on possible future areas of research of relevance to the physics of proton therapy.

Comparison of MCNPX and Geant4 proton energy deposition predictions for clinical use

Several different Monte Carlo codes are currently being used at proton therapy centers to improve upon dose predictions over standard methods using analytical or semi-empirical dose algorithms. There is a need to better ascertain the differences between proton dose predictions from different available Monte Carlo codes. In this investigation Geant4 and MCNPX, the two most-utilized Monte Carlo codes for proton therapy applications, were used to predict energy deposition distributions in a variety of geometries, comprising simple water phantoms, water phantoms with complex inserts and in a voxelized geometry based on clinical CT data. The Gamma analysis was used to evaluate the differences of the predictions between the codes. The results show that in all the cases the agreement was better than clinical acceptance criteria.

Variations in proton scanned beam dose delivery due to uncertainties in magnetic beam steering

The purpose of this work was to develop a method to calculate and study the impact of fluctuations in the magnetic field strengths within the steering magnets in a proton scanning beam treatment nozzle on the dose delivered to the patient during a proton therapy treatment. First, an analytical relationship between magnetic field uncertainties in the steering magnets and the resulting lateral displacements in the position of the delivered scanned beam "dose spot" was established. Next, using a simple 3D dose calculation code and data from a validated Monte Carlo model of the proton scanning beam treatment nozzle, the uniform dose delivery to a 3D treatment volume was calculated. The dose distribution was then recalculated using the calculated lateral displacements due to magnetic field fluctuations to the proton pencil beam position. Using these two calculated dose distributions, the clinical effects of the magnetic field fluctuations were determined. A deliberate displacement of four adjacent spots either toward or away from each other was used to determine the "maximum" dose impact, while a random displacement of all spots was used to establish a more realistic clinical dose impact. Changes in the dose volume histogram (DVH) and the presence of hot and cold spots in the treatment volume were used to quantify the impact of dose-spot displacement. A general analytical relationship between magnetic field uncertainty and final dose-spot position is presented. This analytical relationship was developed such that it can be applied to study magnetic beam steering for any scanned beam nozzle design. Using this relationship the authors found for the example beam steering nozzle used in this study that deliberate lateral displacements of 0.5 mm or random lateral displacements of up to 1.0 mm produced a noticeable dose impact (5% hot spot) in the treatment volume. A noticeable impact (3% decrease in treatment volume coverage) on the DVH was observed for random displacements of up to 1.5 mm. For the scanning nozzle studied in this work, these displacement values correlated with an uncertainty value of 2.04% in the magnetic field values of the nozzle steering magnets. The authors conclude that fluctuations in the dose-spot delivery caused by uncertainty in the magnet fields used for beam steering could have clinically significant effects on the delivered dose distribution. Due to differences in the design and implementation of proton beam scanning nozzles at different treatment facilities, the effects of magnetic field fluctuations of dose delivery should be evaluated and understood for each specific nozzle design during clinical commissioning of the treatment nozzle.

Neutron production from beam-modifying devices in a modern double scattering proton therapy beam delivery system

In this work the neutron production in a passive beam delivery system was investigated. Secondary particles including neutrons are created as the proton beam interacts with beam shaping devices in the treatment head. Stray neutron exposure to the whole body may increase the risk that the patient develops a radiogenic cancer years or decades after radiotherapy. We simulated a passive proton beam delivery system with double scattering technology to determine the neutron production and energy distribution at 200 MeV proton energy. Specifically, we studied the neutron absorbed dose per therapeutic absorbed dose, the neutron absorbed dose per source particle and the neutron energy spectrum at various locations around the nozzle. We also investigated the neutron production along the nozzle's central axis. The absorbed doses and neutron spectra were simulated with the MCNPX Monte Carlo code. The simulations revealed that the range modulation wheel (RMW) is the most intense neutron source of any of the beam spreading devices within the nozzle. This finding suggests that it may be helpful to refine the design of the RMW assembly, e.g., by adding local shielding, to suppress neutron-induced damage to components in the nozzle and to reduce the shielding thickness of the treatment vault. The simulations also revealed that the neutron dose to the patient is predominated by neutrons produced in the field defining collimator assembly, located just upstream of the patient.

Prompt gamma-ray emission from biological tissues during proton irradiation: a preliminary study

In this paper, we present the results of a preliminary study of secondary 'prompt' gamma-ray emission produced by proton-nuclear interactions within tissue during proton radiotherapy. Monte Carlo simulations were performed for mono-energetic proton beams, ranging from 2.5 MeV to 250 MeV, irradiating elemental and tissue targets. Calculations of the emission spectra from different biological tissues and their elemental components were made. Also, prompt gamma rays emitted during delivery of a clinical proton spread-out Bragg peak (SOBP) in a homogeneous water phantom and a water phantom containing heterogeneous tissue inserts were calculated to study the correlation between prompt gamma-ray production and proton dose delivery. The results show that the prompt gamma-ray spectra differ significantly for each type of tissue studied. The relative intensity of the characteristic gamma rays emitted from a given tissue was shown to be proportional to the concentration of each element in that tissue. A strong correlation was found between the delivered SOBP dose distribution and the characteristic prompt gamma-ray production. Based on these results, we discuss the potential use of prompt gamma-ray emission as a method to verify the accuracy and efficacy of doses delivered with proton radiotherapy.

Assessment of the accuracy of an MCNPX-based Monte Carlo simulation model for predicting three-dimensional absorbed dose distributions

In recent years, the Monte Carlo method has been used in a large number of research studies in radiation therapy. For applications such as treatment planning, it is essential to validate the dosimetric accuracy of the Monte Carlo simulations in heterogeneous media. The AAPM Report no 105 addresses issues concerning clinical implementation of Monte Carlo based treatment planning for photon and electron beams, however for proton-therapy planning, such guidance is not yet available. Here we present the results of our validation of the Monte Carlo model of the double scattering system used at our Proton Therapy Center in Houston. In this study, we compared Monte Carlo simulated depth doses and lateral profiles to measured data for a magnitude of beam parameters. We varied simulated proton energies and widths of the spread-out Bragg peaks, and compared them to measurements obtained during the commissioning phase of the Proton Therapy Center in Houston. Of 191 simulated data sets, 189 agreed with measured data sets to within 3% of the maximum dose difference and within 3 mm of the maximum range or penumbra size difference. The two simulated data sets that did not agree with the measured data sets were in the distal falloff of the measured dose distribution, where large dose gradients potentially produce large differences on the basis of minute changes in the beam steering. Hence, the Monte Carlo models of medium- and large-size double scattering proton-therapy nozzles were valid for proton beams in the 100 MeV-250 MeV interval.