Fokker-Planck collision operator
|
|
- Karin Allen
- 5 years ago
- Views:
Transcription
1 DRAFT 1 Fokker-Planck collision operator Felix I. Parra Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford OX1 3NP, UK (This version is of 16 April 18) 1. Introduction In these notes, we study the effect on the plasma of close encounters between charged particles (electrons and ions). These close encounters are the collisions between charge particles that we have ignored in Collisionless Plasma Physics. We can naively estimate how close the particles have to be to notice one another. The characteristic potential energy of a system of two charged particles with charges of the order of the proton charge e separated by a distance b is e /4πɛ b, where ɛ is the vacuum permittivity. For the particles to be affected by their mutual interaction, this potential energy must be of the order of the typical kinetic energy of the particles, that is, of the order of the temperature of the plasma T. Thus, the distance b between two particles that have a significant collision is b = e 4πɛ T. (1.1) The naive estimation in (1.1) is incorrect in weakly coupled plasmas. To define weakly coupled plasmas, we need to define the coupling parameter: the ratio of the characteristic potential energy between two typical particles in the plasma and the characteristic kinetic energy of the plasma. The typical distance between two particles is n 1/3, where n is the particle density (number of particles per unit volume). Then, the potential energy between two typical particles is e n 1/3 /4πɛ, and the coupling parameter is Γ = e n 1/3 4πɛ T. (1.) By definition, in weakly coupled plasmas, the coupling parameter is small, Γ 1. (1.3) As a result, particles are rarely within a distance b of each other, that is, the probablity of finding a particle within a sphere of radius b is very small, 4π 3 nb3 Γ 3 1. (1.4) Importantly, the electric force on a given particle is dominated by particles that are at a distance of the order of the Debye length, ɛ T λ D = e n. (1.5) In weakly coupled plasmas, the Debye length is much larger than the impact parameter b. Indeed, the ratio between the two, also known as plasma parameter, is Λ = λ D b 1 1. (1.6) Γ3/
2 Felix I. Parra Thus, the electric force on a particle exerted by another particle at a Debye length is much smaller than the force exerted by another particle at a distance b. Even though the force due to each individual particle is small, the collective force due to all the particles within a Debye length is comparable to (and as we will see, larger than) the force exerted by the small number of particles that on rare occasions happen to be at a distance b. Indeed, the number of particles within a Debye sphere is very large, The total potential energy due to all these particles is e 4π 3 nλ3 D 1 1. (1.7) Γ3/ 4π 3 nλ3 D 4π 1 4πɛ λ D 3 nλ3 D Λ 4πɛ b e T, (1.8) 4πɛ b Thus, the collective potential energy due to all the particles within the Debye sphere is comparable to the potential energy of a particle that is sufficiently close (a distance b) to modify the kinetic energy significantly. We will show that in fact, the force due to the particles within the Debye sphere dominates over the close encounters. This is in contrast to the collisions considered in the kinetic theory of neutral gases. Collisions in neutral gases are dominated by close encounters, and the collective force of particles far from the particle of interest is negligible. Based on these considerations, we derive the collision operator for charged particles, known as Fokker-Planck collision operator, starting from the Boltzmann operator for binary collisions derived in the Kinetic Theory (Dellar 15). As we have indicated above, the collisional events are not composed of just two particles interacting with each other, but we can still use a binary collision operator. We can consider the individual interaction of the particle of interest with every other particle within a Debye length, and then using the fact that these interactions are small and hence the effects are simply additive, we can sum over all of them. e. Boltzmann collision operator Before we give the Boltzmann collision operator, we discuss binary collisions in detail..1. Binary collisions In a binary collision between particle 1 of species s and particle of species s, we need to determine the final velocities of both particles, v 1f and v f, from their initial velocities, v 1i and v i (see figure 1). For simplicity, we only consider particles that interact via a potential V (r) that only depends on the relative distance r = r 1 r between the position of the two particles. The potential V (r) vanishes for r. The equations of motion for particles 1 and are d r 1 dt = r 1 V = dv dr d r dt = r V = dv dr r 1 r r 1 r, (.1) r r 1 r r 1. (.) Assuming a potential that only depends on the separation between the particles simplifies the problem considerably. There are three important simplifications: (a) Conservation of momentum. Adding equations (.1) and (.), we find that
3 Fokker-Planck collision operator 3 Collision v 1i v 1f v i v f Figure 1. Schematic collision between particles 1 and. the total momentum of both particles is conserved during the collision, that is, d dt (v 1 + v ) =. (.3) This property indicates that the center of mass velocity v CM = v 1 + v + (.4) is constant. This is a convenient result because it means that we can write both velocities v 1 and v as functions of a constant, v CM, and the difference between the two velocities Indeed, g = v 1 v. (.5) v 1 = v CM + g, + v = v CM g. (.6) + Thus, we just need to determine g f = v 1f v f from the initial condition g i = v 1i v i. To find an equation for g(t), we subtract equation (.) divided by from equation (.1) divided by. The result is an equation for the separation between the particles r = r 1 r : the reduced particle equation where d r dt = 1 dv r µ ss dr r, (.7) µ ss = (.8) + is the reduced mass. Once r(t) is calculated using (.7), we can find g = dr/dt by taking a time derivative. (b) Planar motion. The separation between particles r(t) and the velocity difference g(t) remain within a plane during the collision. By taking the cross product of (.7) with µ ss r, we find that the reduced particle s angular momentum L = µ ss r g (.9) is a constant of the motion, dl/dt =. Since L is a constant and L r =, the separation r and the velocity difference g = dr/dt are in the plane perpendicular to L for all time t. The collision characteristics are independent of the orientation of the plane L r = because equation (.7) is the same independently of the plane in which it is projected.
4 4 Felix I. Parra (c) Elastic collisions. Since we have assumed that the force is derived from a potential, the collision is elastic, that is, the total kinetic energy is the same before and after the collision. To show this, we can take the scalar product of (.7) with µ ss g = µ ss (dr/dt) and use (dv/dr)(r/ r ) = r V to obtain that the reduced particle s total energy E = 1 µ ss g + V (r) (.1) is a constant of the motion, de/dt =. Since V (r) for r, the magnitudes of g i and g f are the same, and only the direction is different. The properties described above imply that we just need to determine the change of direction of the velocity difference g due to a collision. Since the plane of the collision is mostly irrelevant, the angle χ between the initial velocity difference g i and final velocity difference g f is the important parameter. We use an orthonormal basis {ˆx i, ŷ i, g i /g i } aligned with the initial velocity difference g i (see figure ). Note that ˆx i and ŷ i are only determined up to a rotation since we only know that they must be in the plane perpendicular to g i. In the basis {ˆx i, ŷ i, g i /g i }, the final velocity difference g f is determined by two angles, χ and ϕ, g f = g i g i sin χ(cos ϕ ˆx i + sin ϕ ŷ i ) + cos χ g i. (.11) The angle ϕ determines the plane in which the collision takes place and χ gives the deflection within that plane. Note that we have defined the operator Q that relates quantities before and after the collision. Instead of the two angles χ and ϕ, we can determine the direction of g f using the unit vector ĉ (see figure ), ĉ = g f g i g f g i = cos ( χ With this unit vector we can define the operator Q as ) ( χ ) gi (cos ϕ ˆx i + sin ϕ ŷ i ) sin. (.1) g i g f = g i (I ĉĉ) g i, (.13) The advantage of using the vector ĉ is that it shows that the operator Q is its own inverse, and hence it gives the conditions before the collision if we know the conditions after, Indeed, g i = g f (I ĉĉ) g f. (.14) g i = (I ĉĉ) (I ĉĉ) g i = [I 4ĉĉ + 4(ĉ ĉ)ĉĉ g i = g i. (.15) Once ĉ or the angles χ and ϕ are given, we can obtain the final g f, and using (.6), determine v 1f and v f, The relation (.16) is its own inverse, v 1f = v 1i v 1i ĉ(ĉ g i ), + v f = v i v i + ĉ(ĉ g i ). (.16) + v 1i = v 1f v 1f ĉ(ĉ g f ), + v i = v f v f + ĉ(ĉ g f ). (.17) +
5 (a) g i Fokker-Planck collision operator 5 (b) ĉ g f (g i ĉ)ĉ g i g f ' ŷ i ˆx i Figure. (a) 3D sketch of the relative position of the final velocity difference g f with respect to the initial velocity difference g i. The angles χ and ϕ and the vector ĉ are sketched. (b) Relative position of the final velocity difference g f with respect to the initial velocity g i in the plane of the collision. The projection (g i ĉ)ĉ is sketched. 1 r 1 g i Figure 3. Impact parameter and its relation to the angle χ. Two impact parameters, ρ 1 (black) and ρ (blue), are given with their respective trajectories and angles, χ 1 and χ. This figure assumes a repelling force that increases when the particles gets closer. We have again used the operator Q that gives the final (initial) velocities for a collision whose vector ĉ and initial (final) velocities are known... Boltzmann collision operator The collision operator C ss gives the rate of change in time of the distribution function of species s, f s (r, v, t), as a result of collisions with species s. The kinetic equation for the distribution function f s (r, v, t) in a collisional plasma is f s t + v f s + Z se (E + v B) v f s = s C ss. (.18) The Boltzmann operator for binary collisions between species s and s was derived using the BBGKY approach in (Dellar 15). Here we use a heuristic derivation. The collision operator is composed of two terms, C ss (r, v, t) d 3 v d 3 r = (collisions/time with v 1f within volume d 3 v around v) (collisions/time with v 1i within volume d 3 v around v). (.19) To estimate the number of collisions per unit time, we use the construction in figure 4. The impact parameter ρ is the distance of closest approach between the two particles if
6 6 Felix I. Parra ˆx i ' g i t d' d g i Location of particle of species s with velocity v 1i ŷ i Figure 4. Sketch of the volume that a particle of species s and velocity within the infinitesimal volume d 3 v i around v i must occupy during the time interval t to pass at a distance [ρ, ρ + dρ and at an angle [ϕ, ϕ + dϕ of a particle of species s with velocity v 1i. The time t is assumed to be sufficiently long that g i t ρ. As a result, we can ignore particles that are a distance ρ ahead of the particle of species s with velocity v 1i even though they interacting with the particle of species s with velocity v 1i. we assume that they do not interact. The impact parameter determines the angle χ that characterizes the collision (see figure 3). The number of collisions characterized by the angles [ϕ, ϕ + dϕ and [χ, χ + dχ of a particle of species s with velocity v 1i in the time interval t with particles of species s with velocity within the infinitesimal volume d 3 v i around v i corresponds to the number of particles of species s with velocity within the infinitesimal volume d 3 v i around v i that pass at a distance [ρ(χ) dρ, ρ(χ) (with dρ = ρ/ χ dχ) and an angle [ϕ, ϕ + dϕ of a particle of species s with velocity v 1i in the time interval t. The number of particles of species s that satisfy these conditions is (density of particles of species s within the volume d 3 v i around v i ) (volume sketched in figure 4) = (f s (r, v i, t) d 3 v i )(g i t ρ dρ dϕ). (.) Dividing by t, multiplying by the number of particles of species s with velocity within the infinitesimal volume d 3 v 1i around v 1i, f s (r, v 1i, t) d 3 v 1i d 3 r, and integrating over all possible initial velocities of particles of species s, and over all possible impact parameters ρ and angles ϕ, we obtain the number of collisions per unit time of particles of species s with initial velocity within an infinitesimal volume d 3 v 1i around v i, collisions/time with initial velocity within an infinitesimal volume d 3 v 1i around v 1i π = d 3 r d 3 v 1i d 3 v i dρ dϕ f s (r, v 1i, t)f s (r, v i, t)g i ρ. (.1) Before continuing, it is convenient to define the differential cross section σ ss (g i, χ). The differential cross section can be calculated if we know the relation between the
7 Fokker-Planck collision operator 7 impact parameter ρ and the angle χ, ρ(χ) (see figure 3), σ ss = ρ ρ sin χ χ. (.) For collisions satisfying equations (.1)-(.), the differential cross section can only depend on the magnitude of the initial velocity difference, g i, and the angle χ due to the symmetry of the equations (see (.7)). With the definition (.), equation (.1) becomes collisions/time with initial velocity within an infinitesimal volume d 3 v 1i of v 1i π π = d 3 r d 3 v 1i d 3 v i dχ dϕ f s (r, v 1i, t)f s (r, v i, t)g i σ ss (g i, χ) sin χ. (.3) Using (.3), we can obtain the two terms in (.19). Since the operator Q gives the initial conditions for some given final conditions, we use it in (.3) to rewrite the first term in (.19) as collisions/time with v 1f within volume d 3 v of v π π = d 3 r d 3 v d 3 v dχ dϕ f s (r, v, t)f s (r, v, t)gσ ss (g, χ) sin χ, (.4) where v i = v, and g = g i = g. Using (.3), the second term in (.19) simply becomes collisions/time with v 1i within volume d 3 v of v π π = d 3 r d 3 v d 3 v dχ dϕ f s (r, v, t)f s (r, v, t)gσ ss (g, χ) sin χ. With these results, equation (.19) can be rewritten as C ss (r, v, t) d 3 v d 3 r = d 3 r d 3 v d 3 r d 3 v π d 3 v dχ π d 3 v dχ π π (.5) dϕ f s (r, v, t)f s (r, v, t)gσ ss (g, χ) sin χ dϕ f s (r, v, t)f s (r, v, t)gσ ss (g, χ) sin χ. (.6) We can simplify further using the fact that the operator Q is its own inverse. The determinant of the Jacobian of the operator Q is then ±1, giving d 3 v d 3 v = d 3 v d 3 v. Using this result, equation (.6) finally becomes π C ss [f s, f s (r, v, t) = d 3 v dχ π [ dϕ gσ ss (g, χ) sin χ f s (r, v, t)f s (r, v, t) f s (r, v, t)f s (r, v, t). (.7) To deduce the Fokker-Planck collision operator, we will not use (.7). Instead, it is convenient to use moments of the collision operator. For any function X(r, v, t), we
8 8 Felix I. Parra calculate the corresponding moment X(r, v, t) C ss [f s, f s (r, v, t) d 3 v = d 3 v X(v) d 3 v X(v) π d 3 v dϕ π π d 3 v dϕ π dχ gσ ss (g, χ) sin χf s (v)f s (v ) dχ gσ ss (g, χ) sin χf s (v)f s (v ), (.8) where we have omitted the dependence on r and t in the right side for brevity. Using the fact that the operator Q is its own inverse (and hence d 3 v d 3 v = d 3 v d 3 v ), and the fact that g = g, the first integral in (.8) can be written as π π d 3 v X(v) d 3 v dϕ dχ gσ ss (g, χ) sin χf s (v)f s (v ) = π d 3 v X(v) d 3 v dϕ π dχ gσ ss (g, χ) sin χf s (v)f s (v ). (.9) Using (.9), and changing the dummy integration variables v and v to v and v, equation (.8) becomes X(r, v, t) C ss [f s, f s (r, v, t) d 3 v = d 3 v f s (v) d 3 v f s (v ) 3. Coulomb collisions π dϕ π dχ gσ ss (g, χ) sin χ [X(v) X(v). (.3) We need to calculate the differential cross section for Coulomb interactions. To do so, we first calculate the relation between the impact parameter ρ and the angle χ. The equations of motion are (.1)-(.) with the potential V (r) = Z sz s e 4πɛ r. (3.1) As we have already explained, the motion throughout the collision is described by equation (.7) for the reduced particle. To solve this equation, we recall that the motion remains within the same plane, and we use the polar coordinates r and θ in figure 3. We also use that the reduced particle angular momentum is conserved, L(t) µ ss r dθ dt = L(t ) µ ss ρg, (3.) and that the reduced particle total energy is constant, E(t) 1 ( ) dr µ ss + 1 ( ) dθ dt µ ss r + Z sz s e = E(t ) 1 dt 4πɛ r µ ss g. (3.3) Solving for dθ/dt from (3.) and for dr/dt from (3.3), we find an equation for r(θ), dr dθ = dr/dt dθ/dt = ±r 1 b ss ρ ρ r r, (3.4) where b ss = Z sz se 4πɛ µ ss g. (3.5)
9 Fokker-Planck collision operator 9 Due to the symmetry of the problem (see figure 3), the solution from θ = π/ χ/ to θ = π χ is the reflection of the solution from θ = to θ = π/ χ/. The sign of dr/dθ determines whether the equation is for the portion of the orbit in which r decreases as θ increases (in figure 3, from r = at θ = to r = r min at θ = π/ χ/) or it is for the section in which r increases as θ increases (in figure 3, from r = r min at θ = π/ χ/ to r = at θ = π χ). We integrate (3.4) from θ = to θ = π/ χ/, This integral gives rmin ρ dr r 1 b ss /r ρ /r = π/ χ/ dθ. (3.6) [ ( ) r=rmin ρ/r + b ss /ρ arcsin = π 1 + b ss /ρ χ. (3.7) The minimum r is given by setting dr/dθ in (3.4) equal to zero, ρ r min = 1 + b ss /ρ b ss /ρ. (3.8) Using this result in (3.7), we find π arcsin ( r= b ss /ρ 1 + b ss /ρ ) = π χ, (3.9) and it can be rewritten as b ss ρ = tan(χ/). (3.1) Using (3.1) in (.), we obtain the differential cross section for Coulomb collisions, also known as Rutherford cross section, σ ss = b ss 4 sin 4 (χ/). (3.11) 4. Fokker-Planck collision operator We have argued in the introduction that the interaction of a given charged particle with other charged particles is dominated by small collisions with all the particles within the Debye sphere around the particle of interest. For a weakly coupled plasma, particles in the Debye sphere are mostly at a distance ρ λ D b ss. As a result, equation (3.1) implies that ( ) bss χ = arctan b ss 1 1. (4.1) ρ ρ Λ Then, we can assume χ 1 and use the approximation σ ss 4b ss χ 4. (4.) for the Rutherford cross section. The fact that χ is small can be used to simplify the operator Q in (.16). For example, the final velocity of a particle of species s that started with velocity v after a collision with a particle of species s with velocity v becomes v v + v + O(χ 3 v), (4.3)
10 1 Felix I. Parra where v = + [χg (cos ϕ ˆx i + sin ϕ ŷ i ) χ g. (4.4) We have kept corrections of order χ because we will need them due to several cancellations. Using equation (4.3) for particles of species s and a similar expression for particles of species s, we can expand the Boltzmann operator in (.7). Performing the expansion directly is tedious, but there is a shortcut that simplifies the calculation. We will use the expression for the moments of the Boltzman collision operator given in (.3). We are going to look for an approximate operator Css FP that satisfies X(r, v, t) Css FP [f s, f s (r, v, t) d 3 v X(r, v, t) C ss [f s, f s (r, v, t) d 3 v (4.5) for any functions X(r, v, t), f s (r, v, t) and f s (r, v, t). This operator must then satisfy Using (4.3), we obtain C FP ss [f s, f s C ss [f s, f s. (4.6) X(v) X(v) = v v X + 1 v v : v v X + O(χ 3 X). (4.7) Using this approximation in (.3), we find X(r, v, t) C ss [f s, f s (r, v, t) d 3 v = d 3 v f s (v) d 3 v f s (v ) π + χ v v : v v X + O(χ 4 X) [ ( χsmall dϕ dχ gσ ss (g, χ) χ v v X ) π + ( ) X(v) X(v) dχ gσ ss (g, χ) sin χ χ } small {{} gb ss X. (4.8) Note that we have split the integral over χ into two intervals: the interval [, χ small, with χ small 1, in which the approximation χ 1 is valid, and the interval [χ small, π that accounts for the collisions with large deflection angles. Using (4.) and (4.4), we take the integrals over the angles ϕ and χ in the region [, χ small, finding π χsmall dϕ dχ gσ ss (g, χ) χ v = 4πm s gb χsmall ss dχ g (4.9) + χ and π χsmall dϕ dχ gσ ss (g, χ) χ v v = πm s g3 b ss ( + ) (ˆx iˆx i + ŷ i ŷ i ) χsmall ( ) 1 χ + O(1) dχ. (4.1) Both of these integrals diverge due to the lower limit χ =. The Coulomb interaction is long range, and it may seem that as a result, every charged particle interacts with all the other charged particles in the system. Fortunately, there is a scale, the Debye length,
11 Fokker-Planck collision operator 11 λ D, beyond which the Coulomb potential is shielded. To obtain the Boltzmann collision operator, we had to assume that we can consider the interaction between charged particles as many uncorrelated binary collisions. Particles beyond the Debye length respond to the motion of the charged particle of interest in a correlated manner. We can only use the Boltzmann collision operator for ρ λ D. Then, according to equation (4.1), the Botzmann operator can only be used for collisions that satisfy χ b ss /λ D. For χ of order b ss /λ D or smaller, we need to take into account the correlations between particles. It can be done rigorously using the BBGKY formalism, and the result is the Balescu- Lenard collision operator (Balescu 196; Lenard 196; Hazeltine & Waelbroeck 4). The Boltzmann collision operator and the Balescu-Lenard operator can be combined into a single collision operator (Frieman & Book 1963) that is rarely used. Using the fact that the Botzmann operator can only be used for collisions that satisfy χ b ss /λ D, the integral χ small χ 1 dχ that appears in (4.9) and (4.1) should be χsmall b ss /λ D dχ χ = ln ( ) λd χ small = ln Λ ss 1, (4.11) b ss where we taken χ small / to be larger than b ss /λ D. The function ln Λ ss is the Coulomb logarithm. We have used the notation ln Λ ss even though, as we will see shortly, the quantity λ D χ small /b ss will not be the plasma parameter Λ defined in (1.6). For ln Λ ss 1, the collision operator for charged particles is greatly simplified. Using (4.11) in (4.9) and (4.1), we obtain and π dϕ χsmall dχ gσ ss (g, χ) χ v = 4π gb ln Λ ss ss b ss /λ D + π χsmall dϕ dχ gσ ss (g, χ) χ b ss /λ D v v g (4.1) πm s g3 b ss ln Λ ss ( + ) (ˆx iˆx i + ŷ i ŷ i ). (4.13) The term due to χ [χ small, π in (4.8) can be neglected as small in 1/ ln Λ ss 1 compared to the terms in (4.1) and (4.13). The contribution from collisions with χ b ss /λ D, included in the Balescu-Lenard collision operator, turns out to be small in 1/ ln Λ ss 1 as well. Thus, as long as 1/ ln Λ ss 1, we can neglect the collisions that do not satisfy b ss /λ D χ χ small. With the results in (4.9) and (4.1), and using ˆx iˆx i + ŷ i ŷ i = I gg/g, equation (4.8) becomes X(r, v, t) C ss [f s, f s (r, v, t) d 3 v where γ ss + 1 g I gg g 3 d 3 v f s (v) : v v X d 3 v f s (v ) [ ( ) g g 3 vx, (4.14) γ ss = πz s Z s e4 ln Λ ss (4πɛ ). (4.15)
12 1 Felix I. Parra To write the right side of (4.14) as X(r, v, t) C FP ss [f s, f s (r, v, t) d 3 v, we integrate by parts in velocity space to find f s (v) g g 3 vx d 3 v = X v (f s (v) gg ) 3 d 3 v (4.16) and f s (v) g I gg g 3 : v v X d 3 v = With these results, equation (4.14) becomes (4.5) with { f s (v )[ C FP ss [f s, f s = γ ss v This is the Fokker-Planck collision operator. [ ( ) X v v f s (v) g I gg g 3 d 3 v. (4.17) 1 v ( g I gg g 3 ( ) g g 3 f s(v) ) f s (v) d 3 v }. (4.18) 5. The Coulomb logarithm The Fokker-Planck collision operator is based on the expansion in 1/ ln Λ ss 1. To find ln Λ ss, we need equations for b ss /λ D and χ small /. Fortunately, both b ss /λ D and χ small / appear inside logarithm. As long as ln Λ ss 1, we do not need to worry about factors of order unity in b ss /λ D and χ small /. (a) Conventionally (Huba 13), the Debye length is defined to be 1/ Zs λ D = e n s, (5.1) ɛ T s Ts / g where the summation is done only over species whose thermal v ts = T s / is faster or of the order of the typical velocity difference g T s / + T s /. The idea is that species with thermal speeds much slower than g cannot respond in the time scale of the collision and will not be able to shield the electric field. Note that the formula in (5.1) is defined such that having more species reduces the Debye length due to additional shielding. (b) The impact parameter b ss depends on the relative velocity g. Since b ss appears within the logarithm, we can replace g by an approximate average value, 3Ts 3Ts g g = +. (5.) The approximate average value of b ss calculated using g is usually employed (Huba 13), b ss b ss = Z sz s e ( + ) 1πɛ ( T s + T s ). (5.3) (c) The limit χ small was imposed in (4.8) so that the χ 1 approximations were valid, and it would be natural then to take χ small / = 1. However, the approximations that we used could fail for reasons different than χ being large. If the plasma is sufficiently hot,
13 the characteristic quantum wavelength, Fokker-Planck collision operator 13 l ss = µ ss g, (5.4) with = h/π and h Planck s constant, becomes larger than the impact parameter b ss, defined in (3.5), and according to (4.1), quantum effects start to matter for small angles χ b ss /l ss 1. If we take χ small / to be ( χ small = min 1, b ) ss, (5.5) l ss and ln Λ ss = ln(λ D χ small /b ss ) is still larger than one, we can ignore the small angle collisions with quantum corrections by saying that the are small in 1/ ln Λ ss 1. The angle χ small depends on g, but it appears inside the logarithm. Thus, we can use the approximation in (5.) to find χ small χsmall ( = min 1, b ss ), (5.6) l ss where l ss = ( + ) 3 ( T s + T s ). (5.7) Thus, conventionally (Huba 13), we define the Coulomb logarithm to be ( ( ) ( )) λd λd ln Λ ss = min ln, ln, (5.8) b ss l ss where λ D, b ss and l ss are defined in (5.1), (5.3) and (5.7), respectively. As an example, we can calculate ln Λ ss for electron-electron, electron-ion and ion-ion collisions. For electron-electron collisions, since m i m e, only the electrons are sufficiently fast to satisfy the condition v te g v ti, giving λ D = ɛ T e /e n e. The impact parameter is b ee = e /1πɛ T e, and the characteristic quantum wavelength is l ee = / 6m e T e. Then, { 18 (1/) ln(ne [1 ln Λ ee = m 3 ) + (3/) ln(t e [1 kev) for T e ev 16 (1/) ln(n e [1 m 3 ) + ln(t e [1 kev) for T e ev. (5.9) For electron-ion collisions, the electrons are the only species that has a thermal speed v te g v ti, giving λ D = ɛ T e /e n e. Neglecting m e m i, we find that the impact parameter is b ei = Ze /1πɛ T e, and the characteristic quantum wavelength is l ei = / 3m e T e. Then, ln Λ ei = { 18 (1/) ln(ne [1 m 3 ) ln Z + (3/) ln(t e [1 kev) for T e 4Z ev 16 (1/) ln(n e [1 m 3 ) + ln(t e [1 kev) for T e 4Z ev. (5.1) Finally, for ion-ion collisions, both electron and ions satisfy v ts g, giving the Debye length λ D = (e n e /ɛ T e + i Z i e n i /ɛ T i ) 1/. Then, lnλ ii = [ 18 ln ln A it i +A i T i Z iz i (A i+a i ) [ AiA i (A it i +A i T i) A i+a i ( n e T e + ) Z 1/ i i n i T i ( n e T e + Z i i n i T i ) 1/ for Ti A i + T i A 7Z i i Z i kev, for Ti A i + T i A 7Z i i Z i kev (5.11)
14 14 Felix I. Parra where the densities n s must be given in 1 m 3, the temperatures T s in kev, and A s = /m p is the mass of species s divided by the mass of the proton. 6. Landau form of the Fokker-Planck collision operator The Fokker-Planck collision operator can be written in a form that makes the symmetry between species s and s explicit: the Landau form. This form is very convenient to show that the Fokker-Planck collision operator conserves particles, momentum and energy and that it has an H-theorem. To obtain the Landau form, we first need to prove two useful relationships, and g g g = g I gg g 3 (6.1) g g g = g g 3. (6.) We proceed to prove them. In an orthonormal basis {ˆx, ŷ, ẑ}, g = g xˆx + g y ŷ + g z ẑ and the gradient g is defined to be Then, for g = g = g f = ˆx f + ŷ f + ẑ f. (6.3) g x g y g z gx + gy + gz, we obtain g g = g xˆx + g y ŷ + g z ẑ gx + gy + gz = g g. (6.4) Taking a second gradient of this expression, we find ( ) g g g g = g = ( ) gg 1 + g g = I g g g g ( gg)g g = I g gg g 3. (6.5) This expression proves (6.1). Taking the trace of (6.1), we obtain ( g ) I gg gg = Trace( g g g) = Trace g 3 = 3 g g g g 3 = g. (6.6) Taking a gradient of this expression, we find ( ) g gg = g = gg g g This expression proves (6.). From equation (6.), we deduce the useful relationships = g g 3. (6.7) v ( g g g) = v ( g g g) = g g 3. (6.8) We proceed to prove them. Using Einstein s repeated index notation and the chain rule, we obtain [ v ( g g g) i = ( ) g = g k 3 g (6.9) v j g j g i v j g k g j g i and [ v ( g g g) i = ( ) g v j = g k 3 g g j g i v j. (6.1) g k g j g i
15 Fokker-Planck collision operator 15 Since g = v v, we find that g k / v j = v k / v j = δ jk and g k / v j = v k / v j = δ jk, where δ jk is Kronecker s delta. Using these expressions in equations (6.9) and (6.1), we obtain and [ v ( g g g) i = 3 g g j g j g i = ( g g g) i (6.11) [ v ( g g g) i = 3 g g j g j g i = ( g g g) i. (6.1) These expressions and (6.) give (6.8). Using (6.1) and (6.8), we can rewrite (4.18) as { f s (v )[ C ss [f s, f s = γ ss v Here we can simplify 1 v ( g g gf s (v)) 1 v ( g g g)f s (v) + 1 v ( g g g)f s (v) d 3 v }, (6.13) v ( g g gf s (v)) v ( g g g)f s (v) = g g g v f s (v). (6.14) We can also integrate by parts to find f s (v ) v ( g g g) d 3 v = g g g v f s (v ) d 3 v. (6.15) Using (6.14) and (6.15) in (6.13), we obtain the Landau form of the Fokker-Planck collision operator { [ } C ss [f s, f s = γ ss f s (v ) v g g g v f s (v) f s(v) v f s (v ) d 3 v. (6.16) 7. Conservation properties of the Fokker-Planck collision operator The Fokker-Planck collision operator conserves particles, momentum and energy. We proceed to show these properties using the Landau form in (6.16) Conservation of particles The Landau collision operator satisfies C ss [f s, f s (r, v, t) d 3 v =. (7.1) This property is a consequence of the fact that the collision operator in (6.16) is a divergence and the fact that the distribution functions f s and f s vanish for large velocities. 7.. Conservation of momentum The friction force on species s due to collision with s is given by F ss = v C ss [f s, f s (r, v, t) d 3 v. (7.)
16 16 Felix I. Parra The Fokker-Planck collision operator conserves momentum, that is, the force of species s on species s is equal and opposite to the force exerted by species s on s, F ss + F s s =. (7.3) For collisions between particle of the same species s, this property implies that F ss =, (7.4) that is, there is no net collisional force of species s on itself. We proceed to prove (7.3). We use the Landau form in (7.) that can be written as where [ Γ ss = g g g C ss [f s, f s = γ ss v Γ ss, (7.5) f s (v ) v f s (v) f s(v) v f s (v ) d 3 v. (7.6) Substituting (7.5) into (7.) and integrating by parts, we obtain F ss = γ ss v( v Γ ss ) d 3 v = γ ss Γ ss v v d 3 v = γ ss Γ ss d 3 v [ = γ ss d 3 v d 3 v f s (v ) g g g v f s (v) f s(v) v f s (v ). (7.7) Exchanging s and s and recalling that γ ss = γ s s, we find [ F s s = γ ss d 3 v d 3 v f s (v ) g g g v f s (v) f s (v) v f s (v ). (7.8) In this equation, the integration variables v and v are dummy variables. We can exchange their names, finding [ F s s = γ ss d 3 v d 3 v f s (v) g g g v f s (v ) f s (v ) v f s (v). (7.9) Summing (7.7) and (7.9), we obtain the momentum conservation equation in (7.3) Conservation of energy The energy gained or lost per unit time by species s due to collision with s is 1 W ss = v C ss [f s, f s (r, v, t) d 3 v. (7.1) The Fokker-Planck collision operator conserves energy, W ss + W s s =. (7.11) For collisions between particle of the same species s, this property implies that W ss =, (7.1) that is, there is no net collisional energy gain or loss due to collision of species s with itself. We proceed to prove (7.11). Substituting (7.5) into (7.1) and integrating by parts, we
17 Fokker-Planck collision operator 17 obtain W ss = γ ss v ( v Γ ss ) d 3 v = γ ss Γ ss v v d 3 v = γ ss v Γ ss d 3 v [ = γ ss d 3 v d 3 v f s (v ) v g g g v f s (v) f s(v) v f s (v ). (7.13) Exchanging s and s and the dummy integration variables v and v, we find [ W s s = γ ss d 3 v d 3 v v f s (v) g g g v f s (v ) f s (v ) v f s (v). (7.14) Summing (7.13) and (7.14), we obtain W ss + W s s = γ ss d 3 v From equation (6.1), we obtain d 3 v g g g g g g g g = g g I gg g 3 proving the energy conservation equation in (7.11). [ f s (v ) v f s (v) f s(v) v f s (v ). (7.15) = g g g g g =, (7.16) g3 8. H-theorem for the Fokker-Planck collision operator The Fokker-Planck collision operator satisfies its own H-theorem, that is, it always produces entropy. The species s entropy production due to collisions with species s is σ ss = ln f s (r, v, t) C ss [f s, f s (r, v, t) d 3 v. (8.1) The H-theorem for the Fokker-Planck collision operator states that σ ss + σ s s, (8.) and σ ss + σ s s is equal to zero only when both f s and f s are Maxwellians with the same average velocity u and temperature T, ( ) 3/ ms f s (r, v, t) = n s (r, t) exp ( v u(r, t) ), πt (r, t) T (r, t) ( ) 3/ ms f s (r, v, t) = n s (r, t) exp ( v u(r, t) ). (8.3) πt (r, t) T (r, t) Note that densities n s and n s can be different, and that densities, average flow and temperature can be general functions of position r and time t. For collisions between particles of the same species s, the H-theoreimplifies to σ ss, (8.4) and σ ss is equal to zero only when f s is a Maxwellian, ( ) 3/ f s (r, v, t) = n s (r, t) exp ( v u s (r, t) ). (8.5) πt s (r, t) T s (r, t)
18 18 Felix I. Parra In this case, the temperature T s and the average velocity u s do not have to be equal to the temperature or average velocity of any of the other species. We first prove (8.), and we then show that σ ss + σ s s = only for Maxwellian distribution functions. Substituting (7.5) into (8.1) and integrating by parts, we obtain σ ss = γ ss ln f s ( v Γ ss ) d 3 v = γ ss Γ ss v ln f s d 3 v = γ ss d 3 v d 3 v v ln f s (v) g g g [ f s (v ) v f s (v) f s(v) v f s (v ). (8.6) Exchanging s and s and the dummy integration variables v and v, we find σ s s = γ [ ss d 3 v d 3 v v ln f s (v f s (v) ) g g g v f s (v ) f s (v ) v f s (v). (8.7) Summing (8.6) and (8.7), we obtain σ ss + σ s s = γ ss d 3 v d 3 v f s (v)f s (v ) a(v, v ) g g g a(v, v ), (8.8) where a(v, v ) = 1 v ln f s (v) 1 v ln f s (v ). (8.9) From equation (6.1), we obtain a g g g a = 1 g a a g g g, (8.1) proving that the entropy production in (8.8) is always positive. We proceed to consider the case in which the entropy production in (8.8) vanishes. Given that the Fokker-Planck collision operator in (6.16) contains second derivatives with respect to velocity of the distribution functions f s and f s, we assume that the distribution functions are continuous and that their derivatives with respect to the velocities are also continuous. According to (8.1), σ ss + σ s s = only if a = (a g/g )g for all v and v (recall that the derivatives with respect to velocity of f s and f s are continuous, and hence a(v, v ) is continuous). Then, for σ ss + σ s s =, the vector a, defined in (8.9) must be parallel to g, giving 1 v ln f s (v) 1 v ln f s (v ) = K(v, v )g, (8.11) where the function K(v, v ) is unknown. Taking v = v in this equation, we obtain that Using this result in (8.11), we find Taking v = in (8.13), we find 1 v ln f s (v) = 1 v ln f s (v). (8.1) v ln f s (v) v ln f s (v ) = K(v, v )g. (8.13) v ln f s (v) = K(v, )v + v ln f s (). (8.14)
19 Substituting equation (8.14) into (8.13), we obtain Fokker-Planck collision operator 19 K(v, )v K(v, )v = K(v, v )(v v ). (8.15) When v and v are linearly independent, this equation implies that K(v, ) = K(v, v ) and that K(v, ) = K(v, v ), leading to K(v, ) = K(v, ). Thus, K(v, ) and hence K(v, v ) are constants. This argument only breaks when v and v are colinear, but due to continuity, K(v, v ) must be a constant in this case as well. The constant K must be negative because v ln f s in (8.14) must be negative for large v to ensure that distribution functions vanish at v. We name the constant 1/T, 1 T = K(v, ) = K(v, v ). (8.16) Using this result and naming the value of v ln f s at v = u T = v ln f s (), (8.17) equation (8.14) becomes Using (8.1), we then obtain v ln f s (v) = (v u). (8.18) T v ln f s (v) = (v u). (8.19) T Integrating equations (8.18) and (8.19) in velocity space gives (8.3). 9. Solutions to the Fokker-Planck collision operator The H-theorem of the Fokker-Planck collision operator implies that the Maxwellians given in (8.3) and in (8.5) are the only solutions to the Fokker-Planck collision operator. The only solutions to the system of equations C ss [f s, f s =, (9.1) C s s[f s, f s =, (9.) are the Maxwellians in (8.3). To prove this, multiply equation (9.1) by ln f s and integrate over velocity space, and multiply equation (9.) by ln f s and integrate over velocity space. These equations give σ ss = = σ s s. Then, according to the H-theorem, f s and f s must be the Maxwellians in (8.3). Similarly, we can show that the only solution to C ss [f s, f s = (9.3) is the Maxwellian in (8.5). The H-theorem is even more powerful. It implies that the steady-state solution of a closed system can only be Maxwellians with the same temperature and velocity for all species. The kinetic equation is f s t + v f s + Z se (E + v B) v f s = s C ss [f s, f s. (9.4) Multiplying this equation by ln f s, we obtain [ t + v + Z se (E + v B) v ( f s ln f s + f s ) = ln f s C ss [f s, f s. (9.5) s
20 Felix I. Parra Using phase space volume conservation, v + v [(Z s e/ )(E + v B) =, we can rewrite the equation as [ t ( f Zs e s ln f s + f s ) + [v( f s ln f s + f s ) + v (E + v B)( f s ln f s + f s ) Integrating over velocity space and the volume V of the plasma, we find = s ln f s C ss [f s, f s. (9.6) closed system d 3 r d 3 v ( f s ln f s + f s ) + d t V S d 3 v v ˆn( f s ln f s + f s ) = d 3 r σ ss, V s V (9.7) where V is the surface that limits the volume, and ˆn is the normal to that surface. Summing over species, we finally obtain Thus, in steady state, s t V d 3 r d 3 v ( f s ln f s + f s ) = s,s s,s V V d 3 r σ ss. (9.8) d 3 r σ ss =. (9.9) Since σ ss + σ s s, the sum of all these entropy productions can only vanish if each one is independently zero. This implies that the only allowed solution are Maxwellians with the same temperature and average velocity. In most cases, the system is not closed, and it is driven by sources and sinks that prevent the distribution functions from becoming exact Maxwellians. REFERENCES Balescu, R. 196 Irreversible processes in ionized gases. Phys. Fluids 3, 5. Dellar, P. 15 Lecture notes for Part 1 of MMathPhys Kinetic Theory. Frieman, E.A. & Book, D.L Convergent classical kinetic equation for a plasma. Phys. Fluids 6, 17. Hazeltine, R.D. & Waelbroeck, F.L. 4 The Framework of Plasma Physics. Westview Press. Huba, J.D. 13 NRL Plasma Formulary. Lenard, A. 196 On Bogoliubov s kinetic equation for a spatially homogeneous plasma. Annals of Physics 3, 39.
Lectures on basic plasma physics: Introduction
Lectures on basic plasma physics: Introduction Department of applied physics, Aalto University Compiled: January 13, 2016 Definition of a plasma Layout 1 Definition of a plasma 2 Basic plasma parameters
More informationTransport coefficients in plasmas spanning weak to strong correlation
Transport coefficients in plasmas spanning weak to strong correlation Scott D. Baalrud 1,2 and Jerome Daligault 1 1 Theoretical Division, Los Alamos National Laboratory 2 Department of Physics and Astronomy,
More informationWave propagation in an inhomogeneous plasma
DRAFT Wave propagation in an inhomogeneous plasma Felix I. Parra Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford OX NP, UK This version is of 7 February 208. Introduction In
More informationChapter 3. Coulomb collisions
Chapter 3 Coulomb collisions Coulomb collisions are long-range scattering events between charged particles due to the mutual exchange of the Coulomb force. Where do they occur, and why they are of interest?
More informationPLASMA: WHAT IT IS, HOW TO MAKE IT AND HOW TO HOLD IT. Felix I. Parra Rudolf Peierls Centre for Theoretical Physics, University of Oxford
1 PLASMA: WHAT IT IS, HOW TO MAKE IT AND HOW TO HOLD IT Felix I. Parra Rudolf Peierls Centre for Theoretical Physics, University of Oxford 2 Overview Why do we need plasmas? For fusion, among other things
More informationLong range interaction --- between ions/electrons and ions/electrons; Coulomb 1/r Intermediate range interaction --- between
Collisional Processes Long range interaction --- between ions/electrons and ions/electrons; Coulomb 1/r Intermediate range interaction --- between ions/electrons and neutral atoms/molecules; Induced dipole
More information8.333: Statistical Mechanics I Fall 2007 Test 2 Review Problems
8.333: Statistical Mechanics I Fall 007 Test Review Problems The second in-class test will take place on Wednesday 10/4/07 from :30 to 4:00 pm. There will be a recitation with test review on Monday 10//07.
More informationFigure 1.1: Ionization and Recombination
Chapter 1 Introduction 1.1 What is a Plasma? 1.1.1 An ionized gas A plasma is a gas in which an important fraction of the atoms is ionized, so that the electrons and ions are separately free. When does
More information12. MHD Approximation.
Phys780: Plasma Physics Lecture 12. MHD approximation. 1 12. MHD Approximation. ([3], p. 169-183) The kinetic equation for the distribution function f( v, r, t) provides the most complete and universal
More informationAPPENDIX Z. USEFUL FORMULAS 1. Appendix Z. Useful Formulas. DRAFT 13:41 June 30, 2006 c J.D Callen, Fundamentals of Plasma Physics
APPENDIX Z. USEFUL FORMULAS 1 Appendix Z Useful Formulas APPENDIX Z. USEFUL FORMULAS 2 Key Vector Relations A B = B A, A B = B A, A A = 0, A B C) = A B) C A B C) = B A C) C A B), bac-cab rule A B) C D)
More informationIntroduction. Chapter Plasma: definitions
Chapter 1 Introduction 1.1 Plasma: definitions A plasma is a quasi-neutral gas of charged and neutral particles which exhibits collective behaviour. An equivalent, alternative definition: A plasma is a
More informationKinetic theory of ions in the magnetic presheath
Kinetic theory of ions in the magnetic presheath Alessandro Geraldini 1,2, Felix I. Parra 1,2, Fulvio Militello 2 1. Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford, Oxford
More informationPlasma collisions and conductivity
e ion conductivity Plasma collisions and conductivity Collisions in weakly and fully ionized plasmas Electric conductivity in non-magnetized and magnetized plasmas Collision frequencies In weakly ionized
More informationResistive MHD, reconnection and resistive tearing modes
DRAFT 1 Resistive MHD, reconnection and resistive tearing modes Felix I. Parra Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford OX1 3NP, UK (This version is of 6 May 18 1. Introduction
More informationCandidacy Exam Department of Physics February 6, 2010 Part I
Candidacy Exam Department of Physics February 6, 2010 Part I Instructions: ˆ The following problems are intended to probe your understanding of basic physical principles. When answering each question,
More information14. Energy transport.
Phys780: Plasma Physics Lecture 14. Energy transport. 1 14. Energy transport. Chapman-Enskog theory. ([8], p.51-75) We derive macroscopic properties of plasma by calculating moments of the kinetic equation
More informationCollision Processes. n n The solution is 0 exp x/ mfp
Collision Processes Collisions mediate the transfer of energy and momentum between various species in a plasma, and as we shall see later, allow a treatment of highly ionized plasma as a single conducting
More informationPlasmas as fluids. S.M.Lea. January 2007
Plasmas as fluids S.M.Lea January 2007 So far we have considered a plasma as a set of non intereacting particles, each following its own path in the electric and magnetic fields. Now we want to consider
More informationThe 3 dimensional Schrödinger Equation
Chapter 6 The 3 dimensional Schrödinger Equation 6.1 Angular Momentum To study how angular momentum is represented in quantum mechanics we start by reviewing the classical vector of orbital angular momentum
More informationPlasma Astrophysics Chapter 1: Basic Concepts of Plasma. Yosuke Mizuno Institute of Astronomy National Tsing-Hua University
Plasma Astrophysics Chapter 1: Basic Concepts of Plasma Yosuke Mizuno Institute of Astronomy National Tsing-Hua University What is a Plasma? A plasma is a quasi-neutral gas consisting of positive and negative
More informationPlasma Physics Prof. V. K. Tripathi Department of Physics Indian Institute of Technology, Delhi
Plasma Physics Prof. V. K. Tripathi Department of Physics Indian Institute of Technology, Delhi Module No. # 01 Lecture No. # 02 Plasma Response to Fields: Fluid Equations Well, friends today I am going
More informationElectromagnetic Field Theory (EMT)
Electromagnetic Field Theory (EMT) Lecture # 9 1) Coulomb s Law and Field Intensity 2) Electric Fields Due to Continuous Charge Distributions Line Charge Surface Charge Volume Charge Coulomb's Law Coulomb's
More informationd 1 µ 2 Θ = 0. (4.1) consider first the case of m = 0 where there is no azimuthal dependence on the angle φ.
4 Legendre Functions In order to investigate the solutions of Legendre s differential equation d ( µ ) dθ ] ] + l(l + ) m dµ dµ µ Θ = 0. (4.) consider first the case of m = 0 where there is no azimuthal
More informationSpin Stability of Aysmmetrically Charged Plasma Dust. I.H. Hutchinson. February 2004
PSFC/JA-04-3 Spin Stability of Aysmmetrically Charged Plasma Dust I.H. Hutchinson February 2004 Plasma Science and Fusion Center Massachusetts Institute of Technology Cambridge, MA 02139 USA This work
More informationFluid equations, magnetohydrodynamics
Fluid equations, magnetohydrodynamics Multi-fluid theory Equation of state Single-fluid theory Generalised Ohm s law Magnetic tension and plasma beta Stationarity and equilibria Validity of magnetohydrodynamics
More informationMultiple Integrals and Vector Calculus (Oxford Physics) Synopsis and Problem Sets; Hilary 2015
Multiple Integrals and Vector Calculus (Oxford Physics) Ramin Golestanian Synopsis and Problem Sets; Hilary 215 The outline of the material, which will be covered in 14 lectures, is as follows: 1. Introduction
More informationLecture 14 (11/1/06) Charged-Particle Interactions: Stopping Power, Collisions and Ionization
22.101 Applied Nuclear Physics (Fall 2006) Lecture 14 (11/1/06) Charged-Particle Interactions: Stopping Power, Collisions and Ionization References: R. D. Evans, The Atomic Nucleus (McGraw-Hill, New York,
More informationElastic Collisions. Chapter Center of Mass Frame
Chapter 11 Elastic Collisions 11.1 Center of Mass Frame A collision or scattering event is said to be elastic if it results in no change in the internal state of any of the particles involved. Thus, no
More informationTHE UNIVERSITY OF SYDNEY FACULTY OF SCIENCE INTERMEDIATE PHYSICS PHYS 2912 PHYSICS 2B (ADVANCED) ALL QUESTIONS HAVE THE VALUE SHOWN
CC0936 THE UNIVERSITY OF SYDNEY FACULTY OF SCIENCE INTERMEDIATE PHYSICS PHYS 91 PHYSICS B (ADVANCED) SEMESTER, 014 TIME ALLOWED: 3 HOURS ALL QUESTIONS HAVE THE VALUE SHOWN INSTRUCTIONS: This paper consists
More informationCharged particle motion in external fields
Chapter 2 Charged particle motion in external fields A (fully ionized) plasma contains a very large number of particles. In general, their motion can only be studied statistically, taking appropriate averages.
More informationTHE UNIVERSITY OF SYDNEY FACULTY OF SCIENCE INTERMEDIATE PHYSICS PHYS 2912 PHYSICS 2B (ADVANCED) ALL QUESTIONS HAVE THE VALUE SHOWN
CC0936 THE UNIVERSITY OF SYDNEY FACULTY OF SCIENCE INTERMEDIATE PHYSICS PHYS 91 PHYSICS B (ADVANCED SEMESTER, 015 TIME ALLOWED: 3 HOURS ALL QUESTIONS HAVE THE VALUE SHOWN INSTRUCTIONS: This paper consists
More informationPhysics 20 Homework 3 SIMS 2016
Physics 20 Homework 3 SIMS 2016 Due: Thursday, August 25 th Special thanks to Sebastian Fischetti for problems 1, 5, and 6. Edits in red made by Keith Fratus. 1. The ballistic pendulum is a device used
More informationDielectrics. Lecture 20: Electromagnetic Theory. Professor D. K. Ghosh, Physics Department, I.I.T., Bombay
What are dielectrics? Dielectrics Lecture 20: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay So far we have been discussing electrostatics in either vacuum or in a conductor.
More informationMonte Carlo Collisions in Particle in Cell simulations
Monte Carlo Collisions in Particle in Cell simulations Konstantin Matyash, Ralf Schneider HGF-Junior research group COMAS : Study of effects on materials in contact with plasma, either with fusion or low-temperature
More informationSingle Particle Motion
Single Particle Motion C ontents Uniform E and B E = - guiding centers Definition of guiding center E gravitation Non Uniform B 'grad B' drift, B B Curvature drift Grad -B drift, B B invariance of µ. Magnetic
More information1 The Kompaneets equation
Introduction to Particle Cosmology, Tutorial notes: The Sunyaev-Zeldovich effect 1 The Kompaneets equation Consider evolution of photon phase space density n(, t) in the presence of an electron gas. Assume
More informationPreliminary Quantum Questions
Preliminary Quantum Questions Thomas Ouldridge October 01 1. Certain quantities that appear in the theory of hydrogen have wider application in atomic physics: the Bohr radius a 0, the Rydberg constant
More informationKinetic Theory of the Presheath and the Bohm Criterion
UW-CPTC 10-3 Kinetic Theory of the Presheath and the Bohm Criterion S D Baalrud and C C Hegna Department of Engineering Physics, University of Wisconsin-Madison, 1500 Engineering Drive, Madison, WI 53706,
More informationMultiple Integrals and Vector Calculus: Synopsis
Multiple Integrals and Vector Calculus: Synopsis Hilary Term 28: 14 lectures. Steve Rawlings. 1. Vectors - recap of basic principles. Things which are (and are not) vectors. Differentiation and integration
More informationProblem Set #5: 5.7,5.9,5.13 (Due Monday, April 8th)
Chapter 5 Magnetostatics Problem Set #5: 5.7,5.9,5.13 (Due Monday, April 8th) 5.1 Biot-Savart Law So far we were concerned with static configuration of charges known as electrostatics. We now switch to
More informationPhysics 212: Statistical mechanics II Lecture IV
Physics 22: Statistical mechanics II Lecture IV Our program for kinetic theory in the last lecture and this lecture can be expressed in the following series of steps, from most exact and general to most
More informationLectures on basic plasma physics: Kinetic approach
Lectures on basic plasma physics: Kinetic approach Department of applied physics, Aalto University April 30, 2014 Motivation Layout 1 Motivation 2 Boltzmann equation (a nasty bastard) 3 Vlasov equation
More informationKinetic Solvers with Adaptive Mesh in Phase Space for Low- Temperature Plasmas
Kinetic Solvers with Adaptive Mesh in Phase Space for Low- Temperature Plasmas Vladimir Kolobov, a,b,1 Robert Arslanbekov a and Dmitry Levko a a CFD Research Corporation, Huntsville, AL 35806, USA b The
More informationElectrolytes. Chapter Basics = = 131 2[ ]. (c) From both of the above = = 120 8[
Chapter 1 Electrolytes 1.1 Basics Here we consider species that dissociate into positively and negatively charged species in solution. 1. Consider: 1 H (g) + 1 Cl (g) + ()+ () = { } = (+ )+ ( ) = 167[
More informationLecture 6: Ideal gas ensembles
Introduction Lecture 6: Ideal gas ensembles A simple, instructive and practical application of the equilibrium ensemble formalisms of the previous lecture concerns an ideal gas. Such a physical system
More information1. (3) Write Gauss Law in differential form. Explain the physical meaning.
Electrodynamics I Midterm Exam - Part A - Closed Book KSU 204/0/23 Name Electro Dynamic Instructions: Use SI units. Where appropriate, define all variables or symbols you use, in words. Try to tell about
More informationThe Distribution Function
The Distribution Function As we have seen before the distribution function (or phase-space density) f( x, v, t) d 3 x d 3 v gives a full description of the state of any collisionless system. Here f( x,
More informationChanges from 2013 NRL Plasma Formulary to the 2016 NRL Plasma Formulary. frictional force and electrical conductivities formulae changed to
Changes from 2013 NRL Plasma Formulary to the 2016 NRL Plasma Formulary p. 32 from to = 6.8 10 8 1/2 = 6.8 10 8 1/2 1+ 1+ ) 1/2 T 3/2 ) T 3/2 p. 37 frictional force and electrical conductivities formulae
More informationParticle collisions and decays
Particle collisions and decays Sourendu Gupta TIFR, Mumbai, India Classical Mechanics 2012 September 21, 2012 A particle decaying into two particles Assume that there is a reaction A B+C in which external
More informationρ c (2.1) = 0 (2.3) B = 0. (2.4) E + B
Chapter 2 Basic Plasma Properties 2.1 First Principles 2.1.1 Maxwell s Equations In general magnetic and electric fields are determined by Maxwell s equations, corresponding boundary conditions and the
More informationLecture Notes for PHY 405 Classical Mechanics
Lecture Notes for PHY 405 Classical Mechanics From Thorton & Marion s Classical Mechanics Prepared by Dr. Joseph M. Hahn Saint Mary s University Department of Astronomy & Physics September 1, 2005 Chapter
More information(a) Write down the total Hamiltonian of this system, including the spin degree of freedom of the electron, but neglecting spin-orbit interactions.
1. Quantum Mechanics (Spring 2007) Consider a hydrogen atom in a weak uniform magnetic field B = Bê z. (a) Write down the total Hamiltonian of this system, including the spin degree of freedom of the electron,
More informationSolution of time-dependent Boltzmann equation for electrons in non-thermal plasma
Solution of time-dependent Boltzmann equation for electrons in non-thermal plasma Z. Bonaventura, D. Trunec Department of Physical Electronics Faculty of Science Masaryk University Kotlářská 2, Brno, CZ-61137,
More informationErrata for Conquering the Physics GRE, edition 2
Errata for Conquering the Physics GRE, edition 2 March 12, 217 This document contains corrections to Conquering the Physics GRE, edition 2. It includes errata that are present in both all printings of
More informationIdeal gases. Asaf Pe er Classical ideal gas
Ideal gases Asaf Pe er 1 November 2, 213 1. Classical ideal gas A classical gas is generally referred to as a gas in which its molecules move freely in space; namely, the mean separation between the molecules
More informationElastic Scattering. R = m 1r 1 + m 2 r 2 m 1 + m 2. is the center of mass which is known to move with a constant velocity (see previous lectures):
Elastic Scattering In this section we will consider a problem of scattering of two particles obeying Newtonian mechanics. The problem of scattering can be viewed as a truncated version of dynamic problem
More informationCHAPTER 20. Collisions & Encounters of Collisionless Systems
CHAPTER 20 Collisions & Encounters of Collisionless Systems Consider an encounter between two collisionless N-body systems (i.e., dark matter halos or galaxies): a perturber P and a system S. Let q denote
More informationarxiv: v1 [cond-mat.stat-mech] 5 Mar 2013
Kinetic theory of spatially homogeneous systems with long-range interactions: III. Application to power-law potentials, plasmas, stellar systems, and to the HMF model arxiv:33.4v [cond-mat.stat-mech] 5
More informationPhysics Dec The Maxwell Velocity Distribution
Physics 301 7-Dec-2005 29-1 The Maxwell Velocity Distribution The beginning of chapter 14 covers some things we ve already discussed. Way back in lecture 6, we calculated the pressure for an ideal gas
More informationCHARGED PARTICLE MOTION IN CONSTANT AND UNIFORM ELECTROMAGNETIC FIELDS
CHARGED PARTICLE MOTION IN CONSTANT AND UNIFORM ELECTROMAGNETIC FIELDS In this and in the following two chapters we investigate the motion of charged particles in the presence of electric and magnetic
More informationOverview spherical accretion
Spherical accretion - AGN generates energy by accretion, i.e., capture of ambient matter in gravitational potential of black hole -Potential energy can be released as radiation, and (some of) this can
More informationPhysical models for plasmas II
Physical models for plasmas II Dr. L. Conde Dr. José M. Donoso Departamento de Física Aplicada. E.T.S. Ingenieros Aeronáuticos Universidad Politécnica de Madrid Physical models,... Plasma Kinetic Theory
More informationLecture: Scattering theory
Lecture: Scattering theory 30.05.2012 SS2012: Introduction to Nuclear and Particle Physics, Part 2 2 1 Part I: Scattering theory: Classical trajectoriest and cross-sections Quantum Scattering 2 I. Scattering
More informationGraduate Written Examination Spring 2014 Part I Thursday, January 16th, :00am to 1:00pm
Graduate Written Examination Spring 2014 Part I Thursday, January 16th, 2014 9:00am to 1:00pm University of Minnesota School of Physics and Astronomy Examination Instructions Part 1 of this exam consists
More informationH atom solution. 1 Introduction 2. 2 Coordinate system 2. 3 Variable separation 4
H atom solution Contents 1 Introduction 2 2 Coordinate system 2 3 Variable separation 4 4 Wavefunction solutions 6 4.1 Solution for Φ........................... 6 4.2 Solution for Θ...........................
More informationPhysics 3323, Fall 2016 Problem Set 2 due Sep 9, 2016
Physics 3323, Fall 26 Problem Set 2 due Sep 9, 26. What s my charge? A spherical region of radius R is filled with a charge distribution that gives rise to an electric field inside of the form E E /R 2
More informationChapter 4. Electrostatic Fields in Matter
Chapter 4. Electrostatic Fields in Matter 4.1. Polarization 4.2. The Field of a Polarized Object 4.3. The Electric Displacement 4.4. Linear Dielectrics 4.5. Energy in dielectric systems 4.6. Forces on
More informationAtomic Structure. Chapter 8
Atomic Structure Chapter 8 Overview To understand atomic structure requires understanding a special aspect of the electron - spin and its related magnetism - and properties of a collection of identical
More informationPart VIII. Interaction with Solids
I with Part VIII I with Solids 214 / 273 vs. long pulse is I with Traditional i physics (ICF ns lasers): heating and creation of long scale-length plasmas Laser reflected at critical density surface Fast
More informationIntroduction and Vectors Lecture 1
1 Introduction Introduction and Vectors Lecture 1 This is a course on classical Electromagnetism. It is the foundation for more advanced courses in modern physics. All physics of the modern era, from quantum
More informationPart II. Interaction with Single Atoms. Multiphoton Ionization Tunneling Ionization Ionization- Induced Defocusing High Harmonic Generation in Gases
- Part II 27 / 115 - 2-28 / 115 Bohr model recap. At the Bohr radius - a B = the electric field strength is: 2 me 2 = 5.3 10 9 cm, E a = e ab 2 (cgs) 5.1 10 9 Vm 1. This leads to the atomic intensity:
More informationElectric fields in matter
Electric fields in matter November 2, 25 Suppose we apply a constant electric field to a block of material. Then the charges that make up the matter are no longer in equilibrium: the electrons tend to
More information221B Lecture Notes Scattering Theory II
22B Lecture Notes Scattering Theory II Born Approximation Lippmann Schwinger equation ψ = φ + V ψ, () E H 0 + iɛ is an exact equation for the scattering problem, but it still is an equation to be solved
More informationAlternative method of developing Maxwell s fields
Apeiron, Vol. 14, No. 4, October 2007 464 Alternative method of developing Maxwell s fields Krzysztof Rȩbilas Zak lad Fizyki. Akademia Rolnicza im. Hugona Ko l l ataja Al. Mickiewicza 21, 31-120 Kraków.
More informationChapter 1 Nature of Plasma
Chapter 1 Nature of Plasma Abstract Charge neutrality is one of fundamental property of plasma. Section 1.2 explains Debye length λ D in (1.2), a measure of shielding distance of electrostatic potential,
More informationHeating and current drive: Radio Frequency
Heating and current drive: Radio Frequency Dr Ben Dudson Department of Physics, University of York Heslington, York YO10 5DD, UK 13 th February 2012 Dr Ben Dudson Magnetic Confinement Fusion (1 of 26)
More informationFokker-Planck Simulation Study of Hot-Tail Effect on Runaway Electron Generation in ITER Disruptions )
Fokker-Planck Simulation Study of Hot-Tail Effect on Runaway Electron Generation in ITER Disruptions ) Hideo NUGA, Akinobu MATSUYAMA 1), Masatoshi YAGI 1) and Atsushi FUKUYAMA Kyoto University, Kyoto-Daigaku-Katsura,
More informationProblem Set #3: 2.11, 2.15, 2.21, 2.26, 2.40, 2.42, 2.43, 2.46 (Due Thursday Feb. 27th)
Chapter Electrostatics Problem Set #3:.,.5,.,.6,.40,.4,.43,.46 (Due Thursday Feb. 7th). Coulomb s Law Coulomb showed experimentally that for two point charges the force is - proportional to each of the
More informationP3317 HW from Lecture 7+8 and Recitation 4
P3317 HW from Lecture 7+8 and Recitation 4 Due Friday Tuesday September 25 Problem 1. In class we argued that an ammonia atom in an electric field can be modeled by a two-level system, described by a Schrodinger
More informationConductors and Insulators
Conductors and Insulators Lecture 11: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay Self Energy of a Charge Distribution : In Lecture 1 we briefly discussed what we called
More informationCHAPTER 22. Astrophysical Gases
CHAPTER 22 Astrophysical Gases Most of the baryonic matter in the Universe is in a gaseous state, made up of 75% Hydrogen (H), 25% Helium (He) and only small amounts of other elements (called metals ).
More informationHomework Assignment 4 Solution Set
Homework Assignment 4 Solution Set PHYCS 442 7 February, 24 Problem (Griffiths 2.37 If the plates are sufficiently large the field near them does not depend on d. The field between the plates is zero (the
More informationSolid State Physics FREE ELECTRON MODEL. Lecture 17. A.H. Harker. Physics and Astronomy UCL
Solid State Physics FREE ELECTRON MODEL Lecture 17 A.H. Harker Physics and Astronomy UCL Magnetic Effects 6.7 Plasma Oscillations The picture of a free electron gas and a positive charge background offers
More information[variable] = units (or dimension) of variable.
Dimensional Analysis Zoe Wyatt wyatt.zoe@gmail.com with help from Emanuel Malek Understanding units usually makes physics much easier to understand. It also gives a good method of checking if an answer
More informationPHYS852 Quantum Mechanics II, Spring 2010 HOMEWORK ASSIGNMENT 8: Solutions. Topics covered: hydrogen fine structure
PHYS85 Quantum Mechanics II, Spring HOMEWORK ASSIGNMENT 8: Solutions Topics covered: hydrogen fine structure. [ pts] Let the Hamiltonian H depend on the parameter λ, so that H = H(λ). The eigenstates and
More informationComputing the Universe, Dynamics IV: Liouville, Boltzmann, Jeans [DRAFT]
Computing the Universe, Dynamics IV: Liouville, Boltzmann, Jeans [DRAFT] Martin Weinberg June 24, 1998 1 Liouville s Theorem Liouville s theorem is a fairly straightforward consequence of Hamiltonian dynamics.
More informationMathematical Preliminaries
Mathematical Preliminaries Introductory Course on Multiphysics Modelling TOMAZ G. ZIELIŃKI bluebox.ippt.pan.pl/ tzielins/ Table of Contents Vectors, tensors, and index notation. Generalization of the concept
More informationPhysics 142 Energy in Mechanics Page 1. Energy in Mechanics
Physics 4 Energy in Mechanics Page Energy in Mechanics This set of notes contains a brief review of the laws and theorems of Newtonian mechanics, and a longer section on energy, because of its central
More informationChapter 2. Dielectric Theories
Chapter Dielectric Theories . Dielectric Theories 1.1. Introduction Measurements of dielectric properties of materials is very important because it provide vital information regarding the material characteristics,
More informationLeast Squares Regression
Least Squares Regression Chemical Engineering 2450 - Numerical Methods Given N data points x i, y i, i 1 N, and a function that we wish to fit to these data points, fx, we define S as the sum of the squared
More informationSensors Plasma Diagnostics
Sensors Plasma Diagnostics Ken Gentle Physics Department Kenneth Gentle RLM 12.330 k.gentle@mail.utexas.edu NRL Formulary MIT Formulary www.psfc.mit.edu/library1/catalog/ reports/2010/11rr/11rr013/11rr013_full.pdf
More informationSommerfeld (1920) noted energy levels of Li deduced from spectroscopy looked like H, with slight adjustment of principal quantum number:
Spin. Historical Spectroscopy of Alkali atoms First expt. to suggest need for electron spin: observation of splitting of expected spectral lines for alkali atoms: i.e. expect one line based on analogy
More informationTheoretische Physik 2: Elektrodynamik (Prof. A-S. Smith) Home assignment 11
WiSe 22..23 Prof. Dr. A-S. Smith Dipl.-Phys. Matthias Saba am Lehrstuhl für Theoretische Physik I Department für Physik Friedrich-Alexander-Universität Erlangen-Nürnberg Problem. Theoretische Physik 2:
More informationThe Time Arrow of Spacetime Geometry
5 The Time Arrow of Spacetime Geometry In the framework of general relativity, gravity is a consequence of spacetime curvature. Its dynamical laws (Einstein s field equations) are again symmetric under
More informationFundamentals of Plasma Physics Transport in weakly ionized plasmas
Fundamentals of Plasma Physics Transport in weakly ionized plasmas APPLAuSE Instituto Superior Técnico Instituto de Plasmas e Fusão Nuclear Luís L Alves (based on Vasco Guerra s original slides) 1 As perguntas
More informationAccurate representation of velocity space using truncated Hermite expansions.
Accurate representation of velocity space using truncated Hermite expansions. Joseph Parker Oxford Centre for Collaborative Applied Mathematics Mathematical Institute, University of Oxford Wolfgang Pauli
More informationFluid Dynamics from Kinetic Equations
Fluid Dynamics from Kinetic Equations François Golse Université Paris 7 & IUF, Laboratoire J.-L. Lions golse@math.jussieu.fr & C. David Levermore University of Maryland, Dept. of Mathematics & IPST lvrmr@math.umd.edu
More informationChapter 18 Thermal Properties of Matter
Chapter 18 Thermal Properties of Matter In this section we define the thermodynamic state variables and their relationship to each other, called the equation of state. The system of interest (most of the
More informationFundamentals of Plasma Physics
Fundamentals of Plasma Physics Definition of Plasma: A gas with an ionized fraction (n i + + e ). Depending on density, E and B fields, there can be many regimes. Collisions and the Mean Free Path (mfp)
More information1 Fundamentals. 1.1 Overview. 1.2 Units: Physics 704 Spring 2018
Physics 704 Spring 2018 1 Fundamentals 1.1 Overview The objective of this course is: to determine and fields in various physical systems and the forces and/or torques resulting from them. The domain of
More information