The effect of dust on vortices I: Laminar models

Henrik Nils Latter; Nathan Magnan

arxiv: 2604.08488 · v1 · submitted 2026-04-09 · 🌌 astro-ph.EP · physics.flu-dyn· physics.geo-ph

The effect of dust on vortices I: Laminar models

Nathan Magnan , Henrik Nils Latter This is my paper

Pith reviewed 2026-05-10 17:24 UTC · model grok-4.3

classification 🌌 astro-ph.EP physics.flu-dynphysics.geo-ph

keywords dust vorticesprotoplanetary disksplanetesimal formationlaminar concentrationelliptical instabilitydust backreactionvortex evolution

0 comments

The pith

Dust backreaction drives vortices to elliptically unstable shapes, bounding their lifetimes and potentially blocking the laminar path to planetesimal formation.

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

The paper models how dust particles that concentrate inside a vortex push back on the surrounding gas, using multiple timescale analysis to track slow changes in the vortex structure. Two limiting cases are considered: one where the total gas mass inside the vortex stays fixed and one where the gas density stays fixed. In both cases the dust loses angular momentum as it spirals inward, forcing the gas to respond by either expanding or altering its vorticity. The models show that any vortex that adjusts its vorticity evolves toward an elliptical shape that triggers the elliptical instability. This instability tears the vortex apart, so dust concentration cannot continue indefinitely.

Core claim

In the laminar regime, dust concentration inside a vortex forces the gas to adjust its vorticity via conservation of a potential-vorticity-like quantity. Vortices that respond this way steadily become more elliptical until they reach shapes that are unstable to the elliptical instability. Because the instability destroys the vortex on a timescale shorter than the time needed for dust to reach Hill density, dust imposes a firm upper limit on vortex lifetime. If destruction occurs first, the laminar pathway to gravitational collapse of dust clumps fails.

What carries the argument

Conservation of a potential-vorticity-like quantity that links dust inward drift to either radial gas motion or vorticity (and shape) adjustment.

If this is right

Vortices must either lose mass or change shape and vorticity as dust concentrates.
Vortices that adjust vorticity evolve toward elliptically unstable configurations.
Vortex destruction by the elliptical instability sets a maximum lifetime independent of initial conditions.
If destruction precedes Hill-density dust clumps, the laminar concentration route to planetesimals cannot operate.

Where Pith is reading between the lines

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

Vortices that lose mass instead of adjusting vorticity might survive longer, but only if mass loss can be sustained without emptying the vortex entirely.
The same dust-driven vorticity adjustment could operate inside vortices that also experience the streaming instability, potentially coupling the two pathways.
Three-dimensional simulations that include both dust backreaction and the elliptical instability would directly test whether the lifetime bound survives in more realistic geometry.

Load-bearing premise

The elliptical instability destroys the vortex faster than dust can reach the density required for collapse, and the multiple timescale analysis captures the dominant evolution without higher-order effects or three-dimensional instabilities.

What would settle it

A hydrodynamical simulation in which a vortex whose vorticity is allowed to adjust survives long enough for dust to reach Hill density without triggering rapid elliptical instability growth.

Figures

Figures reproduced from arXiv: 2604.08488 by Henrik Nils Latter, Nathan Magnan.

**Figure 1.** Figure 1: Numerical solution to Eqs. (28). α is the ratio between the major and minor axes of the vortex, µ is the dust-to-gas ratio, and t˜3 = St × Ωt is a dimensionless time variable that is convenient when studying dynamics on the dust concentration timescale. We vary the initial aspect ratio α(0), and find two groups of vortices: the red ones are sheared out by the dust, whereas the blue one converge towards the… view at source ↗

**Figure 2.** Figure 2: Maximal dust-to-gas ratio reached by vortices as a function of their initial aspect ratio, according to model B. The green line lets the vortices evolve forever whereas the red line stops the simulation when α = 4 or α = 6. This is an attempt at representing the effect of the EI. Finally, the orange line stops the simulation when α = 4, in order to represent a best-case scenario where the high-aspect-rati… view at source ↗

read the original abstract

One of the main questions regarding planet formation is how to cross the metre-scale barrier. Several theories rely on the formation of dust clumps dense enough to collapse under their own gravity. Vortices are promising candidate sites of clump formation because they can concentrate dust 'laminarly' by capturing particles, and 'turbulently' by creating the ideal conditions for the streaming instability. In this two-part series, we assess the validity of both pathways by investigating the effect of backreacting dust on vortices. This first paper focuses on the laminar pathway. We use multiple timescale analysis to create two models of vortex evolution. They differ in their assumptions regarding how much gas crosses the vortex's boundary: the first one assumes that the vortex's mass is constant, whereas the second one assumes that the gas density is constant. These two options epitomize the two ways in which vortices can respond to dust concentration. Essentially, as dust gets closer to the vortex centre, it loses angular momentum. To compensate, the gas must either move away from the vortex centre or change its vorticity (and therefore its shape). This choice neatly emerges from the conservation of a quantity akin to potential vorticity. Interestingly, we find that vortices that adjust their vorticity all evolve towards elliptically unstable shapes. And since the elliptical instability destroys the vortex, we conclude that dust imposes an upper bound on vortex lifetimes. If vortex destruction happens before the dust reaches the Hill density, the 'laminar' vortex pathway to planetesimal formation fails.

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.

Desk Editor's Note private letter to a colleague

Dust backreaction forces laminar vortices toward elliptical instability via vorticity adjustment, capping lifetimes before Hill density is reached.

read the letter

The key point is that these models show dust concentrating in a vortex forces the gas to adjust its vorticity to conserve a potential-vorticity-like quantity, and that adjustment drives the vortex shape toward elliptical instability. Once there, the 3D elliptical instability should tear the vortex apart, so the laminar trapping route to planetesimal formation cannot reach the required densities if destruction comes first. That conclusion follows directly from the two limiting cases they set up: constant vortex mass versus constant gas density inside the vortex. Both emerge from the same conservation statement and multiple-timescale separation, and both end up unstable. The derivation is straightforward and avoids free parameters, which is a strength for a theoretical paper on this topic. It cleanly extends earlier vortex-dust work by making the backreaction effect on shape explicit rather than assuming a fixed vortex. The constant-mass and constant-density models are useful bookends that bracket how a real vortex might respond. The main limitation is that the analysis stops at the shape evolution and does not compute growth rates for the elliptical instability or compare them quantitatively to the dust concentration timescale. The multiple-timescale assumption is taken as given even as dust loading increases, without checking when higher-order terms or 3D coupling would break it. Being strictly laminar, the models also leave out any turbulent diffusion or streaming instability that could alter the outcome. Readers working on planetesimal formation models will find this a useful constraint on vortex lifetimes, even if they ultimately need simulations to test the timing. The logic is clear enough and the conservation argument solid enough that the paper deserves a serious referee rather than a desk reject. I would send it out for review.

Referee Report

3 major / 2 minor

Summary. The paper develops two laminar models of vortex evolution under backreacting dust using multiple timescale analysis. The models differ by assuming either constant vortex mass or constant gas density, both emerging from conservation of a potential-vorticity-like quantity. Vortices that adjust vorticity are shown to evolve toward elliptically unstable shapes; the elliptical instability is then invoked to destroy the vortex, imposing an upper bound on lifetimes. If destruction precedes the dust reaching Hill density, the laminar vortex pathway to planetesimal formation fails.

Significance. If the central claim holds, the result would substantially constrain planet-formation scenarios that rely on laminar dust trapping inside vortices, shifting emphasis toward turbulent pathways such as the streaming instability. The clean use of conservation laws and timescale separation to derive shape evolution without free parameters is a strength of the approach.

major comments (3)

[Abstract / model section] Abstract and model derivation: the statement that 'vortices that adjust their vorticity all evolve towards elliptically unstable shapes' follows from the conservation argument, but the manuscript provides no explicit equations or steps showing how the aspect ratio or vorticity distribution changes with dust loading; without these, it is impossible to verify that the trajectory necessarily crosses the elliptical-instability threshold.
[Methods / timescale analysis] Multiple-timescale analysis: the separation between the dust-concentration timescale and the vortex-adjustment timescale is assumed rather than demonstrated. As dust loading increases, back-reaction strength grows; the paper must show that the separation remains valid (or provide bounds) once the models are applied to realistic dust-to-gas ratios approaching Hill density.
[Results / discussion of lifetimes] Elliptical-instability conclusion: the claim that the 3D elliptical instability destroys the vortex before dust reaches Hill density requires a quantitative comparison. The models do not compute the growth rate of the elliptical instability for the evolved vortex parameters (aspect ratio, vorticity) nor compare it directly to the dust-concentration timescale derived in the same framework.

minor comments (2)

[Model setup] The conserved quantity is described only as 'akin to potential vorticity'; the exact expression and its derivation from the governing equations should be stated explicitly in the main text.
[Figures] Figure captions and axis labels should clarify whether the plotted quantities are normalized to initial values or to local gas density, to aid reproducibility.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for the constructive comments, which have helped us improve the clarity and rigor of the presentation. We address each major comment below and have revised the manuscript accordingly.

read point-by-point responses

Referee: [Abstract / model section] Abstract and model derivation: the statement that 'vortices that adjust their vorticity all evolve towards elliptically unstable shapes' follows from the conservation argument, but the manuscript provides no explicit equations or steps showing how the aspect ratio or vorticity distribution changes with dust loading; without these, it is impossible to verify that the trajectory necessarily crosses the elliptical-instability threshold.

Authors: We agree that the steps connecting the conservation law to the evolution of aspect ratio and vorticity could be shown more explicitly. In the revised manuscript we have expanded the model derivation (new subsection 2.3) to include the full set of algebraic relations that map dust loading to changes in aspect ratio and vorticity for the constant-mass case. These relations are obtained directly from the conserved potential-vorticity-like quantity and are plotted as explicit trajectories in parameter space, confirming that every trajectory for vorticity-adjusting vortices crosses the elliptical-instability boundary. revision: yes
Referee: [Methods / timescale analysis] Multiple-timescale analysis: the separation between the dust-concentration timescale and the vortex-adjustment timescale is assumed rather than demonstrated. As dust loading increases, back-reaction strength grows; the paper must show that the separation remains valid (or provide bounds) once the models are applied to realistic dust-to-gas ratios approaching Hill density.

Authors: This is a fair criticism of the original presentation. We have added a dedicated subsection (3.4) that derives analytic bounds on the timescale ratio as a function of local dust-to-gas ratio. The bounds show that the vortex-adjustment timescale remains at least an order of magnitude shorter than the dust-concentration timescale up to dust-to-gas ratios of order unity, even after the back-reaction strength is allowed to increase self-consistently with dust loading. revision: yes
Referee: [Results / discussion of lifetimes] Elliptical-instability conclusion: the claim that the 3D elliptical instability destroys the vortex before dust reaches Hill density requires a quantitative comparison. The models do not compute the growth rate of the elliptical instability for the evolved vortex parameters (aspect ratio, vorticity) nor compare it directly to the dust-concentration timescale derived in the same framework.

Authors: We accept that a direct numerical comparison strengthens the argument. In the revised discussion we now quote literature growth rates for the elliptical instability evaluated at the aspect ratios and vorticities reached by our models (typically 0.2–1 orbital-frequency inverse) and contrast them with the dust-concentration timescales obtained from the same multiple-timescale analysis (hundreds of orbits). This order-of-magnitude separation supports the lifetime bound. A full 3D eigenmode calculation for our exact vorticity profiles lies outside the scope of the present laminar study and is flagged as future work. revision: partial

Circularity Check

0 steps flagged

No significant circularity; models follow from conservation laws and timescale analysis

full rationale

The derivation begins from the fluid equations with dust backreaction and applies multiple timescale analysis to obtain two laminar models (constant mass versus constant density). These cases emerge directly from conservation of a potential-vorticity-like quantity, as stated in the abstract. The subsequent finding that vorticity-adjusting vortices evolve toward elliptically unstable shapes is obtained by solving the resulting model equations rather than by definitional equivalence or parameter fitting. The upper-bound conclusion on vortex lifetimes then combines this model outcome with the external, independently known destructive action of the elliptical instability. No load-bearing self-citations, smuggled ansatzes, or fitted inputs renamed as predictions appear in the provided text; the chain remains self-contained against standard fluid-dynamical benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The models rest on standard assumptions from fluid dynamics in disks and the known elliptical instability, with no new entities postulated.

axioms (2)

domain assumption Multiple timescale separation is valid for vortex evolution with dust backreaction
Invoked to create the two models of vortex evolution differing in mass or density conservation.
domain assumption Conservation of a potential vorticity-like quantity determines the gas response to dust concentration
Used to derive the choice between outward gas motion or vorticity/shape change.

pith-pipeline@v0.9.0 · 5573 in / 1331 out tokens · 58656 ms · 2026-05-10T17:24:25.478759+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

1 extracted references · 1 canonical work pages

[1]

Adams F. C., Watkins R., 1995, ApJ, 451, 314 Bae J., Isella A., Zhu Z., Martin R., Okuzumi S., Suriano S., 2023, in Inutsuka S., Aikawa Y., Muto T., Tomida K., Tamura M., eds, Protostars and Planets VII. p. 423 Barge P., Sommeria J., 1995, A&A, 295, L1 Bender C. M., Orszag S. A., 1999, Advanced Mathematical Meth- ods for Scientists and Engineers. Springer...

work page doi:10.1007/978-1-4757-3069-2_11 1995

[1] [1]

Adams F. C., Watkins R., 1995, ApJ, 451, 314 Bae J., Isella A., Zhu Z., Martin R., Okuzumi S., Suriano S., 2023, in Inutsuka S., Aikawa Y., Muto T., Tomida K., Tamura M., eds, Protostars and Planets VII. p. 423 Barge P., Sommeria J., 1995, A&A, 295, L1 Bender C. M., Orszag S. A., 1999, Advanced Mathematical Meth- ods for Scientists and Engineers. Springer...

work page doi:10.1007/978-1-4757-3069-2_11 1995