pith. sign in

arxiv: 2604.08020 · v2 · submitted 2026-04-09 · 🌌 astro-ph.SR · physics.space-ph

Chromospheric turbulence as a regulator of stellar wind mass flux

Munehito Shoda , Tom Van Doorsselaere , Allan Sacha Brun This is my paper

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

classification 🌌 astro-ph.SR physics.space-ph
keywords chromospheric turbulencestellar wind mass fluxAlfvén wave driven windscoronal magnetic field scalingwave energy dissipationsolar and stellar windsmass loss regulation
0
0 comments X

The pith

Suppressing chromospheric turbulence increases stellar wind mass flux by up to an order of magnitude.

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

The paper examines why wave-driven models of stellar winds underpredict observed mass-loss rates compared to X-ray flux measurements. It runs one-dimensional simulations that turn chromospheric turbulence dissipation on or off via a simple switch. Suppressing turbulence allows more wave energy to reach the corona, raising particle density and wind speed and thereby lifting the overall mass flux, especially at moderate magnetic field strengths. This change alone brings the simulated mass flux into line with the observed scaling against coronal magnetic field strength. The result matters because accurate wind mass flux governs long-term stellar spin-down, planetary atmosphere erosion, and space weather around other stars.

Core claim

In one-dimensional wave-driven wind simulations comparing cases with and without chromospheric turbulence suppression, suppressing the turbulence leads to a systematic increase in the coronal particle flux and the wind mass flux by up to an order of magnitude, particularly in regions of moderately strong magnetic field. This arises from changes in the Poynting flux at the coronal base and in the asymptotic wind speed. The model with turbulence suppression reproduces the observed empirical scaling between coronal magnetic field strength and mass flux without invoking additional energy input mechanisms such as interchange reconnection.

What carries the argument

The turbulence-suppression switch inside one-dimensional Alfvén-wave-driven wind models that prevents chromospheric dissipation and thereby increases energy transmission to the corona.

If this is right

  • Coronal particle flux rises systematically when chromospheric turbulence is suppressed.
  • Wind mass flux increases by up to an order of magnitude especially in moderately strong magnetic field regions.
  • The simulated scaling of mass flux with coronal magnetic field strength matches observations.
  • Both the Poynting flux delivered to the coronal base and the terminal wind speed change.
  • No additional energy input such as interchange reconnection is required to recover the observed relation.

Where Pith is reading between the lines

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

  • Multi-dimensional simulations could test whether the simple on-off switch overestimates the mass-flux boost once realistic turbulent cascades are resolved.
  • Stellar evolution codes that adopt this regulation would predict slower spin-down for active stars than current wave-only models imply.
  • Spectroscopic measurements of chromospheric line broadening in stars with known mass-loss rates could directly test whether higher turbulence reduces wind flux as predicted.

Load-bearing premise

The one-dimensional wave-driven simulations with a simple turbulence-suppression switch accurately capture the three-dimensional, multi-scale dissipation physics that occurs in real stellar chromospheres.

What would settle it

A measurement or simulation showing that chromospheric turbulence levels do not inversely correlate with observed wind mass flux once magnetic field strength is fixed, or that the empirical magnetic-field-to-mass-flux scaling still requires extra coronal heating even after turbulence suppression.

Figures

Figures reproduced from arXiv: 2604.08020 by Allan Sacha Brun, Munehito Shoda, Tom Van Doorsselaere.

Figure 2
Figure 2. Figure 2: Scatter plot of the coronal magnetic field strength (𝐵𝑟,cb, horizontal axis) versus the coronal flux-tube expansion factor ( 𝑓exp = lim𝑟→∞ 𝑓cor (𝑟 ), vertical axis) used as simulation input in this study. Each symbol represents an individual flux tube. The three lines indicate contours of constant 𝐵𝑟,cb/ 𝑓exp, with the solid, dashed, and dotted lines corresponding to 𝐵𝑟,cb/ 𝑓exp = 1 G, 0.25 G, and 4 G, res… view at source ↗
Figure 1
Figure 1. Figure 1: Top panel: Comparison between the radial magnetic field 𝐵𝑟 distri￾bution derived from the PFSS extrapolation (diamonds) and that used as input in the simulation (solid lines). Different colors represent different magnetic field lines. Bottom panel: Radial distribution of the corresponding flux-tube expansion factor. The color coding is the same as in the top panel. which shows the consistency with the radi… view at source ↗
Figure 3
Figure 3. Figure 3: Radial profiles of time-averaged physical quantities with chromospheric turbulence suppression (solid lines) and without suppression (dashed lines). Panels (a) and (b) show the mass density ⟨𝜌⟩ and temperature ⟨𝑇⟩, respectively. Panel (c) displays the radial outflow velocity ⟨𝑣𝑟 ⟩ (black lines) together with the Alfvén speed (blue lines) and sound speed (red lines). Panels (d)–(f) present the Elsässer vari… view at source ↗
Figure 4
Figure 4. Figure 4: Coronal magnetic field–coronal particle flux relation obtained from the simulations (red diamonds), compared with the observational scaling relation of Wang (2020) (black solid line) and its half-scaled version (black dashed line). Observational data from Stansby et al. (2021) are also shown as blue plus symbols. Top: case with chromospheric turbulence suppression. Bottom: case without suppression. As show… view at source ↗
Figure 6
Figure 6. Figure 6: Comparison of the Alfvén-wave damping length 𝐿D (defined in Eq. (39)) between models with and without chromospheric turbulence sup￾pression. The horizontal axis shows the values obtained with suppression, and the vertical axis those without suppression. The colour of each symbol denotes the magnetic-field strength at the coronal base. The black solid line marks the one-to-one relation. modifying the distri… view at source ↗
Figure 5
Figure 5. Figure 5: Comparison of simulation outputs obtained with and without chro￾mospheric turbulence suppression. The horizontal axis represents the values from the simulations with turbulence suppression, and the vertical axis rep￾resents those from the simulations without suppression. The top, middle, and bottom panels show the coronal particle flux, the Poynting flux at the coronal base, and the asymptotic wind speed, … view at source ↗
Figure 8
Figure 8. Figure 8: Schematic illustration of the toy model for energy transport from the photosphere to the corona. The horizontal axis represents the height above the stellar surface, and the vertical axis shows the unsigned Poynt￾ing flux. The vertical dashed line indicates the location of the transition region. The blue line shows the height profile of the Poynting flux carried by upward-propagating Alfvén waves (𝐹 + A ),… view at source ↗
read the original abstract

The mass flux of solar and stellar winds is a key quantity for stellar evolution and space weather, yet its physical regulation mechanism remains an unsolved problem. In particular, conventional Alfv\'en wave--driven models that self-consistently connect the stellar surface to the stellar wind fail to reproduce the observed scaling between stellar X-ray flux and mass-loss rate, a discrepancy that can be largely attributed to the dissipation of a substantial fraction of the wave energy by chromospheric turbulence. To address this issue, we aim to clarify the role of chromospheric turbulence in regulating the stellar wind mass flux. We perform one-dimensional wave-driven wind simulations, comparing cases with and without chromospheric turbulence suppression to assess its impact on coronal and wind properties. We find that suppressing chromospheric turbulence leads to a systematic increase in the coronal particle flux, and hence the wind mass flux, by up to an order of magnitude, particularly in regions of moderately strong magnetic field. This behavior arises from a combination of changes in the Poynting flux at the coronal base and in the asymptotic wind speed. Furthermore, the model with chromospheric turbulence suppression reproduces the observed empirical scaling between coronal magnetic field strength and mass flux without invoking additional energy input mechanisms such as interchange reconnection. These results identify the chromospheric turbulence as a key factor in regulating stellar wind mass flux and highlight the importance of incorporating its effects in models that connect the stellar surface and the stellar wind.

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

2 major / 2 minor

Summary. The paper claims that one-dimensional Alfvén wave-driven stellar wind simulations show that suppressing chromospheric turbulence increases coronal particle flux and wind mass flux by up to an order of magnitude (especially at moderate magnetic field strengths), arising from changes in base Poynting flux and asymptotic wind speed; the turbulence-suppressed case reproduces the observed empirical scaling between coronal magnetic field strength and mass flux without requiring additional energy input such as interchange reconnection.

Significance. If the quantitative difference between the two simulation cases holds under more realistic modeling, the result would be significant for stellar wind theory: it identifies chromospheric turbulence dissipation as a primary regulator of mass flux, potentially resolving the long-standing mismatch between standard wave-driven models and observed X-ray–mass-loss scalings. The explicit with/without turbulence comparison and the direct attempt to match an empirical B_cor–mass-flux relation are clear strengths that could guide future multi-dimensional modeling.

major comments (2)
  1. [Methods (simulation setup and turbulence suppression)] The central quantitative result—an order-of-magnitude increase in particle flux when turbulence is suppressed—rests entirely on the difference produced by a phenomenological on/off switch for turbulence dissipation inside a strictly 1-D radial wave-driven code (described in the methods section on simulation setup). Real chromospheric turbulence involves 3-D geometry, multi-scale cascades, and spatially intermittent dissipation; the 1-D switch therefore cannot be assumed to capture the relevant physics, making the reported flux increase and the reproduced scaling load-bearing on an untested modeling assumption.
  2. [Results (particle flux and scaling figures)] No error bars, convergence tests with respect to grid resolution or wave spectrum, or direct validation against multi-dimensional chromospheric simulations are reported for the particle-flux differences or the B_cor–mass-flux scaling (results section). Without these, it is impossible to assess whether the claimed systematic increase and scaling match are robust or artifacts of the 1-D discretization.
minor comments (2)
  1. [Abstract and Results] The abstract and results text refer to 'systematic changes' and 'reproduction' of the empirical scaling, but no quantitative measure (e.g., fitted slope, R², or residual) of the agreement with observations is provided.
  2. [Methods] Notation for the turbulence suppression parameter and the precise form of the dissipation term should be defined explicitly with an equation number rather than described only in prose.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their detailed and constructive report. We address each of the major comments below and outline the revisions we will make to the manuscript.

read point-by-point responses

  1. Referee: [Methods (simulation setup and turbulence suppression)] The central quantitative result—an order-of-magnitude increase in particle flux when turbulence is suppressed—rests entirely on the difference produced by a phenomenological on/off switch for turbulence dissipation inside a strictly 1-D radial wave-driven code (described in the methods section on simulation setup). Real chromospheric turbulence involves 3-D geometry, multi-scale cascades, and spatially intermittent dissipation; the 1-D switch therefore cannot be assumed to capture the relevant physics, making the reported flux increase and the reproduced scaling load-bearing on an untested modeling assumption.

    Authors: We agree that the 1D phenomenological approach has limitations and does not fully capture the complexities of 3D chromospheric turbulence. Our model uses a simplified on/off switch to isolate the effect of turbulence dissipation on the wave energy available for coronal heating and wind acceleration, following approaches in prior 1D Alfvén wave-driven wind studies. This allows us to demonstrate the sensitivity of the mass flux to this process. We will revise the methods and discussion sections to more explicitly state the assumptions and limitations of this approach, and to clarify that the results highlight the potential regulatory role of chromospheric turbulence rather than providing a precise quantitative prediction for real stars. We note that addressing full 3D effects would require a different modeling framework, which is beyond the scope of this work but motivates future studies. revision: partial

  2. Referee: [Results (particle flux and scaling figures)] No error bars, convergence tests with respect to grid resolution or wave spectrum, or direct validation against multi-dimensional chromospheric simulations are reported for the particle-flux differences or the B_cor–mass-flux scaling (results section). Without these, it is impossible to assess whether the claimed systematic increase and scaling match are robust or artifacts of the 1-D discretization.

    Authors: We appreciate this observation. Although we conducted internal tests for numerical convergence during the development of the simulations, these were not reported in the manuscript. In the revised version, we will include a dedicated subsection on numerical convergence, presenting tests with varying grid resolutions and wave spectrum parameters, along with error bars on the key figures derived from these variations. For validation against multi-dimensional simulations, we will add references to existing 3D chromospheric turbulence studies and discuss how our 1D results align qualitatively with them, while acknowledging that direct quantitative validation is not feasible within the current 1D framework. revision: yes

Circularity Check

0 steps flagged

No circularity: simulation outcomes are independent of inputs

full rationale

The paper reports results from one-dimensional wave-driven wind simulations that compare two explicit cases (chromospheric turbulence active vs. suppressed via a switch). The reported increase in coronal particle flux and reproduction of the observed B_cor–mass-flux scaling are direct numerical outputs of those runs, not quantities fitted to data and then relabeled as predictions, nor definitions that presuppose the result. No load-bearing self-citation chain, uniqueness theorem, or ansatz smuggling is invoked to force the central claim; the derivation chain therefore remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract supplies insufficient detail to enumerate specific free parameters or axioms; the central comparison rests on the unstated numerical implementation of turbulence suppression and the standard assumptions of 1D Alfvén-wave wind theory.

pith-pipeline@v0.9.0 · 5560 in / 1060 out tokens · 75332 ms · 2026-05-10T17:41:24.309477+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]

    Twelfth International Solar Wind Conference , year = 2010, editor =

    AirapetianV.S.,GlocerA.,KhazanovG.V.,LoydR.O.P.,FranceK.,Sojka J., Danchi W. C., Liemohn M. W., 2017, ApJ, 836, L3 Airapetian V. S., Jin M., Lüftinger T., Boro Saikia S., Kochukhov O., Güdel M., Van Der Holst B., Manchester W. I., 2021, ApJ, 916, 96 Alazraki G., Couturier P., 1971, A&A, 13, 380 Altschuler M. D., Newkirk G., 1969, Sol. Phys., 9, 131 Alvara...

    work page doi:10.1063/1.3395870 2017

  2. [2]

    D., et al., 2023, Nature, 618, 252 Belcher J

    AIP, pp 11–14, doi:10.1063/1.4810977 Bale S. D., et al., 2023, Nature, 618, 252 Belcher J. W., 1971, ApJ, 168, 509 Berger T. E., Title A. M., 1996, ApJ, 463, 365 Berger T. E., Title A. M., 2001, ApJ, 553, 449 Boro Saikia S., Jin M., Johnstone C. P., Lüftinger T., Güdel M., Airapetian V. S., Kislyakova K. G., Folsom C. P., 2020, A&A, 635, A178 Breu C., Pet...

    work page doi:10.1063/1.4810977 2023