The effect of dust on vortices II: Streaming instabilities

This paper demonstrates that a cousin of the streaming instability remains active within vortices by extending resonant drag instability theory to non-modal waves in a vortex-following shearing box, thereby strengthening the case for vortex-induced planetesimal formation while clarifying the nature of dust-driven vortex instabilities.

The Gist

Imagine a giant, swirling whirlpool in a cosmic ocean. This is a vortex in a protoplanetary disk—a swirling cloud of gas and dust where new planets are born. For a long time, scientists thought these whirlpools were like cosmic vacuum cleaners, sucking up dust and clumping it together to form the building blocks of planets (called planetesimals).

But there was a problem. The dust didn't just sit there; it pushed back on the gas. This "pushback" made the whirlpool wobble and eventually fall apart before it could gather enough dust to make a planet. It was like trying to build a sandcastle while the ocean waves kept knocking it down.

This paper is the second part of a two-part study by Nathan Magnan and Henrik Latter. They asked a crucial question: Is there a way for the dust and the whirlpool to work together instead of fighting each other?

Here is the breakdown of their discovery, explained simply:

1. The Problem: The "Tug-of-War"

In the first paper, they showed that if dust just sits in a whirlpool, it eventually destroys the structure. But in this paper, they looked at a specific, chaotic type of turbulence called the Streaming Instability (SI).

Think of the Streaming Instability like a crowd surge. If you have a crowd of people (dust) and a moving walkway (gas), and the people start running slightly faster or slower than the walkway, they can bunch up into dense clusters. Usually, this requires very specific conditions (like a lot of dust and very little wind turbulence) to happen.

2. The New Idea: The "Vortex Shearing Box"

The authors realized that studying a giant, swirling whirlpool is mathematically messy. It's like trying to predict the weather by standing on a spinning carousel.

To solve this, they invented a clever mental trick they call the "Vortex Shearing Box."

• The Analogy: Imagine you are riding a bicycle inside a giant, rotating elliptical racetrack. Instead of trying to map the whole track, you build a small, transparent box around yourself that moves with the track's curves.
• The Result: From inside this box, the complex, swirling motion looks like a simple, rhythmic back-and-forth sway. This allowed them to do the math on how tiny ripples in the dust and gas would behave.

3. The Discovery: The "Cosmic Dance"

When they ran the numbers, they found something amazing. The dust and the gas can dance together in a way that creates instability, but it's a very specific kind of dance.

• The Gas: The gas in the vortex creates invisible "lanes" or channels (like lanes on a highway) that flow back and forth.
• The Dust: The dust particles drift along these lanes.
• The Resonance: The authors found that if the dust drifts at just the right speed, it hits a "sweet spot" (resonance). It's like pushing a child on a swing. If you push at the exact right moment in the swing's cycle, the swing goes higher and higher.

In this cosmic dance, the dust pushes the gas, and the gas pushes the dust back, creating a positive feedback loop. The dust clumps get denser, which pushes the gas harder, which makes the clumps even denser.

4. The Big Surprise: It Works in 2D!

Usually, this kind of "clumping instability" (Streaming Instability) needs a third dimension (up and down) to work. It's like a 3D puzzle that won't fit on a flat table.

However, the authors found that inside a vortex, this instability works even in 2D (flat).

• Why? In a normal disk, the "lanes" of gas are straight. In a vortex, the lanes curve and twist. The authors discovered that the curvature of the vortex acts like a special gear that allows the dust to stay in sync with the gas even without a third dimension.
• The Metaphor: Imagine a 2D video game where a character usually needs to jump (3D) to collect coins. But in this specific level (the vortex), the floor itself bounces up and down, letting the character collect coins without ever leaving the ground.

5. What Does This Mean for Planet Formation?

This is a huge deal for our understanding of how planets form.

• The "Metre-Scale Barrier": One of the biggest mysteries in astronomy is how dust grains grow from the size of sand to the size of boulders. They usually get stuck at the size of a meter because they get blown away by gas.
• The Solution: This paper suggests that vortices are the perfect "incubators." They can trap dust, and thanks to this new "Vortex Streaming Instability," they can rapidly clump that dust together into massive, dense piles.
• The Result: These dense piles can then collapse under their own gravity to form planetesimals (the seeds of planets).

6. The Caveats (The "But...")

The authors are careful to note that this is a theoretical model (math on paper), not a full computer simulation of a real, messy universe.

• The "Baby" Phase: They only studied the very beginning of the instability (the linear phase). They don't know yet if the clumps survive long enough to become planets or if they get shredded by turbulence.
• The Scale: The instability happens on very small scales. It's possible that in a real, giant vortex, other forces might interfere.

Summary

In simple terms, Magnan and Latter discovered that whirlpools in space aren't just dust traps; they are dust accelerators.

By creating a special mathematical "box" to follow the swirl, they proved that dust and gas can lock into a rhythmic, self-reinforcing dance. This dance creates dense clumps of dust much faster than previously thought, potentially solving the mystery of how the first building blocks of our solar system were formed. It's a strong argument that vortices are the missing link in the story of how planets are born.

Technical Summary

1. Problem Statement

The formation of planetesimals in protoplanetary discs (PPDs) faces the "metre-scale barrier," where aerodynamic drag causes pebbles to drift into the star before they can grow. While large-scale anticyclonic vortices are known to concentrate dust, previous work (Paper I) suggested that laminar concentration is too slow to trigger gravitational collapse before the vortex is destroyed by the Elliptical Instability (EI).

The Streaming Instability (SI) is a promising mechanism to rapidly concentrate dust, but it typically requires specific conditions (high metallicity, narrow particle size distribution, weak pressure gradients) that are hard to satisfy in standard disc models. Vortices naturally provide high dust density and narrow size distributions, leading to the hypothesis of a "turbulent" vortex pathway where SI operates within vortices.

However, numerical simulations of dusty vortices have shown that dust backreaction triggers an unknown instability that destroys vortices. The nature of this instability remains unclear:

• Is it a form of the Streaming Instability?
• Why does it appear in 2D simulations, given that the classical SI requires a vertical dimension?
• What are the physical mechanisms driving it?

This paper aims to answer these questions through analytical linear stability analysis.

2. Methodology

The authors employ a linear perturbation analysis of dust-laden vortices, extending the framework of Resonant Drag Instability (RDI) theory.

• Background Flow: They utilize the analytical dusty vortex models derived in Paper I, based on Kida's gaseous vortex model. They assume the dust is well-coupled (small Stokes number, $St \ll 1$) and dilute ($\mu = \rho_d/\rho_g \ll 1$).
• Vortex Shearing Box: To handle the spatial complexity of the elliptical vortex flow, the authors construct a novel coordinate system called the "vortex shearing box." Unlike the standard shearing box that follows Keplerian orbits, this box follows a specific streamline of the vortex. This allows for the study of small-wavelength perturbations by treating the background flow as locally uniform within the box.
• Linear Perturbation Theory: The governing two-fluid equations (gas and dust) are linearized. The authors perform a Fourier decomposition in space with a time-dependent wavevector $\mathbf{k}(t)$ to cancel the advection terms caused by the background flow.
• Wave Analysis: Before solving for instability, they characterize the waves propagating in the vortex:
  • Gas Waves: Inertial waves (non-sinusoidal due to time-dependent rotation) and Zonal flows (non-propagating, purely azimuthal velocity perturbations).
  • Dust Waves: Dust density waves (advected by the background drift) and damped velocity waves.
• Resonant Drag Instability (RDI) Generalization: The authors extend standard RDI theory (which relies on resonant frequencies between gas and dust waves) to non-modal waves (Floquet theory). They derive a resonance condition where the Floquet frequencies of gas and dust waves must satisfy $\omega_g = \omega_d + 2n\pi/T_v$.

3. Key Contributions

1. Development of the Vortex Shearing Box: A new analytical framework to study local perturbations in elliptical vortices, allowing for Fourier analysis in a non-uniform background flow.
2. Extension of RDI Theory: Generalizing the resonance condition for instabilities driven by non-sinusoidal, time-periodic waves (Floquet modes) rather than simple harmonic waves.
3. Mechanism Explanation for 2D Instability: Providing a physical explanation for how the Streaming Instability can operate in 2D within vortices, a phenomenon previously unexplained in the literature.

4. Key Results

A. The Nature of the Instability

The analysis confirms that the instability triggered by dust backreaction in vortices is indeed a form of Streaming Instability (specifically, a Resonant Drag Instability).

• 2D Axisymmetric Mode: The dominant instability in 2D is a "Zonal Flow RDI." It arises from the resonance between the non-propagating zonal flow (gas wave) and the dust density wave.
• 3D Modes: In 3D, the instability includes bands corresponding to the standard SI (resonance between inertial waves and dust density waves) and higher-order harmonics ("pimples").

B. The 2D Mechanism (Why it works in 2D)

The paper resolves the paradox of 2D SI activity:

• In standard discs, zonal flows do not propagate and cannot resonate with radial dust drift.
• In elongated vortices, the Coriolis parameter ($d_t\phi - \Omega$) changes sign as the fluid parcel moves from the major to the minor axis of the ellipse.
• This sign change causes the pressure perturbation associated with the zonal flow to alternate between a bump and a trough.
• If the dust drift timescale matches the vortex half-turnover time (resonance condition), a dust parcel is always inside a pressure bump as it drifts, leading to unbounded density growth. This "forward action" survives in 2D.
• Similarly, the "backward reaction" (dust density amplifying the gas flow) survives because the azimuthal drift direction of the dust changes sign in sync with the flow channels.

C. Parameter Space Exploration

• Dust-to-Gas Ratio ($\mu$): Growth rates scale as $\sqrt{\mu}$, consistent with RDI theory.
• Stokes Number ($St$): Unlike the standard SI (where growth scales with $St$ for small particles), the vortex SI growth rate is independent of $St$ in the dilute limit. This is because the azimuthal drift in vortices scales linearly with $St$ (rather than $St^2$ in discs), making the "slow" backward mechanism as fast as the "fast" radial mechanism.
• Aspect Ratio ($\alpha$): The instability is active in the elliptically stable band ($4 < \alpha < 6$). The growth rates and mode structures vary complexly with $\alpha$, but the "Zonal Flow RDI" (vertical bands in $K_x$ space) remains dominant.

D. Comparison with Simulations

• The authors identify their instability as the "unknown instability" seen in 2D and 3D simulations.
• Discrepancies: The analytical growth rates are slower (e-folding time ~40 orbits) than some simulations (orbital timescale), and the preferred wavelengths are smaller.
• Possible Causes: The simulations may use non-Kida vortex profiles, or the simulated instability might be a different, faster mode (potentially involving compressibility/sound waves) that the current incompressible model misses.

5. Significance

• Validation of Vortex Planet Formation: The results strongly support the "turbulent vortex" pathway. Vortices can naturally concentrate dust and trigger an SI-like instability, potentially bridging the metre-scale barrier.
• Resolution of the 2D Paradox: The paper provides the first physical mechanism explaining how SI-like clumping can occur in 2D simulations, attributing it to the unique geometry and time-dependent Coriolis forces of elliptical vortices.
• Theoretical Advancement: By extending RDI theory to non-modal, periodic flows, the authors demonstrate that streaming instabilities are robust to complex background flows, not just simple Keplerian shear.
• Future Directions: The authors note that while the linear phase is understood, non-linear simulations are required to confirm if this instability leads to gravitational collapse and planetesimal formation. They suggest using the "vortex shearing box" to run local simulations to overcome resolution limits.

In conclusion, Magnan and Latter demonstrate that dust-laden vortices are inherently unstable due to a resonant drag instability closely related to the Streaming Instability, offering a robust mechanism for planetesimal formation that functions even in two-dimensional approximations.