**1. Introduction**

The concern of the amelioration of the plasma reactor is always a domain important in the development technology, among these fields we find glow discharge. A plasma technology in a gas mixture has been studied by several authors. Ono et al. [1] have been studied oxygen-nitrogen gas mixture glow discharge plasma by intervene many chemical reactions in their model. Khomich et al. [2] have been treated the problem of the atomic deposition in the metal surface modification by nitrogen-argon mixture glow discharge in abnormal regime. Ponduri et al. [3] have been analyzed the dissociation of CO2 by dielectric barrier glow discharge, as a consequence of utilization of CO2 gas a lot of kinds of species intervene in the phenomena discharge. Baadj et al. have been [4] investigated Xe-Cl2 gas mixture for the formation of XeCl\* exciplex lamp by means of zerodimensional model. Li et al. [5] have been studied the plasma jet length in Ne, Ar, He and Kr in atmospheric pressure when the excimer molecule formatted from

metastable state of the atomic gases, and they have been identified three modes versus of the gas flow rate, its about laminar, transition, and turbulent jet modes. for electrons, and to give the effect of the electric field on the particles charged, the

*Calculation of the Metastable Atom Densities in Argon and Neon Abnormal Glow Discharges*

The chemical reactions intervene in the discharge are indicated in the **Table 1**. After that, the model in the 1D Cartesian geometry, writes by partial differential

> *<sup>o</sup>* <sup>þ</sup> *nmKio m*

*<sup>o</sup>* <sup>þ</sup> *nmKio m*

*<sup>m</sup>* � *nmKio m* � � � <sup>2</sup>*nmnmKci* � *nm*

� � <sup>þ</sup> *nmnmKci,* (1)

� � <sup>þ</sup> *nmnmKci,* (2)

*<sup>m</sup>* <sup>þ</sup> *<sup>ε</sup>cinmnmKci*

*<sup>o</sup>* <sup>þ</sup> *<sup>ε</sup>io* � *<sup>ε</sup><sup>m</sup>* � �*nmKio*

�<sup>1</sup> [31] is the rate coefficient of chemo-ionization

� � *,* (4)

*τm*

*<sup>o</sup>* is the rate coefficient

*<sup>m</sup>* is the rate coeffi-

*<sup>m</sup>* is the rate coefficient of de-

*m*

ð Þ *n*<sup>þ</sup> � *ne* (5)

*,* (3)

*<sup>∂</sup><sup>x</sup>* <sup>¼</sup> *ne noKio*

*<sup>∂</sup><sup>x</sup>* <sup>¼</sup> *ne noKio*

*∂*2 *V <sup>∂</sup>x*<sup>2</sup> ¼ � *<sup>e</sup>*

of the metastable atoms, electrons and ions, respectively. *K<sup>m</sup>*

of excited atoms, no is the constant background gas density, *Kio*

*<sup>o</sup>* � *nmK<sup>o</sup>*

*<sup>∂</sup><sup>x</sup>* ¼ �*eφeE* <sup>þ</sup> *<sup>ε</sup>mnenmKo*

*<sup>o</sup>* <sup>þ</sup> *<sup>ε</sup>ionoKio*

*εo*

Here ne, n+, nm, *φm*, *φ<sup>e</sup>* and *φ*<sup>þ</sup> are number densities, densities of transport flux

of electron-impact excitation of ground state atoms, *<sup>ε</sup><sup>m</sup>* <sup>¼</sup> <sup>11</sup>*:*55 eV is the energy loss

cient of electron-impact ionization of excited atoms with the energy loss *<sup>ε</sup>io* � *<sup>ε</sup><sup>m</sup>* � �,

*<sup>o</sup>* is the rate coefficient of electron-impact ionization of ground state atoms and

excitation of excited atoms by electron collisions. *<sup>E</sup>* ¼ �*∂V=∂<sup>x</sup>* is the electric field strength. *ε<sup>o</sup>* and e are the permittivity of free space and elementary charge, respec-

processes with the energy gain *<sup>ε</sup>ci* <sup>¼</sup> <sup>2</sup>*ε<sup>m</sup>* � *<sup>ε</sup>io*. *<sup>τ</sup><sup>m</sup>* <sup>¼</sup> <sup>1</sup> <sup>μ</sup>s is the metastable lifetime. *<sup>ε</sup><sup>e</sup>* is the mean electron energy, *φe<sup>ε</sup>* is the electron energy flux. V is the electrostatic potential. *Pec* is the energy loss per electron due to elastic collision of electrons with

Momentum transfer equations for metastable atoms, electrons, ions and electron

Poisson equation is included in the model [29, 30].

*∂ne ∂t* þ *∂φe*

*DOI: http://dx.doi.org/10.5772/intechopen.84293*

*<sup>∂</sup>n*<sup>þ</sup> *∂t* þ *<sup>∂</sup>φ*<sup>þ</sup>

*∂εene ∂t* þ *∂φe<sup>ε</sup>*

*<sup>∂</sup><sup>x</sup>* <sup>¼</sup> *ne noK<sup>m</sup>*

�*nePec* � *ne <sup>ε</sup>mnoKm*

*<sup>ε</sup>io* <sup>¼</sup> <sup>15</sup>*:*76 eV is the energy loss of ionized atoms, *<sup>K</sup><sup>o</sup>*

*∂nm ∂t* þ *∂φ<sup>m</sup>*

tively. *<sup>K</sup>ci* <sup>¼</sup> <sup>8</sup>*:*<sup>1</sup> � <sup>10</sup>�<sup>10</sup> cm3 <sup>s</sup>

the background gas [32].

energy [33–35] are:

**Table 1.**

**77**

*Kinetic scheme of processes.*

equations:

*Kio*

In this chapter, an research is through concerning the role of neon and argon metastable atoms in the discharge. Metastable atoms have been considered by several authors, both theoretically and experimentally. Metastable densities can experimentally be measured by optical absorption method. In theoretical means, a poise equation, including different production and loss terms is assembled to compute the metastable densities. Experimental measurements were executed, for example, for Neon gas in an RF glow discharge by Eckstein et al. [6], for Ar gas in a microwave boosted glow discharge by Uzelac and Leis [7], and for He gas as a function of discharge conditions by Browne and Dunn [8]. In the research of Smith et al. [9] purely relative absorption signals of Argon metastable atoms as a function of current and pressure were experimented. In the research of Strauss et al. [10] and Ferreira et al. [11] Argon metastable densities were experimented in a afterglow of a pulsed discharge and in the Grimm-type glow discharge, respectively, and some mechanisms for the structure of Argon metastable atoms were recommended. Argon metastable densities have also been experimented by Ferreira et al. [12], and by Ferreira and Ricard [13].

A comparison was complete with Argon metastable densities determined from a coupled-electron-metastable-atom model. Hardy and Sheldon [14] have examined Argon, Helium and Neon gases. A comparison between calculated and measured metastable densities was also realized by Kubota et al. [15] for Helium in a RF and DC glow discharge. Lymberopoulos and Economou [16] have established a combined fluid model for the electrons, Ar metastable atoms, and Ar ions in order to study the effect of metastable atoms in the discharge. In the works [17–22] rate constants of a number of collision processes in control of the demolition of metastable atoms were got by combining equilibrium equations with the experimented time-dependent variation of the metastable densities or by analyzing the reliance of the decay constants upon pressure. Den Hartog et al. [23, 24] have investigated Helium gas.

Last-mentioned, Fedoseev and Sukhinin [25] have investigated the influence of metastable Ar atoms on gas discharge plasma with dust particles. Shumova et al. [26] have investigated the effect of metastable Ne atoms and dust particles in a positive column of glow discharge.

The aim of this work is to present, influence of the discharge characteristic in the case of input data obtained by BOLSIG+ code, and in the case of input data taken from the Boltzmann equation in multi term approximation (BMA). We note that both these approaches are widely used. For simulations of positive column of glow discharge, Vasilyak et al. [27] used the first approach, while Sukhinin et al. [28] used the second approach. In Section 2, the mathematical model is delineated; it comprises the boundary and initial conditions as well as the numerical method. In Section 3, the results are discussed for Ar discharge. In Section 4, the test scheme is given. In Section 5, Influence of the input data got from BOLSIG+ on the argon discharge has been shown. In Section 6, the characteristics of the Ne discharge with input data taken from BOLSIG+ are presented. Finally, the conclusion of the chapter is given in Section 7.

### **2. Discharge modeling**

Our mathematical model builds on the first three moments of the Boltzmann equation. The Continuity equations and momentum transfer equations of metastable atom, electrons and positively charged ions. The energy equation is known only *Calculation of the Metastable Atom Densities in Argon and Neon Abnormal Glow Discharges DOI: http://dx.doi.org/10.5772/intechopen.84293*

for electrons, and to give the effect of the electric field on the particles charged, the Poisson equation is included in the model [29, 30].

The chemical reactions intervene in the discharge are indicated in the **Table 1**.

After that, the model in the 1D Cartesian geometry, writes by partial differential equations:

$$\frac{\partial n\_{\epsilon}}{\partial t} + \frac{\partial \rho\_{\epsilon}}{\partial \mathbf{x}} = n\_{\epsilon} \left( n\_{o} K\_{o}^{io} + n\_{m} K\_{m}^{io} \right) + n\_{m} n\_{m} K^{ci},\tag{1}$$

$$\frac{\partial n\_{+}}{\partial t} + \frac{\partial \rho\_{+}}{\partial \mathbf{x}} = n\_{\epsilon} \left( n\_{o} \mathbf{K}\_{o}^{io} + n\_{m} \mathbf{K}\_{m}^{io} \right) + n\_{m} n\_{m} \mathbf{K}^{ci},\tag{2}$$

$$\frac{\partial n\_m}{\partial t} + \frac{\partial \rho\_m}{\partial \mathbf{x}} = n\_\varepsilon \left( n\_o K\_o^m - n\_m K\_m^o - n\_m K\_m^{io} \right) - 2 n\_m n\_m K^{ci} - \frac{n\_m}{\tau\_m},\tag{3}$$

$$\begin{aligned} \frac{\partial \varepsilon\_{\varepsilon} n\_{\varepsilon}}{\partial t} + \frac{\partial \rho\_{\varepsilon \varepsilon}}{\partial \mathbf{x}} &= -\varepsilon \rho\_{\varepsilon} E + \varepsilon^{m} n\_{\varepsilon} n\_{m} K\_{m}^{o} + \varepsilon^{ci} n\_{m} n\_{m} K^{ci} \\ -n\_{\varepsilon} P^{\varepsilon} &- n\_{\varepsilon} \left( \varepsilon^{m} n\_{o} K\_{o}^{m} + \varepsilon^{io} n\_{o} K\_{o}^{io} + \left( \varepsilon^{io} - \varepsilon^{m} \right) n\_{m} K\_{m}^{io} \right) \end{aligned} , \tag{4}$$

$$(\varepsilon^{m}n\_{o}K\_{o}^{m}+\varepsilon^{io}n\_{o}K\_{o}^{io}+\left(\varepsilon^{jo}-\varepsilon^{m}\right)n\_{m}K\_{m}^{io})$$

$$\frac{\partial^{2}V}{\partial\mathbf{x}^{2}}=-\frac{e}{\varepsilon\_{o}}\left(n\_{+}-n\_{e}\right)\tag{5}$$

Here ne, n+, nm, *φm*, *φ<sup>e</sup>* and *φ*<sup>þ</sup> are number densities, densities of transport flux of the metastable atoms, electrons and ions, respectively. *K<sup>m</sup> <sup>o</sup>* is the rate coefficient of electron-impact excitation of ground state atoms, *<sup>ε</sup><sup>m</sup>* <sup>¼</sup> <sup>11</sup>*:*55 eV is the energy loss of excited atoms, no is the constant background gas density, *Kio <sup>m</sup>* is the rate coefficient of electron-impact ionization of excited atoms with the energy loss *<sup>ε</sup>io* � *<sup>ε</sup><sup>m</sup>* � �, *Kio <sup>o</sup>* is the rate coefficient of electron-impact ionization of ground state atoms and *<sup>ε</sup>io* <sup>¼</sup> <sup>15</sup>*:*76 eV is the energy loss of ionized atoms, *<sup>K</sup><sup>o</sup> <sup>m</sup>* is the rate coefficient of deexcitation of excited atoms by electron collisions. *<sup>E</sup>* ¼ �*∂V=∂<sup>x</sup>* is the electric field strength. *ε<sup>o</sup>* and e are the permittivity of free space and elementary charge, respectively. *<sup>K</sup>ci* <sup>¼</sup> <sup>8</sup>*:*<sup>1</sup> � <sup>10</sup>�<sup>10</sup> cm3 <sup>s</sup> �<sup>1</sup> [31] is the rate coefficient of chemo-ionization processes with the energy gain *<sup>ε</sup>ci* <sup>¼</sup> <sup>2</sup>*ε<sup>m</sup>* � *<sup>ε</sup>io*. *<sup>τ</sup><sup>m</sup>* <sup>¼</sup> <sup>1</sup> <sup>μ</sup>s is the metastable lifetime. *<sup>ε</sup><sup>e</sup>* is the mean electron energy, *φe<sup>ε</sup>* is the electron energy flux. V is the electrostatic potential. *Pec* is the energy loss per electron due to elastic collision of electrons with the background gas [32].

Momentum transfer equations for metastable atoms, electrons, ions and electron energy [33–35] are:


**Table 1.** *Kinetic scheme of processes.*

metastable state of the atomic gases, and they have been identified three modes versus of the gas flow rate, its about laminar, transition, and turbulent jet modes. In this chapter, an research is through concerning the role of neon and argon metastable atoms in the discharge. Metastable atoms have been considered by several authors, both theoretically and experimentally. Metastable densities can experimentally be measured by optical absorption method. In theoretical means, a poise equation, including different production and loss terms is assembled to compute the metastable densities. Experimental measurements were executed, for example, for Neon gas in an RF glow discharge by Eckstein et al. [6], for Ar gas in a microwave boosted glow discharge by Uzelac and Leis [7], and for He gas as a function of discharge conditions by Browne and Dunn [8]. In the research of Smith et al. [9] purely relative absorption signals of Argon metastable atoms as a function of current and pressure were experimented. In the research of Strauss et al. [10] and Ferreira et al. [11] Argon metastable densities were experimented in a afterglow of a pulsed discharge and in the Grimm-type glow discharge, respectively, and some mechanisms for the structure of Argon metastable atoms were recommended. Argon metastable densities have also been experimented by Ferreira et al. [12], and

*Solid State Physics - Metastable, Spintronics Materials and Mechanics of Deformable…*

A comparison was complete with Argon metastable densities determined from a coupled-electron-metastable-atom model. Hardy and Sheldon [14] have examined Argon, Helium and Neon gases. A comparison between calculated and measured metastable densities was also realized by Kubota et al. [15] for Helium in a RF and DC glow discharge. Lymberopoulos and Economou [16] have established a combined fluid model for the electrons, Ar metastable atoms, and Ar ions in order to study the effect of metastable atoms in the discharge. In the works [17–22] rate constants of a number of collision processes in control of the demolition of metastable atoms were got by combining equilibrium equations with the experimented time-dependent variation of the metastable densities or by analyzing the reliance of the decay constants upon pressure. Den Hartog et al. [23, 24] have investigated

Last-mentioned, Fedoseev and Sukhinin [25] have investigated the influence of metastable Ar atoms on gas discharge plasma with dust particles. Shumova et al. [26] have investigated the effect of metastable Ne atoms and dust particles in a

The aim of this work is to present, influence of the discharge characteristic in the case of input data obtained by BOLSIG+ code, and in the case of input data taken from the Boltzmann equation in multi term approximation (BMA). We note that both these approaches are widely used. For simulations of positive column of glow discharge, Vasilyak et al. [27] used the first approach, while Sukhinin et al. [28] used the second approach. In Section 2, the mathematical model is delineated; it comprises the boundary and initial conditions as well as the numerical method. In Section 3, the results are discussed for Ar discharge. In Section 4, the test scheme is given. In Section 5, Influence of the input data got from BOLSIG+ on the argon discharge has been shown. In Section 6, the characteristics of the Ne discharge with input data taken from BOLSIG+ are presented. Finally, the conclusion of the chap-

Our mathematical model builds on the first three moments of the Boltzmann equation. The Continuity equations and momentum transfer equations of metastable atom, electrons and positively charged ions. The energy equation is known only

by Ferreira and Ricard [13].

positive column of glow discharge.

ter is given in Section 7.

**2. Discharge modeling**

**76**

Helium gas.

$$
\rho\_{\epsilon} = -n\_{\epsilon}\mu\_{\epsilon}E - \frac{\partial D\_{\epsilon}n\_{\epsilon}}{\partial \mathbf{x}},\tag{6}
$$

The applied potential at the cathode is 400 V. The uniform subdivision of the space interval in 250 elements and a time step Δt = 10 ps, have been utilized.

*Calculation of the Metastable Atom Densities in Argon and Neon Abnormal Glow Discharges*

zone busied the time simulation amidst in 3 <sup>10</sup><sup>5</sup> and 3.6 <sup>10</sup><sup>5</sup> s.

*DOI: http://dx.doi.org/10.5772/intechopen.84293*

**Figure 1.**

**79**

**Figure 1** shows the temporal progression of the potential (**Figure 1(a)**), concentration of electrons (**Figure 1(b)**), concentration of ions (**Figure 1(c)**), concentration of metastable atoms (**Figure 1(d)**), the electric field (**Figure 1(e)**) and the mean electron energy (**Figure 1(f)**). We remark that the discharge is distinguished by three zones: the first one busied the time simulation amidst in 10<sup>11</sup> and 9 <sup>10</sup><sup>6</sup> s, the second one took the time simulation inter-time of 9 <sup>10</sup><sup>6</sup> and 3 <sup>10</sup><sup>5</sup> s, the final

For the first zone, we comment that the metastable atom, electrons and ion concentrations are almost identical. Therefore, the net space charge concentration is unimportant. The electric potential is distinguished by the Laplace form due to the net space charge concentration that is exist. Consequently, the electric field is seemingly constant, besides that the mean electron energy is constant.

In the second zone, we remark a pseudo emergence of the cathodic region, this is characterized by a significance of ion concentration and unimportant of the electron

*Argon spatiotemporal distributions of electric potential (a), electrons volume number (b), ions volume number*

*(c), metastable atom density (d), electric field (e), and mean electron energy (f) at 400 V.*

$$
\rho\_{+} = n\_{+} \mu\_{+} E - \frac{\partial D\_{+} n\_{+}}{\partial \mathbf{x}},
\tag{7}
$$

$$
\rho\_m = -D\_m \frac{\partial n\_m}{\partial \mathbf{x}},\tag{8}
$$

$$
\rho\_{\epsilon\varepsilon} = -n\_{\epsilon}E\mu\_{\epsilon\varepsilon} - \frac{\partial n\_{\epsilon}D\_{\epsilon\varepsilon}}{\partial \mathbf{x}} \tag{9}
$$

Here *μe*, *μ*þ, De and D+ are the electron, ion for mobilities and diffusion coefficients, respectively. *De<sup>ε</sup>* and *μe<sup>ε</sup>* are the diffusivity and mobility of electron energy transport. The ion mobility has brought from Phelps and Petrović [36] and their ion diffusivity is computed agree with Einstein's relation [37]. The coefficients for electrons [38] in argon as dependences on the mean electron energy are got as of INP Greifswald for direct evaluation with the results obtained by Becker et al. [30]. Dm is the metastable atom diffusivity someplace *noDm* <sup>¼</sup> <sup>1</sup>*:*7*x*10<sup>18</sup> cm�<sup>1</sup> <sup>s</sup> �<sup>1</sup> [19].

#### **2.1 Boundary and initial conditions**

The discharge is affected between two parallel plate electrodes and the radius of the electrode is presumed to be higher than the electrode gap and the physical characteristic distributions are approximately uniform along the radial direction. The grounded electrode has been put at *x* ¼ 1 cm, play the role as the anode (*Vanode* ¼ 0). The powered electrode has been put at *x* ¼ 0 cm, which initiates the model discretization, play the role as the cathode (*Vcathode* ¼ �*VDC*).

At time *t* ¼ 0, the metastable atom, electron, and ion densities are supposed constant and equal to 10<sup>3</sup> cm�<sup>3</sup> , and the mean electron energy equal to 1 eV.

Presuming disappearance of the metastable atom density, i.e., *nm* ¼ 0 has been arranged at the cathode, whilst a predominant field-driven flux shut to the cathode the expression *<sup>∂</sup>D*þ*n*þ*=∂<sup>x</sup>* <sup>¼</sup> 0 has been used for the positive ion density <sup>∀</sup>*<sup>t</sup>* > 0. At the anode, the metastable atoms and the electron density are supposed to be zero. The electron flux separating the cathode is calculated by the expression *φe*ð Þ¼� *x* ¼ 0*; t γφ*þð Þ *x* ¼ 0*; t* ∀*t* > 0, the mean electron energy is assumed to be 5 eV at the cathode [36] and the gas temperature is equal to 273 K in the discharge.

#### **2.2 Numerical method**

For the metastable atom and Poisson equations a finite difference method has been employed. The transport equations of the electron energy, electron and ion are also discretized spatially with the finite difference technique. In this method the exponential scheme has been employed into account [37, 39–42]. The discretization of the time terms by the right position of the finite difference technique has been used. Consequently, every discretized equation is defined as a tridiagonal matrix, which is solved by Thomas's algorithm.

### **3. Results and discussion of argon discharge**

In this part, we will analysis the spatiotemporal evolution of the abnormal glow discharge in the existence of metastable atom density. The gas pressure is 133.32 Pa. The neutral species density is computed from the temperature and gas pressure with the ideal gas law. The constant value for the secondary electron yield is 0.06 [36].

#### *Calculation of the Metastable Atom Densities in Argon and Neon Abnormal Glow Discharges DOI: http://dx.doi.org/10.5772/intechopen.84293*

The applied potential at the cathode is 400 V. The uniform subdivision of the space interval in 250 elements and a time step Δt = 10 ps, have been utilized.

**Figure 1** shows the temporal progression of the potential (**Figure 1(a)**), concentration of electrons (**Figure 1(b)**), concentration of ions (**Figure 1(c)**), concentration of metastable atoms (**Figure 1(d)**), the electric field (**Figure 1(e)**) and the mean electron energy (**Figure 1(f)**). We remark that the discharge is distinguished by three zones: the first one busied the time simulation amidst in 10<sup>11</sup> and 9 <sup>10</sup><sup>6</sup> s, the second one took the time simulation inter-time of 9 <sup>10</sup><sup>6</sup> and 3 <sup>10</sup><sup>5</sup> s, the final zone busied the time simulation amidst in 3 <sup>10</sup><sup>5</sup> and 3.6 <sup>10</sup><sup>5</sup> s.

For the first zone, we comment that the metastable atom, electrons and ion concentrations are almost identical. Therefore, the net space charge concentration is unimportant. The electric potential is distinguished by the Laplace form due to the net space charge concentration that is exist. Consequently, the electric field is seemingly constant, besides that the mean electron energy is constant.

In the second zone, we remark a pseudo emergence of the cathodic region, this is characterized by a significance of ion concentration and unimportant of the electron

**Figure 1.**

*Argon spatiotemporal distributions of electric potential (a), electrons volume number (b), ions volume number (c), metastable atom density (d), electric field (e), and mean electron energy (f) at 400 V.*

*<sup>φ</sup><sup>e</sup>* ¼ �*neμeE* � *<sup>∂</sup>Dene*

*Solid State Physics - Metastable, Spintronics Materials and Mechanics of Deformable…*

*<sup>φ</sup>*<sup>þ</sup> <sup>¼</sup> *<sup>n</sup>*þ*μ*þ*<sup>E</sup>* � *<sup>∂</sup>D*þ*n*<sup>þ</sup>

*<sup>φ</sup>e<sup>ε</sup>* ¼ �*neEμe<sup>ε</sup>* � *<sup>∂</sup>neDe<sup>ε</sup>*

Here *μe*, *μ*þ, De and D+ are the electron, ion for mobilities and diffusion coefficients, respectively. *De<sup>ε</sup>* and *μe<sup>ε</sup>* are the diffusivity and mobility of electron energy transport. The ion mobility has brought from Phelps and Petrović [36] and their ion diffusivity is computed agree with Einstein's relation [37]. The coefficients for electrons [38] in argon as dependences on the mean electron energy are got as of INP Greifswald for direct evaluation with the results obtained by Becker et al. [30].

The discharge is affected between two parallel plate electrodes and the radius of

the electrode is presumed to be higher than the electrode gap and the physical characteristic distributions are approximately uniform along the radial direction. The grounded electrode has been put at *x* ¼ 1 cm, play the role as the anode (*Vanode* ¼ 0). The powered electrode has been put at *x* ¼ 0 cm, which initiates the

At time *t* ¼ 0, the metastable atom, electron, and ion densities are supposed

Presuming disappearance of the metastable atom density, i.e., *nm* ¼ 0 has been arranged at the cathode, whilst a predominant field-driven flux shut to the cathode the expression *<sup>∂</sup>D*þ*n*þ*=∂<sup>x</sup>* <sup>¼</sup> 0 has been used for the positive ion density <sup>∀</sup>*<sup>t</sup>* > 0. At the anode, the metastable atoms and the electron density are supposed to be zero.

*φe*ð Þ¼� *x* ¼ 0*; t γφ*þð Þ *x* ¼ 0*; t* ∀*t* > 0, the mean electron energy is assumed to be 5 eV at the cathode [36] and the gas temperature is equal to 273 K in the discharge.

For the metastable atom and Poisson equations a finite difference method has been employed. The transport equations of the electron energy, electron and ion are also discretized spatially with the finite difference technique. In this method the exponential scheme has been employed into account [37, 39–42]. The discretization of the time terms by the right position of the finite difference technique has been used. Consequently, every discretized equation is defined as a tridiagonal matrix,

In this part, we will analysis the spatiotemporal evolution of the abnormal glow discharge in the existence of metastable atom density. The gas pressure is 133.32 Pa. The neutral species density is computed from the temperature and gas pressure with the ideal gas law. The constant value for the secondary electron yield is 0.06 [36].

, and the mean electron energy equal to 1 eV.

model discretization, play the role as the cathode (*Vcathode* ¼ �*VDC*).

The electron flux separating the cathode is calculated by the expression

*∂nm*

*∂x*

*φ<sup>m</sup>* ¼ �*Dm*

Dm is the metastable atom diffusivity someplace *noDm* <sup>¼</sup> <sup>1</sup>*:*7*x*10<sup>18</sup> cm�<sup>1</sup> <sup>s</sup>

**2.1 Boundary and initial conditions**

constant and equal to 10<sup>3</sup> cm�<sup>3</sup>

**2.2 Numerical method**

**78**

which is solved by Thomas's algorithm.

**3. Results and discussion of argon discharge**

*<sup>∂</sup><sup>x</sup> ,* (6)

*<sup>∂</sup><sup>x</sup> ,* (7)

*<sup>∂</sup><sup>x</sup> ,* (8)

(9)

�<sup>1</sup> [19].

concentration. This is explicated by the velocity of the electron species, which speed a more than the ion species and displaced quickly starting at the cathodic region. Consequently, the amount (ne-n+) is considerable which influence directly the electric potential as consequence the important chute that is present. Inevitably, the electric field is intense. The last earns the electron species an important energy. We comment that the metastable atom concentration is important. This discharge is sustained by the secondary emission coefficient as well as the existence of metastable atom concentration. Sooner than t = 3 <sup>10</sup><sup>5</sup> s we remark a pseudo emergence of the negative glow region, where it is typified by the similar electron and ion concentrations. Consequently the net space charge concentration is slight. Hence, the electric potential and the electric field are constants. Automatically, the metastable atoms concentration is diminished.

In the last zone, we remark three dissimilar regions: the anode and the plasma regions and the cathode region. The anode region is characterized by the ion concentration that is less important compared to the electron concentration. In this zone we remark that the convergence of each physical characteristics of the dis-

*Calculation of the Metastable Atom Densities in Argon and Neon Abnormal Glow Discharges*

In this part, we will study the influence of the voltage and gas pressure on the argon discharge. So, the potential at the cathode is taking of 250 V and we will alter the gas pressure. For the influence of the applied potential on the discharge, we take

**Figure 2(a)** shows the metastable atoms concentration plots depending on the pressure in the study state. The metastable atoms concentration augments with increasing pressure. For elevated pressure the gas density augments, which the electron diffusion coefficient turns out to be fewer and the bulk of the plasma rises which the both sheaths of the anode and cathode turn out to be small. These circumstances of the charged particle manipulate on the metastable atoms behavior in the study stated, i.e., the cathodic region is overflowing with the electron and metastable atom concentrations, which go faster the ion species in the existence of

**Figure 2(b)** shows the metastable atoms concentration graphs depending on the

. We judged these results with calculated values established in

. Consequently, our calculated

potential in the study state. For elevated potential the excitation and ionization processes increase, and the charged particle turns out to be raised in the stationary state. Consequently, the metastable atom concentration graphs become growing. The greatest of the metastable atom concentration varies amidst in 2.47 <sup>10</sup><sup>10</sup>

the literature [9, 10, 12, 15]. Depending on the discharge circumstances, all these

values something like in the exact range of the size order. We find that the values of

*Comparison between the results obtained by our calculation using a database of BMA for pd = 1 Torr cm and*

charge are reached at the time 3.6 <sup>10</sup><sup>5</sup> s.

*DOI: http://dx.doi.org/10.5772/intechopen.84293*

the gas pressure at 133.32 Pa. bulk.

the electric field.

and 6.63 <sup>10</sup><sup>11</sup> cm<sup>3</sup>

**Figure 3.**

**81**

**3.1 Influence of the voltage and gas pressure**

value varies amidst in 2 <sup>10</sup><sup>10</sup> and 5 <sup>10</sup><sup>13</sup> cm<sup>3</sup>

*those given by experimental for pd = 0.5 Torr cm in Argon gas.*

**Figure 2.** *Metastable atom densities as a function of pressure (a) and as a function of voltage (b) in the stationary state in Argon gas.*

In the last zone, we remark three dissimilar regions: the anode and the plasma regions and the cathode region. The anode region is characterized by the ion concentration that is less important compared to the electron concentration. In this zone we remark that the convergence of each physical characteristics of the discharge are reached at the time 3.6 <sup>10</sup><sup>5</sup> s.

### **3.1 Influence of the voltage and gas pressure**

concentration. This is explicated by the velocity of the electron species, which speed a more than the ion species and displaced quickly starting at the cathodic region. Consequently, the amount (ne-n+) is considerable which influence directly the electric potential as consequence the important chute that is present. Inevitably, the electric field is intense. The last earns the electron species an important energy. We comment that the metastable atom concentration is important. This discharge is sustained by the secondary emission coefficient as well as the existence of metastable atom concentration. Sooner than t = 3 <sup>10</sup><sup>5</sup> s we remark a pseudo emergence of the negative glow region, where it is typified by the similar electron and ion concentrations. Consequently the net space charge concentration is slight. Hence, the electric potential and the electric field are constants. Automatically, the meta-

*Solid State Physics - Metastable, Spintronics Materials and Mechanics of Deformable…*

*Metastable atom densities as a function of pressure (a) and as a function of voltage (b) in the stationary state in*

stable atoms concentration is diminished.

**Figure 2.**

*Argon gas.*

**80**

In this part, we will study the influence of the voltage and gas pressure on the argon discharge. So, the potential at the cathode is taking of 250 V and we will alter the gas pressure. For the influence of the applied potential on the discharge, we take the gas pressure at 133.32 Pa. bulk.

**Figure 2(a)** shows the metastable atoms concentration plots depending on the pressure in the study state. The metastable atoms concentration augments with increasing pressure. For elevated pressure the gas density augments, which the electron diffusion coefficient turns out to be fewer and the bulk of the plasma rises which the both sheaths of the anode and cathode turn out to be small. These circumstances of the charged particle manipulate on the metastable atoms behavior in the study stated, i.e., the cathodic region is overflowing with the electron and metastable atom concentrations, which go faster the ion species in the existence of the electric field.

**Figure 2(b)** shows the metastable atoms concentration graphs depending on the potential in the study state. For elevated potential the excitation and ionization processes increase, and the charged particle turns out to be raised in the stationary state. Consequently, the metastable atom concentration graphs become growing.

The greatest of the metastable atom concentration varies amidst in 2.47 <sup>10</sup><sup>10</sup> and 6.63 <sup>10</sup><sup>11</sup> cm<sup>3</sup> . We judged these results with calculated values established in the literature [9, 10, 12, 15]. Depending on the discharge circumstances, all these value varies amidst in 2 <sup>10</sup><sup>10</sup> and 5 <sup>10</sup><sup>13</sup> cm<sup>3</sup> . Consequently, our calculated values something like in the exact range of the size order. We find that the values of

#### **Figure 3.**

*Comparison between the results obtained by our calculation using a database of BMA for pd = 1 Torr cm and those given by experimental for pd = 0.5 Torr cm in Argon gas.*

the current densities are 0*:*137, 0*:*508 and 0*:*843 mA/cm<sup>2</sup> related to the pressures of 133*:*32, 250*:*32 and 330*:*32 Pa, in that order.

**Figure 3** represents the current-voltage properties in the stationary state. The results attained beginning of the database of BMA evaluated to those obtained by experimental method [43, 44]. We find that the results get beginning of the database of the BMA are in excellent conformity judgment against to those experimental results [43, 44].
