<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">JMP</journal-id><journal-title-group><journal-title>Journal of Modern Physics</journal-title></journal-title-group><issn pub-type="epub">2153-1196</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/jmp.2019.107050</article-id><article-id pub-id-type="publisher-id">JMP-93013</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Comparison between Theoretical and Experimental Radial Electron Temperature Profiles in a Low Density Weakly Ionized Plasma
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Ahmed</surname><given-names>Rida Galaly</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Guido</surname><given-names>Van Oost</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Department of Physics, Faculty of Science, Beni-Suef University, Beni Suef, Egypt</addr-line></aff><aff id="aff2"><addr-line>Department of Engineering Science, Faculty of Community, Umm Al_Qura University, Mecca, Mecca Province, KSA</addr-line></aff><pub-date pub-type="epub"><day>05</day><month>06</month><year>2019</year></pub-date><volume>10</volume><issue>07</issue><fpage>699</fpage><lpage>716</lpage><history><date date-type="received"><day>15,</day>	<month>May</month>	<year>2019</year></date><date date-type="rev-recd"><day>11,</day>	<month>June</month>	<year>2019</year>	</date><date date-type="accepted"><day>14,</day>	<month>June</month>	<year>2019</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  Experimental and theoretical studies of the radial distribution function of the electron temperature (RDFT) in a low-density plasma and weakly ionized gas for the abnormal glow region are presented. Experimentally, the electron temperatures and densities are measured by a Langmuir probe moved radially from the center to the edge of the cathode electrode for helium gas at different pressures in the low-pressure glow discharge. The comparison of the final experimental data for the radial distribution of electron temperatures and densities for different low pressures ranging from 0.2 to 1.2 torr, with the final proved equation of RDFT confirms that the electron temperatures decrease with increasing product of radial distance and gas pressures, showing a radial decrement dependence of the electron temperature from the center to the edge of the electrode. This is attributed to the increase of the number of electron-atom collisions at higher gas pressures and consequently of the rate of ionization. For the axial distance (
  <em>L</em>) from the tip of the probe to cathode electrode and the cathode electrode radius (
  <em>R</em>), a theoretical and experimental comparison for the two conditions 
  <em>L</em> &lt; 
  <em>R</em> and 
  <em>L </em>&gt; 
  <em>R</em>, for both cases the produced plasma temperatures decrease and densities increase. It is concluded that the RDFT accurately shows a dramatic decrease for 
  <em>L</em> &lt; 
  <em>R</em> by 60% less than RDFT values for 
  <em>L</em> &gt; 
  <em>R</em> similar as for conditions of magnetized and unmagnetized effect for DC plasma. This means that the rate of plasma loss by diffusion decreased for 
  <em>L</em> &lt; 
  <em>R</em>, agrees well with the applied of magnetic field behavior
 
</p></abstract><kwd-group><kwd>Radial Distribution</kwd><kwd> Electron Temperature</kwd><kwd> Ambipolar Diffusion</kwd><kwd> Low Pressure</kwd><kwd> Weakly Ionized Plasma</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>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<sub>ri</sub>) [<xref ref-type="bibr" rid="scirp.93013-ref1">1</xref>] [<xref ref-type="bibr" rid="scirp.93013-ref2">2</xref>] [<xref ref-type="bibr" rid="scirp.93013-ref3">3</xref>] .</p><p>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 [<xref ref-type="bibr" rid="scirp.93013-ref4">4</xref>] , the plasma parameters especially electron temperature and electron density can be determined from the current-voltage characteristic curves of the probe [<xref ref-type="bibr" rid="scirp.93013-ref5">5</xref>] . Diffusion phenomena occur in a plasma when a spatial gradient of the charged species is present [<xref ref-type="bibr" rid="scirp.93013-ref6">6</xref>] . Furthermore, it is caused by the differently charged species having different diffusivities, hence there is a loss of neutrality [<xref ref-type="bibr" rid="scirp.93013-ref7">7</xref>] 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 [<xref ref-type="bibr" rid="scirp.93013-ref8">8</xref>] .</p><p>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 [<xref ref-type="bibr" rid="scirp.93013-ref9">9</xref>] [<xref ref-type="bibr" rid="scirp.93013-ref10">10</xref>] , 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.</p><p>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.</p></sec><sec id="s2"><title>2. Experimental Study</title><sec id="s2_1"><title>2.1. Experimental Set-Up</title><p><xref ref-type="fig" rid="fig1">Figure 1</xref> 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 &gt; L and R &lt; 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 <xref ref-type="table" rid="table1">Table 1</xref>.</p><p><xref ref-type="fig" rid="fig1">Figure 1</xref> 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<sub>P</sub> is applied to the probe.</p></sec><sec id="s2_2"><title>2.2. Experimental Study of Axial Potential and Electric Field Distribution Measurements</title><p>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 <xref ref-type="fig" rid="fig2">Figure 2</xref>, shows that the potential</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Shows the values of operating parameters</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Parameters</th><th align="center" valign="middle" >Values</th></tr></thead><tr><td align="center" valign="middle" >Discharge Current (I<sub>a</sub>)</td><td align="center" valign="middle" >4 - 30 mA</td></tr><tr><td align="center" valign="middle" >Gas Pressure (P)</td><td align="center" valign="middle" >0.1 - 1.5 torr</td></tr><tr><td align="center" valign="middle" >Discharge Voltage (V<sub>a</sub>)</td><td align="center" valign="middle" >200 - 1200 V DC</td></tr><tr><td align="center" valign="middle" >Current Density</td><td align="center" valign="middle" >2 - 15 mA/m<sup>2</sup></td></tr><tr><td align="center" valign="middle" >Working Gas</td><td align="center" valign="middle" >He</td></tr></tbody></table></table-wrap><p>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.</p><p>Values of the electric field distribution are obtained by differentiating the measured potential distribution of <xref ref-type="fig" rid="fig2">Figure 2</xref>. <xref ref-type="fig" rid="fig3">Figure 3</xref> and <xref ref-type="fig" rid="fig4">Figure 4</xref> 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</p><p>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.</p></sec><sec id="s2_3"><title>2.3. Experimental Study of Radial Dependence of the Electron Temperature</title><p>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.</p><p>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 [<xref ref-type="bibr" rid="scirp.93013-ref10">10</xref>] [<xref ref-type="bibr" rid="scirp.93013-ref11">11</xref>] .</p><p><xref ref-type="fig" rid="fig3">Figure 3</xref> and <xref ref-type="fig" rid="fig4">Figure 4</xref> 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 [<xref ref-type="bibr" rid="scirp.93013-ref12">12</xref>] . 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. <xref ref-type="fig" rid="fig5">Figure 5</xref> and <xref ref-type="fig" rid="fig6">Figure 6</xref> show the radial temperatures and densities distributions at different pressures from center to the edge. T<sub>e</sub> is decreased and n<sub>e</sub> increased due to the general trend that values of T<sub>e</sub> and n<sub>e</sub><sub> </sub>are inversely proportional [<xref ref-type="bibr" rid="scirp.93013-ref13">13</xref>] .</p><p>Increasing the helium working pressures from 0.2 to 1.2 torr, the temperature T<sub>e</sub> 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<sub>e</sub> increased at center n<sub>e</sub> from 0.6 &#215; 10<sup>9</sup> cm<sup>−3</sup> to 1.3 &#215; 10<sup>9</sup> cm<sup>−3</sup> and at the edge from 2.5 to 3.95 &#215; 10<sup>9</sup> cm<sup>−3</sup>.</p></sec></sec><sec id="s3"><title>3. Theoretical Consideration of the Radial Dependence of the Electron Temperature</title><p>Conductivity at the tube axis, and charge density of plasma according to the following expression [<xref ref-type="bibr" rid="scirp.93013-ref14">14</xref>] [<xref ref-type="bibr" rid="scirp.93013-ref15">15</xref>]</p><p>D s = D a [ 1 − μ e ρ σ ] (1)</p><p>where D<sub>a</sub> is the ambipolar diffusion coefficient as follows:</p><p>The flow of ions and electrons are the same, hence</p><p>Γ = Γ i = Γ e</p><p>(Congruence approximation). If there is no external electric field, the fluxes of ions and electrons (drift-diffusion model) can be expressed as:</p><p>Γ i = − D i ∇ n i + μ i n i E D</p><p>Γ e = − D e ∇ n e + μ e n e E D</p><p>where E<sub>D</sub> follows Poisson’s equation. If we multiply the first equation by μ e n e and the second one by μ i n i and then subtract them, we get:</p><p>μ i n i Γ e − μ e n e Γ i = μ e n e D i ∇ n i − μ e n e μ i n i E D − μ i n i D e ∇ n e + μ i n i μ e n e E D</p><p>Since Γ = Γ i = Γ e , the terms with E<sub>D</sub> cancel each other out and using n<sub>e</sub> = n<sub>i</sub>, the expression for the flow of particles becomes:</p><p>Γ = − ( D e μ i − D i μ e ) μ i − μ e ∇ n i</p><p>With the ambipolar diffusion coefficient being</p><p>D a = D e μ + − D + μ e μ + − μ e (2)</p><p>and ρ is the charge density in the plasma given by:</p><p>ρ = e [ n + − n e ] (3)</p><p>and σ is the conductivity at the tube axis given by:</p><p>σ = e [ μ + n + − μ e n e ] (4)</p><p>Considering that diffusion coefficient D given by:</p><p>D = k T m ν   m / s 2 (5)</p><p>where μ ε , + are the mobility of electrons (e) and ions (+), kT (ev) is the particle temperature, ν 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</p><p>D s = D e μ + − D + μ e μ + − μ e [ 1 − μ e n + − n e μ + n + − μ e n e ] (6)</p><p>But by substituting by the Einstein relation [<xref ref-type="bibr" rid="scirp.93013-ref16">16</xref>] [<xref ref-type="bibr" rid="scirp.93013-ref17">17</xref>] :</p><p>μ = ⌈ q ⌉ D k T (7)</p><p>into Equation (6) by the value of D, then we get</p><p>D s = μ e μ + K T e − μ e μ + K T + q ( μ + − μ e ) [ 1 − μ e n + − n e μ + n + − μ e n e ] (8)</p><p>or</p><p>D s = μ e μ + ( K T e − K T + ) e ( μ + − μ e ) [ μ + n + − μ e n e − μ e n + + μ e n e μ + n + − μ e n e ] (9)</p><p>Neglecting k T + , where k T e ≫ k T + finally results in</p><p>Then</p><p>D s = n + μ e μ + ( k T e ) e ( μ + − μ e ) [ μ + − μ e μ + n + − μ e n e ] (10)</p><p>D s = n + μ e μ + ( k T e ) e [ 1 μ + n + − μ e n e ]</p><p>Then</p><p>D s = k T e e [ μ e μ + μ + − μ e n e n + ] (11)</p><p>or</p><p>D s = k T e e [ μ e 1 − μ e μ + n e n + ] (12)</p><p>Substituting μ e = e m e ν e and μ i = e m i ν i from Equations (5) and (7) into Equation (12) results in</p><p>D S = k T e e [ ( e m e ν e ) 1 − ( e m e ν e ) n e ( e m i ν i ) n + ] (13)</p><p>or</p><p>D S = k T e [ 1 m e ν e 1 − m i ν i n e m e ν e n + ] (14)</p><p>or</p><p>D S = k T e [ n + m e ν e n + − m i ν i n e ] (15)</p><p>D S ν i = k T e m e ν e ν i − m i ν i 2 n e n + (16)</p><p>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 [<xref ref-type="bibr" rid="scirp.93013-ref18">18</xref>] [<xref ref-type="bibr" rid="scirp.93013-ref19">19</xref>] stating that:</p><sec id="s3_1"><title>3.1. In the Case of L (Axial Distance) &gt; R (Radial Distance)</title><p>D S ν i = ( R 2.405 ) 2 (17)</p><p>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</p><p>( R 2.4 ) 2 = k T e m e ν e ν i − m i ν i 2 n e n + (18)</p><p>Substituting the following equation into Equation (18)</p><p>ν e = N n 〈 Q e - n υ e - n 〉 (19a)</p><p>ν e = ( 3.55 &#215; 10 16 P ) 〈 Q e - n 〉 ( 2 k T e m e ) 1 / 2 (19b)</p><p>ν i = ( 3.55 &#215; 10 16 P ) 〈 Q i - n 〉 ( 2 k T i m i ) 1 / 2 (19c)</p><p>where P in torr, Q e - n and Q i - n are the cross sections of electron-neutral and ion-neutral collisions, respectively results in</p><p>( R 2.4 ) 2 = k T e [ m e ( 2 k T e m e ) 1 / 2 ( 2 k T i m i ) 1 / 2 ( 3.55 &#215; 10 16 P ) 2 〈 Q i 〉 〈 Q e 〉 ] − [ n e n + m i [ ( 3.55 &#215; 10 16 P ) 2 〈 Q i 〉 ] 2 ( 2 k T i m i ) ] (20)</p><p>then</p><p>( R P 2.4 ) 2 = 1 [ ( 4 m e k T i m i k T e ) 1 / 2 ( 3.55 &#215; 10 16 ) 2 〈 Q i 〉 〈 Q e 〉 ] − [ n e n + ( 3.55 &#215; 10 16 〈 Q i 〉 ) 2 ( 2 k T i m i k T e ) ] (21)</p><p>But from [<xref ref-type="bibr" rid="scirp.93013-ref20">20</xref>] [<xref ref-type="bibr" rid="scirp.93013-ref21">21</xref>] μ e μ i ≅ 7.64 ( m i ) 0.5</p><p>Or</p><p>m i ≅ ( 1 7.64 μ e μ i ) 2 (22a)</p><p>Then substituting (7) into (22a) gives</p><p>1 m i ≅ ( 7.64 k T i k T e ) 2 (22b)</p><p>Moreover, substituting (22b), m i m e ≅ 2000 and ν i ν e ≅ 10 − 2 into (21) gives</p><p>( R P 2.4 ) 2 = 1 [ ( 7.64 k T i k T e ) ( 4 m e k T i k T e ) 1 / 2 ( 3.55 &#215; 10 16 ) 2 〈 Q i 〉 〈 Q e 〉 ] − [ n e n + ( 3.55 &#215; 10 16 〈 Q i 〉 ) 2 ( 7.64 k T i k T e ) 2 ( 2 k T i k T e ) ] (23)</p><p>( R P ) 2 = 2.4 2 [ ( 7.64 k T i k T e ) ( 4 m e k T i k T e ) 1 / 2 ( 3.55 &#215; 10 16 ) 2 〈 Q i 〉 〈 Q e 〉 ] − [ n e n + ( 3.55 &#215; 10 16 〈 Q i 〉 ) 2 ( 7.64 k T i k T e ) 2 ( 2 k T i k T e ) ] (24)</p><p>Substituting the values of the cross sections, masses of ions and electrons into (24) gives</p><p>( R P ) 2 = 2.4 2 ( 3.55 ) 2 ( 7.64 ) ( 2 ) [ [ ( k T i k T e ) ( m e k T i k T e ) 1 / 2 ] − [ n e n + 7.64 ( k T i k T e ) 2 ( k T i k T e ) ] ] (25)</p><p>or</p><p>( R P ) 2 = 2.4 2 ( 3.55 ) 2 ( 7.64 ) ( 2 ) ( k T i k T e ) 3 / 2 [ m e 1 / 2 − 7.64 n e n + ( k T i k T e ) 3 / 2 ] (26)</p><p>Then finally</p><p>R P = 0.1729 1 ( k T i k T e ) 3 / 4 [ m e 1 / 2 − 7.64 n e n + ( k T i k T e ) 3 / 2 ] 1 / 2 (27)</p><p>Using the data of the electron temperatures and densities shown in <xref ref-type="fig" rid="fig4">Figure 4</xref> and <xref ref-type="fig" rid="fig5">Figure 5</xref> respectively, into the Equation (27), taking into account that T e ≫ T i ( T i = 0.1 T e ), m e m i ≅ 2000 [<xref ref-type="bibr" rid="scirp.93013-ref22">22</xref>] and for the spherical probe, the ion density is given by [<xref ref-type="bibr" rid="scirp.93013-ref23">23</xref>] :</p><p>n + = I + 0.6 ( k T e m i ) 1 / 2 A P (28)</p><p>Knowing T<sub>e</sub>, ( A p = 4 π r p 2 ) area of spherical probe) and calculating positive ion current I<sub>+</sub> from the I-V characteristic curve of the spherical single probe [<xref ref-type="bibr" rid="scirp.93013-ref24">24</xref>] , n<sub>+ </sub>can be determined.</p><p><xref ref-type="fig" rid="fig7">Figure 7</xref> 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<sub>e</sub> at the edge of the electrode are more pronounced at the center due to edge effect [<xref ref-type="bibr" rid="scirp.93013-ref25">25</xref>] .</p></sec><sec id="s3_2"><title>3.2. In the Case of L (Axial Distance) &lt; R (Radial Distance)</title><p>D S ν i = ( R 2.405 ) 2 + ( L π ) 2 (29)</p><p>and by substituting Equation (29) into (16), then</p><p>( R 2.4 ) 2 + ( L π ) 2 = k T e m e ν e ν i − m i ν i 2 n e n + (30)</p><p>Substituting by Equations (19a, b and c) into Equation (30), respectively results in</p><p>( R 2.4 ) 2 + ( L π ) 2 = k T e [ m e ( 2 k T e m e ) 1 / 2 ( 2 k T i m i ) 1 / 2 ( 3.55 &#215; 10 16 P ) 2 〈 Q i 〉 〈 Q e 〉 ] − [ n e n + m i [ ( 3.55 &#215; 10 16 P ) 2 〈 Q i 〉 ] 2 ( 2 k T i m i ) ] (31)</p><p>then</p><p>( R 2.4 ) 2 + ( L π ) 2 = 1 [ ( 4 m e k T i m i k T e ) 1 / 2 ( 3.55 &#215; 10 16 ) 2 〈 Q i 〉 〈 Q e 〉 ] − [ n e n + ( 3.55 &#215; 10 16 〈 Q i 〉 ) 2 ( 2 k T i m i k T e ) ] (32)</p><p>Substitution from (22a and b) into (32) gives</p><p>( R P 2.4 ) 2 + ( L π ) 2 = 1 [ ( 7.64 k T i k T e ) ( 4 m e k T i k T e ) 1 / 2 ( 3.55 &#215; 10 16 ) 2 〈 Q i 〉 〈 Q e 〉 ] − [ n e n + ( 3.55 &#215; 10 16 〈 Q i 〉 ) 2 ( 7.64 k T i k T e ) 2 ( 2 k T i k T e ) ] (33)</p><p>( R P ) 2 = 2.4 2 [ ( 7.64 k T i k T e ) ( 4 m e k T i k T e ) 1 / 2 ( 3.55 &#215; 10 16 ) 2 〈 Q i 〉 〈 Q e 〉 ] − [ n e n + ( 3.55 &#215; 10 16 〈 Q i 〉 ) 2 ( 7.64 k T i k T e ) 2 ( 2 k T i k T e ) ] − 2.4 2 ( L π ) 2 (34)</p><p>Substituting the values of the cross sections, masses of ions and electrons into (34) gives finally</p><p>R P = ( 0.0299 ( k T i k T e ) 3 / 2 [ m e 1 / 2 − 7.64 n e n + ( k T i k T e ) 3 / 2 ] − 5.76 ( L π ) 2 ) (35)</p><p>For 1.2 torr of applied pressure, <xref ref-type="fig" rid="fig8">Figure 8</xref> shows comparison between the radial distribution of the electron temperature (RDFT) for L &lt; R using Equation (35) and for L &gt; 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.</p><p>The reasons for the radial distribution of the electron temperatures for both cases.</p><p>L &lt; R and L &gt; R are the following:</p><p>a) The temperature decreases for increasing pressure may be due to the relation between the temperature and the pressure given by [<xref ref-type="bibr" rid="scirp.93013-ref26">26</xref>] :</p><p>| T e | = 6.5 | P − 0.066 | (36)</p><p>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 [<xref ref-type="bibr" rid="scirp.93013-ref27">27</xref>] .</p><p>c) The present theoretical and experimental results are compatible and agree fairly well with the behavior of the data of von Engel [<xref ref-type="bibr" rid="scirp.93013-ref28">28</xref>] [<xref ref-type="bibr" rid="scirp.93013-ref29">29</xref>] [<xref ref-type="bibr" rid="scirp.93013-ref30">30</xref>] for the radial distribution of the electron temperature.</p><p>For L &lt; R, RDFT accurately shows a dramatic decrease by 60% less than RDFT values for L &gt; R. This may be attributed to:</p><p>The cathode fall, negative glow and positive column regions are compressed when L &lt; 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 [<xref ref-type="bibr" rid="scirp.93013-ref23">23</xref>] , therefore the current and current density are increased.</p></sec></sec><sec id="s4"><title>4. Conclusion</title><p>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.</p></sec><sec id="s5"><title>Acknowledgements</title><p>GVO acknowledges for the partial financial support from MEPhI and NRU MPEI in the framework of the Russian Academic Excellence Project.</p></sec><sec id="s6"><title>Conflicts of Interest</title><p>The authors declare no conflicts of interest regarding the publication of this paper.</p></sec><sec id="s7"><title>Cite this paper</title><p>Galaly, A.R. and Van Oost, G. (2019) Comparison between Theoretical and Experimental Radial Electron Temperature Profiles in a Low Density Weakly Ionized Plasma. Journal of Modern Physics, 10, 699-716. https://doi.org/10.4236/jmp.2019.107050</p></sec><sec id="s8"><title>Appendix: Proofing of the Schottky Condition</title><p>There are two classes of Schottky condition models using.</p><p>Appendix A: the recurrence relation for the Bessel condition</p><p>Appendix B: the boundary conditions</p><p>D S ν i = K T e m i ν i 2 n e n + + m e ν e ν i (36)</p><p>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</p><p>D S ν i = ( R 2.4 ) 2 can be demonstrated as follows:</p>Appendix A: From the Recurrence Relation for Bessel Condition [<xref ref-type="bibr" rid="scirp.93013-ref31">31</xref>] [<xref ref-type="bibr" rid="scirp.93013-ref32">32</xref>] [<xref ref-type="bibr" rid="scirp.93013-ref33">33</xref>]<p>2 J ′ n ( x ) = J n − 1 ( x ) − J n + 1 ( x ) (37)</p><p>and</p><p>2 n x J n ( x ) = J n − 1 ( x ) + J n + 1 ( x ) (38)</p><p>Subtracting equation (37) from (38), gives</p><p>2 n x J n ( x ) − 2 J ′ n ( x ) = J n − 1 ( x ) + J n + 1 ( x ) − J n − 1 ( x ) + J n + 1 ( x ) (39)</p><p>Or</p><p>n x J n ( x ) − J ′ n ( x ) = J n + 1 ( x )</p><p>Or</p><p>n J n ( x ) − x J ′ n ( x ) = x J n + 1 ( x )</p><p>This results in</p><p>x J ′ n ( x ) = n J n ( x ) − x J n + 1 ( x ) (40)</p><p>Differentiating (40) with respect to x, gives:</p><p>x J ″ n ( x ) + J ′ n ( x ) = n J ′ n ( x ) − x J ′ n + 1 ( x ) − J n + 1 ( x ) (41)</p><p>Multiplying Equation (40) in n/x gives:</p><p>n J ′ n ( x ) = n 2 x J n ( x ) − n J n + 1 ( x ) (42)</p><p>Adding Equations (37) and (38) gives</p><p>x J ′ n ( x ) = x J n − 1 ( x ) − n J n ( x ) (43)</p><p>Putting n = n + 1 in Equation (43) gives</p><p>x J ′ n + 1 ( x ) = x J n ( x ) − ( n + 1 ) J n + 1 ( x ) (44)</p><p>Substituting (39) and (40) into (38), results in</p><p>x J ″ n ( x ) + J ′ n ( x ) = n 2 x J n ( x ) − x J n ( x )</p><p>Then</p><p>J ″ n ( x ) + 1 x J ′ n ( x ) + ( 1 − n 2 x 2 ) J n ( x ) = 0 (45)</p><p>Finally</p><p>J ″ + 1 x J ′ + ( λ − n 2 x 2 ) J = 0 (46)</p><p>For n = 0</p><p>J ″ + 1 x J ′ + ( λ ) J = 0 (47)</p><p>For the Bessel condition in cylindrical coordinate</p><p>x = K r = ( 1 D a τ ) 1 / 2 R = ( λ D a ) 1 / 2 R</p><p>with the solution J n ( x ) = 2 π cos ( x − n π 2 − π 4 ) , and J 0 = 0 for n = 0, gives</p><p>0 = cos ( x − π 4 ) with the solution x = 3 π 4 ≈ 2.4</p><p>x = K r = R ( D τ ) 1 / 2 ≈ 2.4 or R 2.4 = ( D τ ) 1 / 2</p><p>Then</p><p>R 2.4 = ( D S ν i ) 1 / 2</p><p>or</p><p>D S ν i = ( R 2.4 ) 2</p>Appendix B: From the Boundary Conditions [<xref ref-type="bibr" rid="scirp.93013-ref34">34</xref>] [<xref ref-type="bibr" rid="scirp.93013-ref35">35</xref>]<p>n e υ e = − D S ∇ n e</p><p>Then</p><p>D S = − n e υ e ∇ n e (48)</p><p>But</p><p>ν i = ∇ ⋅ n e υ e n e</p><p>Then</p><p>D S ν i = − ( n e ∇ n e ) 2 (49)</p><p>For boundary conditions r = R , n e = n e 0 , and J 0 = 2.4 , then</p><p>∇ n e = − 2.4 n e R (50)</p><p>Substituting from (50) into (49) then we get</p><p>D S ν i = ( R 2.4 ) 2</p></sec></body><back><ref-list><title>References</title><ref id="scirp.93013-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Chen, F.F. (1965) Plasma Diagnostic Techniques. Academic Press, New York.</mixed-citation></ref><ref id="scirp.93013-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">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</mixed-citation></ref><ref id="scirp.93013-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Kagan, Y. and Paskalev, K.K. (1975) Soviet Physics. Technical Physics, 19, 1604.</mixed-citation></ref><ref id="scirp.93013-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Kagan, Y. and Perel, V.I. (1964) Soviet Physics Uspekhi, 6, 767. https://doi.org/10.1070/PU1964v006n06ABEH003611</mixed-citation></ref><ref id="scirp.93013-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Mott-Smith and Langmuir (1926) Physical Review, 28, 727. https://doi.org/10.1103/PhysRev.28.727</mixed-citation></ref><ref id="scirp.93013-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">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</mixed-citation></ref><ref id="scirp.93013-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Holt, E.H. and Haskell, R.E. (1965) Foundations of Plasma Dynamics. 2nd Edition, MacMillan Company, New-York.</mixed-citation></ref><ref id="scirp.93013-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Huang, P.G., Shang, J.S. and Stanfield, S.A. (2011) AIAA Journal, 49, 119. https://doi.org/10.2514/1.J050463</mixed-citation></ref><ref id="scirp.93013-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Bisek, N.J., Boyd, I.D. and Poggie (2009) Journal of Spacecraft and Rockets, 46, 568. https://doi.org/10.2514/1.39032</mixed-citation></ref><ref id="scirp.93013-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Galaly, A.R. and Khedr, M.A. (2015) British Journal of Applied Science &amp; Technology, 11, 1-14. https://doi.org/10.9734/BJAST/2015/17663</mixed-citation></ref><ref id="scirp.93013-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Galaly, A.R. (2016) British Journal of Applied Science &amp; Technology, 13, 1.</mixed-citation></ref><ref id="scirp.93013-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">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</mixed-citation></ref><ref id="scirp.93013-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Yasuda, H.K., Tao, W.H. and Prelas, M.A. (1996) Journal of Vacuum Science &amp; Technology A: Vacuum, Surfaces, and Films, 14, 2113. https://doi.org/10.1116/1.580089</mixed-citation></ref><ref id="scirp.93013-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">Nasser, E. (1971) Fundamentals of Gaseous Ionization and Plasma Electronics. Wiley-Interscience, New York.</mixed-citation></ref><ref id="scirp.93013-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Raizer, Y.P. (1999) Gas Discharge Physics. Springer, New York.</mixed-citation></ref><ref id="scirp.93013-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">Von Engel, A.H. (1964) Ionized Gases. Oxford University Press, Oxford.</mixed-citation></ref><ref id="scirp.93013-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">Franklin, R.N. (1977) Plasma Phenomena in Gas Discharges. Oxford University Press, Oxford.</mixed-citation></ref><ref id="scirp.93013-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">Lieberman, M.A. and Lichtenberg, A.J. (1994) Principles of Plasma Discharges for Materials Processing. Wiley Interscience, New York.</mixed-citation></ref><ref id="scirp.93013-ref19"><label>19</label><mixed-citation publication-type="other" xlink:type="simple">Braithwaite, N.St.J. (2000) Plasma Sources Science and Technology, 9, 517. https://doi.org/10.1088/0963-0252/9/4/307</mixed-citation></ref><ref id="scirp.93013-ref20"><label>20</label><mixed-citation publication-type="other" xlink:type="simple">Piel, A. (2011) Plasma Physics. Springer, Berlin.</mixed-citation></ref><ref id="scirp.93013-ref21"><label>21</label><mixed-citation publication-type="other" xlink:type="simple">Kolobov, V.I. (2013) Physics of Plasmas, 20, Article ID: 101610. https://doi.org/10.1063/1.4823472</mixed-citation></ref><ref id="scirp.93013-ref22"><label>22</label><mixed-citation publication-type="other" xlink:type="simple">Fruchtman, A. (2009) Plasma Sources Science and Technology, 18, Article ID: 025033. https://doi.org/10.1088/0963-0252/18/2/025033</mixed-citation></ref><ref id="scirp.93013-ref23"><label>23</label><mixed-citation publication-type="other" xlink:type="simple">Simon, A. (1955) Physical Review, 98, 317. https://doi.org/10.1103/PhysRev.98.317</mixed-citation></ref><ref id="scirp.93013-ref24"><label>24</label><mixed-citation publication-type="other" xlink:type="simple">Galaly, A.R. (2014) Physical Science International Journal, 4, 930-939. https://doi.org/10.9734/PSIJ/2014/9604</mixed-citation></ref><ref id="scirp.93013-ref25"><label>25</label><mixed-citation publication-type="other" xlink:type="simple">Hutchinson (1990) Principles of Plasma Diagnostics. Cambridge University Press, Cambridge.</mixed-citation></ref><ref id="scirp.93013-ref26"><label>26</label><mixed-citation publication-type="other" xlink:type="simple">Huba, J.D. (2011) Plasma Formulary. Naval Research Laboratory Report NRL/PU/ 6790-94-265.</mixed-citation></ref><ref id="scirp.93013-ref27"><label>27</label><mixed-citation publication-type="other" xlink:type="simple">Druyvesteyn, M.J. and Penning, F.M. (1940) Reviews of Modern Physics, 12, 87. https://doi.org/10.1103/RevModPhys.12.87</mixed-citation></ref><ref id="scirp.93013-ref28"><label>28</label><mixed-citation publication-type="other" xlink:type="simple">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</mixed-citation></ref><ref id="scirp.93013-ref29"><label>29</label><mixed-citation publication-type="other" xlink:type="simple">Von Engel, A. and Cozens, J.R. (1965) International Journal of Electronics, 19, 61-68. https://doi.org/10.1080/00207216508937799</mixed-citation></ref><ref id="scirp.93013-ref30"><label>30</label><mixed-citation publication-type="other" xlink:type="simple">Von Engel, A. (1956) Proceedings of the Physical Society, 59, 468.</mixed-citation></ref><ref id="scirp.93013-ref31"><label>31</label><mixed-citation publication-type="other" xlink:type="simple">Bowman, F. (1958) Introduction to Bessel Conditions. Dover, New York.</mixed-citation></ref><ref id="scirp.93013-ref32"><label>32</label><mixed-citation publication-type="other" xlink:type="simple">Korenev, B.G. (2002) Bessel Condition and Their Applications. CRC Press Company, Boca Raton, London, New York, Washington DC.</mixed-citation></ref><ref id="scirp.93013-ref33"><label>33</label><mixed-citation publication-type="other" xlink:type="simple">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</mixed-citation></ref><ref id="scirp.93013-ref34"><label>34</label><mixed-citation publication-type="other" xlink:type="simple">Gray, A. and Mathews, G.B. (1966) A Treatise on Bessel Conditions and Their Applications to Physics. 2nd Edition, Dover, New York.</mixed-citation></ref><ref id="scirp.93013-ref35"><label>35</label><mixed-citation publication-type="other" xlink:type="simple">Watson, G.N. (1966) A Treatise on the Theory of Bessel conditions. 2nd Edition, Cambridge University Press, Cambridge.</mixed-citation></ref></ref-list></back></article>