Problems For The Course F5170 - Introduction To Plasma Physics
Problems For The Course F5170 - Introduction To Plasma Physics
Problems For The Course F5170 - Introduction To Plasma Physics
2014
Contents
1 Introduction 5
1.1 Theory . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.2 Problems . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.2.1 Derivation of the plasma frequency . . . . . . . . . . . 6
1.2.2 Plasma frequency and Debye length . . . . . . . . . . 7
1.2.3 Debye-Hückel potential . . . . . . . . . . . . . . . . . 8
1
CONTENTS 2
4
Preface
This document contains exercise problems for the course F5170 – Intro-
duction to plasma physics. This work was supported by the Ministry of Re-
search and Education of the Czech Republic, project no. FRVŠ 12/2013/G6.
The most valuable source of information for this document was the book
Fundamentals of Plasma Physics by J. A. Bittencourt [4]. The authors
would be grateful for any notification about eventual errors.
The complete and up-to-date version of this document can be found at
http://physics.muni.cz/~sperka/exercises.html.
Contacts
Jiřı́ Šperka [email protected]
Jan Voráč [email protected]
Lenka Zajı́čková [email protected]
Physical constants
Proton rest mass mp 1, 67 · 10−27 kg
Electron rest mass me 9.109 · 10−31 kg
Elementary charge e 1.602 · 10−19 C
Boltzmann’s constant k 1.38 · 10−23 J K−1
Vacuum permittivity ε0 8.854 · 10−12 A2 s4 kg−1 m−3
Used symbols
Vector quantities are typed in bold face (v), scalar quantities, including
magnitudes of vectors are in italic (v). Tensors are usually in upper-case
calligraphic typeface (P).
5
Operators
scalar product a·b
vector product a×b
di
ith derivative with respect to x dxi
∂
partial derivative ∂x
∂ ∂ ∂
nabla operator ∇ = ( ∂x , ∂y , ∂z )
Laplace operator ∆ = ∇2
D ∂
total time derivative Dt = ∂t +u·∇
Physical quantities
electron concentration ne
electron temperature Te
electron plasma frequency ωpe
Debye length λD
Larmor radius rc
Larmor frequency Ωc
magnetic moment m
force F
electric field intensity E
magnetic field induction B
arb. quantity for one type of particles χα
distribution function f (χα )
mean velocity u
charge denisty ρ
mass density ρm
collision frequency ν
source term due to collisions Sα
scalar pressure p
tensor of kinetic pressure P
mobility of particles Mα
6
Chapter 1
Introduction
1.1 Theory
Electron plasma frequency
s
ne e2 √
ωp = = const ne (1.1)
ε0 m e
7
CHAPTER 1. INTRODUCTION 8
1.2 Problems
1.2.1 Derivation of the plasma frequency
Consider a steady initial state with a uniform number density of electrons
and an equal number of ions such that the total electrical charge is neutral.
Neglect the thermal motion of the particles and assume that the ions are
stationary. Show that a small displacement of a group of electrons leads to
oscillations with the plasma frequency according to the equation (1.1).
Solution The situation is sketched in the figure 1.1. Assume that the
electric field in the plane perpendicular to the x-axis is zero (just like in
the case of an infinitely large charged plane or capacitor). Let us apply the
Gauss’s law to a closed cylindrical surface (only contour of which is sketched
in the figure): I
Q Sne e
E · dS = = x, (1.3)
S 0 0
CHAPTER 1. INTRODUCTION 9
where S is the area of the base of the cylinder. The resulting electric field is
n0 e
Ex = x. (1.4)
0
Inserting this electric field into the equation of motion of a single electron
yields
d2 x n0 e 2
+ x = 0, (1.5)
dt2 me 0
Which is an equation of a harmonic oscillator with the frequency
1/2
n0 e 2
ωpe = . (1.6)
me 0
(a) Earth’s ionosphere with electron concentration ne = 106 cm−3 and elec-
tron temperature k Te = 0.2 eV.
[ωp = 5, 6 × 107 rad · s−1 = 3, 5 × 108 Hz, λD = 3, 3 mm]
(b) A cell of a typical plasma display with electron concentration of 1013 cm−3
and electron temperature of 1 eV. The cell dimension is about 100 µm.
Is the condition that the system dimension should be much greater than
the Debye length fulfilled?
[ωpe = 2, 3 × 1013 rad · s−1 = 3, 6 THz, λD = 21 nm]
(c) A welding arc with electron concentration of 1, 6×1017 cm−3 and electron
temperature of 1, 3 eV
[ωpe = 2, 3 × 1013 rad · s−1 = 3, 6 THz, λD = 21 nm]
(d) A fluorescent lamp with electron concentration of 1010 cm−3 and electron
temperature of 1 eV
[ωpe = 5, 6 × 109 rad · s−1 = 0, 90 GHz, λD = 74 µm]
CHAPTER 1. INTRODUCTION 10
ϕ(r) ne e 2
∇2 ϕ(r) = 2 = ϕ(r) (1.8)
rD ε0 kTe
Solution Put simply the Debye-Hückel potential into the equation (1.8)
and calculate Laplace operator in spherical coordinates
∂2f
1 ∂ 2 ∂f 1 ∂ ∂f 1
∆f = 2 r + 2 sin θ + 2 2 (1.9)
r ∂r ∂r r sin θ ∂θ ∂θ r sin θ ∂ϕ2
Chapter 2
Motion of particles in
electromagnetic fields
2.1 Theory
The Lorentz force is the combination of electric and magnetic force on
a point charge due to electromagnetic fields. If a particle of charge q moves
with velocity v in the presence of an electric field E and a magnetic field B,
then it will experience the Lorentz force
F = q (E + v × B) . (2.1)
The gyroradius (also known as Larmor radius or cyclotron radius) is
the radius of the circular motion of a charged particle in the presence of a
uniform magnetic field:
mv⊥
rg = , (2.2)
|q|B
where rg is the gyroradius, m is the mass of the charged particle, v⊥ is the
velocity component perpendicular to the direction of the magnetic field, q
is the charge of the particle, and B is magnitude of the constant magnetic
field.
Similarly, the frequency of this circular motion is known as the gyrofre-
quency or cyclotron frequency, and is given by:
|q|B
ωg = . (2.3)
m
Note: A cyclotron is a type of particle accelerator in which charged particles
accelerate due to high-frequency electric field. The cyclotron was invented
and patented by Ernest Lawrence of the University of California, Berkeley,
where it was first operated in 1932.
11
CHAPTER 2. MOTION OF PARTICLES IN ELECTROMAGNETIC FIELDS12
2.2 Problems
2.2.1 Magnetic mirror
Magnetic mirrors are used to confine charged particles in a limited vol-
ume. The gradient of magnetic field induction can result in reversing the
direction of drift of a charged particle.
Suppose we have an electron located at z = 0 with initial velocity v0 and
an initial pitch angle ϑ. The magnetic field induction is given by
B(z) = B0 1 + (γ z)2 .
(2.4)
Solution We start with the conservation of kinetic energy and the mag-
netic moment. The kinetic energy conservation condition yields
The z-component of the velocity at the turning point must be zero, which we
immediately use in the equation describing the conservation of the magnetic
moment
me v02 sin2 ϑ me vt2
=
2 B0 2 B0 (1 + (γ zt )2 )
v0 sin ϑ 1 + (γ zt ) ) = vt2
2 2 2
= v02
1 − sin2 ϑ
γ 2 zt2 =
sin2 ϑ
1
zt = . (2.6)
γ tan ϑ
We see that the position of the point of reflection depends only on the
gradient of the magnetic field and on the initial pitch angle.
B(z) = B0 1 + (γ z)4 .
(2.7)
The
initial pitch angle is ϑ.
1/2
1
zt = γ tan ϑ
CHAPTER 2. MOTION OF PARTICLES IN ELECTROMAGNETIC FIELDS13
(b) After that, the electron, having the velocity vx , enters the second part,
→
−
where a transverse electric field E of the magnitude 10−10 V m−1 is
applied. This field is generated by the plates of a capacitor with the
length d = 1 m. What is the vertical displacement of the electron with
respect to the starting position at the end of the second part? First
derive the general solution.
[35.2 m]
→
−
(c) Finally, the electron enters a homogeneous magnetic field B of the mag-
nitude 20.6 µT (this is the magnitude of the horizontal component of the
geomagnetic field induction in Brno). Calculate the Larmor radius, cy-
clotron frequency and the magnitude of the magnetic moment of the
rotating electron.
[rc = 9.72 × 10−6 m, Ωc = 3.6 × 106 rad · s−1 , |m| = 2.7 × 10−23 A · m2 ]
(d) What would be the result for proton, neutron and positron? For illus-
tration see 2.1.
2.2.4 E × B drift
Suppose we have a vacuum chamber with electric field E = 1 kV m−1
perpendicular to magnetic field B = 1 mT. Calculate E × B drift speed for
an
electron inside the chamber.
E
B
CHAPTER 2. MOTION OF PARTICLES IN ELECTROMAGNETIC FIELDS14
Solution Gyroradius:
γβm0 c
r= (2.8)
qB
Angular gyrofrequency:
r
|q|B Ωc p v 2
Ωrel
c = = = Ωc 1 − β 2 = Ωc 1− (2.9)
γm0 γ c
Energy:
mc2
Ek = mγc2 − mc2 = p − mc2 (2.10)
1 − v 2 /c2
h Derive continuity
i equation from Maxwell’s Equations.
∂ρ
∂t +∇·J=0
3.1 Theory
• Phase space is defined by six coordinations (x, y, z, vx , vy , vz ).
16
CHAPTER 3. ELEMENTS OF PLASMA KINETIC THEORY 17
3.2 Problems
3.2.1 Uniform distribution function
Suppose we have system of particles uniformly distributed in space with con-
stant particle number density n, which is characterised by one dimensional
distribution function of speeds F (v):
F (v) = C for v ≤ v0
F (v) = 0 otherwise,
4.1 Theory
• The macroscopic variables, such as number density, flow velocity, ki-
netic pressure or thermal energy flux can be considered as average
values of physical quantities, involving the collective behaviour of a
large number of particles. These macroscopic variables are related to
the various moments of the distribution function.
• The average value of the property χα (~r, ~v , t) for the particles of type
α is defined by
Z
1
hχα (~r, ~v , t)i = χα (~r, ~v , t)fα (~r, ~v , t) d3 v. (4.1)
nα (~r, t) ~v
• For example, the average velocity (or flow velocity) ~uα (~r, t) for the
particles of type α is defined by
Z
1
~uα (~r, t) = hvα (~r, t)i = ~v fα (~r, ~v , t) d3 v. (4.2)
nα (~r, t) ~v
18
CHAPTER 4. AVERAGE VALUES AND MACROSCOPIC VARIABLES19
4.2 Problems
4.2.1 RMS speed
What is the rms speed of the following three electrons (|v1 | = 1, |v2 | = 2
√ |v3 | = 5)?
and
[ 10]
0.002
a)
b)
c)
0.0015
0.001
n
0.0005
0
0 1e+06 2e+06 3e+06 4e+06 5e+06
v [m/s]
0.002
a)
b)
c)
0.0015
0.001
n
0.0005
0
0 1e+06 2e+06 3e+06 4e+06 5e+06
v [m/s]
Figure 4.2: Diagram to the problem of the highest particle density 4.2.5.
CHAPTER 4. AVERAGE VALUES AND MACROSCOPIC VARIABLES21
5.1 Theory
• The equilibrium distribution function fαEq (~r, ~v , t) is the time-independent
solution of the Boltzmann equation in the absence of external forces.
The gamma function Γ(n) is an extension of the factorial function, with its
argument shifted down by 1, to real and complex numbers. That is, if n is
a positive integer:
Γ(n) = (n − 1)! (5.2)
Other important formulas:
∞ Γ( (n+1) √
2 )
Z
n −a x2 1
x e dx = (n+1)
; Γ = π. (5.3)
0 2a 2 2
22
CHAPTER 5. THE EQUILIBRIUM STATE 23
5.2 Problems
5.2.1 Gamma function
Starting from the definition of a Gamma function show that, if n is a positive
integer, then
Γ(n + 1) = n!
R∞
Recipe: First, using integration by parts of Γ(n+1) = 0 xn e−x dx demon-
strate that
Γ(a + 1) = aΓ(a).
Next, it remains to show, that Γ(1) = 1.
(e) Derive the relation for the number of particles passing through a unit
of length in a unit of time from one side (the flux of particles from one
side).
Solution
(a) Integrate the distribution function over the whole velocity space. The
condition that the integral equals the concentration of particles n yields
Z∞ r
m vx2
2kT π
n=C exp − dvx = C . (5.5)
2kT m
−∞
r
m
C=n (5.6)
2kT π
CHAPTER 5. THE EQUILIBRIUM STATE 24
(b) Distribution of particle speeds F (v) from the summation over the both
possible directions is
m v2
r
m
F (v) = 2 n exp − (5.7)
2kT π 2kT
(c) From the condition that the derivation of the distribution F (v) must
equal zero
m v2
0 = v exp − (5.8)
2kT
we will get that the most probable speed is zero.
(d)
Z∞ r
2kT
hvi = v F (v) dv = (5.9)
πm
0
(e)
Z∞ r
kT
Γ= vx f (vx )dvx = n (5.10)
2πm
0
Results:
mn
(a) C = 2πkT
h i
nm v2
(b) F (v) = 2 π v f (v) = kT v exp − m
2kT
q
kT
(c) Most probable speed v = m .
q
kT π
(d) Mean speed hvi = 2m .
q
kT
(e) Γ = n 2mπ
CHAPTER 5. THE EQUILIBRIUM STATE 25
Results:
m
3/2
(a) C = n 2πkT
3/2 2 h i
m v2
(b) F (v) = 4π n 2 πkT v exp − m
2kT
q
2kT
(c) Most probable speed v = m .
q
8kT
(d) Mean speed hvi = πm .
q
kT
(e) Γ = n 2mπ
Particle interactions in
plasmas
6.1 Theory
Collisional phenomena can be divided into two categories:
• inelastic - the internal states of some or all of the particles involved are
changed and particles may be created as well as destroyed. A charged
particle may recombine with another to form a neutral particle or it
can attach itself to a neutral particle to form a heavier charged particle.
The energy state of an electron in an atom may be raised and electrons
can be removed from their atoms resulting in ionization.
The total scattering cross section can be obtained by integrating σ(χ, ε)dΩ
over the entire solid angle:
Z
σt = σ(χ, ε)dΩ. (6.1)
Ω
In the special case, when the interaction potential is isotropic (e.g. Coulomb
potential), we can get the total scattering cross section using the formula
Zπ
σt = 2π σ(χ) sin χdχ. (6.2)
0
For the same case, when the interaction potential is isotropic, we can get
the momentum transfer cross section using the formula:
26
CHAPTER 6. PARTICLE INTERACTIONS IN PLASMAS 27
Zπ
σm = 2π (1 − cos χ)σ(χ) sin χdχ. (6.3)
0
6.2 Problems
6.2.1 Mean free path of Xe ions
Scattering cross section σ for elastic collisions of Xe+ ions with Xe atoms
is approximately independent on their energy with cross section value of
σ = 10−14 cm2 .
A) Calculate mean free path l of Xe+ ions for elastic collisions in a weakly
ionized plasma in xenon atmosphere at room temperature (20 ◦ C) at the
pressure:
a) 1000 Pa
b) 10 Pa
c) 0.1 Pa
B) How long is the time period between two subsequent collisions, if the
mean temperature of Xe ions is T = 1000 K?
Solution:
kT
λ= .
pσ
So the final results for given pressures are:
velocity: r
m
τ =λ .
3kT
So the results are: a) 17 · 10−9 s b) 17 · 10−7 s c) 17 ·10−5 s.
Macroscopic transport
equations
7.1 Theory
From different moments of Boltzmann equation, the following macro-
scopic transport equations can be derived:
• From the condition of conservation of mass the continuity equation
∂ρmα
+ ∇ · (ρmα uα ) = Sα , (7.1)
∂t
where ρmα is the mass density of type-α particles and Sα describes
the creation or destruction of particles due to collisions (ionization,
recombination, etc.).
is the collision term with ναβ being the collision frequency for the
momentum transfer between the particles of type α and particles of
type β. From conservation of momentum during a collision follows
29
CHAPTER 7. MACROSCOPIC TRANSPORT EQUATIONS 30
7.2 Problems
7.2.1 Afterglow
Consider a homogeneous plasma afterglow consisting of electrons and one
type of singly charged positive ions. In this case, the continuity equation is
∂ne
= −kr ne ni , (7.6)
∂t
where kr is the rate coefficient for recombination. The spatial derivatives
vanish because of the spatial uniformity. The concentration of electrons at
ht = 0 is n0 . Calculate
i ne (t > 0). Remember the quasineutrality condition.
n0
ne (t) = n0 kr t+1
This is obviously
uαx,p (t) = C e−ναβ t (7.10)
We now take the parameter C to be time-dependent C = C(t) and calculate
the deriative
duαx (t) dC(t) −ναβ t
= e − C(t) ναβ e−ναβ t (7.11)
dt dt
inserting this into the original equation (7.8) yields
dC(t) −ναβ t
e = ναβ uβx
dt
from which we obtain by integrating
And similarly for all three spatial components. The velocity uα will expo-
nentially approach to the velocity uβ with the rate given by the collision
frequency for momentum transfer ναβ .
Results:
h i
ρmα
(a) u(t) = uα (0) · exp 1+ ρmβ t
h i
uα (0) ρmα
(b) uα (t) = ρmα +ρmβ ρmβ · exp −ναβ 1 + ρmβ t + ρmα
uβ (t) = u(t) + uα (t)
CHAPTER 7. MACROSCOPIC TRANSPORT EQUATIONS 32
and the ideal gas law pe = ne k Te . Show that the heat flow equation can be
written as
Ωce
(qe × B) = −K0 ∇Te + (fe − fe0 ), (7.17)
ν
where
5 k pe
K0 = (7.18)
2 me ν
is the thermal conductivity.
Chapter 8
8.1 Theory
The equations governing the important physical properties of the plasma
as a whole can be obtained by summing the terms for the particular species.
If also several simplifying assumptions are made, the following set of so
called magnetohydrodynamic equations can be derived:
• The continuity equation
∂ρm
+ ∇ · (ρm u) = 0 (8.1)
∂t
• The momentum equation
Du
ρm = J × B − ∇p (8.2)
Dt
• Generalised Ohm’s law
σ0
J = σ0 (E + u × B) − J × B. (8.3)
ne
The electric and magnetic fields are further bound by the Maxwell equations.
In these equations, viscosity and thermal conductivity are neglected.
8.2 Problems
8.2.1 Electric current density
The mean velocity of plasma u is defined as a weighted average of the
mean velocities of the particular species
X ρmα
u= uα (8.4)
α
ρm
33
CHAPTER 8. MACROSCOPIC EQUATIONS FOR A CONDUCTING FLUID34
where ρm is the total mass density of the plasma. Each species has con-
centration nα , charge qα and the so called diffusion velocity wα = uα − u.
Calculate the total electric current density J in terms of the total electric
charge density ρ and the particular densities, charges and diffusion veloci-
ties. Note, that due to the definition of the mean velocity of plasma, the
result is not simply J =
ρ u.
P
J = ρu + nα qα wα
α
and form the equation for the mean velocity of the plasma as a whole
1
u= (ρme ue + ρmi ui ) (8.6)
ρm
derive
h the
drift velocities
ui and ue. i
µ ρm u J µ ρm u J me mi
ui = ρmi me + e , u e = ρme mi − e , µ= me +mi
Dux −σ0 B 2
= ux (8.12)
Dt ρm
This has a form of a simple decay problem. The solution is
σ0 B 2
ux (t) = ux (0) exp − t . (8.13)
ρm
The same holds for uy , so the time dependence ofq the component of the
velocity perpendicular to the magnetic field u⊥ = u2x + u2y is
where
ρm
τ= (8.15)
σ0 B 2
is the characteristic time for diffusion across the magnetic field lines.
Chapter 9
9.1 Theory
In weakly ionised cold plasma, the equation of motion for electrons takes
the simple form of the so called Langevin equation
Due
me = −e (E + ue × B) − νc me ue , (9.1)
Dt
where νc is the collision frequency for momentum transfer between the elec-
trons and the heavy particles.
In the absence of a magnetic field, the current produced by moving
electrons is
J = −e ne ue (9.2)
and the DC conductivity is
ne e2
σ0 = (9.3)
me νc
and the electron mobility
e σ0
Me = − =− . (9.4)
me νc ne e
When an external magnetic field is present, the plasma becomes anisotropic,
and the DC conductivity and electron mobility are described by tensors (see
problem 9.2.2).
In weakly ionised plasma with relatively high density of neutrals, the
diffusion equation for charged species α is
∂nα
= D ∇2 nα . (9.5)
∂t
36
CHAPTER 9. PLASMA CONDUCTIVITY AND DIFFUSION 37
9.2 Problems
9.2.1 DC plasma conductivity
From the Langevin equation for electrons in the absence of magnetic
field and in the steady state
− e E − me νc ue = 0 (9.9)
J = −e ne ue (9.10)
Inserting this into the Langevin equation (9.9), we obtain the expression for
the current density J
ne e2
J= E (9.11)
me νc
The Ohm’s law states
J = σ0 E (9.12)
the DC conductivity is thus
ne e2
σ0 = . (9.13)
me νc
CHAPTER 9. PLASMA CONDUCTIVITY AND DIFFUSION 38
J=S ·E (9.14)
Jx σ⊥ −σH 0 Ex
J y = σH σ⊥ 0 Ey ,
Jz 0 0 σk Ez
where the components of the DC conductivity tensor S are
νc2
σ⊥ = σ0 (9.15)
νc2 + Ω2ce
νc Ωce
σH = 2 σ0 (9.16)
νc + Ω2ce
ne e 2
σk = σ0 = , (9.17)
me νc
where νc is the collision frequency for momentum transfer between electrons
and heavy particles and Ωce is the electron cyclotron frequency due to the
external magnetic field. Find the components of the mobility tensor Me
defined as
ue = Me · E. (9.18)
Results:
M⊥ −MH 0
Me = MH M⊥ 0
0 0 Mk
νc e
M⊥ = − (9.19)
me (νc2 + Ω2ce )
Ωce e
MH = − (9.20)
me (νc2 + Ω2ce )
e
Mk = − (9.21)
me νc
∂n(x, t) ∂ 2 n(x, t)
=D (9.22)
∂t ∂x2
by separation of variables, assuming
Results:
• Tk (t) = T0 exp(−D k 2 t)
10.1 Theory
In a paper from 1923, an American chemist and physicist wrote, that
electrons are repelaed from negative electrode, whereas positive ions are
attracted towards it. Langmuir concluded, that around every negative elec-
trode, a sheath of defined thickness containing only positive ions and neutral
atoms exists. Morever, Langmuir observed, that also the glass wall of the
discharge chamber is negatively charged and repels (or reflects) almost all
electrons [16].
The fact, that insulated objects inside plasma are negatively charged
(in respect to plasma) to floating potential, is caused by higher mobility of
electrons than ions. The thermal velocity of electrons (kB Te /me )1/2 is at
least 100 times higher than the thermal velocity of ions (kB Ti /Mi )1/2 [17].
The first reason for different mobility is higher mass of ions. If we consider
only proton (the lightest ion that can appear in a plasma), than the ratio
between the mass of proton and electron mp /me is 1836. This ratio cor-
responds approximately to the ratio of masses of heavy bowling ball (5 kg)
and ping pong ball (2,7 g). Another reason for higher thermal velocity of
electrons in low-temperature plasma is their higher temperature in respect
to the ions.
The slowest possible speed of ions at the sheath edge is called the Bohm
speed uB . Ionts are accelerated to this speed in a quasi-neutral pre-sheath,
where small electric field exists. The Bohm criterion of plasma sheath is
described by following equation
r
kB Te
us (0) ≥ uB = . (10.1)
Mi
40
CHAPTER 10. SOME BASIC PLASMA PHENOMENA 41
10.2 Problems
10.2.1 Waves in non-magnetized plasma
Plasma of so called E layer of Earth’s ionosphere has electron density ap-
proximately 105 cm−3 and is at altitude of approximatelyo 100 km.
a) Which electromagnetic waves can be reflected from this layer?
b) Calculate the dielectric constant of plasma for the waves with frequencies
of 100 Mhz and 1000 Hz.
b) Calculate the skin depth of the wave with frequency of 1000 Hz.
Solution:
a) All electromagnetic waves with frequency lower than the plasma frequency
(2 839 725 Hz) will be reflected.
ωp2
ε=1− 2
ω
For 100 MHz ε = 0.9991 (positive value, elmag. waves propagate), for
1000 Hz ε = −8064037 (negative value, imaginary refraction index, reflec-
tion). q
c) Skin depth δ appropriately equals to c/ ωp2 − ω 2 , where c is speed of
light. The skin depth for 1000 Hz is 16.8 m.
Boltzmann and
Fokker-Planck collision terms
11.1 Theory
Under several simplifying assumptions (mainly homogeneous and isotropic
distribution function of electronic velocities, molecular chaos, considering
only binary collisions and ignoring external forces), so called Boltzmann
collision integral can be derived
ZZ
∂f
= g σ(g, Ω) [fe (v0 ) f1 (v0 1 ) − fe (v) f1 (v1 )] dΩ d3 v. (11.1)
∂t coll
43
CHAPTER 11. BOLTZMANN AND FOKKER-PLANCK COLLISION TERMS44
where Z Z
h∆vi iav = ∆vi g σ(Ω)dΩfβ1 d3 v1 (11.3)
Ω v1
Z Z
h∆vi ∆vj iav = ∆vi ∆vj g σ(Ω)dΩfβ1 d3 v1 (11.4)
Ω v1
11.2 Problems
11.2.1 Collisions for Maxwell-Boltzmann distribution func-
tion
Consider a plasma in which the electrons and the ions are characterised,
respectively, by the following distribution functions
3/2
me (v − ue )2
me
fe = n0 exp − (11.5)
2 π k Te 2 k Te
3/2
mi (v − ui )2
mi
fi = n0 exp − (11.6)
2 π k Ti 2 k Ti
(a) Calculate the difference (fe (v0 ) fi (v10 ) − fe (v) fi (v1 )).
(b) Show that this plasma of electrons and ions will be in the equilibrium
state, that is, the difference (fe (v0 ) fi (v10 ) − fe (v) fi (v1 )) will vanish if
and only if ue = ui and Te = Ti .
Solution
(a)
3
me mi 3/2
0 1
(fe (v ) fi (v10 ) − fe (v) fi (v1 )) = n20 ×
2πk Te Ti
me (v0 − ue )2 mi (v10 − ui )2
× exp − − −
2 k Te 2 k Ti
me (v − ue )2 mi (v1 − ui )2
− exp − − (11.7)
2 k Te 2 k Ti
(b) For the difference to vanish, the term in the parentheses must equal zero.
This will happen if the arguments of the exponentials will be equal. Let
us rewrite the arguments, omitting the factor −(2 k)−1 :
me 02 mi 02
(v − 2 v0 · ue + u2e ) + (v − 2 v10 · ui + u2i ) (11.8)
Te Ti 1
CHAPTER 11. BOLTZMANN AND FOKKER-PLANCK COLLISION TERMS45
me 2 mi 2
(v − 2 v · ue + u2e ) + (v − 2 v1 · ui + u2i ) (11.9)
Te Ti 1
From the derivation of the Boltzmann collision term follows, that the
pairs of velocities v, v1 and v0 , v10 can be considered as pairs of velocities
before and after an elastic two-body collision. Thus, they are bound by
the conservation laws:
me v 2 + mi v12 me v 02 + mi v102
= (11.10)
2 2
me v + mi v1 = me v0 + mi v10 (11.11)
It is now obvious from the last four equations, that the collision term
will vanish if and only if Te = Ti and ue = ui . In other words, the
distribution function fe will be changed by the collisions only if the
plasma is out of equilibrium – the collisions tend to bring the plasma to
the state of equilibrium.
mi (v10 − ui )2
fi = Ci exp − (11.13)
2 k Ti
Will the difference (fe (v0 ) fi (v10 ) − fe (v) fi (v1 )) be zero for ue = ui ?
[8] Eun Ha Choi, Jeong Chull Ahn, Min Wook Moon, Jin Goo Kim,
Myung Chul Choi, Choon Gon Ryu, Sung Hyuk Choi, Tae Seung Cho,
Yoon Jung, Guang Sup Cho, et al. Electron temperature and plasma
density in surface-discharged alternating-current plasma display panels.
Plasma Science, IEEE Transactions on, 30(6):2160–2164, 2002.
46
BIBLIOGRAPHY 47
pressure helium and argon rf plasma jets. Plasma Sources Science and
Technology, 20(6):065010, 2011.
[11] J Hubert, M Moisan, and A Ricard. A new microwave plasma at at-
mospheric pressure. Spectrochimica Acta Part B: Atomic Spectroscopy,
34(1):1–10, 1979.
[12] Zdeněk Hubička. The low temperature plasma jet sputtering systems
applied for the deposition of thin films. 2012.
[13] Umran S Inan and Marek Golkowski. Principles of plasma physics for
engineers and scientists. Cambridge University Press, 2010.
[14] Jae Duk Kim, Young Ho Na, Young June Hong, Han Sup Uhm,
Eun Ha Choi, et al. Microwave plasma jet system development at
atmospheric pressure using a 2.45 ghz gan hemt devices. In Plasma
Science (ICOPS), 2011 Abstracts IEEE International Conference on,
pages 1–1. IEEE, 2011.
[15] SB Krupanidhi and M Sayer. Position and pressure effects in rf mag-
netron reactive sputter deposition of piezoelectric zinc oxide. Journal
of applied physics, 56(11):3308–3318, 1984.
[16] I. Langmuir. Positive Ion Currents from the Positive Column of Mer-
cury Arcs. Science, 58:290–291, October 1923.
[17] M.A. Lieberman and A.J. Lichtenberg. Principles of plasma discharges
and materials processing. Published by A Wiley-Interscience Publica-
tion, page 388, 1994.
[18] Liming Liu, Ruisheng Huang, Gang Song, and Xinfeng Hao. Be-
havior and spectrum analysis of welding arc in low-power yag-laser–
mag hybrid-welding process. Plasma Science, IEEE Transactions on,
36(4):1937–1943, 2008.
[19] Yasuyuki Noguchi, Akira Matsuoka, Kiichiro Uchino, and Katsunori
Muraoka. Direct measurement of electron density and temperature
distributions in a micro-discharge plasma for a plasma display panel.
Journal of applied physics, 91(2):613–616, 2002.
[20] Kunio Okimura, Akira Shibata, Naohiro Maeda, Kunihide Tachibana,
Youichiro Noguchi, and Kouzou Tsuchida. Preparation of rutile tio2
films by rf magnetron sputtering. JAPANESE JOURNAL OF AP-
PLIED PHYSICS PART 1 REGULAR PAPERS SHORT NOTES
AND REVIEW PAPERS, 34:4950–4950, 1995.
[21] Cheng-gang PAN, Xue-ming HUA, Wang ZHANG, Fang LI, and Xiao
XIAO. Calculating the stark broadening of welding arc spectra by
fourier transform method. 32(7), 2012.
BIBLIOGRAPHY 48