Wolf-Rayet stellar evolution models with improved treatment of the atmosphere
Pith reviewed 2026-05-20 00:19 UTC · model grok-4.3
The pith
Using detailed CMFGEN atmospheres in evolution models lowers effective temperatures of Wolf-Rayet stars to better match observations while leaving internal structure unchanged.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
By applying outer boundary conditions interpolated from a CMFGEN atmosphere grid at each time step in the STAREVOL code, the models produce effective temperatures that are greatly reduced during the Wolf-Rayet phase and agree better with observations. The internal structure and chemical profiles remain barely affected. The authors show that this direct inclusion of detailed atmospheres is equivalent to applying post-processing corrections to temperature and gravity on standard Eddington-gray models. Accurate placement of evolved massive stars in the Hertzsprung-Russell diagram therefore requires realistic atmospheres, even though internal and chemical evolution are unaffected.
What carries the argument
Interpolation of effective temperature and gravity from a grid of CMFGEN model atmospheres to supply outer boundary conditions at every evolutionary time step.
Load-bearing premise
The outer boundary conditions taken from the CMFGEN grid can be applied at each time step without iterative feedback that would change the interior structure or chemical profiles.
What would settle it
A set of high-precision effective-temperature measurements for Galactic Wolf-Rayet stars that lie systematically closer to the hotter values produced by Eddington-gray models than to the cooler values from the new calculations would falsify the claim of improved observational agreement.
Figures
read the original abstract
Evolutionary models of massive stars are quasi-exclusively computed using an Eddington gray atmosphere. This approximation does not accurately describe the complex physical phenomena occurring in the atmosphere of massive stars. We aim to include state-of-the-art atmosphere models in the evolution computations of massive stars and test how the Wolf-Rayet phase is impacted. We computed the evolution of Galactic massive stars with the code STAREVOL. During the advanced phases of evolution, we applied outer boundary conditions interpolated within a grid of CMFGEN model atmospheres at each time step. The effective temperature and effective gravity were extracted from the atmosphere models. We then compared the resulting evolutionary tracks with classical calculations assuming Eddington gray atmospheres. We find that including detailed model atmospheres has a significant impact on the effective radius and temperature of the models during the later stages of the evolution. The effective temperatures of the evolution models computed with detailed model atmospheres are greatly reduced and in better agreement with observations of Wolf-Rayet stars. On the other hand, the internal structure of the models is barely affected by the choice of the atmosphere. We show that applying post-processing corrections on effective temperature and gravity is a method equivalent to our direct inclusion of atmosphere models in evolutionary calculations. The inclusion of detailed atmosphere models in the computation of evolutionary models is necessary to correctly reproduce the position of evolved massive stars in the Hertzsprung-Russell diagram. However, this has no impact on the internal and chemical evolution.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents evolutionary calculations of Galactic massive stars using the STAREVOL code, replacing the standard Eddington gray outer boundary condition with interpolated values of effective temperature and effective gravity drawn from a grid of CMFGEN detailed model atmospheres during advanced evolutionary phases. The authors report that this change produces substantially lower effective temperatures and larger radii in the Wolf-Rayet phase, improving agreement with observations, while the internal structure, convective zones, and chemical profiles remain essentially unchanged. They further show that the same surface corrections can be obtained by post-processing classical models rather than by direct inclusion during the run.
Significance. If the central result holds, the work supplies a practical route to place evolved massive stars at more observationally consistent locations in the HR diagram without materially altering their core evolution or nucleosynthetic yields. The direct numerical comparison of two boundary-condition choices, performed without additional free parameters, is a clear methodological strength. The finding that internal structure is insensitive to the atmosphere treatment would, if robustly demonstrated, simplify the use of existing grids for supernova-progenitor and chemical-evolution studies.
major comments (2)
- [Method of boundary-condition application] The claim that 'the internal structure of the models is barely affected' (abstract and results) rests on a one-way interpolation of Teff and geff from the static CMFGEN grid at each time step. The manuscript does not test whether the resulting change in radius and luminosity requires an iterative update of the atmosphere model itself; for Wolf-Rayet stars the wind base and optical-depth structure are known to be sensitive to these quantities, so even modest adjustments could feed back into mass-loss rate or mixing and alter the chemical profiles reported as unchanged.
- [Results comparing internal structure] No quantitative measures (differences in core mass, extent of the convective core, or abundance profiles at key epochs) are supplied to support the statement that internal structure is 'barely affected.' Without such metrics or an assessment of numerical noise, it is impossible to judge whether the reported invariance is physically meaningful or merely below the resolution of the comparison.
minor comments (2)
- [Abstract] The abstract states that effective temperatures are 'greatly reduced' but supplies neither the magnitude of the shift nor the range of initial masses for which the effect is largest.
- [Numerical method] Details of the interpolation scheme (grid spacing in Teff, log g, and composition; handling of time-step changes) are not described, which limits reproducibility.
Simulated Author's Rebuttal
We thank the referee for their careful reading and constructive comments on our manuscript. The points raised help clarify our methodology and strengthen the evidence for our conclusions. We address each major comment below and indicate where revisions will be made.
read point-by-point responses
-
Referee: [Method of boundary-condition application] The claim that 'the internal structure of the models is barely affected' (abstract and results) rests on a one-way interpolation of Teff and geff from the static CMFGEN grid at each time step. The manuscript does not test whether the resulting change in radius and luminosity requires an iterative update of the atmosphere model itself; for Wolf-Rayet stars the wind base and optical-depth structure are known to be sensitive to these quantities, so even modest adjustments could feed back into mass-loss rate or mixing and alter the chemical profiles reported as unchanged.
Authors: We appreciate the referee's point on the nature of the coupling. Our method performs a one-way interpolation: at each time step the current luminosity and mass are used to select Teff and geff from the static CMFGEN grid, which then supplies the outer boundary condition for the stellar structure equations. A fully iterative scheme, in which the atmosphere model is recomputed after each radius adjustment, was not implemented because of the substantial additional computational cost. However, the manuscript already demonstrates that applying the same Teff and geff corrections in post-processing to otherwise identical Eddington-gray models yields surface properties and internal profiles that are indistinguishable from the direct-inclusion runs. This equivalence indicates that any potential feedback through mass-loss or mixing remains below the level that affects the reported chemical evolution. In the revised manuscript we will expand the discussion of the one-way coupling, explicitly note its limitations for Wolf-Rayet winds, and cite the post-processing comparison as supporting evidence that the approximation is adequate for the present study. revision: partial
-
Referee: [Results comparing internal structure] No quantitative measures (differences in core mass, extent of the convective core, or abundance profiles at key epochs) are supplied to support the statement that internal structure is 'barely affected.' Without such metrics or an assessment of numerical noise, it is impossible to judge whether the reported invariance is physically meaningful or merely below the resolution of the comparison.
Authors: We agree that quantitative metrics are needed to substantiate the claim. In the revised manuscript we will add a dedicated subsection (or table) that reports the differences, at several well-defined epochs (end of core hydrogen burning, start of the Wolf-Rayet phase, and terminal age), in (i) helium-core mass, (ii) mass coordinate of the outer edge of the convective core, and (iii) surface mass fractions of helium, nitrogen, and carbon. These differences will be presented together with an estimate of numerical noise obtained from convergence tests performed with the STAREVOL code. Preliminary inspection of the existing model outputs shows that the variations remain below 0.5 % in core mass and well within the numerical uncertainty, but the revised text will document this explicitly. revision: yes
Circularity Check
Direct numerical comparison of boundary conditions yields independent results with no definitional or self-referential reduction
full rationale
The paper computes stellar evolution tracks in STAREVOL under two explicit outer-boundary prescriptions (interpolated CMFGEN grid versus Eddington gray) and reports the resulting differences in Teff, radius, and internal profiles as direct numerical outcomes. No parameter is fitted to the target WR observations, no prediction is constructed from a subset of the same data, and no uniqueness theorem or ansatz is imported via self-citation to force the conclusion. The statement that internal structure is barely affected is an empirical finding from the runs, not an equivalence by construction. The post-processing equivalence is likewise shown by explicit comparison rather than assumed. The derivation chain is therefore self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
axioms (2)
- standard math Standard equations of stellar structure and evolution hold when outer boundary conditions are updated from detailed atmosphere models.
- domain assumption The CMFGEN grid accurately represents the atmospheres of Galactic Wolf-Rayet stars at the metallicities and luminosities reached by the models.
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We applied outer boundary conditions interpolated within a grid of CMFGEN model atmospheres at each time step. The effective temperatures ... are greatly reduced ... internal structure of the models is barely affected
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
-
[1]
Agrawal, P., Stevenson, S., Szécsi, D., & Hurley, J. 2022, A&A, 668, A90
work page 2022
-
[2]
Alongi, M., Bertelli, G., Bressan, A., et al. 1993, A&AS, 97, 851
work page 1993
-
[3]
Amard, L., Palacios, A., Charbonnel, C., et al. 2019, A&A, 631, A77
work page 2019
-
[4]
Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 Böhm-Vitense, E. 1958, Z. Astrophys., 46, 108
work page 2009
- [5]
- [6]
- [7]
-
[8]
H., Watson, L., Lux, T., et al
Chang, T. H., Watson, L., Lux, T., et al. 2020, ACM Trans. Math. Softw., 46, 1
work page 2020
-
[9]
Chang, T. H., Watson, L. T., Lux, T. C. H., et al. 2018, in Proceedings of the 2018 ACM Southeast Conference, ACMSE ’18 (New York, NY , USA: Association for Computing Machinery)
work page 2018
-
[10]
Cheng, S. J., Goldberg, J. A., Cantiello, M., et al. 2024, ApJ, 974, 270
work page 2024
- [11]
-
[12]
Crowther, P. A. 2007, ARA&A, 45, 177
work page 2007
-
[13]
Crowther, P. A., Dessart, L., Hillier, D. J., Abbott, J. B., & Fullerton, A. W. 2002, A&A, 392, 653 de Jager, C., Nieuwenhuijzen, H., & van der Hucht, K. A. 1988, A&AS, 72, 259
work page 2002
-
[14]
Debnath, D., Sundqvist, J. O., Moens, N., et al. 2024, A&A, 684, A177 Degl’Innocenti, S., Prada Moroni, P. G., Marconi, M., & Ruoppo, A. 2008, Ap&SS, 316, 25 Ekström, S., Georgy, C., Eggenberger, P., et al. 2012, A&A, 537, A146
work page 2024
-
[15]
Eldridge, J. J. & Stanway, E. R. 2022, ARA&A, 60, 455
work page 2022
-
[16]
Ferguson, J. W., Alexander, D. R., Allard, F., et al. 2005, ApJ, 623, 585
work page 2005
- [17]
-
[18]
Garcia, M., Herrero, A., Najarro, F., Lennon, D. J., & Alejandro Urbaneja, M. 2014, ApJ, 788, 64 Gräfener, G. & Hamann, W. R. 2008, A&A, 482, 945 Gräfener, G., Koesterke, L., & Hamann, W. R. 2002, A&A, 387, 244 Gräfener, G., Owocki, S. P., & Vink, J. S. 2012, A&A, 538, A40
work page 2014
-
[19]
H., Meynet, G., Ekström, S., & Georgy, C
Groh, J. H., Meynet, G., Ekström, S., & Georgy, C. 2014, A&A, 564
work page 2014
- [20]
-
[21]
R., Gräfener, G., & Liermann, A
Hamann, W. R., Gräfener, G., & Liermann, A. 2006, A&A, 457, 1015
work page 2006
-
[22]
R., Gräfener, G., Liermann, A., et al
Hamann, W. R., Gräfener, G., Liermann, A., et al. 2019, A&A, 625, A57
work page 2019
-
[23]
Henyey, L., Vardya, M. S., & Bodenheimer, P. 1965, ApJ, 142, 841
work page 1965
- [24]
-
[25]
Hillier, D. J., Lanz, T., Heap, S. R., et al. 2003, ApJ, 588, 1039
work page 2003
-
[26]
Hillier, D. J. & Miller, D. L. 1998, ApJ, 496, 407
work page 1998
-
[27]
Hubeny, I. & Mihalas, D. 2014, Theory of Stellar Atmospheres (Princeton Uni- versity Press)
work page 2014
-
[28]
Iglesias, C. A. & Rogers, F. J. 1996, ApJ, 464, 943
work page 1996
- [29]
-
[30]
Jiang, Y .-F., Cantiello, M., Bildsten, L., Quataert, E., & Blaes, O. 2015, ApJ, 813, 74
work page 2015
-
[31]
Johnson, J. A. 2019, Sci, 363, 474
work page 2019
-
[32]
Josiek, J., Sander, A. A. C., Bernini-Peron, M., et al. 2025, A&A, 697, A193
work page 2025
-
[33]
1990, Stellar Structure and Evolution, second edition (Springer) Krishna Swamy, K
Kippenhahn, R., Weigert, A., & Weiss, A. 1990, Stellar Structure and Evolution, second edition (Springer) Krishna Swamy, K. 1966, ApJ, 145, 174
work page 1990
-
[34]
Lefever, R. R., Sander, A. A. C., Bernini-Peron, M., et al. 2026, A&A, 707, A8
work page 2026
- [35]
-
[36]
Manchon, L., Deal, M., Marques, J. P. C., & Lebreton, Y . 2025, A&A, 704, A79
work page 2025
- [37]
- [38]
- [39]
-
[40]
Martins, F., Depagne, E., Russeil, D., & Mahy, L. 2013, A&A, 554, A23
work page 2013
- [41]
-
[42]
Martins, F., Schaerer, D., Hillier, D. J., & Heydari-Malayeri, M. 2004, A&A, 420, 1087
work page 2004
- [43]
-
[44]
Moens, N., Sundqvist, J. O., El Mellah, I., et al. 2022, A&A, 657, A81 Montalbán, J., D’Antona, F., Kupka, F., & Heiter, U. 2004, A&A, 416, 1081
work page 2022
- [45]
- [46]
- [47]
-
[48]
Pauli, D., Langer, N., Schootemeijer, A., et al. 2026, A&A, 707, A11
work page 2026
- [49]
- [50]
-
[51]
Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992, Nu- merical recipes in FORTRAN. The art of scientific computing (Cambridge University Press)
work page 1992
-
[52]
Sana, H., de Mink, S. E., de Koter, A., et al. 2012, Sci, 337, 444
work page 2012
-
[53]
Sander, A. A. C., Hamann, W.-R., Todt, H., et al. 2019, A&A, 621, A92
work page 2019
-
[54]
Sander, A. A. C., Lefever, R. R., Josiek, J., et al. 2026, Nat. Astron, 10, 290
work page 2026
-
[55]
Sander, A. A. C., Lefever, R. R., Poniatowski, L. G., et al. 2023, A&A, 670, A83
work page 2023
-
[56]
Sander, A. A. C. & Vink, J. S. 2020, MNRAS, 499, 873
work page 2020
-
[57]
Santolaya-Rey, A. E., Puls, J., & Herrero, A. 1997, A&A, 323, 488
work page 1997
-
[58]
Sanyal, D., Grassitelli, L., Langer, N., & Bestenlehner, J. M. 2015, A&A, 580, A20
work page 2015
-
[59]
Schaerer, D., de Koter, A., Schmutz, W., & Maeder, A. 1996, A&A, 310, 837
work page 1996
-
[60]
Schaller, G., Schaerer, D., Meynet, G., & Maeder, A. 1992, A&AS, 96, 269
work page 1992
- [61]
- [62]
-
[63]
Vink, J. S., de Koter, A., & Lamers, H. J. G. L. M. 2001, A&A, 369, 574 Article number, page 12 Thomas V oje et al.: Wolf-Rayet stellar evolution models with improved treatment of the atmosphere Appendix A: Grid of atmosphere models for the 50M ⊙ track Figure A.1 shows the grid of CMFGEN atmosphere models we used to compute the 50M ⊙ track with detailed a...
work page 2001
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.