pith. sign in

arxiv: 2604.15689 · v1 · submitted 2026-04-17 · 🌌 astro-ph.SR

Radiation-driven stellar winds at the fast-slow transition: new hydrodynamic solutions

M.C. Fernandez , R.O.J. Venero , L.S. Cidale , I. Araya , M. Cur\'e This is my paper

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

classification 🌌 astro-ph.SR
keywords solutionsdeltadifferenthydrodynamiclineprofilesstationarywinds
0
0 comments X

The pith

New stable wind solutions fill the elusive gap between fast and δ-slow radiation-driven regimes, showing velocity kinks and distinct H I, He I, Si IV line profiles.

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

Massive stars drive powerful winds by radiation pushing on atoms in their outer layers. The standard way to model this uses three numbers called alpha, delta, and k to describe the force. This leads to fast winds or slower delta-slow winds, but there is a gap in between where no steady wind solutions were known before. This work uses a computer program called ZEUS-3D that simulates the wind changing over time until it settles into a steady state. For a model of a B supergiant star, they found new steady solutions in that gap. Some of these have a sharp bend or kink in how fast the wind is moving at a certain distance from the star. They also calculated what the light from hydrogen, helium, and silicon would look like coming from these winds. The authors suggest that small changes in how ionized the wind is could make the wind switch between different types, which might explain why some stars show changing or patchy winds in observations.

Core claim

We found new stationary solutions in the gap region, alongside their corresponding line profiles, for a typical B supergiant star model. In this model, the new solutions are stable, and some of them present a kink in the velocity profile at a fixed distance from the star, depending on the δ value.

Load-bearing premise

That the modified CAK theory parameters can be chosen such that time-dependent hydrodynamics in ZEUS-3D yields stable stationary solutions in the gap region where previous methods failed, and that ionization perturbations can trigger transitions between regimes.

Figures

Figures reproduced from arXiv: 2604.15689 by I. Araya, L.S. Cidale, M.C. Fernandez, M. Cur\'e, R.O.J. Venero.

Figure 1
Figure 1. Figure 1: Temporal evolution obtained with ZEUS-3D (solid lines) [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Hydrodynamic solutions for T19 model (Teff = 19 000 K), adopting values from 0 < δ < 0.4 and ∆δ = 0.01, using different rotational rates, obtained with ZEUS-3D code. Dashed lines represents the solutions in the gap region. sensitive dependence of the solution on the ionisation parameter in this interval. These new solutions form a continuous transition between the fast and δ-slow regimes. Near the photosph… view at source ↗
Figure 3
Figure 3. Figure 3: Terminal velocity (left panel) and mass-loss rate (right panel) as functions of the [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Synthetic line profiles computed for the three hydrodynamic wind regimes using the T19 model: fast solution ( [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Comparison between synthetic line profiles computed using the hydrodynamic solutions obtained with ZEUS-3D (using the [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗
read the original abstract

Radiation-driven winds of massive stars can be described within the modified CAK theory, which parametrises the radiation force through three key quantities: $\alpha$, $\delta$, and $k$. Different combinations of these parameters, together with rotation, result in three types of stationary solutions, namely fast (or classical), $\delta$-slow, and $\Omega$-slow solutions. The primary objective of this work is to model radiation-driven winds inside the gap region between the fast and $\delta$-slow regimes, where stationary solutions have proven elusive. In addition, we compute synthetic line profiles of H I, He I, and Si IV to illustrate the morphology of different wind regimes. We employ the time-dependent hydrodynamic code ZEUS-3D, capable of obtaining stationary solutions by progressing through an initial solution. Then we compute the line profiles solving the transfer equation for an expanding atmosphere, assuming spherical symmetry in the comoving frame, under non-local thermodynamic equilibrium (NLTE) conditions. We found new stationary solutions in the gap region, alongside their corresponding line profiles, for a typical B supergiant star model. In this model, the new solutions are stable, and some of them present a kink in the velocity profile at a fixed distance from the star, depending on the $\delta$ value. Perturbations in the wind ionisation may trigger transitions between different hydrodynamic regimes and offer a plausible explanation for structured and variable winds. A systematic investigation of these effects will be the subject of future work. Furthermore, we investigate the resulting line profiles from different hydrodynamic solutions and compare them with those predicted by a velocity profile given by a $\beta$-law using the same global wind parameters.

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.

Circularity Check

0 steps flagged

Numerical integration yields emergent stationary solutions with no circular reduction

full rationale

The paper obtains new stationary wind solutions inside the fast–δ-slow gap by evolving the time-dependent hydrodynamic equations (continuity and momentum with modified CAK force) in ZEUS-3D from chosen initial conditions until a time-independent state appears. The final (ρ, v) profiles are therefore determined by the dynamics of the PDE system rather than being predefined by the input parameters α, δ, k or by any fitted quantity. Synthetic line profiles are computed afterward from the resulting velocity law and compared to a β-law; this is a post-processing diagnostic, not a prediction that loops back to the inputs. No load-bearing step invokes a self-citation whose validity depends on the present result, nor is any ansatz or uniqueness theorem smuggled in to force the outcome. The derivation chain is therefore self-contained as a numerical experiment.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The claim depends on the numerical method's ability to converge and the validity of the modified CAK parametrization for the gap parameters.

free parameters (1)
  • α, δ, k
    Parameters parametrizing the radiation force in modified CAK theory; combinations are explored to find gap solutions.
axioms (1)
  • domain assumption The time-dependent code ZEUS-3D can obtain stationary solutions by evolving from an initial state in the gap region.
    This is the key assumption allowing discovery of new solutions where stationary methods failed.

pith-pipeline@v0.9.0 · 5623 in / 1315 out tokens · 65362 ms · 2026-05-10T08:24:01.558916+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

36 extracted references · 36 canonical work pages

  1. [1]

    Abbott, D. C. 1982, ApJ, 259, 282

    work page 1982

  2. [2]

    2018, MNRAS, 477, 755

    Araya, I., Curé, M., ud-Doula, A., Santillán, A., & Cidale, L. 2018, MNRAS, 477, 755

    work page 2018

  3. [3]

    2012, MNRAS, 424, 728

    Barai, P., Proga, D., & Nagamine, K. 2012, MNRAS, 424, 728

    work page 2012

  4. [4]

    P., Waldron, W

    Cassinelli, J. P., Waldron, W. L., Sanders, W. T., et al. 1981, ApJ, 250, 677

    work page 1981

  5. [5]

    I., Abbott, D

    Castor, J. I., Abbott, D. C., & Klein, R. I. 1975, ApJ, 195, 157

    work page 1975

  6. [6]

    & Kunasz, P

    Catala, C. & Kunasz, P. B. 1987, A&A, 174, 158

    work page 1987

  7. [7]

    Cidale, L. S. 1993, PhD thesis, National University of La Plata, Argentina

    work page 1993

  8. [8]

    Cidale, L. S. 1998, ApJ, 502, 824

    work page 1998

  9. [9]

    S., Haucke, M., Arias, M

    Cidale, L. S., Haucke, M., Arias, M. L., et al. 2023, A&A, 677, A176

    work page 2023

  10. [10]

    Cidale, L. S. & Ringuelet, A. E. 1993, ApJ, 411, 874

    work page 1993

  11. [11]

    Clarke, D. A. 1996, ApJ, 457, 291

    work page 1996

  12. [12]

    Clarke, D. A. 2010, ApJS, 187, 119 Curé, M. 2004, ApJ, 614, 929 Curé, M. & Araya, I. 2023, Galaxies, 11, 68 Curé, M., Cidale, L., & Granada, A. 2011, ApJ, 737, 18 de Burgos, A., Keszthelyi, Z., Simón-Díaz, S., & Urbaneja, M. A. 2024, A&A, 687, L16

    work page 2010

  13. [13]

    1995, A&A, 299, 523

    Feldmeier, A. 1995, A&A, 299, 523

    work page 1995

  14. [14]

    Friend, D. B. & Abbott, D. C. 1986, ApJ, 311, 701

    work page 1986

  15. [15]

    L., Mac Low, M.-M., & Abel, T

    Fujita, A., Martin, C. L., Mac Low, M.-M., & Abel, T. 2003, ApJ, 599, 50

    work page 2003

  16. [16]

    Gayley, K. G. 1995, ApJ, 454, 410

    work page 1995

  17. [17]

    S., Venero, R

    Haucke, M., Cidale, L. S., Venero, R. O. J., et al. 2018, A&A, 614, A91

    work page 2018

  18. [18]

    F., Gammie, C

    Hawley, J. F., Gammie, C. F., & Balbus, S. A. 1995, ApJ, 440, 742

    work page 1995

  19. [19]

    & Lanz, T

    Hubeny, I. & Lanz, T. 1995, ApJ, 439, 875

    work page 1995

  20. [20]

    & Henrichs, H

    Kaper, L. & Henrichs, H. F. 1994, Ap&SS, 221, 115

    work page 1994

  21. [21]

    F., Nichols, J

    Kaper, L., Henrichs, H. F., Nichols, J. S., & Telting, J. H. 1999, A&A, 344, 231

    work page 1999

  22. [22]

    Kramida, A., Ralchenko, Y ., Reader, J., & Team, N. A. 2024, NIST Atomic Spec- tra Database (ver. 5.11),https://physics.nist.gov/asd, national Insti- tute of Standards and Technology, Gaithersburg, MD Krtiˇcka, J. & Kubát, J. 2016, Advances in Space Research, 58, 710

    work page 2024

  23. [23]

    Kudritzki, R. P. 2002, ApJ, 577, 389

    work page 2002

  24. [24]

    Kurucz, R. L. 1979, ApJS, 40, 1

    work page 1979

  25. [25]

    Lamers, H. J. G. L. M., Snow, T. P., & Lindholm, D. M. 1995, ApJ, 455, 269

    work page 1995

  26. [26]

    & Puls, J

    Markova, N. & Puls, J. 2008, A&A, 478, 823

    work page 2008

  27. [27]

    & Kunasz, P

    Mihalas, D. & Kunasz, P. B. 1978, ApJ, 219, 635

    work page 1978

  28. [28]

    & Pudritz, R

    Ouyed, R. & Pudritz, R. E. 1997, ApJ, 482, 712

    work page 1997

  29. [29]

    Pauldrach, A., Puls, J., & Kudritzki, R. P. 1986, A&A, 164, 86

    work page 1986

  30. [30]

    2000, A&AS, 141, 23

    Puls, J., Springmann, U., & Lennon, M. 2000, A&AS, 141, 23

    work page 2000

  31. [31]

    2001, PhD thesis, National University of La Plata, Argentina Rubio-Díez, M

    Rohrmann, R. 2001, PhD thesis, National University of La Plata, Argentina Rubio-Díez, M. M., Sundqvist, J. O., Najarro, F., et al. 2022, A&A, 658, A61

    work page 2001

  32. [32]

    M., Mihalas, D., & Norman, M

    Stone, J. M., Mihalas, D., & Norman, M. L. 1992, ApJS, 80, 819

    work page 1992

  33. [33]

    Stone, J. M. & Norman, M. L. 1992, ApJS, 80, 753

    work page 1992

  34. [34]

    2000, in Astronomical Society of the Pacific Conference Series, V ol

    Venero, R., Cidale, L., & Ringuelet, A. 2000, in Astronomical Society of the Pacific Conference Series, V ol. 214, IAU Colloq. 175: The Be Phenomenon in Early-Type Stars, ed. M. A. Smith, H. F. Henrichs, & J. Fabregat, 607

    work page 2000

  35. [35]

    Venero, R. O. J., Curé, M., Cidale, L. S., & Araya, I. 2016, ApJ, 822, 28

    work page 2016

  36. [36]

    Venero, R. O. J., Curé, M., Puls, J., et al. 2024, MNRAS, 527, 93 Article number, page 9 of 10 A&A proofs:manuscript no. aa58896-26 Appendix A: Terminal velocity and mass-loss rates for the T19 model (Ω =0.2) Table A.1: Stellar wind parameters for the solutions found in the fast andδ-slow regimes, and in the gap region for the T19 model withΩ =0.2, comput...

    work page 2024