Onboard cone-beam CT-based replan evaluation for head and neck proton therapy.

PURPOSE: Quality assurance computed tomography (QACT) is the current clinical practice in proton therapy to evaluate the needs for replan. QACT could falsely indicate replan because of setup issues that would be solved on the treatment machine. Deforming the treatment planning CT (TPCT) to the pretreatment CBCT may eliminate this issue. We investigated the performance of replan evaluation based on deformed TPCT (TPCTdir) for proton head and neck (H&N) therapy. METHODS AND MATERIALS: Twenty-eight H&N datasets along with pretreatment CBCT and QACT were used to validate the method. The changes in body volume were analyzed between the no-replan and replan groups. The dose on the TPCTdir, the deformed QACT (QACTdir), and the QACT were calculated by applying the clinical plans to these image sets. Dosimetric parameters' changes, including ΔD95, ΔDmean, and ΔD1 for the clinical target volumes (CTVs) were calculated. Receiver operating characteristic curves for replan evaluation based on ΔD95 on QACT and TPCTdir were calculated, using ΔD95 on QACTdir as the reference. A threshold for replan based on ΔD95 on TPCTdir is proposed. The specificities for the proposed method were calculated. RESULTS: The changes in the body contour were 95.8 ± 83.8 cc versus 305.0 ± 235.0 cc (p < 0.01) for the no-replan and replan groups, respectively. The ΔD95, ΔDmean, and ΔD1 are all comparable for all the evaluations. The differences between TPCTdir and QACTdir evaluations were 0.30% ± 0.86%, 0.00 ± 0.22 Gy, and -0.17 ± 0.61 Gy for CTV ΔD95, ΔDmean, and ΔD1, respectively. The corresponding differences between the QACT and QACTdir were 0.12% ± 1.1%, 0.02 ± 0.32 Gy, and -0.01 ± 0.71 Gy. CTV ΔD95 > 2.6% in TPCTdir was chosen as the threshold to trigger QACT/replan. The corresponding specificity was 94% and 98% for the clinical practice and the proposed method, respectively. CONCLUSIONS: The replan evaluation based on TPCTdir provides better specificity than that based on the QACT.

Characterizing devices for validation of dose, dose rate, and LET in ultra high dose rate proton irradiations.

BACKGROUND: Ultra high dose rate (UHDR) radiotherapy using ridge filter is a new treatment modality known as conformal FLASH that, when optimized for dose, dose rate (DR), and linear energy transfer (LET), has the potential to reduce damage to healthy tissue without sacrificing tumor killing efficacy via the FLASH effect. PURPOSE: Clinical implementation of conformal FLASH proton therapy has been limited by quality assurance (QA) challenges, which include direct measurement of UHDR and LET. Voxel DR distributions and LET spectra at planning target margins are paramount to the DR/LET-related sparing of organs at risk. We hereby present a methodology to achieve experimental validation of these parameters. METHODS: Dose, DR, and LET were measured for a conformal FLASH treatment plan involving a 250-MeV proton beam and a 3D-printed ridge filter designed to uniformly irradiate a spherical target. We measured dose and DR simultaneously using a 4D multi-layer strip ionization chamber (MLSIC) under UHDR conditions. Additionally, we developed an "under-sample and recover (USRe)" technique for a high-resolution pixelated semiconductor detector, Timepix3, to avoid event pile-up and to correct measured LET at high-proton-flux locations without undesirable beam modifications. Confirmation of these measurements was done using a MatriXX PT detector and by Monte Carlo (MC) simulations. RESULTS: MC conformal FLASH computed doses had gamma passing rates of >95% (3 mm/3% criteria) when compared to MatriXX PT and MLSIC data. At the lateral margin, DR showed average agreement values within 0.3% of simulation at 100 Gy/s and fluctuations ∼10% at 15 Gy/s. LET spectra in the proximal, lateral, and distal margins had Bhattacharyya distances of <1.3%. CONCLUSION: Our measurements with the MLSIC and Timepix3 detectors shown that the DR distributions for UHDR scenarios and LET spectra using USRe are in agreement with simulations. These results demonstrate that the methodology presented here can be used effectively for the experimental validation and QA of FLASH treatment plans.

Optimization of hexagonal-pattern minibeams for spatially fractionated radiotherapy using proton beam scanning.

PURPOSE: In this study, we investigated computationally and experimentally a hexagonal-pattern array of spatially fractionated proton minibeams produced by proton pencil beam scanning (PBS) technique. Spatial fractionation of dose delivery with millimeter or submillimeter beam size has proven to be a promising approach to significantly increase the normal tissue tolerance. Our goals are to obtain an optimized minibeam design and to show that it is feasible to implement the optimized minibeams at the existing proton clinics. METHODS: An optimized minibeam arrangement is one that would produce high peak-to-valley dose ratios (PVDRs) in normal tissues and a PVDR approaching unity at the Bragg peak. Using Monte Carlo (MC) code TOPAS we simulated proton pencil beams that mimic those available at the existing proton therapy facilities and obtained a hexagonal-pattern array of minibeams by collimating the proton pencil beams through the 1-3 mm diameter pinholes of a collimator. We optimized the minibeam design by considering different combinations of parameters including collimator material and thickness (t), center-to-center (c-t-c) distance, and beam size. The optimized minibeam design was then evaluated for normal tissue sparing against the uniform pencil beam scanning (PBS) by calculating the therapeutic advantage (TA) in terms of cell survival fraction. Verification measurements using radiochromic films were performed at the Emory proton therapy center (EPTC). RESULTS: Optimized hexagonal-pattern minibeams having PVDRs of >10 at phantom surface and of >3 at depths up to 6 cm were achieved with 2 mm diameter modulated proton minibeams (with proton energies between 120 and 140 MeV) corresponding to a spread-out-Bragg-peak (SOBP) over the depth of 10-14 cm. The results of the film measurements agree with the MC results within 10%. The TA of the 2 mm minibeams against the uniform PBS is >3 from phantom surface to the depth of 5 cm and then smoothly drops to ~1.5 as it approaches the proximal edge of the SOBP. For 2 mm minibeams and 6 mm c-t-c distance, we delivered 1.72 Gy at SOBP for 7.2 × 7.2 × 4 cm3 volume in 48 s. CONCLUSIONS: We conclude that it is feasible to implement the optimized hexagonal-pattern 2 mm proton minibeam radiotherapy at the existing proton clinics, because desirable PVDRs and TAs are achievable and the treatment time is reasonable.