Gate/Geant4-based Monte Carlo simulation for calculation of dose distribution of 400 MeV/u carbon ion beam and fragments in water

Gate/Geant4-based Monte Carlo simulation for calculation of dose distribution of 400 MeV/u carbon ion beam and fragments in water 1 2 first-author 1 1 corresp Corresponding author.E-mail address: zhaosj@zzu.edu.cn Corresponding author.E-mail address: zhaosj@zzu.edu.cn zhaosj@zzu.edu.cn Physical Science and Technology College, Zhengzhou University, Zhengzhou 450001, China Physical Science and Technology College, Zhengzhou University, Zhengzhou 450001, China College of Science, Henan University of Technology, Zhengzhou 450001, China College of Science, Henan University of Technology, Zhengzhou 450001, China The applications of carbon ion beam in tumor therapy have attracted more attention in recent years. Monte Carlo simulation is an important approach to obtain accurate radiotherapy parameters. In this work, a 400 MeV/u carbon ion beam incident on water phantom was simulated with Gate/Geant4 tools. In methods, the authors set up a carbon ion beam source according to the experiment parameters of Haettner E, defined the geometries and materials, set up the physics processes, and designed the means of information collection. In results, the authors obtained the longitudinal dose distribution, the lateral dose distribution, and the relative uncertainty of dose. The dose contributions of all kinds of fragments were calculated detailedly, and compared with the Francis Z’s results. This work is helpful for people’s understanding of the dose distributions produced by carbon ion beam and fragments in water. The simulation method is also significative for radiotherapy treatment planning of carbon ion beam, and it is easy to extend. For obtaining a special result, we may change the particle energy, particle type, target material, target geometry, physics process, detector, etc. Gate carbon ion beam fragments dose distribution Monte Carlo simulation 1 INTRODUCTION 1 The applications of proton and heavy ion in tumor therapy have attracted more attention in developed countries. Worldwide, proton and carbon ion treatment devices have been undergoing fast development in recent years. The number of patients treated with proton and carbon ion radiotherapy has been gradually increasing. In the clinic, the accuracy of radiotherapy treatment planning plays a very important role in radiotherapy quality [ 1 1 , 2 2 ] . At present, Monte Carlo simulation is an important approach to obtain accurate radiotherapy parameters. Carbon ion has better biological effectiveness than proton [ 3 3 ] . Water is considered as a kind of matter which is very approximate to body tissue, so we proposed a simulation that ion beam was incident on water phantom. Doing this work can help people to understand better the dose distribution of carbon ion beam and fragments in water. In this work, the Geant4 and Gate were chosen as simulation tools. Geant4 [ 4 4 , 5 5 ] is a Monte Carlo simulation tool to simulate energetic particles through materials. It can be used in high energy physics, nuclear physics [ 5 5 ] , accelerator physics, medical physics [ 7 7 - 10 10 7-10 ] , and space science [ 11 11 ] . The other tool, Gate [ 12 12 - 16 16 12-16 ] , an opensource software, is a toolbox based on Geant4. It is mainly used for numerical simulations in medical imaging and radiotherapy. Some articles [ 17 17 - 19 19 17-19 ] have showed the applications of Gate in radiation therapy. We adopted the Gate6.2 and Geant4.9.5 in this work. 2 SIMULATION SCHEME 1 2.1 Carbon ion beam settings 2 In this work, the settings of carbon ion beam were consistent with the source parameters of Haettner E [ 20 20 ] experiment. The energy of carbon ion was set to monoenergic 400 MeV/u. The beam was set to Gaussian distribution on the vertical and horizontal direction, and <math id="M1"><mrow><msub><mi>δ</mi><mi>x</mi></msub><mo>=</mo><msub><mi>δ</mi><mi>y</mi></msub><mo>=</mo><mn>2.12</mn><mi>m</mi><mi>m</mi></mrow></math> , of which the full width at half maximum(FWHM) was <math id="M2"><mrow><mn>5</mn><mi>m</mi><mi>m</mi></mrow></math> . The angle scattering parameter of beam was <math id="M3"><mrow><msub><mi>δ</mi><mrow><mi>a</mi><mi>n</mi><mi>g</mi></mrow></msub><mo>=</mo><mn>1</mn><mi>m</mi><mi>r</mi><mi>a</mi><mi>d</mi></mrow></math> . The number of carbon ions used for this simulation was <math id="M4"><mrow><mn>1</mn><mo>×</mo><msup><mrow><mn>10</mn></mrow><mn>5</mn></msup></mrow></math> . 2.2 Water phantom settings 2 As shown in Fig. 1 Fig. 1 , the water phantom of this simulation was a cuboid, of which the end surface was 20 cm×20 cm square, the length was 40 cm. The material parameter of water phantom was set to G4_WATER. Carbon ion beam was incident vertically on the center of left end surface. Carbon ion beam incident direction was marked as incident center line. 2.3 Detector geometry 2 In this simulation, three kinds of detector geometries, thin slice A, thin slice B, and thin strip C, were used for collecting information in water phantom, of which the geometric sizes were 200mm×200mm×1mm, 0.5mm×200mm×400mm and 0.5mm×200mm×1mm, respectively. Specially, the thin strip C was in a slice A. A thin slice A can be divided into 400 thin strip C. The thin slice A was used for collecting the longitudinal dose information of water phantom. The thin slice B was used for collecting the lateral dose information of water phantom. The thin strip C was used for collecting the lateral dose information of thin slice A in different depth. The material parameters of three geometries were also set to G4_WATER. As shown in Fig. 1 Fig. 1 , the distance between the center of thin slice A and the incident point was set to <math id="M5"><mi>d</mi></math> . 2.4 Gate method 2 In Gate, all the things of simulation must be written in a script file with mac as suffix. The script file included the several aspects: source and particle management, defining geometry and material, setting up the physics, actor management. The source, particle, geometry and material have been described ahead A. The physics processes of this simulation were set to default settings of particle radiotherapy in Gate. The settings were contained in two files. The first file, physicslist_EM_std.mac, included the standard Electromagnetic processes for leptons. The second, physicslist_HAD_std.mac, included the Electromagnetic processes for hadrons and Hadronic processes for hadrons. Actors are the tools allowed to interact with the simulation, which are used for collecting information during simulating. In this work, DoseActors was used to record the particle dose distribution; PhaseSpaceActors was used to record the particle type, energy, and location. In addition, SecondaryProductionActors was used to record the secondary particles, and particleFilter was used to filter particles. 3 RESULTS AND DISCUSSION 1 3.1 Total dose of carbon ion beam 2 3.1.1 Longitudinal total dose of carbon ion beam 3 The incident carbon ions interacted with water, and produced a variety of fragments. These fragments included ionized electrons, photons, neutrons, and H, He, Li, Be, B, C, N, O, F, Ne elements (including their isotopes). Carbon ions and all fragments transfered in water and produced dose deposits, which formed the total dose of carbon ion beam. In this simulation, the Bragg peak position depth of 400 MeV/u carbon ion beam in water was 274.5 mm (corresponding to the depth value of the position of peak height 80%). This is consistent with the Haettner E’s experimental result <math id="M6"><mrow><mn>274.7</mn><mo>±</mo><mn>1.0</mn></mrow></math> mm [ 20 20 ] . The dose distributions are shown in Fig. 2(a) Fig. 2(a) . In Fig. 2(a) Fig. 2(a) , the solid line is total dose distribution of carbon ion beam, the dashed line is dose distribution of primary 12 C ion, and the dotted line is dose distribution of all fragments. The primary 12 C ion dose has a minimum before the Bragg peak. Fragments dose reaches its maximum at 252mm depth before the Bragg peak. After the Bragg peak, the total dose line and fragments dose line are almost entirely overlapping. This is because the primary 12 C ion contributes almost nothing after the Bragg peak, the tail of total dose is primarily due to the contribution of fragments. Fig. 2(b) Fig. 2(b) is the longitudinal distribution of relative uncertainty, where the solid line indicates total dose, the dashed line indicates primary 12 C ion dose, and the dotted line indicates fragments dose. The relative uncertainty line of fragments dose has a minimum before the Bragg peak. After the Bragg peak, the relative uncertainty lines of total dose and fragments dose are almost overlapping. Fig. 2(c) Fig. 2(c) shows the whole distribution of relative uncertainty of primary 12 C ion dose. 3.1.2 Lateral total dose of carbon ion beam 3 The lateral dose distributions are shown in Fig. 3(a) Fig. 3(a) . In Fig. 3(a) Fig. 3(a) , the transverse coordinate 0 mm is the transverse center of water phantom. The solid line, the dashed line, and the dotted line represent total dose, primary 12 C ion, and all fragments respectively. It can be seen that the lateral dose distribution curve of primary 12 C ion, is narrower than the other two curves, which indicates the scattering of primary 12 C ion is more inconspicuous than fragments during transmission. The lateral total dose and lateral fragments dose are almost overlapping when the absolute value of abscissa is greater than 5mm. The dashed line and dotted line intersect at a point where the lateral coordinate is approximately equal to ±1.6 mm. In the intersection point, the lateral dose of primary 12 C ion and the lateral dose of fragments are equal. These dose distributions are consistent with spatial fragments distributions published by Matsufuji N [ 21 21 ] , Kusano Y [ 22 22 ] , and Braunn B [ 23 23 ] . Fig. 3(b) Fig. 3(b) is the lateral distribution of relative uncertainty. The solid line, the dashed line, and the dotted line represent total dose, primary 12 C ion and all fragments respectively. It can be seen from the graph that the solid line and dotted line are almost overlapping when the absolute value of abscissa is greater than 5mm, and the value of relative uncertainty is approximately equal to 0.01. 3.1.3 Lateral total dose in water slices 3 Four detector slices, such as slice A of Fig. 1 Fig. 1 , were arranged at the longitudinal direction 25 mm, 150 mm, 274 mm, and 300 mm depth of water phantom, respectively. These were used to collect the lateral dose information of carbon ion beam at the four different depths. The simulation results are shown in Fig. 4 Fig. 4 . It can be seen from Fig. 4(a) Fig. 4(a) and Fig. 4(b) Fig. 4(b) that the lateral dose line at 25 mm depth is narrow and high, and the relative dose at the center (0 mm abscissa position) is about 3.03. That is to say that the dose of incident center line at 25 mm depth is 3 times of lateral total dose maximum of carbon ion beam. In addition, the lateral dose line at 25mm depth is the steepest descent when it extends from center to both sides. The other three lines of 150 mm, 274 mm, and 300 mm depth are rising in turn when they extend from center to both sides. This is just because of particle scattering during transmission. More detector slices were arranged to collect detailed information of dose distribution at different depths. The dose distribution of incident center line is shown in Fig. 5 Fig. 5 . The relative doses of incident center line at 5mm, 25mm, 100mm, 150mm, 225mm, 274mm, and 300mm depth are 3.03, 2.90, 1.65, 0.92, 0.52, 1.22, and 0.03 respectively. 3.2 Dose of fragments 2 3.2.1 Longitudinal dose of fragments 3 In this simulation, a variety of fragments were considered. These fragments included electrons, photons, neutrons, and H-, He-, Li-, Be-, B-, C-, N-, O-, F-, Ne elements (including their isotopes). As shown in Fig. 6 Fig. 6 , each fragment forms its own dose distribution in water phantom. The thick solid line indicates total dose distribution of all fragments. Others indicate the dose distributions of H-, He-, Li-, Be-, B-, C isotopes (except 12 C) and electron respectively. Because the doses formed by N-, O-, F-, Ne isotopes and γ are so small, their dose curves are not shown in Fig. 6 Fig. 6 . From Table 1 Table 1 , we know the dose contribution of N-, O-, F-, Ne isotopes and γ in the range (1-400 mm) are between 0.1% to 0.00025%. These values are really very small. Francis Z data description: a : Carbon element dose, contains all the C isotopes dose. b : Contains only the 1 H(proton) dose. c : Contains only 4 He(α) dose. 1 1 Fragments 4 1 Dose contribution（%） 1 1   3 1 This work 1 1 Francis 1 1   1 1 1-400mm 1 1 1-275mm 1 1 276-400mm 1 1 1 to unknown mm 1 1 12 C 1 1 58.39776 1 1 65.85672 1 1 0.83115 1 1 64 a 1 1 H 1 1 8.73729 1 1 6.86052 1 1 23.22173 1 1 14 b 1 1 He 1 1 9.22231 1 1 5.11578 1 1 40.91557 1 1 13 c 1 1 Li 1 1 2.19841 1 1 1.41168 1 1 8.27021 1 1 1.7 1 1 Be 1 1 2.46668 1 1 1.78520 1 1 7.72622 1 1 1.3 1 1 B 1 1 5.81482 1 1 4.83029 1 1 13.41321 1 1 4.2 1 1 C（except 12 C） 1 1 5.08843 1 1 5.73208 1 1 0.12090 1 1 - 1 1 N 1 1 0.10605 1 1 0.10560 1 1 0.10947 1 1 - 1 1 O 1 1 0.09918 1 1 0.09320 1 1 0.14529 1 1 - 1 1 F 1 1 0.00243 1 1 0.00245 1 1 0.00227 1 1 - 1 1 Ne 1 1 0.00055 1 1 0.00058 1 1 0.00030 1 1 - 1 1 γ 1 1 0.00025 1 1 0.00026 1 1 0.00017 1 1 - 1 1 e - 1 1 7.82327 1 1 8.16501 1 1 5.18574 1 1 - 1 1 Neutron 1 1 0 1 1 0 1 1 0 1 1 - As can be seen from Fig. 6 Fig. 6 , the dose curves of H-, He-, Li-, Be-, B isotopes reach their peaks nearby the Bragg peak position. The dose curve of C isotopes (except 12 C) reaches its peak at 252 mm depth, and the dose curve of all fragments reaches its peak at 252 mm depth. Electron dose curve is also very special. It begins to fall sharply at about 250 mm depth, and becomes relatively flat again after the Bragg peak. It also can be seen from Fig. 6 Fig. 6 that the doses of H, He, B, C, and e - maintain a relatively high level before the Bragg peak, while the doses of Li, and Be are relatively low. After the Bragg peak the dose of He isotopes remains high level, while the dose of B isotopes falls steeply. 3.2.2 Lateral dose of fragments 3 The lateral dose distributions of the various fragments in water phantom are shown in Fig.7 Fig.7 . It can be seen from Fig. 7 Fig. 7 that the lateral dose curves of H and He decline relatively flat when they extend from center to both sides. This indicates that the scattering of H and He is obvious. The lateral dose curves of downmost N and O decline steeply when they extend from center to both sides. This indicates that the scattering of N and O is not obvious. Other fragments, such as Li, Be, B and C, their scattering characteristics are between H and He and N and O. It suggests that the element with greater atomic number is more difficult to scatter. In Fig. 7 Fig. 7 , the electron is not in line with the above law. The lateral dose curve is not flat when it extends from center to both sides. This is because the electron is affected mainly by the ionization interaction, and the electron dose distribution is affected by the comprehensive distribution of carbon ions and fragments. 3.3 Dose contributions of fragments 2 The dose contributions of all kinds of fragments in water phantom are shown in TABLE 1 TABLE 1 . The right column of TABLE 1 TABLE 1 shows the data obtained by Francis Z [ 24 24 ] . Francis Z gave the dose contributions of C ( 12 C and all other C isotopes), 1 H, 4 He, Li, Be, and B. The dose contributions of other elements or particles weren’t listed in his article. He only gave that the sum of dose contributions of other elements and particles was not more than 1.7%. In this work, the dose contribution of carbon element ( 12 C and all other C isotopes) in the range (1-400 mm) is 63.5%. It is similar with the Francis Z’s result 64%. Comparing the data of this study and Francis Z, we find that the differences of dose contributions of Li, Be and B are unobvious. While the differences of dose contributions of H, He and e - are obvious. In addition, we can see that the dose contributions of N, O, F, Ne, and γ are very tiny. Especially, the neutron dose contribution is zero. This is because in Geant4 neutral particles almost do not deposit any dose directly but mostly indirectly. In TABLE 1 TABLE 1 , we also list the dose contributions of various fragments in the range (1-275mm) and range (276-400mm). Comparing the data of these two ranges, we can see that some change hugely, while others change slightly. The dose contributions of some fragments of this work are consistent with the Francis Z’s [ 24 24 ] , while others have the differences. In our opinions, the reasons which lead to the differences include the following four aspects: simulation tools, physics processes, calculated particles, and phantom size. In simulation tools, this work used the Gate6.2/Geant4.9.5.1, while Francis Z used the Geant4.9.5. In physics processes, this work used the default physics settings of Gate. The first settings file, physicslist_EM_std.mac, included the standard electromagnetic processes for leptons. The second file, physicslist_HAD_std.mac, included the electromagnetic processes for hadrons and Hadronic processes for hadrons. While in Francis Z’s work, ions interactions are taken into account using the inelastic processes, e.g. ionizations, and the fragmentation process is simulated according to the quantum molecular dynamics (QMD) model. In calculated particles, this work included γ, e - , neutron, and H-, He-, Li-, Be-, B-, C-, N-, O-, F-, Ne isotopes. While Francis Z calculated the H-, He-, Li-, Be-, B-, C isotopes, and whether the dose contribution of e - was considered in his study can’t be judged. In phantom size, the dimension of water phantom in this simulation was 20cm×20cm×40cm. While Francis Z. didn’t describe the special size of water phantom in his calculation; he only said it was an infinite water phantom. 4 SUMMARY 1 As mentioned above, a 400 MeV/u carbon ion beam incident on water phantom was simulated with Gate/Geant4 tools. The carbon ion beam source was set up according to the experiment parameters of Haettner E. The geometries, materials, physics processes and detectors were defined or designed respectively. In this work, the authors obtained the longitudinal total dose distribution, the lateral total dose distribution, and the relative uncertainty of total dose. In a similar way, the longitudinal dose distribution- and the lateral dose distribution of fragments were shown. The dose contributions of all kinds of fragments were calculated detailedly, and compared with the Francis Z’s results. This work is helpful for people’s understanding of the dose distributions produced by carbon ion beam and fragments in water. The simulation method is also significative for radiotherapy treatment planning of carbon ion beam, and it is easy to extend. For obtaining a special result, we may change the particle energy, particle type, target material, target geometry, physics process, detector, etc. References Bache S T , Juang T , Belley M D . Investigating the accuracy of microstereotactic-body-radiotherapy utilizing anatomically accurate 3D printed rodent-morphic dosimeters . Med Phys , 2015 , 42 : 846 - 855 . DOI: 10.1118/1.4905489 http://doi.org/10.1118/1.4905489 Investigating the accuracy of microstereotactic-body-radiotherapy utilizing anatomically accurate 3D printed rodent-morphic dosimeters WANG Jie , PEI Xi , CAO Rui-Fen , et al . A multiphase direct aperture optimization for inverse planning in radiotherapy . Nuclear Science and Techniques , 2015 , 26 : 010502 . DOI: 10.13538/j.1001-8042/nst.26.010502 http://doi.org/10.13538/j.1001-8042/nst.26.010502 A multiphase direct aperture optimization for inverse planning in radiotherapy Frese M C , Yu V K , Stewart R D , et al . A Mechanism-Based Approach to Predict the Relative Biological Effectiveness of Protons and Carbon Ions in Radiation Therapy . Int J Radiation Oncol Biol Phys , 2012 , 83 : 442 - 450 . DOI: 10.1016/j.ijrobp.2011.06.1983 http://doi.org/10.1016/j.ijrobp.2011.06.1983 A Mechanism-Based Approach to Predict the Relative Biological Effectiveness of Protons and Carbon Ions in Radiation Therapy Agostinelli S , Allison J , Amako K , et al . GEANT4—a simulation toolkit . Nucl Instrum Methods Phys Res Sect A , 2003 , 506 : 250 - 303 . DOI: 10.1016/S0168-9002(03)01368-8 http://doi.org/10.1016/S0168-9002(03)01368-8 GEANT4—a simulation toolkit Zahra N , Frisson T , Grevillot L , et al . Influence of Geant4 parameters on dose distribution and computation time for carbon ion therapy simulation . Phys Medica , 2010 , 26 : 202 - 208 . DOI: 10.1016/j.ejmp.2009.12.001 http://doi.org/10.1016/j.ejmp.2009.12.001 Influence of Geant4 parameters on dose distribution and computation time for carbon ion therapy simulation QIN Xue , ZHOU Rong , HAN Ji-Feng , et al . GEANT4 simulation of the characteristic gamma-ray spectrum of TNT under soil induced by DT neutrons . Nuclear Science and Techniques , 2015 , 26 : 010501 . DOI: 10.13538/j.1001-8042/nst.26.010501 http://doi.org/10.13538/j.1001-8042/nst.26.010501 GEANT4 simulation of the characteristic gamma-ray spectrum of TNT under soil induced by DT neutrons Böhlen T T , Cerutti F , Dosanjh M , et al . Benchmarking nuclear models of FLUKA and GEANT4 for carbon ion therapy . Phys Med Biol , 2010 , 55 : 5833 - 5847 . DOI: 10.1088/0031-9155/55/19/014 http://doi.org/10.1088/0031-9155/55/19/014 Benchmarking nuclear models of FLUKA and GEANT4 for carbon ion therapy Pshenichnov I , Mishustin I and Greiner W . Distributions of positron-emitting nuclei in proton and carbon-ion therapy studied with GEANT4 . Phys Med Biol , 2006 , 51 : 6099 - 6112 . DOI: 10.1088/0031-9155/51/23/011 http://doi.org/10.1088/0031-9155/51/23/011 Distributions of positron-emitting nuclei in proton and carbon-ion therapy studied with GEANT4 Napoli M D , Romano F , D'Urso D , et al . Nuclear reaction measurements on tissue-equivalent materials and GEANT4 Monte Carlo simulations for hadrontherapy . Phys Med Biol , 2014 , 59 : 7643 - 7652 . DOI: 10.1088/0031-9155/59/24/7643 http://doi.org/10.1088/0031-9155/59/24/7643 Nuclear reaction measurements on tissue-equivalent materials and GEANT4 Monte Carlo simulations for hadrontherapy Pshenichnov I , Mishustin , Greiner W . Comparative study of depth-dose distributions for beams of light and heavy nuclei in tissue-like media . Nucl Instr and Meth in Phys Res B , 2008 , 266 : 1094 - 1098 . DOI: 10.1016/j.nimb.2008.02.025 http://doi.org/10.1016/j.nimb.2008.02.025 Comparative study of depth-dose distributions for beams of light and heavy nuclei in tissue-like media LIU Yun-Peng , TANG Xiao-Bin , XU Zhi-Heng , et al . Energy Deposition, Parameter Optimization, and Performance Analysis of Space Radiation Voltaic Batteries . Nuclear Science and Techniques , 2014 , 25 : S010402 . DOI: 10.13538/j.1001-8042/nst.25.S010402 http://doi.org/10.13538/j.1001-8042/nst.25.S010402 Energy Deposition, Parameter Optimization, and Performance Analysis of Space Radiation Voltaic Batteries Santin G , Strul D , Lazaro D , et al . GATE: A Geant4-based simulation platform for PET, SPECT integrating movement and time management . IEEE Trans Nucl Sci , 2003 , 50 : 1516 - 1521 . DOI: 10.1109/TNS.2003.817974 http://doi.org/10.1109/TNS.2003.817974 GATE: A Geant4-based simulation platform for PET, SPECT integrating movement and time management Sarrut D , Bardiès M , Boussion N , et al . A review of the use and potential of the GATE Monte Carlo code for radiation therapy and dosimetry applications . Med Phys , 2014 , 41 : 1 - 14 . DOI: 10.1118/1.4871617 http://doi.org/10.1118/1.4871617 A review of the use and potential of the GATE Monte Carlo code for radiation therapy and dosimetry applications Jan S , Benoit D , Becheva E , et al . GATE V6: a major enhancement of the GATE simulation platform enabling modelling of CT and radiotherapy . Phys Med Biol , 2011 , 56 : 881 - 901 . DOI: 10.1088/0031-9155/56/4/001 http://doi.org/10.1088/0031-9155/56/4/001 GATE V6: a major enhancement of the GATE simulation platform enabling modelling of CT and radiotherapy Thiam C O , Breton V , Donnarieix D , et al . Validation of a dose deposited by low-energy photons using GATE/GEANT4 . Phys Med Biol , 2008 , 53 : 3039 - 3056 . DOI: 10.1088/0031-9155/53/11/019 http://doi.org/10.1088/0031-9155/53/11/019 Validation of a dose deposited by low-energy photons using GATE/GEANT4 Grevillot L , Frisson T , Maneval D , et al . Simulation of a 6 MV Elekta Precise Linac photon beam using GATE/GEANT4 . Phys Med Biol , 2011 , 56 : 903 - 918 . DOI: 10.1088/0031-9155/56/4/002 http://doi.org/10.1088/0031-9155/56/4/002 Simulation of a 6 MV Elekta Precise Linac photon beam using GATE/GEANT4 Grevillot L , Bertrand D , Dessy F , et al . GATE as a GEANT4-based Monte Carlo platform for the evaluation of proton pencil beam scanning treatment plans . Phys Med Biol , 2012 , 57 : 4223 - 4244 . DOI: 10.1088/0031-9155/57/13/4223 http://doi.org/10.1088/0031-9155/57/13/4223 GATE as a GEANT4-based Monte Carlo platform for the evaluation of proton pencil beam scanning treatment plans Kurosu K , Takashina M , Koizumi M , et al . Optimization of GATE and PHITS Monte Carlo code parameters for uniform scanning proton beam based on simulation with FLUKA general-purpose code . Nucl Instr and Meth in Phys Res B , 2014 , 336 : 45 - 54 . DOI: 10.1016/j.nimb.2014.06.009 http://doi.org/10.1016/j.nimb.2014.06.009 Optimization of GATE and PHITS Monte Carlo code parameters for uniform scanning proton beam based on simulation with FLUKA general-purpose code Grevillot L , Frisson T , Zahra N , et al . Optimization of GEANT4 settings for Proton Pencil Beam Scanning simulations using GATE . Nucl Instr and Meth in Phys Res B , 2010 , 268 : 3295 - 3305 . DOI: 10.1016/j.nimb.2010.07.011 http://doi.org/10.1016/j.nimb.2010.07.011 Optimization of GEANT4 settings for Proton Pencil Beam Scanning simulations using GATE Haettner E , Iwase H , Krämer M , et al . Experimental study of nuclear fragmentation of 200 and 400 MeV/u 12C ions in water for applications in particle therapy . Phys Med Biol , 2013 , 58 : 8265 - 8279 . DOI: 10.1088/0031-9155/58/23/8265 http://doi.org/10.1088/0031-9155/58/23/8265 Experimental study of nuclear fragmentation of 200 and 400 MeV/u 12C ions in water for applications in particle therapy Matsufuji N , Komori M , Sasaki H , et al . Spatial fragment distribution from a therapeutic pencil-like carbon beam in water . Phys Med Biol , 2005 , 50 : 3393 - 3403 . DOI: 10.1088/0031-9155/50/14/014 http://doi.org/10.1088/0031-9155/50/14/014 Spatial fragment distribution from a therapeutic pencil-like carbon beam in water Kusano Y , Kanai T , Kase Y , et al . Dose contributions from large-angle scattered particles in therapeutic carbon beams . Med Phys , 2007 , 34 : 193 - 198 . DOI: 10.1118/1.2402328 http://doi.org/10.1118/1.2402328 Dose contributions from large-angle scattered particles in therapeutic carbon beams Braunn B , Labalme M , Ban G , et al . Nuclear reaction measurements of 95 MeV/u 12C interactions on PMMA for hadrontherapy . Nucl Instr and Meth in Phys Res B , 2011 , 269 : 2676 - 2684 . DOI: 10.1016/j.nimb.2011.08.010 http://doi.org/10.1016/j.nimb.2011.08.010 Nuclear reaction measurements of 95 MeV/u 12C interactions on PMMA for hadrontherapy Francis Z , Seif E , Incerti S , et al . Carbon ion fragmentation effects on the nanometric level behind the Bragg peak depth . Phys Med Biol , 2014 , 59 : 7691 - 7702 . DOI: 10.1088/0031-9155/59/24/7691 http://doi.org/10.1088/0031-9155/59/24/7691 Carbon ion fragmentation effects on the nanometric level behind the Bragg peak depth 10.1007/s41365-016-0097-3 This work was supported by the National Natural Science Foundation of China (No. H1809) . 2015-08-12 2016-02-02 2016-02-18 2016-07-08 2016-08-20 2021-04-11 © China Science Publishing & Media Ltd. (Science Press), Shanghai Institute of Applied Physics, the Chinese Academy of Sciences, ChineseNuclear Society and Springer Nature Singapore Pte Ltd. 2016 This is a post-peer-review, pre-copyedit version of an article published in Nuclear Science and Techniques. The final authenticated version is available online at: http://dx.doi.org/10.1007/s41365-016-0097-3 http://dx.doi.org/10.1007/s41365-016-0097-3 . 2016 Nuclear Science and Techniques 04 27