advanced31 min studyLesson 292 in Physics

The electron plasma oscillation frequency

Derive the restoring field of displaced electrons and distinguish angular frequency, oscillation period, and thermal corrections.

On this page

What you’ll understand

  • Derive the electron restoring force with the correct charge signs.
  • Convert angular frequency into ordinary frequency and period.
  • State when fixed-ion cold plasma oscillations are a useful approximation.

Displacing an electron population away from its neutralizing ions creates an electric field that pulls it back. The inertia of the electrons makes the response overshoot, producing a collective oscillation. This motion is not an isolated electron circling an ion. Its restoring field is generated by the coherent displacement of an entire population, so its frequency encodes density rather than an assumed microscopic orbital radius.

Derive the sign before naming the frequency

Imagine a broad uniform electron slab displaced a small distance ξ\xi toward positive x relative to a stationary positive background. The positive charge uncovered at the left and excess negative charge at the right produce an electric field toward positive x inside the overlapping region. Its magnitude is n0eξ/ϵ0n_0e\xi/\epsilon_0. An electron has charge −e-e, so its force points toward negative x: the displacement produces a restoring force rather than an instability.

Newton's law then has exactly the harmonic-oscillator form. Density is the equilibrium electron number per cubic metre, and the ions provide a fixed neutralizing background during the fast motion. The slab construction ignores edge fringing and large displacements. An equivalent local derivation combines linear electron continuity, momentum, and Poisson equations. It gives the same cold longitudinal frequency in a uniform unmagnetized plasma without requiring the whole apparatus to translate rigidly.

Ex=n0eϵ0ξ,ξ¨+ωpe2ξ=0,ωpe=n0e2ϵ0meE_x=\frac{n_0e}{\epsilon_0}\xi,\qquad \ddot\xi+\omega_{pe}^2\xi=0,\qquad \omega_{pe}=\sqrt{\frac{n_0e^2}{\epsilon_0m_e}}
These SI equations assume cold nonrelativistic electron motion, immobile ions, small displacement, and no collisions or applied magnetic field. The angular frequency ω_pe is in rad/s; the positive coefficient follows from the two electron-charge factors in the electric restoring force.↔ Scroll wide equations horizontally.

A timescale is not automatically a propagating wave

For initial displacement ξ0\xi_0 and zero initial velocity, displacement is a cosine. The electric field stores energy when charge separation is largest; electron bulk kinetic energy is largest as the slab passes through neutrality. In this ideal model the total quadratic energy is conserved. Zero field at an instant therefore does not mean that the disturbance has disappeared: its energy can temporarily reside entirely in coherent particle motion.

The cold electrostatic dispersion relation has no wave-number dependence: ω=ωpe\omega=\omega_{pe}. Consequently its group velocity vanishes in that idealization. A spatially patterned oscillation can exist without this model carrying a packet at a finite group speed. Finite temperature changes the dispersion; magnetic fields introduce additional branches. One must not use the cold electron formula as the frequency of every disturbance in a plasma, especially slow ion-supported waves.

ξ(t)=ξ0cos⁡(ωpet),fpe=ωpe2π,Tpe=2πωpe,λDe=vteωpe,vte=kBTeme\xi(t)=\xi_0\cos(\omega_{pe}t),\qquad f_{pe}=\frac{\omega_{pe}}{2\pi},\qquad \mathcal T_{pe}=\frac{2\pi}{\omega_{pe}},\qquad \lambda_{De}=\frac{v_{te}}{\omega_{pe}},\quad v_{te}=\sqrt{\frac{k_BT_e}{m_e}}
The cosine uses the specified zero initial velocity. Ordinary frequency f_pe is in Hz and period T_pe in seconds. The thermal speed convention uses no factor of two, which makes the stated electron Debye-length relation exact by definition. It does not add thermal dynamics to the cold oscillator.↔ Scroll wide equations horizontally.

Check the fixed-ion and collisionless windows

For a singly charged electron–ion plasma, the ion plasma frequency is smaller by me/mi\sqrt{m_e/m_i}. On a few electron periods, ion inertia therefore makes the fixed-background approximation reasonable. It is not a statement that ions never move. On longer times, their density response can change the mode structure and screening. If both species are mobile in a cold relative-displacement calculation, the squared relative-mode frequency contains both species' plasma-frequency squares.

Weak electron collisions require νe≪ωpe\nu_e\ll\omega_{pe} for many visible oscillations; neutral collisions can spoil this condition in a partially ionized gas. A warm Maxwellian high-phase-speed Langmuir wave has the approximate correction ω2≃ωpe2+3k2kBTe/me\omega^2\simeq\omega_{pe}^2+3k^2k_BT_e/m_e when kλDe≪1k\lambda_{De}\ll1. Collisionless damping depends on the velocity distribution and is not represented by the cold oscillator. A damping coefficient inserted for illustration must be distinguished from a kinetic derivation of that damping.

Density from an electron oscillation

For n0=1016 m−3n_0=10^{16}\,\mathrm{m^{-3}}, use e=1.60×10−19 Ce=1.60\times10^{-19}\,\mathrm C, me=9.11×10−31 kgm_e=9.11\times10^{-31}\,\mathrm{kg}, and ϵ0=8.85×10−12 F/m\epsilon_0=8.85\times10^{-12}\,\mathrm{F/m}. Find angular frequency, ordinary frequency, and period in the cold fixed-ion model.

  1. First compute n0e2/(ϵ0me)=3.175×1019 s−2n_0e^2/(\epsilon_0m_e)=3.175\times10^{19}\,\mathrm{s^{-2}}. The density appears to the first power in this squared frequency, so quadrupling density will only double the frequency.

  2. Taking the square root gives ωpe=5.63×109 rad/s\omega_{pe}=5.63\times10^9\,\mathrm{rad/s}. Divide by 2π2\pi to obtain fpe=8.97×108 Hzf_{pe}=8.97\times10^8\,\mathrm{Hz}, rather than calling the angular-frequency number a number of cycles per second.

  3. Invert the ordinary frequency: Tpe=1.12×10−9 s\mathcal T_{pe}=1.12\times10^{-9}\,\mathrm s. This is a nanosecond scale, which explains why ions can often be treated as stationary during the first few electron oscillations but laboratory electronics need substantial bandwidth to resolve them.

The result

The model predicts about 5.63 billion radians per second, 897 MHz, and a 1.12 ns period. These quantities describe one mode under its stated assumptions, not a collision rate.

Explore the model

Interactive lab

Interactive lab · oscillator

Isolate a restoring frequency

Reset and set both controls to damping ζ = 0 and natural frequency ω₀ = 1; predict a period of 2π. Reset and set both controls to ζ = 0 and ω₀ = 2; predict half that period, then Step or Play. The unit-mass model starts at displacement 1 with zero velocity and does not compute a plasma density or kinetic damping.

Isolate a restoring frequencyUnit mass, x(0) = 1 and v(0) = 0. The analytic solution includes underdamped, critically damped and overdamped cases; time is measured in seconds. Current readings: Displacement x(t): 1; Velocity v(t): 0; Regime: underdamped. All parameters are adjustable with the labeled controls below.02.557.510−1.15−0.5700.571.15Time · sDisplacementunderdamped
Displacement x(t)Current time
Displacement x(t)
1
Velocity v(t)
0
Regime
underdamped
t = 0 s

Unit mass, x(0) = 1 and v(0) = 0. The analytic solution includes underdamped, critically damped and overdamped cases; time is measured in seconds. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forA doubled restoring angular frequency halves the period; dimensional plasma parameters must be supplied outside this illustrative oscillator.

Check your understanding

If electron density increases by a factor of nine while the cold fixed-ion assumptions remain valid, what happens to the electron plasma period?

Try it yourself

A cold plasma diagnostic measures an electron oscillation frequency of 448.5 MHz under the same assumptions and constants as the worked example. Estimate the density and period. Explain why using 448.5 million directly as an angular frequency would misestimate density.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

Compare ordinary frequencies first; density scales with their squared ratio.

Compare with the solution

The measured frequency is half 897 MHz, giving one quarter of the reference density: n0≃2.50×1015 m−3n_0\simeq2.50\times10^{15}\,\mathrm{m^{-3}}. Its period is 1/(448.5×106)=2.23 ns1/(448.5\times10^6)=2.23\,\mathrm{ns}. In an absolute density calculation the angular frequency must be 2πf2\pi f. Omitting that factor makes the inferred density too small by (2π)2≃39.5(2\pi)^2\simeq39.5, because density is proportional to squared angular frequency. The comparison assumes the measured branch is genuinely the electron mode.

Take it with you

  • The electron charge sign makes the self-consistent field restoring.
  • Density sets squared plasma frequency; radians per second and Hz differ by 2π.
  • Ion response, thermal dispersion, and damping mark distinct limits of the cold model.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.