J Eley

Minibeam Therapy With Protons and Light Ions: Physical Feasibility and Potential to Reduce Radiation Side Effects and to Facilitate Hypofractionation

PURPOSE Despite several advantages of proton therapy over megavoltage x-ray therapy, its lack of proximal tissue sparing is a concern. The method presented here adds proximal tissue sparing to protons and light ions by turning their uniform incident beams into arrays of parallel, small, or thin (0.3-mm) pencil or planar minibeams, which are known to spare tissues. As these minibeams penetrate the tissues, they gradually broaden and merge with each other to produce a solid beam. METHODS AND MATERIALS Broadening of 0.3-mm-diameter, 109-MeV proton pencil minibeams was measured using a stack of radiochromic films with plastic spacers. Monte Carlo simulations were used to evaluate the broadening in water of minibeams of protons and several light ions and the dose from neutron generated by collimator. RESULTS A central parameter was tissue depth, where the beam full width at half maximum (FWHM) reached 0.7 mm, beyond which tissue sparing decreases. This depth was 22 mm for 109-MeV protons in a film stack. It was also found by simulations in water to be 23.5 mm for 109 MeV proton pencil minibeams and 26 mm for 116 MeV proton planar minibeams. For light ions, all with 10 cm range in water, that depth increased with particle size; specifically it was 51 mm for Li-7 ions. The ∼2.7% photon equivalent neutron skin dose from the collimator was reduced 7-fold by introducing a gap between the collimator and the skin. CONCLUSIONS Proton minibeams can be implemented at existing particle therapy centers. Because they spare the shallow tissues, they could augment the efficacy of proton therapy and light particle therapy, particularly in treating tumors that benefit from sparing of proximal tissues such as pediatric brain tumors. They should also allow hypofractionated treatment of all tumors by allowing the use of higher incident doses with less concern about proximal tissue damage.

A 4D-optimization concept for scanned ion beam therapy

BACKGROUND AND PURPOSE Scanned carbon beam therapy offers advantageous dose distributions and an increased biological effect. Treating moving targets is complex due to sensitivity to range changes and interplay. We propose a 4D treatment planning concept that considers motion during particle number optimization. MATERIAL AND METHODS The target was subdivided into sectors, one for each motion phase of a 4D-CT. Each sector was non-rigidly transformed to its motion phase and there targeted by a dedicated raster field (RST). Therefore, the resulting 4D-RST compensated target motion and range changes. A 4D treatment control system (TCS) was needed for synchronized delivery to the measured patient motion. 4D-optimized plans were simulated for 9 NSCLC lung cancer patients and compared to static irradiation at end-exhale. A prototype TCS was implemented and successfully tested in a film experiment. RESULTS The 4D-optimized treatment plan resulted in only slightly lower dose coverage of the target compared to static optimization, with V 95% of 97.9% (median, range 96.5-99.4%) vs. 99.3% (98.5-99.8%), with negligible overdose. The conformity number was comparable at 88.2% (85.1-92.5%) vs. 85.2% (79.9-91.2%) for 4D and static, respectively. CONCLUSION We implemented and tested a 4D treatment plan optimization method resulting in highly conformal dose delivery.

Validation of a track repeating algorithm for intensity modulated proton therapy: Clinical cases study

Monte Carlo (MC) methods are acknowledged as the most accurate technique to calculate dose distributions. However, due its lengthy calculation times, they are difficult to utilize in the clinic or for large retrospective studies. Track-repeating algorithms, based on MC-generated particle track data in water, accelerate dose calculations substantially, while essentially preserving the accuracy of MC. In this study, we present the validation of an efficient dose calculation algorithm for intensity modulated proton therapy, the fast dose calculator (FDC), based on a track-repeating technique. We validated the FDC algorithm for 23 patients, which included 7 brain, 6 head-and-neck, 5 lung, 1 spine, 1 pelvis and 3 prostate cases. For validation, we compared FDC-generated dose distributions with those from a full-fledged Monte Carlo based on GEANT4 (G4). We compared dose-volume-histograms, 3D-gamma-indices and analyzed a series of dosimetric indices. More than 99% of the voxels in the voxelized phantoms describing the patients have a gamma-index smaller than unity for the 2%/2 mm criteria. In addition the difference relative to the prescribed dose between the dosimetric indices calculated with FDC and G4 is less than 1%. FDC reduces the calculation times from 5 ms per proton to around 5 μs.

Implementation of an Analytical Model for Leakage Neutron Equivalent Dose in a Proton Radiotherapy Planning System

Equivalent dose from neutrons produced during proton radiotherapy increases the predicted risk of radiogenic late effects. However, out-of-field neutron dose is not taken into account by commercial proton radiotherapy treatment planning systems. The purpose of this study was to demonstrate the feasibility of implementing an analytical model to calculate leakage neutron equivalent dose in a treatment planning system. Passive scattering proton treatment plans were created for a water phantom and for a patient. For both the phantom and patient, the neutron equivalent doses were small but non-negligible and extended far beyond the therapeutic field. The time required for neutron equivalent dose calculation was 1.6 times longer than that required for proton dose calculation, with a total calculation time of less than 1 h on one processor for both treatment plans. Our results demonstrate that it is feasible to predict neutron equivalent dose distributions using an analytical dose algorithm for individual patients with irregular surfaces and internal tissue heterogeneities. Eventually, personalized estimates of neutron equivalent dose to organs far from the treatment field may guide clinicians to create treatment plans that reduce the risk of late effects.

Robustness of target dose coverage to motion uncertainties for scanned carbon ion beam tracking therapy of moving tumors.

Beam tracking with scanned carbon ion radiotherapy achieves highly conformal target dose by steering carbon pencil beams to follow moving tumors using real-time magnetic deflection and range modulation. The purpose of this study was to evaluate the robustness of target dose coverage from beam tracking in light of positional uncertainties of moving targets and beams. To accomplish this, we simulated beam tracking for moving targets in both water phantoms and a sample of lung cancer patients using a research treatment planning system. We modeled various deviations from perfect tracking that could arise due to uncertainty in organ motion and limited precision of a scanned ion beam tracking system. We also investigated the effects of interfractional changes in organ motion on target dose coverage by simulating a complete course of treatment using serial (weekly) 4DCTs from six lung cancer patients. For perfect tracking of moving targets, we found that target dose coverage was high ([Formula: see text] was 94.8% for phantoms and 94.3% for lung cancer patients, respectively) but sensitive to changes in the phase of respiration at the start of treatment and to the respiratory period. Phase delays in tracking the moving targets led to large degradation of target dose coverage (up to 22% drop for a 15° delay). Sensitivity to technical uncertainties in beam tracking delivery was minimal for a lung cancer case. However, interfractional changes in anatomy and organ motion led to large decreases in target dose coverage (target coverage dropped approximately 8% due to anatomy and motion changes after 1 week). Our findings provide a better understand of the importance of each of these uncertainties for beam tracking with scanned carbon ion therapy and can be used to inform the design of future scanned ion beam tracking systems.

Risk-optimized proton therapy to minimize radiogenic second cancers

Proton therapy confers substantially lower predicted risk of second cancer compared with photon therapy. However, no previous studies have used an algorithmic approach to optimize beam angle or fluence-modulation for proton therapy to minimize those risks. The objectives of this study were to demonstrate the feasibility of risk-optimized proton therapy and to determine the combination of beam angles and fluence weights that minimizes the risk of second cancer in the bladder and rectum for a prostate cancer patient. We used 6 risk models to predict excess relative risk of second cancer. Treatment planning utilized a combination of a commercial treatment planning system and an in-house risk-optimization algorithm. When normal-tissue dose constraints were incorporated in treatment planning, the risk model that incorporated the effects of fractionation, initiation, inactivation, repopulation and promotion selected a combination of anterior and lateral beams, which lowered the relative risk by 21% for the bladder and 30% for the rectum compared to the lateral-opposed beam arrangement. Other results were found for other risk models.

Efficient interplay effect mitigation for proton pencil beam scanning by spot-adapted layered repainting evenly spread out over the full breathing cycle

PURPOSE To develop and implement a practical repainting method for efficient interplay effect mitigation in proton pencil beam scanning (PBS). METHODS AND MATERIALS A new flexible repainting scheme with spot-adapted numbers of repainting evenly spread out over the whole breathing cycle (assumed to be 4 seconds) was developed. Twelve fields from 5 thoracic and upper abdominal PBS plans were delivered 3 times using the new repainting scheme to an ion chamber array on a motion stage. One time was static and 2 used 4-second, 3-cm peak-to-peak sinusoidal motion with delivery started at maximum inhalation and maximum exhalation. For comparison, all dose measurements were repeated with no repainting and with 8 repaintings. For each motion experiment, the 3%/3-mm gamma pass rate was calculated using the motion-convolved static dose as the reference. Simulations were first validated with the experiments and then used to extend the study to 0- to 5-cm motion magnitude, 2- to 6-second motion periods, patient-measured liver tumor motion, and 1- to 6-fraction treatments. The effect of the proposed method was evaluated for the 5 clinical cases using 4-dimensional (4D) dose reconstruction in the planning 4D computed tomography scan. The target homogeneity index, HI = (D2 - D98)/Dmean, of a single-fraction delivery is reported, where D2 and D98 is the dose delivered to 2% and 98% of the target, respectively, and Dmean is the mean dose. RESULTS The gamma pass rates were 59.6% ± 9.7% with no repainting, 76.5% ± 10.8% with 8 repaintings, and 92.4% ± 3.8% with the new repainting scheme. Simulations reproduced the experimental gamma pass rates with a 1.3% root-mean-square error and demonstrated largely improved gamma pass rates with the new repainting scheme for all investigated motion scenarios. One- and two-fraction deliveries with the new repainting scheme had gamma pass rates similar to those of 3-4 and 6-fraction deliveries with 8 repaintings. The mean HI for the 5 clinical cases was 14.2% with no repainting, 13.7% with 8 repaintings, 12.0% with the new repainting scheme, and 11.6% for the 4D dose without interplay effects. CONCLUSIONS A novel repainting strategy for efficient interplay effect mitigation was proposed, implemented, and shown to outperform conventional repainting in experiments, simulations, and dose reconstructions. This strategy could allow for safe and more optimal clinical delivery of thoracic and abdominal proton PBS and better facilitate hypofractionated and stereotactic treatments.

Comparative Risk Predictions of Second Cancers After Carbon-Ion Therapy Versus Proton Therapy

PURPOSE This work proposes a theoretical framework that enables comparative risk predictions for second cancer incidence after particle beam therapy for different ion species for individual patients, accounting for differences in relative biological effectiveness (RBE) for the competing processes of tumor initiation and cell inactivation. Our working hypothesis was that use of carbon-ion therapy instead of proton therapy would show a difference in the predicted risk of second cancer incidence in the breast for a sample of Hodgkin lymphoma (HL) patients. METHODS AND MATERIALS We generated biologic treatment plans and calculated relative predicted risks of second cancer in the breast by using two proposed methods: a full model derived from the linear quadratic model and a simpler linear-no-threshold model. RESULTS For our reference calculation, we found the predicted risk of breast cancer incidence for carbon-ion plans-to-proton plan ratio, , to be 0.75 ± 0.07 but not significantly smaller than 1 (P=.180). CONCLUSIONS Our findings suggest that second cancer risks are, on average, comparable between proton therapy and carbon-ion therapy.

Comparison of multi-institutional Varian ProBeam pencil beam scanning proton beam commissioning data

Purpose Commissioning beam data for proton spot scanning beams are compared for the first two Varian ProBeam sites in the United States, at the Maryland Proton Treatment Center (MPTC) and Scripps Proton Therapy Center (SPTC). In addition, the extent to which beams can be matched between gantry rooms at MPTC is investigated. Method Beam data for the two sites were acquired with independent dosimetry systems and compared. Integrated depth dose curves (IDDs) were acquired with Bragg peak ion chambers in a 3D water tank for pencil beams at both sites. Spot profiles were acquired at different distances from the isocenter at a gantry angle of 0° as well as a function of gantry angles. Absolute dose calibration was compared between SPTC and the gantries at MPTC. Dosimetric verification of test plans, output as a function of gantry angle, monitor unit (MU) linearity, end effects, dose rate dependence, and plan reproducibility were compared for different gantries at MPTC. Results The IDDs for the two sites were similar, except in the plateau region, where the SPTC data were on average 4.5% higher for lower energies. This increase in the plateau region decreased as energy increased, with no marked difference for energies higher than 180 MeV. Range in water coincided for all energies within 0.5 mm. The sigmas of the spot profiles in air were within 10% agreement at isocenter. This difference increased as detector distance from the isocenter increased. Absolute doses for the gantries measured at both sites were within 1% agreement. Test plans, output as function of gantry angle, MU linearity, end effects, dose rate dependence, and plan reproducibility were all within tolerances given by TG142. Conclusion Beam data for the two sites and between different gantry rooms were well matched.

Potential of discrete Gaussian edge feathering method for improving abutment dosimetry in eMLC-delivered segmented-field electron conformal therapy

PURPOSE The purpose of this work was to investigate the potential of discrete Gaussian edge feathering of the higher energy electron fields for improving abutment dosimetry in the planning volume when using an electron multileaf collimator (eMLC) to deliver segmented-field electron conformal therapy (ECT). METHODS A discrete (five-step) Gaussian edge spread function was used to match dose penumbras of differing beam energies (6-20 MeV) at a specified depth in a water phantom. Software was developed to define the leaf eMLC positions of an eMLC that most closely fit each electron field shape. The effect of 1D edge feathering of the higher energy field on dose homogeneity was computed and measured for segmented-field ECT treatment plans for three 2D PTVs in a water phantom, i.e., depth from the water surface to the distal PTV surface varied as a function of the x-axis (parallel to leaf motion) and remained constant along the y-axis (perpendicular to leaf motion). Additionally, the effect of 2D edge feathering was computed and measured for one radially symmetric, 3D PTV in a water phantom, i.e., depth from the water surface to the distal PTV surface varied as a function of both axes. For the 3D PTV, the feathering scheme was evaluated for 0.1-1.0-cm leaf widths. Dose calculations were performed using the pencil beam dose algorithm in the Pinnacle(3) treatment planning system. Dose verification measurements were made using a prototype eMLC (1-cm leaf width). RESULTS 1D discrete Gaussian edge feathering reduced the standard deviation of dose in the 2D PTVs by 34, 34, and 39%. In the 3D PTV, the broad leaf width (1 cm) of the eMLC hindered the 2D application of the feathering solution to the 3D PTV, and the standard deviation of dose increased by 10%. However, 2D discrete Gaussian edge feathering with simulated eMLC leaf widths of 0.1-0.5 cm reduced the standard deviation of dose in the 3D PTV by 33-28%, respectively. CONCLUSIONS A five-step discrete Gaussian edge spread function applied in 2D improves the abutment dosimetry but requires an eMLC leaf resolution better than 1 cm.