Skip to main content

Comparison of photon volumetric modulated arc therapy, intensity-modulated proton therapy, and intensity-modulated carbon ion therapy for delivery of hypo-fractionated thoracic radiotherapy



The aim of the present study was to compare the dose distribution generated from photon volumetric modulated arc therapy (VMAT), intensity modulated proton therapy (IMPT), and intensity modulated carbon ion therapy (IMCIT) in the delivery of hypo-fractionated thoracic radiotherapy.

Methods and materials

Ten selected patients who underwent thoracic particle therapy between 2015 and 2016 were re-planned to receive a relative biological effectiveness (RBE) weighted dose of 60 Gy (i.e., GyE) in 15 fractions delivered with VMAT, IMPT, or IMCIT with the same optimization criteria. Treatment plans were then compared.


There were no significant differences in target volume dose coverage or dose conformity, except improved D95 was found with IMCIT compared with VMAT (p = 0.01), and IMCIT was significantly better than IMPT in all target volume dose parameters. Particle therapy led to more prominent lung sparing at low doses, and this result was most prominent with IMCIT (p < 0.05). Improved sparing of other thoracic organs at risk (OARs) was observed with particle therapy, and IMCIT further lowered the D1cc and D5cc for major blood vessels, as compared with IMPT (p = 0.01).


Although it was comparable to VMAT, IMCIT led to significantly better tumor target dose coverage and conformity than did IMPT. Particle therapy, compared with VMAT, improved thoracic OAR sparing. IMCIT, compared with IMPT, may further improve normal lung and major blood vessel sparing under limited respiratory motion.


Early-stage non-small cell lung cancer (NSCLC) and lung metastases have been treated with stereotactic body radiation therapy (SBRT) with excellent clinical outcome [1,2,3,4]. Locally advanced NSCLC has traditionally been treated with conventionally fractionated concurrent chemo-radiation [5]. However, dose escalation with conventional fractionation has not demonstrated any clinical advantage in any phase III randomized controlled trials [6]. As a result, alternative radio-therapeutic strategies are being sought. Given the clinical efficacy that has been observed with hypo-fractionated dose schedules, such as SBRT, this approach has also been increasingly considered as a treatment strategy for locally advanced and selected stage IV NSCLC in recent years [7,8,9].

Excellent dose distribution can be achieved with image-guided and intensity-modulated radiotherapy (IG-IMRT) delivered with volumetric modulated arc therapy (VMAT) in the thorax [10, 11]. As a result, VMAT has been quickly adopted in clinical settings to deliver thoracic radiotherapy [12]. However, sparing thoracic organs at risk (OARs) while maintaining adequate target volume dose coverage with VMAT remains challenging in selected patients. On the contrary, particle therapy (PT), which includes proton and heavy ion therapies, may have an advantage over VMAT in OAR sparing. This improvement is due to PT’s physical properties, which allow for better normal tissue protection, while heavy ions’ increased radiobiological effectiveness (RBE) increases the tumorcidal effect of radiotherapy over photons [13, 14]. PT’s dosimetric advantages have been shown in multiple studies [15,16,17,18,19,20,21,22]. Given these advantages, PT is increasingly being considered for the delivery of hypo-fractionated thoracic radiotherapy [23, 24]. To date, PT has been mostly delivered with passive scattering systems. With further improvements in technology, intensity modulation of particle beams with active beam scanning has been developed for better dose conformity and control of the OAR dose. However, little is known about how active scanning beams compare with VMAT in dose distribution in the treatment of lung tumors. In this study, we explored the dosimetric characteristics of intensity-modulated proton therapy (IMPT) and carbon ion therapy (IMCIT) with raster scanning beams in the delivery of hypo-fractionated thoracic radiotherapy in comparison to VMAT under limited respiratory motion.


Patients, tumors and treatment characteristics

Ten selected patients who underwent either IMPT or IMCIT between 2015 and 2016 were re-planned. Among them, 9 patients had stage I - III NSCLC, and 1 patient had a single lung metastasis from adenoid cystic carcinoma of the head and neck (Table 1). All patients were simulated head first in prone or supine positions. Each patient was immobilized with Vac-Lok™ cushions (CIVCO Medical Solutions, IA, USA). During CT simulation, free-breathing CT and 4D CT (divided evenly over 10 phases of the respiratory cycle) were obtained with 3-mm slice thickness to account for respiratory motion. All patients were planned for gated treatments with a pressure sensor-based motion management system (AZ-773 V, Anzai Medical, Co., Ltd., Tokyo, Japan). The gating window (30–40% duty cycle) was selected to keep the respiratory motion to no more than 5 mm cranio-caudally (Table 1). This motion limitation was accomplished through visual inspection of the 4D CT's to select the phases from the end of expiration to the beginning of inspiration that were confined within a cranio-caudal tumor motion amplitude of ≤5 mm. The primary gross tumor volume (GTV) was delineated at the lung window level, and any nodal metastases and soft tissue extensions were delineated at the soft tissue window level on the maximum exhalation CT. To account for tumor motion, an iGTV was created with the maximal intensity projection (MIP) generated from CTs of the selected respiratory phases. The clinical target volume (CTV) was the iGTV with an 8-mm expansion. The planning target volume (PTV) was the CTV with a 5-mm expansion to account for set up errors and residual tumor motion. The CTV and PTV were created on the average CT of the selected respiratory phases. Particular attention was paid to avoid overlapping of any target volumes with the OARs. The lungs, esophagus, spinal cord, and the heart were contoured for each patient on the average CT. The major blood vessels and major airway were contoured only when they were adjacent to the PTV. Target and OAR volume delineation was performed with MIM v6.5.9 software (MIM software Inc., Cleveland, OH, USA).

Table 1 Patient, tumor, and treatment characteristics

All plans prescribed 60 Gy/GyE (relative biological effectiveness, RBE, weighted dose for IMPT and IMCIT) delivered in 15 daily fractions to the PTV with heterogeneity corrections. For VMAT, a 1.0-Gy physical dose is identical to 1.0 GyE. For IMPT, a 1.0-Gy physical dose is equal to a constant factor of 1.1 GyE at all positions. However, a 1.0-Gy physical dose for IMCIT was scaled to biological doses according to the LEM-1 model. Raster-continuous spot scanning without stopping between spots was delivered with dynamic intensity control to decrease the time required to deliver large doses (resulting in a maximum of 1–2 ms per spot). In our beam delivery system, the beam intensity could be dynamically changed between spots according to the total number of particles required for each spot. The change in beam intensity could be performed within 100 μs. Because our beam delivery system was calibrated to the number of particles per unit dose Gy, 13 levels of intensity for both carbon-ion and proton beams could be used clinically (1.3 × 106 to 6.5 × 107/s for carbon ions and 5.0 × 107 to 2.6 × 109/s for protons). However, only the final 6 highest levels (within ¼ of the maximum intensities) were used for actual treatment. The average dose rate at the center of each spot during beam delivery was 2.0 to 10.0 Gy/s or 12.0 to 60.0 Gy/min for both carbon-ion and proton beams. The dose delivery time of each energy-layer was typically within 1.0 s. With a typical energy-layer switching time of approximately 4–5 s for both carbon-ion and proton beams, the doses for spots in each energy layer were typically delivered within a 5-s breathing cycle.

The plans were optimized to have ≥ 90% of the PTV receiving 100% of the prescription dose and ≥ 95% of the PTV receiving ≥ 95% of the prescription dose. OAR dose constraints were adopted from previous studies on hypo-fractionated thoracic radiotherapy [8, 23]. PTV coverage took precedence over OAR sparing in all plans, which were optimized under the same set of planning criteria.

Photon planning

VMAT plans were generated with Rapid Arc (RA) using the Eclipse Treatment Planning System V.11 (Varian Medical Systems. Palo Alto, CA, USA) with 6 MV photons. Plans were generated with co-planar partial arcs to spare as much contralateral lung as possible. The arc angle and collimator setting were based on the target size, shape, and location. Then, 2-arc plans were created for each case at dynamic multi-leaf collimator motion, variable dose rate (max: 600 MU/min), and variable gantry rotation speed to generate a modulated dose distribution. The arc length varied from 0° to 320°. All plans were optimized using the Progressive Resolution Optimizer (PRO, v11.0.31), and the volume dose was calculated with the Anisotropic Analytical Algorithm (AAA, V11.0.31) with a 2.5-mm grid size resolution.

PT planning (IMPT & IMCIT)

IMPT and IMCIT plans were generated with multi-field optimization (MFO) by using the Syngo® PT Planning System v10 (Siemens. Erlangen, Germany). The effective dose calculation model (Local Effect Model, LEM) was used for IMCIT dose calculation [25, 26], and IMPT dose calculation was performed with a constant RBE factor of 1.1. Plan optimization was done on average CT's with iGTV density overridden (to the density of muscle) [27]. All plans were implemented in an oblique (45°) beam-line room with a raster scanning delivery technique. With the fixed 45° beam-line, three different angles of the robotic couch were selected to perform non-coplanar optimization. Angle selection depended on the tumor location and size, beam path according to OAR location, and the innate restriction of the robotic couch to avoid any unnecessary dose and to evade strong density heterogeneities upon the beam’s entry. Beam directions traversing patient support and immobilization were avoided. The setting of the focal spot size for the lateral full-width half-maximum (FWHM) of the scanning beam was 8 mm, and the longitudinal beam spot range step was 3 mm. Considering the calculation time and accuracy, we used 3 mm as the dose grid resolution for both proton and carbon ion treatment plans. Finally, the virtual target expansion was specified as 3 mm in lateral, proximal, and distal directions to allow optimal placement of additional raster spots outside the PTV within the expansion volume to satisfy target coverage requirements. The estimated treatment delivery times (median) for carbon-ion and proton plans were 736 (281–1009) seconds and 514.5 (276–703.5) seconds, respectively.

Plan comparison

All plans were transferred into MIM software for direct comparison of the dose parameters obtained from the VMAT, IMPT, and IMCIT plans. OAR dose parameters, including the mean lung dose (MLD), the volume of the normal lung receiving 5 Gy, 10 Gy, and 20 Gy (V5, V10, and V20, respectively) for the total lung (both lungs – iGTV), the ipsilateral lung and the contralateral lung, the heart’s mean dose (HMD), V5, V15, and V30, the maximum dose (Dmax), doses to 0.01 cm3, and 1 cm3 (D0.01cm 3 and D1cc, respectively) of the esophagus, spinal cord, major vessels and the major airways, were compared. For the PTV, the dose covering 95% and 99% of the PTV (D95, D99), the percentage of PTV receiving ≥95% of the prescribed dose (V95), Dmin, Dmax, Dmean, the homogeneity index (HI), and the conformity index (CI) were compared among the 3 treatment techniques. The CI and HI were defined as [20]:

$$ \mathrm{CI}={\mathrm{V}}_{95\%\mathrm{Rx}}/{\mathrm{V}}_{\mathrm{PTV}} $$
$$ \mathrm{HI}={\mathrm{D}}_{\mathrm{max}}/{\mathrm{D}}_{\mathrm{Rx}} $$

where V95%Rx and VPTV represent the volume of tissue  receiving ≥ 95% of the prescribed dose and the volume of the PTV, respectively. Dmax and DRx represent the maximum dose within the PTV and the prescribed dose, respectively. The new Conformity Index (nCI) was also assessed to evaluate the degree to which the prescribed isodose volume conforms to the target volume [28].

$$ \mathrm{nCI}=\left[\left(\mathrm{treatment}\ \mathrm{volume}\right)\times \left(\mathrm{prescription}\ \mathrm{isodose}\ \mathrm{volume}\right)\right]/{\left(\mathrm{volume}\ \mathrm{of}\ \mathrm{the}\ \mathrm{target}\ \mathrm{covered}\ \mathrm{by}\ \mathrm{the}\ \mathrm{prescription}\ \mathrm{isodose}\ \mathrm{volume}\right)}^2 $$

Dose evaluation with recalculation on selected 4D phases within the gating window

Gating has been applied in recognition of the effect of the interplay on dose due to variations in the temporal relationship between dynamic beam scanning and target motion, as described in previous studies [29, 30]. By applying gating to limit the target movement to within 5 mm, dose heterogeneity was found to be within the clinically acceptable range of 3% to 5% at our institution (Additional file 1). Without the ability to fully account for this interplay, we focused mainly on the effect of uncertainties in dose due to anatomical changes within the gating window of 5 mm by retrospectively recalculating the IMCT and IMPT plans using selected 4D CT phases within the gating window. The spot pattern and doses obtained in the 3D plan were used to recalculate the doses on selected CT phases. Dose distributions of all recalculated 4D phases were accumulated to the reference CT by utilizing the deformable registration workflow provided by MIM 6.5.9. Selected dose parameters from the corresponding DVH’s were then extracted and compared with those from the original 3D plans.

Statistical analysis

The primary objective of this study was the dosimetric comparison of VMAT, IMPT and IMCIT for delivery of hypo-fractionated thoracic radiotherapy. The dose parameters were obtained from treatment plans for the three radio-therapeutic modalities for each patient and were analyzed as continuous variables. A sample size of 10 patients had 80% power to detect one standard deviation difference by using a two-sided paired t-test at a 0.05 significance level. Descriptive statistics were used to summarize patient data, including OAR dose parameters and target volume dose coverage parameters. Categorical data were described using contingency tables including counts and percentages; continuous variables were summarized with descriptive statistical measures (i.e., mean (± s.d.)). The Wilcoxon signed rank test was used to assess dose parameters or target volume dose coverage as paired data (within same patients but different radio-therapeutic modalities) without the assumption of a normal distribution. A p-value <0.05 was considered to be statistically significant. Statistical analyses were performed using SAS 9.2 and R software, version R 3.1.3 (SAS Institute Inc., Cary, NC, USA).


Target volume dose distribution

No statistically significant differences in any pertinent PTV dose coverage parameters, with the exception of D95, was observed between VMAT and IMCIT (Table 2). In contrast, statistically significant differences in all parameters favoring fixed-beam IMCIT over fixed-beam IMPT were observed (p < 0.05). Although no statistically significant differences in V95, D95, and D99 were observed between VMAT and fixed-beam IMPT, VMAT led to significantly lower Dmax, CI, HI, and nCI (p < 0.05).

Table 2 Target volume dose distribution

Dose to the normal lungs

The volumes of normal lung tissue receiving low doses were significantly smaller for PT (Additional file 2: Table S1; Fig. 1). However, the V50’s for the total and ipsilateral lung were slightly, but significantly, higher for PT than VMAT (p < 0.05). A trend toward an increased ipsilateral lung V60 was also observed with IMPT compared with VMAT (p = 0.06). PT led to significantly lower V5, V10, and V20 for the total, ipsilateral, and contralateral lung than did photon VMAT (p < 0.05, Table 3). The contralateral lung’s MLD was significantly lower with PT (p < 0.05). The MLD for the total lung was similar between VMAT, IMPT, and IMCIT, whereas a lower MLD was observed for IMCIT compared with VMAT, with a trend toward statistical significance (p = 0.07). In addition, IMCIT further decreased the V5 and V10 for the total and ipsilateral lung from that achieved with IMPT (p < 0.05). Additionally, lower ipsilateral lung V30 for IMCIT compared with IMPT exhibited a trend toward statistical significance (p = 0.09).

Fig. 1
figure 1

Dose distribution for normal lungs. a Dose distribution for the total lung (both lungs – iGTV). b Dose distribution for the ipsilateral lung. c Dose distribution for the contralateral lung

Table 3 Normal lung dose volume parameters

Dose to the heart

When compared with VMAT, IMCIT resulted in significantly lower MHD, the D1cc and the D5cc (Table 4), and heart’s V5 – V20 (p < 0.05), (Fig. 2a, Additional file 2: Table S2). Similarly, IMPT led to significantly lower MHD and V5 – V20 than VMAT (p < 0.05). However, no significant difference in heart sparing was observed between IMCIT and IMPT.

Table 4 Selected dose parameters for other OARs
Fig. 2
figure 2

Dose distribution for the heart and the esophagus. a Dose distribution for the heart. b Dose distribution for the esophagus

Dose to the esophagus

Although no significant differences between IMPT and IMCIT was found, both led to significantly improved esophageal sparing, as compared with VMAT. The mean dose, Dmax, D0.01cc, D1cc, D5cc, and the V5 - V30 were significantly lower with PT (p < 0.05). These findings are shown in Table 4, Additional file 2: Table S2, and Fig. 2b.

Dose to other thoracic organs

Dose reduction for the spinal cord was observed with PT for all dose parameters assessed (Table 4). However, no significant difference in spinal cord sparing was observed between IMPT and IMCIT. Similar findings were observed for the major blood vessels in the vicinity of the PTV. IMCIT further decreased the D1cc and D5cc to levels below those achieved with IMPT (p = 0.01). For major airway sparing, no significant difference between IMPT and IMCIT was observed, whereas PT, compared with VMAT, resulted in significantly lower mean doses, D1cc and D5cc for major airways (Table 4).

Dose evaluation with recalculation on selected 4D phases for carbon ion and proton plans

Selected dose parameters for the PTV and the common OARs after dose recalculation with dose summation from 4D phases within the gating window, in comparison to those obtained from 3D planning, are shown in Additional file 2: Table 3. No statistically significant differences in the OAR dose parameters was observed. Among the 6 dose parameters assessed for the PTV (V95, Dmax, Dmin, CI, HI, nCI), dose recalculation decreased Dmax and HI within the PTV (p < 0.05).


To the best of our knowledge, this study is the first study comparing photon VMAT with IMPT and IMCIT in the definitive treatment of loco-regionally confined lung tumors with hypo-fractionated radiotherapy. Although their clinical significance remains to be determined, our results indicated that fixed-beam IMCIT, compared with fixed-beam IMPT, may improve dose distribution for the tumor target, an effect possibly related to the sharper lateral dose gradient associated with carbon ions. The target volume dose coverage between VMAT and fixed-beam IMCIT was similar. Although the D95 was slightly improved with IMCIT compared with VMAT, the difference may be too small to have any clinical impact. However, fixed-beam IMCIT with 3 fields is greatly limited by the number and angles of beams that can be used to generate the most conformal treatment plan. In contrast, VMAT is one of the most mature and sophisticated forms of intensity-modulated photon therapy, and it can generate excellent dose conformity to the tumor target [11]. This difference indicates the great potential for IMCIT to further improve the dose distribution and the clinical efficacy of thoracic radiotherapy as it evolves in the future, given its known advantages over photons and protons (sharper lateral dose gradient, higher relative biological effectiveness, lower oxygen enhancement ratio, and more densely ionizing tracks that can lead to increased tumor DNA damage, less cell cycle dependence, and a stronger immunological response) [13, 14].

Suboptimal dose conformity due to the limitations of fixed-beam PT appears to be more prominent with IMPT, because protons, compared with carbon ions, are associated with greater lateral penumbra. Improving dose conformity for IMPT is important, because it may decrease the risk of radiation pneumonitis [31]. Although IMPT and CIRT, compared with 3D–CRT, have been shown to improve target volume coverage when multiple fields are used, IMRT has been shown to have better dose conformity and homogeneity than passively scattered PT [17,18,19,20,21,22]. Our study comparing IMPT and VMAT demonstrated similar findings. IMPT arcs may lead to better dose conformity than VMAT at the cost of increased dose heterogeneity, a finding that warrants further investigation [20].

PT has traditionally been shown to decrease the low dose volumes in the lungs [15,16,17,18,19,20,21,22]. IMPT, and especially IMPT arcs, compared with passively scattered PT, may further decrease the normal lung dose [16, 17, 19, 20]. Despite prominent low dose sparing, either passively scattered PT or IMPT may lead to higher volumes of the normal lung receiving doses >50% of the prescribed dose than photon therapy delivered with 3D–CRT or IMRT [17, 18]. However, this unique phenomenon of proton therapy has not been previously observed in comparisons of CIRT with 3D–CRT [21, 22]. Both IMPT and IMCIT significantly lowered the normal lung V5 - V20 in this study. This result was similar to those from previous studies comparing proton and photon therapies [15,16,17,18,19,20]. In contrast, VMAT, compared with either IMPT or IMCIT, appears to be associated with slightly lower but clinically non-significant volumes of normal lungs receiving high doses close to the prescribed dose, probably because of the suboptimal dose conformity associated with fixed-beam particle therapy (Fig. 1). In addition, IMCIT, compared with IMPT, may further decrease the ipsilateral lung volumes receiving low doses, thus leading to reduced total lung volumes receiving low doses. These findings may be associated with both the physical properties of particle beam therapy and the better dose conformity observed with VMAT and IMCIT compared with IMPT. However, improved normal lung sparing with IMCIT compared with VMAT was not found to be as dramatic as previously observed when passively scattered CIRT was compared with 3D–CRT, possibly because of the better control of normal lung dose associated with VMAT compared with 3D–CRT [21, 22]. Given PT’s physical properties, this sparing may be improved with arc-based methods of beam delivery [20].

Significant sparing of the heart, esophagus, and other thoracic OARs was observed with PT compared with VMAT (Table 4, Fig. 2). This finding was consistent with those of previous studies comparing PT and photon therapy [15,16,17,18,19,20,21,22]. Although no significant difference between IMCIT and IMPT was observed in the sparing of the heart, the esophagus, the spinal cord, or the major airway, the sparing of the heart and the major blood vessels in the high dose volumes may be best achieved with IMCIT. Overall, our findings suggested that although they are delivered with a limited number of beams and beams angles, fixed-beam IMPT and IMCIT may have an OAR sparing advantage over photon VMAT, which is more prominent with IMCIT than IMPT in the sparing of normal lungs and major blood vessels, owing to its improved dose conformity. Such an advantage may be further augmented as PT is more developed in the future.

Intensity modulation through active scanning may further improve dose conformity and OAR sparing in the thorax [16, 17, 19, 20]. This advantage is significantly limited by organ motion due to changes in radiological path length as a result of organ motion adjacent to the tumor, inter-field motion, and the interplay of interference between beam and tumor motion, which may result in tumor under-dose and OAR over-dose [32]. The interplay effect has been shown to significantly deteriorate target dose coverage and dose homogeneity when the motion amplitude is greater than 8 mm [33]. The correlation between motion amplitude and dose heterogeneity due to interplay, which leads to increased under-dosing within the target volume, has been further demonstrated in a 4D Monte Carlo simulation [34]. Such under-dosing has been modeled to significantly decrease the 2-year local control. However, this problem may be mitigated through fractionation, which leads to a motion-averaging effect on dose distribution, and the selection of larger spot size for the scanning beam. Together, these strategies can retain dose homogeneity with motion amplitudes of <20 mm. The mitigating effect of fractionation on interplay has also been shown in a series of 11 patients with stage III NSCLC for whom conventional fractionation with 35 fractions and hypo-fractionation with 10 fractions were modeled [35]. These observations further support our utilization of a hypo-fractionated schedule in the current study.

Respiratory motion is mainly managed through rescanning, gating, and beam tracking [32, 36, 37]. Rescanning mitigates the interplay effect by averaging the under-dosing and over-dosing patterns of the dose distribution. This process can be accomplished through volumetric or energy-slice by energy-slice rescanning while motion parameters are kept different between rescans. Rescanning mitigates the interplay effects without mitigating the target motion. Thus, adequate margins need to be maintained for adequate target dose coverage, possibly leading to suboptimal OAR sparing in conjunction with the dose blurring at the field borders. In addition, rescanning may be beneficial only when significant interplay exists under moderate to large respiratory motion [38]. Gating, which is already in clinical use for photon therapy to minimize respiratory motion-related dose degradation and unnecessary OAR irradiation, has been a major approach of respiratory motion management for particle therapy to minimize interplay. Through 4D dose calculation, adequate target volume dose coverage and dose homogeneity have been demonstrated with a gating window (GW) of ≤5 mm, and lung dose can be further decreased within the clinically acceptable range with shorter GWs up to 1 mm for spot scanning proton therapy [39]. This process supports our motion selection criteria of limiting the cranio-caudal motion to ≤5 mm through selecting the motion phases, including the maximal expiratory phase, by visual inspection of each patient’s 4D CT. Although 4D dose calculation to better account for interplay was not performed in this study, dose heterogeneities due to interplay for target motion of ≤5 mm were found to be within the clinically acceptable 5% for both carbon ions and protons at our institution (Additional file 1). How to best account for interplay in thoracic particle therapy will be more thoroughly assessed in future studies. A limitation for evaluating the effects of range uncertainties and interplay in this study is that errors may exist in the correlation between surrogate marker motion and internal tumor motion, and irregularities in the breathing pattern may also be present during real-time treatment. Such errors may be further decreased through periodical stereoscopic imaging intra-fractionally and phase-controlled rescanning (PCR) combined with gating for fast scanning particle beams [40]. Excellent dose conformity and enhanced OAR sparing have been demonstrated for PCR combined with gating [41]. However, rescanning during gating may not be necessary in all cases, especially when fractionated treatments are delivered, as suggested by the findings discussed above [34, 35, 38]. To account for irregular breathing patterns throughout the entire course of fractionated particle therapy, amplitude-based gating based on tumor location observed in 4D CT datasets has been adopted and routinely used clinically [42]. Although active tumor tracking may lead to the smallest high dose volume, this approach places a high demand on scanning speed, which must allow for rapid alteration of beam energy to adjust the Bragg peak in depth in relation to tumor motion [43]. This longitudinal compensation may be achieved through the use of wedges [44, 45]. However, tumor tracking remains an area of active research in thoracic particle therapy.

Owing to the scanning beam’s sensitivity to motion, 4D dose calculation for thoracic particle therapy has been advocated. The temporal density changes due to respiratory motion, which result in range uncertainties, can be incorporated into the treatment planning process. This process leads to decreased interplay and more robust treatment plans that, in contrast to 3D planning, can avoid unexpected under-dosing [46, 47]. Robust 4D planning is preferred for thoracic particle therapy, because dose errors in 3D planning are not always dependent on motion [47, 48]. However, this approach can be labor intensive and technically demanding. Thus, 3D planning with adequate motion management is still an acceptable approach in clinical practice [48]. No clinically significant differences in most of the commonly evaluated dose parameters were observed in our dose recalculation. This finding may be due to the limiting of respiratory motion in this selected group of patients, because the primary goal of the study was to compare the dose distribution among photon VMAT, IMCIT, and IMPT with minimal influence from respiratory motion. Therefore, 4D planning may have a more significant impact when significant respiratory motion is encountered. Our 4D dose evaluation is only a limited approximation of the actual dose, owing to a lack of full consideration of the interplay between respiratory motion and the dynamic beam scanning, which could be achieved only with more robust approaches for 4D dose calculation that are not commercially available currently [29, 30, 46, 47]. As a result, the limited dose variation observed after 4D dose evaluation, especially for the OAR dose parameters, may not fully capture the interplay effect and potentially be missing larger than actually observed dose uncertainties. This issue is also the major limitation of our study. Nevertheless, this study represents a step forward from 3D planning in accounting for range uncertainties that may lead to an overall decreased interplay effect [47]. The best approaches to 4D dose calculation and how to fully account for the interplay effect are beyond the scope of the current study. Robust 4D dose calculation for scanning beam particle therapy should be further assessed, and 4D planning should be considered whenever feasible, owing to the lack of direct correlation between motion amplitude and the interplay [29, 30, 47].

Another major limitation of this study is that fixed-angle particle beams are used. Dose distribution may be further improved with gantry-based systems through augmenting the known advantages of active scanning particle therapy with the increased number of angles for beam delivery. This question remains to be investigated in the future.


In comparison with VMAT, fixed-beam IMCIT led to comparable dose conformity under limited respiratory motion. However, IMCIT had significantly better tumor target dose coverage and conformity than did IMPT. Although both IMPT and IMCIT led to significantly better thoracic OAR sparing than VMAT, IMCIT may further improve normal lung and major blood vessel sparing, as compared with IMPT.



Three-dimensional conformal radiotherapy


Anisotropic analytical algorithm


Conformity index


Carbon ion radiotherapy


Clinical target volume

D95 & D99 :

Dose covering 95% and 99% of the volume, respectively

Dmax :

Maximum dose

Dmean :

Mean dose

Dmin :

Minimum dose

Dncc :

Dose to n cc of the volume

DRx :

Prescribed dose


Full-width half-maximum


Gross tumor volume


Heterogeneity index


Heart’s mean dose


Image-guided and intensity-modulated radiotherapy


Internal GTV


Intensity-modulated carbon ion therapy


Intensity-modulated proton therapy


Intensity-modulated radiotherapy


Local effect model


Multi-field optimization


Maximal intensity projection


Mean lung dose


New conformity index


Non-small cell lung cancer


Organs at risk


Phase controlled rescanning


Progressive resolution optimizer


Particle therapy


Planning target volume


Rapid arc


Radiobiological effectiveness


Stereotactic body radiation therapy

V95 & V99 :

Volume receiving at least 95% and 99% of the prescribed dose, respectively

V95%Rx and VPTV :

The volume of tissue receiving ≥95% of the prescribed dose and the volume of the PTV, respectively


Volumetric modulated arc therapy

Vn :

% of a volume receiving n Gy


  1. Chi A, Liao Z, Nguyen NP, Xu J, Stea B, Komaki R. Systemic review of the patterns of failure following stereotactic body radiation therapy in early-stage non-small-cell lung cancer: clinical implications. Radiother Oncol. 2010;94:1–11.

    Article  PubMed  Google Scholar 

  2. Louie AV, Rodrigues G, Hannouf M, Zaric GS, Palma DA, Cao JQ, et al. Stereotactic body radiotherapy versus surgery for medically operable stage I non-small-cell lung cancer: a Markov model-based decision analysis. Int J Radiat Oncol Biol Phys. 2011;81:964–73.

    Article  PubMed  Google Scholar 

  3. Shirvani SM, Jiang J, Chang JY, Welsh JW, Gomez DR, Swisher S, et al. Comparative effectiveness of 5 treatment strategies for early-stage non-small cell lung cancer in the elderly. Int J Radiat Oncol Biol Phys. 2012;84:1060–70.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Nuyttens JJ, van der Voort van Zyp NC, Verhoef C, Maat A, van Klaveren RJ, van der Holt B, et al. Stereotactic body radiation therapy for oligometastases to the lung: a phase 2 study. Int J Radiat Oncol Biol Phys. 2015;91:337–43.

    Article  PubMed  Google Scholar 

  5. Curran WJ Jr, Paulus R, Langer CJ, Komaki R, Lee JS, Hauser S, et al. Sequential vs. concurrent chemoradiation for stage III non-small cell lung cancer: randomized phase III trial RTOG 9410. J Natl Cancer Inst. 2011;103:1452–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Bradley JD, Paulus R, Komaki R, Masters G, Blumenschein G, Schild S, et al. Standard-dose versus high-dose conformal radiotherapy with concurrent and consolidation carboplatin plus paclitaxel with or without cetuximab for patients with stage IIIA or IIIB non-small-cell lung cancer (RTOG 0617): a randomised, two-by-two factorial phase 3 study. Lancet Oncol. 2015;16:187–99.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Cannon DM, Mehta MP, Adkison JB, Khuntia D, Traynor AM, Tomé WA, et al. Dose-limiting toxicity after hypofractionated dose-escalated radiotherapy in non-small-cell lung cancer. J Clin Oncol. 2013;31:4343–8.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Westover KD, Loo BW, Gerber DE, Iyengar P, Choy H, Diehn M, et al. Precision hypofractionated radiation therapy in poor performing patients with non-small cell lung cancer: phase I dose escalation trial. Int J Radiat Oncol Biol Phys. 2015;93:72–81.

    Article  PubMed  Google Scholar 

  9. Din OS, Harden SV, Hudson E, Mohammed N, Pemberton LS, Lester JF, et al. Accelerated hypo-fractionated radiotherapy for non small cell lung cancer: results from 4 UK centres. Radiother Oncol. 2013;109:8–12.

    Article  PubMed  Google Scholar 

  10. Verbakel WF, Senan S, Cuijpers JP, Slotman BJ, Lagerwaard FJ. Rapid delivery of stereotactic radiotherapy for peripheral lung tumors using volumetric intensity-modulated arcs. Radiother Oncol. 2009;93:122–4.

    Article  PubMed  Google Scholar 

  11. Ong CL, Verbakel WF, Cuijpers JP, Slotman BJ, Lagerwaard FJ, Senan S. Stereotactic radiotherapy for peripheral lung tumors: a comparison of volumetric modulated arc therapy with 3 other delivery techniques. Radiother Oncol. 2010;97:437–42.

    Article  PubMed  Google Scholar 

  12. Bongers EM, Botticella A, Palma DA, Haasbeek CJ, Warner A, Verbakel WF, et al. Predictive parameters of symptomatic radiation pneumonitis following stereotactic or hypofractionated radiotherapy delivered using volumetric modulated arcs. Radiother Oncol. 2013;109:95–9.

    Article  PubMed  Google Scholar 

  13. Liao Z, Lin SH, Cox JD. Status of particle therapy for lung cancer. Acta Oncol. 2011;50:745–56.

    Article  PubMed  Google Scholar 

  14. Kamada T, Tsujii H, Blakely EA, Debus J, De Neve W, Durante M, et al. Carbon ion radiotherapy in Japan: an assessment of 20 years of clinical experience. Lancet Oncol. 2015;16:e93–100.

    Article  PubMed  Google Scholar 

  15. Chang JY, Zhang X, Wang X, Kang Y, Riley B, Bilton S, et al. Significant reduction of normal tissue dose by proton radiotherapy compared with three-dimensional conformal or intensity-modulated radiation therapy in stage I or stage III non-small-cell lung cancer. Int J Radiat Oncol Biol Phys. 2006;65:1087–96.

    Article  PubMed  Google Scholar 

  16. Georg D, Hillbrand M, Stock M, Dieckmann K, Pötter R. Can protons improve SBRT for lung lesions? Dosimetric considerations. Radiother Oncol. 2008;88:368–75.

    Article  CAS  PubMed  Google Scholar 

  17. Macdonald OK, Kruse JJ, Miller JM, Garces YI, Brown PD, Miller RC, et al. Proton beam radiotherapy versus three-dimensional conformal stereotactic body radiotherapy in primary peripheral, early-stage non-small-cell lung carcinoma: a comparative dosimetric analysis. Int J Radiat Oncol Biol Phys. 2009;75:950–8.

    Article  PubMed  Google Scholar 

  18. Hoppe BS, Huh S, Flampouri S, Nichols RC, Oliver KR, Morris CG, et al. Double-scattered proton-based stereotactic body radiotherapy for stage I lung cancer: a dosimetric comparison with photon-based stereotactic body radiotherapy. Radiother Oncol. 2010;97:425–30.

    Article  PubMed  Google Scholar 

  19. Register SP, Zhang X, Mohan R, Chang JY. Proton stereotactic body radiation therapy for clinically challenging cases of centrally and superiorly located stage I non-small-cell lung cancer. Int J Radiat Oncol Biol Phys. 2011;80:1015–22.

    Article  PubMed  Google Scholar 

  20. Seco J, Gu G, Marcelos T, Kooy H, Willers H. Proton arc reduces range uncertainty effects and improves conformity compared with photon volumetric modulated arc therapy in stereotactic body radiation therapy for non-small cell lung cancer. Int J Radiat Oncol Biol Phys. 2013;87:188–94.

    Article  PubMed  Google Scholar 

  21. Ebara T, Shimada H, Kawamura H, Shirai K, Saito J, Kawashima M, et al. Dosimetric analysis between carbon ion radiotherapy and stereotactic body radiotherapy in stage I lung cancer. Anticancer Res. 2014;34:5099–104.

    PubMed  Google Scholar 

  22. Kubo N, Saitoh J, Shimada H, Shirai K, Kawamura H, Ohno T, et al. Dosimetric comparison of carbon ion and x-ray radiotherapy for stage IIIA non-small cell lung cancer. J Radiat Res. 2016;57:1–7.

    Article  Google Scholar 

  23. Gomez DR, Gillin M, Liao Z, Wei C, Lin SH, Swanick C, et al. Phase I study of dose escalation in hypofractionated proton beam therapy for non-small cell lung cancer. Int J Radiat Oncol Biol Phys. 2013;86:665–70.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Takahashi W, Nakajima M, Yamamoto N, Yamashita H, Nakagawa K, Miyamoto T, et al. A prospective nonrandomized phase I/II study of carbon ion radiotherapy in a favorable subset of locally advanced non-small cell lung cancer (NSCLC). Cancer. 2015;121:1321–7.

    Article  CAS  PubMed  Google Scholar 

  25. Scholz M, Kraft G. Track structure and the calculation of biological effects of heavy charged particles. Adv Space Res. 1996;18:5–14.

    Article  CAS  PubMed  Google Scholar 

  26. Elsässer T, Krämer M, Scholz M. Accuracy of the local effect model for the prediction of biologic effects of carbon ion beams in vitro and in vivo. Int J Radiat Oncol Biol Phys. 2008;71:866–72.

    Article  PubMed  Google Scholar 

  27. Kang Y, Zhang X, Chang JY, Wang H, Wei X, Liao Z, et al. 4D proton treatment planning strategy for mobile lung tumors. Int J Radiat Oncol Biol Phys. 2007;67:906–14.

    Article  PubMed  Google Scholar 

  28. Ohtakara K, Hayashi S, Hoshi H. The relation between various conformity indices and the influence of the target coverage difference in prescription isodose surface on these values in intracranial stereotactic radiosurgery. Br J Rad. 2012;85:e223–8.

    Article  CAS  Google Scholar 

  29. Richter D, Saito N, Chaudhri N, Härtig M, Ellerbrock M, Jäkel O, et al. Four-dimensional patient dose reconstruction for scanned ion beam therapy of moving liver tumors. Int J Radiat Oncol Biol Phys. 2014;89:175–81.

    Article  PubMed  Google Scholar 

  30. Bernatowicz K, Peroni M, Perrin R, Weber DC, Lomax A. Four-dimensional dose reconstruction for scanned proton therapy using liver 4DCT-MRI. Int J Radiat Oncol Biol Phys. 2016;95:216–23.

    Article  PubMed  Google Scholar 

  31. Yamashita H, Nakagawa K, Nakamura N, Koyanagi H, Tago M, Igaki H, et al. Exceptionally high incidence of symptomatic grade 2-5 radiation pneumonitis after stereotactic radiation therapy for lung tumors. Radiat Oncol. 2007;2:21.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Bert C, Durante M. Motion in radiotherapy: particle therapy. Phys Med Biol. 2011;56:R113–44.

    Article  CAS  PubMed  Google Scholar 

  33. Bert C, Grӧzinger SO, Rietzel E. Quantification of interplay effects of scanned particle beams and moving targets. Phys Med Biol. 2008;53:2253–65.

    Article  PubMed  Google Scholar 

  34. Grassberger C, Dowdell S, Lomax A, Sharp G, Shackleford J, Choi N, et al. Motion interplay as a function of patient parameters and spot size in spot scanning proton therapy for lung cancer. Int J Radiat Oncol Biol Phys. 2013;86:380–6.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Li Y, Kardar L, Li X, Li H, Cao W, Chang JY, et al. On the interplay effects with proton scanning beams in stage III lung cancer. Med Phys. 2014;41:021721.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Rietzel E. Respiratory motion management in particle therapy. Med Phys. 2010;37:449–60.

    Article  PubMed  Google Scholar 

  37. Kubiak T. Particle therapy of moving targets-the strategies for tumour motion monitoring and moving targets irradiation. Br J Radiol. 2016;89:20150275.

    Article  PubMed  Google Scholar 

  38. Knopf A-C, Hong TS, Lomax A. Scanned proton radiotherapy for mobile targets-the effectiveness of re-scanning in the context of different treatment planning approaches and for different motion characteristics. Phys Med Biol. 2011;56:7257–71.

    Article  PubMed  Google Scholar 

  39. Kanehira T, Matsuura T, Takao S, Matsuzaki Y, Fujii Y, Fujii T, et al. Impact of real-time image gating on spot scanning proton therapy for lung tumors: a simulation study. Int J Radiat Oncol Biol Phys. 2017;97:173–81.

    Article  CAS  PubMed  Google Scholar 

  40. Furukawa T, Inaniwa T, Sato S, Tomitani T, Minohara S, Noda K, et al. Design study of a raster scanning system for moving target irradiation in heavy-ion radiotherapy. Med Phys. 2007;34:1085–97.

    Article  CAS  PubMed  Google Scholar 

  41. Takahashi W, Mori S, Nakajima M, Yamamoto N, Inaniwa T, Furukawa T, et al. Carbon-ion scanning lung treatment planning with respiratory-gated phase-controlled rescanning: simulation study using 4-dimensional CT data. Radiat Oncol. 2014;9:238.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Mori S, Inaniwa T, Furukawa T, Takahashi W, Nakajima M, Shirai T, et al. Amplitude-based gated phase-controlled rescanning in carbon-ion scanning beam treatment planning under irregular breathing conditions using lung and liver 4DCTs. J Radiat Res. 2014;55:948–58.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Saito N, Bert C, Chaudhri N, Gemmel A, Schardt D, Durante M, et al. Speed and accuracy of a beam tracking system for treatment of moving targets with scanned ion beams. Phys Med Biol. 2009;54:4849–62.

    Article  CAS  PubMed  Google Scholar 

  44. Grӧzinger SO, Bert C, Haberer T, Kraft G, Rietzel E. Motion compensation with a scanned ion beam: a technical feasibility study. Radiat Oncol. 2008;3:34.

    Article  Google Scholar 

  45. Bert C, Gemmel A, Saito N, Chaudhri N, Schardt D, Durante M, et al. Dosimetric precision of an ion beam tracking system. Radiat Oncol. 2010;5:61.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Graeff C, Lüchtenborg R, Eley JG, Durante M, Bert C. A 4D-optimization concept for scanned ion beam therapy. Radiother Oncol. 2013;109:419–24.

    Article  PubMed  Google Scholar 

  47. Liu W, Schild SE, Chang JY, Liao Z, Chang YH, Wen Z, et al. Exploratory study of 4D versus 3D robust optimization in intensity modulated proton therapy for lung cancer. Int J Radiat Oncol Biol Phys. 2016;95:523–33.

    Article  PubMed  Google Scholar 

  48. Li H, Liu W, Park P, Matney J, Liao Z, Chang J, et al. Evaluation of the systematic error in using 3D dose calculation in scanning beam proton therapy for lung cancer. J Appl Clin Med Phys. 2014;15:4810.

    PubMed  PubMed Central  Google Scholar 

Download references


None for this manuscript.



Availability of data and materials

All data reported.

Author information

Authors and Affiliations



AC, LL, and WH conceived the study. AC and LL planned the study. SW and HY analyzed the data. AC drafted the manuscript, and all authors read and approved the final manuscript.

Corresponding authors

Correspondence to Alexander Chi or Wen-Chien Hsi.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

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

Additional information

Accepted for poster presentation during the Multidisciplinary Thoracic Cancers Symposium in March 2017

Additional files

Additional file 1:

Measurement of dose heterogeneity for particle therapy. (DOC 27 kb)

Additional file 2: Table S1.

Normal lung dose distribution. Table S2. Dose distribution for the heart and the esophagus. Table S3. Recalculated 4D (doses accumulated from 4D CT phases within the gating window) vs. 3D dose comparison for selected dose parameters for intensity modulated a) carbon ion therapy, and b) proton therapy with the Wilcoxon rank test. (DOC 114 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chi, A., Lin, LC., Wen, S. et al. Comparison of photon volumetric modulated arc therapy, intensity-modulated proton therapy, and intensity-modulated carbon ion therapy for delivery of hypo-fractionated thoracic radiotherapy. Radiat Oncol 12, 132 (2017).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: