Beam Position Projection Algorithms in Proton Pencil Beam Scanning.

Beam position uncertainties along the beam trajectory arise from the accelerator, beamline, and scanning magnets (SMs). They can be monitored in real time, e.g., through strip ionization chambers (ICs), and treatments can be paused if needed. Delivery is more reliable and accurate if the beam position is projected from monitored nozzle parameters to the isocenter, allowing for accurate online corrections to be performed. Beam position projection algorithms are also used in post-delivery log file analyses. In this paper, we investigate the four potential algorithms that can be applied to all pencil beam scanning (PBS) nozzles. For some combinations of nozzle configurations and algorithms, however, the projection uses beam properties determined offline (e.g., through beam tuning or technical commissioning). The best algorithm minimizes either the total uncertainty (i.e., offline and online) or the total offline uncertainty in the projection. Four beam position algorithms are analyzed (A1-A4). Two nozzle lengths are used as examples: a large nozzle (1.5 m length) and a small nozzle (0.4 m length). Three nozzle configurations are considered: IC after SM, IC before SM, and ICs on both sides. Default uncertainties are selected for ion chamber measurements, nozzle entrance beam position and angle, and scanning magnet angle. The results for other uncertainties can be determined by scaling these results or repeating the error propagation. We show the propagation of errors from two locations and the SM angle to the isocenter for all the algorithms. The best choice of algorithm depends on the nozzle length and is A1 and A3 for the large and small nozzles, respectively. If the total offline uncertainty is to be minimized (a better choice if the offline uncertainty is not stable), the best choice of algorithm changes to A1 for the small nozzle for some hardware configurations. Reducing the nozzle length can help to reduce the gantry size and make proton therapy more accessible. This work is important for designing smaller nozzles and, consequently, smaller gantries. This work is also important for log file analyses.

Late Radiation-Induced Carotid Artery Stenosis and Stroke in Pediatric Patient Treated With Proton Radiation Therapy for Skull-Base Chordoma.

Radiation vasculopathy is a well-recognized late complication of radiation therapy. We present a case of a stroke 29 years after high-dose proton radiation therapy for skull-base chordoma due to occlusion of bilateral internal carotid arteries.

Pencil beam scanning dose calibration at reduced source-to-axis distance.

BACKGROUND: Pencil beam scanning (PBS) monitoring chambers use an ionization control signal, monitor units (MUs), or gigaprotons (Gp) to irradiate a pencil beam and normalize dose calculations. The nozzle deflects the beam from the nozzle axis by an angle subtended at the source-to-axis distance (τ) from the isocenter. If the angle is not correctly considered in calibrations or calculations, it can lead to systematic errors. PURPOSE: Aspects to consider for machines of various τs are fourfold. First, for the machine, there is a pathlength change of proton tracks in the monitor chamber. Second, for measurements, a uniform-square irradiation over a plane, with constant Gp per spot, does not deliver uniform dose in a measurement plane. Third, for Monte Carlo (MC) simulations, Gp (and not MU) is proportional to simulating a number of protons. Fourth, for pencil beam algorithms (PBA), MU or Gp may be used for pencil beam weight, but usage needs to be consistent. Another consideration is the beam shape change from circular to oval in the projection onto voxels. METHODS: Coordinate systems for PBS delivery are described. RESULTS: Users of intermediate-τ machines, corresponding to the onset of 1% pathlength corrections within the scanned field size, must not assume that MUs are proportional to the number of particles in MC simulations, and the PBA may need pathlength corrections. For a field size of 24 × 24 cm2 , intermediate-τ machines correspond to 59 cm ≤ τ < 120 cm. For a field size of 40 × 40 cm2 , intermediate-τ machines correspond to 98 cm ≤ τ < 200 cm. Small-τ machines correspond to τ < 59 and 98 cm at these field sizes, respectively, which require corrections in projecting the beam shape onto voxels. CONCLUSIONS: Identifying corrections due to the pencil beam angle and their onset are important for reducing the outer diameter of proton therapy gantries. The use of Gp (or the number of protons) meterset standardizes data interchange and helps to reduce systematic errors due to the angle of the beam.

DICOM-RT Ion interface to utilize MC simulations in routine clinical workflow for proton pencil beam radiotherapy.

To adopt Monte Carlo (MC) simulations as an independent dose calculation method for proton pencil beam radiotherapy, an interface that converts the plan information in DICOM format into MC components such as geometries and beam source is a crucial element. For this purpose, a DICOM-RT Ion interface (https://github.com/topasmc/dicom-interface) has been developed and integrated into the TOPAS MC code to perform such conversions on-the-fly. DICOM-RT objects utilized in this interface include Ion Plan (RTIP), Ion Beams Treatment Record (RTIBTR), CT image, and Dose. Beamline geometries, gantry and patient coordinate systems, and fluence maps are determined from RTIP and/or RTIBTR. In this interface, DICOM information is processed and delivered to a MC engine in two steps. A MC model, which consists of beamline geometries and beam source, to represent a treatment machine is created by a DICOM parser of the interface. The complexities from different DICOM types, various beamline configurations and source models are handled in this step. Next, geometry information and beam source are transferred to TOPAS on-the-fly via the developed TOPAS extensions. This interface with two treatment machines was successfully deployed into our automated MC workflow which provides simulated dose and LET distributions in a patient or a water phantom automatically when a new plan is identified. The developed interface provides novel features such as handling multiple treatment systems based on different DICOM types, DICOM conversions on-the-fly, and flexible sampling methods that significantly reduce the burden of handling DICOM based plan or treatment record information for MC simulations.

MRI-based IMPT planning for prostate cancer.

PURPOSE: Treatment planning for proton therapy requires the relative proton stopping power ratio (RSP) information of the patient for accurate dose calculations. RSP are conventionally obtained after mapping of the Hounsfield units (HU) from a calibrated patient computed tomography (CT). One or multiple CT are needed for a given treatment which represents additional, undesired dose to the patient. For prostate cancer, magnetic resonance imaging (MRI) scans are the gold standard for segmentation while offering dose-less imaging. We here quantify the clinical applicability of converted MR images as a substitute for intensity modulated proton therapy (IMPT) treatment of the prostate. METHODS: MRCAT (Magnetic Resonance for Calculating ATtenuation) is a Philips-developed technology which produces a synthetic CT image consisting of five HU from a specific set of MRI acquisitions. MRCAT and original planning CT data sets were obtained for ten patients. An IMPT plan was generated on the MRCAT for each patient. Plans were produced such that they fulfill the prostate protocol in use at Massachusetts General Hospital (MGH). The plans were then recomputed onto the nominal planning CT for each patient. Robustness analyses (±5 mm setup shifts and ±3.5 % range uncertainties) were also performed. RESULTS: Comparison of MRCAT plans and their recomputation onto the planning CT plan showed excellent agreement. Likewise, dose perturbations due to setup shifts and range uncertainties were well within clinical acceptance demonstrating the clinical viability of the approach. CONCLUSIONS: This work demonstrate the clinical acceptability of substituting MR converted RSP images instead of CT for IMPT planning of prostate cancer. This further translates into higher contouring accuracy along with lesser imaging dose.

Design of a QA method to characterize submillimeter-sized PBS beam properties using a 2D ionization chamber array.

Pencil beam scanning (PBS) periodic quality assurance (QA) programs ensure the beam delivered to patients is within technical specifications. Two critical specifications for PBS delivery are the beam width and position. The aim of this study is to investigate whether a 2D ionization chamber array, such as the MatriXX detector (IBA Dosimetry, Schwarzenbruck, Germany), can be used to characterize submillimeter-sized PBS beam properties. The motivation is to use standard equipment, which may have pixel spacing coarser than the pencil beam size, and simplify QA workflow. The MatriXX pixels are cylindrical in shape with 4.5 mm diameter and are spaced 7.62 mm from center to center. Two major effects limit the ability of using the MatriXX to measure the spot position and width accurately. The first effect is that too few pixels sample the Gaussian shaped pencil beam profile and the second effect is volume averaging of the Gaussian profile over the pixel sensitive volumes. We designed a method that overcomes both limitations and hence enables the use of the MatriXX to characterize sub-millimeter-sized PBS beam properties. This method uses a cross-like irradiation pattern that is designed to increase the number of sampling data points and a modified Gaussian fitting technique to correct for volume averaging effects. Detector signals were calculated in this study and random noise and setup errors were added to simulate measured data. With the techniques developed in this work, the MatriXX detector can be used to characterize the position and width of sub-millimeter, σ = 0.7 mm, sized pencil beams with uncertainty better than 3% relative to σ. With the irradiation only covering 60% of the MatriXX, the position and width of σ = 0.9 mm sized pencil beams can be determined with uncertainty better than 3% relative to σ. If one were to not use a cross-like irradiation pattern, then the position and width of σ = 3.6 mm sized pencil beams can be determined with uncertainty better than 3% relative to σ. If one were to not use a cross-like pattern nor volume averaging corrections, then the position and width of σ = 5.0 mm sized pencil beams can be determined with uncertainty better than 3% relative to σ. This work helps to simplify periodic QA in proton therapy because more routinely used ionization chamber arrays can be used to characterize narrow pencil beam properties.

National Cancer Institute Workshop on Proton Therapy for Children: Considerations Regarding Brainstem Injury.

PURPOSE: Proton therapy can allow for superior avoidance of normal tissues. A widespread consensus has been reached that proton therapy should be used for patients with curable pediatric brain tumor to avoid critical central nervous system structures. Brainstem necrosis is a potentially devastating, but rare, complication of radiation. Recent reports of brainstem necrosis after proton therapy have raised concerns over the potential biological differences among radiation modalities. We have summarized findings from the National Cancer Institute Workshop on Proton Therapy for Children convened in May 2016 to examine brainstem injury. METHODS AND MATERIALS: Twenty-seven physicians, physicists, and researchers from 17 institutions with expertise met to discuss this issue. The definition of brainstem injury, imaging of this entity, clinical experience with photons and photons, and potential biological differences among these radiation modalities were thoroughly discussed and reviewed. The 3 largest US pediatric proton therapy centers collectively summarized the incidence of symptomatic brainstem injury and physics details (planning, dosimetry, delivery) for 671 children with focal posterior fossa tumors treated with protons from 2006 to 2016. RESULTS: The average rate of symptomatic brainstem toxicity from the 3 largest US pediatric proton centers was 2.38%. The actuarial rate of grade ≥2 brainstem toxicity was successfully reduced from 12.7% to 0% at 1 center after adopting modified radiation guidelines. Guidelines for treatment planning and current consensus brainstem constraints for proton therapy are presented. The current knowledge regarding linear energy transfer (LET) and its relationship to relative biological effectiveness (RBE) are defined. We review the current state of LET-based planning. CONCLUSIONS: Brainstem injury is a rare complication of radiation therapy for both photons and protons. Substantial dosimetric data have been collected for brainstem injury after proton therapy, and established guidelines to allow for safe delivery of proton radiation have been defined. Increased capability exists to incorporate LET optimization; however, further research is needed to fully explore the capabilities of LET- and RBE-based planning.

Brainstem Injury in Pediatric Patients With Posterior Fossa Tumors Treated With Proton Beam Therapy and Associated Dosimetric Factors.

PURPOSE: Proton radiation therapy is commonly used in young children with brain tumors for its potential to reduce late effects. However, some proton series report higher rates of brainstem injury (0%-16%) than most photon series (2.2%-8.6%). We report the incidence of brainstem injury and a risk factor analysis in pediatric patients with posterior fossa primary tumors treated with proton radiation therapy at our institution. METHODS AND MATERIALS: The study included 216 consecutive patients treated between 2000 and 2015. Dosimetry was available for all but 4 patients. Grade 2 to 5 late brainstem toxicity was assessed by the National Cancer Institute Common Terminology Criteria for Adverse Events version 4.0. RESULTS: The histologies include medulloblastoma (n=154, 71.3%), ependymoma (n=56, 25.9%), and atypical teratoid rhabdoid tumor (n=6, 2.8%). The median age at irradiation was 6.6 years (range, 0.5-23.1 years); median dose, 54 gray relative biological effectiveness (Gy RBE) (range, 46.8-59.4 Gy RBE); and median follow-up period, 4.2 years (range, 0.1-15.3 years) among 198 survivors. Of the patients, 83.3% received chemotherapy; 70.4% achieved gross total resection. The crude rate of injury was 2.3% in all patients, 1.9% in those with medulloblastoma, 3.6% in those with ependymoma, and 0% in those with atypical teratoid rhabdoid tumor. The 5-year cumulative incidence of injury was 2.0% (95% confidence interval, 0.7%-4.8%). The median brainstem dose (minimum dose received by 50% of brainstem) in the whole cohort was 53.6 Gy RBE (range, 16.5-56.8 Gy RBE); maximum point dose within the brainstem (Dmax), 55.2 Gy RBE (range, 48.4-60.5 Gy RBE); and mean dose, 50.4 Gy RBE (range, 21.1-56.7 Gy RBE). In the 5 patients with injury, the median minimum dose received by 50% of the brainstem was 54.6 Gy RBE (range, 50.2-55.1 Gy RBE); Dmax, 56.2 Gy RBE (range, 55.0-57.1 Gy RBE); mean dose, 51.3 Gy RBE (range, 45.4-54.4 Gy RBE); and median volume of the brainstem receiving ≥55 Gy RBE (V55), 27.4% (range, 0%-59.4%). Of the 5 patients with injury, 4 had a brainstem Dmax in the highest quartile (≥55.8 Gy RBE, P = .016) and a V55 in the highest tertile (>6.0%) of the cohort distribution (P = .047). Of the 5 patients with injury, 3 were aged >6 years (age range, 4.1-22.8 years), and 4 of 5 patients received chemotherapy and achieved gross total resection. CONCLUSIONS: The incidence of injury in pediatric patients with posterior fossa tumors is consistent with previous reports in the photon setting. Our data suggest that when Dmax and V55 are kept <55.8 Gy RBE and ≤6.0%, respectively, the 5-year rate of radiation brainstem injury would be <2%.

Effects of spot parameters in pencil beam scanning treatment planning.

BACKGROUND: Spot size σ (in air at isocenter), interspot spacing d, and spot charge q influence dose delivery efficiency and plan quality in Intensity Modulated Proton Therapy (IMPT) treatment planning. The choice and range of parameters varies among different manufacturers. The goal of this work is to demonstrate the influence of the spot parameters on dose quality and delivery in IMPT treatment plans, to show their interdependence, and to make practitioners aware of the spot parameter values for a certain facility. Our study could help as a guideline to make the trade-off between treatment quality and time in existing PBS centers and in future systems. METHODS: We created plans for seven patients and a phantom, with different tumor sites and volumes, and compared the effect of small-, medium-, and large-spot widths (σ = 2.5, 5, and 10 mm) and interspot distances (1σ, 1.5σ, and 1.75σ) on dose, spot charge, and treatment time. Moreover, we quantified how postplanning charge threshold cuts affect plan quality and the total number of spots to deliver, for different spot widths and interspot distances. We show the effect of a minimum charge (or MU) cutoff value for a given proton delivery system. RESULTS: Spot size had a strong influence on dose: larger spots resulted in more protons delivered outside the target region. We observed dose differences of 2-13 Gy (RBE) between 2.5 mm and 10 mm spots, where the amount of extra dose was due to dose penumbra around the target region. Interspot distance had little influence on dose quality for our patient group. Both parameters strongly influence spot charge in the plans and thus the possible impact of postplanning charge threshold cuts. If such charge thresholds are not included in the treatment planning system (TPS), it is important that the practitioner validates that a given combination of lower charge threshold, interspot spacing, and spot size does not result in a plan degradation. Low average spot charge occurs for small spots, small interspot distances, many beam directions, and low fractional dose values. CONCLUSIONS: The choice of spot parameters values is a trade-off between accelerator and beam line design, plan quality, and treatment efficiency. We recommend the use of small spot sizes for better organ-at-risk sparing and lateral interspot distances of 1.5σ to avoid long treatment times. We note that plan quality is influenced by the charge cutoff. Our results show that the charge cutoff can be sufficiently large (i.e., 106 protons) to accommodate limitations on beam delivery systems. It is, therefore, not necessary per se to include the charge cutoff in the treatment planning optimization such that Pareto navigation (e.g., as practiced at our institution) is not excluded and optimal plans can be obtained without, perhaps, a bias from the charge cutoff. We recommend that the impact of a minimum charge cut impact is carefully verified for the spot sizes and spot distances applied or that it is accommodated in the TPS.

Characterization of proton pencil beam scanning and passive beam using a high spatial resolution solid-state microdosimeter.

PURPOSE: This work aims to characterize a proton pencil beam scanning (PBS) and passive double scattering (DS) systems as well as to measure parameters relevant to the relative biological effectiveness (RBE) of the beam using a silicon on insulator (SOI) microdosimeter with well-defined 3D sensitive volumes (SV). The dose equivalent downstream and laterally outside of a clinical PBS treatment field was assessed and compared to that of a DS beam. METHODS: A novel silicon microdosimeter with well-defined 3D SVs was used in this study. It was connected to low noise electronics, allowing for detection of lineal energies as low as 0.15 keV/µm. The microdosimeter was placed at various depths in a water phantom along the central axis of the proton beam, and at the distal part of the spread-out Bragg peak (SOBP) in 0.5 mm increments. The RBE values of the pristine Bragg peak (BP) and SOBP were derived using the measured microdosimetric lineal energy spectra as inputs to the modified microdosimetric kinetic model (MKM). Geant4 simulations were performed in order to verify the calculated depth-dose distribution from the treatment planning system (TPS) and to compare the simulated dose-mean lineal energy to the experimental results. RESULTS: For a 131 MeV PBS spot (124.6 mm R90 range in water), the measured dose-mean lineal energy yD¯ increased from 2 keV/µm at the entrance to 8 keV/µm in the BP, with a maximum value of 10 keV/µm at the distal edge. The derived RBE distribution for the PBS beam slowly increased from 0.97 ± 0.14 at the entrance to 1.04 ± 0.09 proximal to the BP, then to 1.1 ± 0.08 in the BP, and steeply rose to 1.57 ± 0.19 at the distal part of the BP. The RBE distribution for the DS SOBP beam was approximately 0.96 ± 0.16 to 1.01 ± 0.16 at shallow depths, and 1.01 ± 0.16 to 1.28 ± 0.17 within the SOBP. The RBE significantly increased from 1.29 ± 0.17 to 1.43 ± 0.18 at the distal edge of the SOBP. CONCLUSIONS: The SOI microdosimeter with its well-defined 3D SV has applicability in characterizing proton radiation fields and can measure relevant physical parameters to model the RBE with submillimeter spatial resolution. It has been shown that for a physical dose of 1.82 Gy at the BP, the derived RBE based on the MKM model increased from 1.14 to 1.6 in the BP and its distal part. Good agreement was observed between the experimental and simulation results, confirming the potential application of SOI microdosimeter with 3D SV for quality assurance in proton therapy.

Evaluating Intensity Modulated Proton Therapy Relative to Passive Scattering Proton Therapy for Increased Vertebral Column Sparing in Craniospinal Irradiation in Growing Pediatric Patients.

PURPOSE: At present, proton craniospinal irradiation (CSI) for growing children is delivered to the whole vertebral body (WVB) to avoid asymmetric growth. We aimed to demonstrate the feasibility and potential clinical benefit of delivering vertebral body sparing (VBS) versus WVB CSI with passively scattered (PS) and intensity modulated proton therapy (IMPT) in growing children treated for medulloblastoma. METHODS AND MATERIALS: Five plans were generated for medulloblastoma patients, who had been previously treated with CSI PS proton radiation therapy: (1) single posteroanterior (PA) PS field covering the WVB (PS-PA-WVB); (2) single PA PS field that included only the thecal sac in the target volume (PS-PA-VBS); (3) single PA IMPT field covering the WVB (IMPT-PA-WVB); (4) single PA IMPT field, target volume including thecal sac only (IMPT-PA-VBS); and (5) 2 posterior-oblique (-35°, +35°) IMPT fields, with the target volume including the thecal sac only (IMPT2F-VBS). For all cases, 23.4 Gy (relative biologic effectiveness [RBE]) was prescribed to 95% of the spinal canal. The dose, linear energy transfer, and variable-RBE-weighted dose distributions were calculated for all plans using the tool for particle simulation, version 2, Monte Carlo system. RESULTS: IMPT VBS techniques efficiently spared the anterior vertebral bodies (AVBs), even when accounting for potential higher variable RBE predicted by linear energy transfer distributions. Assuming an RBE of 1.1, the V10 Gy(RBE) decreased from 100% for the WVB techniques to 59.5% to 76.8% for the cervical, 29.9% to 34.6% for the thoracic, and 20.6% to 25.1% for the lumbar AVBs, and the V20 Gy(RBE) decreased from 99.0% to 17.8% to 20.0% for the cervical, 7.2% to 7.6% for the thoracic, and 4.0% to 4.6% for the lumbar AVBs when IMPT VBS techniques were applied. The corresponding percentages for the PS VBS technique were higher. CONCLUSIONS: Advanced proton techniques can sufficiently reduce the dose to the vertebral body and allow for vertebral column growth for children with central nervous system tumors requiring CSI. This was true even when considering variable RBE values. A clinical trial is planned for VBS to the thoracic and lumbosacral spine in growing children.

Impact of spot charge inaccuracies in IMPT treatments.

BACKGROUND: Spot charge is one parameter of pencil-beam scanning dose delivery system whose accuracy is typically high but whose required value has not been investigated. In this work we quantify the dose impact of spot charge inaccuracies on the dose distribution in patients. Knowing the effect of charge errors is relevant for conventional proton machines, as well as for new generation proton machines, where ensuring accurate charge may be challenging. METHODS: Through perturbation of spot charge in treatment plans for seven patients and a phantom, we evaluated the dose impact of absolute (up to 5× 106 protons) and relative (up to 30%) charge errors. We investigated the dependence on beam width by studying scenarios with small, medium and large beam sizes. Treatment plan statistics included the Γ passing rate, dose-volume-histograms and dose differences. RESULTS: The allowable absolute charge error for small spot plans was about 2× 106 protons. Larger limits would be allowed if larger spots were used. For relative errors, the maximum allowable error size for small, medium and large spots was about 13%, 8% and 6% for small, medium and large spots, respectively. CONCLUSIONS: Dose distributions turned out to be surprisingly robust against random spot charge perturbation. Our study suggests that ensuring spot charge errors as small as 1-2% as is commonly aimed at in conventional proton therapy machines, is clinically not strictly needed.

A Greedy reassignment algorithm for the PBS minimum monitor unit constraint.

Proton pencil beam scanning (PBS) treatment plans are made of numerous unique spots of different weights. These weights are optimized by the treatment planning systems, and sometimes fall below the deliverable threshold set by the treatment delivery system. The purpose of this work is to investigate a Greedy reassignment algorithm to mitigate the effects of these low weight pencil beams. The algorithm is applied during post-processing to the optimized plan to generate deliverable plans for the treatment delivery system. The Greedy reassignment method developed in this work deletes the smallest weight spot in the entire field and reassigns its weight to its nearest neighbor(s) and repeats until all spots are above the minimum monitor unit (MU) constraint. Its performance was evaluated using plans collected from 190 patients (496 fields) treated at our facility. The Greedy reassignment method was compared against two other post-processing methods. The evaluation criteria was the γ-index pass rate that compares the pre-processed and post-processed dose distributions. A planning metric was developed to predict the impact of post-processing on treatment plans for various treatment planning, machine, and dose tolerance parameters. For fields with a pass rate of 90 ± 1% the planning metric has a standard deviation equal to 18% of the centroid value showing that the planning metric and γ-index pass rate are correlated for the Greedy reassignment algorithm. Using a 3rd order polynomial fit to the data, the Greedy reassignment method has 1.8 times better planning metric at 90% pass rate compared to other post-processing methods. As the planning metric and pass rate are correlated, the planning metric could provide an aid for implementing parameters during treatment planning, or even during facility design, in order to yield acceptable pass rates. More facilities are starting to implement PBS and some have spot sizes (one standard deviation) smaller than 5 mm, hence would require small spot spacing. While this is not the only parameter that affects the optimized plan, the perturbation due to the minimum MU constraint increases with decreasing spot spacing. This work could help to design the minimum MU threshold with the goal to keep the γ-index pass rate above an acceptable value.

Impact of Spot Size and Beam-Shaping Devices on the Treatment Plan Quality for Pencil Beam Scanning Proton Therapy.

PURPOSE: This study aimed to assess the clinical impact of spot size and the addition of apertures and range compensators on the treatment quality of pencil beam scanning (PBS) proton therapy and to define when PBS could improve on passive scattering proton therapy (PSPT). METHODS AND MATERIALS: The patient cohort included 14 pediatric patients treated with PSPT. Six PBS plans were created and optimized for each patient using 3 spot sizes (∼12-, 5.4-, and 2.5-mm median sigma at isocenter for 90- to 230-MeV range) and adding apertures and compensators to plans with the 2 larger spots. Conformity and homogeneity indices, dose-volume histogram parameters, equivalent uniform dose (EUD), normal tissue complication probability (NTCP), and integral dose were quantified and compared with the respective PSPT plans. RESULTS: The results clearly indicated that PBS with the largest spots does not necessarily offer a dosimetric or clinical advantage over PSPT. With comparable target coverage, the mean dose (Dmean) to healthy organs was on average 6.3% larger than PSPT when using this spot size. However, adding apertures to plans with large spots improved the treatment quality by decreasing the average Dmean and EUD by up to 8.6% and 3.2% of the prescribed dose, respectively. Decreasing the spot size further improved all plans, lowering the average Dmean and EUD by up to 11.6% and 10.9% compared with PSPT, respectively, and eliminated the need for beam-shaping devices. The NTCP decreased with spot size and addition of apertures, with maximum reduction of 5.4% relative to PSPT. CONCLUSIONS: The added benefit of using PBS strongly depends on the delivery configurations. Facilities limited to large spot sizes (>∼8 mm median sigma at isocenter) are recommended to use apertures to reduce treatment-related toxicities, at least for complex and/or small tumors.

PBS machine interlocks using EWMA.

Delivery of pencil beam scanning (PBS) requires the on-line measurement of several beam parameters. If the measurement is outside of specified tolerances and a binary threshold algorithm is used, the beam will be paused. Given instrumentation and statistical noise such a system can lead to many pauses which could increase the treatment time. Statistical quality control methods are typically used on manufacturing lines to monitor a process and give early detection of a gradual problem and stop the process if a deviation is statistically significant. These methods can be used to develop a more intuitive algorithm for (PBS) delivery systems that is robust and safe and leads to decreased treatment times. The Exponentially Weighted Moving Average (EWMA) control scheme monitors deviations in beam properties which are averaged over a specified number of measurements with greater weight applied to the more recent ones. Simulation of an EWMA-style algorithm safely detected shifts in random and systematic delivery errors without false alarms. Binary and EWMA methods can be combined for improved reliability without sacrificing patient safety. In the EWMA method, the mean of a beam property can be related to systematic uncertainties and the standard deviation can be related to random uncertainties. This method allows one to have separate interlock levels for each type of uncertainty and to detect systematic trends.

Shortening delivery times of intensity modulated proton therapy by reducing proton energy layers during treatment plan optimization.

PURPOSE: To shorten delivery times of intensity modulated proton therapy by reducing the number of energy layers in the treatment plan. METHODS AND MATERIALS: We have developed an energy layer reduction method, which was implemented into our in-house-developed multicriteria treatment planning system "Erasmus-iCycle." The method consisted of 2 components: (1) minimizing the logarithm of the total spot weight per energy layer; and (2) iteratively excluding low-weighted energy layers. The method was benchmarked by comparing a robust "time-efficient plan" (with energy layer reduction) with a robust "standard clinical plan" (without energy layer reduction) for 5 oropharyngeal cases and 5 prostate cases. Both plans of each patient had equal robust plan quality, because the worst-case dose parameters of the standard clinical plan were used as dose constraints for the time-efficient plan. Worst-case robust optimization was performed, accounting for setup errors of 3 mm and range errors of 3% + 1 mm. We evaluated the number of energy layers and the expected delivery time per fraction, assuming 30 seconds per beam direction, 10 ms per spot, and 400 Giga-protons per minute. The energy switching time was varied from 0.1 to 5 seconds. RESULTS: The number of energy layers was on average reduced by 45% (range, 30%-56%) for the oropharyngeal cases and by 28% (range, 25%-32%) for the prostate cases. When assuming 1, 2, or 5 seconds energy switching time, the average delivery time was shortened from 3.9 to 3.0 minutes (25%), 6.0 to 4.2 minutes (32%), or 12.3 to 7.7 minutes (38%) for the oropharyngeal cases, and from 3.4 to 2.9 minutes (16%), 5.2 to 4.2 minutes (20%), or 10.6 to 8.0 minutes (24%) for the prostate cases. CONCLUSIONS: Delivery times of intensity modulated proton therapy can be reduced substantially without compromising robust plan quality. Shorter delivery times are likely to reduce treatment uncertainties and costs.

A novel approach to postmastectomy radiation therapy using scanned proton beams.

PURPOSE: Postmastectomy radiation therapy (PMRT), currently offered at Massachusetts General Hospital, uses proton pencil beam scanning (PBS) with intensity modulation, achieving complete target coverage of the chest wall and all nodal regions and reduced dose to the cardiac structures. This work presents the current methodology for such treatment and the ongoing effort for its improvements. METHODS AND MATERIALS: A single PBS field is optimized to ensure appropriate target coverage and heart/lung sparing, using an in-house-developed proton planning system with the capability of multicriteria optimization. The dose to the chest wall skin is controlled as a separate objective in the optimization. Surface imaging is used for setup because it is a suitable surrogate for superficial target volumes. In order to minimize the effect of beam range uncertainties, the relative proton stopping power ratio of the material in breast implants was determined through separate measurements. Phantom measurements were also made to validate the accuracy of skin dose calculation in the treatment planning system. Additionally, the treatment planning robustness was evaluated relative to setup perturbations and patient breathing motion. RESULTS: PBS PMRT planning resulted in appropriate target coverage and organ sparing, comparable to treatments by passive scattering (PS) beams but much improved in nodal coverage and cardiac sparing compared to conventional treatments by photon/electron beams. The overall treatment time was much shorter than PS and also shorter than conventional photon/electron treatment. The accuracy of the skin dose calculation by the planning system was within ±2%. The treatment was shown to be adequately robust relative to both setup uncertainties and patient breathing motion, resulting in clinically satisfying dose distributions. CONCLUSIONS: More than 25 PMRT patients have been successfully treated at Massachusetts General Hospital by using single-PBS fields. The methodology and robustness of both the setup and the treatment have been discussed.

Dose uncertainties in IMPT for oropharyngeal cancer in the presence of anatomical, range, and setup errors.

PURPOSE: Setup, range, and anatomical uncertainties influence the dose delivered with intensity modulated proton therapy (IMPT), but clinical quantification of these errors for oropharyngeal cancer is lacking. We quantified these factors and investigated treatment fidelity, that is, robustness, as influenced by adaptive planning and by applying more beam directions. METHODS AND MATERIALS: We used an in-house treatment planning system with multicriteria optimization of pencil beam energies, directions, and weights to create treatment plans for 3-, 5-, and 7-beam directions for 10 oropharyngeal cancer patients. The dose prescription was a simultaneously integrated boost scheme, prescribing 66 Gy to primary tumor and positive neck levels (clinical target volume-66 Gy; CTV-66 Gy) and 54 Gy to elective neck levels (CTV-54 Gy). Doses were recalculated in 3700 simulations of setup, range, and anatomical uncertainties. Repeat computed tomography (CT) scans were used to evaluate an adaptive planning strategy using nonrigid registration for dose accumulation. RESULTS: For the recalculated 3-beam plans including all treatment uncertainty sources, only 69% (CTV-66 Gy) and 88% (CTV-54 Gy) of the simulations had a dose received by 98% of the target volume (D98%) >95% of the prescription dose. Doses to organs at risk (OARs) showed considerable spread around planned values. Causes for major deviations were mixed. Adaptive planning based on repeat imaging positively affected dose delivery accuracy: in the presence of the other errors, percentages of treatments with D98% >95% increased to 96% (CTV-66 Gy) and 100% (CTV-54 Gy). Plans with more beam directions were not more robust. CONCLUSIONS: For oropharyngeal cancer patients, treatment uncertainties can result in significant differences between planned and delivered IMPT doses. Given the mixed causes for major deviations, we advise repeat diagnostic CT scans during treatment, recalculation of the dose, and if required, adaptive planning to improve adequate IMPT dose delivery.

Increasing maximum tumor dose to manage range uncertainties in IMPT treatment planning.

The accuracy of intensity modulated proton therapy (IMPT) is sensitive to range uncertainties. Geometric margins, as dosimetric surrogates, are ineffective and robust optimization strategies are needed. These, however, lead to increased normal tissue dose. We explore here how this dose increase can be reduced by increasing the maximum tumor dose instead. We focus on range uncertainties, modeled by scaling the stopping powers 5% up (undershoot) or down (overshoot) compared to the nominal scenario. Robust optimization optimizes for target dose conformity in the most likely scenario, not the worst, while constraining target coverage for the worst-case scenario. Non-robust plans are also generated. Different maximum target doses are applied (105% versus 120% versus 140%) to investigate the effect on normal tissue dose reduction. The method is tested on a homogeneous and a lung phantom and on a liver patient. Target D99 of the robust plans equals the prescription dose of 60 GyE for all scenarios, but decreases to 36 GyE for the non-robust plans. The mean normal tissue dose in a 2 cm ring around the target is 11% to 31% higher for the robust plans. This increase can be reduced to -8% and 3% (compared to the non-robust plan) by allowing a maximum tumor dose of 120% instead of 105%. Thus robustness leads to more normal tissue dose, but it can be compensated by allowing a higher maximum tumor dose.

Proton arc reduces range uncertainty effects and improves conformality compared with photon volumetric modulated arc therapy in stereotactic body radiation therapy for non-small cell lung cancer.

PURPOSE: To describe, in a setting of non-small cell lung cancer (NSCLC), the theoretical dosimetric advantages of proton arc stereotactic body radiation therapy (SBRT) in which the beam penumbra of a rotating beam is used to reduce the impact of range uncertainties. METHODS AND MATERIALS: Thirteen patients with early-stage NSCLC treated with proton SBRT underwent repeat planning with photon volumetric modulated arc therapy (Photon-VMAT) and an in-house-developed arc planning approach for both proton passive scattering (Passive-Arc) and intensity modulated proton therapy (IMPT-Arc). An arc was mimicked with a series of beams placed at 10° increments. Tumor and organ at risk doses were compared in the context of high- and low-dose regions, represented by volumes receiving >50% and <50% of the prescription dose, respectively. RESULTS: In the high-dose region, conformality index values are 2.56, 1.91, 1.31, and 1.74, and homogeneity index values are 1.29, 1.22, 1.52, and 1.18, respectively, for 3 proton passive scattered beams, Passive-Arc, IMPT-Arc, and Photon-VMAT. Therefore, proton arc leads to a 30% reduction in the 95% isodose line volume to 3-beam proton plan, sparing surrounding organs, such as lung and chest wall. For chest wall, V30 is reduced from 21 cm(3) (3 proton beams) to 11.5 cm(3), 12.9 cm(3), and 8.63 cm(3) (P=.005) for Passive-Arc, IMPT-Arc, and Photon-VMAT, respectively. In the low-dose region, the mean lung dose and V20 of the ipsilateral lung are 5.01 Gy(relative biological effectiveness [RBE]), 4.38 Gy(RBE), 4.91 Gy(RBE), and 5.99 Gy(RBE) and 9.5%, 7.5%, 9.0%, and 10.0%, respectively, for 3-beam, Passive-Arc, IMPT-Arc, and Photon-VMAT, respectively. CONCLUSIONS: Stereotactic body radiation therapy with proton arc and Photon-VMAT generate significantly more conformal high-dose volumes than standard proton SBRT, without loss of coverage of the tumor and with significant sparing of nearby organs, such as chest wall. In addition, both proton arc approaches spare the healthy lung from low-dose radiation relative to photon VMAT. Our data suggest that IMPT-Arc should be developed for clinical use.