Glow Discharge in a High-Velocity Air Flow: The Role of the Associative Ionization Reactions Involving Excited Atoms

A kinetic scheme for non-equilibrium regimes of atmospheric pressure air discharges is developed. A distinctive feature of this model is that it includes associative ionization with the participation of N(2D, 2P) atoms. The thermal dissociation of vibrationally excited nitrogen molecules and the electronic excitation from all the vibrational levels of the nitrogen molecules are also accounted for. The model is used to simulate the parameters of a glow discharge ignited in a fast longitudinal flow of preheated (T0 = 1800–2900 K) air. The results adequately describe the dependence of the electric field in the glow discharge on the initial gas temperature. For T0 = 1800 K, a substantial acceleration in the ionization kinetics of the discharge is found at current densities larger than 3 A/cm2, mainly due to the N(2P) + O(3P) → NO+ + e process; being the N(2P) atoms produced via quenching of N2(A3∑u+) molecules by N(4S) atoms. Correspondingly, the reduced electric field noticeably falls because the electron energy (6.2 eV) required for the excitation of the N2(A3∑u+) state is considerably lower than the ionization energy (9.27 eV) of the NO molecules. For higher values of T0, the associative ionization N(2D) + O(3P) → NO+ + e process (with a low–activation barrier of 0.38 eV) becomes also important in the production of charged particles. The N(2D) atoms being mainly produced via quenching of N2(A3∑u+) molecules by O(3P) atoms.


Introduction
A number of experiments have been reported on non-equilibrium regimes of discharges in atmospheric-pressure air, in particular, glow-type discharges in open ambient air (e.g., [1][2][3][4][5][6][7][8][9][10][11][12]), and in fast longitudinal flows of air [13,14]. In [13,14], the parameters of a low-current glow discharge in a high-velocity flow of preliminary heated air (T 0 = 1800-2900 K) at atmospheric pressure, is investigated. The discharge is ignited between two needle electrodes oriented along the axis of the gas flow directed from the cathode to the anode. High-speed gas flows have been used to provide sufficient cooling of discharges at high pressures [15,16]. If the gas residence time in the discharge is small as compared with the vibrational-translational (V-T) relaxation time, the gas heating is almost suppressed and the gas is in a strongly non-equilibrium state (i.e., characterized by a level of vibrational energy which considerable exceeds its equilibrium value). These regimes, which correspond to relatively high electron number densities, 10 18 -10 19 m −3 , and relatively low gas temperatures, 2000-3000 K, are of interest for many practical applications, including plasma decontamination and sterilization, material processing, modification of electromagnetic waves propagation, and plasma aerodynamics (see the

Ion-Ion Recombination
T e and T g units are in Kelvin.
Among the processes of charged particle production, the model incorporates the associative ionization reaction with the participation of excited atoms (R18), which have no activation barrier; and the near-threshold reaction (R16), with a low-activation barrier of 0.38 eV. The efficiency of these processes is determined by the decrease in ionization threshold resulting from the chemical bond energy of the compound molecular ion. The reaction between the N( 2 P) and O( 3 P) atoms was observed in experiments [54]; the rate constant of the decay of the N( 2 P) atoms in the reaction with O( 3 P) was found to be (1.7 ± 0.4) 10 -17 m 3 /s, with the formation of NO + as a significant reaction channel. This value agrees with estimate [36]. Cross section data for the ionization reaction between the N( 2 D) and O( 3 P) were measured in [39]. Following [25], the rate coefficient of the reaction (R13) was taken from [35]. The calculations are performed using rates for reactions (R22) and (R23) taken from [44]. The coefficients for electron-impact processes are calculated for the mixture N 2 -20%O 2 (the rate constants for electrons depend on the mixture composition, i.e. on the degree of dissociation of oxygen, only slightly [20,25]). The influence of the deformation of the EEDF (due to superelastic collisions) on the high-threshold electron-impact processes (electronic excitation, dissociation, ionization and dissociative attachment), is accounted for as in [15,20,34]. The increase of the electronic excitation rates (including the electron-impact dissociation) due to the inclusion of excitation from all the vibrational levels of N 2 (X 1 g + , v) (v = 0 to 67, according to [55,56]) is considered as [57][58][59], is the constant for the jth excitation process from vibrational ground-state N 2 (X 1 g + , v = 0) molecules calculated with the BOLSIG+ code, ω is the nitrogen vibrational quantum and k b the Boltzmann's constant. T v is the vibrational temperature of the nitrogen molecules related to its mean-vibrational energy as (1) results from the simple assumption of reducing the threshold of the cross sections for the corresponding N 2 (X 1 g process according to the excitation energy of the reactants [60][61][62] and of a Boltzmann distribution for the populations of N 2 (X 1 g [58]. The conditions of vibrational excitation also facilitate the thermal dissociation of the nitrogen. The dependence of the rate of the thermal dissociation of N 2 (X 1 g + , v) on the vibrational temperature T v is take into account through the Macheret-Fridman model [63,64], which was found to be the most accurate for nitrogen dissociation under non-equilibrium conditions [65]. The vibrational temperatures of other molecular air components, such as O 2 and NO, are close to the gas temperature due to a fast V-T relaxation of these species (e.g., [15]).
For the conditions of [13,14], the preheated airflow enters the discharge region with a mean-velocity u 0 ≈ 450 m/s, the discharge length along the flow L = 3.5 cm, the discharge radius R ≈ 1.6 mm, the gas pressure p = 1 atm. The Reynolds number, based on the nozzle diameter (D = 1 cm), is 1× 10 4 -2 × 10 4 . For such conditions, the radial transport of particles (and also heat) due to radial inhomogeneities is negligible small as compared to the longitudinal (along the axis) convective transport caused by axial inhomogeneities. The radial diffusion timescale for neutral and charged species may be estimated as 15]) accounts for the gas turbulence effects as the gas flow is likely fully turbulent [66]. (Note that D t is several times greater than that corresponding to laminar flows [15].) It follows that the radial diffusion timescale R 2 /D t is about 1.3 × 10 −4 s under the analyzed conditions. This value can be compared to the mean-convective transport time of the gas in the discharge region, τ ≈ u 0 /(2d) ≈ 3.9 × 10 −5 s. The factor 2 in the estimation of τ, considers that, on average, the gas is transported along half the length of the discharge [15]. The inequality R 2 /D t >> τ indicates that the transport of particles in the radial direction is negligible as compared to the longitudinal transport caused by convection. See also estimates in [23]. Note that spatially resolved optical measurements of N 2 C state emission in the conditions of [13,14] show that the radius of the discharge is nearly constant along the discharge axis (i.e., the radial diffusion of electrons from the discharge core is small); thus, supporting the above estimation. In addition, the diffusion timescale of heat may be estimated as N c p R 2 /λ, where c p is the mean-specific heat at constant pressure per heavy particle [67] and λ the translational thermal conductivity of heavy particles. For T g = 2000-3000 K, λ is about 0.1 W/(s m) [68]. The corresponding ratio N c p R 2 /λ/τ is of the order of 100.
To reveal the influence of the associative ionization with the participation of exited atoms, on the air ionization kinetics in the conditions of [13,14], the balance equations are solved in a local (volume-averaged) approximation [69]; but the longitudinal inhomogeneities caused by convection is also accounted for. This is done by introducing a term (Y i0 − Y i )/τ [15,70] in the balance equation for each Y i plasma quantity; where Y i0 is the value of the Y i quantity at the section where the preheated air enters the discharge column (just after the cathodic part of the discharge).
The balance equations for the plasma particles under the local approximation, where [Y i ] is the number density and S i j is the rate of production of the Y i species in the jth reaction (negative if the species is destroyed), are solved for the species N 2 (A 3 u + ), N 2 (a' 1 u − ), NO, N( 4 S), -, and electrons. The rate to produce the N 2 (A 3 u + ) state by cascading is assumed to be equal to the sum of the rates for the production by electron impact of the N 2 (B 3 Π g ) and N 2 (C 3 Π u ) states [71,72]. Estimates show [71] that this simplification causes a negligible error in the rate of production of the N 2 (A 3 u + ) state. By making this assumption, the populations of the N 2 (B 3 Π g ) and N 2 (C 3 Π u ) states do not need to be computed directly. The density of the dominating sort of positive ions NO + (all primary ions formed in ionizing collisions of electrons with air components other than NO convert quickly to NO + under the conditions considered) is obtained from the condition of quasi-neutrality. The density of the dominant species N 2 (X 1 g + , v) is obtained from the constancy of the pressure, while the density of O 2 is calculated by the condition of conservation of N and O nuclei. The balance equations for the plasma particles are coupled with the balance equations describing the mean-vibrational energy of the nitrogen molecules and the mean-kinetic energy of the gas, η V and η T are the fractions of electron energy transferred to gas heating and to the vibrational excitation of N 2 (X 1 g + , v) molecules, and σ is the electrical conductivity of the plasma. ε v (T g ) is the equilibrium value of the mean-vibrational energy of the nitrogen molecule. The fraction η V of electron energy going to vibrational excitation of N 2 (X 1 g + , v) is determined from the corresponding electron energy loss coefficients given by the BOLSIG+ software. The fraction η T is determined in the same manner but considering the electron energy loss coefficients corresponding to the excitation of the translational and rotational modes of molecules in air, and the vibrational excitation of O 2 followed by a fast V-T relaxation. τ VT is the timescale of the V-T energy relaxation by molecules and atoms collisions [34]. The V-T relaxation of N 2 (X 1 g + , v) molecules on O( 3 P) atoms is described by using the experimental data presented in [73][74][75]. c p is the specific heat of the heavy particle of the Y i species calculated by assuming that the translational and rotational energies stored per particle are equal to their classical values (= 7/2 k b for a diatomic molecule because temperatures are equilibrated in air discharges [34]). Q R is the 'fast' gas heating rate term due to the electronic-translational (E-T) relaxation energy suggested in [76], to describe observations at moderate values of reduced electric fields. Reactions where energy release is considered for the gas heating term Q R in Equation (4) are accompanied by the exothermic energy ε R in the right side of the equations in Table 1. The current density j of the discharge can be related with the electric field strength in the discharge column by using the Ohm's law. The balance equations are solved numerically by a finite-difference explicit method with the second-order approximation in time. Because of the stiffness of the equations (there is a wide range of timescales related with different plasma processes, which reduces the time-step needed for accurate numerical integration), a short time-step for integration (= 1.0 × 10 -10 s) is used. At the initial time-instant (t = 0), it is assumed that the reduced electric field has a value in the range 100-200 Td, intermediate between relatively low E/N in the column and high E/N in the cathode part. (A variation of this value inside the above interval has a negligible effect on the calculated discharge parameters [23]). The electron number density for a given current density value is determined by the Ohm's law. The initial densities of negative ions are estimated using the local balance equations, and the density of positive ions is obtained from the condition of quasi-neutrality. According to [13,14], the preheated air enters the discharge column in a state close to the thermodynamic equilibrium at T 0 , hence, the densities of neutral species at t = 0 are taken to be close to the equilibrium values at T 0 . The initial conditions for Equations (3) and (4) are as follows: T g (0) = T 0 , ε v (0) = ε v (T 0 ). The balance equations are integrated from these initial values up to times of about 10 −3 s (it was sufficient for the density of each species to converge within an error of about 10 −3 to its equilibrium value). The gas flow velocity u is assumed to vary proportionally to T g , to consider the increase of gas temperature due to gas heating inside the discharge [49]. In [13,14], the current density values are measured by emission spectroscopy of the N 2 C state, and thus are representative of the discharge region with high electron number density. That is, the model output gives averaged values of the plasma quantities over the discharge axis, but corresponding to the central region of the current-carrying area of the discharge.

Results and Discussion
Gas and vibrational temperatures calculated at various values of T 0 versus the discharge current density is given in Figure 1. The gas and vibrational temperatures increase with j. For j > 1 A/cm 2 the values of T v are larger at smaller T 0 , thus increasing the non-equilibrium state of the discharge. This is expected since the rate of V-T energy relaxation increases strongly with the gas temperature. The same trend is shown in [23]. Over the mean-convective transport time (τ ≈ 33 µs) the gas at T 0 = 1800 K is heated by about 300 K for j = 5 A/cm 2 , due to an incomplete V-T energy relaxation. Note under these specific conditions, the timescale of the V-T energy relaxation of N 2 (X 1 g + , v) molecules on O( 3 P) atoms collisions is τ VT ≈ 86 µs. The major source of gas heating, however, is the V-T relaxation of nitrogen molecules. The E-T energy relaxation does not play a relevant role under these conditions (i.e., η T + η V ≈ 1 [15,34]) The T e values (not shown) softly decreases while the discharge current density (and gas temperature) increases; with values around 9000 K, typical of this kind of low-current discharges in molecular gases [20][21][22][23][24].   Figure 2 shows a comparison between the electron-impact dissociation from all the vibrational levels of N2(X 1 ∑g + , v), with the thermal dissociation for vibrational excited molecules for T0 = 1800 K. As it can be seen, the thermal dissociation, stimulated by vibrational non-equilibrium, dominates over the electron-impact dissociation for j > 2 A/cm 2 ; being the dominant mechanism in the production of N( 4 S) atoms under the conditions considered.  Figure 2 shows a comparison between the electron-impact dissociation from all the vibrational levels of N 2 (X 1 g + , v), with the thermal dissociation for vibrational excited molecules for T 0 = 1800 K.
As it can be seen, the thermal dissociation, stimulated by vibrational non-equilibrium, dominates over the electron-impact dissociation for j > 2 A/cm 2 ; being the dominant mechanism in the production of N( 4 S) atoms under the conditions considered.  Figure 2 shows a comparison between the electron-impact dissociation from all the vibrational levels of N2(X 1 ∑g + , v), with the thermal dissociation for vibrational excited molecules for T0 = 1800 K. As it can be seen, the thermal dissociation, stimulated by vibrational non-equilibrium, dominates over the electron-impact dissociation for j > 2 A/cm 2 ; being the dominant mechanism in the production of N( 4 S) atoms under the conditions considered. The number density of several neutral species versus the discharge current density for T0 = 1800 K is shown in Figure 3. The density of NO molecules is fairly high under the conditions The number density of several neutral species versus the discharge current density for T 0 = 1800 K is shown in Figure 3. The density of NO molecules is fairly high under the conditions considered, due to the high-initial gas temperature value. The rather-high vibrational non-equilibrium state of the discharge for high-values of j promotes the production of N( 4 S) atoms due to thermal dissociation of N 2 (X 1 g + , v) molecules (the N( 4 S) atoms number density being orders of magnitude higher than the corresponding to the local thermal equilibrium at T g [77]), and also leads to a significantly speeds up of the production of N 2 (A 3 u + ) molecules by electron-impact excitation from all the vibrational  Figure 4. The rates of electron production and loss via various mechanisms versus the discharge current density for T 0 = 1800 K is shown in Figure 4. At this low-gas temperature value, the electrons are efficiently lost in reactions of dissociative (R84) and three-body (R83) attachment. These processes, however, are balanced to a great extent by the rapid destruction of negative ions in their interaction with NO molecules via the reaction (R89), and with O( 3 P) atoms for current density values higher than Note that the convective charge transport does not significantly perturbs the local particle balance for j > 1 A/cm 2 , when the time of the electron-ion recombination (which is inversely proportional to the electron number density) is smaller than the convective transport time of the gas in the discharge. The role of the processes involving negative ions is less important for higher values of T 0 , due to the accumulation of components such as NO and O( 3 P), resulting in an increase in the detachment rate.
due to thermal dissociation of N2(X 1 ∑g + , v) molecules (the N( 4 S) atoms number density being orders of magnitude higher than the corresponding to the local thermal equilibrium at Tg [77]), and also leads to a significantly speeds up of the production of N2(A 3 ∑u + ) molecules by electron-impact excitation from all the vibrational levels of N2(X 1 ∑g + , v). The increase of the density of N2(A 3 ∑u + ) molecules, cause in turn, a rise in the concentration of O( 3 P) atoms because O2 molecules are intensively dissociated in the quenching reaction (R44) of N2(A 3 ∑u + ) molecules. Subsequent quenching of N2(A 3 ∑u + ) molecules by N( 4 S) leads in turn to the generation of N( 2 P) metastable atoms trough reactions (R50), while the quenching of N2(A 3 ∑u + ) by O( 3 P) atoms produce both O( 1 S) metastable atoms trough reactions (R48) and N( 2 D) metastable atoms trough reaction (R49). These exited atoms may participate in reactions of associative ionization (R14), (R16) and (R18), as shown in Figure 4. The rates of electron production and loss via various mechanisms versus the discharge current density for T0 = 1800 K is shown in Figure 4. At this low-gas temperature value, the electrons are efficiently lost in reactions of dissociative (R84) and three-body (R83) attachment. These processes, however, are balanced to a great extent by the rapid destruction of negative ions in their interaction with NO molecules via the reaction (R89), and with O( 3 P) atoms for current density values higher than 3 A/cm 2 (see Figure 3), via the reactions (R86) and (R88). The curve labeled as 'effective' attachment represents the difference between the rates of attachment (R83) and (R84), and detachment (R85)-(R90) reactions. The convection due to gas flow as well as the electron attachment are the main channels for the loss of electrons at current density values lower than 1 A/cm 2 . Note that the gas flow has practically no direct effect on the motion of electrons. However, due to the coupling between the electrons and ions through the ambipolar field, the electron trajectories correspond to a field distribution in the discharge that takes into account the removal of ions by the gas flow (note that the flow velocity of the gas (450 m/s) and the drift velocity of ions in the electric field are comparable [15,16]). For larger j values, however, the loss of electrons is dominated by the fast electron-ion recombination reaction (R23). Note that the convective charge transport does not significantly perturbs the local particle balance for j > 1 A/cm 2 , when the time of the electron-ion recombination (which is inversely proportional to the electron number density) is smaller than the convective transport time of the gas in the discharge. The role of the processes involving negative ions is less important for higher values of T0, due to the accumulation of components such as NO and O( 3 P), resulting in an increase in the detachment rate.  It is seen in Figure 4 that the contribution of NO molecules to ionization processes via reaction (R3) is significant because of their relatively low ionization energy (9.27 eV). This result agrees with the inferences made in [23][24][25]. However, the accumulation of N( 2 P) metastable atoms for current density values higher than 3 A/cm 2 , significantly speeds up the ionization kinetics of the discharge; mainly via the following reactions: the rate coefficient of reaction (R18) is weakly gas temperature dependent (it has no activation barrier), and independent of the reduced electric field. In particular, the charged particle generation It is seen in Figure 4 that the contribution of NO molecules to ionization processes via reaction (R3) is significant because of their relatively low ionization energy (9.27 eV). This result agrees with the inferences made in [23][24][25]. However, the accumulation of N( 2 P) metastable atoms for current density values higher than 3 A/cm 2 , significantly speeds up the ionization kinetics of the discharge; mainly via the following reactions: e + N 2 (X) → e + N 2 (A) (R9) N 2 (A) + N 4 S → N 2 (X) + N 2 P (R50) the rate coefficient of reaction (R18) is weakly gas temperature dependent (it has no activation barrier), and independent of the reduced electric field. In particular, the charged particle generation is controlled predominantly by reaction (R18) for j > 5 A/cm 2 . For higher values of the initial gas temperature of the discharge T 0 (not shown), the following reactions become also important in the gas ionization: because the rate coefficient of reaction (R16) is strongly gas temperature dependent (it has a low-activation barrier of 0.38 eV), and independent of the reduced electric field. The associative ionization reaction (R13) (with a high-activation barrier of 2.76 eV), does not play a significant role under the considered conditions due to the low-gas temperature of the discharge. Figure 5 shows the number density of the charged particles in the discharge versus the discharge current density for T 0 = 1800 K. For j = 2 A/cm 2 the electron number density is N e ≈ 2 × 10 18 m -3 , in good agreement with the estimate in [23,24]. It is also observed that the electron number density increases by one order of magnitude when the applied current density is increased by one order of magnitude. This trend is consistent with the estimations in [13]. The total number density of negative ions decreases with j, due to the accumulation of O( 3 P) atoms resulting in an increase in the detachment rate.  The reduced electric field (1 Td  10 −21 V m 2 ) versus the discharge current density is shown in Figure 6, both measured in [9,13] and obtained by simulation for T0 = 1800 K. Curves are the results of calculations, with (solid line) and without (dashed line) associative ionization reactions (R14)-(R16) and (R18). The reduced electric field given by the solid line follows the changes in the ionization mechanisms shown in Figure 4. It begins to appreciably decrease at current density values higher than 3 A/cm 2 (i.e., when the reaction (R18) begins to significantly contribute to the ionization kinetics of the discharge) because the electron energy (6.2 eV) required for the excitation of the N2(A 3 ∑u + ) state is considerably lower than the ionization energy (9.27 eV) of the NO molecules. The calculated reduced electric field values considering the ionization reactions in metastable atomic collisions are seen to be in fairly good agreement with the experimental data over one order of magnitude in current density. The reduced electric field (1 Td ≡ 10 −21 V m 2 ) versus the discharge current density is shown in Figure 6, both measured in [9,13] and obtained by simulation for T 0 = 1800 K. Curves are the results of calculations, with (solid line) and without (dashed line) associative ionization reactions (R14)-(R16) and (R18). The reduced electric field given by the solid line follows the changes in the ionization mechanisms shown in Figure 4. It begins to appreciably decrease at current density values higher than 3 A/cm 2 (i.e., when the reaction (R18) begins to significantly contribute to the ionization kinetics of the discharge) because the electron energy (6.2 eV) required for the excitation of the N 2 (A 3 u + ) state is considerably lower than the ionization energy (9.27 eV) of the NO molecules. The calculated reduced electric field values considering the ionization reactions in metastable atomic collisions are seen to be in fairly good agreement with the experimental data over one order of magnitude in current density. Figure 7 shows the average electric field and the reduced electric field versus the initial gas temperature of the discharge T 0 = 1800-2900 K, for a current density j = 1 A/cm 2 . The results of calculations are in good agreement with the experimental data [13]. Note that the reduced electric field depends weakly on the gas temperature: as T 0 increases from 1800 to 2900 K, the reduced field decreases from 46 to 41 Td. The decrease is caused by the formation of NO molecules with low-ionization energy (note from Figure 4 that at j = 1 A/cm 2 , the main ionization process is via the reaction (R3)). This agrees with the results of calculations performed in [24,25].
(R16) and (R18). The reduced electric field given by the solid line follows the changes in the ionization mechanisms shown in Figure 4. It begins to appreciably decrease at current density values higher than 3 A/cm 2 (i.e., when the reaction (R18) begins to significantly contribute to the ionization kinetics of the discharge) because the electron energy (6.2 eV) required for the excitation of the N2(A 3 ∑u + ) state is considerably lower than the ionization energy (9.27 eV) of the NO molecules. The calculated reduced electric field values considering the ionization reactions in metastable atomic collisions are seen to be in fairly good agreement with the experimental data over one order of magnitude in current density.  Figure 7 shows the average electric field and the reduced electric field versus the initial gas temperature of the discharge T0 = 1800-2900 K, for a current density j = 1 A/cm 2 . The results of calculations are in good agreement with the experimental data [13]. Note that the reduced electric field depends weakly on the gas temperature: as T0 increases from 1800 to 2900 K, the reduced field

1.
A kinetic scheme for non-equilibrium regimes of atmospheric pressure air discharges is developed. An improvement of the model is that it considers associative ionization with the participation of N( 2 D, 2 P) exited atoms. 2. The model is used to simulate the parameters of a glow discharge ignited in a fast longitudinal flow of preheated (T0 = 1800-2900 K) air. The results adequately describe the dependence of the electric field in the glow discharge on the initial gas temperature. 3. The rather-high vibrational non-equilibrium state of the discharge for high current density values, promotes the production of N( 4 S) atoms due to thermal dissociation of N2(X 1 ∑g + , v) molecules, and also leads to a significantly speeds up of the production of N2(A 3 ∑u + ) molecules by electron-impact excitation from all the vibrational levels of N2(X 1 ∑g + , v). 4. For T0 = 1800 K, the accumulation of N( 2 P) metastable atoms at current density values higher than 3 A/cm 2 , significantly speed up the ionization kinetics of the discharge; mainly via the following reactions: Figure 7. Electric field strength (solid line) and reduced electric field (dashed line) as functions of the initial gas temperature of the discharge T 0 for j = 1 A/cm 2 .

1.
A kinetic scheme for non-equilibrium regimes of atmospheric pressure air discharges is developed.
An improvement of the model is that it considers associative ionization with the participation of N( 2 D, 2 P) exited atoms. 2.
The model is used to simulate the parameters of a glow discharge ignited in a fast longitudinal flow of preheated (T 0 = 1800-2900 K) air. The results adequately describe the dependence of the electric field in the glow discharge on the initial gas temperature. 3.
The rather-high vibrational non-equilibrium state of the discharge for high current density values, promotes the production of N( 4 S) atoms due to thermal dissociation of N 2 (X 1 g + , v) molecules, and also leads to a significantly speeds up of the production of N 2 (A 3 u + ) molecules by electron-impact excitation from all the vibrational levels of N 2 (X 1 g + , v).

4.
For T 0 = 1800 K, the accumulation of N( 2 P) metastable atoms at current density values higher than 3 A/cm 2 , significantly speed up the ionization kinetics of the discharge; mainly via the following reactions: e + N 2 X 1 Σ + g , v → e + N 2 A 3 Σ + u , N 2 A 3 Σ + u + N 4 S → N 2 X 1 Σ + g , v + N 2 P , O 3 P + N 2 P → NO + + e, accordingly; the reduced electric field begins to appreciably decrease at current densities higher than 3 A/cm 2 , because the electron energy (6.2 eV) required for the excitation of the N 2 (A 3 u + ) state is considerably lower than the ionization energy (9.27 eV) of the NO molecules. For higher values of T 0 , the following reactions become also important in the charged particles production: Funding: This research was funded by Universidad Tecnológica Nacional, grants 5418 and 5447, and Agencia Nacional de Promoción Científica y Tecnológica, grant PICT 2015-1 553.