Dust cloud lifetimes of Scallop-shell stars

Luke Bouma; Moira M. Jardine; Simon Daley-Yates

arxiv: 2604.27116 · v1 · submitted 2026-04-29 · 🌌 astro-ph.SR

Dust cloud lifetimes of Scallop-shell stars

Simon Daley-Yates , Moira M. Jardine , Luke Bouma This is my paper

Pith reviewed 2026-05-07 09:53 UTC · model grok-4.3

classification 🌌 astro-ph.SR

keywords M-dwarfsstellar prominencesdust cloudsscallop-shell starslight curve dipsMHD simulationscentrifugal breakoutsTESS photometry

0 comments

The pith

Dust trapped in stellar prominences around fast-rotating M-dwarfs can last for tens of rotations via repeated slingshot ejections.

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

The paper runs a two-dimensional magnetohydrodynamic simulation of a rapidly rotating M-dwarf to track dust grains injected once into cool, magnetically confined gas clouds. A passive tracer shows that centrifugal breakouts repeatedly eject material while chromospheric evaporation adds fresh gas from below, producing an exponential decay in dust content with a half-life of roughly six stellar rotations. This decay is presented as a lower limit under the assumption of perfect dust-gas coupling. The resulting synthetic light-curve signatures reproduce the gradual fading of dips seen in TESS and K2 observations of scallop-shell stars. The work therefore connects the long-lived photometric features directly to the cyclic slingshot dynamics of the prominences.

Core claim

For a star with radius 0.6 solar radii, mass 0.3 solar masses and rotation period 0.32 days, the dust content in a prominence decays exponentially with a minimum half-life of approximately six stellar rotations. Synthetic velocity-phase maps show a single, phase-locked feature that fades steadily, matching observed river plots. The simulation distinguishes three classes of features: persistent non-decaying dips from quiescent prominences below co-rotation, gradually fading dips from slingshot prominences near co-rotation, and abrupt disappearances from reconnection events.

What carries the argument

A passive tracer field advected with the MHD flow that records the single injection and subsequent removal of collisionally charged dust grains during recurrent centrifugal breakouts.

If this is right

Dust-bearing prominences near co-rotation produce gradually fading photometric features that can still be detected after tens of rotations.
Quiescent prominences below co-rotation yield stable, non-decaying dips.
Magnetic reconnection events produce sudden, complete loss of the photometric signal.
The observed longevity of scallop-shell morphology is a natural outcome of the repeated slingshot cycle rather than a static structure.

Where Pith is reading between the lines

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

If real dust decouples from the gas during ejection, actual lifetimes could be shorter than the simulated lower limit.
Three-dimensional geometry or different stellar parameters might change the ejection frequency and therefore the observed decay rate.
The same slingshot mechanism may operate on other rapidly rotating cool stars that show transient absorption features.

Load-bearing premise

Dust and gas remain fully coupled at all times, so that any gas loss immediately removes the associated dust and sets a lower bound on cloud lifetime.

What would settle it

High-cadence photometric monitoring of an individual scallop-shell star that shows whether the amplitude of the phase-locked dips declines exponentially with a half-life near six rotations or remains constant over many tens of rotations.

Figures

Figures reproduced from arXiv: 2604.27116 by Luke Bouma, Moira M. Jardine, Simon Daley-Yates.

**Figure 1.** Figure 1: Development of cyclic ejection of gas clouds. From left to right, the five stages are: (i) heating concentrated at low heights evaporates plasma from the chromosphere (shown as the lower blue horizontal line); (ii) at larger heights there is insufficient heating to prevent cooling; (iii) condensations form and rain out of the cloud either downwards (below the co-rotation radius which is shown as the black … view at source ↗

**Figure 2.** Figure 2: Two dimensional plots of the dust structure around our simulated star. The top is the initial state with uniform tracers in the prominence region. Bottom is the final state. The orange contours depict the density distribution and are equidistant in log space. The black lines show the magnetic field geometry and the colour map shows the tracer, and therefore dust, concentration. The dust bearing prominence … view at source ↗

**Figure 3.** Figure 3: Tracer decay. Here we integrate the tracer concentration over the whole simulation grid, for all outputs. We then normalise each output to the total integrated tracers from the first time step. This way we can see how the action of cyclic slingshot ejections remove tracers, and therefore dust from the prominence (solid line). We then fit a exponential (dashed line) to this decay and work out the dust prom… view at source ↗

**Figure 4.** Figure 4: Dynamic spectrum of the prominence gas, which is illustrative of mock Hα tracks. The colour map is normalised between zero and one, as the absolute value is not important, the rate at which the feature travels through the line centre gives the distance of from the rotational axis. The faster the feature crosses, the further it is from the axis (see Jardine et al. (2020) for examples of this method). Each c… view at source ↗

**Figure 5.** Figure 5: River plot from the dynamic spectrum view at source ↗

read the original abstract

We investigate the survival of dust trapped in magnetically confined cool gas clouds (or {\it prominences}) around rapidly rotating M-dwarfs exhibiting the ``scallop-shell'' light-curve morphology. Using a two-dimensional magnetohydrodynamic simulation, we extend previous coronal prominence models to include a passive tracer field to allow for a single injection of collisionally charged dust grains. The tracer evolution reveals how recurrent centrifugal breakouts--the slingshot process--remove dust and gas from the prominence while chromospheric evaporation replenishes gas from below. For our simulated star, which has $R_{\ast} = 0.6 R_{\odot}$, $M_{\ast} = 0.3 M_{\odot}$, and $P_{\ast} = 0.32$ days, the resulting dust content decays exponentially with a minimum half-life of approximately 6 stellar rotations, representing a lower limit set by our assumption of fully coupled dust and gas dynamics. Synthetic velocity-phase diagnostics show a single, phase-locked feature that fades steadily, reproducing the behaviour of dips seen in TESS and K2 light curves. Comparison with observed river plots suggests a natural classification: (i) persistent, non-decaying features formed by quiescent prominences below co-rotation; (ii) gradually fading features produced by slingshot prominences near co-rotation; and (iii) abrupt disappearances linked to magnetic reconnection and flare-driven ejections. These results demonstrate that dust-bearing prominences--undergoing repeated slingshots--can persist for tens of rotations, linking the observed longevity of the scallop-shell photometric features with the dynamic cycle of prominence slingshot ejections.

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

This 2D MHD run with a passive dust tracer shows exponential decay in prominence dust with a minimum half-life of 6 rotations under full coupling, giving a concrete mechanism for scallop-shell feature longevity but only as a lower limit.

read the letter

The paper adds a passive dust tracer to existing 2D coronal prominence models and integrates it forward under MHD. For the chosen M-dwarf parameters the dust content falls exponentially because repeated centrifugal breakouts remove material while chromospheric evaporation resupplies gas. The resulting half-life is stated as a minimum of roughly six rotations, and the synthetic phase-velocity signals produce a steadily fading, phase-locked dip that resembles some TESS and K2 river plots. They also sort observed features into three groups: persistent ones below co-rotation, gradually fading ones near co-rotation, and abrupt ones tied to reconnection events. That quantitative decay time and the explicit three-way split are the new pieces relative to the earlier prominence literature they cite.

Referee Report

2 major / 2 minor

Summary. The paper investigates dust survival in magnetically confined prominences around rapidly rotating M-dwarfs showing scallop-shell light curves. Using a 2D MHD simulation with a passive tracer for collisionally charged dust under the assumption of full dust-gas coupling, it reports exponential decay of dust content with a minimum half-life of ~6 stellar rotations for parameters R*=0.6 R_sun, M*=0.3 M_sun, P*=0.32 d. This is presented as a lower limit, with recurrent centrifugal breakouts (slingshots) removing material while chromospheric evaporation replenishes gas. Synthetic diagnostics reproduce phase-locked fading, and comparison to TESS/K2 river plots yields a three-category classification of features: persistent (quiescent prominences), gradually fading (slingshot prominences), and abrupt (reconnection/flare ejections). The central claim is that dust-bearing prominences undergoing repeated slingshots can persist for tens of rotations, linking the observed longevity to the dynamic ejection cycle.

Significance. If the result holds, the work supplies a dynamical, parameter-free forward model connecting prominence slingshot cycles to the multi-rotation persistence of scallop-shell photometric dips, offering a physical basis for feature classification without invoking static structures. The forward integration of the MHD equations with an external passive tracer (no inverted parameters) is a methodological strength.

major comments (2)

[simulation results and discussion of assumptions] The reported minimum half-life of ~6 rotations and the claim of persistence over tens of rotations rest on a single 2D MHD run with perfect dust-gas coupling and fixed stellar parameters; the manuscript does not present tests varying the coupling timescale or exploring 3D field topology/reconnection, leaving the sensitivity of the removal rate (and thus whether the value is truly a lower limit) unquantified.
[comparison with observations] The three-category classification of observed river plots (persistent non-decaying, gradually fading, abrupt disappearances) is introduced via post-hoc comparison with TESS/K2 data rather than as a direct prediction from the simulation, which models only the slingshot case; this weakens the claimed link between the simulated dynamics and the full range of observed morphologies.

minor comments (2)

[abstract] The abstract states the half-life as 'approximately 6 stellar rotations' but does not specify the exact fitting procedure or time window used to extract the exponential decay constant from the tracer evolution.
[methods] Notation for the passive tracer field and its injection is introduced without an explicit equation defining its evolution equation or normalization.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for their constructive comments and positive evaluation of the methodological approach. We address each major comment below, with revisions made to clarify limitations and strengthen the presentation of results.

read point-by-point responses

Referee: The reported minimum half-life of ~6 rotations and the claim of persistence over tens of rotations rest on a single 2D MHD run with perfect dust-gas coupling and fixed stellar parameters; the manuscript does not present tests varying the coupling timescale or exploring 3D field topology/reconnection, leaving the sensitivity of the removal rate (and thus whether the value is truly a lower limit) unquantified.

Authors: We agree that the half-life and persistence claims derive from a single 2D simulation assuming perfect dust-gas coupling. This choice isolates the slingshot mechanism to establish a baseline lower limit under ideal coupling. We have added a dedicated paragraph in the discussion section explaining the rationale for full coupling (as a conservative case where dust follows gas removal exactly) and qualitatively addressing how weaker coupling or 3D reconnection might alter rates, while explicitly stating that a full parameter exploration lies beyond the present scope due to computational cost. The lower-limit framing is retained with these caveats. revision: partial
Referee: The three-category classification of observed river plots (persistent non-decaying, gradually fading, abrupt disappearances) is introduced via post-hoc comparison with TESS/K2 data rather than as a direct prediction from the simulation, which models only the slingshot case; this weakens the claimed link between the simulated dynamics and the full range of observed morphologies.

Authors: The simulation models only the slingshot case, which we associate with the gradually fading morphology. The classification is an observational framework motivated by contrasting the simulated fading behavior against TESS/K2 river plots, where persistent features are interpreted as stable sub-co-rotation prominences and abrupt ones as reconnection events outside the model scope. We have revised the text in Section 4 to explicitly distinguish the simulation's direct prediction (gradual fading via slingshots) from the broader interpretive scheme, avoiding any implication of direct predictions for all categories. revision: yes

standing simulated objections not resolved

The sensitivity of the dust removal rate to variations in coupling timescale and 3D field topology/reconnection remains unquantified, as no additional simulations were performed.

Circularity Check

0 steps flagged

Forward 2D MHD simulation with passive tracer produces independent dust lifetime output

full rationale

The central result is the exponential decay of the dust tracer (half-life ~6 rotations) obtained by integrating the MHD equations forward in time for fixed stellar parameters (R*=0.6 R_sun, M*=0.3 M_sun, P*=0.32 d) under the explicit assumption of perfect dust-gas coupling. This is an emergent numerical output, not a fitted parameter, not a redefinition of inputs, and not dependent on a self-citation chain for its value. The paper states the result is a lower limit due to the coupling assumption and 2D geometry; no load-bearing step reduces by construction to the inputs or to prior work by the same authors. The classification of observed features is interpretive comparison, not a derivation. The chain is therefore self-contained as a controlled numerical experiment.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard MHD plasma equations plus one key modeling choice (perfect dust-gas coupling) that is introduced to obtain the reported lower-limit lifetime; no new particles or forces are postulated.

axioms (2)

standard math The equations of ideal magnetohydrodynamics govern the evolution of the stellar corona and prominences.
Invoked as the basis for the 2D simulation framework extending previous models.
ad hoc to paper Dust grains remain fully coupled to the gas flow at all times.
Explicitly stated as the assumption that sets the reported half-life as a lower limit.

pith-pipeline@v0.9.0 · 5604 in / 1446 out tokens · 79517 ms · 2026-05-07T09:53:52.510311+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

1 extracted references · 1 canonical work pages

[1]

A., 2004, ApJ, 612, 519 Bouma L

Ansdell M., et al., 2016, ApJ, 816, 69 Aulanier G., Demoulin P., 1998, A&A, 329, 1125 Aulanier G., Demoulin P., van Driel-Gesztelyi L., Mein P., Deforest C., 1998, A&A, 335, 309 van Ballegooijen A. A., 2004, ApJ, 612, 519 Bouma L. G., Jardine M. M., 2025, ApJ, 988, L3 Bouma L. G., Jayaraman R., Rappaport S., Rebull L. M., Hillen- brand L. A., Winn J. N., ...

work page 2016

[1] [1]

A., 2004, ApJ, 612, 519 Bouma L

Ansdell M., et al., 2016, ApJ, 816, 69 Aulanier G., Demoulin P., 1998, A&A, 329, 1125 Aulanier G., Demoulin P., van Driel-Gesztelyi L., Mein P., Deforest C., 1998, A&A, 335, 309 van Ballegooijen A. A., 2004, ApJ, 612, 519 Bouma L. G., Jardine M. M., 2025, ApJ, 988, L3 Bouma L. G., Jayaraman R., Rappaport S., Rebull L. M., Hillen- brand L. A., Winn J. N., ...

work page 2016