research36 min studyLesson 305 in Physics

Two-stream and bump-on-tail instabilities

Derive an exact cold two-beam instability band and distinguish reactive coupling from resonant growth driven by a positive velocity-space slope.

On this page

What you’ll understand

  • Calculate the unstable branch of a symmetric cold two-stream model.
  • Distinguish reactive beam coupling from weak resonant bump-on-tail growth.
  • Explain why a local positive slope alone is not a complete stability test.

A plasma away from a single Maxwellian can store free energy in relative drifts or a secondary velocity-space bump. Small electrostatic disturbances may then grow by drawing energy from the particles. Two-stream and bump-on-tail instabilities describe related but distinct limits. Cold beams can couple reactively without a smooth resonant slope, whereas a weak dilute bump often modifies a Langmuir wave through resonant energy exchange. Stating the equilibrium distribution is more informative than assigning every growing electrostatic wave the same name.

Cold beam susceptibilities can be combined exactly

Consider two equal-density electron beams with velocities plus and minus U, their total density defining the electron plasma frequency. Immobile ions neutralize the total electron charge. Each beam contributes half the total squared plasma frequency in its Doppler-shifted cold susceptibility. This factor of one half matters: assigning the full density to both beams double counts the electrons and shifts the instability threshold. The electrostatic model neglects transverse electromagnetic instabilities of the same drifting system.

Combining the two rational terms gives a quadratic equation in frequency squared. Introduce a dimensionless squared wavenumber parameter so that the lower branch can be inspected directly. For a nonzero parameter below one, that branch is negative, giving one exponentially growing root and one decaying root. At exactly zero wavenumber the original rational equation has degeneracies, so the limiting low-frequency branch should not be mistaken for an additional finite-frequency physical mode at k = 0.

ϵ=1−ωp2/2(ω−kU)2−ωp2/2(ω+kU)2=0,Q=k2U2ωp2\epsilon=1-\frac{\omega_p^2/2}{(\omega-kU)^2}-\frac{\omega_p^2/2}{(\omega+kU)^2}=0,\qquad Q=\frac{k^2U^2}{\omega_p^2}
ωp uses total beam density, U is in m/s, and k in m⁻¹. Both beams are cold delta distributions and have equal densities. The equilibrium is homogeneous with a stationary neutralizing ion background.↔ Scroll wide equations horizontally.

An imaginary root is a reactive instability

Solving the quadratic yields two values of normalized frequency squared. For the unstable branch choose the square root with positive imaginary frequency under exp(ikx − iωt), which gives temporal growth. This is an exact result within the cold two-beam model, not a weak-damping perturbation of a smooth Maxwellian. Finite beam spread changes the dispersion and may suppress the instability, so the cold threshold cannot be used unchanged for an arbitrary warm distribution.

A warm distribution can have two maxima and a minimum between them without being unstable at every wavenumber. The full dielectric response or a suitable Penrose/Nyquist test determines whether a growing root exists. A necessary shape feature and a sufficient spectral condition are different statements. This distinction is especially useful when a velocity histogram has noisy positive slopes: not every small local rise in sampled data proves an exponentially growing collective mode.

ω±2ωp2=Q+12±121+8Q,0<Q<1⟹ω−2<0\frac{\omega_\pm^2}{\omega_p^2}=Q+\frac12\pm\frac12\sqrt{1+8Q},\qquad 0<Q<1\Longrightarrow\omega_-^2<0
The ± labels denote branches of frequency squared, not propagation direction. A negative lower value gives growth rate γ = √(−ω₋²) in s⁻¹. The instability band assumes the ideal cold symmetric electrostatic model.↔ Scroll wide equations horizontally.

A dilute bump can reverse resonant energy transfer

For a weak bump on a broad background, a Langmuir branch can remain close to its original real frequency while the distribution slope at phase velocity changes sign. On a positive-energy branch with positive real dielectric frequency derivative, a positive total resonant slope gives growth instead of damping. Use the total distribution: a rising beam contribution can be outweighed by a falling bulk contribution. Other species and collisions can also supply competing damping.

Linear growth cannot continue indefinitely without changing its own assumptions. Growing fields trap resonant particles, modify the velocity distribution, and can couple to other modes. A plateau can reduce resonant drive, but the saturation amplitude and late-time state depend on coherence, spectrum, boundaries, and collisions. A linear growth rate is therefore a prediction for the early interval of a specified equilibrium, not a universal estimate of the final turbulent transport or a proof of complete relaxation to a Maxwellian.

γ≃πωp2k2∂ωϵr g0′(ωr/k)(weak resonant correction)\gamma\simeq\frac{\pi\omega_p^2}{k^2\partial_\omega\epsilon_r}\,g_0^{\prime}(\omega_r/k)\quad\text{(weak resonant correction)}
This expression assumes a simple weak-growth root, k > 0, the causal exp(ikx − iωt) convention, and one normalized electron distribution with fixed ions. Inferring growth from a positive slope additionally requires ∂ωεr > 0 for the branch.↔ Scroll wide equations horizontally.

An unstable cold branch with a simple radical

In the equal cold two-beam model, choose Q = 3/8 and ωp = 8 × 10⁵ rad/s. Find the lower squared frequency and its temporal growth rate.

  1. The discriminant factor is 1+8Q=1+3=41+8Q=1+3=4, so its square root is exactly 2.

  2. The lower normalized squared frequency is 3/8+1/2−1=−1/83/8+1/2-1=-1/8. Therefore ω−2=−ωp2/8\omega_-^2=-\omega_p^2/8, which is negative.

  3. Choose the growing root ω=+iωp/8\omega=+i\omega_p/\sqrt8. The growth rate is γ=(8×105)/8≃2.83×105 s−1\gamma=(8\times10^5)/\sqrt8\simeq2.83\times10^5\,\mathrm{s^{-1}} under the stated time convention.

The result

The lower branch has ω² = −ωp²/8 and grows with gamma approximately 2.83 × 10⁵ s⁻¹. This is not a weak-growth approximation.

Explore the model

Interactive lab

Interactive lab · function

A slope changes sign locally

Reset and set BOTH Coefficient a to 1 and Position x to −1; predict tangent slope −2. Then Reset and set BOTH a to 1 and x to 1; predict slope +2. The exact curve y = ax² illustrates a derivative sign only. It is not a normalized nonnegative velocity distribution on the full line and does not calculate a Penrose criterion, dielectric root, or growth rate.

A slope changes sign locallyThis lab uses only y = ax². Playback traces the curve within the fixed integration interval from 0 to x; it does not estimate numerical error or change the chosen point. Current readings: y at chosen x: 0.608; Tangent slope · 2ax: 1.35; Signed area · ax³/3: 0.182. All parameters are adjustable with the labeled controls below.−2.4−1.201.22.4−9−4.504.59xy(x, ax²)
y = ax²TangentSigned integral
y at chosen x
0.608
Tangent slope · 2ax
1.35
Signed area · ax³/3
0.182
t = 0 s

This lab uses only y = ax². Playback traces the curve within the fixed integration interval from 0 to x; it does not estimate numerical error or change the chosen point. 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 local slope is one ingredient of resonant exchange; collective stability also requires the full response.

Check your understanding

In the symmetric cold two-beam model using total-density ωp, which nonzero Q range gives a negative lower ω² branch?

Try it yourself

For the same cold model set Q = 3. Find both normalized squared frequencies and decide whether it has temporal growth at this wavenumber.

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

Give me a hint

Here the square-root factor is √25. Inspect the sign of each squared-frequency branch.

Compare with the solution

The factor is 1+8(3)=5\sqrt{1+8(3)}=5. Thus ω+2/ωp2=3+1/2+5/2=6\omega_+^2/\omega_p^2=3+1/2+5/2=6 and ω−2/ωp2=3+1/2−5/2=1\omega_-^2/\omega_p^2=3+1/2-5/2=1. Both are positive, giving real roots ω=±6ωp\omega=\pm\sqrt6\omega_p and ω=±ωp\omega=\pm\omega_p. There is no exponential growth at this wavenumber in the ideal cold model, although the same equilibrium is unstable at other wavenumbers within 0 < Q < 1.

Take it with you

  • A symmetric cold beam instability follows from an exact coupled susceptibility calculation.
  • The total resonant slope and wave-energy sign determine weak bump-driven growth.
  • Linear instability predicts an early growth stage, not a universal nonlinear saturation state.

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.