pith. sign in

arxiv: 2604.11262 · v1 · submitted 2026-04-13 · 🌌 astro-ph.EP

Streaming instabilities in weakly ionized protoplanetary discs: the Ambipolar Streaming Instability (AmSI)

Arnaud Pierens , Min-Kai Lin This is my paper

Pith reviewed 2026-05-10 16:26 UTC · model grok-4.3

classification 🌌 astro-ph.EP
keywords ambipolar diffusionresonant drag instabilityprotoplanetary discsstreaming instabilityAlfven wavesdust dynamicsplanet formationMHD
0
0 comments X

The pith

Ambipolar diffusion triggers a resonant drag instability that clumps dust even in low-density regions of protoplanetary discs.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

In weakly ionized protoplanetary discs, ambipolar diffusion modifies the frequency of Alfvén waves. This creates a wide resonance with the drift between gas and dust, destabilizing the waves and driving a strong resonant drag instability. The instability grows at significant rates even when dust is scarce and particles are tightly coupled to the gas. Dust feedback also damps the oblique modes of the magnetorotational instability. If correct, the process supplies a route from dust coagulation to planetesimal formation across a wider range of disc conditions.

Core claim

The paper claims that ambipolar diffusion in a dusty, magnetized disc leads to the Ambipolar Streaming Instability (AmSI), a resonant drag instability in which an Alfvén wave is destabilized by the relative gas-dust drift. The main effect of ambipolar diffusion is to alter the Alfvén wave frequency, producing a large resonance width. This yields significant growth rates for dust-to-gas ratios as low as 0.01 and for particles with stopping times much less than the orbital period.

What carries the argument

The resonant drag instability (RDI) of an Alfvén wave whose frequency is shifted by ambipolar diffusion, creating a broad resonance with the gas-dust drift velocity.

If this is right

  • The instability maintains high growth rates in dust-poor discs.
  • It operates efficiently for tightly coupled dust particles.
  • Dust feedback suppresses the oblique MRI modes.
  • The mechanism may connect grain coagulation directly to planetesimal formation.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • This could allow planetesimal formation in disc regions previously considered too dust-poor for other streaming instabilities.
  • The same ambipolar modification of wave frequencies might influence the launching of MHD disc winds.
  • Nonlinear evolution of the instability could be tested with global disc simulations that include both ambipolar diffusion and dust dynamics.

Load-bearing premise

The analysis assumes a specific weakly ionized disc model that combines ohmic resistivity and ambipolar diffusion, with linear stability analysis and dust feedback on MRI modes capturing the dominant dynamics.

What would settle it

A direct numerical simulation or observation showing either rapid dust clumping at ambipolar-diffusion-dominated radii with low dust-to-gas ratio, or the complete absence of such clumping despite the predicted growth rates.

Figures

Figures reproduced from arXiv: 2604.11262 by Arnaud Pierens, Min-Kai Lin.

Figure 1
Figure 1. Figure 1: Growth rates of unstable modes for different amplitudes of the azimuthal magnetic field, for a dust-free disc with ϵ = 0 (top panel), a dust-poor disc with ϵ = 0.2 (middle panel), a dust-rich disc with ϵ = 3. Here, the Stokes number is fixed to to St = 0.1 and the ambipolar Elsasser number to ΛAD = 100. The solid line corresponds to the resonant drag instability (RDI) condition for the streaming instabilit… view at source ↗
Figure 2
Figure 2. Figure 2: Growth rates of unstable modes for different amplitudes of the azimuthal magnetic field. Here, the dust-to-gas ratio is fixed to ϵ = 3, the Stokes number to St = 0.1 and the ambipolar Elsasser number to ΛAD = 10. 10 3 10 0 10 3 Kx 10 3 10 1 10 1 10 3 10 5 Kz St = 10 4 10 3 10 0 10 3 Kx St = 10 3 10 3 10 0 10 3 Kx St = 10 2 10 3 10 0 10 3 Kx St = 10 1 10 5 10 4 10 3 10 2 10 1 10 0 / = 0.01, By/Bz = 1, z = 1… view at source ↗
Figure 3
Figure 3. Figure 3: Growth rates of unstable modes for different values of the Stokes number, for ϵ = 0.01 (top) and ϵ = 0.2 (bottom). Here, the ambipolar Elsasser number to ΛAD = 10 and the azimuthal field to By/Bz = −1. 10 ≲ K˜ x, K˜ z ≲ 103 , however, we find that ambipolar diffusion, on the contrary, tends to increase the growth rates of the unstable modes located slightly away from the straight line corresponding to the … view at source ↗
Figure 6
Figure 6. Figure 6: For ϵ = 0.01 and St = 0.001, ratio between oscillation frequen￾cies ℑ(σ) found by solving the stability problem (Eq. 23) and frequen￾cies expected from the resonance condition kx(wx − vx). a higher ambipolar diffusivity. The smaller Elsasser number also causes the classic SI modes to appear in the blue-green rectangu￾lar region corresponding to 1 ≲ K˜ x ≲ 102 and 102 ≲ K˜ z ≲ 105 , with growth rates that i… view at source ↗
Figure 5
Figure 5. Figure 5: Velocity and magnetic field components along nˆ = kˆ × bˆ as a function of vertical wavenumber K˜ z , for ϵ = 3, St = 0.1, ΛAD = 100 (upper panel) and ϵ = 0.2, St = 0.001, ΛAD = 10 (lower panel). Here, the azimuthal field is fixed to By/Bz = −1 and the radial wavenumber to K˜ x = 103 . 10 0 10 1 10 2 10 3 10 4 Kx 10 0 10 1 10 2 10 3 10 4 10 5 Kz = 0.01, St=0.001, AD = 10, By/BZ = -1.0 0.95 1.00 1.05 ( ) kx… view at source ↗
Figure 7
Figure 7. Figure 7: Nonlinear evolution of the AmSI with St = 0.01 (blue) and St = 0.1 (orange). Left: the maximum gas velocity perturbation. Right: the maximum dust-to-gas ratio [PITH_FULL_IMAGE:figures/full_fig_p010_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Snapshots of the dust-to-gas ratios at the end of the AmSI simulations shown in [PITH_FULL_IMAGE:figures/full_fig_p010_8.png] view at source ↗
read the original abstract

The regions of protoplanetary discs where planets can form are believed to be weakly ionised, suggesting thereby that non-ideal magneto-hydrodynamics (MHD) effects play an important role in the disc dynamics and in the planet formation process. In particular, the combined effect of ohmic resistivity and ambipolar diffusion can be responsible for launching MHD-driven disc winds. In this context, we focus on the effect of ambipolar diffusion (AD) and examine the stability of a dusty, magnetized disc by employing both linear stability analyses and numerical simulations. We show that dust feedback tends to stabilize the MRI oblique modes involved in the ambipolar-shear instability. We also find that ambipolar diffusion leads to the onset of a strong resonant drag instability (RDI), in which an Alfv\'en wave is destabilized by the relative drift between the gas and dust components. The main impact of AD is to modify the Alfv\'en wave frequency, resulting in a large resonance width. The instability is found to have significant growth rates even in dust-poor discs and for tightly coupled particles, which may help to bridge the gap between growth of dust grains through coagulation and planetesimal formation.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 3 minor

Summary. The manuscript investigates the stability of dusty, weakly ionized protoplanetary discs subject to non-ideal MHD effects, focusing on the combined influence of ohmic resistivity and ambipolar diffusion (AD). Employing linear stability analysis and numerical simulations, the authors show that dust feedback stabilizes the oblique modes of the magnetorotational instability (MRI) associated with the ambipolar-shear instability. They further demonstrate that AD induces a strong resonant drag instability (RDI), termed the Ambipolar Streaming Instability (AmSI), in which an Alfvén wave is destabilized by the relative drift between gas and dust. The primary effect of AD is to modify the Alfvén wave frequency, producing a large resonance width. The instability exhibits significant growth rates even at low dust-to-gas ratios and for tightly coupled particles (small Stokes numbers).

Significance. If the central results hold, this work identifies a potentially important mechanism for dust clumping and planetesimal formation in regions of protoplanetary discs that are dust-poor and where particles remain tightly coupled to the gas, thereby helping to bridge the gap between grain coagulation and gravitational collapse. The combination of linear analysis (including explicit treatment of dust feedback on MRI modes) and numerical simulations provides a solid foundation for the claims. The identification of a broad resonance arising from AD-modified wave frequencies is a clear strength, as is the demonstration that the instability operates across a wider parameter space than classical streaming instabilities.

minor comments (3)
  1. The abstract states that the instability has 'significant growth rates' but provides no quantitative values or explored parameter ranges (e.g., specific dust-to-gas ratios or Stokes numbers); adding a brief quantification or reference to a figure/table would improve clarity.
  2. In the linear stability analysis, the dispersion relation for the modified Alfvén wave under AD should be presented with all terms explicitly defined and numbered for easy reference when discussing the resonance width.
  3. The numerical simulation section would benefit from a short statement on the grid resolution, box size, and boundary conditions used to confirm that the reported growth rates are converged and not affected by numerical diffusion.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their positive assessment of our manuscript, including the recognition of the Ambipolar Streaming Instability (AmSI) as a potentially important mechanism for dust clumping in dust-poor, tightly coupled regimes. We appreciate the recommendation for minor revision and the acknowledgment of the strengths in our linear analysis and numerical simulations.

Circularity Check

0 steps flagged

No significant circularity

full rationale

The paper derives the Ambipolar Streaming Instability (AmSI) and its growth rates directly from linear stability analysis of the governing MHD equations that incorporate ambipolar diffusion, ohmic resistivity, and dust feedback. The resonance width and destabilization of Alfvén waves follow from solving the dispersion relation without parameter fitting to the target instability or load-bearing self-citations that presuppose the result. Numerical simulations provide independent verification rather than circular confirmation. The approach is self-contained against the stated model assumptions.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 1 invented entities

Based solely on the abstract; the central claim rests on standard assumptions of weakly ionized disc MHD and linear perturbation theory rather than new free parameters or invented entities.

free parameters (2)
  • ambipolar diffusion coefficient
    Controls the modification of Alfvén wave frequency and resonance width; value not specified in abstract.
  • dust-to-gas ratio
    Instability growth rates depend on this parameter; abstract claims significance even at low values but does not report specific fitted values.
axioms (2)
  • domain assumption Protoplanetary disc regions where planets form are weakly ionized, so non-ideal MHD effects (ohmic resistivity and ambipolar diffusion) dominate.
    Stated as the physical context motivating the study.
  • domain assumption Linear stability analysis of the coupled gas-dust-MHD equations accurately identifies the dominant unstable modes.
    Basis for the reported growth rates and resonance properties.
invented entities (1)
  • Ambipolar Streaming Instability (AmSI) no independent evidence
    purpose: Name for the resonant drag instability arising from ambipolar diffusion in dusty discs.
    Newly introduced term in the abstract; no independent observational or experimental evidence provided.

pith-pipeline@v0.9.0 · 5513 in / 1641 out tokens · 85475 ms · 2026-05-10T16:26:38.056366+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

2 extracted references · 2 canonical work pages

  1. [1]

    P., Simon, J

    Abod, C. P., Simon, J. B., Li, R., & et al. 2019, ApJ, 883, 192 Bai, X.-N. 2011, ApJ, 739, 50 Bai, X.-N. 2017, ApJ, 845, 75 Bai, X.-N. & Stone, J. M. 2013, ApJ, 769, 76 Balbus, S. A. & Hawley, J. F. 1991, ApJ, 376, 214 Balsara, D. S., Tilley, D. A., Rettig, T., & et al. 2009, MNRAS, 397, 24 Béthune, W., Lesur, G., & Ferreira, J. 2017, A&A, 600, A75 Blum, ...

    work page 2019

  2. [2]

    Right: St=0.1

    Left: St=0.01. Right: St=0.1. Laibe, G. & Price, D. J. 2014, MNRAS, 440, 2136 Latter, H. N. & Papaloizou, J. 2017, MNRAS, 472, 1432 Lesur, G. 2021, Journal of Plasma Physics, 87, 205870101 Lesur, G. R. J., Baghdadi, S., Wafflard-Fernandez, G., et al. 2023, A&A, 677, A9 Lim, J., Baronett, S. A., Simon, J. B., et al. 2025, ApJ, 993, 12 Lin, M.-K. & Hsu, C.-...

    work page 2014