Energy loss of degrader in SC200 proton therapy facility

Energy loss of degrader in SC200 proton therapy facility 1 2 corresp first-author Corresponding author: jfeng@ipp.ac.cn Corresponding author: jfeng@ipp.ac.cn jfeng@ipp.ac.cn 1 2 1 2 1 2 1 2 1 2 1 Institute of Plasma Physics, Chinese Academy of Sciences, Hefei 230031, China Institute of Plasma Physics, Chinese Academy of Sciences, Hefei 230031, China University of Science and Technology of China, Hefei 230026, China University of Science and Technology of China, Hefei 230026, China The proton beam energy determines the range of particles and thus where the dose is deposited. According to the depth of tumors, an energy degrader is needed to modulate the proton beam energy in proton therapy facilities based on cyclotrons, because the energy of beam extracted from the cyclotron is fixed. The energy loss was simulated for the graphite degrader used in the beamline at the Superconducting Cyclotron of 200 MeV in Hefei (SC200). After adjusting the mean excitation energy of the graphite used in the degrader to 76 eV, we observed an accurate match between the simulations and measurements. We also simulated the energy spread of the degraded beam and the transmission of the degrader using theoretical formulae. The results agree well with the Monte Carlo simulation. Degrader Energy loss Mean excitation energy Energy spread Transmission 1. Introduction 1 Since 1946, when R. Wilson revealed the advantageous dose distributions of protons, it has been shown that energetic protons can be useful for cancer therapy applications. Proton depth dose distributions are characterized by a relatively low dose at shallow depths, a peak near the end of the proton range, and then a rapid falloff. The peak is called the Bragg peak, as it was discovered by William Henry Bragg in 1903 [ 1 1 ] . Proton therapy facilities can deliver a high dose of ionizing radiation to a deep-seated tumor while not exceeding the tolerance dose of the intervening normal tissue and delivering nearly no dose to normal tissues beyond the tumor [ 2 2 , 3 3 ] . The SC200 proton therapy facility (Superconducting Cyclotron of 200 MeV in Hefei, China) consists of a 200 MeV superconducting cyclotron, whose maximum beam intensity is approximately 400 nA, and a beamline that guides the beam to two treatment rooms, Gantry treatment room and Fixed beam treatment room [ 4 4 , 5 5 ] . The beam energy is modulated using a degrader, which is inserted into the beam trajectory. The degrader can reduce the beam energy from 200 MeV to a value specified by the treatment planning, which varies at SC200 in the range of 70–190 MeV. The beam intensity of the degraded proton beam at the entrances of the treatment rooms depends strongly on the beam degrading process. The minimum transmission is obtained at 70 MeV, owing to the increase in the beam emittance caused by the multiple scattering that occurs when the proton beam is degraded. The increase in the energy spread, which must be limited in order to reduce the distal falloff of the Bragg peak, affects the transmission. Some proton therapy facilities with a high primary beam energy degrade the beam energy in two steps for higher transmission. First, the beam energy is degraded to approximately 150 MeV in the beamline. Then, another degrader in the treatment room degrades the energy to the required energy. However, this is disadvantageous for treatments that require a sharp distal falloff of the Bragg peak. On the other hand, this method is very attractive for pencil beam scanning facilities. The larger distal falloff can reduce the number of scanning layers and the treatment time [ 6 6 ] . In SC200, the beam must be controlled over the full available range. Thus, there is only one degrader in SC200, and detailed information about the degrading process is required. We do not aim for a universal theory that can solve all problems regarding protons passing through materials but rather for a set of equations allowing a degrader to be designed for a proton therapy facility. Moreover, we want to predict the transmission and properties of the beam, similar to the result of Monte Carlo program. In this study, we investigated the energy loss and beam loss of the degrader using a theoretical calculation method. All of the programs were developed and executed in MATLAB. Then, we used several Monte Carlo codes to simulate the process of protons passing through the degrader. Herein, we present the experimental results to verify the rationality of the physical model. 2. Experimental Section 1 2.1 Degrader description 2 The Energy Selection System (ESS) is an important subsystem of proton therapy facilities based on cyclotrons. This subsystem is necessary in proton therapy, as the beam energy extracted from the cyclotron is fixed and the required energy of treatment plans varies according to the depth of tumors. The ESS can adjust the energy and other properties of the extracted beam in accordance with the requirement of the treatment plan [ 7 7 ] . A typical ESS often includes a degrader, collimators, slits, quadrupole magnets, bending magnets, correction magnets, vacuum equipment, and beam diagnostic equipment [ 8 8 - 10 10 8-10 ] . The beam extracted from the SC200 cyclotron is focused at the degrader, by passing it through four quadrupole magnets. A double waist, in the horizontal and vertical directions, is formed in the degrader by a special design of beam optics, for minimizing the emittance after the degrader and enhancing the transmission of the ESS. Before the degrader, there is a "control point" to monitor the parameters of the extracted beam, including a current monitor placed permanently in the beam trajectory, a profile monitor, and a beam stopper, which is inserted into the beam when necessary. After the "control point", collimator 1 is mainly used for reducing the influence of stray particles. The layout of ESS in SC200 was shown in Fig. 1 Fig. 1 . The degrader consists of two parts, which are placed at the left and right sides of the beam. Each part has several full and half wedges, and the mid-plane of the degrader is the same as beam. This design of wedges can significantly reduce the time of degrader movement, ensuring a change time of < 2 s/layer. The wedges are made of high-density graphite (ρ = 1.84 ± 0.005 g/cm 3 ) and fixed on an aluminum base platform. On the base platform, we designed a water-cooled brass connection to remove heat from the wedges, which is caused by energy loss. Two servomotors are used to drive the wedges, which adjust the overlap of the degrader. These servomotors have very high control precision, along with a grating scale to verify their correct setting. To obtain a uniform thickness of the degrader, the entrance and exit of the degrader are parallel to each other and perpendicular to the axis of beam. After the degrader, the beam passes through three collimators in sequence. The first and the third collimators are made of brass, have different cone apertures, and define the acceptable emittance for the downstream beamline. According to the requirement of the gantry and design of the downstream beamline, all of these apertures are circular for obtaining a round-beam at the "control point" before the gantry. The first collimator defines the cross-sectional size of the degraded beam. The second collimator, which is very close to the first one, is made of graphite for absorbing stray particles. The third collimator defines the maximum scattering angle. This system of collimators controls the emittance in the range of <16π mm·mrad, owing to the acceptance of the gantry. Next are six quadruple magnets and two bending magnets, forming an achromatic system. The beam waist in the horizontal direction is designed at the center of the two bending magnets. At this location, the function between the envelope curve and the momentum spread is maximized (28 mm / %ΔP), and the momentum slit is placed here. The momentum spread of the transported beam can be selected by a slit with two servomotors in the horizontal directions. The maximum momentum spread is ±1.2%. A momentum spread of ±0.6% is used for the gantry, to obtain a spot whose full width at half maximum is 4–10 mm at the isocenter. 2.2 Ionization energy-loss calculation 2 To ensure that all of the beam is degraded, the width of the degrader overlap should be greater than the diameter of the beam. This method can avoid the uncertainty of the degraded distribution, owing to the statistical fluctuation that arises when a proton passes through a slice. On the other hand, this affects the available range of the beam energy. The kinematic maximum possible electron recoil kinetic energy is given as <math display="block" id="M1"><mrow><msub><mi>T</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>=</mo><mfrac><mrow><mn>2</mn><msub><mi>m</mi><mtext>e</mtext></msub><msup><mi>c</mi><mn>2</mn></msup><msup><mi>β</mi><mn>2</mn></msup><msup><mi>γ</mi><mn>2</mn></msup></mrow><mrow><mn>1</mn><mo>+</mo><msup><mrow><mrow><mo>(</mo><mrow><msub><mi>m</mi><mtext>e</mtext></msub><mo>/</mo><mi>M</mi></mrow><mo>)</mo></mrow></mrow><mn>2</mn></msup><mo>+</mo><mn>2</mn><mrow><mo>(</mo><mrow><msub><mi>m</mi><mtext>e</mtext></msub><mo>/</mo><mi>M</mi></mrow><mo>)</mo></mrow><mo>⋅</mo><mi>γ</mi></mrow></mfrac><mo>,</mo></mrow></math> where m e is the mass of the electron, c is the speed of light, β is the ratio of the proton speed to light speed, γ = (E+E 0 )/E 0 , E 0 is the particle’s rest energy, E is the kinetic energy, and M is the mass of the projectile. According to the definition of a "thick" medium, the average ionization energy loss in thickness δx should be defined by <math display="block" id="M2"><mrow><mrow><mo>(</mo><mrow><mfrac><mrow><mtext>d</mtext><mi>E</mi></mrow><mrow><mtext>d</mtext><mi>x</mi></mrow></mfrac></mrow><mo>)</mo></mrow><mi>δ</mi><mi>x</mi><mo>></mo><mo>></mo><msub><mi>T</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>.</mo></mrow></math> The distribution of the degraded energy passed through a "thick" medium should be similar to a Gaussian distribution, and that for a "thin" medium should be similar to a Landau distribution. For a 200-MeV beam, the threshold between "thick" and "thin" is approximately 2.04 g/cm 2 . In SC200, the minimum width of the overlap is approximately 13 mm, and the thickness of graphite is approximately 12 mm, allowing a degradation of an incident 200-MeV proton beam to the range of 70–191 MeV. We use the Bethe–Bloch formula for the energy loss and range calculation [ 11 11 - 14 14 11-14 ] . The energy loss due to only ionization is given by <math display="block" id="M3"><mrow><mo>−</mo><mfrac><mrow><mtext>d</mtext><mi>E</mi></mrow><mrow><mtext>d</mtext><mi>x</mi></mrow></mfrac><mo>=</mo><mo>−</mo><mi>K</mi><msubsup><mi>Z</mi><mn>1</mn><mn>2</mn></msubsup><mi>ρ</mi><mfrac><mrow><msub><mi>Z</mi><mn>2</mn></msub></mrow><mrow><msub><mi>A</mi><mn>2</mn></msub></mrow></mfrac><mfrac><mn>1</mn><mrow><msup><mi>β</mi><mn>2</mn></msup></mrow></mfrac><mrow><mo>(</mo><mrow><mfrac><mn>1</mn><mn>2</mn></mfrac><mi>l</mi><mi>n</mi><mfrac><mrow><mn>2</mn><msub><mi>m</mi><mtext>e</mtext></msub><msup><mi>c</mi><mn>2</mn></msup><msup><mi>β</mi><mn>2</mn></msup><msup><mi>γ</mi><mn>2</mn></msup></mrow><mrow><msup><mi>I</mi><mn>2</mn></msup></mrow></mfrac><msub><mi>T</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>−</mo><msup><mi>β</mi><mn>2</mn></msup><mo>−</mo><mfrac><mi>C</mi><mi>Z</mi></mfrac><mo>−</mo><mfrac><mi>δ</mi><mn>2</mn></mfrac></mrow><mo>)</mo></mrow><mo>,</mo></mrow></math> where K is a constant equal to 0.3071; Z 1 is the charge of the incident particle (in units of e ); ρ is the density of the degrader material; Z 2 and A 2 are the atomic number and mass of the target atom, respectively; I is the mean excitation energy; δ is the density-effect correction; and C is the shell correction. If the degrader is made of a compound material, the energy loss should be calculated for each component with relevant density: <math display="block" id="M4"><mrow><mo>−</mo><mfrac><mrow><mtext>d</mtext><mi>E</mi></mrow><mrow><mi>d</mi><mi>x</mi></mrow></mfrac><mo>=</mo><mo>−</mo><mi>K</mi><msubsup><mi>Z</mi><mn>1</mn><mn>2</mn></msubsup><mfrac><mn>1</mn><mrow><msup><mi>β</mi><mn>2</mn></msup></mrow></mfrac><mstyle displaystyle="true"><munder><mo>∑</mo><mi>k</mi></munder><mrow><msub><mi>ρ</mi><mi>k</mi></msub><mfrac><mrow><msub><mi>Z</mi><mrow><mn>2</mn><mi>k</mi></mrow></msub></mrow><mrow><msub><mi>A</mi><mrow><mn>2</mn><mi>k</mi></mrow></msub></mrow></mfrac><mrow><mo>(</mo><mrow><mfrac><mn>1</mn><mn>2</mn></mfrac><mi>l</mi><mi>n</mi><mfrac><mrow><mn>2</mn><msub><mi>m</mi><mtext>e</mtext></msub><msup><mi>c</mi><mn>2</mn></msup><msup><mi>β</mi><mn>2</mn></msup><msup><mi>γ</mi><mn>2</mn></msup></mrow><mrow><msub><mi>I</mi><mi>k</mi></msub><msup><mrow/><mn>2</mn></msup></mrow></mfrac><msub><mi>T</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>−</mo><msup><mi>β</mi><mn>2</mn></msup><mo>−</mo><mfrac><mrow><msub><mi>C</mi><mi>k</mi></msub></mrow><mi>Z</mi></mfrac><mo>−</mo><mfrac><mrow><msub><mi>δ</mi><mi>k</mi></msub></mrow><mn>2</mn></mfrac></mrow><mo>)</mo></mrow></mrow></mstyle><mo>.</mo></mrow></math> The range of protons with the incident kinetic energy E 0 in materials is defined by <math display="block" id="M5"><mrow><mi>R</mi><mo stretchy="false">(</mo><msub><mi>E</mi><mn>0</mn></msub><mo stretchy="false">)</mo><mo>=</mo><mstyle displaystyle="true"><mrow><msubsup><mo>∫</mo><mn>0</mn><mrow><msub><mi>E</mi><mn>0</mn></msub></mrow></msubsup><mrow><mfrac><mrow><mtext>d</mtext><mi>E</mi></mrow><mrow><mrow><mo>|</mo><mrow><mtext>d</mtext><mi>E</mi><mo>/</mo><mtext>d</mtext><mi>x</mi></mrow><mo>|</mo></mrow></mrow></mfrac></mrow></mrow></mstyle><mo>.</mo></mrow></math> The resulting range was obtained by integrating Eq. (3) via the Runge–Kutta method. We calculated the proton range for several materials, including pure and compound materials. 2.3 Energy-spread calculation 2 The process of passing a beam through the degrader material has obvious randomness, which incurs fluctuation of the energy loss. This means that a monoenergetic beam passing through degrader material will yield an energy spread (or momentum spread). This spread and the final energy are defined by the initial energy of the incident beam and the thickness of the degrader material, which are nonlinear [ 15 15 - 17 17 15-17 ] . In combination with the stochastic process, the energy spread should be defined by <math display="block" id="M6"><mrow><msubsup><mi>σ</mi><mrow><msub><mi>E</mi><mrow><mtext>final</mtext></mrow></msub></mrow><mn>2</mn></msubsup><mo>=</mo><msup><mrow><mrow><mo>(</mo><mrow><mfrac><mrow><mtext>d</mtext><msub><mi>E</mi><mrow><mtext>final</mtext></mrow></msub><mo>/</mo><mtext>d</mtext><mi>x</mi></mrow><mrow><mtext>d</mtext><msub><mi>E</mi><mrow><mtext>initial</mtext></mrow></msub><mo>/</mo><mtext>d</mtext><mi>x</mi></mrow></mfrac></mrow><mo>)</mo></mrow></mrow><mn>2</mn></msup><mrow><mo>[</mo><mrow><msubsup><mi>σ</mi><mrow><msub><mi>E</mi><mrow><mtext>initial</mtext></mrow></msub></mrow><mn>2</mn></msubsup><mo>+</mo><msup><mrow><mrow><mo>(</mo><mrow><mfrac><mrow><mtext>d</mtext><msub><mi>E</mi><mrow><mtext>initial</mtext></mrow></msub></mrow><mrow><mtext>d</mtext><mi>x</mi></mrow></mfrac></mrow><mo>)</mo></mrow></mrow><mn>2</mn></msup><mstyle displaystyle="true"><mrow><msubsup><mo>∫</mo><mrow><msub><mi>E</mi><mrow><mtext>initial</mtext></mrow></msub></mrow><mrow><msub><mi>E</mi><mrow><mtext>final</mtext></mrow></msub></mrow></msubsup><mrow><mfrac><mrow><mi>N</mi><mrow><mo>(</mo><mi>E</mi><mo>)</mo></mrow></mrow><mrow><mi>M</mi><msup><mrow><mrow><mo>(</mo><mi>E</mi><mo>)</mo></mrow></mrow><mn>3</mn></msup></mrow></mfrac><mtext>d</mtext><mi>E</mi></mrow></mrow></mstyle></mrow><mo>]</mo></mrow><mo>,</mo></mrow></math> <math display="block" id="M7"><mrow><mi>N</mi><mrow><mo>(</mo><mi>E</mi><mo>)</mo></mrow><mo>=</mo><mi>K</mi><msub><mi>m</mi><mtext>e</mtext></msub><msup><mi>c</mi><mn>2</mn></msup><mi>ρ</mi><mfrac><mrow><msub><mi>Z</mi><mn>2</mn></msub></mrow><mrow><msub><mi>A</mi><mn>2</mn></msub></mrow></mfrac><mrow><mo>(</mo><mrow><mn>1</mn><mo>−</mo><msup><mi>β</mi><mn>2</mn></msup><mo>/</mo><mn>2</mn></mrow><mo>)</mo></mrow><mtext> </mtext><msup><mi>γ</mi><mn>2</mn></msup><mo>,</mo></mrow></math> where M ( E ) is the mean ionization energy loss, and N ( E ) is the stochastic energy spread for thin targets, which depends on the energy. According to the aforementioned theories and methods, we wrote the calculation programs using MATLAB. 2.4 Proton-nucleus collision model 2 When the proton beam passes through the degrader, the number of emerging protons undergoes exponential decay owing to the inelastic collision between the proton and nucleus. The proton–nucleus collision is a special case of the nucleus–nucleus collision. These parameters of the collision system are associated with the physics process. In general, Coulomb interaction plays a leading role at lower energies, and Pauli blocking has more influence at higher energies. For two colliding ions, the total cross section σ T equals the reaction cross section σ R plus the elastic cross section σ e . The empirical models often approximate the Bradt–Peters form. In addition, the effect of the Coulomb interaction should be considered at low energies. R. K. Tripathi [ 18 18 ] proposed the following model for calculating the reaction cross section: <math display="block" id="M8"><mrow><msub><mi>σ</mi><mtext>R</mtext></msub><mo>=</mo><mi>π</mi><mtext> </mtext><msubsup><mi>r</mi><mn>0</mn><mn>2</mn></msubsup><msup><mrow><mrow><mo>(</mo><mrow><msubsup><mi>A</mi><mn>1</mn><mrow><mn>1</mn><mo>/</mo><mn>3</mn></mrow></msubsup><mo>+</mo><msubsup><mi>A</mi><mn>2</mn><mrow><mn>1</mn><mo>/</mo><mn>3</mn></mrow></msubsup><mo>+</mo><msub><mi>δ</mi><mtext>E</mtext></msub></mrow><mo>)</mo></mrow></mrow><mn>2</mn></msup><mrow><mo>(</mo><mrow><mn>1</mn><mo>−</mo><mfrac><mi>B</mi><mrow><msub><mi>E</mi><mrow><mtext>cm</mtext></mrow></msub></mrow></mfrac></mrow><mo>)</mo></mrow><mo>,</mo></mrow></math> where r 0 is a constant not dependent on the energy, A 1 is the mass of the incident particle, E cm is the colliding center of mass energy, and B is the Coulomb interaction barrier, which depends on the energy. More importantly, the study provided a good approximation for proton–nucleus collisions. Because of the small compression effect for proton–nucleus collisions, a coefficient D , that is related to the density of the colliding system is defined as a constant to be a very good match at all energies. The proportion of the proton beam in the entire decay process is determined by <math display="block" id="M9"><mrow><mi>N</mi><mo>=</mo><msub><mi>N</mi><mn>0</mn></msub><msup><mi>e</mi><mrow><mo>−</mo><mstyle displaystyle="true"><mo>∑</mo><mrow><mi>ρ</mi><mo>/</mo><mi>λ</mi></mrow></mstyle></mrow></msup><mo>,</mo></mrow></math> where N and N 0 represent the numbers of protons in the degraded and incident beams, respectively, and λ is the mean free path. 2.5 Experiment 2 Two graphite triple wedge blocks ( Fig. 2a Fig. 2a ) were produced by the Institute of Solid State Physics. The average grain size of this isostatic graphite is only 7 μm, and the bulk density is approximately 1.84 g/cm 3 . These two blocks were mounted on an aluminum frame with two servomotors for driving them. The thickness of graphite can reach 130 mm, which is suitable for SC200. (a) Graphite blocks (b) MLIC As the transmission depends strongly on the energy, it is important to set the degraded beam energy with high accuracy. We placed the graphite blocks on the treatment bed of Gantry 2 at the Paul Scherrer Institute (PSI), and aligned them carefully in order to make the beam entrance perpendicular to the beam. The energy of the "initial beam", before these blocks, was divided into 116 points from 70 to 230 MeV by adjusting the beamline and was measured by monitors in the nozzle. The thickness of the graphite blocks was set at 15 positions (intervals of 10 mm from 6 mm to 126 mm, 0 mm, and 130 mm). We used a multi-layer ionization chamber (MLIC, Fig. 2b Fig. 2b ) to measure the beam energy after the graphite blocks. The measurement range of this MLIC was approximately 277 mm water phantom, with an error of approximately 1 mm. The beam currents before and after the graphite blocks were measured to obtain the transmission by using monitors. 3. Results and Discussion 1 3.1 Mean excitation energy of graphite 2 In the calculation of the energy loss in the degrader, the most important parameter is the mean excitation energy. In earlier literature, this parameter is often given by the experimental formula I = I 0 Z . The coefficient I 0 differs among various sources. It is often assumed to be 13 eV, and the coefficient of light atom is slightly lower—closer to 10 eV. Except for hydrogen, the mean excitation energy of each material should be measured experimentally. In the web database PSTAR of the National Institute of Standards and Technology (NIST), which calculates stopping powers, ranges, and related quantities for protons, the material properties of two kinds of graphite are given. One is amorphous carbon, which has a density close to 2.00 g/cm 3 and mean excitation energy of 81 eV. The other is graphite, which has a density of approximately 1.70 g/cm 3 and mean excitation energy of 78 eV. We set up an experiment to study the mean excitation energy of graphite, which is chosen to make the degrader. The results in Fig. 3 Fig. 3 show that the functions of the energy loss and graphite thickness under different "initial beam" energies agree with the linear rule. This reveals that the graphite blocks are uniform. In addition, we used different I values to simulate the energy loss of the degrader, from 60 to 100 eV. According to the percentage difference between the degraded energy from the measurement and simulation, we obtained a function of the degrader length for different mean excitation energies of graphite, as shown in Fig. 4 Fig. 4 . The experimental results indicate that 76 eV is the most suitable. 3.2 Ionization energy loss of degrader 2 We calculated the ionization energy loss of a proton beam passed through the degrader using Eq. (3). For the interaction between incident protons and carbon atoms, three kinds of corrections are considered: relativistic correction, shell correction, and density–effect correction. The relativistic correction begins to take effect at approximately 100 MeV and plays an important role after several GeV [ 19 19 ] . At 200 MeV, the relativistic correction is < 1%. This is equal to the density–effect correction, which arises around 300 MeV and increases to 5% at 1.5 GeV. However, the shell correction differs. It is more important at a lower energy. When the incident energy is < 4 MeV, the shell correction is approximately 2%. It is almost not working above 20 MeV, as shown in Fig. 5 Fig. 5 . As shown in Fig. 6 Fig. 6 , the stopping power of different materials for the proton beam is similar. At the zone of βγ = 3–4, the stopping power is minimized, called minimum ionizing particles (MIPs). In SC200, the maximum βγ is < 0.7. That is, the stopping power depends linearly on 1/ β 2 . The energy loss in the degrader under a 200–MeV "initial beam" is shown in Fig. 7 Fig. 7 . We compared the results of Monte Carlo simulation [ 20 20 , 21 21 ] with theoretical results and measurement data. The error is < 1.354% between the two simulation methods and is < 2.263% between the calculation and experiment. According to the results of the error analysis, the errors between the calculation results and the experiment data mainly result from the equivalent thickness of the monitors in the MLIC (each monitor is equivalent to 2.2 mm water). 3.3 Energy spread of degraded beam 2 We used Eq. (6) to calculate the energy spread after the degrader. First, the initial beam with the energy spread in the Gaussian distribution was divided into monochrome beams with different weights. Second, every monochrome beam was calculated for the distribution of the degraded beam. The final result was obtained using multi-Gaussian superimposition. As shown in Fig. 8 Fig. 8 , the calculation and simulation results were compared and agreed well. When the thickness is relatively small, the formula results are slightly lower than the Monte Carlo results, owing to the several δ electrons with high energy. The distribution has a long tail in the high energy loss zone and is more similar to the Landau distribution than to the Gaussian distribution. Fig. 9 Fig. 9 shows the standard deviation of the degraded beam with an initial beam in the 0%–1.0% energy spread. Obviously, the minus error bars are longer than plus error bars, which is caused by the tail at a low energy. A typical energy spread of the degraded beam is shown in Fig. 10 Fig. 10 . 3.4 Transmission of degrader 2 According to the empirical formula for the proton–nucleus collision, we calculated the reaction cross section of two materials that are typically used to make degraders: beryllium and carbon. We compared the results with experimental data reported by R. K. Tripathi [ 18 18 ] . The comparison is presented in Fig. 11 Fig. 11 . The agreement between the empirical formula and the experiment is very good. Obviously, the cross sections are affected by the Coulomb interaction at the lower energy end. Pauli blocking effects become increasingly important at a high energy. For confirming the loss of beam current when protons pass through degrader, we calculated the transmission of the graphite degrader using the cross–section function and compared the result with Monte Carlo simulation and experimental data, as shown in Fig. 12 Fig. 12 . The transmission error is < 1.593% between the two simulation methods and < 1.719% between the calculation and experiment. 4. Conclusion 1 We developed a degrader simulation program to calculate the energy loss and beam loss for protons passing through a degrader, according to theoretical formulae and empirical formulae. The results of the simulation of the degrader in the SC200 proton therapy facility exhibit very good agreement with experimental measurements and Monte Carlo simulation results. Thus, the rationality of the physical model is verified. After the mean excitation energy of the graphite degrader used in SC200 is adjusted, the energy–loss calculation for the degrader reproduces the measured value within 2.263% at 70–200 MeV. This means that the standard material properties are not accurate enough for high-accuracy simulations. All of the properties used in nuclear stopping simulation should be measured by experiments. According to an analysis of the main component in the ionization energy–loss analysis, we defined the scope of each component and calculated the stopping power of several common materials. Finally, we constructed the physical model of the proton–nucleus collision using the approximate formulae. The transmission of the degrader closely approximated that of the Monte Carlo simulation, with an error of <1.593%. References J.M. Slater , O.J. Archambeau , D.W. Miller , et al ., The proton treatment center at Loma Linda University Medical Center: rationale for and description of its development . Int. J. Radiat. Oncol. Biol. Phys . 22 : 383 - 389 ( 1992 ). DOI: 10.1016/0360-3016(92)90058-p http://doi.org/10.1016/0360-3016(92)90058-p The proton treatment center at Loma Linda University Medical Center: rationale for and description of its development W.P. Jones , G.P.A. Berg . Design of a beam transport system for a proton radiation therapy facility . Particle Accelerator Conference, 1999. Proceedings of the . IEEE , 1999 : 2519 - 2521 vol. 4 . DOI: 10.1109/pac.1999.792752 http://doi.org/10.1109/pac.1999.792752 Design of a beam transport system for a proton radiation therapy facility J.B. Flanz , F. Gerardi , E.L. Hubbard , Design considerations for a proton therapy beamline with an energy degrader, Fourteenth International Conference on the Application of Accelerators in Research & Industry . American Institute of Physics , 1997 : 1257 - 1260 . DOI: 10.1063/1.52693 http://doi.org/10.1063/1.52693 Design considerations for a proton therapy beamline with an energy degrader, Fourteenth International Conference on the Application of Accelerators in Research & Industry G. Karamysheva , Y. Bi , G. Chen , et al ., Compact superconducting cyclotron SC200 for proton therapy . Proceedings of Cyclotrons 2016 , Zurich, Switzerland , Sep. 2016 . Compact superconducting cyclotron SC200 for proton therapy X.H. Zeng , J.X. Zheng , Y.T. Song , et al . Beam optics study for energy selection system of SC200 superconducting proton cyclotron . Nucl. Sci. Tech ., 2018 , 29 : 134 . DOI: 10.1007/s41365-018-0462-5 http://doi.org/10.1007/s41365-018-0462-5 Beam optics study for energy selection system of SC200 superconducting proton cyclotron E. Pedroni , R. Bearpark , T. Böhringer , et al ., The PSI Gantry 2: a second generation proton scanning gantry . Med. Phys . 14 , 25 - 34 ( 2004 ). DOI: 10.1078/0939-3889-00194 http://doi.org/10.1078/0939-3889-00194 The PSI Gantry 2: a second generation proton scanning gantry M.J. van Goethem , R. van der Meer , H.W. Reist , et al ., Geant4 simulations of proton beam transport through a carbon or beryllium degrader and following a beam line . Phys. Med. Biol . 54 , 5831 - 5846 ( 2009 ). DOI: 10.1088/0031-9155/54/19/011 http://doi.org/10.1088/0031-9155/54/19/011 Geant4 simulations of proton beam transport through a carbon or beryllium degrader and following a beam line F. Stichelbaut , Y. Jongen . Properties of an energy degrader for light ions . 2014 , 4 : 272 - 275 . DOI: 10.15669/pnst.4.272 http://doi.org/10.15669/pnst.4.272 Properties of an energy degrader for light ions E.W. Cascio , S. Sarkar , A Continuously Variable Water Beam Degrader For The Radiation Test Beamline at The Francis H. Burr Proton Therapy Center . 2007 IEEE Radiation Effects Data Workshop , 2007 : 30 - 33 . DOI: 10.1109/redw.2007.4342536 http://doi.org/10.1109/redw.2007.4342536 A Continuously Variable Water Beam Degrader For The Radiation Test Beamline at The Francis H. Burr Proton Therapy Center J.A. Brennsæter , The Influence of the Energy Degrader Material for a Therapeutical Proton Beam . 2015 . The Influence of the Energy Degrader Material for a Therapeutical Proton Beam Stopping Powers and Ranges for Protons and Alpha Particles , International Commission on Radiation Units and Measurements , 1993 . Stopping Powers and Ranges for Protons and Alpha Particles D.E. Groom , Energy loss in matter by heavy particles. Particle Data Group Notes, PDG-93-06 . Energy loss in matter by heavy particles. Particle Data Group Notes, PDG-93-06 D.E. Groom , S.R. Klein , Passage of particles through matter . Eur. Phys. J. C , 15 : 163 - 173 ( 2000 ). DOI: 10.1007/bf02683419 http://doi.org/10.1007/bf02683419 Passage of particles through matter M. Behar , C.D. Denton , R.C. Fadanelli , et al ., Experimental and theoretical determination of the stopping power of ZrO2 films for protons and α-particles . Eur. Phys. J. D , 2010 , 59 : 209 - 213 . DOI: 10.1140/epjd/e2010-00164-x http://doi.org/10.1140/epjd/e2010-00164-x Experimental and theoretical determination of the stopping power of ZrO2 films for protons and α-particles C.C. Montanari , J.E. Miraglia , S. Heredia-Avalos , et al ., Calculation of energy-loss straggling of C, Al, Si, and Cu for fast H, He, and Li ions . Phys. Rev. A , 2007 , 75 : 441 - 445 . DOI: 10.1103/PhysRevA.75.022903 http://doi.org/10.1103/PhysRevA.75.022903 Calculation of energy-loss straggling of C, Al, Si, and Cu for fast H, He, and Li ions Y. Kido , T. Hioki . Measurements of energy loss and straggling for fast H+ in metals and their compounds by means of a nuclear resonant reaction . Phys. Rev. B , 1983 , 27 : 2667 - 2673 . DOI: 10.1103/physrevb.27.2667 http://doi.org/10.1103/physrevb.27.2667 Measurements of energy loss and straggling for fast H+ in metals and their compounds by means of a nuclear resonant reaction M. Behar , R.C. Fadanelli , I. Abril , et al ., Energy-loss straggling study of proton and alpha-particle beams incident onto ZrO2, and Al2O3 films . Eur. Phys. J. D , 2011 , 64 : 297 - 301 . DOI: 10.1140/epjd/e2011-20352-4 http://doi.org/10.1140/epjd/e2011-20352-4 Energy-loss straggling study of proton and alpha-particle beams incident onto ZrO2, and Al2O3 films R.K. Tripathi , F.A. Cucinotta , J.W. Wilson . Universal Parameterization of Absorption Cross Sections . Nucl. Instrum. Methods B . 117 : 347 ( 1996 ). DOI: 10.1016/0168-583x(96)00331-x http://doi.org/10.1016/0168-583x(96)00331-x Universal Parameterization of Absorption Cross Sections H. Paul , A comparison of recent stopping power tables for light and medium-heavy ions with experimental data, and applications to radiotherapy dosimetry . Nucl. Instrum. Methods 247 : 166 - 172 ( 2006 ). DOI: 10.1016/j.nimb.2006.01.059 http://doi.org/10.1016/j.nimb.2006.01.059 A comparison of recent stopping power tables for light and medium-heavy ions with experimental data, and applications to radiotherapy dosimetry T. Aso , A. Kimura , S. Tanaka , et al . Verification of the dose distributions with GEANT4 simulation for proton therapy . IEEE T. Nucl. Sci ., 52 : 896 - 901 ( 2005 ),. DOI: 10.1109/nssmic.2004.1462685 http://doi.org/10.1109/nssmic.2004.1462685 Verification of the dose distributions with GEANT4 simulation for proton therapy H. Paganetti , H. Jiang , S.Y. Lee , et al ., Accurate Monte Carlo simulations for nozzle design, commissioning and quality assurance for a proton radiation therapy facility . Med. Phys ., 31 : 2107 ( 2004 ). DOI: 10.1118/1.1762792 http://doi.org/10.1118/1.1762792 Accurate Monte Carlo simulations for nozzle design, commissioning and quality assurance for a proton radiation therapy facility The authors are grateful for the financial support from Hefei CAS Ion Medical and Technical Devices Co., Ltd. 10.1007/s41365-018-0526-6 This work was supported in part by the National Natural Science Foundation of China (No. 51525703). 2017-12-25 2018-04-21 2018-05-13 2019-01-02 2019-01-01 2021-04-10 © 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. 2019 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-018-0526-6 http://dx.doi.org/10.1007/s41365-018-0526-6 . 2019 Nuclear Science and Techniques 01 30