- Open Access
A step towards international prospective trials in carbon ion radiotherapy: investigation of factors influencing dose distribution in the facilities in operation based on a case of skull base chordoma
Radiation Oncology volume 14, Article number: 24 (2019)
Carbon ion radiotherapy (CIRT) has been delivered to more than 20,000 patients worldwide. International trials have been recommended in order to emphasize the actual benefits. The ULICE program (Union of Light Ion Centers in Europe) addressed the need for harmonization of CIRT practices. A comparative knowledge of the sources and magnitudes of uncertainties altering dose distribution and clinical effects during the whole CIRT procedure is required in that aim.
As part of ULICE WP2 task group, we sent a centrally reviewed questionnaire exploring candidate sources of uncertainties in dose deposition to the ten CIRT facilities in operation by February 2017. We aimed to explore native beam characterization, immobilization, anatomic data acquisition, target volumes and organs at risks delineation, treatment planning, dose delivery, quality assurance prior and during treatment. The responders had to consider the clinical case of a clival chordoma eligible for postoperative CIRT according to their clinical practice. With the results, our task group discussed ways to harmonize CIRT practices.
We received 5 surveys from facilities that have treated 77% of the patients worldwide per November 2017. We pointed out the singularity of the facilities and beam delivery systems, a divergent definition of target volumes, the multiplicity of TPS and equieffective dose calculation approximations.
Multiple uncertainties affect equieffective dose definition, deposition and calculation in CIRT. Although it is not possible to harmonize all the steps of the CIRT planning between the centers, our working group proposed counter-measures addressing the improvable limitations.
CIRT should be reported similarly to initiate international trials
Different uncertainties affect dose-effect relationship throughout CIRT planning
Equieffective dose and target volume definition and application require particular attention
We propose some suggestions to harmonize the practices
Carbon-ion radiation therapy (CIRT) is available in 10 centers worldwide . While Japanese centers have studied a broad spectrum of indications with dose escalation trials, hypofractionation and passive beam shaping (PS), the German facilities have treated smaller but homogeneous cohorts, developed fully active beam shaping (ABS), and a specific mapping of RBE-weighted dose. Per end of 2017, more than 20,000 patients have been treated worldwide with CIRT, mainly in phase I-II trials [4,5,6,7].
Due to the technical and financial investments required to expand CIRT, the decisions-makers are requesting phase III studies. These trials are particularly complex to lead in CIRT due to the rarity of indications and available resources, but also due to the different practices and environments of the dedicated centers [8, 9].
To face these challenges, it is critical that the CIRT treatments be reported similarly. Experts from IAEA and the International Commission on Radiation Units and measurements (ICRU) met in 2004 and 2006 to address these issues and published dedicated proceedings . However, there is still no official recommendation.
The ULICE project
In this framework, ULICE was a FP7 project (part of ENLIGHT) initiated in 2009 to address specific issues such as particle radiobiology, intra-fractional moving targets, adaptive treatment planning, database organization and gantry engineering . In particular, ULICE-WP2 aimed to adapt joint concepts for dose volume and outcome assessment, standard operating procedures for clinical trial design and a clinical research infrastructure.
The objective of our study was to describe and compare each stage of the treatment process in the existing CIRT centers worldwide illustrated in a case of skull base chordoma (SBC).
SBC is a rare and slow growing malignancy arising from the remnants of the notochord, optimally managed with surgery and postoperative radiation . SBC is one of the most consolidated indications of CIRT [13,14,15,16,17].
ULICE-WP2 included the following members: radiation oncologists involved in CIRT either in treating (n = 3) or in planned facilities (n = 6) or in ICRU (n = 1), physicists (n = 6), computer scientist expert in imaging (n = 1) or machine learning (n = 1), radiation biologists (n = 2) and quality manager (n = 1). Based on ten working parties, we produced a set of consensual reports released to the European Commission after a formal approval procedure involving the Steering Committee, the Technical Project Board and an External Scientific Advisory Board [18,19,20,21,22,23].
We drafted a survey addressing the potential sources of uncertainties altering clinical effects as far we may hypothesize, based on ICRU and our own reports [4, 10, 19,20,21,22,23,24,25,26,27,28] (Table 1). The questionnaire was centrally reviewed within the group and sent to the expert radiation oncologist(s) involved in SBC management in the 10 CIRT facilities in operation by February 2017.
The responders had to consider the following clinical case: a 50 y.o. man with a clival chordoma previously operated with a macroscopic residual tumor (R2) and eligible to one the following state-of-art CIRT schedules:
A) 45 Gy EQD in 3 Gy EQD/d, 5 days a week, 15 fractions on the initial PTV (PTV1) including the GTV and suspected subclinical disease – followed by a boost to a cumulative dose of 63–66 Gy EQD in 3 Gy EQD/d, 5 days a week, by 6–7 additional fractions, on a restricted PTV2 including the GTV visible on the MRI  (Fig. 1).
B) 34.2–39.6 Gy EQD in 3.8–4.4 Gy EQD/d, 4 days a week, 9 fractions on the PTV1 – followed by a boost to a cumulative dose of 60.8–70.4 Gy EQD in 3.8–4.4 Gy EQD/d, 4 days a week, by 7 additional fractions on the PTV2 .
We received 5 surveys from seven facilities representing 77% of the patients treated worldwide per November 2017 (16,704 pts. out of 21,675, Table 2A). HIT and MIT use exactly the same equipment and procedures and we distinguished HIMAC “old” and “new” facilities. We analyzed the various hypothesized sources of uncertainties.
Native beam characterization
According to Table 2B, different equipment has been set up in the successive facilities since 1994. In Japan, the two synchrotrons of the pilot Heavy Ion Medical Accelerator in Chiba (HIMAC) produce ion beams from 4He to 54Xe up to a maximum energy of 800 MeV/u. PS was first implemented with a beam-wobbling and ridge filter delivery method. The same type of facility has been reproduced in Hyogo, Gunma and Tosu, making the PS delivery method the most widely used in the world so far (Table 2). At the same time, GSI reproduced and improved the ABS system originally outlined in the pioneer facility in Berkeley - standardizing the raster scanning delivery system in European then Chinese facilities. More recently the new HIMAC facility implemented a high dose rate ABS system incorporating superconducting technology, which offers fast 3D repainting as well as a reduced gantry.
Anatomic data acquisition
Patient setup may interfere on intra/inter fraction tumor motion and the related CTV-PTV margin to apply (Table 2C).
Patients are usually treated lying in supine position - with rolling up to 20° when gantries are not available. They can also be irradiated in prone or lateral position.
In the Japanese centers, the patients are immobilized with a 3 mm-thick non-perforated thermoplastic mask over a MOLDCARE™ pillow. The mask is wrapped around the patient-pillow-cradle block and fixed with straps. Vacuum bags are used to maintain the rest of the body.
The German centers use preferentially 3-point soft continuous thermoplastic masks with individual dental splints and a set of standardized head supports.
CNAO resort to customized rigid non-perforated thermoplastic-masks, mouth bites and head-rests or moldable body pillows.
Simulation CT scan
All facilities acquire anatomic data with a CT scan with (Japan, Germany) or without (CNAO) intravenous contrast. The CT slice thickness is less than 3 mm. Hounsfield units from the simulation CT scan as well as from the mask need to be converted in relative stopping powers to accurately define water equivalent tissue thickness regarding the range of the particles.
GTV definition - imaging - fusion
All the centers perform image fusion with MRI sequences and non-deformable algorithms based on anatomical landmarks. Based on specific imaging findings, axial T1 post gadolinium and T2 fat saturated or Flair are mostly used to delineate SBC .
Procedure for CTV(s) delineation
Although the definition of the GTV appears consensual between the participating centers, there are large variations in defining the CTVs and especially the low-risk CTV (Table 2D). Subjective determinants such as clinical experience come into account.
In the treatment of cranial lesions, the major source of geometrical deviation is the relative motion of skin and bone anatomy, significantly affecting the repositioning accuracy [33,34,35]. In addition, the particle range variation is reported to be around 1 mm, representing a further source of systematic error .
Most facilities apply a 2–3 mm margin to define the PTV according to two methods: i) isotropic expansion around the related CTV; ii) portal-specific PTV.
CNAO applies the most hypofractionated dose schedule by default (B); however, whenever the PTV cannot be properly covered, schedule (A) is applied.
The European and Japanese teams apply different units and different formulae to specify and calculate equieffective dose (EQD) in relation with their radiobiological modelling.
In Mizoe et al. report, EQD (named “clinical dose”) was expressed in photon-equivalent doses (gray equivalent dose [GyE]). Planning aims to cover the PTV with at least 90% of the prescribed dose .
At GSI/HIT/MIT, EQD is rather specified to each scan spot covering the PTV. Throughout the dose distribution, at least 90% of the prescribed dose aims to cover more than 95% of the PTV unless OAR constraints limit the dose. Hot and cold spots never exceed the maximum diameter of 15 mm .
In its ULICE-report, CNAO specified dose in Gy [RBE] .
At NIRS, the calculation of “clinical dose” is based on the relationships between RBE and LET of human salivary gland tumor (HSG) cells for flattening the SOBP and on the clinical experience with fast neutrons. The RBE of carbon ions is assumed to be 3.0 at the distal region of SOBP . This method was implemented in the initial HiPLAN treatment planning system (TPS) and disseminated to the other Japanese centers . With the implementation of the active scanning delivery method at the NIRS new facility, the modified microdosimetric kinetic model (MKM) was integrated into a novel TPS .
In Europe, the RBE is determined for each single spot through the local effect model (LEM) which assumes that equal local energy deposit should lead to equal local effects in the cell nucleus, independently of the radiation quality .
Dose constraints on OAR
Dose constraints applied to OAR are summarized in Table 2E.
At HIT-MIT, equivalent EQD in 2Gy fractions are retrieved to cope with the corresponding QUANTEC recommendations for optic structures as well as brain stem surface, assuming alpha/beta values of 2Gy.
Treatment planning system
Overall, there are no major differences on the basic functionalities of the TPS (Table 2F).
The Japanese centers are equipped with XiO-N (ELEKTA, Stockholm, Sweden and Mitsubishi Electric, Tokyo, Japan) with a divergent pencil beam method algorithm based on the modified MKM EQD calculation model since 2012.
In the European centers, the inverse TPS - Syngo RT Planning (Siemens, Erlangen, Germany) is presently used. Syngo RT Planning includes a detailed radiobiological modelling of the effects of high LET radiation based on the LEM model and adopts pencil beam algorithms too.
Noteworthy, CNAO recently equipped with the Carbon Planning module on Ray Station (RaySearch Laboratories AB, Sweden) - including modules for carbon pencil beam scanning with robust biological optimization (using both the MKM and LEM models).
Two delivery systems
Since the beginning of this decade, most CIRT facilities have adopted the latter in order to limit the reduction of the beam fluence and fragmentation. Among the participating centers, HIMAC pilot facility and Gunma still use PS.
In PS, beam parameters such as MLC positions, SOBP size, range shifter thickness, range compensator shape, etc. are adapted to satisfy the tolerance doses for OARs.
With ABS, treatment planning is either possible by single field uniform dose optimization (SFUD) or by multiple field optimization (IMPT, Intensity Modulated Particle Therapy). Two to four irradiation fields are chosen at HIT-MIT.
Methods to ensure positioning reproducibility
Positioning accuracy is controlled using orthogonal x-rays or cone-beam-CTs. Set-up deviations are corrected prior to irradiation by correction with the vector of the robotic table. A 6 degrees-of-freedom correction is permitted in the most advanced facilities.
In CNAO optical tracking and in-room imaging systems are integrated .
Discussion and suggestions from ULICE WP2
Van Herk reviewed the most common uncertainties affecting dose distribution in external beam radiotherapy (EBRT) . Particle therapy raises additional issues and uncertainties in dose deposition mainly linked to the finite range as demonstrated in protontherapy .
Analyzing the whole CIRT planning process in the participating centers, we pointed out the singularity of the facilities and beam delivery systems, a divergent definition of target volumes, the multiplicity of TPS and EQD calculation approximations. Our multidisciplinary working group proposed the following counter-measures addressing the improvable limitations.
Beam production and delivery
Differences in equipment between centers may influence the components of the physical dose at the patient’s entrance.
It appears possible to accurately measure the physical dose and its composition in terms of individual particle energy, direction and LET spectrum at the exit of the delivery system. Such parameters in addition to their spatial distribution in relation with cross-section anatomy could be stored to feed more robust modeling systems in the future.
Immobilization and anatomic data acquisition
devices and the CIRT beam has to be emphasized as reported in our clinical model [45, 46]. Although they cannot easily be standardized, positioning systems and anatomic data acquisition (contrast medium use) must be perfectly described in the trial protocol as this could affect particle range calculation.
Target volume definition and delineation
The GTV and CTV are oncological volumes defined in the dedicated ICRU Reports [4, 24,25,26,27,28]. Contrarily to GTV, CTV delineation varies substantially between the centers. We could hold a joint symposium to adopt consensus delineation guidelines with standardized imaging procedures as performed in other tumor locations  - keeping in mind the natural history of the disease [48, 49].
Conversely to EBRT, adding internal and set-up margins quadratically to construct PTV  appears inadequate in CIRT since depth margins are to be incorporated to account for range uncertainties in addition to the typical lateral geometric uncertainties. When using a single beam, the different margins could theoretically combined as an anisotropic shell around the CTV . Portal-specific PTV should be theoretically employed in PS. However, when two (or more) non-parallel opposed particle beams are applied, delineation of a (unique) PTV would generate complex computation issues for the TPS . ABS raise additional issues. The dose is not prescribed to a reference point but to each voxel. Moreover, organ motion, artificial material and tissue heterogeneities may have a strong impact on dose distribution [51, 52]. Overall, specific attention and complex irradiation techniques such as gating, repainting and tracking are needed to alleviate the motion effect in CIRT although this factor is not in the foreground in SBC [53,54,55].
Dose specification and calculation
The ICRU recommends that the absorbed dose distribution (3D) be systematically reported, to allow interpretation of the effects and to reconstruct the irradiation conditions if needed [27, 56]. However, specification of these quantities alone is not usually sufficient to predict biological effects when comparing treatments in different conditions. In CIRT the particle energy spectrum varies along the ion path within the treated volume. Each voxel receives a mixed field of particles of different LET and nature due to fragmentation, so the same physical dose at a given point may correspond to different particle spectra, each of which could lead to a particular biological and clinical effect .
Despite different specifications, EQD should be reported in addition to absorbed dose [10, 57,58,59]. These two quantities are expressed in Gy. In order to avoid confusion, the symbol for EQD should be written “D (EQD)” with a space between “D” and “EQD”.
EQD is the product of absorbed dose and a transformation function (no dimension) which includes all the factors that influence the biological/medical effect(s) – including RBE  (Fig. 2). It is likely that the transformation functions are distinct between the centers and their determinants need thus to be reported.
In addition the calculation of EQD is not equivalent between the active facilities and comparing different TPS is a very difficult task in CIRT due to the many interconnected components and the different biological models .
Each individual TPS generate dose calculation uncertainties as soon the beam penetrates into the patient due to the modeling of mixed radiation field and tissue heterogeneities  (Fig. 3). Monte Carlo simulations are probably the most accurate and reliable to estimate those uncertainties . NIRS estimated that EQDs produced at NIRS and GSI are biologically identical using HSG and in vivo assays . CNAO thus recalculated their treatment plans - minimizing differences in physical dose distributions between the two treatment plans. FLUKA has been interfaced with LEM I to retrieve physical doses. Overall, the NIRS model estimates a lower RBE of carbon ions when compared to LEM. They proposed conversion tables to minimize the physical dose variations in the PTV [64,65,66].
So far, research in treatment planning of CIRT is usually based on a fixed RBE or α/β value. Since fractionation regimens differ from 3D-CRT, α/β ratios for dose-limiting toxicity and the specific tumor cell type are expected to differ substantially . Radiobiological data for both human tumors and healthy tissues are incomplete. Ideally, prospective experiments should be performed. Overlapping of PTV and OAR raises additional issues .
Quality assurance procedures
Each center has developed its own quality assurance procedures based on local equipment and the available “homemade” measuring tools. When international clinical trials will be initiated, it will be critical to have a uniform quality assurance procedure.
ICRU/IAEA Committee Report 76 - which some of the members of our committee belonged to - proposed recommendations for harmonized dosimetry guidelines and accurate beam calibration .
Thereupon, ETOILE is the first transnational prospective randomized trial comparing definitive carbon ion therapy versus photon or combined photon and protontherapy as standard treatment for unresectable or macroscopically uncompleted resected radioresistant tumors. Eligible tumors are axial chordoma (except of base of skull), adenoid cystic carcinoma of head and neck (except of trachea) and sarcomas of any site (except chondrosarcoma of the skull base), non-previously irradiated and without pre-planned surgery or chemotherapy after the clinical trial procedure. Randomization is balanced 1 for 1. Patients of the experimental arm are treated in carbon ions centers in Europe and patients of the standard arm are treated in France in their closest participating radiotherapy center. An accrual of 250 patients is needed and an absolute difference of 20% of relapse free survival at five years is awaited. The main endpoint is the progression free survival at five years (ClinicalTrials.gov Identifier: NCT02838602) We applied the aforementioned proposals when we drafted the protocol .
Multiple uncertainties affect EQD definition, deposition and calculation in CIRT. Although it is not possible to harmonize all the steps of the treatment planning between the active centers, we proposed some of the following suggestions in order to further pool the results in the frame of international trials. Reporting only EQD appears clearly insufficient. Registering and storing raw individual particle number, energy and direction prior entry into the patient and LET in each voxel (for all fragments) according to the simulation CT scan, could allow for centralized EQD reconstruction with existing or novel radiobiological models. On the other hand, it is theoretically possible to standardize the definition and delineation of structures.
Active beam shaping
Carbon ion radiotherapy
Centro Nazionale di Adroterapia Oncologica
Gesellschaft für Schwerionenforschung
Heavy Ion Medical Accelerator in Chiba
Heidelberg Ion Beam Therapy Center
International Atomic Energy Agency
International Commission on Radiation Units and Measurements
local effect model
linear energy transfer
National Institute of Radiological Sciences
Organ(s) at risk
Passive beam shaping
Particle Therapy Co-Operative Group
Relative biological effectiveness
Skull base chordoma(s)
Spread-out Bragg peak
Union of light ion centers in Europe
Barendsen GW, Koot CJ, Van Kersen GR, Bewley DK, Field SB, Parnell CJ. The effect of oxygen on impairment of the proliferative capacity of human cells in culture by ionizing radiations of different LET. Int J Radiat Biol Relat Stud Phys Chem Med. 1966;10(4):317–27.
Furusawa Y, Fukutsu K, Aoki M, Itsukaichi H, Eguchi-Kasai K, Ohara H, et al. Inactivation of aerobic and hypoxic cells from three different cell lines by accelerated (3) he-, (12) C- and (20) ne-ion beams. Radiat Res. 2000;154(5):485–96.
Particle Therapy Co-Operative Group. Particle therapy facilities in operation (last update: April 2018) [Available from: https://www.ptcog.ch/index.php/facilities-in-operation.
ICRU. Prescribing, Recording, and Reporting Brachytherapy for Cancer of the Cervix (Report 89). J ICRU. 2013;13(1–2):NP.
Combs SE, Ellerbrock M, Haberer T, Habermehl D, Hoess A, Jakel O, et al. Heidelberg Ion Therapy Center (HIT): Initial clinical experience in the first 80 patients. Acta Oncol. 2010;49(7):1132–40.
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(2):e93–e100.
Pommier P, Lievens Y, Feschet F, Borras JM, Baron MH, Shtiliyanova A, et al. Simulating demand for innovative radiotherapies: an illustrative model based on carbon ion and proton radiotherapy. Radiother Oncol. 2010;96(2):243–9.
Vanderstraeten B, Verstraete J, De Croock R, De Neve W, Lievens Y. In search of the economic sustainability of hadron therapy: the real cost of setting up and operating a hadron facility. Int J Radiat Oncol Biol Phys. 2014;89(1):152–60.
Peeters A, Grutters JP, Pijls-Johannesma M, Reimoser S, De Ruysscher D, Severens JL, et al. How costly is particle therapy? Cost analysis of external beam radiotherapy with carbon-ions, protons and photons. Radiother Oncol. 2010;95(1):45–53.
IAEA. Dose reporting in ion-beam therapy. Vienna: IAEA; 2007.
Combs SE, Djosanjh M, Potter R, Orrechia R, Haberer T, Durante M, et al. Towards clinical evidence in particle therapy: ENLIGHT, PARTNER, ULICE and beyond. J Radiat Res. 2013;54(Suppl 1):i6–12.
Di Maio S, Temkin N, Ramanathan D, Sekhar LN. Current comprehensive management of cranial base chordomas: 10-year meta-analysis of observational studies. J Neurosurg. 2011;115(6):1094–105.
Fujisawa H, Genik PC, Kitamura H, Fujimori A, Uesaka M, Kato TA. Comparison of human chordoma cell-kill for 290 MeV/n carbon ions versus 70 MeV protons in vitro. Radiat Oncol. 2013;8:91.
Schulz-Ertner D, Karger CP, Feuerhake A, Nikoghosyan A, Combs SE, Jakel O, et al. Effectiveness of carbon ion radiotherapy in the treatment of skull-base chordomas. Int J Radiat Oncol Biol Phys. 2007;68(2):449–57.
Uhl M, Mattke M, Welzel T, Roeder F, Oelmann J, Habl G, et al. Highly effective treatment of skull base chordoma with carbon ion irradiation using a raster scan technique in 155 patients: first long-term results. Cancer. 2014;120(21):3410–7.
Aydemir E, Bayrak OF, Sahin F, Atalay B, Kose GT, Ozen M, et al. Characterization of cancer stem-like cells in chordoma. J Neurosurg. 2012;116(4):810–20.
Mizoe JE. Review of carbon ion radiotherapy for skull base tumors (especially chordomas). Rep Pract Oncol Radiother. 2016;21(4):356–60.
Vogin G, Wambersie A, Pötter R, Beuve M, Combs S, Magrin G, et al. Concepts and terms for dose volume parameters in carbon-ion radiotherapy: recommendations of the ULICE taskforce. Cancer Radiother. 2018;22(8):802–9 (In press).
ULICE WP2 working group. DJRA 2.1: Harmonisation of concepts and terms for volume and dose parameters in photon, proton and carbon-ion therapy. 2011.
ULICE WP2 working group. DJRA 2.4: Joint dosimetry protocol structure enabling intercomparison between centres, including modern dosimetric and microdosimetric approaches (e.g. selection of stopping powers). 2011.
ULICE WP2 working group. DJRA 2.5: Harmonisation and recommendations for prescribing and reporting absorbed doses and dose volume histogrammes based on 3D and 4D concepts. 2012.
ULICE WP2 working group. DJRA 2.7: Document on joint outcome assessment: disease control, recurrence assessment, morbidity assessment, Quality of Life. 2012.
ULICE WP2 working group. DJRA 2.9: Integrated concept of 3D/4D absorbed dose and variations in biological effects with RBE, fractionation, overall time. 2014.
ICRU. ICRU Report 62: Prescribing, Recording and Reporting Photon Beam Therapy (Supplement to ICRU Report 50). J ICRU. 1999;os32(1):1–52.
ICRU. Prescribing, recording, and reporting Electron beam therapy. J ICRU. 2004;4(1):5–9.
ICRU. Prescribing, Recording, and Reporting Proton-Beam Therapy (Report 78). J ICRU. 2007;7(2).
ICRU. Prescribing, Recording, and Reporting Photon-Beam Intensity-Modulated Radiation Therapy (IMRT): Contents (report 83). J ICRU. 2010;10(1):NP.
Jones D. ICRU report 50—prescribing, recording and reporting photon beam therapy. Med Phys. 1994;21(6):833–4.
Nikoghosyan AV, Karapanagiotou-Schenkel I, Munter MW, Jensen AD, Combs SE, Debus J. Randomised trial of proton vs. carbon ion radiation therapy in patients with chordoma of the skull base, clinical phase III study HIT-1-Study. BMC Cancer. 2010;10:607.
Mizoe JE, Hasegawa A, Takagi R, Bessho H, Onda T, Tsujii H. Carbon ion radiotherapy for skull base chordoma. Skull Base. 2009;19(3):219–24.
Erdem E, Angtuaco EC, Van Hemert R, Park JS, Al-Mefty O. Comprehensive review of intracranial chordoma. Radiographics. 2003;23(4):995–1009.
ULICE WP10 working group. DNA 10.3: Imaging and contouring procedures for chordoma of the skull base at the Italian National Center for Hadrontherapy Oncology (CNAO). 2012.
Engelsman M, Rosenthal SJ, Michaud SL, Adams JA, Schneider RJ, Bradley SG, et al. Intra- and interfractional patient motion for a variety of immobilization devices. Med Phys. 2005;32(11):3468–74.
Karger CP, Jakel O, Debus J, Kuhn S, Hartmann GH. Three-dimensional accuracy and interfractional reproducibility of patient fixation and positioning using a stereotactic head mask system. Int J Radiat Oncol Biol Phys. 2001;49(5):1493–504.
Jensen AD, Winter M, Kuhn SP, Debus J, Nairz O, Munter MW. Robotic-based carbon ion therapy and patient positioning in 6 degrees of freedom: setup accuracy of two standard immobilization devices used in carbon ion therapy and IMRT. Radiat Oncol. 2012;7:51.
Rietzel E, Schardt D, Haberer T. Range accuracy in carbon ion treatment planning based on CT-calibration with real tissue samples. Radiat Oncol. 2007;2:14.
Matsufuji N, Kanai T, Kanematsu N, Miyamoto T, Baba M, Kamada T, et al. Specification of Carbon Ion Dose at the National Institute of Radiological Sciences (NIRS). J Radiat Res. 2007;48(Suppl A):A81–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(1):201–10.
Inaniwa T, Furukawa T, Kase Y, Matsufuji N, Toshito T, Matsumoto Y, et al. Treatment planning for a scanned carbon beam with a modified microdosimetric kinetic model. Phys Med Biol. 2010;55(22):6721–37.
Elsasser T, Kramer 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(3):866–72.
Kitagawa A, Fujita T, Muramatsu M, Biri S, Drentje AG. Review on heavy ion radiotherapy facilities and related ion sources (invited). Rev Sci Instrum. 2010;81(2):02B909.
Baroni G, Riboldi M, Spadea MF, Tagaste B, Garibaldi C, Orecchia R, et al. Integration of Enhanced Optical Tracking Techniques and Imaging in IGRT. J Radiat Res. 2007;48(Suppl A):A61–74.
van Herk M. Errors and margins in radiotherapy. Semin Radiat Oncol. 2004;14(1):52–64.
Cheng C, Zhao L, Zhao Q, Moskvin V, Buchsbaum J, Das I. SU-E-T-295: Factors Affecting Accuracy in Proton Therapy. Med Phys. 2012;39(6Part14):3771.
Fattori G, Riboldi M, Scifoni E, Kramer M, Pella A, Durante M, et al. Dosimetric effects of residual uncertainties in carbon ion treatment of head chordoma. Radiother Oncol. 2014;113(1):66–71.
ULICE WP10 working group. DNA 10.4: Immobilization Procedures in the Active Centers. 2013.
Lee AW, Ng WT, Pan JJ, Poh SS, Ahn YC, AlHussain H, et al. International guideline for the delineation of the clinical target volumes (CTV) for nasopharyngeal carcinoma. Radiother Oncol. 2018;126(1):25–36.
Fagundes MA, Hug EB, Liebsch NJ, Daly W, Efird J, Munzenrider JE. Radiation therapy for chordomas of the base of skull and cervical spine: patterns of failure and outcome after relapse. Int J Radiat Oncol Biol Phys. 1995;33(3):579–84.
Vogin G, Calugaru V, Bolle S, George B, Oldrini G, Habrand JL, et al. Investigation of ectopic recurrent skull base and cervical chordomas: the Institut Curie's proton therapy center experience. Head Neck. 2016;38(Suppl 1):E1238–46.
Cabal GA, Jakel O. Dynamic target definition: a novel approach for PTV definition in ion beam therapy. Radiother Oncol. 2013;107(2):227–33.
Ammazzalorso F, Graef S, Weber U, Wittig A, Engenhart-Cabillic R, Jelen U. Dosimetric consequences of intrafraction prostate motion in scanned ion beam radiotherapy. Radiother Oncol. 2014;112(1):100–5.
Bert C, Durante M. Motion in radiotherapy: particle therapy. Phys Med Biol. 2011;56(16):R113–44.
Lomax AJ. Intensity modulated proton therapy and its sensitivity to treatment uncertainties 2: the potential effects of inter-fraction and inter-field motions. Phys Med Biol. 2008;53(4):1043–56.
Unkelbach J, Oelfke U. Incorporating organ movements in inverse planning: assessing dose uncertainties by Bayesian inference. Phys Med Biol. 2005;50(1):121–39.
Maleike D, Unkelbach J, Oelfke U. Simulation and visualization of dose uncertainties due to interfractional organ motion. Phys Med Biol. 2006;51(9):2237–52.
ICRU. Fundamental quantities and units for ionizing radiation (report 85). J ICRU. 2011;11(1):1–31.
Wambersie A, Menzel HG, Andreo P, DeLuca PM Jr, Gahbauer R, Hendry JH, et al. Isoeffective dose: a concept for biological weighting of absorbed dose in proton and heavier-ion therapies. Radiat Prot Dosim. 2011;143(2–4):481–6.
BIPM (Bureau international des Poids et Mesures). Système international d’unités (SI), 7ème édition. France: Sèvres; 1998.
Bentzen SM, Dorr W, Gahbauer R, Howell RW, Joiner MC, Jones B, et al. Bioeffect modeling and equieffective dose concepts in radiation oncology--terminology, quantities and units. Radiother Oncol. 2012;105(2):266–8.
Doolan PJ, Alshaikhi J, Rosenberg I, Ainsley CG, Gibson A, D'Souza D, et al. A comparison of the dose distributions from three proton treatment planning systems in the planning of meningioma patients with single-field uniform dose pencil beam scanning. J Appl Clin Med Phys. 2015;16(1):4996.
Gueulette J, Wambersie A. Comparison of the Methods of Specifying Carbon Ion Doses at NIRS and GSI. J Radiat Res. 2007;48(Suppl A):A97–A102.
Mairani A, Brons S, Cerutti F, Fasso A, Ferrari A, Kramer M, et al. The FLUKA Monte Carlo code coupled with the local effect model for biological calculations in carbon ion therapy. Phys Med Biol. 2010;55(15):4273–89.
Uzawa A, Ando K, Koike S, Furusawa Y, Matsumoto Y, Takai N, et al. Comparison of biological effectiveness of carbon-ion beams in Japan and Germany. Int J Radiat Oncol Biol Phys. 2009;73(5):1545–51.
Fossati P, Molinelli S, Matsufuji N, Ciocca M, Mirandola A, Mairani A, et al. Dose prescription in carbon ion radiotherapy: a planning study to compare NIRS and LEM approaches with a clinically-oriented strategy. Phys Med Biol. 2012;57(22):7543–54.
Molinelli S, Magro G, Mairani A, Matsufuji N, Kanematsu N, Inaniwa T, et al. Dose prescription in carbon ion radiotherapy: how to compare two different RBE-weighted dose calculation systems. Radiother Oncol. 2016;120(2):307–12.
Steinstrater O, Grun R, Scholz U, Friedrich T, Durante M, Scholz M. Mapping of RBE-weighted doses between HIMAC- and LEM-based treatment planning systems for carbon ion therapy. Int J Radiat Oncol Biol Phys. 2012;84(3):854–60.
Jones B. A simpler energy transfer efficiency model to predict relative biological effect for protons and heavier ions. Front Oncol. 2015;5:184.
Steinstrater O, Scholz U, Friedrich T, Kramer M, Grun R, Durante M, et al. Integration of a model-independent interface for RBE predictions in a treatment planning system for active particle beam scanning. Phys Med Biol. 2015;60(17):6811–31.
Jones B, Underwood TS, Carabe-Fernandez A, Timlin C, Dale RG. Fast neutron relative biological effects and implications for charged particle therapy. Br J Radiol. 2011;84(Spec 1):S11–8.
Dreher C, Scholz C, Pommer M, Brons S, Prokesch H, Ecker S, et al. Optimization of carbon ion treatment plans by integrating tissue specific alpha/beta-values for patients with non-Resectable pancreatic Cancer. PLoS One. 2016;11(10):e0164473.
IAEA. TRS no. 398; absorbed dose determination in external beam radiotherapy. Vienna: IAEA; 2000.
ULICE WP2 working group. DJRA 2.10: Adaptation of the so far existing protocols (that will be provided by WP 10) within the facilities according to the needs as defined in this SOP for clinical trial design for hadron-therapy – in cooperation with WP 11. 2013.
We are grateful to Richard Pötter, Michael Beuve, Stephanie E Combs, Giulio Magrin, Ramona Mayer, Ulrike Mock, David Sarrut and Thomas Schreiner who regularly participated to ULICE WP2 meetings.
ULICE is co-funded by the European Commission within the Framework Program 7 Capacities Specific Program, under Grant Agreement Number 228436.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.