research36 min studyLesson 304 in Physics

Ion-acoustic waves and electron response

Combine fast electron shielding with ion inertia, derive finite-Debye-length dispersion, and state when kinetic damping invalidates an acoustic fluid picture.

On this page

What you’ll understand

  • Derive Boltzmann electron density with the correct potential sign.
  • Obtain ion-acoustic dispersion using cold ions and finite charge separation.
  • Identify the temperature and phase-speed ordering needed for weak damping.

An ion-acoustic wave is a collective density disturbance in which ion inertia works against electron pressure. Electrons respond rapidly along the disturbance and approximately maintain a Boltzmann density, while ions supply most of the oscillating mass. This separation can yield an acoustic branch without assuming frequent collisions. Its fluid-looking dispersion hides kinetic restrictions: the wave should outrun most ions while remaining slow compared with electrons if resonant damping is to stay weak.

Fast electrons need not be collisional

For an isothermal electron response with negligible inertia along the wave, electric force balances the electron pressure gradient. Integrating gives electron density proportional to the exponential of positive e times potential divided by electron thermal energy. The positive sign is physical: a positive potential attracts electrons and raises their density. Linearization requires potential energy small compared with thermal energy. The response is an approximation along connected electron trajectories, not a universal law in a sheath or trapped-electron population.

A collisionless kinetic electron population can produce this leading real response when phase speed is much less than its thermal speed. Its small imaginary response still carries electron Landau damping and is discarded by the simplest Boltzmann closure. Thus Boltzmann electrons do not imply that all electrons have no inertia in every situation, nor do they establish a collisional Maxwellian at every instant. The equilibrium distribution, connectivity, frequency, and perturbation amplitude determine when the closure is appropriate.

ne=n0exp⁡(eϕkBTe)≃n0(1+eϕkBTe),λDe2=ϵ0kBTen0e2n_e=n_0\exp\left(\frac{e\phi}{k_BT_e}\right)\simeq n_0\left(1+\frac{e\phi}{k_BT_e}\right),\qquad \lambda_{De}^2=\frac{\epsilon_0k_BT_e}{n_0e^2}
e > 0, φ is potential in volts, Te in kelvins, and electron Debye length λDe in metres. The linear form requires |eφ| ≪ kBTe. It describes a specified isothermal electron response, with kinetic corrections omitted.↔ Scroll wide equations horizontally.

Ion inertia closes the elementary dispersion

Take singly charged cold ions with no equilibrium flow. Linear continuity connects ion density to velocity, and the ion momentum equation connects velocity to electric acceleration. For a Fourier disturbance these give density proportional to potential times wavenumber squared divided by frequency squared. Combining that ion response with the Boltzmann electron response in Poisson yields an acoustic frequency divided by a Debye-length correction. Quasineutrality is therefore a controlled long-wavelength limit, rather than exact equality imposed before all field dynamics are considered.

At wavelengths long compared with the electron Debye length, phase speed approaches the electron-pressure sound speed. At shorter wavelengths, charge separation reduces the phase speed in this cold-ion model and frequency approaches the ion plasma frequency. That mathematical limit does not ensure that neglected kinetic or electron-inertia effects remain small for a chosen experiment. A warm-ion pressure term adds a closure-dependent correction; the familiar factor three corresponds to one-dimensional adiabatic compression in an appropriate high-phase-speed limit.

ωr2=k2cs21+k2λDe2,cs2=kBTemi,cs,warm2=kBTe+γikBTimi\omega_r^2=\frac{k^2c_s^2}{1+k^2\lambda_{De}^2},\qquad c_s^2=\frac{k_BT_e}{m_i},\qquad c_{s,\mathrm{warm}}^2=\frac{k_BT_e+\gamma_i k_BT_i}{m_i}
The first formula uses cold ions, Boltzmann electrons, and negligible damping. mi is in kilograms, cs in m/s, and k in m⁻¹. The warm expression additionally assumes a specified ion pressure closure; γi is not fixed without that assumption.↔ Scroll wide equations horizontally.

The kinetic ordering tests the closure

Define each species' one-component thermal speed as square root of thermal energy divided by mass. Weakly damped ion-acoustic propagation typically needs phase speed much greater than the ion thermal speed and much smaller than the electron thermal speed. The first requirement suppresses resonant ions, while the second permits the nearly Boltzmann electron response. A sufficiently large electron-to-ion temperature ratio creates this window in the ordinary electron–ion mass hierarchy.

If electron and ion temperatures are comparable, the acoustic phase speed is near the ion thermal population and ion Landau damping can be substantial. A real fluid sound frequency then does not prove that a well-defined weakly damped mode exists. Electron drift relative to ions can instead provide free energy and destabilize the branch, requiring the full shifted distributions. Distinguish a closure calculation, a weak-damping kinetic expansion, and a measured propagating signal before assigning the same label to their results.

vTi≪ωrk≪vTe,vTs=kBTsms,ωr≪ωpev_{Ti}\ll\frac{\omega_r}{k}\ll v_{Te},\qquad v_{Ts}=\sqrt{\frac{k_BT_s}{m_s}},\qquad \omega_r\ll\omega_{pe}
These are sufficient scale-separation conditions for the elementary ion-acoustic interpretation, with Te ≫ Ti usually needed to make ion damping small. They are orderings, not exact universal instability thresholds. The thermal speeds use the same convention as the dielectric lesson.↔ Scroll wide equations horizontally.

A finite-shielding correction

A cold-ion, Boltzmann-electron model has kBTe/mi = 9.0 × 10⁸ m²/s². At k = 20 m⁻¹ the electron Debye length is 0.0375 m. Find the phase speed and angular frequency, keeping the Debye correction.

  1. The long-wavelength sound speed is cs=9.0×108=3.0×104 m/sc_s=\sqrt{9.0\times10^8}=3.0\times10^4\,\mathrm{m/s}. The dimensionless shielding parameter is kλDe = 0.75.

  2. The denominator for phase speed is 1+0.752=1.5625=1.25\sqrt{1+0.75^2}=\sqrt{1.5625}=1.25. Thus phase speed is 30,000/1.25 = 24,000 m/s.

  3. Multiply by wavenumber to obtain ωr=20(24000)=4.8×105 rad/s\omega_r=20(24000)=4.8\times10^5\,\mathrm{rad/s}. The quasineutral long-wavelength estimate would instead give 6.0 × 10⁵ rad/s.

The result

The model predicts phase speed 24 km/s and angular frequency 4.8 × 10⁵ rad/s. Since kλDe is not very small, the finite-Debye correction matters.

Explore the model

Interactive lab

Interactive lab · waves

A closure relates perturbations with a phase

Reset and set BOTH Phase difference to 0 degrees and Amplitude ratio to 1; predict aligned peaks. Then Reset and set BOTH Phase difference to 90 degrees and Amplitude ratio to 1; predict a quarter-cycle offset. The fixed-frequency co-propagating waves illustrate relative phase only. They do not calculate ion-acoustic dispersion, electron shielding, a species density ratio, or Landau damping.

A closure relates perturbations with a phaseTwo waves share the same frequency and wave number. Playback advances t in sin(2x − t) + a sin(2x − t + phase). Current readings: Sum at x = 0: 0.866; Resultant amplitude: 1.732; Phase difference: 60°. All parameters are adjustable with the labeled controls below.−3.14−1.5701.573.14−3.4−1.701.73.4x · radiansAmplitude
SumFirst componentSecond component
Sum at x = 0
0.866
Resultant amplitude
1.732
Phase difference
60°
t = 0 s

Two waves share the same frequency and wave number. Playback advances t in sin(2x − t) + a sin(2x − t + phase). 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 real leading density response and a phase-shifted kinetic correction are different pieces of a response model.

Check your understanding

Why can a fluid calculation predict an acoustic frequency while a plasma with Te ≈ Ti has no weakly damped ion-acoustic signal?

Try it yourself

In the long-wavelength cold-ion limit, quadruple electron temperature at fixed density, ion mass, and wavenumber. Initially kλDe = 0.02. Find the changes in cs, λDe, and leading acoustic frequency, and check the wavelength ordering.

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

Give me a hint

Both sound speed and Debye length scale with the square root of electron temperature under these fixed conditions.

Compare with the solution

Quadrupling Te doubles cs and doubles λDe. The new kλDe is 0.04, still much smaller than one. Hence the leading acoustic frequency kcs doubles. The finite-Debye formula gives a ratio 21+0.022/1+0.0422\sqrt{1+0.02^2}/\sqrt{1+0.04^2}, just below two. This comparison holds within the stipulated cold-ion response and does not by itself determine a damping rate.

Take it with you

  • Boltzmann electron density rises with positive electrostatic potential.
  • Ion inertia and electron pressure produce an acoustic branch with a finite-Debye correction.
  • A real fluid frequency does not guarantee weak damping when resonant ions are abundant.

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.