Impact of Initial Electron Development on Sodium Plasma Generation Using LIBORS Technique : Numerical Modeling

The present work reports an investigation on the role played by Na ions formed through triatomic associative ionization collision of Na(4d) atoms with Na2 ground state molecules during the early phase of sodium plasma generation by laser ionization based on resonance saturation (LIBORS). Such ionization mechanism is observed experimentally for the first time by Tapalian and Smith (1993) [1]. In their experiment, stepwise atomic excitations are created using two CW dye lasers; one laser is tuned to 589 nm to excite the Na(3s) to Na(3p) D2 transition of sodium and the other laser is tuned 569 nm to excite the Na(3p) to Na(4d) transition. The analysis is grounded on a numerical study of the role of seed electron processes on the temporal evolution of sodium plasma formation by laser irradiation. A previously developed numerical model based on LIBORS technique is modified and adopted. In the present study, the sodium atom is treated as an atom comprises 22 levels namely: a ground state, 18 excited states and three ionic states (atomic, molecular and tri-atomic). The model tackled various collisional and radiative processes that act to enhance and deplete the free electrons generated in the interaction region. The contribution of these processes is signified by solving numerically a system of time-dependent rate equations, which couple the generated atomic and ionic species with the laser fields. Meanwhile, it solves the time-dependent Boltzmann equation for the electron energy distribution function (EEDF) of the generated electrons. The computed values of the EEDF, time evolution of both excited states population and the formed ionic species considering the individual effect of associative ionization, Penning, and photo-ionization and triatomic associative ionization justified the important effect of each of these ionizing processes in creating the early stage electrons. These seed electrons are assumed to rapidly gain energy through superelastic collisions leading eventually to plasma development. How to cite this paper: Hamam, K.A. and Gamal, Y.E.E.-D. (2018) Impact of Initial Electron Development on Sodium Plasma Generation Using LIBORS Technique: Numerical Modeling. Journal of Modern Physics, 9, 669-684. https://doi.org/10.4236/jmp.2018.94046 Received: February 12, 2018 Accepted: March 27, 2018 Published: March 30, 2018 Copyright © 2018 by authors and Scientific Research Publishing Inc. This work is licensed under the Creative Commons Attribution International License (CC BY 4.0). http://creativecommons.org/licenses/by/4.0/ Open Access DOI: 10.4236/jmp.2018.94046 Mar. 30, 2018 669 Journal of Modern Physics K. A. Hamam, Y. E. E.-D. Gamal


Introduction
Study of laser-produced plasmas is one of the rapidly growing fields of present-day physics.The plasma generated by intense laser pulses in gases is a very well recognized phenomenon.The productivity of this process was found to be increased by many orders of magnitude when the laser wavelength is tuned to a selected transition of alkali vapors leading to fast ionization and plasma generation.Such phenomenon involves many excitation and ionization collisions besides radiative transfer processes.
The pioneering observation of this phenomenon was given by McIlrath and Lucatorto [2] when measurements are performed with laser radiation tuned to the first resonance line in dense sodium vapor.These measurements included some experimental study of the collisional mechanisms responsible for complete ionization in a sodium vapor illuminated by a CW laser radiation of wavelength λ = 589.0nm tuned to the 3s-3p resonance line.Ever since efficient ionization by resonance radiation has been tested in different alkali-metal vapors [3] [4] [5] [6].The innovative work given in Ref. 2, has been extended by various experimental [7]- [12], and theoretical [13]- [19] studies.All aimed to investigate the excitation energy transfer and collisional ionization processes in a laser-excited alkali vapor.The investigation offered good evidence for their importance in radiative transfers involved different subjects such as: in gas lasers, astrophysics, controlled thermonuclear plasma, plasma physics, photochemistry [20] [21] as well as in creating plasma guide channels [22].
The results of McIlrath and Lucatorto [2] were interpreted by Measures et al., [22] [23], on the basis of superelastic heating of electrons in colliding with the excited atoms (collisions of the second kind).The theoretical model that describes the formation of resonance plasma is referred to as LIBORS.It proceeds as follows: initially, the CW laser field saturates the atomic resonance transition.
Then free electrons are produced in the medium due to associative and Penning ionization processes [12].These electrons gain energy in superelastic collisions by deactivating excited atoms, with the subsequent production of an electron avalanche resulting in complete ionization of the gas.
On the other hand, ionization of alkali atoms by nanosecond resonant laser pulses has played a strong part in coupling energy into vapor.Thus the feedback of resonant irradiation of dense sodium vapor, by nanosecond pulsed laser, has been the subject of some investigations by different researchers [24] [25] [26].
These studies revealed that enormous atomic, molecular, and photon interactions could take place in the exciting medium.Of specific importance are the io-nization mechanisms and energy transfer processes which are acting to enhance the electrons density leading eventually to plasma formation.Theoretical modeling of ionization of sodium vapor by nanosecond resonant laser pulses tuned to the 3s-4p transition, is presented by Stacewicz and Latek [27].In this study, the laser pulse saturates this selected optical transition and produces a small number of initial electrons.Then the avalanche ionization develops without the presence of the laser radiation, using the energy stored in the excited atoms.
On the other hand, in (1993) Tapalian and Smith [1] observed experimentally for the first time the formation of Na 3+ ions.These ions are produced through triatomic associative ionization interaction between Na(4d) atoms and Na 2 ground state molecules (which are assumed to be present at ≈ 0.5% percent in the atomic beam), during sodium plasma generation by laser ionization based on resonance saturation (LIBORS), according to the reaction ( ) Na 4d Na Na e The manifestation of this process requires a high population of the Na (4d) state.Therefore, in these measurements, stepwise atomic excitations are created using two CW dye lasers.The first laser is tuned to 589 nm to excite the Na(3s) to Na(3p) D2 transition of sodium and the second laser is tuned to 569 nm to excite the formed Na(3p) to the Na(4d) state.Their results showed the interesting importance of this process where they suggested that the collision cross section for this process is of the order of the cross-section for the atomic associative ionization collision of Na(4d) atoms with Na(3s) atoms.
This associative ionization reaction mechanism presents several possible applications.One of these applications is the study of cluster formation in atomic beams, where the formed  [28] to interpret the experimental observation which is executed to study the associative ionization, and Penning ionization processes (atom-atom collisions).
These are considered as main sources for seed electrons development in the ionization of sodium vapor irradiated by a CW laser source tuned to the transition 3s-3p prearranged by Caree el al. [29].Then, this model is modified [30] and could successfully investigate the finding of the experimental observation given in [1] that carried out to measure to the first time the formation of tri-atomic ions in laser ionization of sodium vapor.
As an extension to this work, the present study examines the ionization of sodium vapor irradiated by CW laser radiation taken into account the occurrence of various seed electron generation processes.The ionization mechanism proceeds as follows: 1) heating of the generated seed electrons through superelastic collisions, 2) initiation of electron impact ionization by the heated electrons leading to the development of new electrons.These electrons are then combined with steps (1) and (2) lead to ionization burnout which results in an almost complete ionization.These processes are compensated by loss of electrons from the interaction region through radiative and three body recombination processes.

Basic Formalism
A detailed study of the theoretical background and modeling of ionization based on resonance saturation using tuned CW laser source is presented in our previous paper [28].Sodium vapor has ionization energy, E ion , greater than twice E 12 , which is the energy between the ground state (3s) and the first excited state (the saturation level, 3p), and less than the three times that energy, i.e., E ion < 3.0 x E 12 .This is a property of the alkali metal elements.Here, we recall that the modified LIBORS model [30] is applied to analyze plasma formation induced by pulsed laser radiation.The sodium atom is presumed to comprising 22 levels namely; a ground state, 18 excited states and three ionic states (atomic, molecular, and tri-atomic).Taken these assumption into consideration the following subsequent chain of ionization events take place in the interaction region: 1) two-photon ionization of the 3p state, 2) single photoionization of the highly excited nl states 3) associative ionization of the 3p-3p states, 4) associative ionization of the 3p-3d states, 5) Molnar-Hornbeck ionization in a collision between ground state and any of the highly excited states, 6) Penning ionization in a collision between the 3p state and any of the highly excited states, 7) laser-induced Penning ionization a collision between two 3p states in the presence of the laser photon, 8) electron impact excitation of electrons having energy ε > E x (E x is the excitation energy from the ground or any of the excited states), 9) super-elastic collisions between highly excited states and low energy electrons ε 0 < ε, 10) Collisional ionization by electrons of energy ε > E ion , (where ε is the total energy of the electron and is assumed to equal ε = (ε 1 ) + (ε 2 ) + E ion where ε 1 and ε 2 are the energies of the impinging and ejected electrons respectively), 11) energy pooling collisions between two 3p states, 12) radiative recombination, 13) three-body recombination.To simplify the calculations, we have assumed that in the electron impact ionization and its inverse processes three-body recombination, the impinging and ejected electrons are assumed to carry the same amount of energy (i.e.) ε 1 = ε 2 .

Rate Equations
The above mentioned physical processes are encountered in the applied model given in Ref. 30 This model solves iteratively a set of rate equations that represent the time variation of the instantaneous population density of the various species present in the formed plasma.Besides, the model solves the Boltzmann equation for the electron energy distribution for the given normalization conditions.These are written as follows: The rate of growth of the molecular, atomic, and the tri-atomic ions is given by: The time evolution of the electron energy distribution function is expressed by Boltzmann equation including all the processes which involve electrons interactions and written as: The normalization conditions considered in this analysis are taken as: ( ) ( ) N(3s), N(3p) and N(n) represent the population density of the ground, first excited (saturated), and highly excited states respectively.n e (ε) refers to the free electron density of energy ε.N 0 is the density of Na vapor.F represents the photon flux density.v is average velocity of atoms, in cm/sec.Note that the factor 1/2 appeared in the terms express the energy pooling and associative ionization mechanisms are added to correct for possible double counting of each colliding pair of identical particles [31].
The definition and source of the rate coefficients and cross sections of the collisional and photoionization processes appeared in the rate Equations (2) to (8) are defined in Table 1.

Results and Discussions
Equations ( 2) to (10) are incorporated into a computer program together with the relations which describe the rate coefficients and cross-sections of the physical processes encountered into the model and solved numerically applying Runge-Kutta fourth order technique.The calculations are performed under the experimental conditions given by Taplain and Smith [1].In their experiment, the generation of sodium plasma induced by tuning a CW broadband dye laser to 589 nm to excite the Na3s ( 2 S 1/2 ) to Na3p ( 2 P 3/2 ) transition and a CW frequency-stabilized narrow-band dye laser is tuned to 569 nm to excite the Na3p ( 2 P 3/2 ) to Na4d ( 2 D 5/2 ).The applied laser beams of intensities up to 100 W/cm 2 are collimated and directed to the Na vapor chamber.The analysis considered atomic vapor density as 4.0 × 10 12 cm −3 at temperature 500 -600 K for exposure time vary from 100 up to 500 ns.The irradiated medium is found to evolve free electrons as well as atomic, molecular and trimer ions.Consequently, the study is devoted to examining the individual contribution of each process responsible for generating the seed electron and hence positive ions namely: associative ionization, Penning ionization, photo-ionization, and tri-atomic associative Table 1.Definition of the rate coefficients and cross sections appeared in the rate equations.R 21 (sec −1 ) Stimulated emission rate coefficient for transition from level 2 to 1 [24] R 26 (sec −1 ) Stimulated emission rate coefficient for transition from level 2 to 6 [32] K EP Energy pooling collisions [11] K mn (cm 3 sec −1 ) Electron-atom collisional excitation or de-excitation between the states m and n [17] K nc (cm 3 sec −1 ) Electron collisional ionization rate coefficient from any energy level (n) [6] K RD (cm 3 sec −1 ) Radiative recombination rate coefficient to any energy level (n) [33] K cn (cm 6 sec −1 ) Three body recombination rate coefficient [33] A 21 , A nm (n > m) Einstein coefficient for the spontaneous emission from the resonance level (2 → 1) and transition from any level n > m (n → m) [34] K AI1 and K AI2 , Associative ionization for the reactions (3p-3p) [15], and (3p-3d) [16] K MHI Hornbeck-Molnar ionization [35] K PI Penning ionization [29] ( ) Single-photon ionization cross section from any high energy level n [36] ( ) σ Two-photon ionization cross section from the saturated energy level (3P) [36] PL σ Laser-induced Penning ionization cross section [20] 2 4 Na d K + Triatomic associative ionization rate coefficient [1] ionization to plasma formation.This investigation is obtained via a comparative study of the EEDF as well as the temporal evolution of electrons and ions density during the exposure time in the presence and absence of each of these ionization processes.

Associative Ionization
To examine the exact contribution of the associative ionization (AI) to the growth rate of the seed electrons calculations are carried out to obtain the EEDF in the presence and absence of this process for sodium atomic density 4 × 10 15 cm −3 and exposure time 500 ns.This is displayed in Figure 1 as curves (1) and (2) respectively.This figure demonstrates that the structure of the electron spectra has a non-Maxwellian distribution.Moreover, the effect of introducing the AI process (curve 1) appeared only in the low energy region (0.0 -0.6 eV) where an increase is observed in the peak labeled (a) at an electron energy of about 0.44 eV.This increase reflects the role played by the AI process (produced by a collision of two Na3p atoms) in enriching the density of the low energy electrons.
Beyond this energy range (~0.6 eV up to 5.25 eV) however, the structure of the electron spectra represented by both curves (1) and ( 2) coincides.
To clarify this result the temporal variation of the density of electrons corresponding to the low energy region (0.0 < ε < 0.6 eV) in the presence (curve 1) and absence (curve 2) of AI is calculated and shown in Figure 2(a).It is noticed here that curve (1) lies much above the curve (2) over the whole exposure time duration.Also, the electron density grows at a constant rate.The obtained growth  rate is found to be ~5.17

Penning Ionization and Photoionization
In studying the effect of Penning ionization (PI) on the production of seed electrons, calculations showed that this process interplays with the photo-ionization process in line with the applied experimental conditions.This result could be attributed to the fact that the occurrence of both processes depends on the population of atoms in the highly excited states (5s, 4d/4f, 6s, and 5d).These states are mostly populated through energy pooling collisions of two Na(3p) atoms.
While the 4p and 3d states are rather populated by radiative decay from these upper states which are produced by collisions [21] [37].So, to account for the actual contribution of the PI process, a study is performed to compare the structure of the EEDF in the absence and presence of photo-ionization.This is displayed in Figure 3(a), and Figure 3(b) respectively.In this figure, the peaks appeared in the structure of the electron energy spectrum reflect the contribution of a specified physical process as shown in Table 2. ground state atoms or ionize the populated excited states.These processes lead to the generation of low energy electrons.This, in turn, could elucidate the increase of the EEDF corresponds to the peak (a) in the presence of PI process.
On the other hand, to clarify the role played by the PHI process on the plasma development Figure 3(b) illustrates the calculated EEDF in the presence of this process under the action and omission of the PI process.This figure showed an increase of the values of the peaks (a to f) with a shift of the whole electron energy spectrum towards the high energy regime in the presence and absence of PI as displayed by curves (1), and (2) respectively.This result is depicted in Table 3.Moreover, It is shown from this figure that, despite the growth of the values of the minima appeared at energies ~0.75 eV and ~2.75 eV, there is a clear similarity in the structure of the EEDF represented by the curves (1 and 2).This result is investigated following to the fact that, PI acts only to enhance the generation of the electrons density carrying energy in the range 0.75 eV < ε < 1.25 eV.These electrons are subsequently heated up through super-elastic collision (SEC) with Na(3p) atoms.Those excited atoms boost the kinetic energy of the created electrons by 2.1 eV.This, in turn, enhances the population of electrons having the energy of about ~2.75 e V.Over the remaining energy range (2.75 eV < ε < 5.0 eV) both curves (1, and 2) are almost found to coincide.The behavior displayed in this figure is demonstrated in Table 3 which undertaking the pronounced effect of the PHI process on the rate of electron production and its interplay with the PI process.Bearing in mind that, the occurrence of both mechanisms mainly depends on the population of the same highly excited states.This result is assured from the fast growth of the EEDF obtained at the first peak (a) under the effect of PI and absence of PHI given in Table 3.Consequently, in dealing with PI process, one should take into account the competing role between PI and photo-ionization processes in relation to the experimental condition under investigation.
As we have seen above, the mutual effect of photo-ionization and PI processes leads to a change in the electron energy structure.Accordingly, it was necessary to have a detailed study on the real impact of the photoionization process on the

Tri-Atomic Associative Ionization
According to what mentioned above, this process results in the interaction between atoms in the 4d state and ground state molecules and leads to the formation of triatomic ions and free electrons.So, for this process to take place the Na(4d) state should be highly populated.This is the role played by the 569 nm laser source which is tuned to the transition 3p-4d in the experimental measurement under investigation [1]. Figure 6

Conclusions
In -Photoionization showed a pronounced effect on the enhancement of the atomic ions density during the examined exposure time range.Besides, the peaks of the EEDF showed an observable increase, and their position undergoes slight shifting towards the high-energy region.This shift points out the increase in the density of the high-energy electrons.
-The formation of Tri-atomic associative ionization is found to play a vital role in enhancing the seed electron population which in turn accelerates the plasma generation (less exposure time).Besides, it increases the possibility of studying cluster formation in atomic beams, for example, the formed These clusters may found several applications.

Figure 1 .
Figure 1.Comparison between the EEDF calculated in the presence (curve 1) and absence (curve 2) of associative ionization for the Na density N 0 = 4 × 10 15 cm −3 and exposure time 500 ns.

Figure 2 .
Figure 2. Temporal evolution of (a) density of electrons generated in the energy range 0.0 < ε < 0.6 eV, and (b) molecular ions density for Na atomic density 4.0 × 10 +15 cm −3 in the presence (curve 1) and absence (curve 2) of AI.
× 10 3 at t = 500 ns, which corresponds to an increase of the low energy electron density from ~1.4 × 10 7 cm −3 to ~8.0 × 10 10 cm −3 .Figure 2(b) illustrates similar result but for the time variation of the formed molecular ions through AI process.Here curve (1) represents the presence of AI, and curve (2) is calculated in its absence.The similarity is shown in Figure 2(a) and Figure 2(b) for the time evolution of; electrons, and molecular ions respectively assured that AI is mainly responsible for the generation of electrons in the low energy range (0.0 -0.6 eV).Moreover, the high growth rate of these low energy electrons represented by curve (1) elucidated the significant role played by AI process in accelerating plasma generation.

Figure 3 (Figure 3 .
Figure 3(a) (absence of photoionization) shows that the structure of the calculated EEDF under the effect of PI (curve 1) is summarized as follows: decrease at the minimum position corresponds to electron energies 1.7 eV and 3.6 eV (PI and the first SEC of electrons generated by PI process respectively).An increase is shown at the peak (a) at energy ~0.28 eV (AI) followed by an increase in electron energies 0.5 eV, 2.4 eV, and 4.4 eV.This performance of the EEDF can be explained based on the fact that electrons generated by AI (0.28 eV) are heated up through the first (2.4 eV) and second (4.45 eV) SECs.Hence, they can undergo inelastic collisions and lose their energies in either exciting

Figure 4 andFigure 4 .
Figure 4. Time evolution of electron density in the presence (curve 1) and absence (curve 2) of the photoionization process.

Figure 5 .
Figure 5. Same as in Figure 4 but for the atomic ions density.

Figure 6 .
Figure 6.Electron energy distribution function calculated for atomic sodium density of 4 × 10 15 cm −3 and exposure time 500 ns in the presence (curve 1) and absence (curve 2) of tri-atomic associative ionization process.
the present work, using a modified model based on LIBORS technique, the results of calculations of the EEDF in the presence and absence of each of the seed electron generation processes namely: Associative ionization, Penning ionization, photoionization, and Tri-atomic associative ionization showed the non-Maxwellian distribution of the electrons.Besides, the growth rate of electrons which reflects the ions formation is achieved by comparing the electron (ions) density in the presence and absence of each of the considered ionization processes.This leads to the following conclusion: -Associative ionization results only in an increase of the electron density by an order of ≈ 10 4 on the energy range 0.0 -0.25 eV.This, in turn, represents a fair contribution to plasma development.The similar growth rate is obtained for the corresponding molecular ions density.-Penningionization and photoionization are found to signify two strongly competing processes both sharing the same highly excited states which are populated via energy pooling process.The later process diminishes the contribution of the Penning ionization to the plasma generation.Consequently, in studying the role of Penning ionization process, photoionization process is ignored.

3
Na + , may react with another Na atom (ground state or excited) to form 4 Na + .

Table 3 .
Peak label, peak energy position, EEDF calculated under the effect and ignorance of PI for the absence (a) and presence (b) of photoionization process.