Transport characteristics of space charge dominated multi- species deuterium beam in electrostatic low energy beam line

LOW ENERGY ACCELERATOR, RAY AND APPLICATIONS

Transport characteristics of space charge dominated multi- species deuterium beam in electrostatic low energy beam line

Xiao-Long Lu，

Yu Zhang，

Jun-Run Wang，

Zhi-Wu Huang，

Zhan-Wen Ma，

Yao Yang，

Ze-En Yao

Nuclear Science and Techniques

Vol.29, No.4

Article number 51

Published in print 01 Apr 2018

Available online 14 Mar 2018

DOI：10.1007/s41365-018-0384-2

40200

The transport characteristics of a space charge dominated multi-species deuterium beam consisting of $D_{1}^{+}$ , $D_{2}^{+}$ , and $D_{3}^{+}$ particles in an electrostatic low energy beam line are studied. First, the envelope equations of the primary $D_{1}^{+}$ beam are derived considering the space charge effects caused by all particles. Second, the evolution of the envelope of the multi-species deuterium beam is simulated using the PIC code TRACK, with the results showing a significant effect of the unwanted beam on the transport of the primary beam. Finally, different injected beam parameters are used to study beam matching, and a new beam extraction system for the existing duoplasmatron source is designed to obtain the ideal injected beam parameters that allow a $D_{1}^{+}$ beam of up to 50 mA to pass unobstructed through the electrostatic low energy beam transport line in the presence of an unwanted ( $D_{2}^{+}$ , $D_{3}^{+}$ ) beam of 20 mA; at the same time, distortions of the beam emittance and particle distributions are observed.

Envelope equationSpace charge effectsMulti-species beamElectrostatic LEBT

I. INTRODUCTION

High intensity accelerators such as the High Intensity Heavy Ion Accelerator Facility (HIAF), Accelerator Driven Sub-critical System (ADS), Intense Neutron Generator, etc., are required for scientific research and industrial applications. High intensity low energy beam transport (LEBT) lines are crucial for the development of high intensity accelerators; however, the associated strong space charge effects present challenges to the transport of high intensity beams in the LEBTs^[1,2]. The LEBT design is normally based on electrostatic or magnetostatic focusing. The advantage of a magnetostatic LEBT is that the beam can be fully neutralized by the residual gas present in the line, while in an electrostatic LEBT the beam is fully un-neutralized, resulting in strong beam divergence. Therefore, magnetostatic LEBTs are usually adopted. However, electrostatic LEBTs can be very compact because focusing of a high intensity low energy beam using an electric field is more effective than using a magnetic field^[3]. Hence, electrostatic LEBTs are very attractive for some high intensity accelerators used in education laboratories and industrial applications, such as high intensity ion implanters, intense neutron generators based on accelerators, and others^[4].

In our laboratory at Lanzhou University, an intense DT/DD neutron generator based on a Cockcroft–Walton accelerator is under construction, and the electrostatic LEBT design is adopted. As shown in figure 1, the designed electrostatic LEBT can be viewed as a system of two electrostatic lenses (A1 and A2) combined with a drift tube^[5,6]. The multi-species beam contains $D_{1}^{+}$ , $D_{2}^{+}$ , and $D_{3}^{+}$ particles directly injected into the electrostatic LEBT; therefore, the transport of the primary particles $D_{1}^{+}$ is affected by the space charge effects not only from the $D_{1}^{+}$ particles themselves but also from the unwanted $D_{2}^{+}$ and $D_{3}^{+}$ particles^[7].

Fig. 1

(Color online) Electrostatic LEBT for a multi-species deuterium beam; the beam extraction system of the duoplasmatron source is indicated by the black rectangle; the voltage values are shown; the pink lines represent equipotential surfaces

In this paper, the envelope equations of the primary beam in an electrostatic LEBT are deduced considering the space charge effects of both the primary and the unwanted beams. In addition, PIC simulations are performed to study the influences of the unwanted beam on the transport of the primary beam. Finally, the influence of the parameters of the beam injected from the duoplasmatron source on the transport of the multi-species beam are studied by the PIC simulations, and a new beam extraction system for the duoplasmatron source is designed.

II. THEORETICAL ANALYSIS

A. Beam space charge field

We consider an axisymmetric high intensity multi-species continuous beam propagating through an axisymmetric electrostatic LEBT. We will use a cylindrical coordinate system $(r, z, θ)$ . We consider a uniform charge density profile, which can be expressed as^[8]

ρ_{i} (r, s) = \frac{λ_{i}}{π r_{b i}^{2} (s)},

(1)

where $λ_{i}$ is the linear charge density, $r_{b i}$ is the beam radius, and $s = z$ is the axial coordinate along the beam direction.

According to Gauss’s theorem, $\nabla \cdot E = \frac{ρ}{ε_{0}}$ , the x-component of the space charge electric field can be expressed as^[9]

E_{x i} (x) = {_{\frac{I_{i}}{2 π ε_{0} β_{i} c} \frac{x}{r^{2}}, x > r_{b i}}^{\frac{I_{i}}{2 π ε_{0} β_{i} c r_{b i}^{2}} x, x \leq r_{b i}},

(2)

where $I_{i}$ is the beam current of the species $i$ , $β_{i} c \approx v_{z i}$ is the particle axial velocity in the paraxial approximation, $r = \sqrt{x^{2} + y^{2}}$ is the distance from the beam axis in the radial direction, and $x$ and $y$ are the transverse components in the horizontal and vertical directions, respectively. The constants $ε_{0}$ , $β_{i}$ , and $c$ are the permittivity, the relativistic parameter, and the speed of light in vacuum, respectively.

B. Envelope equations of a multi-species beam

Assume that the multi-species beam contains particles of two species, a $D_{1}^{+}$ beam with a radius of $r_{1}$ and $D_{2}^{+}$ beam with a radius of $r_{2}$ . The non-relativistic transverse equation of motion describing a $D_{1}^{+}$ particle evolving in an applied electric accelerating field and the beam space charge field can be written as^[10]

m_{1} \ddot{x} = q (E_{x}^{appl} + E_{D_{1} x}^{self} + E_{D_{2} x}^{self})

(3)

Here, $\ddot{x} = v_{_{z 1}}^{2} x'' + v_{z 1} v_{_{z 1}}^{'} x'$ ; primes denote derivatives with respect to the axial coordinate $z$ ; $v_{z 1} \approx β_{1} c$ is the axial velocity of the $D_{1}^{+}$ particle; $m_{1}$ and $q$ are the mass and charge of the $D_{1}^{+}$ particle, respectively; $E_{x}^{appl}$ , $E_{D_{1} x}^{self}$ , and $E_{D_{2} x}^{self}$ are the applied electric field, the space charge field of the $D_{1}^{+}$ beam, and the space charge field of the $D_{2}^{+}$ beam along the x-direction, respectively.

In the paraxial approximation, $E_{x}^{appl} = - \frac{1}{2} \frac{\partial^{2} φ}{\partial z^{2}} x = - \frac{φ''}{2} x$ , $ϕ$ is the normalized potential corresponding to the applied on-axis electric field, $q ϕ = \frac{1}{2} m_{1} β_{1}^{2} c^{2}$ , while $E_{D_{1} x}^{self}$ and $E_{D_{2} x}^{self}$ can be calculated with Eq.(2). Then, the equations of motion of the $D_{1}^{+}$ particle in the x-direction can be written as^[11]

{_{x^{″} + \frac{ϕ^{'}}{2 ϕ} x^{'} + \frac{ϕ^{″}}{4 ϕ} x - \frac{I_{1}}{4 ϕ π ε_{0} β_{1} c} \frac{x}{r_{1}^{2}} - \frac{I_{2}}{4 ϕ π ε_{0} β_{2} c} \frac{x}{r^{2}} = 0 \begin{matrix} (r_{1} > r_{2}) \end{matrix}}^{x^{″} + \frac{ϕ^{'}}{2 ϕ} x^{'} + \frac{ϕ^{″}}{4 ϕ} x - \frac{I_{1}}{4 ϕ π ε_{0} β_{1} c r_{1}^{2}} x - \frac{I_{2}}{4 ϕ π ε_{0} β_{2} c r_{2}^{2}} x = 0 \begin{matrix} (r_{1} \leq r_{2}) \end{matrix}} .

(4)

In order to derive the envelope equations, we let $R_{x}^{2} = X_{b}^{2} = 2 〈 x^{2} 〉$ , where $X_{b}$ is the beam envelope in the x-direction, and $〈 x^{2} 〉 = \frac{\int_{0}^{r_{b i}} \int_{0}^{2 π} x^{2} ρ (r, s) r d r d θ}{\int_{0}^{r_{b i}} \int_{0}^{2 π} ρ (r, s) r d r d θ}$ is the transverse average of the particles^[8-11]. Then,

{X^{″}}_{b} = \frac{2 〈 x x^{″} 〉}{X_{b}} + \frac{4 〈 x^{2} 〉 〈 {x^{'}}^{2} 〉 - 4 {〈 x x^{'} 〉}^{2}}{X_{b}^{3}}

(5)

For $r_{1} \leq r_{2}$ , equation (4) can be transformed into equation (6) by multiplying by $x$

and averaging over the distribution of the particles,

〈 x x^{″} 〉 + \frac{ϕ^{'}}{2 ϕ} 〈 x x^{'} 〉 + \frac{ϕ^{″}}{4 ϕ} 〈 x^{2} 〉 - \frac{I_{1}}{4 ϕ π ε_{0} β_{1} c r_{1}^{2}} 〈 x^{2} 〉 - \frac{I_{2}}{4 ϕ π ε_{0} β_{2} c r_{2}^{2}} 〈 x^{2} 〉 = 0

(6)

Thus, combining Eq.(5) with Eq.(6), the envelope for the $D_{1}^{+}$ species of the multi-species beam can be described as

{X^{″}}_{b} + \frac{ϕ^{'}}{2 ϕ} {X^{'}}_{b} + \frac{ϕ^{″}}{4 ϕ} X_{b} - \frac{I_{1}}{4 ϕ π ε_{0} β_{1} c} \frac{1}{X_{b}} - \frac{I_{2}}{4 ϕ π ε_{0} β_{2} c} \frac{1}{X_{bd 2}^{2}} X_{b} - \frac{ε_{x d 1}^{2}}{X_{b}^{3}} = 0,

(7)

where $X_{bd 2}$ is the beam envelope of the $D_{2}^{+}$ species in the x-direction, $ε_{x d 1}$ is the beam emittance of the $D_{1}^{+}$ species in the x-direction, and $ε_{x d 1}^{2} = 4 〈 x^{2} 〉 〈 {x^{'}}^{2} 〉 - 4 {〈 x x^{'} 〉}^{2}$ .

Similarly, for $r_{1} > r_{2}$ , equation (5) can be written as

〈 x x^{″} 〉 + \frac{ϕ^{'}}{2 ϕ} 〈 x x^{'} 〉 + \frac{ϕ^{″}}{4 ϕ} 〈 x^{2} 〉 - \frac{I_{1}}{4 ϕ π ε_{0} β_{1} c} \frac{1}{r_{1}^{2}} 〈 x^{2} 〉 - \frac{I_{2}}{4 ϕ π ε_{0} β_{2} c} 〈 \frac{x^{2}}{r^{2}} 〉 = 0,

(8)

where 〈 \frac{x^{2}}{r^{2}} 〉 = \frac{\int_{r_{2}}^{r_{1}} \int_{0}^{2 π} \frac{x^{2}}{r^{2}} ρ_{1} (r, s) r d r d θ}{\int_{0}^{r_{1}} \int_{0}^{2 π} ρ_{1} (r, s) r d r d θ} = \frac{1}{2} (1 - \frac{r_{2}^{2}}{r_{1}^{2}}) = \frac{1}{2} (1 - \frac{X_{bd 2}^{2}}{X_{b}^{2}})

(9)

Thus, the envelope for the $D_{1}^{+}$ species of the multi-species beam can be described as

{X^{″}}_{b} + \frac{ϕ^{'}}{2 ϕ} {X^{'}}_{b} + \frac{ϕ^{″}}{4 ϕ} X_{b} - \frac{I_{1}}{4 ϕ π ε_{0} β_{1} c} \frac{1}{X_{b}} - \frac{I_{2}}{4 ϕ π ε_{0} β_{2} c} \frac{1}{2} (1 - \frac{X_{bd 2}^{2}}{X_{b}^{2}}) - \frac{ε_{x d 1}^{2}}{X_{b}^{3}} = 0

(10)

From the above, the x-direction envelope for the $D_{1}^{+}$ species of the multi-species beam consisting of $D_{1}^{+}$ and $D_{2}^{+}$ can be expressed as

{_{{X^{″}}_{b} + \frac{ϕ^{'}}{2 ϕ} {X^{'}}_{b} + \frac{ϕ^{″}}{4 ϕ} X_{b} - \frac{I_{1}}{4 ϕ π ε_{0} β_{1} c} \frac{1}{X_{b}} - \frac{I_{2}}{4 ϕ π ε_{0} β_{2} c} \frac{1}{2} (1 - \frac{X_{bd 2}^{2}}{X_{b}^{2}}) - \frac{ε_{x d 1}^{2}}{X_{b}^{3}} = 0 \begin{matrix} (r_{1} > r_{2}) \end{matrix}}^{{X^{″}}_{b} + \frac{ϕ^{'}}{2 ϕ} {X^{'}}_{b} + \frac{ϕ^{″}}{4 ϕ} X_{b} - \frac{I_{1}}{4 ϕ π ε_{0} β_{1} c} \frac{1}{X_{b}} - \frac{I_{2}}{4 ϕ π ε_{0} β_{2} c} \frac{1}{X_{bd 2}^{2}} X_{b} - \frac{ε_{x d 1}^{2}}{X_{b}^{3}} = 0 \begin{matrix} (r_{1} \leq r_{2}) \end{matrix}}

(11)

In the y-direction, the envelope equations have a similar form.

For the multi-species beam containing $D_{1}^{+}$ , $D_{2}^{+}$ , and $D_{3}^{+}$ , equation (3) can be written as

m_{1} \ddot{x} = q (E_{x}^{appl} + E_{D_{1} x}^{self} + E_{D_{2} x}^{self} + E_{D_{3} x}^{self}),

(12)

where $E_{D_{3} x}^{self}$ is the space charge field of the $D_{3}^{+}$ beam along the x-direction.

Then, the equations of motion of $D_{1}^{+}$ in x-direction can be described as

(13)

Using the same method as that used for the derivation of Eq.(11), the x-direction envelope for the $D_{1}^{+}$ species of the multi-species beam consisting of $D_{1}^{+}$ , $D_{2}^{+}$ , and $D_{3}^{+}$ can be expressed as

(14)

In the y-direction, the envelope equations for the particles of the three species have a similar form.

III. SIMULATION

A. Beam extraction system of the duoplasmatron source

1. Current operational beam extraction system

The current operational duoplasmatron ion source has been used for decades, and it is very reliable. The particle distributions in the phase space obtained in previous experiments^[12] are shown in Fig.2. For the LEBT simulation, the following parameters of the old beam extraction system are assumed^[5]: the normalized beam emittance is $ε_{x} = ε_{y} = 0.71$ mm·mrad, and the Twiss parameters are $α_{x} = α_{y} =$ −4 and $β_{x} = β_{y} =$ 80 cm/rad, respectively. The fraction of $D_{1}^{+}$ in the multi-species beam is approximately 75%, and $D_{2}^{+}$ and $D_{3}^{+}$ comprise approximately 18% and approximately 7% of the total beam, respectively.

Fig. 2

(Color online)Particle distributions in the phase space for the old beam extraction system^[12]

2. New designed beam extraction system

As shown in figure 3, a new beam extraction system of the duoplasmatron source is designed by the PIC simulation method. Compared to the old beam extraction system, the size of the extraction electrode, the dip angle of the anode, and the distance between the anode and the extraction electrode are changed. In addition, two small permanent magnet lenses with a magnetic strength of 1000 Gs are added. The simulated particle distributions in the phase space are shown in Fig. 4, where the normalized beam emittance is $ε_{x} = ε_{y} = 2.5$ mm·mrad, and the Twiss parameters are $α_{x} = α_{y} =$ −3.5 and $β_{x} = β_{y} = 57$ cm/rad, respectively.

Fig. 3

(Color online)New designed beam extraction system: (a) axial magnetic field; (b) configuration of the system, where the red lines represent the trajectories of the 50 mA

D_{1}^{+}

beam, and the green lines represent equipotential surfaces

Fig. 4

(Color online)Particle distributions in the phase space for the new designed beam extraction system

B. Transport characteristics of the multi-species beam

1. Effects of the unwanted beam composition on the primary beam transport

It is clear that using a pure beam is beneficial for improving the beam transport efficiency and decreasing the beam power. Therefore, various beam separation elements, such as dipole-bending magnets, Wien filters, or apertures combined with solenoids, are usually utilized to obtain a pure beam in the LEBT^[3]. However, using such beam separation elements will not only have the disadvantage of significantly increasing the length of the LEBT but is also not conducive to industrial application marketing due to the increased cost. In our scheme, the multi-species beam will be injected into the electrostatic LEBT directly.

In order to understand the influence of the unwanted beam on the primary beam, the evolution of the $D^{+}$ beam envelope for various ratios of species in the multi-species beam is simulated using the PIC code TRACK^[14]. The simulation results are illustrated in figure 5. They show that the envelope radius of the $D_{1}^{+}$ beam increases with the increasing unwanted beam intensity. The four envelopes corresponding to Beam (40 mA $D_{1}^{+}$ ), Beam (30 mA $D_{1}^{+}$ , 10 mA $D_{2}^{+}$ ), Beam (30 mA $D_{1}^{+}$ , 10 mA $D_{3}^{+}$ ), and Beam (30 mA $D_{1}^{+}$ , 5 mA $D_{2}^{+}$ , 5 mA $D_{3}^{+}$ ) almost fully overlap. The two envelopes corresponding to Beam (30 mA $D_{1}^{+}$ , 5 mA $D_{2}^{+}$ ) and Beam (30 mA $D_{1}^{+}$ , 5 mA $D_{3}^{+}$ ) almost fully overlap as well. This illustrates that the total intensity of the unwanted beams has a significant influence on the primary beam divergence. However, the influence of each species on the primary beam envelope is not obvious.

Fig. 5

(Color online)Evolution of the

D_{1}^{+}

beam envelopes of the multi-species beam with different ratios of the species

2. Effects of the injection beam parameters on the multi-species beam transport

Matching the extracted beam to the subsequent electrostatic LEBT beam is important for high intensity beam transport^[15]. Therefore, two different sets of the injected beam parameters from the old beam extraction system and the new designed beam extraction system, as described above, are used to simulate the effects of these parameters on the transport of the multi-species beam by the PIC code TRACK. The results are described below.

a. Transport efficiency

For the LEBT simulation, we use two different sets of the initial beam parameters obtained from the new and old beam extraction systems. For the new beam extraction system, the beam envelope radius is $X_{D 1} = X_{D 2} = X_{D 3} = 1.7$ cm, and the beam divergence angle is ${X^{'}}_{D 1} = {X^{'}}_{D 2} = {X^{'}}_{D 3} = 110$ mrad. For the old beam extraction system, the beam envelope radius is $X_{D 1} = X_{D 2} = X_{D 3} = 1.1$ cm, and the beam divergence angle is ${X^{'}}_{D 1} = {X^{'}}_{D 2} = {X^{'}}_{D 3} = 55$ mrad. The intensity of each species and the simulation results are shown in figure 6. It can be seen that using the initial beam parameters obtained from the new beam extraction system reduces the beam divergence for the beam transport in the electrostatic LEBT.

Fig. 6

(Color online) Evolution of the

D_{1}^{+}

beam envelopes of the multi-species beam with different ratios of the species for different injected beam parameters

The beam transport efficiency of each species in the electrostatic LEBT is simulated by the PIC code. The intensities of the $D_{1}^{+}$ , $D_{2}^{+}$ , and $D_{3}^{+}$ beams are 50 mA, 10 mA, and 10 mA, respectively. The beam envelope radius and divergence angle are the same with those described in the previous paragraph. The simulation results are shown in figure 7. When the initial beam parameters obtained from the new extraction system are used for the simulation, the transport efficiency of the $D_{1}^{+}$ , $D_{2}^{+}$ , and $D_{3}^{+}$ beams are 100%, 94%, and 39%, respectively. For the old beam extraction system, the corresponding values are 100%, 58%, and 52%. This implies that in the new beam extraction system, the $D_{3}^{+}$ beam will quickly hit the LEBT pipe and lose the majority of its particles.

Fig. 7

(Color online) Transport efficiency of each species beam for the new and old beam extraction systems

b. Emittance distortion and particle distributions

In the electrostatic LEBT, due to the presence of high electrostatic fields, the high intensity beam is fully un-neutralized. The space charge effect of the un-neutralized beam will cause a significant increase in emittance and strong beam filamentation. In addition, the nonlinear space charge field will cause emittance distortion^[16]. At the same time, the beam parameters at the output of the electrostatic LEBT are important for matching the beam to the subsequent beam transport elements. Therefore, we simulate the beam emittance and particle distributions. In the simulation, the fraction of each species beam is set based on the current operational duoplasmatron source, and the intensities of the $D_{1}^{+}$ , $D_{2}^{+}$ , and $D_{3}^{+}$ beams are set to 30 mA, 7 mA, and 3 mA, respectively. For the new beam extraction system, the beam envelope radius is $X_{D 1} = X_{D 2} = X_{D 3} = 1.7$ cm, and the beam divergence angle is ${X^{'}}_{D 1} = {X^{'}}_{D 2} = {X^{'}}_{D 3} = 110$ mrad. For the old beam extraction system, the beam envelope radius is $X_{D 1} = X_{D 2} = X_{D 3} = 1.1$ cm, and the beam divergence angle is ${X^{'}}_{D 1} = {X^{'}}_{D 2} = {X^{'}}_{D 3} = 55$ mrad. A water-bag distribution is used. The simulation results are shown in figures 8, 9 and 10.

Fig. 8

(Color online) Evolution of the envelopes of the multi-species beam for different injected beam parameters

Fig. 9

(Color online) Beam emittance at the output of the electrostatic LEBT;

D_{1}^{+}

D_{2}^{+}

, and

D_{3}^{+}

are represented by yellow, orange, and blue points, respectively; particle x–x’ (a) and particle y–y’ (b) with the injected beam parameters for the old beam extraction system; particle x–x’ (c) and particle y–y’ (d) with the injected beam parameters for the new designed beam extraction system

Fig. 10

(Color online) Particle distributions in real space at the output of the electrostatic LEBT;

D_{1}^{+}

D_{2}^{+}

, and

D_{3}^{+}

are represented by yellow, orange, and blue points, respectively; of the two graphs correspond to the old (a) and new designed (b) beam extraction systems, respectively

Figure 8 shows the simulated envelopes of the multi-species beam with different injected beam parameters, confirming that the injected beam parameters have a significant influence on the transport of the multi-species beam, and therefore it is necessary to match the beam to the electrostatic LEBT by optimizing the injected beam parameters.

Figure 9 shows the simulated beam emittance at the output of the electrostatic LEBT for different injected beam parameters. It can be seen that, for the new designed beam extraction system (Figs. 9(a) and 9(b)), the emittance distortion of the primary $D_{1}^{+}$ beam at the output of the electrostatic LEBT is more severe than that for the old beam extraction system (Figs. 9(c) and 9(d)); however, the particle distributions of the primary beam are more focused, as shown in figure 10. For the new designed beam extraction system, the emittance of the primary beam at the output is approximately 10.1 mm·mrad, ~4 times greater than that at the input, which is approximately 2.5 mm·mrad. However, for the old beam extraction system, the emittance of the primary beam at the output is approximately 9.1 mm·mrad, ~12 times greater than that at the input, which is approximately 0.7 mm·mrad.

From Fig. 10, it can be seen that hollow beams are formed for all species ( $D_{1}^{+}$ , $D_{2}^{+}$ , and $D_{3}^{+}$ ) when the beam is extracted from the old beam extraction system. From Fig. 10(a), the hollow structure of the $D_{1}^{+}$ beam is obvious, and exhibits a spot size of approximately 5 cm. However, in Fig.10(b), this effect is reduced significantly, and the spot size of the beam for the new extraction system is approximately 3 cm, which is obtained considering 98% of all $D_{1}^{+}$ particles and neglecting the distorted portion. It is found that the appropriate injected beam parameters play an important role in both the beam transport and the beam spot size in the electrostatic LEBT. This can be helpful for developing compact ion implanters and intense neutron generators.

IV. CONCLUSION

The envelope equations of the multi-species deuterium beam in an electrostatic LEBT are derived considering the space charge effects caused by the particles of all species. The evolution of the envelopes of the multi-species beam is simulated in the presence of different ratios of the components $D_{1}^{+}$ , $D_{2}^{+}$ , and $D_{3}^{+}$ , and for two sets of the injected beam parameters corresponding to the old and new beam extraction systems of the duoplasmatron source. The simulation results show that the envelope radius of the $D_{1}^{+}$ beam increases with the increase of the unwanted beam intensity. The new injected beam parameters improve the transport of the primary beam in the electrostatic LEBT. However, the corresponding emittance is found to be distorted at the output of the electrostatic LEBT. The results indicate that, with the appropriate injected beam parameters, an electrostatic LEBT is a promising choice for developing compact ion implanters and intense neutron generators.