Assume dilute, nondegenerate ideal electrons, protons and neutral atoms in thermal and chemical equilibrium, with charge neutrality np=ne. Set μH=μp+μe, using each ideal species' translational factor nλth3/g and hydrogen ground-state binding energy −R. Neglect excited states and take the heavy proton/atom thermal-wavelength ratio as one. With electron spin and neutral-atom degeneracies treated in the convention of the source, the remaining factor isDifferent partition-degeneracy conventions must be propagated consistently. Screening, dense-plasma level dissolution and departures from local thermodynamic equilibrium lie outside this derivation.
nHne2=(h22πmekBT)3/2e−R/(kBT).
For total nuclear density n, ne=nζ and nH=n(1−ζ). Saha's equation becomes nζ2/(1−ζ)=CT3/2e−R/(kBT). At ζ≪1,The factor two in the exponent follows from taking a square root; using the full ionization temperature would double the predicted logarithmic slope.
ζ≃C/nT3/4e−TH/T,m=43,TH=2kBR≃7.89×104 K.
Write u=1/T. Then logζ=constant−(3/4)logu−THu, giving local slope dlog10ζ/du=−(TH+3T/4)/log10. At 5000–6000 K its magnitude is about 3.6×104 K. The low-ionization tail in the original graph drops roughly one decade per 2.8×10−5 K−1, consistent with this value.
At u=2.0×10−4 K−1, the graph gives T≃5000 K and ζ≃0.004. Invert the full equation to get n=CT3/2e−R/(kBT)(1−ζ)/ζ2, yielding n∼1018 m−3. The neighboring tail markers give the same order; plotting precision and the squared dependence on the read ionization fraction preclude claiming many significant digits.
For weak potentials, Boltzmann electrons give δne=neeϕ/(kBT) against fixed ions. Poisson's equation becomes ∇2ϕ=ϕ/λD2 outside a test charge, giving a Yukawa potential with λD=ϵ0kBT/(nee2). For a small electron displacement x relative to fixed ions, charge separation produces restoring field E=neex/ϵ0, so mx¨=−nee2x/ϵ0. Hence ωp=nee2/(ϵ0me). The two formulas respectively encode thermal screening and electron inertia.
Choose ND=(4π/3)neλD3. The stipulated estimate then givesThus density cancels and conductivity grows as T3/2 at fixed ion charge. Physically faster electrons are less effectively deflected by Coulomb encounters. Numerical factors depend on the Debye-volume convention, and a realistic collision calculation introduces a Coulomb logarithm and evolving ionization; the cancellation is approximate rather than exact for every real plasma.
τ≃3Znee44πϵ02me(kBT)3/2,σ≃3Ze2me4πϵ02(kBT)3/2.
In the cold solid, phonon scattering makes τ∝T−1 and therefore conductivity falls as temperature rises, at roughly fixed carrier density. After melting and ionization, the character of scattering changes; a hot Coulomb plasma eventually follows an increasing T3/2 branch. A minimum in the crossover region is natural, but a weakly coupled plasma formula is not quantitatively justified throughout the warm dense regime.
Room temperature corresponds to about 0.026 eV, or 2.6×10−5 keV. Read the left-hand branch there and draw a horizontal line to the rising hot branch. The original graph gives a return temperature of order kBT′∼1–2 keV, with graph-reading uncertainty, not the temperature of minimum conductivity.
At full ionization ne=13ρ/(27mu)≃7.8×1029 m−3. Inserting the inferred temperature into ND=(4π/3)ne[ϵ0kBT′/(nee2)]3/2 gives a Debye-sphere occupancy of order 60–170. This is much greater than one and supports a collective weak-coupling description more than the conductivity-minimum region does. The approximate collision-time formula and omitted Coulomb logarithm do not warrant forcing the entire plotted conductivity curve to have its exact prefactor.
Take electrostatic perturbations proportional to ei(kz−ωt), with ions immobile on the fast electron time scale. Linearized Vlasov gives f1=−i(qE1/m)(∂f0/∂vz)/(ω−kvz). Integrating over transverse velocities gives F0(vz)=n0m/(2πkBT)e−mvz2/(2kBT). Combine the perturbed charge density with Poisson's equation:Since F0′=−mvzF0/(kBT), this is preciselyThe pole selects particles moving at phase speed. A causal initial-value prescription specifies how the contour passes it and produces resonant energy exchange, including Landau damping for a decreasing Maxwellian. A bare real improper integral without that prescription is incomplete. The unclosed brace in the printed Maxwellian exponent is typographical; the normalized Gaussian above fixes its meaning.
1+ϵ0mke2∫ω−kvzF0′(vz)dvz=0.
1−πϵ0mk2n0e2(2kBTm)3/2∫ω−kvzvze−mvz2/(2kBT)dvz=0.
When ∣kvz/ω∣≪1 over the populated velocities, expand the denominator in powers of kvz/ω. Odd Gaussian moments vanish. The first two nonzero terms use ⟨vz2⟩=vth2 and ⟨vz4⟩=3vth4, yielding 1−ωp2/ω2−3ωp2k2vth2/ω4+⋯=0. Solving consistently to order k2 givesThis is a long-wavelength/high-phase-speed approximation, not an exact Maxwellian dispersion relation at arbitrary k.
ω2≃ωp2+3k2vth2.
Differentiate the approximate dispersion: vg=3kvth2/ω, while vϕ=ω/k. Therefore vgvϕ=3vth2. In the high-phase-speed regime that justified the approximation, the group speed is correspondingly small. Trying to make vg much larger than vth would put the phase speed into the thermally populated distribution, where resonant Landau damping is strong and the undamped-fluid estimate becomes unreliable. This explains the stated tendency for these electron electrostatic waves, not a universal prohibition on fast electromagnetic waves in plasma.
Ancestors (9)
- 2016
- Past exam of the short options paper shared by Prelims and Parts A and B of the physics course of the University of Oxford
- Short options paper shared by Prelims and Parts A and B of the physics course of the University of Oxford
- Physics course of the University of Oxford
- Undergraduate course of the University of Oxford
- Course of the University of Oxford
- University of Oxford
- List of universities
- Home
Discussion (0)
New discussionThere are no discussions about this article yet.
Articles by others on the same topic (0)
There are currently no matching articles.