REVIEW 3 major objections 7 minor 48 references
Hydrodynamical mass-loss rates for very massive stars II. New theoretical mass-loss predictions at solar metallicity (Z = 0.02)
T0 review · 3 major / 7 minor · reviewed 2026-07-31 · grok-4.5
Pith's one-line read A single continuous mass-loss recipe captures the kink from optically thin O-star winds to optically thick WNh winds at solar metallicity.
desk verdict Solid, usable solar-Z mass-loss recipe from a large PoWR-HD grid; kink and bistability confirmed across L/X, anchored at Arches and Pauli, with the main systematic being prescribed clumping/turbulence. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
A log-sum-exponential bridge between low- and high-Eddington power laws (the paper’s Eq. 5), fed by hydrodynamically consistent PoWR-HD solutions that predict mass-loss rate and velocity structure instead of assuming a prescribed beta-law.
What would settle it
Obtain independent mass-loss rates and Eddington parameters for Of/WNh transition stars in other young massive clusters; if the rates at the spectral thin-to-thick transition systematically miss log(Mdot) ≈ −5.15 at the matching luminosity, the absolute calibration fails.
Extended reading notes
Core claim
Across 178 PoWR-HD models at Z=0.02, mass loss versus classical Eddington parameter shows a kink: a shallow scaling near 2.8 for optically thin O-star winds becomes a steep scaling near 10 once winds are optically thick. The kink occurs where the flux-weighted optical depth at the sonic point is order unity and the wind efficiency is about 0.4. One continuous fit captures that kink, two iron ionisation bistability dips, and explicit luminosity and hydrogen-abundance scalings, and it recovers the Arches transition rate of about log(Mdot) = −5.2.
Load-bearing premise
Clumping and turbulent velocity are fixed by hand in every model rather than predicted, so the entire fitted recipe inherits those prescribed choices.
Editorial extensions
If this is right
- Stellar evolution codes can replace separate O-star and VMS wind recipes with one continuous function over ~20–500 solar masses at Z=0.02.
- Above the kink, steeper high-Eddington mass loss will change retained mass, black-hole masses, and pair-instability outcomes relative to older shallow recipes.
- Below the transition the new rates lie below classic Vink et al. (2001) rates, reducing main-sequence mass loss for ordinary O stars.
- Auxiliary relations convert surface Teff and high-T Eddington estimates into the inner-boundary quantities the fit needs, enabling direct use in structure codes.
- The same grid supplies terminal velocities and ionising photon rates as secondary outputs across the parameter space.
Reading between the lines
- Because the kink tracks sonic-point optical depth of order unity, lower-metallicity hydro grids should shift the kink in Eddington space and change the initial-mass threshold for WNh-type winds.
- Once depth-dependent turbulence from multi-D envelope simulations is solved inside the same hydro code, absolute rates may move by tenths of a dex, especially at cool temperatures.
- Agreement on the ZAMS with empirical Mdot–Γe slopes implies that, at solar metallicity, the dominant remaining uncertainty for very massive star endpoints is less the main-sequence wind law than binary mass transfer and post-main-sequence mass loss.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The authors present a continuous theoretical mass-loss recipe for massive and very massive stars at Z=0.02, derived from a grid of 178 hydrodynamically consistent PoWR-HD non-LTE wind models spanning M⋆≈16–500 M⊙, log(L⋆/L⊙)=5.5–6.8, T⋆=12–50 kK, and X=0.01–0.9. They confirm a kink in the Ṁ–Γe relation (shallow slope ∼2.8 below the kink; steep ∼10 above), linked to τF,sonic∼1 and η∼0.4, and two Fe-driven bistability features near T⋆≈25 and 17 kK. These behaviours, plus explicit L⋆ and X scalings, are captured in a single fitting formula (Eq. 5) with auxiliary T⋆–Teff(τR=2/3) and incomplete-ionisation Γe corrections for evolutionary implementation. Absolute scale is tested against the model-independent Arches O-to-WNh transition mass-loss rate (predicted log Ṁ≈−5.15 vs ≈−5.2) and against the Pauli et al. (2025) empirical Ṁ–Γe trend on the ZAMS.
Significance. Mass loss is the dominant uncertainty for VMS evolution, pair-instability boundaries, and heavy black-hole formation. A hydrodynamically self-consistent, continuous recipe that joins optically thin O-star and optically thick WNh regimes—and is anchored to the Arches transition point rather than only to internal model normalisations—is of clear practical value for stellar evolution. Strengths include the large documented grid, explicit coefficients with uncertainties and quoted RMSE (0.12 dex on Ṁ), auxiliary implementation relations, and two external anchors (Arches transition; Pauli et al. 2025) not used to set the hydro solutions. Terminal-velocity and ionising-flux predictions add further utility. The work is a natural and useful extension of Paper I and of the η-framework of Sabhahit et al. (2022).
major comments (3)
- [§2, Appendix A, Eq. (5)] §2 and Appendix A: Every hydrodynamic solution (hence every coefficient in Eq. 5) is computed with prescribed micro-clumping (Dcl=1→10, onset τcl=0.1) and a fixed radially constant vturb=70.71 km s−1. Appendix A already shows that vturb variations can shift Ṁ at the ∼0.2 dex level (especially at cool T⋆). The Arches τF,sonic=1 anchor (Fig. 8) remains persuasive and is partly robust to vturb, but the manuscript does not quantify how changing the clumping law (factor or onset) would move the absolute Ṁ(Γe) curve or the high-Γe slope. For a recipe intended for evolution codes, please add a short, explicit systematic-uncertainty floor (or a limited clumping sensitivity test) so users know the precision claimed beyond the 0.12 dex fit RMSE.
- [§3.1, §3.4, Fig. 6] §3.3–3.4 and Fig. 6: The models systematically under-predict terminal velocities relative to Arches transition objects (η≈0.4 at the kink vs the η≈0.6 used in the Vink & Gräfener 2012 argument). The paper correctly notes that the recipe is parametrized in Γe rather than η, so Ṁ at the kink is still well matched. However, evolutionary applications often need both Ṁ and v∞ (wind momentum, mechanical feedback). Please state more clearly whether v∞ from the grid (or a recommended scaling) should be used with Eq. 5, and flag the v∞ under-prediction as a separate limitation rather than only as an explanation for low η.
- [§3.3, Table 2, §4.3] §3.3 and §4.3: Implementation requires iterating Eqs. (5), (6), and (9) from structure-code Teff(τR=2/3) and Γe,high-T. The text describes the procedure but does not demonstrate numerical stability or provide a minimal worked example (e.g., one ZAMS row from Table 2 with intermediate T⋆, Γe, and Ṁ iterates). Given the steep high-Γe branch, a 10% Γe error is stated to cost ∼0.5 dex in Ṁ; a short convergence note or pseudocode would make the “suitable for stellar evolution” claim load-bearing rather than aspirational.
minor comments (7)
- [Abstract, §5] Abstract and Conclusions quote suitability for 20–500 M⊙ (abstract) and 30–500 M⊙ (conclusions). Align the stated mass range.
- [Fig. 1] Fig. 1 caption refers to “best-fit relations described in Sect. 3.3” before that section; a forward reference is fine but ensure all solid-line sequences are uniquely identifiable against Table 1 sequence indices (as in Fig. 2).
- [§3.3, Eq. (5)] Eq. (5): T2 is written as depending on Γe/Γe,ref while T1 is fixed; a one-sentence physical motivation for allowing only the cooler dip centre to shift would help readers.
- [Fig. 9, §4.2] Fig. 9 vs Fig. 1: the apparent high-Γe slope difference at fixed Teff versus fixed T⋆ is important for users; consider a brief callout in the Fig. 9 caption pointing to the T⋆–Teff correction in §3.3.
- [Appendix C] Table C.1 is very long; if journal length is an issue, consider moving the bulk to CDS/Zenodo and retaining a representative subset plus a machine-readable full table statement in the text.
- [Throughout] Minor typography: “3 turb”, “3∞”, “PoWR hd/PoWRhd” appear in several spellings; unify notation for vturb, v∞, and PoWR-HD across text and tables.
- [§3.2, §4.2] References to Vink et al. (2026, submitted) and related in-press items should be updated or marked consistently at proof stage.
Circularity Check
No significant circularity: the recipe is an explicit fit to independent hydro models, externally checked against Arches and Pauli et al.
full rationale
The central product is a multi-parameter fit (Eq. 5) to 178 PoWR-HD hydrodynamic solutions in which Mdot and the velocity law are outputs of the equation of motion, not prescribed inputs. That fit describing its own training grid is ordinary recipe construction, not a circular derivation. Absolute scale is not forced by the Arches transition criterion: the models independently place τF,sonic≈1 and η≈0.4 at the kink and yield log Mdot(τF,sonic=1)≈−5.15 in the Arches L–T window (Fig. 8), which is then compared to the externally estimated transition rate ≈−5.2; the paper even reports a mild η mismatch (∼0.4 vs ∼0.6) rather than tuning to the Vink & Gräfener factor. ZAMS application uses MESA structure L and Teff and is compared to the independent Pauli et al. (2025) empirical Mdot–Γe relation without refitting. Clumping and vturb are fixed priors that set a systematic floor (App. A), which is a modeling assumption, not a definitional loop. Overlap with Vink & Gräfener (2012) and Paper I is normal sequential self-citation and does not make the hydro solutions or the external numerical anchors reduce to inputs by construction.
Assumptions & free parameters
free parameters (4)
- log Mdot_0 and kink/bistability coefficients in Eq. 5 =
log Mdot_0=-5.448±0.021; fhigh slope=9.243±0.277; Γe,ref base 0.43; etc.
- Clumping factor Dcl and onset τcl =
Dcl,max=10, τcl=0.1
- Turbulent velocity vturb =
70.71 km/s (main); α_min=0.062, α_diff=0.713, T0=47.8 kK
- T*–Teff and Γe ionization correction coefficients =
e.g. 0.542±0.012 in Γe correction; a,b coefficients in Eq. 6
assumptions (6)
- domain assumption Stationary, spherically symmetric, non-LTE CMF radiative transfer with hydro equation of motion solved self-consistently (PoWR-HD framework).
- domain assumption Micro-clumping formalism adequately represents wind inhomogeneity for mass-loss prediction.
- domain assumption Iron-group line driving at Z=0.02 with Grevesse & Sauval (1998) scaled abundances dominates the wind force; recent solar O revision is negligible.
- domain assumption Classical Eddington parameter Γe (electron-scattering) is the primary organizing variable for Mdot, with secondary explicit L and X terms.
- domain assumption Transition mass-loss diagnostic η≈0.6 τF,sonic with τF,sonic~1 gives a nearly model-independent Mdot at the O/WNh boundary (Vink & Gräfener 2012).
- ad hoc to paper Log-sum-exponential connection plus two Gaussian-like bistability dips in T* is a sufficient phenomenological form for the grid.
Cite this review
Pith. "Pith review of Hydrodynamical mass-loss rates for very massive stars II. New theoretical mass-loss predictions at solar metallicity (Z = 0.02)." pith.science (2026). https://pith.science/paper/Y2LMQKJ7
@misc{pith2026260728012,
author = {Pith},
title = {Pith review of: Hydrodynamical mass-loss rates for very massive stars II. New theoretical mass-loss predictions at solar metallicity (Z = 0.02)},
year = {2026},
howpublished = {\url{https://pith.science/paper/Y2LMQKJ7}},
note = {Machine review of arXiv:2607.28012}
}
abstract
The evolutionary pathways and ultimate fates of very massive stars are governed primarily by mass loss through radiatively-driven winds. We present a new theoretical mass-loss prescription for (very) massive stars, capturing the complex dependence on the Eddington parameter $\Gamma_e$, luminosity, temperature, and hydrogen abundance. We calculated an extensive grid of 178 hydrodynamically consistent wind-atmosphere models in non-local thermodynamic equilibrium using the PoWR-HD code, predicting wind properties such as the mass-loss rate and terminal velocity self-consistently. The grid spans masses $M_*$ = 16-500 Msun, luminosities $\log(L_*/L_\odot) = 5.5-6.8$, inner boundary temperatures $T_* = 12-50$ kK, and hydrogen mass fractions X = 0.01-0.9, at a fixed metallicity of Z=0.02. We confirm the presence of a mass-loss kink in the $\dot{M}-\Gamma_e$ relation across the explored parameter space. The kink marks the transition from a shallow scaling ($\sim 2.8$) at low $\Gamma_\mathrm{e}$ for optically thin O-star winds to a steeper scaling ($\sim 10$) for optically thick winds at high $\Gamma_e$. We derive comprehensive fitting relations capturing both the kink behaviour and two bistability jumps arising from iron ionisation changes, and provide auxiliary relations for implementation into stellar evolutionary calculations. Our prescription correctly reproduces the model-independent transition mass-loss rate in the Arches Cluster, confirming the accuracy of our predicted rates at the O-to-WNh transition. Application of our recipe to the Zero Age Main Sequence provides excellent agreement with recent empirical $\dot{M}-\Gamma_e$ relation obtained for a wide range of temperatures and Eddington parameters. We provide a physically motivated, continuous, and empirically anchored mass-loss recipe for (very) massive stars, suitable for stellar evolution calculations in the 20-500 Msun range.
Figures
Figures from the paper (6 more)
Reference graph
Works this paper leans on
-
[1]
J., & Scott, P
Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481
2009
-
[2]
Bernini-Peron, M., Sander, A. A. C., Sabhahit, G. N., et al. 2026, A&A, 708, A206
2026
-
[3]
M., Gräfener, G., Vink, J
Bestenlehner, J. M., Gräfener, G., Vink, J. S., et al. 2014, A&A, 570, A38 Björklund, R., Sundqvist, J. O., Puls, J., & Najarro, F. 2021, A&A, 648, A36 Björklund, R., Sundqvist, J. O., Singh, S. M., Puls, J., & Najarro, F. 2023, A&A, 676, A109
2014
-
[4]
A., de Koter, A., Bestenlehner, J
Brands, S. A., de Koter, A., Bestenlehner, J. M., et al. 2022, A&A, 663, A36
2022
-
[5]
I., Abbott, D
Castor, J. I., Abbott, D. C., & Klein, R. I. 1975, ApJ, 195, 157
1975
-
[6]
2015, MNRAS, 452, 1068
Chen, Y ., Bressan, A., Girardi, L., et al. 2015, MNRAS, 452, 1068
2015
-
[7]
A., Schnurr, O., Hirschi, R., et al
Crowther, P. A., Schnurr, O., Hirschi, R., et al. 2010, MNRAS, 408, 731
2010
-
[8]
Crowther, P. A. & Walborn, N. R. 2011, MNRAS, 416, 1311
2011
Show all 48 references
-
[9]
O., Moens, N., et al
Debnath, D., Sundqvist, J. O., Moens, N., et al. 2024, A&A, 684, A177
2024
-
[10]
E., Marchant, P., & Justham, S
Farmer, R., Renzo, M., de Mink, S. E., Marchant, P., & Justham, S. 2019, ApJ, 887, 53 Gräfener, G. 2021, A&A, 647, A13 Gräfener, G. & Hamann, W. R. 2005, A&A, 432, 633 Gräfener, G. & Hamann, W. R. 2008, A&A, 482, 945 Gräfener, G., Koesterke, L., & Hamann, W. R. 2002, A&A, 387, 244
2019
-
[11]
& Sauval, A
Grevesse, N. & Sauval, A. J. 1998, Space Sci. Rev., 85, 161
1998
-
[12]
Hamann, W. R. & Gräfener, G. 2003, A&A, 410, 993
2003
-
[13]
Hamann, W. R. & Koesterke, L. 1998, A&A, 335, 1003
1998
-
[14]
R., Vink, J
Higgins, E. R., Vink, J. S., Hirschi, R., Laird, A. M., & Sabhahit, G. N. 2023, MNRAS, 526, 534 Krtiˇcka, J., Kubát, J., & Krtiˇcková, I. 2021, A&A, 647, A28
2023
-
[15]
Krumholz, M. R. & Matzner, C. D. 2009, ApJ, 703, 1352
2009
-
[16]
Lamers, H. J. G. L. M., Snow, T. P., & Lindholm, D. M. 1995, ApJ, 455, 269
1995
-
[17]
2012, ARA&A, 50, 107
Langer, N. 2012, ARA&A, 50, 107
2012
-
[18]
R., Sander, A
Lefever, R. R., Sander, A. A. C., Bernini-Peron, M., et al. 2025, A&A, 700, A2
2025
-
[19]
J., Paumard, T., et al
Martins, F., Hillier, D. J., Paumard, T., et al. 2008, A&A, 478, 219
2008
-
[20]
2025, A&A, 704, A121 Müller-Horn, J., Ramachandran, V ., El-Badry, K., et al
Moens, N., Debnath, D., Verhamme, O., et al. 2025, A&A, 704, A121 Müller-Horn, J., Ramachandran, V ., El-Badry, K., et al. 2026, A&A, 708, A187
2025
-
[21]
M., Hamann, W.-R., et al
Pauli, D., Oskinova, L. M., Hamann, W.-R., et al. 2025, A&A, 697, A114
2025
-
[22]
S., & Gräfener, G
Petrov, B., Vink, J. S., & Gräfener, G. 2016, MNRAS, 458, 1999
2016
-
[23]
S., & Najarro, F
Puls, J., Vink, J. S., & Najarro, F. 2008, The Astronomy and Astrophysics Re- view, 16, 209
2008
-
[24]
Ramachandran, V ., Sander, A. A. C., Pauli, D., et al. 2024, arXiv e-prints, arXiv:2406.17678
2024 arXiv
-
[25]
S., Krumholz, M
Roy, A., Sutherland, R. S., Krumholz, M. R., Heger, A., & Dopita, M. A. 2020, MNRAS, 494, 3861
2020
-
[26]
N., Vink, J
Sabhahit, G. N., Vink, J. S., Higgins, E. R., & Sander, A. A. C. 2022, MNRAS, 514, 3736
2022
-
[27]
N., Vink, J
Sabhahit, G. N., Vink, J. S., & Sander, A. A. C. 2026, A&A, 706, A97
2026
-
[28]
N., Vink, J
Sabhahit, G. N., Vink, J. S., Sander, A. A. C., et al. 2025, A&A, 696, A200
2025
-
[29]
N., Vink, J
Sabhahit, G. N., Vink, J. S., Sander, A. A. C., & Higgins, E. R. 2023, MNRAS, 524, 1529
2023
-
[30]
2015, A&A, 577, A13
Sander, A., Shenar, T., Hainich, R., et al. 2015, A&A, 577, A13
2015
-
[31]
Sander, A. A. C. 2015, PhD thesis, University of Potsdam, Germany
2015
-
[32]
Sander, A. A. C., Hamann, W. R., Todt, H., Hainich, R., & Shenar, T. 2017, A&A, 603, A86
2017
-
[33]
Sander, A. A. C. & Vink, J. S. 2020, MNRAS, 499, 873
2020
-
[34]
Schneider, F. R. N., Ramírez-Agudelo, O. H., Tramper, F., et al. 2018, A&A, 618, A73
2018
-
[35]
S., & de Koter, A
Smith, N., Vink, J. S., & de Koter, A. 2004, ApJ, 615, 475
2004
-
[36]
2024, A&A, 692, A91
Verhamme, O., Sundqvist, J., de Koter, A., et al. 2024, A&A, 692, A91
2024
-
[37]
Vink, J. S. 2018, A&A, 615, A119
2018
-
[38]
Vink, J. S. 2022, ARA&A, 60, 203
2022
-
[39]
S., de Koter, A., & Lamers, H
Vink, J. S., de Koter, A., & Lamers, H. J. G. L. M. 1999, A&A, 350, 181
1999
-
[40]
S., de Koter, A., & Lamers, H
Vink, J. S., de Koter, A., & Lamers, H. J. G. L. M. 2001, A&A, 369, 574
2001
-
[41]
Vink, J. S. & Gräfener, G. 2012, ApJ, 751, L34
2012
-
[42]
S., Heger, A., Krumholz, M
Vink, J. S., Heger, A., Krumholz, M. R., et al. 2015, Highlights of Astronomy, 16, 51
2015
-
[43]
S., Muijres, L
Vink, J. S., Muijres, L. E., Anthonisse, B., et al. 2011, A&A, 531, A132
2011
-
[44]
S., Sabhahit, G
Vink, J. S., Sabhahit, G. N., & Sander, A. A. C. 2026, arXiv e-prints, arXiv:2606.17130
2026 arXiv
-
[45]
Vink, J. S. & Sander, A. A. C. 2021, MNRAS, 504, 2051
2021
-
[46]
Winch, E. R. J., Vink, J. S., Higgins, E. R., & Sabhahitf, G. N. 2024, MNRAS, 529, 2980
2024
-
[47]
E., Blinnikov, S., & Heger, A
Woosley, S. E., Blinnikov, S., & Heger, A. 2007, Nature, 450, 390
2007
-
[48]
2013, MNRAS, 433, 1114 Article number, page 12 of 18 Gautham N
Yusof, N., Hirschi, R., Meynet, G., et al. 2013, MNRAS, 433, 1114 Article number, page 12 of 18 Gautham N. Sabhahit,⋆ et al.: Hydrodynamical mass-loss rates for very massive stars Appendix A: Turbulent velocity correction for mass-loss rate The mass-loss recipe in Eq. (5) is b...
2013
Reviewed July 31, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.