To quantitatively assess target and organs-at-risk (OAR) dose rate based on three proposed proton PBS dose rate metrics and study FLASH intensity-modulated proton therapy (IMPT) treatment planning using transmission beams. An in-house FLASH planning platform was developed to optimize transmission (shoot-through) plans for nine consecutive lung cancer patients previously planned with proton SBRT. Dose and dose rate calculation codes were developed to quantify three types of dose rate calculation methods (dose-averaged dose rate (DADR), average dose rate (ADR), and dose-threshold dose rate (DTDR)) based on both phantom and patient treatment plans. Two different minimum MU/spot settings were used to optimize two different dose regimes, 34-Gy in one fraction and 45-Gy in three fractions. The OAR sparing and target coverage can be optimized with good uniformity (hotspot < 110% of prescription dose). ADR, accounting for the spot dwelling and scanning time, gives the lowest dose rate; DTDR, not considering this time but a dose-threshold, gives an intermediate dose rate, whereas DADR gives the highest dose rate without considering any time or dose-threshold. All three dose rates attenuate along the beam direction, and the highest dose rate regions often occur on the field edge for ADR and DTDR, whereas DADR has a better dose rate uniformity. The differences in dose rate metrics have led a large variation for OARs dose rate assessment, posing challenges to FLASH clinical implementation. This is the first attempt to study the impact of the dose rate models, and more investigations and evidence for the details of proton PBS FLASH parameters are needed to explore the correlation between FLASH efficacy and the dose rate metrics.

Purpose: While transmission proton beams have been demonstrated to achieve ultra-high dose rate FLASH therapy delivery, they are unable to spare normal tissues distal to the target. This study aims to compare FLASH treatment planning using single energy Bragg peak proton beams versus transmission proton beams in lung tumors and to evaluate Bragg peak plan optimization, characterize plan quality, and quantify organ-at-risk (OAR) sparing. Materials and Methods: Both Bragg peak and transmission plans were optimized using an in-house platform for 10 consecutive lung patients previously treated with proton stereotactic body radiation therapy (SBRT). To bring the dose rate up to the FLASH-RT threshold, Bragg peak plans with a minimum MU/spot of 1200 and transmission plans with a minimum MU/spot of 400 were developed. Two common prescriptions, 34 Gy in 1 fraction and 54 Gy in 3 fractions, were studied with the same beam arrangement for both Bragg peak and transmission plans (n = 40 plans). RTOG 0915 dosimetry metrics and dose rate metrics based on different dose rate calculations, including average dose rate (ADR), dose-averaged dose rate (DADR), and dose threshold dose rate (DTDR), were investigated. We then evaluated the effect of beam angular optimization on the Bragg peak plans to explore the potential for superior OAR sparing. Results: Bragg peak plans significantly reduced doses to several OAR dose parameters, including lung V_7.4Gy and V_7Gy by 32.0% (p < 0.01) and 30.4% (p < 0.01) for 34Gy/fx plans, respectively; and by 40.8% (p < 0.01) and 41.2% (p < 0.01) for 18Gy/fx plans, respectively, compared with transmission plans. Bragg peak plans have ~3% less in DADR and ~10% differences in mean OARs in DTDR and DADR relative to transmission plans due to the larger portion of lower dose regions of Bragg peak plans. With angular optimization, optimized Bragg peak plans can further reduce the lung V_7Gy by 20.7% (p < 0.01) and V_7.4Gy by 19.7% (p < 0.01) compared with Bragg peak plans without angular optimization while achieving a similar 3D dose rate distribution. Conclusion: The single-energy Bragg peak plans achieve superior dosimetry performances in OARs to transmission plans with comparable dose rate performances for lung cancer FLASH therapy. Beam angle optimization can further improve the OAR dosimetry parameters with similar 3D FLASH dose rate coverage.

The purpose of this study is to determine whether organ sparing and target coverage can be simultaneously maintained for pencil beam scanning (PBS) proton therapy treatment of thoracic tumors in the presence of motion, stopping power uncertainties, and patient setup variations. Ten consecutive patients that were previously treated with proton therapy to 66.6/1.8 Gy (RBE) using double scattering (DS) were replanned with PBS. Minimum and maximum intensity images from 4D CT were used to introduce flexible smearing in the determination of the beam specific PTV (BSPTV). Datasets from eight 4D CT phases, using ?3\% uncertainty in stopping power and ?3?\text{mm} uncertainty in patient setup in each direction, were used to create 8?12?10=960 PBS plans for the evaluation of 10 patients. Plans were normalized to provide identical coverage between DS and PBS. The average lung V20, V5, and mean doses were reduced from 29.0%, 35.0%, and 16.4 Gy with DS to 24.6%, 30.6%, and 14.1 Gy with PBS, respectively. The average heart V30 and V45 were reduced from 10.4% and 7.5% in DS to 8.1% and 5.4% for PBS, respectively. Furthermore, the maximum spinal cord, esophagus, and heart doses were decreased from 37.1 Gy, 71.7 Gy, and 69.2 Gy with DS to 31.3 Gy, 67.9 Gy, and 64.6 Gy with PBS. The conformity index (CI), homogeneity index (HI), and global maximal dose were improved from 3.2, 0.08, 77.4 Gy with DS to 2.8, 0.04, and 72.1 Gy with PBS. All differences are statistically significant, with p-values <0.05 , with the exception of the heart V45 ( p=0.146 ). PBS with BSPTV achieves better organ sparing and improves target coverage using a repainting method for the treatment of thoracic tumors. Incorporating motion-related uncertainties is essential. PACS number: 87.55.D

This article describes the design, construction, and properties of an anthropomorphic thorax phantom with a moving surrogate tumor. This novel phantom permits detection of dose both inside and outside a moving tumor and within the substitute lung tissue material. A 3D printer generated the thorax shell composed of a chest wall, spinal column, and posterior regions of the phantom. Images of a computed tomography scan of the thorax from a patient with lung cancer provided the template for the 3D printing. The plastic phantom is segmented into two materials representing the muscle and bones, and its geometry closely matches a patient. A surrogate spherical plastic tumor controlled by a 3D linear stage simulates a lung tumor’s trajectory during normal breathing. Sawdust emulates the lung tissue in terms of average and distribution in Hounsfield numbers. The sawdust also provides a forgiving medium that permits tumor motion and sandwiching of radiochromic film inside the mobile surrogate plastic tumor for dosimetry. A custom cork casing shields the film and tumor and eliminates film bending during extended scans. The phantom, lung tissue surrogate, and radiochromic film are exposed to a seven field plan based on an ECLIPSE plan for 6 MV photons from a Trilogy machine delivering 230 cGy to the isocenter. The dose collected in a sagittal plane is compared to the calculated plan. Gamma analysis finds 8.8% and 5.5% gamma failure rates for measurements of large amplitude trajectory and static measurements relative to the large amplitude plan, respectively. These particular gamma analysis results were achieved using parameters of 3% dose and 3 mm, for regions receiving doses >150 cGy. The plan assumes a stationary detection grid unlike the moving radiochromic film and tissues. This difference was experimentally observed and motivated calculated dose distributions that incorporated the phase of the tumor periodic motion. These calculations modestly improve agreement between the measured and intended doses.

Background: For lung tumors with large motion amplitudes, the use of proton pencil beam scanning (PBS) can produce large dose errors. In this study, we assess under what circumstances PBS can be used to treat lung cancer patients who exhibit large tumor motion, based on the quantification of tumor motion and the dose interplay. Material and methods: PBS plans were optimized on average 4DCT datasets using a beam-specific PTV method for 10 consecutive patients with locally advanced non-small-cell-lung-cancer (NSCLC) treated with proton therapy to 6660/180 cGy. End inhalation (CT0) and end exhalation (CT50) were selected as the two extreme scenarios to acquire the relative stopping power ratio difference (Δrsp) for a respiration cycle. The water equivalent difference (ΔWET) per radiological path was calculated from the surface of patient to the iCTV by integrating the Δrsp of each voxel. The magnitude of motion of voxels within the target follows a quasi-Gaussian distribution. A motion index (MI (>5mm WET)), defined as the percentage of target voxels with an absolute integral ΔWET larger than 5 mm, was adopted as a metric to characterize interplay. To simulate the treatment process, 4D dose was calculated by accumulating the spot dose on the corresponding respiration phase to the reference phase CT50 by deformable image registration based on spot timing and patient breathing phase. Results: The study indicated that the magnitude of target underdose in a single fraction plan is proportional to the MI (p < .001), with larger motion equating to greater dose degradation and standard deviations. The target homogeneity, minimum, maximum and mean dose in the 4D dose accumulations of 37 fractions varied as a function of MI. Conclusions: This study demonstrated that MI can predict the level of dose degradation, which potentially serves as a clinical decision tool to assess whether lung cancer patients are potentially suitable to receive PBS treatment.