Effect of magnetic field distribution on ECR ion source discharge

Effect of magnetic field distribution on ECR ion source discharge 1 corresp first-author Corresponding author, Yu-Guo Liu, liuyg10@lzu.edu.cn Corresponding author, Yu-Guo Liu, liuyg10@lzu.edu.cn liuyg10@lzu.edu.cn 1 1 1 1 1 Institute of Nuclear Physics and Chemistry, CAEP, Mianyang 621900, China Institute of Nuclear Physics and Chemistry, CAEP, Mianyang 621900, China An electron cyclotron resonance ion source model that use a coaxial antenna to inject the RF power was simulated, and the corresponding hydrogen plasma densities of the ion source were calculated at different magnetic field distributions via a fluid model. The results show that most of the microwave power is absorbed around the resonance surface near the antenna when the plasma density is below the cut off density of 7.4 × 10 16 m -3 for 2.45 GHz microwaves, and the simulation results also show that the plasma density is strongly influenced by the position and shape of the resonance surface, where a larger resonance surface would improve the plasma density. ECR ion source Magnetic field distribution Plasma simulation 1. INTRODUCTION 1 The electron cyclotron resonance (ECR) ion source was first proposed by Geller et al. [ 1 1 ] in the late 1960’s, and high-current ECR ion sources have been rapidly developed for various purposes since the 1990s [ 2 2 - 6 6 2-6 ] . The generation and distribution of plasma in an ECR ion source are strongly influenced by the magnetic field distribution, and use of a suitable magnetic field configuration could improve the discharge efficiency [ 7 7 , 8 8 ] . Previous works showed that the beam current of an ECR ion source is closely related to the axial magnetic field [ 9 9 , 10 10 ] , and the optimal magnet configuration is generally represented by a magnetic field distribution along the center axis of the discharge chamber. However, for ion sources with different magnets or solenoids, the optimal magnetic field distributions along the center axis of the discharge chamber may not be the same [ 2 2 , 6 6 , 10 10 ] . Therefore, the influence of the resonance surface on the plasma density rather than the magnet field strength along the center axis of the discharge chamber was studied. An ECR ion source model with a coaxial antenna that injects the RF power was simulated. The hydrogen plasma density of the source was calculated at different magnetic field distributions via a fluid model [ 11 11 , 12 12 ] . The influence of the shape and position of the resonance surface on the plasma density was analyzed. 2. THE SIMULATION MODEL 1 Figure 1 Figure 1 shows a schematic of the structure of the simulation model. The time evolution of the microwaves and plasma is simulated in azimuthally symmetric cylindrical coordinates, and the direction of the antenna in the calculation model is along the z _axis, as shown in Fig. 1 Fig. 1 . The 2.45 GHz RF power is axially introduced to the discharge chamber via an antenna 8 mm in diameter. The 2D axisymmetric model represents the antenna in the discharge chamber and the permanent ring magnets. The discharge chamber is a cylinder approximately 50 mm in diameter and 50 mm in length. The magnetic field of the ion source is generated by two NdFeB permanent ring magnets, and the direction of magnetization is along the positive z _axis. The parameters of the permanent ring magnets are shown in Table 1 Table 1 . 1 1 Inner diameter (mm) 1 1 Outer diameter (mm) 1 1 Height (mm) 1 1 Remanence (T) 1 1 μ r 1 1 66 1 1 110 1 1 20 1 1 1.23 1 1 1.05 The magnetic field B of the permanent ring magnets is obtained by solving the Maxwell_Ampere equations: <math display="block" id="M1"><mrow><mo>∇</mo><mo>×</mo><mrow><mo>(</mo><mrow><msubsup><mi>μ</mi><mn>0</mn><mrow><mo>−</mo><mn>1</mn></mrow></msubsup><msubsup><mi>μ</mi><mtext>r</mtext><mrow><mtext>-1</mtext></mrow></msubsup><mtext>B</mtext></mrow><mo>)</mo></mrow><mo>−</mo><mi>σ</mi><mtext>v</mtext><mo>×</mo><mtext>B</mtext><mo>=</mo><msub><mtext>J</mtext><mtext>m</mtext></msub></mrow></math> <math display="block" id="M2"><mrow><mtext>B</mtext><mo>=</mo><mo>∇</mo><mo>×</mo><mtext>A</mtext></mrow></math> where μ is the permeability, J m is the azimuthal current at the surface of the ring magnets, and A is the magnetic vector potential. The microwave fields were obtained by solving the wave equation: <math display="block" id="M3"><mrow><mo>∇</mo><mo>×</mo><msubsup><mi>μ</mi><mn>0</mn><mrow><mo>−</mo><mn>1</mn></mrow></msubsup><mrow><mo>(</mo><mrow><mo>∇</mo><mo>×</mo><mtext>E</mtext></mrow><mo>)</mo></mrow><mo>−</mo><msubsup><mi>k</mi><mn>0</mn><mn>2</mn></msubsup><mrow><mo>(</mo><mrow><msub><mi>ε</mi><mi>r</mi></msub><mo>−</mo><mfrac><mrow><mi>j</mi><mi>σ</mi></mrow><mrow><mi>ω</mi><msub><mi>ε</mi><mn>0</mn></msub></mrow></mfrac></mrow><mo>)</mo></mrow><mo>⋅</mo><mtext>E</mtext><mo>=</mo><mn>0</mn><mo>,</mo></mrow></math> where E is the electric field, μ 0 is the permeability of free space, ε 0 is the permittivity of free space, ε r is the relatively permittivity, k 0 is the wave number, ω is the angle frequency of microwave, σ is the full tensor of the plasma conductivity. The electron density n e is described by the follow equations: <math display="block" id="M4"><mrow><mfrac><mo>∂</mo><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><mrow><mo>(</mo><mrow><msub><mi>n</mi><mtext>e</mtext></msub></mrow><mo>)</mo></mrow><mo>+</mo><mo>∇</mo><mo>⋅</mo><mrow><mo>[</mo><mrow><mo>−</mo><msub><mi>n</mi><mtext>e</mtext></msub><mrow><mo>(</mo><mrow><msub><mtext>μ</mtext><mtext>e</mtext></msub><mo>×</mo><mtext>E</mtext></mrow><mo>)</mo></mrow><mo>−</mo><msub><mi>T</mi><mtext>e</mtext></msub><msub><mtext>μ</mtext><mtext>e</mtext></msub><mo>⋅</mo><mo>∇</mo><msub><mi>n</mi><mtext>e</mtext></msub></mrow><mo>]</mo></mrow><mo>=</mo><msub><mi>R</mi><mtext>e</mtext></msub><mo>,</mo></mrow></math> <math display="block" id="M5"><mrow><mfrac><mo>∂</mo><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><mrow><mo>(</mo><mrow><msub><mi>n</mi><mtext>e</mtext></msub></mrow><mo>)</mo></mrow><mo>+</mo><mo>∇</mo><mo>⋅</mo><mrow><mo>[</mo><mrow><mo>−</mo><msub><mi>n</mi><mtext>e</mtext></msub><mrow><mo>(</mo><mrow><mfrac><mn>5</mn><mn>3</mn></mfrac><msub><mtext>μ</mtext><mtext>e</mtext></msub><mo>⋅</mo><mtext>E</mtext></mrow><mo>)</mo></mrow><mo>−</mo><msub><mi>T</mi><mtext>e</mtext></msub><msub><mtext>μ</mtext><mtext>ε</mtext></msub><mo>×</mo><mo>∇</mo><msub><mi>n</mi><mtext>e</mtext></msub></mrow><mo>]</mo></mrow><mo>+</mo><mtext>E</mtext><mo>×</mo><msub><mtext>Γ</mtext><mtext>e</mtext></msub><mo>=</mo><msub><mi>R</mi><mtext>ε</mtext></msub><mo>,</mo></mrow></math> <math display="block" id="M6"><mrow><msub><mi>Γ</mi><mtext>e</mtext></msub><mo>=</mo><mo>−</mo><msub><mi>n</mi><mtext>e</mtext></msub><mtext>E</mtext><mo>−</mo><msub><mi>T</mi><mtext>e</mtext></msub><mo>∇</mo><mo>⋅</mo><msub><mtext>μ</mtext><mtext>e</mtext></msub></mrow></math> where T e is the electron temperature, R e and R ε is the electron source and energy loss due to inelastic collisions, respectively; and <math id="M7"><mrow><mtext>E</mtext><mo>⋅</mo><msub><mi>Γ</mi><mtext>e</mtext></msub></mrow></math> is the heating term representing electron motion and heating of electrons by microwaves. μ ε is the energy mobility, and μ e is the tensor of the electron mobility. The working gas is hydrogen in the simulation, and the reaction processes that are taken into account are listed in Table 2 Table 2 . The cross section of each reaction is shown in Fig. 2 Fig. 2 [ 13 13 - 15 15 13-15 ] . 1 1 Reaction type 1 1 Reaction 1 1   1 1 Elastic collision 1 1 e+H 2 →e+H 2 1 1 (I) 1 1   1 1 e+H→e+ H 1 1 (II) 1 1 Dissociative excitation 1 1 e+H 2 →e+ H+H 1 1 (III) 1 1   1 1 e+ <math id="M8"><mrow><msubsup><mtext>H</mtext><mn>2</mn><mo>+</mo></msubsup></mrow></math> →e+ H+H + 1 1 (IV) 1 1   1 1 e+ <math id="M9"><mrow><msubsup><mtext>H</mtext><mn>3</mn><mo>+</mo></msubsup></mrow></math> →e+2H+H + 1 1 (V) 1 1   1 1 e+H+ →e+H+ <math id="M10"><mrow><msubsup><mtext>H</mtext><mn>2</mn><mo>+</mo></msubsup></mrow></math> 1 1 (VI) 1 1 Dissociative combination 1 1 e+ <math id="M11"><mrow><msubsup><mtext>H</mtext><mn>2</mn><mo>+</mo></msubsup></mrow></math> →H+H 1 1 (VII) 1 1   1 1 e+ <math id="M12"><mrow><msubsup><mtext>H</mtext><mn>3</mn><mo>+</mo></msubsup></mrow></math> →H+H+H 1 1 (VIII) 1 1 Excitation 1 1 e+H→e+H * 1 1 (IX) 1 1 <math id="M13"><mrow><msubsup><mtext>H</mtext><mn>3</mn><mo>+</mo></msubsup></mrow></math> generation 1 1 H 2 + <math id="M14"><mrow><msubsup><mtext>H</mtext><mn>2</mn><mo>+</mo></msubsup></mrow></math> → <math id="M15"><mrow><msubsup><mtext>H</mtext><mn>3</mn><mo>+</mo></msubsup></mrow></math> +H 1 1 (X) 1 1 Ionization 1 1 e+H→2e+H + 1 1 (XI) 1 1   1 1 e+H * →2e+ H + 1 1 (XII) 1 1   1 1 e+H 2 →2e+ <math id="M16"><mrow><msubsup><mtext>H</mtext><mn>2</mn><mo>+</mo></msubsup></mrow></math> 1 1 (XIII) 1 1 Reaction with wall 1 1 H + +wall→1/2H 2 1 1   1 1 <math id="M17"><mrow><msubsup><mtext>H</mtext><mn>2</mn><mo>+</mo></msubsup></mrow></math> +wall→1/2H 2 1 1   1 1 H+wall→1/2H 2 1 1   1 1 <math id="M18"><mrow><msubsup><mtext>H</mtext><mn>3</mn><mo>+</mo></msubsup></mrow></math> +wall→H+H 2 1 1   3. RESULTS AND DISCUSSION 1 The typical solution results are shown in Fig. 3 Fig. 3 , where the gas pressure is 1 Pa, the microwave power is 100 W, and the gap length D between two permanent ring magnets is 24 mm. Figure 3(a) Figure 3(a) shows that most of the microwave power is absorbed in the resonance surface ( B = 875 G) near the top of the antenna because resonance between the electrons and the microwave field occurs in this region; as a result, the microwave field dose most of the work on the electrons to ionize the hydrogen gas near the resonance surface. The corresponding maximum electron density is approximately 5 × 10 16 m -3 , which is below the cut off density of 7.4 × 10 16 m -3 for 2.45 GHz microwaves [ 16 16 ] , and the maximum electron temperature is approximately 10 eV. The magnetic field in the discharge chamber could be changed by moving two permanent ring magnets along the z _axis with a fixed gap length D or changing the gap length D between two permanent ring magnets. The magnetic field strength along the source z _axis for different gap lengths D is shown in Fig. 4 Fig. 4 . The minimum magnetic field strength decreases from 1100 to 550 G when the gap length D increases from 20 to 30 mm, and the corresponding maximum magnetic field strength decrease from 1200 to 1100 G. The distributions of H + ions and maximum H + ion densities are shown in Figs. 5 Figs. 5 and 6 6 , respectively, as the gap length D increase from 20 to 30 mm. Figure 5 Figure 5 shows that the area of the resonance surface in the discharge chamber is increased with increasing permanent ring magnets gap length D , and the minimum distance between the antenna and the resonance surface is decreased at the same time. Figure 6 Figure 6 shows that the maximum density of H + ions first increases with increasing gap length D and then decreases when D is larger than 24 mm, and the maximum H + ion density is between 2.6 × 10 16 m -3 and 2.8 × 10 16 m -3 . In addition, Fig. 6 Fig. 6 also shows that the fluctuation of the maximum densities of H+ 2 ions and H+ 3 ions is smaller than that of H + ions when the magnet gap length D is changed. The distributions of the H + ion densities are also calculated as the ring magnets moved along the z _axis with a fixed gap length D of 24 mm; the midpoint of the ring magnets is placed at the origin of the cylindrical coordinates, as shown in Fig. 1 Fig. 1 , when the corresponding movement distance M is zero. The density distributions and the corresponding maximum densities of H + ions are shown in Figs. 7 Figs. 7 and 8 8 , respectively, as the midpoint of the ring magnets moves along the z _axis from -4 to 7 mm. Figure 7 Figure 7 shows that the minimum distance between the antenna and the resonance surface is decreased as the ring magnets move. Figure 8 Figure 8 shows that the maximum density of H + ions is sensitive to the distance between the resonance surface and the antenna. The maximum density of H + ions increases from 2.3 × 10 16 m -3 to 3.2 × 10 16 m -3 when the midpoint of the ring magnets moves along the z _axis from -4 to 4 mm and then decrease to 2.78 × 10 16 m -3 when the midpoint of the ring magnets move from 4 to 7 mm along the z _axis. The normal electric field component, which is perpendicular to the axial direction is calculated and shown in Fig. 9 Fig. 9 , which shows that the normal electric field strength is attenuated exponentially along the z _axis [ 17 17 ] . The electron density distributions are shown in Fig. 10 Fig. 10 . Figure 10(a) Figure 10(a) shows the electron density distributions along the z _axis as the ring magnets gap length D is changed, and Fig. 10(b) Fig. 10(b) shows the electron density distributions along the z _axis as the ring magnets move along the z _axis with a fixed gap length D of 24 mm. Figure 10 Figure 10 shows that regardless of whether the gap length D of the ring magnets is changed or the ring magnets are moved along the z _axis with a fixed gap length D , the electron density is highest in the central area and then decrease to zero on both sides of the z _axis. Figures 5 Figures 5 and 7 7 show that regardless of whether the gap length D of the ring magnets is increased or the ring magnets are moved toward the antenna along the z _axis with a fixed gap length D , the minimum distance between the resonance surface and the top of the antenna is reduced. Therefore, the normal electric field strength along the z _axis would increase according to figure 9 figure 9 , and the deposited microwave power could be expressed as <math display="block" id="M19"><mrow><mi>P</mi><mo>=</mo><mo>−</mo><mfrac><mrow><mi>e</mi><mi>ω</mi></mrow><mrow><mn>2</mn><mi>π</mi></mrow></mfrac><mstyle displaystyle="true"><mrow><msubsup><mo>∫</mo><mn>0</mn><mrow><mfrac><mrow><mn>2</mn><mi>π</mi></mrow><mi>ω</mi></mfrac></mrow></msubsup><mrow><mtext>E</mtext><mo>×</mo><msub><mtext>v</mtext><mtext>e</mtext></msub><mtext>d</mtext><mi>t</mi></mrow></mrow></mstyle></mrow></math> where E is the electric field and v e is the electron mean velocity. Consequently the deposited microwave power would increase with decreasing the distance between the resonance surface and the antenna; on the other side, the electron density near the resonance surface would first increase, as shown in Fig. 10 Fig. 10 , when the distance between the resonance surface and the antenna is reduced. Thus, the density of H + ions would increase initially. When the minimum distance between the resonance surface and the top of the antenna is further reduced, the deposited microwave power near the resonance surface increases continuously as the electron density near the resonance surface decreases. As a result, the number of electrons used to ionize the hydrogen gas decreases, and the density of H + ions would decrease as a consequence. Therefore, the graphs in Figs. 6 Figs. 6 and 8 8 both have a tendency to rise initially and then fall. However, the shape of the resonance surface also changes with increasing ring magnets gap length D , which may be the main reason for the discrepancy between Figs. 6 Figs. 6 and 8 8 . The maximum density of H + ions under different magnetic field distributions, which was obtained by replacing the abscissas of Figs. 6 Figs. 6 and 8 8 with the minimum distance between the resonance surface and the top of the antenna is shown in Fig. 11 Fig. 11 . Figure 11 Figure 11 shows that the maximum density of H + ions has a tendency to rise at first and then fall when the minimum distance between the resonance surface and the antenna is reduced. Further, the maximum H + ions densities under the same abscissa values are not equal in Fig. 11 Fig. 11 ; the variation of the H + ion density caused by moving the magnets rings with a fixed gap length D is larger than changing the ring magnets gap length D , indicating that the maximum H + ion density is associated with not only the position but also the shape of the resonance surface. The distributions of the H + ion density under different resonance surface shapes are shown in Fig.12 Fig.12 , where the resonance surface shapes were changed by adjusting the size and gap length of the magnets. The minimum distance between the resonance surface and the antenna is 15 mm as the ring magnets are moved along the z _axis. The maximum H + ion density in figure 12.a figure 12.a is 2.8 × 10 16 m -3 , and the Fig.12(b) Fig.12(b) is 1.25 × 10 16 m -3 , whereas that in Fig.12(c) Fig.12(c) is only 3.39 × 10 9 m -3 . Although the minimum distance between the resonance surface and the top of the antenna is the same in each figure, the corresponding H + ion density is quite different. The resonance surface in Fig.12(c) Fig.12(c) is located at the side of the discharge chamber, and the resonance surface area is also smallest compared with those in the other figures; therefore the H + ion density is lowest in Fig. 12(c) Fig. 12(c) . The resonance surfaces in Fig. 12(a) Fig. 12(a) and Fig.12(b) Fig.12(b) are both located in the center of the discharge chamber; however, the resonance surface area in Fig.12(a) Fig.12(a) is larger than that in Fig.12(b) Fig.12(b) . As a result the maximum H + ion density of Fig.12(a) Fig.12(a) is approximately twice that in Fig.12(b) Fig.12(b) . This result shows that the area of the resonance surface has an important effect on the H + ion density; a large flat resonance surface would improve the H + ion density. 4. CONCLUSION 1 A permanent magnet ECR ion source model was simulated via a fluid model using hydrogen as the working gas. The simulation results show that the maximum electron density is approximately 5 × 10 16 m -3 , and the maximum electron temperature is approximately 10 eV. The simulation results also show that the most of the microwave power is deposited around the resonance surface near the antenna when the plasma density is below the cut off density (7.4 × 10 16 m -3 ) for 2.45 GHz microwaves. The H + ion densities under different magnetic field configurations were calculated; the results show that the H + ion density is strongly influenced by the shape and position of the resonance surface. The H + ion density is sensitive to the distance between the antenna and the resonance surface, and a larger resonance surface would improve the H + ion density. REFERENCES R. Geller , B. Jacquot , C. Jacquot , Suppression of instabilities and collision limited plasma accumulation in the pleiade mirror device , Phys. Lett.A . 30 , 80 - 81 ( 1969 ). doi: 10.1016/0375-9601(69)91146-3 http://doi.org/10.1016/0375-9601(69)91146-3 Suppression of instabilities and collision limited plasma accumulation in the pleiade mirror device J. Sherman , A. Arvin , L. Hansborough et al ., Development of a 130-mA, 75-kV high voltage column for high-intensity dc proton injectors , Rev. Sci. Instrum . 69 , 1017 - 1019 ( 1998 ). doi: 10.1063/1.1148707 http://doi.org/10.1063/1.1148707 Development of a 130-mA, 75-kV high voltage column for high-intensity dc proton injectors M. Delaunay , E. Touchais . Electron cyclotron resonance plasma ion source for material depositions , Rev. Sci. Instrum . 69 , 2320 - 2324 ( 1998 ). doi: 10.1063/1.1148938 http://doi.org/10.1063/1.1148938 Electron cyclotron resonance plasma ion source for material depositions R. Gobin , P. Y. Beauvais , D. Bogard , et al ., High intensity ECR ion source (H+, D+, H-) developments at CEA/Saclay , Rev. Sci. Instrum . 73 , 922 - 924 ( 2002 ). doi: 10.1063/1.1428783 http://doi.org/10.1063/1.1428783 High intensity ECR ion source (H+, D+, H-) developments at CEA/Saclay B. Cui , R. Wang , Y. Ma , et al ., A high intensity microwave ion source for high current RFQ , Rev. Sci. Instrum . 75 , 1457 - 1459 ( 2004 ). doi: 10.1063/1.1690460 http://doi.org/10.1063/1.1690460 A high intensity microwave ion source for high current RFQ Z. Z. Song , S. X. Peng , J. X. Yu et al ., Minipermanent magnet high-current microwave ion source , Rev. Sci. Instrum . 77 , 03A305-1-03A305-3 ( 2006 ). doi: 10.1063/1.2150807 http://doi.org/10.1063/1.2150807 Minipermanent magnet high-current microwave ion source O. A. Popov , Effects of magnetic field and microwave power on electron cyclotron resonance-type plasma characteristics , J.Vac.Sci.Technol.A 9 , 711 - 716 ( 1991 ). doi: 10.1116/1.577349 http://doi.org/10.1116/1.577349 Effects of magnetic field and microwave power on electron cyclotron resonance-type plasma characteristics T. Taylor , J.F. Mouris , An advanced high-current low-emittance dc microwave proton source , Nucl. Instrum. Method. A , 336 , 1 - 5 ( 1993 ). doi: 10.1016/0168-9002(93)91074-W http://doi.org/10.1016/0168-9002(93)91074-W An advanced high-current low-emittance dc microwave proton source J. S. C. Wills , R. A. Lewis , J. Diserens , et al ., A compact high-current microwave-driven ion source , Rev. Sci. Instrum . 69 , 65 - 68 ( 1998 ). doi: 10.1063/1.1148479 http://doi.org/10.1063/1.1148479 A compact high-current microwave-driven ion source S. X. Peng , R. Xu , J. Zhao et al ., The influence of magnetic field configuration on an electron cyclotron resonance ion source a , Rev. Sci. Instrum . 79 , 02A310-1-02A310-3 ( 2008 ). doi: 10.1063/1.2812343 http://doi.org/10.1063/1.2812343 The influence of magnetic field configuration on an electron cyclotron resonance ion source a G. J. M. Hagelaar , K. Makasheva , L. Garrigues , Boeuf, Modelling of a dipolar microwave plasma sustained by electron cyclotron resonance , J. Phys. D: Appl. Phys . 42 , 1 - 12 ( 2009 ). doi: 10.1088/0022-3727/42/19/194019 http://doi.org/10.1088/0022-3727/42/19/194019 Boeuf, Modelling of a dipolar microwave plasma sustained by electron cyclotron resonance D. Hemmers , M. David , H. Kempkens , et al ., Self-consistent modelling of overdense plasmas in ECR discharges , J. Physics. D: Appl. Phys . 31 2155 - 2164 ( 1998 ). doi: 10.1088/0022-3727/31/17/013 http://doi.org/10.1088/0022-3727/31/17/013 Self-consistent modelling of overdense plasmas in ECR discharges H. Zhang , C. Wan . Ion Source and High Power Neutral Beam Source , ( Atomic Energy Press , Beijing , 1985 ), pp. 218 - 221 . (in Chinese) Ion Source and High Power Neutral Beam Source http://jila.colorado.edu/~avp/collision_data/electronneutral/ELECTRON.TXT http://jila.colorado.edu/~avp/collision_data/electronneutral/ELECTRON.TXT https://nl.lxcat.net/data/set_type.php https://nl.lxcat.net/data/set_type.php M.A. Lieberman , A.J. Lichtenberg . Principles of Plasma Discharge and Materials Processing , 2 nd edn. ( John Wiley & Sons , New Jersey , 2005 ), pp. 94 . Principles of Plasma Discharge and Materials Processing M. Feix , Excitation level of plasma oscillations , Phys.Lett . 9 , 123 - 125 ( 1964 ). doi: 10.1016/0031-9163(64)90106-4 http://doi.org/10.1016/0031-9163(64)90106-4 Excitation level of plasma oscillations 10.1007/s41365-018-0464-3 2017-10-20 2018-01-05 2018-01-12 2018-07-25 2018-09-01 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. 2018 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-0464-3 http://dx.doi.org/10.1007/s41365-018-0464-3 . 2018 Nuclear Science and Techniques 09 29