Session 14: Plasmas in Equilibrium

Session : Plasmas in Equilibrium Ionization and Conduction in a High-pressure Plasma A normal gas at T < 3000◦ K is a good electrical insulator, because there are almost no free electrons in it. For pressures > 0.1 atm, collision among molecule and other particles are frequent enough that we can assume Local Thermodynamic Equilibrium (LTE), and in particular, ionization-recombination reactions are governed by the law of mass action. Consider neutral atoms (n) which ionize singly to ions (i) and electrons (e): n ↔ e+i (1) One form of the law of mass action (in terms of number densities nj = pj /kT , where T is the same for all species) is, n e ni = S(T ) nn (2) Where the “Saha function” S is given (according to Statistical Mechanics) as, S(T ) = 2 qi qn 2πme kT h2 3/2 exp − eVi kT (3) − Ground state degeneracy of ion (qi = 1 for H+) − Ground state degeneracy of neutral (qn = 2 for H) − Mass of electron (me = 9.11 × 10−31 kg) − Boltzmann constant (k = 1.38 × 10−23 J/K) − Plank’s constant (h = 6.62 × 10−34 Js) − Ionization potential of the atom (Vi = 13.6 V for H) Except for very narrow “sheaths” near walls, plasmas are quasi-neutral: n e = ni (4) n2e = S(T ) nn (5) So that, can be used. The ionization fraction is defined as α = ne . ne + nn Given T , this relates ne to nn . A second relation is needed and very often it is a specification of the overall pressure, 1 p = (ne + ni + nn )kT = (2ne + nn )kT (6) Combining (5) and (6), n2e � p � = S(T ) − 2ne = S(T ) (n − 2ne ) kT Where n = p/kT is the total number density of all particles. We then have, n2e + 2Sne − Sn = 0 and, ne = −S + √ S 2 + Sn = n _ n 1+ 1+ S (7) Since S increases very rapidly with T , the limits of (7) are, √ Sn n ne (T → ∞) → 2 ne (T ≈ 0) → Weak ionization Full ionization E will drive a current density Ej through Once an electron population exists, an electric field E the plasma. To understand this quantitatively, consider the momentum balance of a “typical” electron. It sees an electrostatic force, E FEe = −eE (8) It also sees a “frictional” force due to transfer of momentum each time it collides with some other particle (neutral or ion). Collisions with other electrons are not counted, because the momentum transfer is in that case internal to the electron species. The ions and neutrals are almost at rest compared to the fast-moving electrons, and we define an effective collision as one in which the electron’s directed momentum is fully given up. Suppose there are νe of these collisions per second (νe is the collision frequency per electron). The electron loses momentum at a rate −meEuνe , where Eu is the mean directed velocity of electrons, and so, FEf riction = −meEu νe On average, FEe + FEf riction = 0 E so that, or meEuνe = −eE, 2 (9)  ~u = − e m e νe  ~ E (10) e is the electron “mobility” [(m/s)/(volt/m)]. The current density is me νe the flux of charge due to motion of all charges. If only the electron motion is counted (it dominates in this case), The group µ = ~j = −ene~u =  ne e2 me νe  ~ E (11) The group, σ = ne e 2 me νe (12) is the electrical conductivity of the plasma (Si/m). Let us consider the collision frequency. Suppose a neutral is regarded as a sphere with a cross-section area Qen . Electrons moving at random with (thermal) velocity ce intercept the area Qen at a rate equal to their flux ne ce Qen . Since a whole range of speeds ce exists, we use the average value c̄e for all electrons. But this is for all electrons colliding with one neutral. We are interested in the reverse (all neutrals, one electron), so the part of νe due to neutrals should be nn c̄e Qen . Adding the part due to ions, νe = nn c̄e Qen + ni c̄e Qei (13) c̄e is very different (usually much larger) than ue . Most of the thermal motion is fast, but in random directions, so that on average it nearly cancels out. The non-canceling remainder is ue . Think of a swarm of bees moving furiously back and forth, but moving (as a swarm) slowly. The number of electrons per unit volume that have a velocity vector ~ce ending in a “box”, d3 ce = dcex dcey dcez in velocity space is defined as fe (~ce , ~x)d3 ce , where fe (~ce , ~x) is the distribution function of the electrons, which depends (for a given location ~x and time t) on the three components of ~ce . In an equilibrium situation all directions are equally likely (isotropic), so fe = fe (ce ) only, and one can show that the form is Maxwellian,  f e = ne me 2πkTe 3/2   me c2e exp − 2kTe where With the normalization, ZZZ f e d 3 ce = n e The mean velocity is then, 3 c2e = c2ex + ce2y + c2ez (14) 1 c̄e = ne ZZZ r ce fe d3 ce = 8k Te πme (15) For electrons, c̄e = 6210 p Te (m/s with Te in ◦ K) (16) Note that if there is current, the distribution cannot be strictly Maxwellian (or even isotropic). But since ue  ce , the mean thermal velocity is very close to Eq. (15) anyway. Regarding the cross sections Qen and Qei , they depend on the collision velocity c̄e , especially Qei . This is because the e − i coulombic interaction is “soft”, so a very fast electron can pass nearly undeflected near an ion, whereas a slow one will be strongly deflected. The complete theory yields an expression, Qei ≈ 2.95 × 10−10 ln Λ Te2 (m2 with Te in ◦ K) (17) where 1 ln Λ ≈ −11.35 + 2 ln Te − 2 ln p 2 (Te in ◦ K, p in atm) (18) so that ln Λ is usually around 6 − 12, and can even be taken as a constant (∼ 8) in rough calculations. For the neutral hydrogen atoms, the collisions are fairly “hard”, and one can use the approximation, QeH ≈ 2 × 10−19 m2 4 (19)