Ion-Ion Acoustic Instability (IIAI)

• Ion-Ion Acoustic Instability (IIAI) is a beam-driven, kinetic electrostatic instability arising from relative drifts between ion populations in plasma.
• The instability requires hot electrons relative to ions, with growth rates sensitive to beam density, drift speed, and temperature ratios.
• Observations from the Parker Solar Probe and FPC diagnostics link IIAI to narrowband ion-acoustic bursts and energy transfer from ion beams to cores.

Searching arXiv for recent and foundational papers directly relevant to ion-ion acoustic instability. Ion-Ion Acoustic Instability (IIAI) is an electrostatic instability driven by relative drift between two ion populations, typically a proton core and a drifting proton beam, in a plasma where the electron response permits an ion-acoustic-like mode to grow rather than be suppressed by electron Landau damping. In the near-Sun solar-wind literature, IIAI is discussed as a beam-driven ion-acoustic variant operating in a proton-electron plasma and motivated in part by Parker Solar Probe observations of narrowband electrostatic emissions interpreted as ion-acoustic waves produced most likely by the ion-ion beam instability (Mozer et al., 2021). More generally, the term is used most precisely for a kinetic instability in which the ion-acoustic phase speed lies within the positive-slope region of a drifting ion distribution, so that beam ions transfer energy to the wave, while the core and electrons receive or dissipate energy according to the local resonance structure (Afify et al., 13 Jan 2026).

1. Definition and basic physical picture

The defining feature of IIAI is that the free energy source is the relative drift between two ion populations, not an electron current. In the solar-wind-motivated formulation, the plasma contains three Maxwellian species: a stationary proton core, a drifting proton beam, and background electrons, and the instability is treated as a kinetic, electrostatic ion-acoustic mode driven by drifting ion populations (Afify et al., 13 Jan 2026). In a closely related 1D electrostatic formulation, the system is described as a proton core, a proton beam, and electrons, with quasi-neutrality $n_e=n_c+n_b$ and zero current enforced by an electron drift chosen to balance the ion beam current (Afify et al., 2024).

The physical mechanism is resonant wave-particle interaction. Ion-acoustic waves satisfy approximately

$\omega_r \simeq k_\parallel c_s,$

with ion-acoustic speed

$$c_s = \sqrt{\frac{3k_B T_{\parallel i}+k_B T_{\parallel e}}{m_i}},$$

and resonant interaction occurs when

$v_\parallel = v_{\mathrm{res}} = \omega/k.$

The sign of the velocity-space gradient at the resonance determines the energy flow: if

$\frac{\partial f}{\partial v_\parallel} > 0$

at the resonant velocity, particles lose energy to the wave and drive instability; if

$\frac{\partial f}{\partial v_\parallel} < 0,$

the particles gain energy and the wave is Landau damped (Afify et al., 13 Jan 2026). In this sense, IIAI is a beam-driven ion-acoustic instability whose onset depends on where the ion-acoustic resonance falls relative to the velocity-space slope of the drifting beam.

A recurrent condition in the literature is that electrons must be sufficiently hot relative to ions for electron Landau damping to remain weak enough not to suppress growth. In the near-Sun observational case, $T_e/T_i \sim 5$ is described as favorable for ion-acoustic waves because ion-acoustic damping is weak when electrons are sufficiently hotter than ions (Mozer et al., 2021). In the PSP-motivated kinetic study of electrostatic bursts, the instability is likewise said to require electron-to-core and beam-to-core temperature ratios slightly different from reported values during electrostatic burst detection, and the summary states that instability requires roughly $T_e/T_c \gtrsim 8-10$ together with a narrow interval in beam drift and modest $T_b/T_c$ (Afify et al., 2024).

2. Linear kinetic formulation and instability threshold

In the fully kinetic solar-wind treatment, the underlying model is a 1D-1V Vlasov-Poisson system for species $s\in\{c,b,e\}$: $\omega_r \simeq k_\parallel c_s,$0 with Poisson’s equation

$\omega_r \simeq k_\parallel c_s,$1

and $\omega_r \simeq k_\parallel c_s,$2 (Afify et al., 13 Jan 2026). The corresponding linear electrostatic dispersion relation is written as

$\omega_r \simeq k_\parallel c_s,$3

where $\omega_r \simeq k_\parallel c_s,$4,

$\omega_r \simeq k_\parallel c_s,$5

This relation is used to map how the instability growth rate depends on electron temperature, beam density, and beam drift (Afify et al., 13 Jan 2026).

A second linear kinetic representation, written for drifting Maxwellian species, is

$\omega_r \simeq k_\parallel c_s,$6

with

$\omega_r \simeq k_\parallel c_s,$7

and

$\omega_r \simeq k_\parallel c_s,$8

In this formulation the control parameters are $\omega_r \simeq k_\parallel c_s,$9, the beam-core drift $$c_s = \sqrt{\frac{3k_B T_{\parallel i}+k_B T_{\parallel e}}{m_i}},$$0, and the temperature ratios $$c_s = \sqrt{\frac{3k_B T_{\parallel i}+k_B T_{\parallel e}}{m_i}},$$1 and $$c_s = \sqrt{\frac{3k_B T_{\parallel i}+k_B T_{\parallel e}}{m_i}},$$2 (Afify et al., 2024).

The threshold is not a single universal value but a boundary in parameter space. The 2025 study of nonthermal electrons defines the threshold in the standard linear-kinetic way by solving for $$c_s = \sqrt{\frac{3k_B T_{\parallel i}+k_B T_{\parallel e}}{m_i}},$$3, with instability for $$c_s = \sqrt{\frac{3k_B T_{\parallel i}+k_B T_{\parallel e}}{m_i}},$$4 and threshold at $$c_s = \sqrt{\frac{3k_B T_{\parallel i}+k_B T_{\parallel e}}{m_i}},$$5; the threshold is explored chiefly in the space of $$c_s = \sqrt{\frac{3k_B T_{\parallel i}+k_B T_{\parallel e}}{m_i}},$$6, $$c_s = \sqrt{\frac{3k_B T_{\parallel i}+k_B T_{\parallel e}}{m_i}},$$7, $$c_s = \sqrt{\frac{3k_B T_{\parallel i}+k_B T_{\parallel e}}{m_i}},$$8, and electron-distribution shape and temperature (&&&10&&&). This is consistent with the 2024 PSP-motivated analysis, which finds that the IIAI only exists in a restricted window in the plane of $$c_s = \sqrt{\frac{3k_B T_{\parallel i}+k_B T_{\parallel e}}{m_i}},$$9 versus $v_\parallel = v_{\mathrm{res}} = \omega/k.$0, with both minimum and maximum drift for fixed beam fraction (Afify et al., 2024).

Several threshold trends recur across the studies. Increasing $v_\parallel = v_{\mathrm{res}} = \omega/k.$1 or $v_\parallel = v_{\mathrm{res}} = \omega/k.$2 strengthens growth, while varying the drift speed changes growth nonmonotonically as the resonance moves relative to the beam gradient (Afify et al., 13 Jan 2026). Increasing $v_\parallel = v_{\mathrm{res}} = \omega/k.$3 stabilizes the system, and decreasing $v_\parallel = v_{\mathrm{res}} = \omega/k.$4 also stabilizes it (Afify et al., 2024). For the dilute PSP-like beam with $v_\parallel = v_{\mathrm{res}} = \omega/k.$5 and $v_\parallel = v_{\mathrm{res}} = \omega/k.$6, the 2024 study finds the electron-to-core temperature threshold is approximately

$v_\parallel = v_{\mathrm{res}} = \omega/k.$7

With the observed $v_\parallel = v_{\mathrm{res}} = \omega/k.$8, the nominal PSP parameters are therefore slightly below threshold in the simplest Maxwellian model (Afify et al., 2024).

3. Plasma conditions and near-Sun observational evidence

A major observational stimulus for IIAI studies is the Parker Solar Probe event on 2021-01-18 to 2021-01-19 near the spacecraft’s 20 solar radius perihelion, where FIELDS and SWEAP measured continuous narrowband electrostatic emissions for about 12 hours (Mozer et al., 2021). The emissions were observed at $v_\parallel = v_{\mathrm{res}} = \omega/k.$9–$\frac{\partial f}{\partial v_\parallel} > 0$0 Hz in the spacecraft frame, with a later interval around $\frac{\partial f}{\partial v_\parallel} > 0$1 Hz, and they were below the local ion plasma frequency and without a magnetic-field counterpart, supporting an electrostatic interpretation (Mozer et al., 2021).

The waves appeared as wave packets with shock-like envelopes repeating at about $\frac{\partial f}{\partial v_\parallel} > 0$2 Hz, and the repetition was phase correlated with a few-Hz electromagnetic fluctuation in $\frac{\partial f}{\partial v_\parallel} > 0$3 and $\frac{\partial f}{\partial v_\parallel} > 0$4 (Mozer et al., 2021). This produced the interpretive picture of ion-acoustic packets triggered in synchrony with a lower-frequency electromagnetic oscillation rather than random isolated bursts (Mozer et al., 2021).

The measured plasma state was favorable for ion-acoustic excitation. SWEAP core+beam Maxwellian fits gave a core density of $\frac{\partial f}{\partial v_\parallel} > 0$5, a beam density of $\frac{\partial f}{\partial v_\parallel} > 0$6, beam-core drift speed about $\frac{\partial f}{\partial v_\parallel} > 0$7, core perpendicular temperature $\frac{\partial f}{\partial v_\parallel} > 0$8, beam perpendicular temperature $\frac{\partial f}{\partial v_\parallel} > 0$9, core anisotropy $\frac{\partial f}{\partial v_\parallel} < 0,$0, beam anisotropy $\frac{\partial f}{\partial v_\parallel} < 0,$1, and electron temperature $\frac{\partial f}{\partial v_\parallel} < 0,$2, implying $\frac{\partial f}{\partial v_\parallel} < 0,$3 (Mozer et al., 2021). The background solar-wind speed was about $\frac{\partial f}{\partial v_\parallel} < 0,$4, described as slow wind, and the beam was anti-sunward and faster than the core (Mozer et al., 2021).

The observational identification as ion-acoustic waves rests on four explicit arguments: the emissions are electrostatic, their phase speed is of order ion thermal or ion-acoustic speed, they occur under $\frac{\partial f}{\partial v_\parallel} < 0,$5 with an ion beam present, and their properties match a marginal ion-ion acoustic instability regime (Mozer et al., 2021). The authors state that inspection of Gary and Omidi (1987) shows the plasma is at the marginal stability threshold of the ion-ion acoustic instability, although the paper does not derive a full growth-rate calculation from first principles (Mozer et al., 2021).

The same observational paper is careful about its scope. It states explicitly that it does not present a classic, fully developed IIAI calculation, but instead provides a strong observational case that the detected waves are ion-acoustic waves most likely triggered by an ion-beam or ion-ion acoustic instability operating near marginal stability in the near-Sun solar wind (Mozer et al., 2021). The contribution is therefore phenomenological and interpretive rather than a new dispersion calculation.

4. Nonlinear evolution, bursts, and phase-space structures

The 2024 PSP-motivated Vlasov study addresses whether IIAI can explain the high-frequency electrostatic bursts observed between about 15 and 25 solar radii (Afify et al., 2024). Its main claim is that a proton core plus drifting proton beam embedded in a sufficiently hot electron background can drive IIAI, and that the resulting nonlinear electrostatic structures resemble the PSP bursts in frequency, duration, and amplitude (Afify et al., 2024).

For the simulation cases with $\frac{\partial f}{\partial v_\parallel} < 0,$6, $\frac{\partial f}{\partial v_\parallel} < 0,$7, $\frac{\partial f}{\partial v_\parallel} < 0,$8 and $\frac{\partial f}{\partial v_\parallel} < 0,$9, and $T_e/T_i \sim 5$0 or $T_e/T_i \sim 5$1, the measured exponential growth rates are

$T_e/T_i \sim 5$2

slightly larger than the linear predictions but consistent with them (Afify et al., 2024). These simulations reproduce the theoretical trends that more beam density yields faster growth, a hotter beam yields slower growth, and the fastest growth occurs near the wavelength predicted by linear theory (Afify et al., 2024).

The nonlinear saturated state contains traveling vortices or islands in the beam distribution and strong phase-space trapping of beam protons (Afify et al., 2024). In the weakest case, density and temperature changes remain small, beam density varies by about $T_e/T_i \sim 5$3, and no substantial heating is observed (Afify et al., 2024). The normalized saturated electric-field amplitude is about

$T_e/T_i \sim 5$4

converted to roughly $T_e/T_i \sim 5$5 using PSP-like parameters, which is said to be compatible with the observed burst amplitudes within uncertainties (Afify et al., 2024).

The same study estimates a growth time

$T_e/T_i \sim 5$6

using the weakest-case growth rate and $T_e/T_i \sim 5$7, and states that this is in good agreement with the observed burst timescales (Afify et al., 2024). A plausible implication is that near-threshold IIAI can account for both the temporal scale and the burst-like morphology of the PSP event, although the paper also notes that the simplest Maxwellian model requires $T_e/T_i \sim 5$8 somewhat larger than the nominal observed value and that burst decay seen by PSP is not reproduced in the periodic 1D model (Afify et al., 2024).

The observational PSP analysis proposes a staged nonlinear evolution in which low-frequency electromagnetic fluctuations may locally destabilize the plasma, producing quasi-monochromatic ion-acoustic emissions, and subsequent nonlinear processes including dispersion, nonlinear steepening, and particle trapping may turn these into shock-like wave packets (Mozer et al., 2021). This combination of narrowband spectrum, repetition, and shock-like packet envelopes is explicitly claimed not to have been previously reported for ion-acoustic dynamics (Mozer et al., 2021).

5. Velocity-space energy transfer and diagnostic signatures

A central development in the recent IIAI literature is the use of field-particle correlation (FPC) to diagnose energy transfer in fully kinetic simulations (Afify et al., 13 Jan 2026). Starting from the Vlasov equation and multiplying by kinetic energy $T_e/T_i \sim 5$9, the phase-space energy transfer rate is written as

$T_e/T_c \gtrsim 8-10$0

with

$T_e/T_c \gtrsim 8-10$1

Integration over space, velocity, and time gives the cumulative nonlinear energy transfer

$T_e/T_c \gtrsim 8-10$2

while direct energy accounting uses

$T_e/T_c \gtrsim 8-10$3

The paper shows that $T_e/T_c \gtrsim 8-10$4 and $T_e/T_c \gtrsim 8-10$5 agree, confirming that the FPC-based expression measures the same physical energy exchange as direct energy accounting (Afify et al., 13 Jan 2026).

The actual single-point FPC diagnostic is

$T_e/T_c \gtrsim 8-10$6

with $T_e/T_c \gtrsim 8-10$7 (Afify et al., 13 Jan 2026). Time averaging over an interval longer than the wave period removes oscillatory reversible exchange and leaves secular net transfer, which the paper describes as separating energy “sloshing” from true dissipation-like transfer (Afify et al., 13 Jan 2026). In the cases studied, $T_e/T_c \gtrsim 8-10$8 is found long enough to reveal secular transfer (Afify et al., 13 Jan 2026).

The identified IIAI signatures are specific. For the proton beam, the dominant signature is negative FPC at the resonant velocity, indicating that beam particles lose energy to the wave (Afify et al., 13 Jan 2026). For the proton core, the signature is positive at the same resonant velocity, showing that the core gains energy (Afify et al., 13 Jan 2026). Electrons generally show much weaker net energization; in the fiducial case, the electron FPC is mostly oscillatory, with little or no clear secular transfer (Afify et al., 13 Jan 2026). The paper stresses that the core is the main recipient of the beam’s lost energy, whereas the electrons receive only a minor share (Afify et al., 13 Jan 2026).

These results supply a velocity-space interpretation of IIAI saturation. The beam is the primary energy source for the instability, the core is the main beneficiary, and electrons are involved in the field-particle interaction but are not the dominant energy sink under the simulated solar-wind-like conditions (Afify et al., 13 Jan 2026). This provides a kinetic complement to the more phenomenological observational picture of triggered ion-acoustic packets in the young solar wind.

6. Electron-distribution effects and near-threshold sensitivity

Because the instability is often close to marginality in PSP-like conditions, the shape of the electron distribution is an important question. The 2025 study of nonthermal electrons compares Maxwellian, kappa, and core-strahl electron models in a PSP-relevant IIAI regime (Afify et al., 22 Sep 2025).

For the kappa study, the adopted reference parameters are

$T_e/T_c \gtrsim 8-10$9

The 1D standard kappa distribution is

$T_b/T_c$0

with

$T_b/T_c$1

The key result is that decreasing $T_b/T_c$2 reduces the IIAI growth rate, meaning that kappa electrons stabilize the instability (Afify et al., 22 Sep 2025). The explanation given is that lower $T_b/T_c$3 increases the electron phase-space density near the resonance and enhances electron Landau damping (Afify et al., 22 Sep 2025).

The same paper validates this with 1D1V Vlasov-Poisson simulations using $T_b/T_c$4, $T_b/T_c$5, periodic boundaries in $T_b/T_c$6, zero-flux boundaries in velocity space, and the mode $T_b/T_c$7 (Afify et al., 22 Sep 2025). Example growth-rate comparisons are: $T_b/T_c$8, $T_b/T_c$9 theory and $s\in\{c,b,e\}$0 simulation; $s\in\{c,b,e\}$1, $s\in\{c,b,e\}$2 theory and $s\in\{c,b,e\}$3 simulation; $s\in\{c,b,e\}$4, $s\in\{c,b,e\}$5 theory and $s\in\{c,b,e\}$6 simulation (Afify et al., 22 Sep 2025).

For core-strahl electrons, modeled as two Maxwellians with fixed total electron density and zero current, the study finds that a hot strahl tends to destabilize the IIAI relative to a single Maxwellian electron population: higher $s\in\{c,b,e\}$7 increases growth, higher $s\in\{c,b,e\}$8 increases growth, and the unstable range in $s\in\{c,b,e\}$9 broadens with increasing strahl density (Afify et al., 22 Sep 2025). However, the effect is still described as modest (Afify et al., 22 Sep 2025).

A major conclusion is the confirmation of the Jones et al. (1975) effective temperature

$\omega_r \simeq k_\parallel c_s,$00

for a core-strahl electron population (Afify et al., 22 Sep 2025). The authors report that replacing the core-strahl distribution with a single Maxwellian at $\omega_r \simeq k_\parallel c_s,$01 gives theoretical dispersion relations that are identical, simulation results that are very close, and wave-particle interaction patterns that are the same (Afify et al., 22 Sep 2025). The representative values in their table are $\omega_r \simeq k_\parallel c_s,$02 to $\omega_r \simeq k_\parallel c_s,$03 (Afify et al., 22 Sep 2025).

The overall conclusion of that study is explicit: kappa electrons stabilize the IIAI somewhat, core-strahl electrons destabilize it somewhat, but neither effect is large enough to dramatically lower the threshold or fully explain the PSP-observed IIAI event (Afify et al., 22 Sep 2025). This suggests that electron non-Maxwellianity should be included in stability assessments, but should not be assumed to be the dominant explanation for threshold crossing.

7. Distinctions, related instabilities, and limitations of the term

The term IIAI is sometimes used broadly for beam- or flow-driven ion-acoustic destabilization, but the mechanisms grouped under that label are not all equivalent. The most direct form is the kinetic beam-driven ion-acoustic mode described above, in which a drifting ion beam resonates with an electrostatic wave and drives growth (Afify et al., 13 Jan 2026, Afify et al., 2024). The near-Sun observational interpretation likewise uses the phrase “ion beam instability, also known as the ion-ion acoustic instability,” and concludes that the observed waves were “produced most likely by the ion-ion beam instability” (Mozer et al., 2021).

By contrast, some related papers describe mechanisms that are ion-acoustic in outcome but not classic IIAI. The electron-hole oscillatory velocity instability is explicitly distinguished from classic ion-ion acoustic instability: it is a hole instability mediated by passing ions, can persist when ion-ion type instability and Buneman instability are ruled out, and emits ion-acoustic-like density waves as a nonlinear consequence (Zhou et al., 2017). Its instability criterion concerns the electron-hole speed in the ion frame, with one example giving $\omega_r \simeq k_\parallel c_s,$04, and its emitted perturbations propagate in the ion frame with the ion sound speed, mainly in the opposite direction to the electron-hole velocity (Zhou et al., 2017). This is therefore ion-driven and ion-acoustic in radiation, but not an ion-ion acoustic instability in the standard sense (Zhou et al., 2017).

A different non-kinetic mechanism appears in finite-length plasma systems. “Ion sound instability driven by ion beam” studies a finite-length ion-flow-driven ion-sound instability in which stationary ion flow creates positive- and negative-energy modes via Doppler shift, and instability develops through coupling of these modes mediated by boundary reflection rather than by kinetic beam resonance (Koshkarov et al., 2014). The quasineutral limit is stable, and finite Debye-length dispersion is essential for instability (Koshkarov et al., 2014). In the language of IIAI, this is a beam- or flow-driven ion-acoustic instability in a bounded plasma, but its distinctive mechanism is reactive and boundary mediated rather than resonant (Koshkarov et al., 2014).

Further caution is needed because some papers focus on the ordinary electron-ion ion-acoustic instability rather than IIAI proper. The reconnection study reports ion-acoustic instability driven by electron-ion drift/current in the diffusion region and explicitly does not discuss a separate, distinct ion-ion acoustic branch (Li et al., 13 May 2025). The helicon-discharge studies likewise concern electron-driven high-frequency ion-acoustic instability in cylindrical, radially inhomogeneous plasmas, not ion-ion coupling (Mikhailenko et al., 2022, Mikhailenko et al., 2022). These works are relevant for ion-acoustic turbulence and thresholds but should not be conflated with the proton core-beam IIAI of solar-wind studies.

A final limitation concerns present-day in situ diagnosis. The FPC study concludes that the identified IIAI signatures are physically well suited to single-spacecraft diagnosis, but the timescale over which they develop is too fast for current missions: Parker Solar Probe’s fastest SPAN-i cadence is about $\omega_r \simeq k_\parallel c_s,$05 s, whereas the IIAI growth times inferred there are of order tens of milliseconds (Afify et al., 13 Jan 2026). This means that direct resolution of the full temporal development of proton and electron velocity-space signatures is generally beyond current inner-heliosphere particle sampling capabilities (Afify et al., 13 Jan 2026).

Taken together, the recent literature supports a precise use of IIAI as a beam-driven, kinetic, electrostatic ion-acoustic instability in a multi-ion proton-electron plasma, with resonance-controlled onset, strong sensitivity to electron damping, and nonlinear evolution that can produce narrowband bursts, beam trapping, ion holes, and phase-space energy transfer from the beam to the core (Afify et al., 2024, Afify et al., 13 Jan 2026). The PSP observations strengthen the case that such conditions are realized near the Sun and may occur close to marginal stability, where small changes in local plasma parameters or low-frequency triggering can repeatedly excite oblique ion-acoustic wave packets with shock-like envelopes (Mozer et al., 2021).