Diffusive Instabilities in Dusty Disks: Linear Growth and Nonlinear Breakdown

The gist The basic criterion for diffusion-slope-driven instability remains approximately βdiff ≲ −2 for small dust stopping times, with gas feedback providing only modest quantitative changes for parameters motivated by streaming-instability turbulence

Model Formulation

The paper revisits the diffusive instability in dusty disks by including an incompressible, viscous gas that responds azimuthally and couples to the dust through drag. The model starts with an advection-diffusion equation for a dust fluid where the diffusion flux F is interpreted as arising from correlations between density and velocity fluctuations, defined as vdiff ≡ −F/ρd. The pressurized dust momentum equation is derived, which can be written in terms of the effective dust velocity v = v + vdiff. The key relationship linking effective dust pressure and diffusion is that c2d = Dts (9) and vdiff = -Dρd∇ρd. Further allowing for a power law density dependence in dust diffusion, βD ≡ ∂(Dρd)/∂ρd = D(1 + βdiff) and dust viscosity, βν ≡ ∂(νρd)/∂ρd = ν(1 + βvisc). The final momentum equation is expressed in terms of the dimensionless diffusion slope as ∂v/∂t + v∇v = 1/ρd∇ · (ρdvdiffvdiff) + u − vts − 1/ρd Dts (2 + βdiff)∇ρd + 1/ρd∇ · T.

Linear Analysis and Dispersion Relation

The instability is a constitutive instability driven by the density-dependence of the transport coefficients themselves, rather than background gradients. The system is linearized around a constant background state, leading to a fourth-order dispersion relation n4 + a3n3 + a2n2 + a1n + a0 = 0. A sufficient criterion for an unstable static mode is that the coefficient of the leading term, a0, must be less than zero. This immediately implies that βdiff < −2 and βvisc < 1 are approximate criteria for diffusive instability driven by the diffusion slope and viscosity slope respectively. The growth rates for the fastest growing mode are analyzed using dimensionless quantities like τs ≡ tsomega, σ ≡ n/omega, K ≡ kH, δ0 ≡ D0/(csH), α0 ≡ ν0/(csH), and αg ≡ νg/(csH). The growth rate γ ≈ -(2 + βdiff) + 4Scg/3(1 + ε)ξ for σ << 1, confirming that βdiff ≲ −2 leads to instability.

Nonlinear Evolution and Saturation Issues

The nonlinear evolution of the model is investigated numerically, and the authors find that for power-law closures D ∝ Σ βdiff with βdiff < 0, the nonlinear evolution does not saturate. Instead, steepening gradients amplify the nonlinear dust-pressure term and drive finite-time collapse into increasingly sharp spikes. To address this, a simple piecewise closure is tested where the negative diffusion slope operates only over a finite density interval to eliminate blowup and produce peak densities controlled by the imposed saturation scale. The paper concludes that realistic nonlinear saturation requires additional physics beyond the present closure.

Numerical Tests and Comparison to Streaming Instability

The numerical tests employ a one-dimensional pseudo-spectral solver using an implicit-explicit (IMEX) time-integration strategy, which treats the linear, constant coefficient part of the viscous operator implicitly while evaluating all remaining terms explicitly. The results show that for the diffusion slope driven instability, the dominant Fourier mode in noise initial conditions converges toward the fastest-growing eigenmode predicted by linear theory. However, in nonlinear evolution, the model does not qualitatively reproduce multidimensional streaming-instability simulations where dense filaments saturate at finite amplitude. The authors point out that their model lacks a radial gas pressure gradient, which is generally included in simulations of nonlinear streaming instability to set the background equilibrium drift.

Proposed Modifications and Conclusions

The authors propose a piecewise diffusion law D(Σ) = D0fpw(Σ) and viscosity law ν(Σ) = ν0fpw(Σ) with fpw(Σ) defined piecewise to enforce saturation at high dust densities. This modification eliminates blowup, leading to narrow filaments whose peak densities are controlled by the saturation scale Σsat. The paper supports the robustness of diffusive instabilities as a linear organizing mechanism whenever dust diffusion decreases sufficiently steeply with density, while reiterating that realistic nonlinear outcomes depend sensitively on how turbulent diffusion behaves at high dust loading and on multidimensional saturation effects not captured by their one-dimensional framework.

How it works

The model is built upon an advection-diffusion equation where the flux F is related to density and velocity fluctuations, defining vdiff ≡ −F/ρd. The momentum equation for the dust fluid incorporates terms arising from dust viscosity T, drag forces, and a pressure term derived from the diffusion relation. The effective dust velocity v = v + vdiff is used to simplify the momentum equation into a form where transport coefficients are explicitly dependent on density.

The system is analyzed by linearizing around a constant background state to find the dispersion relation n4 + a3n3 + a2n2 + a1n + a0 = 0. The instability is driven by the density-dependence of transport coefficients, meaning regions of enhanced density experience reduced effective diffusion, creating positive feedback.

Key Findings and Limitations

The paper identifies several key points regarding the instability's behavior. These include:

1. Linear Criterion: The basic criterion for diffusion-slope-driven growth is βdiff ≲ −2 for small dust stopping times.

2. Nonlinear Behavior: For power-law closures D ∝ Σ βdiff with βdiff < 0, the nonlinear evolution does not saturate but instead drives finite-time collapse into increasingly sharp spikes.

3. Stabilization Attempts: A piecewise closure where the negative diffusion slope operates only over a finite density interval eliminates blowup and produces peak densities controlled by the imposed saturation scale Σsat.

4. Model Mismatch: The model does not qualitatively reproduce nonlinear outcomes of multidimensional streaming-instability simulations because it lacks a radial gas pressure gradient and other multidimensional physics like Kelvin-Helmholtz instabilities.

Future Directions

The authors suggest two complementary interpretations for the mismatch. First, the instability may still be relevant as a linear organizing mechanism to set filament spacing during early organization. Second, the closure itself may become invalid at high densities because drift reduction and self-generated turbulence limit further sharpening, mechanisms absent in their 1D framework. Future work should focus on calibrating density-dependent transport closures from streaming instability simulations and extending the model to include a background pressure gradient.

Conclusion

The results support the robustness of diffusive instabilities as a linear organizing mechanism in dusty protoplanetary disks whenever dust diffusion decreases sufficiently steeply with density, while emphasizing that realistic nonlinear outcomes depend sensitively on how turbulent diffusion behaves at high dust loading and on multidimensional saturation effects not captured by their model.

--- Page 1 ---

Draft version February 19, 2026 Typeset using LATEX twocolumn style in AASTeX7.0.1

Diffusive Instabilities in Dusty Disks: Linear Growth and Nonlinear Breakdown Konstantin Gerbig 1, 2, ∗ and Min-Kai Lin (林明楷) 3, 4 Department of Astronomy & Astrophysics, University of Chicago, Chicago, IL 60637, USA Department of Astronomy, Yale University, New Haven, CT 06511, USA Institute of Astronomy and Astrophysics, Academia Sinica, Taipei 10617, Taiwan Physics Division National Center for Theoretical Sciences Department of Astronomy & Astrophysics at the University of Chicago

ABSTRACT We revisit the diffusive instability in dusty disks that arises when the dust mass diffusivity and/or viscosity decreases sufficiently steeply with increasing dust density. Our updated model includes an incompressible, viscous gas that responds azimuthally and couples to the dust through drag. We show that the basic criterion for diffusion-slope-driven instability remains approximately βdiff ≲ −2 for small dust stopping times, with gas feedback providing only modest quantitative changes for parameters motivated by streaming-instability turbulence. We perform nonlinear numerical calculations and confirm linear growth and mode selection toward the fastest-growing wavenumber. However, for power-law closures D ∝ Σ βdiff with βdiff < 0, the nonlinear evolution does not saturate. Instead, steepening gradients amplify the nonlinear dust-pressure term and drive finite-time collapse into increasingly sharp spikes.

1. Improved system can predict filament width and peak density self-consistently by incorporating a physical saturation scale Sigma sat into transport closures like D(Sigma) = D0 fpw(Sigma). This allows for a more accurate representation of nonlinear outcomes compared to the current model which does not saturate at high amplitudes but instead collapses into sharp spikes.

2. Improved system can distinguish between linear organizing mechanisms and nonlinear saturation effects in dusty turbulence. By running noise initial condition simulations, the system can confirm that the instability could help set filament spacing during early organization of a turbulent dust field while acknowledging that realistic nonlinear saturation requires additional physics beyond the present closure.

3. Improved system can identify when current one-dimensional models are inadequate for simulating multidimensional astrophysical phenomena. The AI can flag situations where instabilities like Kelvin-Helmholtz instabilities or eddy diffusion in nonlinear streaming turbulence become relevant, which are not captured by the current framework, as noted in Multidimensional simulations show that parasitic instabilities also contribute to saturation.

We revisit the diffusive instability in dusty disks that arises when the dust mass diffusivity and/or viscosity decreases sufficiently steeply with increasing dust density. Our updated model includes an incompressible, viscous gas that responds azimuthally and couples to the dust through drag. We show that the basic criterion for diffusion-slope-driven instability remains approximately β diff-2 for small dust stopping times, with gas feedback providing only modest quantitative changes for parameters motivated by streaming-instability turbulence. We perform nonlinear numerical calculations and confirm linear growth and mode selection toward the fastest-growing wavenumber. However, for power-law closures D proportional toΣ β diff with β diff<0, the nonlinear evolution does not saturate. Instead, steepening gradients amplify the nonlinear dust-pressure term and drive finite-time collapse into increasingly sharp spikes. Motivated by the absence of multidimensional saturation channels in our 1D framework, we test a simple piecewise closure in which the negative diffusion slope operates only over a finite density interval. This modification eliminates blowup and produces peak densities controlled by the imposed saturation scale. Our results support diffusive instabilities as a linear organizing mechanism in dusty turbulence, while highlighting that realistic nonlinear saturation requires additional physics beyond the present closure.