Skip to main content

Adaptive planning based on single beam optimization in passive scattering carbon ion radiotherapy for patients with pancreatic cancer



Daily anatomical deviations may distort the dose distribution in carbon ion radiotherapy (CIRT), which may cause treatment failure. Therefore, this study aimed to perform re-planning to maintain the dose coverage in patients with pancreatic cancer with passive scattering CIRT.


Eight patients with pancreatic cancer and 95 daily computed tomography (CT) sets were examined. Two types of adaptive plans based on new range compensators (RCs) (AP-1) and initial RCs (AP-2) were generated. In AP-2, each beam was optimized by manually adjusting the range shifter thickness and spread-out Bragg peak size to make dose reduction by < 3% of the original plan. Doses of the original plan with bone matching (BM) and tumor matching (TM) were examined for comparison. We calculated the accumulated dose using the contour and intensity-based deformable image registration algorithm. The dosimetric differences in respect to the original plan were compared between methods.


Using TM and BM, mean ± standard deviations of daily CTV V95 (%) difference from the original plan was − 5.1 ± 6.2 and − 8.8 ± 8.8, respectively, but 1.2 ± 3.4 in AP-1 and − 0.5 ± 2.1 in AP-2 (P < 0.001). AP-1 and AP-2 enabled to maintain a satisfactory accumulated dose in all patients. The dose difference was 1.2 ± 2.8, − 2,1 ± 1.7, − 7.1 ± 5.2, and − 16.5 ± 15.0 for AP-1, AP-2, TM, and BM, respectively. However, AP-2 caused a dose increase in the duodenum, especially in the left–right beam.


The possible dose deterioration should be considered when performing the BM, even TM. Re-planning based on single beam optimization in passive scattering CIRT seems an effective and safe method of ensuring the treatment robustness in pancreatic cancer. Further study is necessary to spare healthy tissues, especially the duodenum.


Carbon ion radiotherapy (CIRT) has emerged as a viable treatment option for unresectable pancreatic cancer in the last decades [1,2,3,4]. A multi-institutional study recently reported promising clinical outcomes in the overall survival and local control, which are superior to those obtained with conventional chemoradiation therapy or chemotherapy alone [1]. These are mainly attributable to its ability to provide a concentrated dose distribution to the target and better biological effect induced by a higher linear energy transfer [5, 6]. However, the carbon ion beam is highly sensitive to anatomical changes during treatment, which may cause a prominent dose distortion especially in pancreatic cancer.

Several studies have found that large position variation of pancreatic tumor and random gastrointestinal (GI) tract deformation exit between treatment fractions [7,8,9,10], unlike lung and liver tumors where good tumor reproducibility via tumor matching (TM) could account for interfractional deviations ensuring a robust treatment in CIRT [11,12,13,14]. Anatomical changes like GI deformation (gas volume changes) along the beam path may cause dose degradation, indicating that a poor dose distribution may be obtained even with TM [15]. Houweling et al. reported that using TM, the dose coverage was reduced by 10% on average and 17% for bone matching (BM) with CIRT [9]. Another study used the diaphragm structural matching to reduce positioning errors; however, this method seems to be helpful to estimate the tumor position, limiting the expected effects on maintaining the dose distribution [16]. Therefore, the current image-guided RT repositioning is limited in providing a robust CIRT treatment. Some studies attempted to improve the dose distribution by modeling daily gastrointestinal organs variations [17] or optimizing the original plan (OP) (e.g., worst-case optimization) [18]. These efforts help maintain the robustness of the dose coverage; however, the single plan applied for the entire treatment course cannot fully account for the patient-specific inter-fractional organ deformation in pancreatic cancer [19]. The full-scale re-optimization based on up-to-date imaging data seems necessary.

Adaptive radiotherapy has been considered as an effective technology for the treatment of pancreatic cancer in photon RT [8, 19]. However, only few studies have been conducted on adaptive CIRT due to various limitations [20]. Moreover, the passive scattering CIRT has a bottleneck that a patient-specific range compensator (RC) is needed at pretreatment making the performance of adaptive CIRT challenging. In addition, due to the use of single field irradiation on each treatment day, the evaluation methods for adaptive planning in photon RT seem not applicable to CIRT. There has been no study to address this issue. Therefore, this study aimed to explore the feasibility and possible benefits of employing adaptive CIRT based on single beam optimization. We proposed two adaptive plan protocols based on passive scattering CIRT for eight patients with pancreatic cancer. Dose deviations in respect to the OP were evaluated in both daily and accumulated dose distributions, as well as those investigated with BM and TM were used for comparison.


Patient and CT imaging acquisition

We studied eight consecutive patients with pancreatic cancer who were treated with passive scattering CIRT at our facility between March 2018 and February 2019. Patient characteristics can be found in Table 1. To obtain a satisfactory dose distribution by multiple beams, planning CT imaging was performed at the end of expiration in both the supine and prone positions because the fixed beam line is used at our facility (Aquilion LB, Self-Propelled, Canon Medical Systems, Japan) [15]. On the treatment day, daily CT sets were acquired by a separate CT under the same settings and conditions as the planning CT in the treatment room. Thus, a total of 16 planning and 95 daily CT sets were obtained (one CT set was missing because of the equipment failure). We performed gated irradiation within 30% of gating level.

Table 1 Patient characteristics

Original treatment plan

The treatment planning was established using a XiO-N system that employs a pencil beam algorithm (Elekta Sweden, Mitsubishi Electric, Japan) [21]. Four broad beams (angles: 0° (anterior–posterior), 90° (left–right), 270° (right-left), and 180° (posterior-anterior) based on a passive scattering technique were applied (180° beam was not applied for patient 5). Only one beam was delivered per treatment day as follows: 0°, 90°, 270°, 0°, 90°, 270°, 0°, 90°, 270°, 180°, 180°, and 180°. Our facility uses Gy (RBE) as the unit of the clinical dose, which was calculated based on the physical dose and relative biological effectiveness (RBE) [5]. The prescription dose was 55.2 Gy (RBE) in 12 fractions. The gross tumor volume (GTV) was determined for the primary tumor using a contrast-enhanced CT acquired at pretreatment. To protect the GI organs, a planning organ at risk volume (PRV)-GI defined as (GI boundary + 2 mm)—GTV was applied at our facility. In the first nine fractions, the clinical target volume 1 (CTV1) was defined as the GTV plus a 5-mm margin, including locoregional lymph nodes and neural plexus regions minus the PRV-GI margin; CTV2 for the last three boost fractions was generated by adding a 5-mm margin to the GTV, excluding the PRV-GI margin. Planning target volume (PTV) was established by adding a 3-mm margin to the CTV, excluding the PRV-GI margin. Treatment plans were performed to deliver at least 95% of the prescribed dose to 95% of the CTV and GTV (V95 > 95%). The dose coverage was slightly compromised when the distance between PTV and organs at risk (OARs) were close or overlapping. The dose to the most exposed 2 cc (D2cc) of the duodenum and stomach should receive < 44 Gy (RBE). A volume of the duodenum and stomach receiving 30 Gy (RBE) (V30 [Gy (RBE)]) should be < 10 cc as possible.

Adaptive plan

The contours were delineated on each daily CT sets by two experienced radiation oncologists. We proposed two simulation adaptive plans as AP-1 and AP-2. AP-1 was generated under an idealized condition of no RC limitation, i.e., a new RC was applied for each beam according to daily anatomical changes; AP-2 used the same RC applied in each OP beam. The same beam angle arrangement as the OP was applied in both adaptive plans. In AP-1, optimal spread-out Bragg peak (SOBP) size and range shifter were chosen in XiO-N calculation ensuring a conformal dose at both proximal and distal edges. In AP-2, both SOBP and the range shifter were manually adjusted to satisfy the criteria ensuring the target dose reduction within 3% of the OP beam. Normally, it is necessary to reduce the OAR dose while maintaining the target dose as much as possible; however, it is difficult to generate a plan that meets these requirements in pancreatic cancer. In our proposed method, the priority strategy was to ensure target coverage. Multi-leaf collimator was optimized for both adaptive plans.

Dose accumulation and evaluation parameters

To examine the dosimetric difference between the conventional positioning and adaptive plans, dose distributions with BM and TM were calculated. First, daily CT images were transferred to XiO-N for dose calculation (with the same parameters as the OP) after a rigid registration with the corresponding planning CT images with BM and TM using MIM Maestro (MIM Software, Cleveland, OH, USA). Then, dose distributions were warped and transferred to planning CT images in supine position using the DIR matrix (for one missing CT set, images acquired on the most adjacent treatment day were substituted). To reduce the DIR uncertainties, contour and intensity-based DIR (hybrid-DIR) was used [22, 23]. Lastly, accumulated doses were calculated with addition of twelve dose distributions. The accumulated doses for adaptive plans were obtained with the same process above after re-planning on each daily CT image in XiO-N.

The volume percentage receiving 95%, 90% of prescribed dose (V95, V90), the dose covered by 98% of volume (D98), and the mean dose were used to evaluate dose distributions for CTV and GTV. D2cc, mean dose, V10 [Gy (RBE)], V20 [Gy (RBE)], and V30 [Gy (RBE)] were evaluated for the duodenum and stomach. The dosimetric differences between proposed methods and OP were examined. Moreover, the dose reduction of > 5% in CTV V95 was defined as unacceptable. The significance of statistical difference between methods was analyzed using the Wilcoxson signed-rank test and the Bonferoni correction was used to adjust the P value. The Friedman test was used for multiple comparisons. P < 0.05 was considered significant.


Dosimetric difference in daily dose

An example for dose distributions in four beam angles with AP-1, AP-2, TM, and BM is shown in Fig. 1. Using BM and TM, 56.8% and 47.5% of 96 fractional dose reductions were unacceptable in CTV, respectively, in contrast with 1% in AP-1 and AP-2. A significant difference was observed between methods (P < 0.001) (Fig. 2a). The dosimetric difference of mean ± standard deviation in CTV V95 was 1.2 ± 3.4, − 0.5 ± 2.1, − 5.1 ± 6.2, and − 8.8 ± 8.8 for AP-1, AP-2, TM and BM, respectively. In AP-2, four fractional CTV doses (one in 270° beam, three in 90° beam) failed to be meet the dose constraint because of severe dose distortion caused by anatomical changes; two fractional doses (both in 90° beam) were decreased by > 3% in AP-1 due to an extreme narrow distance between tumor and OARs. Adaptive plans had the best performance in average of CTV V95 at 180° and worst at 90°, but superior to BM and TM at all beam angles (Fig. 3).

Fig. 1

Example of dose distributions using AP-1, AP-2, TM, and BM (patient 1). The left panel shows dose distributions of the original plan at 0°, 90°, 270°, 180° beam angle, and all beams. The right panels from the top to the bottom show dose distributions with AP-1, AP-2, TM, and BM in the first, second, third, tenth, and all fractions, respectively. TM, tumor matching; BM, bone matching. Filled red, gross tumor volume; filled green, clinical target volume; filled blue, kidney; filled purple, stomach; filled white, duodenum; filled yellow, spinal cord

Fig. 2

The mean dose difference of all fractions in CTV V95 (a) and GTV V95 (b) with AP-1, AP-2, TM, and BM. These differences are in respect to the original plan. GTV, gross tumor volume; CTV, clinical target volume; ns, P > 0.05. A significant difference was obtained between methods without ns (P < 0.05)

Fig. 3

The mean dose difference of all fractions in CTV V95 at each beam angle. CTV1 doses were based on the first nine fractions. These differences are in respect to the original plan. CTV, clinical target volume; ns, P > 0.05. A significant difference was obtained between methods without ns (P < 0.05)

Dosimetric difference in accumulated dose

For accumulated dose in CTV2 V95, the dose reduction of > 5% was observed in 3 and 6 patients with TM and BM, respectively; however, all patient doses were acceptable with AP-1 and AP-2 (Fig. 4). The dosimetric change in CTV2 V95 (%) was 1.2 ± 2.8, − 2.1 ± 1.7, − 7.1 ± 5.2, and − 16.5 ± 15.0 for AP-1, AP-2, TM, and BM, respectively. A better accumulated dose was obtained in CTV1 for adaptive plans. All dose reductions in CTV1 V95 (%) were < 3% in both AP-1 and AP-2. However, accumulated doses of 6 and 7 patients were unacceptable with TM and BM, respectively. The dosimetric difference in CTV1 V95 (%) was 0.1 ± 1.1, − 1.3 ± 1.4, − 11.3 ± 8.0, and − 14.0 ± 9.5 for AP-1, AP-2, TM, and BM, respectively. Figure 5 shows dosimetric differences in both daily and accumulated doses for each patient.

Fig. 4

The mean dose difference of accumulated dose in CTV1 V95 (a), CTV2 V95 (b), and GTV V95 (c) with AP-1, AP-2, TM, and BM. CTV1 doses were based on the first nine fractions. These differences are in respect to the original plan. TM, tumor matching; BM, bone matching; GTV, gross tumor volume; CTV, clinical target volume

Fig. 5

The dose difference of CTV V95 in each patient. CTV1 doses were based on the first nine fractions. These differences are in respect to the original plan. CTV, clinical target volume; ns, P > 0.05. A significant difference was obtained between methods without ns (P < 0.05)

Dosimetric difference in OARs

Compared with the OP, the median duodenum doses in D2cc and mean doses were increased in all the methods, especially AP-2 (Table 2). In one case (patient 5), the dose exceeded the tolerance in D2cc when TM, BM, and AP-2 were applied. Lower median doses were obtained in V10 [Gy (RBE)] but higher doses were observed in V20 [Gy (RBE)] and V30 [Gy (RBE)]. The stomach doses in all the cases were within the tolerance in all the methods, and no significant differences were obtained as compared with the OP (P < 0.05). The average differences of the daily dose-volume histograms for the duodenum and stomach at each beam angle between the adaptive methods are shown in Fig. 6. Compared with AP-1, AP-2 caused a higher dose in the duodenum, especially at 90° beam angle; however, the delivered dose to the stomach was comparable at all angles between the two adaptive plans. The details of the dose-volume parameters are shown in Table 2.

Table 2 The dose-volume parameters of accumulated dose
Fig. 6

The average dose-volume histogram (solid line) of all fractions at each beam angle in the duodenum (upper) and stomach (below). Dashed lines show the standard deviation


In the present study, even performing TM, the dose distribution was decreased by 7% with respect to the OP, which was consistent with that of the previous studies [9, 15]. Although performing adaptive CIRT with a new RC for each fraction based on passive irradiation technology is impractical, AP-1 was proposed to examine the distance between doses with ideal condition and actual conditions. Only AP-1 allowed for a dose increase in average fractional dose changes in respect to the OP. AP-2 had a slight reduction, but prominently superior to BM and TM, in which approximately half of fractional dose reductions were > 5% in CTV V95 (Fig. 2). This finding indicates that adaption planning might be an effective strategy to improve the target dose, especially when employing the scanning technology where the doses would be similar to those with AP-1. However, in view of the actual treatment condition, AP-2 will be discussed in particular below.

Moriya et al. verified the effectiveness of the adaptive plan with range optimization method for three patients with lung cancer and one with abdominal lymph node metastases receiving passive scattering proton [24, 25]; one highlight of their study was that the adaptive plans were calculated using a developed range-optimization system, which can reduce labor burden. However, approximately 30 min are required for optimization calculation (for six beams), which seems still too long for online re-planning. In addition, only performing range optimization seems not enough in cases with large anatomical changes, e.g., the SOBP size proposed in OP may be too large or small to account for these changes. Therefore, full-scale re-optimization is crucial. Moreover, the same dose constraint as that of OP was used in their study to evaluate the adaptive plan. This strategy is similar to that used in photon RT, where doses based on all beams can be evaluated using the rotating gantry application [8, 26]. This is a conservative method that ensures a safe dose in OARs assuming that dose distribution in each fraction is similar. However, in single beam optimization, this strategy may not always work because the positional relationship between the target and OAR changes largely, especially in case of pancreatic cancer. Hence, we developed a method using the dose parameters of each beam in OP as a reference, and a dose reduction of ≤ 3% was used as the dose constraint for the target. Although XiO-N cannot automatically select the parameters that satisfy the dose constraint because it is not an inverse planning system, the algorithm to realize this system is simple. In fact, we manually performed the aforementioned optimization in only 10–15 min.

Interestingly, although only the target dose constraint was prioritized in this study, our strategy seems effective to balance the dose distribution between the target and OARs. The accumulated dose reductions in all patients were < 5% in CTV2 and < 3% in CTV1 in both adaptive plans; only one duodenal dose was clinically unacceptable in AP-2 (Fig. 4 and Table 2). Compared with BM and TM, re-planning proposal in this study maintained superior target dose coverage while ensuring a safe OAR dose in most patients (P < 0.05) (Table 2). One thing worth mentioning is that the extent of dose changes in CTV2 and CTV1 may be different. For example, in patient 5, an acceptable dose reduction (− 3.0%) in CTV2 was obtained with TM; however, a prominent dose degradation in CTV1 was observed (− 16.6%) (Fig. 5). Therefore, just evaluating the CTV1 or CTV2 dose was inadequate.

Comparing with AP-1 and OP, AP-2 had a comparable delivered dose in the stomach but a higher dose in the duodenum (Table 2). The fact is that the dose distortion caused by anatomical changes in BM and TM also existed in AP-2, especially in the 90° beam where large intestinal deformation usually occurred (Fig. 6). This is a limitation of AP-2. To address this problem, other beam angles such as posterior oblique direction seems a reasonable option. Beam angles from the 135° to 210° are considered to be effective to obtain a robust plan accounting for interfractional deviations by avoiding passing through the bowel gas [27,28,29]. However, considering the poor flexibility of the beam-angle arrangement (limited angles are available) in the fixed-beam port system, oblique beam selection remains challenging but merits close attention in the future study.

This study has some limitations. First, the contours were delineated manually, which is time-consuming and impractical for online adaptive RT. Developing a fast and accurate contouring system seems necessary, which will be considered in our future work. Second, this study was conducted in an idealized manner that daily adaption was performed to evaluate the accumulated dose. This is quite different from other adaptive protocols used in photon RT where re-planning is only performed for unacceptable cases examined prior to treatment [30]. However, considering the high ratio of unacceptable fractional dose reduction (approximately 50% in TM and BM in this study) in CIRT, daily re-planning seems necessary to ensure a safe treatment.

Third, the impact of intrafractional organ motion on dose distribution was not considered in this study [10]. The National Institute of Radiological Sciences (NIRS) has proposed a field-specific target volume (FTV) to account for the intrafractional variation in pancreatic cancer [27]. However, the FTV is just effective with the condition that the respiratory motion and organ deformation are reproducible during treatment, and possible interfraction deviations are not considered. Magnetic resonance (MR)-guided plan adaption seems a promising and ideal technology to consider both inter- and intrafractional deviations for pancreatic cancer [31,32,33], which is expected to be available in the particle therapy in a few years.

Another limitation is that a dose-constraint OAR value for each fraction is lacking, which may cause over irradiation to OARs when the whole treatment is completed. In this study, although dose constraints of OARs were assigned a lower priority than the target coverage, to minimize the OAR dose, such adaptive plan would be selected for AP-2: the one with an acceptable dose coverage in the target but minimal dose distribution in OARs. However, one patient’s duodenum dose was not tolerated in this study. Therefore, further studies considering both target and OARs is warranted to improve the adaptive CIRT. Lastly, DIR may cause uncertainties in dose accumulation. Kubota et al. [22] quantified the errors in hDIR and showed satisfactory dice similarity coefficients (DSCs) in CTV (0.94 ± 0.05); worse DSCs in the stomach and duodenum (0.85 ± 0.09 and 0.81 ± 0.06, respectively). These correspond to a 1.50% difference in median CTV V95 and 2.46% and 0.68% differences in median stomach and duodenum V50, respectively. Although, these errors seem within tolerance, a validation study with large samples is needed and improving the accuracy of DIR is necessary in further study.


Improving the dose distribution for pancreatic cancer in CIRT with current image-guided positioning strategies is challenging. The possible dose reduction should be carefully considered when employing BM, even TM. The dose coverage in the target was significantly improved by re-planning in this study. Although AP-2 has uncertainties in OARs sparing due to RC limitation, the dose reduction constraint proposed for target seems effective to minimize OARs dose to some extent. However, considering a prominent dose increase in the 90° beam in the duodenum, other beam angle arrangements should be examined to improve this technology in the future. In fact, many issues still need to be addressed including the software and hardware aspects before performing adaptive CIRT in clinical practice, such as adjusting the range shifter and/or the ridge filter online is not yet possible at our facility, but it is not technically difficult. However, adaptive methods based on single beam optimization proposed in this study provide important evidence for the feasibility of implementing plan adaption with passive scattering CIRT in patients with pancreatic cancer.

Availability of data and materials

The datasets analyzed during the current study are available from the corresponding author on reasonable request.



carbon ion radiotherapy


computed tomography


range compensator


tumor matching


bone matching




original plan


planning target volume


clinical target volume


gross tumor volume


relative biological effectiveness


planning organ at risk volume


spread out Bragg peak


organs at risk


field-specific target volume


satisfactory dice similarity coefficients


  1. 1.

    Kawashiro S, Yamada S, Okamoto M, Ohno T, Nakano T, Shinoto M, et al. Multi-institutional study of carbon-ion radiotherapy for locally advanced pancreatic cancer: Japan Carbon-ion Radiation Oncology Study Group (J-CROS) Study 1403 Pancreas. Int J Radiat Oncol Biol Phys. 2018;101:1212–21.

    Article  PubMed  Google Scholar 

  2. 2.

    Shinoto M, Yamada S, Terashima K, Yasuda S, Shioyama Y, Honda H, et al. Working Group for Pancreas Cancer. Carbon ion radiation therapy with concurrent gemcitabine for patients with locally advanced pancreatic cancer. Int J Radiat Oncol Biol Phys. 2016;95:498–504.

    CAS  Article  PubMed  Google Scholar 

  3. 3.

    Kawashiro S, Mori S, Yamada S, Miki K, Nemoto K, Tsuji H, et al. Dose escalation study with respiratory-gated carbon-ion scanning radiotherapy using a simultaneous integrated boost for pancreatic cancer: simulation with four-dimensional computed tomography. Br J Radiol. 2017;90:20160790.

    Article  PubMed  PubMed Central  Google Scholar 

  4. 4.

    Okamoto M, Shiba S, Okazaki S, Miyasaka Y, Shibuya K, Kiyohara H, et al. Feasibility and safety of repeated carbon ion radiotherapy for locally advanced unresectable pancreatic cancer. Cancers. 2021;13:665.

    Article  PubMed  PubMed Central  Google Scholar 

  5. 5.

    Ohno T. Particle radiotherapy with carbon ion beams. EPMA J. 2013;4:1–7.

    Article  Google Scholar 

  6. 6.

    Kanai T, Endo M, Minohara S, Miyahara N, Koyama-ito H, Tomura H, et al. Biophysical characteristics of HIMAC clinical irradiation system for heavy-ion radiation therapy. Int J Radiat Oncol Biol Phys. 1999;44:201–10.

    CAS  Article  PubMed  Google Scholar 

  7. 7.

    van der Horst A, Wognum S, Dávila Fajardo R, de Jong R, van Hooft JE, Fockens P, et al. In-terfractional position variation of pancreatic tumors quantified using intratumoral fiducial markers and daily cone beam computed tomography. Int J Radiat Oncol Biol Phys. 2013;87:202–8.

    Article  PubMed  Google Scholar 

  8. 8.

    Liu F, Erickson B, Peng C, Li XA. Characterization and management of interfractional anatomic changes for pancreatic cancer radiotherapy. Int J Radiat Oncol Biol Phys. 2012;83:e423–9.

    Article  PubMed  Google Scholar 

  9. 9.

    Houweling AC, Fukata K, Kubota Y, Shimada H, Rasch CR, Ohno T, et al. The impact of inter-fractional anatomical changes on the accumulated dose in carbon ion therapy of pancreatic cancer patients. Radiother Oncol. 2016;119:319–25.

    Article  PubMed  Google Scholar 

  10. 10.

    Kumagai M, Hara R, Mori S, Yanagi T, Asakura H, Kishimoto R, et al. Impact of intrafractional bowel gas movement on carbon ion beam dose distribution in pancreatic radiotherapy. Int J Radiat Oncol Biol Phys. 2009;73:1276–81.

    CAS  Article  PubMed  Google Scholar 

  11. 11.

    Li Y, Kubota Y, Kubo N, Mizukami T, Sakai M, Kawamura H, et al. Dose assessment for patients with stage I non-small cell lung cancer receiving passive scattering carbon-ion radiotherapy using daily computed tomographic images: a prospective study. Radiother Oncol. 2020;144:224–30.

    CAS  Article  PubMed  Google Scholar 

  12. 12.

    Sakai M, Kubota Y, Saitoh JI, Irie D, Shirai K, Okada R, et al. Robustness of patient positioning for interfractional error in carbon ion radiotherapy for stage I lung cancer: bone matching versus tumor matching. Radiother Oncol. 2017;129:95–100.

    Article  Google Scholar 

  13. 13.

    Kubota Y, Katoh H, Shibuya K, Shiba S, Abe S, Sakai M, et al. Comparison between bone matching and marker matching for evaluation of intra- and inter-fractional changes in accumulated of carbon ion radiotherapy for hepatocellular carcinoma. Radiother Oncol. 2019;137:77–82.

    Article  PubMed  Google Scholar 

  14. 14.

    Abe S, Kubota Y, Shibuya K, Koyama Y, Abe T, Ohno T, et al. Fiducial marker matching versus vertebral body matching: dosimetric impact of patient positioning in carbon ion radiotherapy for primary hepatic cancer. Phys Med. 2017;33:114–20.

    Article  PubMed  Google Scholar 

  15. 15.

    Kubota Y, Okamoto M, Shiba S, Okazaki S, Matsui T, Li Y, et al. Robustness of daily dose for each beam angle and accumulated dose for inter-fractional anatomical changes in passive carbon-ion radiotherapy for pancreatic cancer: bone matching versus tumor matching. Radiother Oncol. 2021;157:85–92.

    CAS  Article  PubMed  Google Scholar 

  16. 16.

    Itabashi Y, Kubota Y, Okamoto M, Tsuda K, Shiba S, Hoshino Y, et al. Potential pitfalls of diaphragm structural matching in carbon-ion radiotherapy for pancreatic cancer. Anticancer Res. 2019;39:4351–6.

    CAS  Article  PubMed  Google Scholar 

  17. 17.

    Magallon-Baro A, Loi M, Milder MTW, Granton PV, Zolnay AG, Nuyttens JJ, et al. Modeling daily changes in organ-at-risk anatomy in a cohort of pancreatic cancer patients. Radiother Oncol. 2019;134:127–34.

    Article  PubMed  Google Scholar 

  18. 18.

    Steitz J, Naumann P, Ulrich S, Haefner MF, Sterzing F, Oelfke U, et al. Worst case optimization for interfractional motion mitigation in carbon ion therapy of pancreatic cancer. Radiat Oncol. 2016;11:134.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  19. 19.

    Ahunbay EE, Kimura B, Liu F, Erickson BA, Li XA. Comparison of various online strategies to account for inter-fractional variations for pancreatic cancer. Int J Radiat Oncol Biol Phys. 2013;86:914–21.

    Article  PubMed  Google Scholar 

  20. 20.

    Li Y, Kubota Y, Tashiro M, Ohno T. Value of three-dimensional imaging systems for image-guided carbon ion radiotherapy. Cancers (Basel). 2019;11:297.

    CAS  Article  Google Scholar 

  21. 21.

    Kanematsu N. Dose calculation algorithm of fast fine-heterogeneity correction for heavy charged particle radiotherapy. Phys Med. 2011;27:97–102.

    Article  PubMed  Google Scholar 

  22. 22.

    Kubota Y, Okamoto M, Li Y, Shiba S, Okazaki S, Komatsu S, et al. Evaluation of intensity- and contour-based deformable image registration accuracy in pancreatic cancer patients. Cancers (Basel). 2019;11:1447.

    CAS  Article  Google Scholar 

  23. 23.

    Li Y, Kubota Y, Okamoto M, Shiba S, Okazaki S, Matsui T, et al. Determination of deformable image registration algorithms for accumulating dose in carbon-ion radiotherapy for pancreatic cancer. Anticancer Res. 2021;41:835–43.

    CAS  Article  PubMed  Google Scholar 

  24. 24.

    Moriya S, Tachibana H, Hotta K, Nakamura N, Sakae T, Akimoto T. Range optimization for target and organs at risk in dynamic adaptive passive scattering proton beam therapy—a proof of concept. Phys Med. 2018;56:66–73.

    Article  PubMed  Google Scholar 

  25. 25.

    Moriya S, Tachibana H, Hotta K, Nakamura N, Sakae T, Akimoto T. Feasibility of dynamic adaptive passive scattering proton therapy with computed tomography image guidance in the lung. Med Phys. 2017;44:4474–81.

    Article  PubMed  Google Scholar 

  26. 26.

    Boldrini L, Cusumano D, Cellini F, Azario L, Mattiucci GC, Valentini V. Online adaptive magnetic resonance guided radiotherapy for pancreatic cancer: state of the art, pearls and pitfalls. Radiat Oncol. 2019;14:71.

    Article  PubMed  PubMed Central  Google Scholar 

  27. 27.

    Koom WS, Mori S, Furuich W, Yamada S. Beam direction arrangement using a superconducting rotating gantry in carbon ion treatment for pancreatic cancer. Br J Radiol. 2019;92:20190101.

    Article  PubMed  PubMed Central  Google Scholar 

  28. 28.

    Shiomi M, Mori S, Shinoto M, Nakayama Y, Kamada T, Yamada S. Comparison of carbon-ion passive and scanning irradiation for pancreatic cancer. Radiother Oncol. 2016;119:326–30.

    Article  PubMed  Google Scholar 

  29. 29.

    Batista V, Richter D, Combs SE, Jäkel O. Planning strategies for inter-fractional robustness in pancreatic patients treated with scanned carbon therapy. Radiat Oncol. 2017;12:94.

    Article  PubMed  PubMed Central  Google Scholar 

  30. 30.

    Lim-Reinders S, Keller BM, Al-Ward S, Sahgal A, Kim A. Online adaptive radiation therapy. Int J Radiat Oncol Biol Phys. 2017;99:994–1003.

    Article  PubMed  Google Scholar 

  31. 31.

    Lagerwaard F, Bohoudi O, Tetar S, Admiraal MA, Rosario TS. Combined inter- and intrafractional plan adaptation using fraction partitioning in magnetic resonance-guided radiotherapy delivery. Cureus. 2018;10: e2434.

    Article  PubMed  PubMed Central  Google Scholar 

  32. 32.

    Bohoudi O, Bruynzeel AME, Senan S, Cuijpers JP, Slotman BJ, Lagerwaard FJ, et al. Fast and robust online adaptive planning in stereotactic MR-guided adaptive radiation therapy (SMART) for pancreatic cancer. Radiother Oncol. 2017;125:439–44.

    CAS  Article  PubMed  Google Scholar 

  33. 33.

    Corradini S, Alongi F, Andratschke N, Belka C, Boldrini L, Cellini F, et al. MR-guidance in clinical reality: current treatment challenges and future perspectives. Radiat Oncol. 2019;14:92.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

Download references


The authors thank the radiology technologists, medical doctors, nurses, and medical physicists at GHMC for their valuable insights.


This research received no external funding.

Author information




Conceptualization, YK; methodology, YK, YL, MT and MO; software, YK, YL; validation, YK and YL; formal analysis, YL; investigation, YK, YL, MO, SS, SO, and TM; resources, YK, YL; data curation, YK and YL; writing—original draft preparation, YL; writing—review and editing, YL, YK, SS, SO, TM, TO and MT; visualization, YL and YK; supervision, YK, MO, TN and TO; project administration, YK, MO, and TO. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Yoshiki Kubota.

Ethics declarations

Ethics approval and consent to participate

This study was approved by our facility’s institutional review board and was registered at the University Hospital Medical Information Network Clinical Trials Registry (UMIN-CTR trial number: 000029495). All patients gave written informed consent for inclusion before they participated in the study.

Consent for publication

Informed consent for publication was collected from all the patients.

Competing interests

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Li, Y., Kubota, Y., Okamoto, M. et al. Adaptive planning based on single beam optimization in passive scattering carbon ion radiotherapy for patients with pancreatic cancer. Radiat Oncol 16, 111 (2021).

Download citation


  • Carbon-ion radiotherapy
  • Adaptive planning
  • Pancreatic cancer
  • Accumulated dose assessment
  • Robustness of treatment