Slow Massive-Star Formation and the O-Star ZAMS Deficit

arXiv:2610.31482 · astro-ph.SR, astro-ph.GA · Submitted 2026-10-22 · Read on arXiv

Listen

Transcript

Introduction to the show: ident: Astrophysics Radio. Generated commentary on the latest astrophysics papers.

Vera: I'm Vera.

Jocelyn: And I'm Jocelyn, here with Subrahmanyan.

Vera: We're looking at the paper: "Slow Massive-Star Formation and the O-Star ZAMS Deficit".

Jocelyn: Here's the gist of the paper. The gist The scarcity of optically visible O stars near the theoretical zero-age main sequence (ZAMS) reflects prolonged accretion, suggesting that massive stars assemble over Myr timescales and their birthline is displaced from the standard ZAMS The Problem and Hypothesis Optically visible O stars are scarce near the theoretical zero-age main sequence (ZAMS), a deficit that has puzzled astronomers for over four decades Holgado et al.

Vera: (2020) found a pronounced deficit near the ZAMS for initial masses of approximately 30–70 M⊙, while Schootemeijer et al.

Jocelyn: (2021) found that Small Magellanic Cloud stars above ≃ 20 M⊙ are similarly rare during the first half of core-H burning A complementary discrepancy appears during the embedded phase, where compact H ii regions are far more abundant relative to optically visible O stars than expected if they are freely expanding nebulae around stars that have completed their formation Both discrepancies may measure the duration of massive-star growth The paper tests whether the Galactic O-star ZAMS deficit can provide a second, independent clock for this growth The Inertial-Inflow Model and Growth Law The inertial-inflow model (IIM) of Padoan et al.

Vera: (2020) provides a physically different picture, where massive-star growth is controlled by the turnover time of the turbulent flow on the scale from which the mass is collected The IIM predicts that formation time increases approximately as tform ∝ m f(1/2) This growth law is tested using accreting stellar-evolution models computed with MESA The model begins with the same pre-main-sequence seed model of mass mseed = 2 M⊙ and grows at the constant rate m acc (m f) = m f − mseed / tform(m f) Constraints from Observational Data The analysis compares the end-of-accretion points of accreting models with the full observed envelope of Galactic O stars in the Kiel diagram The hot-side boundary of this distribution is used to constrain growth-law parameters (alpha, tau0) The mismatch between the birthline and the fixed observational envelope is quantified by Q(alpha, tau0) Joint Constraint and Conclusion The Kiel-envelope analysis favors (alpha, tau0) = (0.52, 4.3 Myr), while the luminosity-function analysis gives posterior median values of (0.48, 4.15 Myr) The two measurements are mutually consistent because they probe the same underlying Galactic population and constrain the same parameters through different observables The combination using Fisher’s method yields projected 90% intervals 0.43 < alpha < 0.60 and 3.9 < tau0/Myr < 5.0, which is much tighter than either measurement alone The jointly preferred parameters define the effective Galactic growth law tform ≃ 4.3 Myr(m f/60 M⊙) 0.5 This work blurs the conventional boundary between massive-star formation and stellar evolution, showing that massive stars can spend a substantial fraction of their core-hydrogen-burning lifetime still assembling their mass How it works The analysis constructs birthlines by defining the formation time as a function of m f using the growth law tform(m f) = τ0 (m f/60 M⊙) α Non-rotating, solar-metallicity accreting models are computed with MESA, employing specific parameters for overshooting and mixing The resulting 24 distinct growth laws are sampled by computing 17 or 18 accreting tracks for different values of m f Key Findings The paper demonstrates that the accreting models naturally reproduce both the displacement and turnover of the hot-side envelope of Galactic O stars in the Kiel diagram The resulting Kiel-envelope constraint is shown in Fig.

Paper summary: Vera: So, to wrap up the "Slow Massive-Star Formation and the O-Star ZAMS Deficit" paper, what does this mean for us in terms of understanding massive star birth? It suggests that these stars are still actively growing when we observe them near their zero-age main sequence.

Jocelyn: It really shifts the view from instantaneous collapse to a more gradual process where accretion lasts for millions of years, leading to a displaced birthline rather than an exact match with the ZAMS.

Subrahmanyan: The title itself points to this idea—that massive star formation is slow enough that we see the effects of prolonged accretion on their observable properties. This supports models like the inertial-inflow model which predict mass supply from turbulent flows over time, rather than just a single rapid collapse event.

Vera: It means when we study these young populations, we can’t just assume they are born exactly where the ZAMS is; we have to account for this growth phase quantified by those tau zero and alpha parameters <ref:2610.31482#pg1>.

Jocelyn: And from a radio astronomy perspective, if formation is slow, it affects when and how these stars interact with their environment as they evolve. It gives us a longer window to study them before they settle onto the main sequence.

Subrahmanyan: Ultimately, this work suggests that incorporating stellar growth into models that currently start from fully formed massive stars will have consequences well beyond just solving the ZAMS deficit; it affects how we understand the entire evolutionary path of massive stars.

Conclusion: Vera: So we’re wrapping up on "Slow Massive-Star Formation and the O-Star ZAMS Deficit." Basically, this paper shows that the lack of bright, visible O stars near where they should be born is because they take a long time to assemble their mass.

Jocelyn: Yeah, so instead of thinking massive stars are born instantly right on the zero-age main sequence line, it suggests prolonged accretion over millions of years. That means their birthline is shifted away from that standard ZAMS point.

Subrahmanyan: If we follow the inertial-inflow model they tested, this implies massive star assembly happens over Myr timescales. Core hydrogen burning starts before all the gas is actually accreted onto the star's core, and then it becomes visible later when accretion finally stops.

Vera: The authors used data from Galactic O stars in a Kiel diagram to test this idea against accreting stellar models. They found that their specific growth law parameters fit the observed envelope displacement and the high-mass turnover well.

Jocelyn: And they combined that with constraints from luminosity functions of compact H II regions, which gave them a tighter range for those growth law parameters. The final results point to a formation time roughly four point three Myr for stars around sixty times the mass of the sun.

Subrahmanyan: That means we have to rethink how we model the entire life cycle of massive stars because they spend a significant chunk of their core hydrogen burning lifetime still adding mass.

Vera: It really blurs that line between formation and evolution. The next thing to consider is how this growth affects the feedback mechanisms those stars have on the surrounding galaxy.

Paolo Padoan, Samuel Shipkoff, Mark Gieles, Fabrice Martins

Department of Physics and Astronomy, Dartmouth College · Institut de Ci¢ns del Cosmos, Universitat de Barcelona · Waterville Valley Academy · ICREA Institute of Science and Technology for Research in Astronomy and Astrophysics · Institut d’Estudis Espacials de Catalunya

astro-ph.SR, astro-ph.GA

Submitted: 2026-10-22

Updated: 2026-10-22

Comments: Submitted to MNRAS Letters

License: http://creativecommons.org/licenses/by/4.0/

The gist: The gist The scarcity of optically visible O stars near the theoretical zero-age main sequence (ZAMS) reflects prolonged accretion, suggesting that massive stars assemble over Myr timescales and

Key concepts

Zero-Age Main Sequence (ZAMS): The theoretical point on a stellar diagram where a star officially begins its core hydrogen burning phase after all gravitational collapse and accretion is complete. Astronomers observe a deficit of stars near this point, which the paper suggests is due to ongoing mass accumulation.
Inertial-Inflow Model (IIM): A physical model suggesting that massive star growth is governed by the turnover time of turbulent flows as material falls onto the star. This model predicts that a star's formation time increases with its final mass, providing a mechanism for prolonged accretion.
Growth Law ($t_{\mathrm{form}}$): A mathematical relationship describing how long it takes for a massive star to reach its final mass ($m_f$). The paper tests the specific law $t_{\mathrm{form}}=τ_0(m_{\mathrm{f}}/60 M_\odot)^α$, using observational data to determine the constants $\tau_0$ and $\alpha$ that best fit the observed stellar populations.
Birthline Displacement: The actual location of a star on a diagram when it is finally considered 'formed' based on accretion models, rather than its theoretical ZAMS position. The paper shows that this birthline is displaced from the standard ZAMS, reflecting the time spent accreting mass before core burning starts.

Summary

The gist The scarcity of optically visible O stars near the theoretical zero-age main sequence (ZAMS) reflects prolonged accretion, suggesting that massive stars assemble over Myr timescales and their birthline is displaced from the standard ZAMS

The Problem and Hypothesis

Optically visible O stars are scarce near the theoretical zero-age main sequence (ZAMS), a deficit that has puzzled astronomers for over four decades Holgado et al. (2020) found a pronounced deficit near the ZAMS for initial masses of approximately 30–70 M⊙, while Schootemeijer et al. (2021) found that Small Magellanic Cloud stars above ≃ 20 M⊙ are similarly rare during the first half of core-H burning A complementary discrepancy appears during the embedded phase, where compact H ii regions are far more abundant relative to optically visible O stars than expected if they are freely expanding nebulae around stars that have completed their formation Both discrepancies may measure the duration of massive-star growth The paper tests whether the Galactic O-star ZAMS deficit can provide a second, independent clock for this growth

The Inertial-Inflow Model and Growth Law

The inertial-inflow model (IIM) of Padoan et al. (2020) provides a physically different picture, where massive-star growth is controlled by the turnover time of the turbulent flow on the scale from which the mass is collected The IIM predicts that formation time increases approximately as tform ∝ m f(1/2) This growth law is tested using accreting stellar-evolution models computed with MESA The model begins with the same pre-main-sequence seed model of mass mseed = 2 M⊙ and grows at the constant rate m acc (m f) = m f − mseed / tform(m f)

Constraints from Observational Data

The analysis compares the end-of-accretion points of accreting models with the full observed envelope of Galactic O stars in the Kiel diagram The hot-side boundary of this distribution is used to constrain growth-law parameters (alpha, tau0) The mismatch between the birthline and the fixed observational envelope is quantified by Q(alpha, tau0)

Joint Constraint and Conclusion

The Kiel-envelope analysis favors (alpha, tau0) = (0.52, 4.3 Myr), while the luminosity-function analysis gives posterior median values of (0.48, 4.15 Myr) The two measurements are mutually consistent because they probe the same underlying Galactic population and constrain the same parameters through different observables The combination using Fisher’s method yields projected 90% intervals 0.43 < alpha < 0.60 and 3.9 < tau0/Myr < 5.0, which is much tighter than either measurement alone The jointly preferred parameters define the effective Galactic growth law tform ≃ 4.3 Myr(m f/60 M⊙) 0.5 This work blurs the conventional boundary between massive-star formation and stellar evolution, showing that massive stars can spend a substantial fraction of their core-hydrogen-burning lifetime still assembling their mass

How it works

The analysis constructs birthlines by defining the formation time as a function of m f using the growth law tform(m f) = τ0 (m f/60 M⊙) α Non-rotating, solar-metallicity accreting models are computed with MESA, employing specific parameters for overshooting and mixing The resulting 24 distinct growth laws are sampled by computing 17 or 18 accreting tracks for different values of m f

Key Findings

The paper demonstrates that the accreting models naturally reproduce both the displacement and turnover of the hot-side envelope of Galactic O stars in the Kiel diagram The resulting Kiel-envelope constraint is shown in Fig. 2, with best-fitting parameters (alpha, tau0) = (0.52, 4.3 Myr) The combination of the Kiel-envelope and luminosity-function constraints substantially tightens the allowed parameter range The final joint 90% confidence intervals are projected as 0.43 < alpha < 0.60 and 3.9 < tau0/Myr < 5.0, based on the Fisher minimum at (alpha, tau0) = (5)

Limitations and Implications

The numerical confidence regions remain conditional on the adopted stellar models and observational inputs A separate systematic concern is that neglecting turbulent pressure in O-star atmospheres can lead to underestimated gravities, which would primarily move the higher-mass envelope downward in the Kiel diagram The value τ0 ≃ 4.3 Myr inferred here should be interpreted as an effective normalization for the Galactic O-star population rather than as a universal value for individual star-forming regions Future population modelling should explicitly include this environmental distribution The results suggest that incorporating stellar growth into models that presently begin from fully formed massive stars will have consequences well beyond the ZAMS deficit addressed here

References

Baraffe I., et al., 2023, MNRAS, 519, 5333 Beech M., Mitalas R., 1994, ApJS, 95, 517 Berlanas S. R., et al., 2020, A&A, 642, A168 Bernasconi P. A., Maeder A., 1996, A&A, 307, 829 Bonnell I. A., Bate M. R., 2006, MNRAS, 370, 488 Bonnell I. A., Bate M. R., Clarke C. J., Pringle J. E., 2001, MNRAS, 323, 785 Brott I., et al., 2011, A&A, 530 Carretero-Castrillo M., Benaglia P., Paredes J. M., Ribó M., 2025, Astronomy & Astrophysics, 694, A250 Castro N., Fossati L., Langer N., Simón-Díaz S., Schneider F. R. N., Izzard R. G., 2014, A&A, 570, L13 Choi J., Dotter A., Conroy C., Cantiello M., Paxton B., Johnson B. D., 2016, ApJ, 823, 102 Churchwell E., 2002, ARA&A, 40, 27 Dotter A., 2016, ApJS, 222, 8 Gieles M.

Improvements for AI systems

Improved stellar population modeling for massive stars by incorporating prolonged accretion histories into evolutionary tracks rather than assuming a standard zero-age main sequence (ZAMS). This allows AI systems to predict observed offsets in effective temperature and surface gravity based on accretion timescales, which is described as the end of massive-star growth does not coincide with the standard non-accreting ZAMS.
Enhanced constraint determination on massive star formation laws by combining constraints from different observational clocks, specifically the Kiel diagram analysis and luminosity function measurements. This allows AI to find transverse degeneracy directions and substantially tighten the allowed parameter range for growth laws like tform = tau0 (mf/60 Msun) alpha.
Development of a joint fitting routine using Fisher's method to combine constraints from multiple observables, such as the Kiel-envelope analysis and luminosity functions. This system can yield projected 90 percent confidence intervals (e g., 0.43 < alpha < 0.60 and 3.9 < tau0/Myr < 5.0) by minimizing a combined statistic like XF, fw = -2 ln pKiel, fw + ln pLF.
Modeling the environmental dependence of star formation rates by incorporating cloud-scale simulation predictions into the characteristic Galactic formation time (tau0). This allows AI to interpret tau0 as an effective normalization for the Galactic O-star population rather than as a universal value for individual star-forming regions, accounting for environmental spread in surface density.
Systematic correction module for spectroscopic gravity measurements that accounts for turbulent pressure support in O-star atmospheres. This module can adjust the inferred log gtrue values to better match theoretical models, which the paper suggests would primarily move the higher-mass envelope downward in the Kiel diagram.

Abstract

The scarcity of optically visible O stars near the theoretical zero-age main sequence (ZAMS) has puzzled astronomers for over four decades. We argue that it reflects prolonged accretion. If massive stars assemble over Myr timescales, as the inertial-inflow model predicts, core-hydrogen burning begins before accretion ends, and the stars become optically visible at the end of accretion on a birthline displaced from the standard ZAMS. We test this interpretation by comparing the hot-side envelope of Galactic O stars in the Kiel diagram with accreting stellar-evolution models, adopting a growth law t form=τ 0(m f/60 M) α with formation time t form and final mass m f. The resulting birthlines reproduce both the displacement and the high-mass turnover of the observed envelope, with best-fitting parameters (α,τ 0)=(0.52,4.3, Myr), consistent with independent constraints from the luminosity functions of compact H II regions and OB stars. The two constraints have orthogonal degeneracy directions, so combining them narrows the allowed parameter space, yielding projected 90% confidence intervals 0.43<α<0.60 and 3.9<τ 0/Myr<5.0. The inferred growth law agrees with cloud-scale simulations and the inertial-inflow model, and implies that mass assembly and main-sequence evolution overlap substantially in massive stars, with broader implications for stellar evolution and feedback.

Sources

You're Gonna Need a Bigger Core: Calibrating Massive Star Models against Galactic OB-type Stars