Compact HII Regions as Clocks of Massive-Star Formation: Evidence for Long Formation Timescales

We revisit the luminosity function (LF) of compact HII regions in the context of the inertial–inflow model, in which massive stars assemble over extended, mass-dependent timescales. The comparison of the compact-HII-region LF with that of OB stars is traditionally used to estimate the compact-HII-phase lifetime and is often cited as evidence for the classical ``lifetime problem" of HII regions. We show that once stellar growth during the ionizing phase is included, the LF comparison instead constrains massive-star formation timescales, so the lifetime problem turns into evidence for prolonged growth. We illustrate the principle with a simple analytic model and then forward-model the two LFs with Monte Carlo realizations. We also derive revised Galactic LFs for compact HII regions and OB stars from the Red MSX Source survey and the Alma Luminous Star catalogue. The joint LF constraints imply a growth law where the formation time is about 2 Myr for a $60,M_\odot$ star, with a square-root dependence on mass. The revised OB-star LF exhibits a statistically significant knee at $\log_{10}(L_{\rm k}/L_\odot)=5.0$, while the HII-region LF knee occurs at lower luminosity, as expected in the interpretation that HII regions are powered by stars that are still growing in mass. We conclude that massive stars in the Milky Way form over Myr timescales that increase with their final mass.

💡 Research Summary

In this paper the authors revisit the luminosity functions (LFs) of compact H II regions and OB stars in the Milky Way, interpreting them within the framework of the inertial‑inflow model (IIM) for massive‑star formation. The classical “lifetime problem” – the apparent discrepancy between the short dynamical expansion time of a compact H II region (∼10⁴ yr) and the much longer lifetime inferred from the ratio of compact H II regions to OB stars (∼3 × 10⁵ yr) – rests on the assumption that the ionizing star has already reached its final main‑sequence mass when the H II region becomes observable. The IIM instead posits that massive stars continue to accrete significant mass while they are already ionizing their surroundings; thus the ionizing luminosity L_HII is systematically lower than the final main‑sequence luminosity L_OB(m_f). Consequently, the appropriate LF comparison is not at a fixed luminosity but between ϕ_HII(L) and the OB‑star LF at a higher luminosity L′ > L, which naturally yields a longer effective age for the compact‑H II phase.

The authors formalize the IIM with two key relations. First, the formation time scales with final mass as t_form(m_f) = τ₀ (m_f/m₀)^α, where α reflects the turbulent velocity scaling and τ₀ is an environment‑dependent normalization. Second, the mass growth is approximately linear in time, m(t; m_f) ≈ m_f · t/t_form(m_f) for 0 ≤ t ≤ t_form. They adopt a broken‑power‑law initial mass function (IMF) with low‑mass slope s (≈2.35), high‑mass slope s_f, and a knee mass m_k set by the maximum cloud mass. By assuming power‑law relations L ∝ m^γ and main‑sequence lifetime t_MS ∝ m^−δ, they derive analytic expressions for the LF slopes: β_OB = 1 − s − δγ for the OB‑star LF and β_HII = 1 + α − sγ for the compact‑H II LF. The difference Δβ = β_HII − β_OB = α + δγ is independent of the IMF slope, meaning that the observed LF slope difference directly constrains the growth exponent α.

Using modern catalogs, the authors construct updated Galactic LFs. For OB stars they employ the third release of the ALMA Luminous Star catalog (ALS III), cross‑matched with Gaia DR3 and the Galactic O‑Star catalog, yielding 7,191 stars within 6 kpc. For compact H II regions they use the Red MSX Source (RMS) survey, extending to 18 kpc. Both samples are corrected for vertical exponential distribution (scale height h = 39 pc) and for luminosity‑dependent completeness via an empirical effective‑volume method. The resulting OB‑star LF shows two power‑law segments with slopes β_OB,1 = −0.91 ± 0.01 (below log L/L⊙ ≈ 5.0) and β_OB,2 = −1.63 ± 0.14 (above), indicating a statistically significant knee at log L_k,OB ≈ 5.0 (≈25 M⊙). The compact‑H II LF exhibits shallower slopes β_HII,1 = −0.26 ± 0.07 and β_HII,2 = −0.99 ± 0.09, with its knee at a lower luminosity log L_k,HII ≈ 4.6 (≈18 M⊙). The offset between the two knees is precisely what the IIM predicts: H II regions are powered by stars that have not yet reached their final mass, so at a given luminosity they correspond to a higher‑mass, later‑stage OB star.

To quantify the model parameters, the authors perform forward Monte‑Carlo simulations. They generate synthetic populations by drawing final masses from the broken‑power‑law IMF, assigning formation times via t_form(m_f) = τ₀ (m_f/m₀)^α, and computing the time spent in each luminosity bin during the ionizing phase. By fitting the simulated LFs to the observed ones, they find best‑fit values α ≈ 0.5 and τ₀ ≈ 2 Myr for a reference mass m₀ = 60 M⊙. This implies that a 60 M⊙ star takes on average ~2 Myr to assemble, and that formation time scales as the square root of the final mass (t_form ∝ m_f^0.5). Such timescales are an order of magnitude longer than the core‑collapse timescales (∼10⁵ yr) traditionally assumed, supporting the view that massive‑star formation is regulated by large‑scale turbulent inflows rather than rapid, isolated collapse.

The paper concludes that the “lifetime problem” of compact H II regions is in fact a manifestation of prolonged stellar growth. The LF shapes and the positions of their knees provide independent, observational constraints on the mass‑dependent formation timescales predicted by the inertial‑inflow model. Massive stars in the Milky Way therefore appear to form over several Myr, with more massive stars taking proportionally longer. This result has broad implications for interpreting the UV output of star‑forming galaxies, the feedback budget in molecular clouds, and the origin of the high‑mass end of the IMF. Future work should aim at refining distance and completeness corrections, extending the analysis to extragalactic samples, and directly probing the accretion signatures of ionizing massive protostars to further test the IIM framework.