[en] Purpose: The spot characteristics for proton pencil beam scanning (PBS) were measured and analyzed over a 16 month period, which included one major site configuration update and six cyclotron interventions. The results provide a reference to establish the quality assurance (QA) frequency and tolerance for proton pencil beam scanning. Methods: A simple treatment plan was generated to produce an asymmetric 9-spot pattern distributed throughout a field of 16 × 18 cm for each of 18 proton energies (100.0–226.0 MeV). The delivered fluence distribution in air was measured using a phosphor screen based CCD camera at three planes perpendicular to the beam line axis (x-ray imaging isocenter and up/down stream 15.0 cm). The measured fluence distributions for each energy were analyzed using in-house programs which calculated the spot sizes and positional deviations of the Gaussian shaped spots. Results: Compared to the spot characteristic data installed into the treatment planning system, the 16-month averaged deviations of the measured spot sizes at the isocenter plane were 2.30% and 1.38% in the IEC gantry x and y directions, respectively. The maximum deviation was 12.87% while the minimum deviation was 0.003%, both at the upstream plane. After the collinearity of the proton and x-ray imaging system isocenters was optimized, the positional deviations of the spots were all within 1.5 mm for all three planes. During the site configuration update, spot positions were found to deviate by 6 mm until the tuning parameters file was properly restored. Conclusions: For this beam delivery system, it is recommended to perform a spot size and position check at least monthly and any time after a database update or cyclotron intervention occurs. A spot size deviation tolerance of <15% can be easily met with this delivery system. Deviations of spot positions were <2 mm at any plane up/down stream 15 cm from the isocenter

[en] Purpose: To introduce a method to model the 3D dose distribution of laterally asymmetric proton beamlets resulting from collimation. The model enables rapid beamlet calculation for spot scanning (SS) delivery using a novel penumbra-reducing dynamic collimation system (DCS) with two pairs of trimmers oriented perpendicular to each other. Methods: Trimmed beamlet dose distributions in water were simulated with MCNPX and the collimating effects noted in the simulations were validated by experimental measurement. The simulated beamlets were modeled analytically using integral depth dose curves along with an asymmetric Gaussian function to represent fluence in the beam’s eye view (BEV). The BEV parameters consisted of Gaussian standard deviations (sigmas) along each primary axis (σ_x1,σ_x2,σ_y1,σ_y2) together with the spatial location of the maximum dose (μ_x,μ_y). Percent depth dose variation with trimmer position was accounted for with a depth-dependent correction function. Beamlet growth with depth was accounted for by combining the in-air divergence with Hong’s fit of the Highland approximation along each axis in the BEV. Results: The beamlet model showed excellent agreement with the Monte Carlo simulation data used as a benchmark. The overall passing rate for a 3D gamma test with 3%/3 mm passing criteria was 96.1% between the analytical model and Monte Carlo data in an example treatment plan. Conclusions: The analytical model is capable of accurately representing individual asymmetric beamlets resulting from use of the DCS. This method enables integration of the DCS into a treatment planning system to perform dose computation in patient datasets. The method could be generalized for use with any SS collimation system in which blades, leaves, or trimmers are used to laterally sharpen beamlets

[en] Purpose: To compare the dose distribution between customized planning (CP) and adopting a single plan (SP) in multifractionated high-dose-rate brachytherapy and to establish predictors for the necessity of CP in a given patient. Methods and Materials: A total of 50 computed tomography-based plans for 10 patients were evaluated. Each patient had received 6 Gy for five fractions. The clinical target volume and organs at risk (i.e., rectum, bladder, sigmoid, and small bowel) were delineated on each computed tomography scan. For the SP approach, the same dwell position and time was used for all fractions. For the CP approach, the dwell position and time were reoptimized for each fraction. Applicator position variation was determined by measuring the distance between the posterior bladder wall and the tandem at the level of the vaginal fornices. Results: The organs at risk D_2cc (dose to 2 cc volume) was increased with the SP approach. The dose variation was statistically similar between the tandem and ring and tandem and ovoid groups. The bladder D_2cc dose was 81.95-105.42 Gy₂ for CP and 82.11-122.49 Gy₂ for SP. In 5 of the 10 patients, the bladder would have been significantly overdosed with the SP approach. The variation of the posterior bladder wall distance from that in the first fraction was correlated with the increase in the bladder D_2cc (SP/CP), with a correlation coefficient of -0.59. Conclusion: Our results support the use of CP instead of the SP approach to help avoid a significant overdose to the bladder. This is especially true for a decrease in the posterior wall distance of ≥0.5 cm compared with that in the first fraction.

[en] Purpose: Accelerated tumor repopulation has significant implications in low–dose rate (LDR) brachytherapy. Repopulation onset time remains undetermined for cervical cancer. The purpose of this study was to determine the onset time of accelerated repopulation in cervical cancer, using clinical data. Methods and Materials: The linear quadratic (LQ) model extended for tumor repopulation was used to analyze clinical data and magnetic resonance imaging-based three-dimensional tumor volumetric regression data from 80 cervical cancer patients who received external beam radiotherapy (EBRT) and LDR brachytherapy. The LDR dose was converted to EBRT dose in 1.8-Gy fractions by using the LQ formula, and the total dose ranged from 61.4 to 99.7 Gy. Patients were divided into 11 groups according to total dose and treatment time. The tumor control probability (TCP) was calculated for each group. The least χ² method was used to fit the TCP data with two free parameters: onset time (T_k) of accelerated repopulation and number of clonogens (K), while other LQ model parameters were adopted from the literature, due to the limited patient data. Results: Among the 11 patient groups, TCP varied from 33% to 100% as a function of radiation dose and overall treatment time. Higher dose and shorter treatment duration were associated with higher TCP. Using the LQ model, we achieved the best fit with onset time T_k of 19 days and K of 139, with uncertainty ranges of (11, 22) days for T_k and (48, 1822) for K, respectively. Conclusion: This is the first report of accelerated repopulation onset time in cervical cancer, derived directly from clinical data by using the LQ model. Our study verifies the fact that accelerated repopulation does exist in cervical cancer and has a relatively short onset time. Dose escalation may be required to compensate for the effects of tumor repopulation if the radiation therapy course is protracted.

[en] Purpose: To study various standardized uptake value (SUV)-based approaches to ascertain the best strategy for delineating metabolic tumor volumes (MTV). Methods and Materials: Twenty-two consecutive previously treated esophageal cancer patients with positron emission tomography (PET) imaging and computed tomography (CT)-based radiotherapy plans were studied. At the level of the tumor epicenter, MTVs were delineated at 11 different thresholds: SUV ≥2, ≥2.5, ≥3, ≥3.5 (SUV_n); ≥40%, ≥45%, and ≥50% of the maximum (SUV_n%); and mean liver SUV + 1, 2, 3, and 4 standard deviations (SUV_Lnσ). The volume ratio and conformality index were determined between MTVs, and the corresponding CT/endoscopic ultrasound-based gross tumor volume (GTV) at the epicenter. Means were analyzed by one-way analysis of variance for repeated measures and further compared using a paired t test for repeated measures. Results: The mean conformality indices ranged from 0.33 to 0.48, being significantly (p < 0.05) closest to 1 at SUV_2.5 (0.47 ± 0.03) and SUV_L4σ (0.48 ± 0.03). The mean volume ratios ranged from 0.39 to 2.82, being significantly closest to 1 at SUV_2.5 (1.18 ± 0.36) and SUV_L4σ (1.09 ± 0.15). The mean value of the SUVs calculated using the SUV_L4σ approach was 2.4. Conclusions: Regardless of the SUV thresholding method used (i.e., absolute or relative to liver mean), a threshold of approximately 2.5 yields the highest conformality index and best approximates the CT-based GTV at the epicenter. These findings may ultimately aid radiation oncologists in the delineation of the entire GTV in esophageal cancer patients.