Commissioning and Validation of CT Number to SPR Calibration in Carbon Ion Therapy Facility.

Purpose: We performed computed tomography (CT)-stopping power ratio (SPR) calibration in a carbon-ion therapy facility and evaluated SPR estimation accuracy. Materials and Methods: A polybinary tissue model method was used for the calibration of CT numbers and SPR. As a verification by dose calculation, we created a virtual phantom to which the CT-SPR calibration table was applied. Then, SPR was calculated from the change in the range of the treatment planning beam when changing to 19 different CT numbers, and the accuracy of the treatment planning system (TPS) calculation of SPR values from the CT-SPR calibration table was validated. As a verification by measurement, 5 materials (water, milk, olive oil, ethanol, 40% K2HPO4) were placed in a container, and the SPR was obtained by measurement from the change in the range of the beam that passed through the materials. Results: The results of the dose calculations of the TPS showed that the results agreed within 1% for the lower CT numbers up to 1000 HU, but there was a difference of 3.0% in the higher CT number volume. The difference between the SPR calculated by TPS and the SPR caused by the difference in the energy of the incident particles agreed within 0.51%. The accuracy of SPR estimation was measured, and the error was within 2% for all materials tested. Conclusion: These results indicate that the SPR estimation errors are within the range of errors that can be expected in particle therapy. From commissioning and verification results, the CT-SPR calibration table obtained during this commissioning process is clinically applicable.

Evaluation of the availability of single-position treatment with a rotating gantry and the validity of deformable image registration dose assessment for pancreatic cancer carbon-ion radiotherapy.

BACKGROUND: This study aimed to evaluate the clinical acceptability of rotational gantry-based single-position carbon-ion radiotherapy (CIRT) to reduce the gastrointestinal (GI) dose in pancreatic cancer. We also evaluated the usefulness of the deformable image registration (DIR)-based dosimetry method for CIRT. MATERIAL AND METHODS: Fifteen patients with pancreatic cancer were analyzed. The treatment plans were developed for four beam angles in the supine (SP plan) and prone (PR plan) positions. In the case of using multiple positions, the treatment plan was created with two angles for each of the supine and prone position (SP + PR plan). Dose evaluation for multiple positions was performed in two ways: by directly adding the values of the DVH parameters for each position treatment plan (DVH sum), and by calculating the DVH parameters from the accumulative dose distribution created using DIR (DIR sum). The D2cc and D6cc of the stomach and duodenum were recorded for each treatment plan and dosimetry method and compared. RESULTS: There were no significant differences among any of the treatment planning and dosimetry methods (p > 0.05). The DVH parameters for the stomach and duodenum were higher in the PR plan and SP plan, respectively, and DVH sum tended to be between the SP and PR plans. DVH sum and DIR sum, DVH sum tended to be higher for D2cc and DIR sum tended to be higher for D6cc. CONCLUSION: There were no significant differences in the GI dose, which suggests that treatment with a simple workflow performed in one position should be clinically acceptable. In CIRT, DIR-based dosimetry should be carefully considered because of the potential for increased uncertainty due to the steep dose distributions.

A robust treatment planning approach for chest motion in postmastectomy chest wall intensity modulated radiation therapy.

PURPOSE: Chest wall postmastectomy radiation therapy (PMRT) should consider the effects of chest wall respiratory motion. The purpose of this study is to evaluate the effectiveness of robustness planning intensity modulated radiation therapy (IMRT) for respiratory movement, considering respiratory motion as a setup error. MATERIAL AND METHODS: This study analyzed 20 patients who underwent PMRT (10 left and 10 right chest walls). The following three treatment plans were created for each case and compared. The treatment plans are a planning target volume (PTV) plan (PP) that covers the PTV within the body contour with the prescribed dose, a virtual bolus plan (VP) that sets a virtual bolus in contact with the body surface and prescribing the dose that includes the PTV outside the body contour, and a robust plan (RP) that considers respiratory movement as a setup uncertainty and performs robust optimization. The isocenter was shifted to reproduce the chest wall motion pattern and the doses were recalculated for comparison for each treatment plan. RESULT: No significant difference was found between the PP and the RP in terms of the tumor dose in the treatment plan. In contrast, VP had 3.5% higher PTV Dmax and 5.5% lower PTV V95% than RP (p < 0.001). The RP demonstrated significantly higher lung V20Gy and Dmean by 1.4% and 0.4 Gy, respectively, than the PP. The RP showed smaller changes in dose distribution affected by chest wall motion and significantly higher tumor dose coverage than the PP and VP. CONCLUSION: We revealed that the RP demonstrated comparable tumor doses to the PP in treatment planning and was robust for respiratory motion compared to both the PP and the VP. However, the organ at risk dose in the RP was slightly higher; therefore, its clinical use should be carefully considered.

Error on the stopping power ratio of ERKODENT's mouthpiece for head and neck carbon ion radiotherapy treatment.

The errors on the stopping power ratio (SPR) of mouthpiece samples from ERKODENT were evaluated. Erkoflex and Erkoloc-pro from ERKODENT and samples that combined Erkoflex and Erkoloc-pro were computed tomography (CT)-scanned using head and neck (HN) protocol at the East Japan Heavy Ion Center (EJHIC), and the values were averaged to obtain the CT number. The integral depth dose of the Bragg curve with and without these samples was measured for 292.1, 180.9, and 118.8 MeV/u of the carbon-ion pencil beam using an ionization chamber with concentric electrodes at the horizontal port of the EJHIC. The average value of the water equivalent length (WEL) of each sample was obtained from the difference between the range of the Bragg curve and the thickness of the sample. To calculate the difference between the theoretical and measured values, the theoretical CT number and SPR value of the sample were calculated using the stoichiometric calibration method. Compared with the Hounsfield unit (HU)-SPR calibration curve used at the EJHIC, the SPR error on each measured and theoretical value was calculated. The WEL value of the mouthpiece sample had an error of approximately 3.5% in the HU-SPR calibration curve. From this error, it was evaluated that for a mouthpiece with a thickness of 10 mm, a beam range error of approximately 0.4 mm can occur, and for a mouthpiece with a thickness of 30 mm, a beam range error of approximately 1 mm can occur. For a beam passing through the mouthpiece in HN treatment, it would be practical to consider a mouthpiece margin of 1 mm to avoid beam range errors if ions pass through the mouthpiece.

Development of a quantitative analysis method for assessing patient body surface deformation using an optical surface tracking system.

This study aimed to develop a new method to quantitatively analyze body shape changes in patients during radiotherapy without additional radiation exposure using an optical surface tracking system. This method's accuracy was evaluated using a cubic phantom with a known shift. Surface images of three-dimensionally printed phantoms, which simulated the head and neck shapes of real patients before and after treatment, were used to create a deformation surface area histogram. The near-maximum deformation value covering an area of 2 cm2 in the surface image (Def-2cm2) was calculated. A volumetric modulated arc therapy (VMAT) plan was also created on the pre-treatment phantom, and the dose distribution was recalculated on the post-treatment phantom to compare the dose indices. Surface images of four patients were analyzed to evaluate Def-2cm2 and examine whether this method can be used in clinical cases. Experiments with the cubic phantom resulted in a mean deformation error of 0.08 mm. With head and neck phantoms, the Def-2cm2 value was 17.5 mm, and the dose that covered 95% of the planning target volume in the VMAT plan decreased by 11.7%, indicating that deformation of the body surface may affect the dose distribution. Although analysis of the clinical data showed no clinically relevant deformation in any of the cases, slight skin sagging and respiratory changes in the body surface were observed. The proposed method can quantitatively and accurately evaluate the deformation of a body surface. This method is expected to be used to make decisions regarding modifications to treatment plans.

Comparison of predictive performance for toxicity by accumulative dose of DVH parameter addition and DIR addition for cervical cancer patients.

We compared predictive performance between dose volume histogram (DVH) parameter addition and deformable image registration (DIR) addition for gastrointestinal (GI) toxicity in cervical cancer patients. A total of 59 patients receiving brachytherapy and external beam radiotherapy were analyzed retrospectively. The accumulative dose was calculated by three methods: conventional DVH parameter addition, full DIR addition and partial DIR addition. ${D}_{2{cm}^3}$, ${D}_{1{cm}^3}$ and ${D}_{0.1{cm}^3}$ (minimum doses to the most exposed 2 cm3, 1cm3 and 0.1 cm3 of tissue, respectively) of the rectum and sigmoid were calculated by each method. V50, V60 and V70 Gy (volume irradiated over 50, 60 and 70 Gy, respectively) were calculated in full DIR addition. The DVH parameters were compared between toxicity (≥grade1) and non-toxicity groups. The area under the curve (AUC) of the receiver operating characteristic (ROC) curves were compared to evaluate the predictive performance of each method. The differences between toxicity and non-toxicity groups in ${D}_{2{cm}^3}$ were 0.2, 5.7 and 3.1 Gy for the DVH parameter addition, full DIR addition and partial DIR addition, respectively. The AUCs of ${D}_{2{cm}^3}$ were 0.51, 0.67 and 0.57 for DVH parameter addition, full DIR addition and partial DIR addition, respectively. In full DIR addition, the difference in dose between toxicity and non-toxicity was the largest and AUC was the highest. AUCs of V50, V60 and V70 Gy were 0.51, 0.63 and 0.62, respectively, and V60 and V70 were high values close to the value of ${D}_{2{cm}^3}$ of the full DIR addition. Our results suggested that the full DIR addition may have the potential to predict toxicity more accurately than the conventional DVH parameter addition, and that it could be more effective to accumulate to all pelvic irradiation by DIR.

The impact of 4DCT-ventilation imaging-guided proton therapy on stereotactic body radiotherapy for lung cancer.

Functional lung avoidance during radiotherapy can help reduce pulmonary toxicity. This study assessed the potential impact of four-dimensional computed tomography (4DCT)-ventilation imaging-guided proton radiotherapy (PT) on stereotactic body radiotherapy (SBRT) by comparing it with three-dimensional conformal radiotherapy (3D-CRT) and volumetric modulated arc therapy (VMAT), which employ photon beams. Thirteen lung cancer patients who received SBRT with 3D-CRT were included in the study. 4DCT ventilation was calculated using the patients' 4DCT data, deformable image registration, and a density-change-based algorithm. Three functional treatment plans sparing the functional lung regions were developed for each patient using 3D-CRT, VMAT, and PT. The prescribed doses and dose constraints were based on the Radiation Therapy Oncology Group 0618 protocol. We evaluated the region of interest (ROI) and functional map-based dose-function metrics for 4DCT ventilation and the irradiated dose. Using 3D-CRT, VMAT, and PT, the percentages of the functional lung regions that received ≥ 5 Gy (fV5) were 26.0%, 21.9%, and 10.7%, respectively; the fV10 were 14.4%, 11.4%, and 9.0%, respectively; and fV20 were 6.5%, 6.4%, and 6.6%, respectively, and the functional mean lung doses (fMLD) were 5.6 Gy, 5.2 Gy, and 3.8 Gy, respectively. These results indicated that PT resulted in a significant reduction in fMLD, fV5, and fV10, but not fV20. The use of PT reduced the radiation to highly functional lung regions compared with those for 3D-CRT and VMAT while meeting all dose constraints.

Quantitative analysis of intra-fractional variation in CT-based image guided brachytherapy for cervical cancer patients.

We quantified intra-fractional dose variation and organ movement during CT-based 3D-image guided brachytherapy (3D-IGBT) in cervical cancer patients. Fifteen patients who underwent CT-based 3D-IGBT were studied. For all patients, pre-delivery CT for treatment planning after applicator insertion and post-delivery CT after dose delivery without changing the applicator position were acquired. Pre- and post-delivery CT were rigidly fused by matching the inserted applicator and planned dose on pre-delivery CT (pre-delivery dose) was mapped on post-delivery CT (post-delivery dose). D2, D1, and D0.1 cm3 of the rectum and bladder were compared between pre- and post-delivery doses with contours on each CT image. Organ movement and deformation was evaluated using deformation vector fields calculated by deformable image registration between pre- and post-delivery CT. We also evaluated dose variation and DVF between with and without a catheter to control filling. Differences in all DVH parameters were <±3% in physical dose and ± 5% in EQD2. However, a > 15% dose difference was found in 13.8% of the fractions in rectum D2 cm3 and in 11.1% of those in bladder D2 cm3. The mean value of DVF for bladder was larger than that of rectum, especially for the superior-inferior (S-I) direction. Insertion catheters in bladder reduced mean dose and DVF variation compared with that of without catheters. In fraction groups with large dose increasing, DVF in the S-I direction was significantly larger than that of other fraction groups. Our results indicated that preparation is needed to reduce changes in the S-I direction affect dose variation.

Measurement of incident position of hypervelocity particles on piezoelectric lead zirconate titanate detector.

A cosmic dust detector for use onboard a satellite is currently being developed by using piezoelectric lead zirconate titanate (PZT). The characteristics of the PZT detector have been studied by bombarding it with hypervelocity iron (Fe) particles supplied by a Van de Graaff accelerator. One central electrode and four peripheral electrodes were placed on the front surface of the PZT detector to measure the impact positions of the incident Fe particles. It was demonstrated that the point of impact on the PZT detector could be identified by using information on the time at which the first peak of the output signal obtained from each electrode appeared.