Minglei Kang

3D printer generated thorax phantom with mobile tumor for radiation dosimetry.

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 tumors 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.

Beam-specific planning target volumes incorporating 4D CT for pencil beam scanning proton therapy of thoracic tumors.

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 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.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 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

Use of a novel two-dimensional ionization chamber array for pencil beam scanning proton therapy beam quality assurance

The need to accurately and efficiently verify both output and dose profiles creates significant challenges in quality assurance of pencil beam scanning (PBS) proton delivery. A system for PBS QA has been developed that combines a new two-dimensional ionization chamber array in a waterproof housing that is scanned in a water phantom. The MatriXX PT has the same detector array arrangement as the standard MatriXXEvolution but utilizes a smaller 2 mm plate spacing instead of 5 mm. Because the bias voltage of the MatriXX PT and Evolution cannot be changed, PPC40 and FC65-G ionization chambers were used to assess recombination effects. The PPC40 is a parallel plate chamber with an electrode spacing of 2 mm, while the FC65-G is a Farmer chamber FC65-G with an electrode spacing of 2.8 mm. Three bias voltages (500, 200, and 100 V) were used for both detectors to determine which radiation type (continuous, pulse or pulse-scanned beam) could closely estimate Pion from the ratios of charges collected. In comparison with the MatriXXEvolution, a significant improvement in measurement of absolute dose with the MatriXX PT was observed. While dose uncertainty of the MatriXXEvolution can be up to 4%, it is <1% for the MatriXX PT. Therefore the MatriXXEvolution should not be used for QA of PBS for conditions in which ion recombination is not negligible. Farmer chambers should be used with caution for measuring the absolute dose of PBS beams, as the uncertainty of Pion can be <1%; chambers with an electrode spacing of 2 mm or smaller are recommended. PACS number: 87.53.Qc.The need to accurately and efficiently verify both output and dose profiles creates significant challenges in quality assurance of pencil beam scanning (PBS) proton delivery. A system for PBS QA has been developed that combines a new two‐dimensional ionization chamber array in a waterproof housing that is scanned in a water phantom. The MatriXX PT has the same detector array arrangement as the standard MatriXXEvolution but utilizes a smaller 2 mm plate spacing instead of 5 mm. Because the bias voltage of the MatriXX PT and Evolution cannot be changed, PPC40 and FC65‐G ionization chambers were used to assess recombination effects. The PPC40 is a parallel plate chamber with an electrode spacing of 2 mm, while the FC65‐G is a Farmer chamber FC65‐G with an electrode spacing of 2.8 mm. Three bias voltages (500, 200, and 100 V) were used for both detectors to determine which radiation type (continuous, pulse or pulse‐scanned beam) could closely estimate Pion from the ratios of charges collected. In comparison with the MatriXXEvolution, a significant improvement in measurement of absolute dose with the MatriXX PT was observed. While dose uncertainty of the MatriXXEvolution can be up to 4%, it is <1% for the MatriXX PT. Therefore the MatriXXEvolution should not be used for QA of PBS for conditions in which ion recombination is not negligible. Farmer chambers should be used with caution for measuring the absolute dose of PBS beams, as the uncertainty of Pion can be <1%; chambers with an electrode spacing of 2 mm or smaller are recommended. PACS number: 87.53.Qc

Evaluation of Motion Mitigation using Abdominal Compression in the Clinical Implementation of Pencil Beam Scanning Proton Therapy of Liver Tumors

Purpose: To determine whether individual liver tumor patients can be safely treated with pencil beam scanning proton therapy. This study reports a planning preparation workflow that can be used for beam angle selection and the evaluation of the efficacy of abdominal compression (AC) to mitigate motion. Methods: Four‐dimensional computed tomography scans (4DCT) with and without AC were available from 10 liver tumor patients with fluoroscopy‐proven motion reduction by AC, previously treated using photons. For each scan, the motion amplitudes and the motion‐induced variation of water‐equivalent thickness (ΔWET) in each voxel of the target volume were evaluated during treatment plan preparation. Optimal proton beam angles were selected after volume analysis of the respective beam‐specific planning target volume (BSPTV). M⊥80 and ΔWET80 derived from the 80th percentiles of perpendicular motion amplitude (M⊥) and ΔWET were compared with and without AC. Proton plans were created on the average CT to achieve target coverage similar to that of the conventional photon treatments. 4D dynamic dose calculation was performed postplan by synchronizing proton beam delivery timing patterns to the 4DCT phases to assess interplay and fractionation effects, and to determine motion criteria for subsequent patient treatment. Results: Selected coplanar beam angles ranged between 180° and 39°, primarily from right lateral (oblique) and posterior (oblique) directions. While AC produced a significant reduction in mean Liver‐GTV dose, any reduction in mean heart dose was patient dependent and not significant. Similarly, AC resulted in reductions in M⊥, ΔWET, and BSPTV volumes and improved dose degradation (ΔD95 and ΔD1) within the CTV. For small motion (M⊥80 < 7 mm and ΔWET80 < 5 mm), motion mitigation was not needed. For moderate motion (M⊥80 7–10 mm or ΔWET80 5–7 mm), AC produced a modest improvement. For large motion (M⊥80 > 10 mm or ΔWET80 > 7 mm), AC and/or some other form of mitigation strategies were required. Conclusion: A workflow for screening patients’ motion characteristics and optimizing beam angle selection was established for the pencil beam scanning proton therapy treatment of liver tumors. Abdominal compression was found to be useful at mitigation of moderate and large motion.

A study of the beam-specific interplay effect in proton pencil beam scanning delivery in lung cancer

Abstract 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.

A benchmarking method to evaluate the accuracy of a commercial proton monte carlo pencil beam scanning treatment planning system

&NA; AcurosPT is a Monte Carlo algorithm in the Eclipse 13.7 treatment planning system, which is designed to provide rapid and accurate dose calculations for proton therapy. Computational run‐time in minimized by simplifying or eliminating less significant physics processes. In this article, the accuracy of AcurosPT was benchmarked against both measurement and an independent MC calculation, TOPAS. Such a method can be applied to any new MC calculation for the detection of potential inaccuracies. To validate multiple Coulomb scattering (MCS) which affects primary beam broadening, single spot profiles in a Solidwater® phantom were compared for beams of five selected proton energies between AcurosPT, measurement and TOPAS. The spot Gaussian sigma in AcurosPT was found to increase faster with depth than both measurement and TOPAS, suggesting that the MCS algorithm in AcurosPT overestimates the scattering effect. To validate AcurosPT modeling of the halo component beyond primary beam broadening, field size factors (FSF) were compared for multi‐spot profiles measured in a water phantom. The FSF for small field sizes were found to disagree with measurement, with the disagreement increasing with depth. Conversely, TOPAS simulations of the same FSF consistently agreed with measurement to within 1.5%. The disagreement in absolute dose between AcurosPT and measurement was smaller than 2% at the mid‐range depth of multi‐energy beams. While AcurosPT calculates acceptable dose distributions for typical clinical beams, users are cautioned of potentially larger errors at distal depths due to overestimated MCS and halo implementation.

Technical Note: Validation of halo modeling for proton pencil beam spot scanning using a quality assurance test pattern.

PURPOSE The purpose of this paper is to demonstrate the utility of a comprehensive test pattern in validating calculation models that include the halo component (low-dose tails) of proton pencil beam scanning (PBS) spots. Such a pattern has been used previously for quality assurance purposes to assess spot shape, position, and dose. METHODS In this study, a scintillation detector was used to measure the test pattern in air at isocenter for two proton beam energies (115 and 225 MeV) of two IBA universal nozzles (UN #1 and UN #2). Planar measurements were compared with calculated dose distributions based on the weighted superposition of location-independent (UN #1) or location-dependent (UN #2) spot profiles, previously measured using a pair-magnification method and between two nozzles. RESULTS Including the halo component below 1% of the central dose is shown to improve the gamma-map comparison between calculation and measurement from 94.9% to 98.4% using 2 mm/2% criteria for the 115 MeV proton beam of UN #1. In contrast, including the halo component below 1% of the central dose does not improve the gamma agreement for the 115 MeV proton beam of UN #2, due to the cutoff of the halo component at off-axis locations. When location-dependent spot profiles are used for calculation instead of spot profiles at central axis, the gamma agreement is improved from 98.0% to 99.5% using 2 mm/2% criteria. The two nozzles clearly have different characteristics, as a direct comparison of measured data shows a passing rate of 89.7% for the 115 MeV proton beam. At 225 MeV, the corresponding gamma comparisons agree better between measurement and calculation, and between measurements in the two nozzles. CONCLUSIONS In addition to confirming the primary component of individual PBS spot profiles, a comprehensive test pattern is useful for the validation of the halo component at off-axis locations, especially for low energy protons.

Thoracic radiotherapy (TRT) improved survival in both oligo- and polymetastatic extensive stage small cell lung cancer

There has been no previous study on the efficacy of the thoracic radiotherapy (TRT) in oligometastatic or polymetastatic extensive stage small-cell lung cancer (ES-SCLC) to the overall survival (OS). In a group of 270 ES-SCLC cases retrospective study, 78 patients (28.9%) had oligometastases and 192 (71.1%) had polymetastases, among which 51 oligometastatic patients (65.4%) and 93 polymetastatic patients (51.6%) received TRT. Propensity score matching (PSM) was utilized. The 2-year OS, progression free survival (PFS) and local control (LC) in oligometastatic and polymetastatic patients were 22.8% and 4.5% (p < 0.001), 12.0% and 3.8% (p < 0.001), and 36.7% and 6.1% (p < 0.001), respectively. The 2-year OS in oligometastatic patients with the chemotherapy + radiotherapy and chemotherapy alone were 25.2% and 12.7% (p = 0.002), in contrast to 10.0% and 6.8% (p = 0.030) in polymetastatic patients. The estimated hazard ratios for survival were 2.9 and 1.7 for both oligometastatic and polymetastatic patients with radiotherapy. The polymetastatic group has a lower LC (6.1% v.s. 36.7%, (p < 0.001)), due to polymetastases patients receiving involved-sites radiotherapy with low dose schemas. TRT improved OS of patients with oligometastases and polymetastases. Our study demonstrated that aggressive TRT might be a suitable addition of chemotherapy when treating ES-SCLC patients with oligometastases and polymetastases.

Validation and clinical implementation of an accurate Monte Carlo code for pencil beam scanning proton therapy

Abstract Monte Carlo (MC)‐based dose calculations are generally superior to analytical dose calculations (ADC) in modeling the dose distribution for proton pencil beam scanning (PBS) treatments. The purpose of this paper is to present a methodology for commissioning and validating an accurate MC code for PBS utilizing a parameterized source model, including an implementation of a range shifter, that can independently check the ADC in commercial treatment planning system (TPS) and fast Monte Carlo dose calculation in opensource platform (MCsquare). The source model parameters (including beam size, angular divergence and energy spread) and protons per MU were extracted and tuned at the nozzle exit by comparing Tool for Particle Simulation (TOPAS) simulations with a series of commissioning measurements using scintillation screen/CCD camera detector and ionization chambers. The range shifter was simulated as an independent object with geometric and material information. The MC calculation platform was validated through comprehensive measurements of single spots, field size factors (FSF) and three‐dimensional dose distributions of spread‐out Bragg peaks (SOBPs), both without and with the range shifter. Differences in field size factors and absolute output at various depths of SOBPs between measurement and simulation were within 2.2%, with and without a range shifter, indicating an accurate source model. TOPAS was also validated against anthropomorphic lung phantom measurements. Comparison of dose distributions and DVHs for representative liver and lung cases between independent MC and analytical dose calculations from a commercial TPS further highlights the limitations of the ADC in situations of highly heterogeneous geometries. The fast MC platform has been implemented within our clinical practice to provide additional independent dose validation/QA of the commercial ADC for patient plans. Using the independent MC, we can more efficiently commission ADC by reducing the amount of measured data required for low dose “halo” modeling, especially when a range shifter is employed.

Receipt of thoracic radiation therapy and radiotherapy dose are correlated with outcomes in a retrospective study of three hundred and six patients with extensive stage small-cell lung cancer

BACKGROUND The importance of the thoracic radiation therapy (TRT) dose has not been clearly defined in extensive stage small-cell lung cancer (ES-SCLC) and it is unclear whether improved TRT dose translates into a survival benefit. METHODS 306 patients with ES-SCLC were retrospectively reviewed, of which 170 received IMRT/CRT fractionation RT after ChT, and 136 received chemotherapy (ChT) alone. We adopted the time-adjusted BED (tBED) for effective dose fractionation calculation. Due to the nonrandomized nature of this study, we compared the ChT+RT with ChT groups that matched on possible confounding variables. RESULTS Patients achieved 2-year OS, PFS and LC rates of 19.7%, 10.7% and 28.4%, respectively. After propensity score matching, (113 cases for each group), the rates of OS, PFS and LC at 2 years were 21.4%, 7.7% and 34.5% for ChT+TRT, and 10.3% (p<0.001), 4.6% (p<0.001) and 6.3% for ChT only (p<0.001), respectively. Among propensity score matching patients, 56 cases for each group received the high dose (tBED>50 Gy) TRT and received low dose (tBED≤50 Gy) TRT. Two-year OS, PFS and LC rates were 32.3%, 15.3% and 47.1% for the high dose compared with 17.0% (p<0.001), 12.9% (p=0.097) and 34.7% (p=0.029) for low dose radiotherapy. CONCLUSIONS TRT added to ChT improved ES-SCLC patient OS. High dose TRT improved OS over lower doses. Our results suggest that high-dose thoracic radiation therapy may be a reasonable consideration in select patients with ES-SCLC.