Comparison between Theoretical and Experimental Radial Electron Temperature Profiles in a Low Density Weakly Ionized Plasma

Show more

1. Introduction

Plasma is the fourth state of matter, considered to be a quasi-neutral medium. However, when a diagnostic such as probe (a small metallic electrode) is inserted into a weakly ionized plasma, a very thin sheath is formed around the conducting surface of the probe due to the redistribution of charges. The amplitude and the polarity of the probe potential control the motion of electrons and ions near the probe. When the probe potential is sufficiently negative, only the ions can reach the probe surface. The probe current is thus equal to the ion random current (I_{ri}) [1] [2] [3] .

Ambipolar diffusion is the most important process taking place within a weakly ionized plasma and considered as a vital process for the distribution of the plasma parameters. Assuming a Maxwellian velocity distribution [4] , the plasma parameters especially electron temperature and electron density can be determined from the current-voltage characteristic curves of the probe [5] . Diffusion phenomena occur in a plasma when a spatial gradient of the charged species is present [6] . Furthermore, it is caused by the differently charged species having different diffusivities, hence there is a loss of neutrality [7] due to plasma rapidly diffusion for some of charged species more than the others. Furthermore, a minor loss of neutrality, however, induces an ambipolar electric field which, if the Debye length is sufficiently small, slows down the fast-diffusing species and speeds up the slow-diffusing species in such a way that the plasma remains quasi-neutral [8] .

The experimental and numerical radial distribution of electron temperature and density from the axis of the tube up to the tube wall have been amply investigated [9] [10] , but not the theoretical derivation in terms of the effect of Schottky condition and ambipolar diffusion due to the recurrence relation of the Bessel condition.

In the present work, an experimental study of the radial dependence of the electron temperature in a low-density plasma using a single probe as diagnostic technique was performed. Theoretically, it was proven that the ambipolar diffusion is an important factor in the weakly-ionized gas. Furthermore, a theoretical physical model for the radial distribution of the temperature was derived. The theoretical considerations are compatible with the experimental data of the temperatures at different pressures in a low-pressure glow discharge of a DC (cold cathode) magnetron sputtering unit.

2. Experimental Study

2.1. Experimental Set-Up

Figure 1 shows the experimental setup of DC (cold cathode) magnetron sputtering unit to generate a glow discharge in a glass tube between two circular, parallel and movable metallic discs acting as electrodes. Two parallel electrodes made of aluminum, are enclosed in the discharge cell, one of the two electrodes is grounded represented cathode electrode and the other is movable to change the axial distance represented the anode electrode, where the axial distance (L) between tip of the probe and the cathode electrode, 5 cm and with 3 cm in radius (R) taking into account that we will deal with two conditions R > L and R < L. The discharge unit is evacuated using a rotary pump to a base pressure of 7 mtorr. A pressure gauge is connected to the discharge tube to measure the inside gas pressure. A stationary DC-glow discharge was generated between two both electrodes, for different parameters as shown in Table 1.

Figure 1 also shows the schematic diagram of the spherical single Langmuir probe circuit. The probe made of molybdenum wire (diameter 3.0-mm and length 0.5 mm) and the tip of a probe is inside the glow discharge plasma. The single probe (between the cathode fixed at the ground potential, and the anode) move axially to measure the axial potential distribution between the two electrodes. Also, it moves radially from the center to the edge of the cathode electrode, and a potential V_{P} is applied to the probe.

2.2. Experimental Study of Axial Potential and Electric Field Distribution Measurements

The potential distribution was measured at discharge currents of 10 mA and gas pressure of (0.2 - 1.2 torr) for He as shown in Figure 2, shows that the potential

Figure 1. Schematic diagram of the experimental setup and the single Langmuir probe circuit.

Figure 2. Axial potential distribution at constant current (I = 10 mA).

Table 1. Shows the values of operating parameters.

distribution can be divided into three regions. In region I (AB) (cathode fall), the potential increases sharply within a small discharge length. Whenever, gas breakdown takes place in the tube, a rapid growth in the rate of ionization, near the cathode is detected. Meanwhile the electrons transfer in the electric field much more faster than positive ions (due to their masses), electrons are swept rapidly towards the anode leaving a dense positive space charge near the cathode. Thus, the electric field is distorted and most of the applied potential is dropped across a narrow space in front of the cathode. Region II (BC) (negative glow) indicates that the potential decreases slightly and hence the electric field will be weak, since this region contains many free electrons. In region III (CD) (positive column), the potential distribution is nearly constant and linear since the positive and negative carriers densities are closely, equal. The positive column can be extended to any length to fill the remaining space between the end of the negative glow and the anode.

Values of the electric field distribution are obtained by differentiating the measured potential distribution of Figure 2. Figure 3 and Figure 4 show the electric field distribution for He discharge. In the cathode fall region at edge and center, high electric field is observed which is decreased sharply away from the

Figure 3. I-V curves of the single probe for cathode fall region at different He pressures at the edge.

Figure 4. Measured I-V curves of the single probe for cathode fall region at different He pressures at the center.

cathode. This is related to the intense positive space charge which lies in front of the cathode. This acts as an accelerator for the electrons towards anode. Thus, the electrons emitted from the cathode are then accelerated away, until they reach the negative glow region where the electric field becomes weak (zero and sometimes it reaches a negative value). In this region, the gained kinetic energy of the electrons is dissipated in collisions with the atoms of the gas and thus secondary electrons would be produced. In the positive column region constant and linear electric field is needed to maintain discharge along the large length of the column which is required to carry the discharge current.

2.3. Experimental Study of Radial Dependence of the Electron Temperature

The region investigated in the present work is the abnormal region. The radial dependence of the electron temperature in low-density plasma at edge and center of the cathode electrode has been studied for different He pressures.

The plasma parameters like the electron temperature and electron density can be determined from the current-voltage (I-V) characteristic curve of the single Langmuir probe, based on the theory and the fundamental technique discussed in detail in many articles [10] [11] .

Figure 3 and Figure 4 show I-v characteristic curves of the single probe at different He pressures at edge and center of the cathode electrode for the abnormal glow region [12] . The probe was moved to investigate the radial distribution from the center to the edge of the cathode, in the direction perpendicular to the direction of the electric field lines. Figure 5 and Figure 6 show the radial temperatures and densities distributions at different pressures from center to the edge. T_{e} is decreased and n_{e} increased due to the general trend that values of T_{e} and n_{e}_{ }are inversely proportional [13] .

Increasing the helium working pressures from 0.2 to 1.2 torr, the temperature T_{e} decreases at the center from 9.7 eV to 5.5 eV and at the edge from 4.1 to 1.15 eV. Values of densities n_{e} increased at center n_{e} from 0.6 × 10^{9} cm^{−3} to 1.3 × 10^{9} cm^{−3} and at the edge from 2.5 to 3.95 × 10^{9} cm^{−3}.

3. Theoretical Consideration of the Radial Dependence of the Electron Temperature

Conductivity at the tube axis, and charge density of plasma according to the following expression [14] [15]

${D}_{s}={D}_{a}\left[1-{\mu}_{e}\frac{\rho}{\sigma}\right]$ (1)

where D_{a} is the ambipolar diffusion coefficient as follows:

The flow of ions and electrons are the same, hence

Figure 5. The Experimental radial temperature distribution of abnormal region for helium at different pressures.

Figure 6. Measured radial density distribution of abnormal region for helium at different pressures.

$\Gamma ={\Gamma}_{i}={\Gamma}_{e}$

(Congruence approximation). If there is no external electric field, the fluxes of ions and electrons (drift-diffusion model) can be expressed as:

${\Gamma}_{i}=-{D}_{i}\nabla {n}_{i}+{\mu}_{i}{n}_{i}{E}_{D}$

${\Gamma}_{e}=-{D}_{e}\nabla {n}_{e}+{\mu}_{e}{n}_{e}{E}_{D}$

where E_{D} follows Poisson’s equation. If we multiply the first equation by
${\mu}_{e}{n}_{e}$ and the second one by
${\mu}_{i}{n}_{i}$ and then subtract them, we get:

${\mu}_{i}{n}_{i}{\Gamma}_{e}-{\mu}_{e}{n}_{e}{\Gamma}_{i}={\mu}_{e}{n}_{e}{D}_{i}\nabla {n}_{i}-{\mu}_{e}{n}_{e}{\mu}_{i}{n}_{i}{E}_{D}-{\mu}_{i}{n}_{i}{D}_{e}\nabla {n}_{e}+{\mu}_{i}{n}_{i}{\mu}_{e}{n}_{e}{E}_{D}$

Since
$\Gamma ={\Gamma}_{i}={\Gamma}_{e}$ , the terms with E_{D} cancel each other out and using n_{e} = n_{i}, the expression for the flow of particles becomes:

$\Gamma =\frac{-\left({D}_{e}{\mu}_{i}-{D}_{i}{\mu}_{e}\right)}{{\mu}_{i}-{\mu}_{e}}\nabla {n}_{i}$

With the ambipolar diffusion coefficient being

${D}_{a}=\frac{{D}_{e}{\mu}_{+}-{D}_{+}{\mu}_{e}}{{\mu}_{+}-{\mu}_{e}}$ (2)

and $\rho $ is the charge density in the plasma given by:

$\rho =e\left[{n}_{+}-{n}_{e}\right]$ (3)

and $\sigma $ is the conductivity at the tube axis given by:

$\sigma =e\left[{\mu}_{+}{n}_{+}-{\mu}_{e}{n}_{e}\right]$ (4)

Considering that diffusion coefficient D given by:

$D=\frac{kT}{m\nu}\text{\hspace{0.17em}}\text{m}/{\text{s}}^{\text{2}}$ (5)

where ${\mu}_{\epsilon ,+}$ are the mobility of electrons (e) and ions (+), kT (ev) is the particle temperature, $\nu $ is the collision frequency between electrons and neutral atoms in Hz. By substituting with Equations (2, 3 and 4) into Equation (1), then we get

${D}_{s}=\frac{{D}_{e}{\mu}_{+}-{D}_{+}{\mu}_{e}}{{\mu}_{+}-{\mu}_{e}}\left[1-{\mu}_{e}\frac{{n}_{+}-{n}_{e}}{{\mu}_{+}{n}_{+}-{\mu}_{e}{n}_{e}}\right]$ (6)

But by substituting by the Einstein relation [16] [17] :

$\mu =\frac{\lceil q\rceil D}{kT}$ (7)

into Equation (6) by the value of D, then we get

${D}_{s}=\frac{{\mu}_{e}{\mu}_{+}K{T}_{e}-{\mu}_{e}{\mu}_{+}K{T}_{+}}{q\left({\mu}_{+}-{\mu}_{e}\right)}\left[1-{\mu}_{e}\frac{{n}_{+}-{n}_{e}}{{\mu}_{+}{n}_{+}-{\mu}_{e}{n}_{e}}\right]$ (8)

or

${D}_{s}=\frac{{\mu}_{e}{\mu}_{+}\left(K{T}_{e}-K{T}_{+}\right)}{e\left({\mu}_{+}-{\mu}_{e}\right)}\left[\frac{{\mu}_{+}{n}_{+}-{\mu}_{e}{n}_{e}-{\mu}_{e}{n}_{+}+{\mu}_{e}{n}_{e}}{{\mu}_{+}{n}_{+}-{\mu}_{e}{n}_{e}}\right]$ (9)

Neglecting $k{T}_{+}$ , where $k{T}_{e}\gg k{T}_{+}$ finally results in

Then

${D}_{s}=\frac{{n}_{+}{\mu}_{e}{\mu}_{+}\left(k{T}_{e}\right)}{e\left({\mu}_{+}-{\mu}_{e}\right)}\left[\frac{{\mu}_{+}-{\mu}_{e}}{{\mu}_{+}{n}_{+}-{\mu}_{e}{n}_{e}}\right]$ (10)

${D}_{s}=\frac{{n}_{+}{\mu}_{e}{\mu}_{+}\left(k{T}_{e}\right)}{e}\left[\frac{1}{{\mu}_{+}{n}_{+}-{\mu}_{e}{n}_{e}}\right]$

Then

${D}_{s}=\frac{k{T}_{e}}{e}\left[\frac{{\mu}_{e}{\mu}_{+}}{{\mu}_{+}-{\mu}_{e}\frac{{n}_{e}}{{n}_{+}}}\right]$ (11)

or

${D}_{s}=\frac{k{T}_{e}}{e}\left[\frac{{\mu}_{e}}{1-\frac{{\mu}_{e}}{{\mu}_{+}}\frac{{n}_{e}}{{n}_{+}}}\right]$ (12)

Substituting ${\mu}_{e}=\frac{e}{{m}_{e}{\nu}_{e}}$ and ${\mu}_{i}=\frac{e}{{m}_{i}{\nu}_{i}}$ from Equations (5) and (7) into Equation (12) results in

${D}_{S}=\frac{k{T}_{e}}{e}\left[\frac{\left(\frac{e}{{m}_{e}{\nu}_{e}}\right)}{1-\frac{\left(\frac{e}{{m}_{e}{\nu}_{e}}\right){n}_{e}}{\left(\frac{e}{{m}_{i}{\nu}_{i}}\right){n}_{+}}}\right]$ (13)

or

${D}_{S}=k{T}_{e}\left[\frac{\frac{1}{{m}_{e}{\nu}_{e}}}{1-\frac{{m}_{i}{\nu}_{i}{n}_{e}}{{m}_{e}{\nu}_{e}{n}_{+}}}\right]$ (14)

or

${D}_{S}=k{T}_{e}\left[\frac{{n}_{+}}{{m}_{e}{\nu}_{e}{n}_{+}-{m}_{i}{\nu}_{i}{n}_{e}}\right]$ (15)

$\frac{{D}_{S}}{{\nu}_{i}}=\frac{k{T}_{e}}{{m}_{e}{\nu}_{e}{\nu}_{i}-\frac{{m}_{i}{\nu}_{i}^{2}{n}_{e}}{{n}_{+}}}$ (16)

Experimentally, due to the diffusion process during plasma formation an interesting process occur in the plasma formation stage of the basil discharge. Theoretically, basil discharge process related with the recurrence relation for Bessel condition as shown in appendix A, where from the recurrence relation for the Bessel condition, the Schottky condition can be derived and proofed [18] [19] stating that:

3.1. In the Case of L (Axial Distance) > R (Radial Distance)

$\frac{{D}_{S}}{{\nu}_{i}}={\left(\frac{R}{2.405}\right)}^{2}$ (17)

A proof of the Schottky condition will be discussed briefly by two methods in Appendix A, where (R cm) represents the radial distance from the center to the edge through the cathode electrode measured by the Langmuir probe moved radially, and by substituting Equation (17) into (16), then

${\left(\frac{R}{2.4}\right)}^{2}=\frac{k{T}_{e}}{{m}_{e}{\nu}_{e}{\nu}_{i}-\frac{{m}_{i}{\nu}_{i}^{2}{n}_{e}}{{n}_{+}}}$ (18)

Substituting the following equation into Equation (18)

${\nu}_{e}={N}_{n}\langle {Q}_{e\text{-}n}{\upsilon}_{e\text{-}n}\rangle $ (19a)

${\nu}_{e}=\left(3.55\times {10}^{16}P\right)\langle {Q}_{e\text{-}n}\rangle {\left(\frac{2k{T}_{e}}{{m}_{e}}\right)}^{1/2}$ (19b)

${\nu}_{i}=\left(3.55\times {10}^{16}P\right)\langle {Q}_{i\text{-}n}\rangle {\left(\frac{2k{T}_{i}}{{m}_{i}}\right)}^{1/2}$ (19c)

where P in torr, ${Q}_{e\text{-}n}$ and ${Q}_{i\text{-}n}$ are the cross sections of electron-neutral and ion-neutral collisions, respectively results in

${\left(\frac{R}{2.4}\right)}^{2}=\frac{k{T}_{e}}{\left[{m}_{e}{\left(\frac{2k{T}_{e}}{{m}_{e}}\right)}^{1/2}{\left(\frac{2k{T}_{i}}{{m}_{i}}\right)}^{1/2}{\left(3.55\times {10}^{16}P\right)}^{2}\langle {Q}_{i}\rangle \langle {Q}_{e}\rangle \right]-\left[\frac{{n}_{e}}{{n}_{+}}{m}_{i}{\left[{\left(3.55\times {10}^{16}P\right)}^{2}\langle {Q}_{i}\rangle \right]}^{2}\left(\frac{2k{T}_{i}}{{m}_{i}}\right)\right]}$ (20)

then

${\left(\frac{RP}{2.4}\right)}^{2}=\frac{1}{\left[{\left(\frac{4{m}_{e}k{T}_{i}}{{m}_{i}k{T}_{e}}\right)}^{1/2}{\left(3.55\times {10}^{16}\right)}^{2}\langle {Q}_{i}\rangle \langle {Q}_{e}\rangle \right]-\left[\frac{{n}_{e}}{{n}_{+}}{\left(3.55\times {10}^{16}\langle {Q}_{i}\rangle \right)}^{2}\left(\frac{2k{T}_{i}}{{m}_{i}k{T}_{e}}\right)\right]}$ (21)

But from [20] [21] $\frac{{\mu}_{e}}{{\mu}_{i}}\cong 7.64{\left({m}_{i}\right)}^{0.5}$

Or

${m}_{i}\cong {\left(\frac{1}{7.64}\frac{{\mu}_{e}}{{\mu}_{i}}\right)}^{2}$ (22a)

Then substituting (7) into (22a) gives

$\frac{1}{{m}_{i}}\cong {\left(7.64\frac{k{T}_{i}}{k{T}_{e}}\right)}^{2}$ (22b)

Moreover, substituting (22b), $\frac{{m}_{i}}{{m}_{e}}\cong 2000$ and $\frac{{\nu}_{i}}{{\nu}_{e}}\cong {10}^{-2}$ into (21) gives

${\left(\frac{RP}{2.4}\right)}^{2}=\frac{1}{\left[\left(7.64\frac{k{T}_{i}}{k{T}_{e}}\right){\left(\frac{4{m}_{e}k{T}_{i}}{k{T}_{e}}\right)}^{1/2}{\left(3.55\times {10}^{16}\right)}^{2}\langle {Q}_{i}\rangle \langle {Q}_{e}\rangle \right]-\left[\frac{{n}_{e}}{{n}_{+}}{\left(3.55\times {10}^{16}\langle {Q}_{i}\rangle \right)}^{2}{\left(7.64\frac{k{T}_{i}}{k{T}_{e}}\right)}^{2}\left(\frac{2k{T}_{i}}{k{T}_{e}}\right)\right]}$ (23)

${\left(RP\right)}^{2}=\frac{{2.4}^{2}}{\left[\left(7.64\frac{k{T}_{i}}{k{T}_{e}}\right){\left(\frac{4{m}_{e}k{T}_{i}}{k{T}_{e}}\right)}^{1/2}{\left(3.55\times {10}^{16}\right)}^{2}\langle {Q}_{i}\rangle \langle {Q}_{e}\rangle \right]-\left[\frac{{n}_{e}}{{n}_{+}}{\left(3.55\times {10}^{16}\langle {Q}_{i}\rangle \right)}^{2}{\left(7.64\frac{k{T}_{i}}{k{T}_{e}}\right)}^{2}\left(\frac{2k{T}_{i}}{k{T}_{e}}\right)\right]}$ (24)

Substituting the values of the cross sections, masses of ions and electrons into (24) gives

${\left(RP\right)}^{2}=\frac{{2.4}^{2}}{{\left(3.55\right)}^{2}\left(7.64\right)\left(2\right)\left[\left[\left(\frac{k{T}_{i}}{k{T}_{e}}\right){\left(\frac{{m}_{e}k{T}_{i}}{k{T}_{e}}\right)}^{1/2}\right]-\left[\frac{{n}_{e}}{{n}_{+}}7.64{\left(\frac{k{T}_{i}}{k{T}_{e}}\right)}^{2}\left(\frac{k{T}_{i}}{k{T}_{e}}\right)\right]\right]}$ (25)

or

${\left(RP\right)}^{2}=\frac{{2.4}^{2}}{{\left(3.55\right)}^{2}\left(7.64\right)\left(2\right){\left(\frac{k{T}_{i}}{k{T}_{e}}\right)}^{3/2}\left[{m}_{e}^{1/2}-7.64\frac{{n}_{e}}{{n}_{+}}{\left(\frac{k{T}_{i}}{k{T}_{e}}\right)}^{3/2}\right]}$ (26)

Then finally

$RP=0.1729\frac{1}{{\left(\frac{k{T}_{i}}{k{T}_{e}}\right)}^{3/4}{\left[{m}_{e}^{1/2}-7.64\frac{{n}_{e}}{{n}_{+}}{\left(\frac{k{T}_{i}}{k{T}_{e}}\right)}^{3/2}\right]}^{1/2}}$ (27)

Using the data of the electron temperatures and densities shown in Figure 4 and Figure 5 respectively, into the Equation (27), taking into account that ${T}_{e}\gg {T}_{i}$ ( ${T}_{i}=0.1{T}_{e}$ ), $\frac{{m}_{e}}{{m}_{i}}\cong 2000$ [22] and for the spherical probe, the ion density is given by [23] :

${n}_{+}=\frac{{I}_{+}}{0.6}{\left(\frac{k{T}_{e}}{{m}_{i}}\right)}^{1/2}{A}_{{}_{P}}$ (28)

Knowing T_{e}, (
${A}_{p}=4\text{\pi}{r}_{p}^{2}$ ) area of spherical probe) and calculating positive ion current I_{+} from the I-V characteristic curve of the spherical single probe [24] , n_{+ }can be determined.

Figure 7 shows the radial distribution of the electron temperature (RDFT) theoretically using Equation (27), in the abnormal cathode fall region from the center to the edge of the electrode as a function of He as an inert gas. The RDFT from the theoretical method presented above accurately shows a dramatic radial reduction of the electron temperature for any region (cathode fall or negative glow or positive column) and for any applied pressure. Furthermore, the decrements of T_{e} at the edge of the electrode are more pronounced at the center due to edge effect [25] .

3.2. In the Case of L (Axial Distance) < R (Radial Distance)

$\frac{{D}_{S}}{{\nu}_{i}}={\left(\frac{R}{2.405}\right)}^{2}+{\left(\frac{L}{\text{\pi}}\right)}^{2}$ (29)

and by substituting Equation (29) into (16), then

${\left(\frac{R}{2.4}\right)}^{2}+{\left(\frac{L}{\text{\pi}}\right)}^{2}=\frac{k{T}_{e}}{{m}_{e}{\nu}_{e}{\nu}_{i}-\frac{{m}_{i}{\nu}_{i}^{2}{n}_{e}}{{n}_{+}}}$ (30)

Substituting by Equations (19a, b and c) into Equation (30), respectively results in

Figure 7. Theoretical radial distribution of the electron temperature (RDFT).

${\left(\frac{R}{2.4}\right)}^{2}+{\left(\frac{L}{\text{\pi}}\right)}^{2}=\frac{k{T}_{e}}{\left[{m}_{e}{\left(\frac{2k{T}_{e}}{{m}_{e}}\right)}^{1/2}{\left(\frac{2k{T}_{i}}{{m}_{i}}\right)}^{1/2}{\left(3.55\times {10}^{16}P\right)}^{2}\langle {Q}_{i}\rangle \langle {Q}_{e}\rangle \right]-\left[\frac{{n}_{e}}{{n}_{+}}{m}_{i}{\left[{\left(3.55\times {10}^{16}P\right)}^{2}\langle {Q}_{i}\rangle \right]}^{2}\left(\frac{2k{T}_{i}}{{m}_{i}}\right)\right]}$ (31)

then

${\left(\frac{R}{2.4}\right)}^{2}+{\left(\frac{L}{\text{\pi}}\right)}^{2}=\frac{1}{\left[{\left(\frac{4{m}_{e}k{T}_{i}}{{m}_{i}k{T}_{e}}\right)}^{1/2}{\left(3.55\times {10}^{16}\right)}^{2}\langle {Q}_{i}\rangle \langle {Q}_{e}\rangle \right]-\left[\frac{{n}_{e}}{{n}_{+}}{\left(3.55\times {10}^{16}\langle {Q}_{i}\rangle \right)}^{2}\left(\frac{2k{T}_{i}}{{m}_{i}k{T}_{e}}\right)\right]}$ (32)

Substitution from (22a and b) into (32) gives

${\left(\frac{RP}{2.4}\right)}^{2}+{\left(\frac{L}{\text{\pi}}\right)}^{2}=\frac{1}{\left[\left(7.64\frac{k{T}_{i}}{k{T}_{e}}\right){\left(\frac{4{m}_{e}k{T}_{i}}{k{T}_{e}}\right)}^{1/2}{\left(3.55\times {10}^{16}\right)}^{2}\langle {Q}_{i}\rangle \langle {Q}_{e}\rangle \right]-\left[\frac{{n}_{e}}{{n}_{+}}{\left(3.55\times {10}^{16}\langle {Q}_{i}\rangle \right)}^{2}{\left(7.64\frac{k{T}_{i}}{k{T}_{e}}\right)}^{2}\left(\frac{2k{T}_{i}}{k{T}_{e}}\right)\right]}$ (33)

${\left(RP\right)}^{2}=\frac{{2.4}^{2}}{\left[\left(7.64\frac{k{T}_{i}}{k{T}_{e}}\right){\left(\frac{4{m}_{e}k{T}_{i}}{k{T}_{e}}\right)}^{1/2}{\left(3.55\times {10}^{16}\right)}^{2}\langle {Q}_{i}\rangle \langle {Q}_{e}\rangle \right]-\left[\frac{{n}_{e}}{{n}_{+}}{\left(3.55\times {10}^{16}\langle {Q}_{i}\rangle \right)}^{2}{\left(7.64\frac{k{T}_{i}}{k{T}_{e}}\right)}^{2}\left(\frac{2k{T}_{i}}{k{T}_{e}}\right)\right]}-{2.4}^{2}{\left(\frac{L}{\text{\pi}}\right)}^{2}$ (34)

Substituting the values of the cross sections, masses of ions and electrons into (34) gives finally

$RP=\sqrt{\left(\frac{0.0299}{{\left(\frac{k{T}_{i}}{k{T}_{e}}\right)}^{3/2}\left[{m}_{e}^{1/2}-7.64\frac{{n}_{e}}{{n}_{+}}{\left(\frac{k{T}_{i}}{k{T}_{e}}\right)}^{3/2}\right]}-5.76{\left(\frac{L}{\text{\pi}}\right)}^{2}\right)}$ (35)

For 1.2 torr of applied pressure, Figure 8 shows comparison between the radial distribution of the electron temperature (RDFT) for L < R using Equation (35) and for L > R using Equation (27), in the abnormal cathode fall region from the center to the edge of the electrode as a function of He as an inert gas.

The reasons for the radial distribution of the electron temperatures for both cases.

L < R and L > R are the following:

a) The temperature decreases for increasing pressure may be due to the relation between the temperature and the pressure given by [26] :

$\left|{T}_{e}\right|=6.5\left|{P}^{-0.066}\right|$ (36)

b) As the pressure increase leads to a further increase in the breakdown voltage, a sharp increase in electron density and a higher electron-electron collision frequency, leads to the decrease of the electron temperatures [27] .

c) The present theoretical and experimental results are compatible and agree fairly well with the behavior of the data of von Engel [28] [29] [30] for the radial distribution of the electron temperature.

For L < R, RDFT accurately shows a dramatic decrease by 60% less than RDFT values for L > R. This may be attributed to:

Figure 8. Comparison between the radial distribution of the electron temperature (RDFT) for L < R and for L > R, at P = 1.2 torr.

The cathode fall, negative glow and positive column regions are compressed when L < R, i.e. the cathode fall region is compressed in thickness and henceforth higher potentials are expected. Thus, strong electric field is produced and therefore ions would accelerate, and more efficient sputtering processes take place. The produced plasma temperatures decrease, and densities increases; this means that the rate of plasma loss by diffusion decreased similarly as in the applied of magnetized DC plasma [23] , therefore the current and current density are increased.

4. Conclusion

Theoretically, using the Einstein relation of ambipolar diffusion, charge density of the plasma, conductivity at the axis of the tube, the Schottky condition, and the cross sections of ion-neutral and electron-neutral collisions, the decrease of the radial distribution function of the electron temperature with increasing product of the radial distance and gas pressures can be determined. The theoretical study agrees well with experimental data of electron and ion temperatures and gives the radial decrement dependence of the electron temperatures from the center to the edge of the electrode.

Acknowledgements

GVO acknowledges for the partial financial support from MEPhI and NRU MPEI in the framework of the Russian Academic Excellence Project.

Appendix: Proofing of the Schottky Condition

There are two classes of Schottky condition models using.

Appendix A: the recurrence relation for the Bessel condition

Appendix B: the boundary conditions

$\frac{{D}_{S}}{{\nu}_{i}}=\frac{K{T}_{e}}{\frac{{m}_{i}{\nu}_{i}^{2}{n}_{e}}{{n}_{+}}+{m}_{e}{\nu}_{e}{\nu}_{i}}$ (36)

Due to the diffusion process during plasma formation an interesting process occur in the plasma formation stage of the basil discharge. The Schottky condition stating that

$\frac{{D}_{S}}{{\nu}_{i}}={\left(\frac{R}{2.4}\right)}^{2}$ can be demonstrated as follows:

Appendix A: From the Recurrence Relation for Bessel Condition [31] [32] [33]

$2{{J}^{\prime}}_{n}\left(x\right)={J}_{n-1}\left(x\right)-{J}_{n+1}\left(x\right)$ (37)

and

$\frac{2n}{x}{J}_{n}\left(x\right)={J}_{n-1}\left(x\right)+{J}_{n+1}\left(x\right)$ (38)

Subtracting equation (37) from (38), gives

$\frac{2n}{x}{J}_{n}\left(x\right)-2{{J}^{\prime}}_{n}\left(x\right)={J}_{n-1}\left(x\right)+{J}_{n+1}\left(x\right)-{J}_{n-1}\left(x\right)+{J}_{n+1}\left(x\right)$ (39)

Or

$\frac{n}{x}{J}_{n}\left(x\right)-{{J}^{\prime}}_{n}\left(x\right)={J}_{n+1}(\; x\; )$

Or

$n{J}_{n}\left(x\right)-x{{J}^{\prime}}_{n}\left(x\right)=x{J}_{n+1}(\; x\; )$

This results in

$x{{J}^{\prime}}_{n}\left(x\right)=n{J}_{n}\left(x\right)-x{J}_{n+1}\left(x\right)$ (40)

Differentiating (40) with respect to x, gives:

$x{{J}^{\u2033}}_{n}\left(x\right)+{{J}^{\prime}}_{n}\left(x\right)=n{{J}^{\prime}}_{n}\left(x\right)-x{{J}^{\prime}}_{n+1}\left(x\right)-{J}_{n+1}\left(x\right)$ (41)

Multiplying Equation (40) in n/x gives:

$n{{J}^{\prime}}_{n}\left(x\right)=\frac{{n}^{2}}{x}{J}_{n}\left(x\right)-n{J}_{n+1}\left(x\right)$ (42)

Adding Equations (37) and (38) gives

$x{{J}^{\prime}}_{n}\left(x\right)=x{J}_{n-1}\left(x\right)-n{J}_{n}\left(x\right)$ (43)

Putting n = n + 1 in Equation (43) gives

$x{{J}^{\prime}}_{n+1}\left(x\right)=x{J}_{n}\left(x\right)-\left(n+1\right){J}_{n+1}\left(x\right)$ (44)

Substituting (39) and (40) into (38), results in

$x{{J}^{\u2033}}_{n}\left(x\right)+{{J}^{\prime}}_{n}\left(x\right)=\frac{{n}^{2}}{x}{J}_{n}\left(x\right)-x{J}_{n}(\; x\; )$

Then

${{J}^{\u2033}}_{n}\left(x\right)+\frac{1}{x}{{J}^{\prime}}_{n}\left(x\right)+\left(1-\frac{{n}^{2}}{{x}^{2}}\right){J}_{n}\left(x\right)=0$ (45)

Finally

${J}^{\u2033}+\frac{1}{x}{J}^{\prime}+\left(\lambda -\frac{{n}^{2}}{{x}^{2}}\right)J=0$ (46)

For n = 0

${J}^{\u2033}+\frac{1}{x}{J}^{\prime}+\left(\lambda \right)J=0$ (47)

For the Bessel condition in cylindrical coordinate

$x=Kr={\left(\frac{1}{{D}_{a}\tau}\right)}^{1/2}R={\left(\frac{\lambda}{{D}_{a}}\right)}^{1/2}R$

with the solution ${J}_{n}\left(x\right)=\sqrt{\frac{2}{\text{\pi}}}\mathrm{cos}\left(x-\frac{n\text{\pi}}{2}-\frac{\text{\pi}}{4}\right)$ , and ${J}_{0}=0$ for n = 0, gives

$0=\mathrm{cos}\left(x-\frac{\text{\pi}}{4}\right)$ with the solution $x=\frac{3\text{\pi}}{4}\approx 2.4$

$x=Kr=\frac{R}{{\left(D\tau \right)}^{1/2}}\approx 2.4$ or $\frac{R}{2.4}={\left(D\tau \right)}^{1/2}$

Then

$\frac{R}{2.4}={\left(\frac{{D}_{S}}{{\nu}_{i}}\right)}^{1/2}$

or

$\frac{{D}_{S}}{{\nu}_{i}}={\left(\frac{R}{2.4}\right)}^{2}$

Appendix B: From the Boundary Conditions [34] [35]

${n}_{e}{\upsilon}_{e}=-{D}_{S}\nabla {n}_{e}$

Then

${D}_{S}=-\frac{{n}_{e}{\upsilon}_{e}}{\nabla {n}_{e}}$ (48)

But

${\nu}_{i}=\frac{\nabla \cdot {n}_{e}{\upsilon}_{e}}{{n}_{e}}$

Then

$\frac{{D}_{S}}{{\nu}_{i}}=-{\left(\frac{{n}_{e}}{\nabla {n}_{e}}\right)}^{2}$ (49)

For boundary conditions $r=R$ , ${n}_{e}={n}_{e0}$ , and ${J}_{0}=2.4$ , then

$\nabla {n}_{e}=-\frac{2.4{n}_{e}}{R}$ (50)

Substituting from (50) into (49) then we get

$\frac{{D}_{S}}{{\nu}_{i}}={\left(\frac{R}{2.4}\right)}^{2}$

References

[1] Chen, F.F. (1965) Plasma Diagnostic Techniques. Academic Press, New York.

[2] Radmilovic-Radjenovic, M. and Radjenovic, B. (2006) Journal of Physics D: Applied Physics, 39, 3002.

https://doi.org/10.1088/0022-3727/39/14/019

[3] Kagan, Y. and Paskalev, K.K. (1975) Soviet Physics. Technical Physics, 19, 1604.

[4] Kagan, Y. and Perel, V.I. (1964) Soviet Physics Uspekhi, 6, 767.

https://doi.org/10.1070/PU1964v006n06ABEH003611

[5] Mott-Smith and Langmuir (1926) Physical Review, 28, 727.

https://doi.org/10.1103/PhysRev.28.727

[6] Parent, B., Macheret, S.O. and Shneider, M.N. (2011) Journal of Computational Physics, 230, 8010-8027.

https://doi.org/10.1016/j.jcp.2011.07.006

[7] Holt, E.H. and Haskell, R.E. (1965) Foundations of Plasma Dynamics. 2nd Edition, MacMillan Company, New-York.

[8] Huang, P.G., Shang, J.S. and Stanfield, S.A. (2011) AIAA Journal, 49, 119.

https://doi.org/10.2514/1.J050463

[9] Bisek, N.J., Boyd, I.D. and Poggie (2009) Journal of Spacecraft and Rockets, 46, 568.

https://doi.org/10.2514/1.39032

[10] Galaly, A.R. and Khedr, M.A. (2015) British Journal of Applied Science & Technology, 11, 1-14.

https://doi.org/10.9734/BJAST/2015/17663

[11] Galaly, A.R. (2016) British Journal of Applied Science & Technology, 13, 1.

[12] Galaly, A.R. and El Akshar, F.F. (2013) Physica Scripta, 88, Article ID: 065503.

https://doi.org/10.1088/0031-8949/88/06/065503

[13] Yasuda, H.K., Tao, W.H. and Prelas, M.A. (1996) Journal of Vacuum Science & Technology A: Vacuum, Surfaces, and Films, 14, 2113.

https://doi.org/10.1116/1.580089

[14] Nasser, E. (1971) Fundamentals of Gaseous Ionization and Plasma Electronics. Wiley-Interscience, New York.

[15] Raizer, Y.P. (1999) Gas Discharge Physics. Springer, New York.

[16] Von Engel, A.H. (1964) Ionized Gases. Oxford University Press, Oxford.

[17] Franklin, R.N. (1977) Plasma Phenomena in Gas Discharges. Oxford University Press, Oxford.

[18] Lieberman, M.A. and Lichtenberg, A.J. (1994) Principles of Plasma Discharges for Materials Processing. Wiley Interscience, New York.

[19] Braithwaite, N.St.J. (2000) Plasma Sources Science and Technology, 9, 517.

https://doi.org/10.1088/0963-0252/9/4/307

[20] Piel, A. (2011) Plasma Physics. Springer, Berlin.

[21] Kolobov, V.I. (2013) Physics of Plasmas, 20, Article ID: 101610.

https://doi.org/10.1063/1.4823472

[22] Fruchtman, A. (2009) Plasma Sources Science and Technology, 18, Article ID: 025033.

https://doi.org/10.1088/0963-0252/18/2/025033

[23] Simon, A. (1955) Physical Review, 98, 317.

https://doi.org/10.1103/PhysRev.98.317

[24] Galaly, A.R. (2014) Physical Science International Journal, 4, 930-939.

https://doi.org/10.9734/PSIJ/2014/9604

[25] Hutchinson (1990) Principles of Plasma Diagnostics. Cambridge University Press, Cambridge.

[26] Huba, J.D. (2011) Plasma Formulary. Naval Research Laboratory Report NRL/PU/ 6790-94-265.

[27] Druyvesteyn, M.J. and Penning, F.M. (1940) Reviews of Modern Physics, 12, 87.

https://doi.org/10.1103/RevModPhys.12.87

[28] Von Engel, A. and Corrigan, S.J.B. (1958) Proceedings of the Physical Society, 72, 786.

https://doi.org/10.1088/0370-1328/72/5/314

[29] Von Engel, A. and Cozens, J.R. (1965) International Journal of Electronics, 19, 61-68.

https://doi.org/10.1080/00207216508937799

[30] Von Engel, A. (1956) Proceedings of the Physical Society, 59, 468.

[31] Bowman, F. (1958) Introduction to Bessel Conditions. Dover, New York.

[32] Korenev, B.G. (2002) Bessel Condition and Their Applications. CRC Press Company, Boca Raton, London, New York, Washington DC.

[33] Arfken, G. (1985) Bessel Conditions. In: Mathematical Methods for Physicists, 3rd Edition, Academic Press, Orlando, Ch. 11, 573-636.

https://doi.org/10.1016/B978-0-12-059820-5.50019-7

[34] Gray, A. and Mathews, G.B. (1966) A Treatise on Bessel Conditions and Their Applications to Physics. 2nd Edition, Dover, New York.

[35] Watson, G.N. (1966) A Treatise on the Theory of Bessel conditions. 2nd Edition, Cambridge University Press, Cambridge.