Probing the evolving massive star population in Orion with kinematic and radioactive tracers

We assemble a census of the most massive stars in Orion, then use stellar isochrones to estimate their masses and ages, and use these results to establish the stellar content of Orion’s individual OB associations. From this, our new population synthesis code is utilized to derive the history of the emission of UV radiation and kinetic energy of the material ejected by the massive stars, and also follow the ejection of the long-lived radioactive isotopes 26Al and 60Fe. In order to estimate the precision of our method, we compare and contrast three distinct representations of the massive stars. We compare the expected outputs with observations of 26Al gamma-ray signal and the extent of the Eridanus cavity. We find an integrated kinetic energy emitted by the massive stars of 1.8(+1.5-0.4)times 10^52 erg. This number is consistent with the energy thought to be required to create the Eridanus superbubble. We also find good agreement between our model and the observed 26Al signal, estimating a mass of 5.8(+2.7-2.5) times 10^-4 Msol of 26Al in the Orion region. Our population synthesis approach is demonstrated for the Orion region to reproduce three different kinds of observable outputs from massive stars in a consistent manner: Kinetic energy as manifested in ISM excavation, ionization as manifested in free-free emission, and nucleosynthesis ejecta as manifested in radioactivity gamma-rays. The good match between our model and the observables does not argue for considerable modifications of mass loss. If clumping effects turn out to be strong, other processes would need to be identified to compensate for their impact on massive-star outputs. Our population synthesis analysis jointly treats kinematic output and the return of radioactive isotopes, which proves a powerful extension of the methodology that constrains feedback from massive stars.

💡 Research Summary

The paper presents a comprehensive population‑synthesis study of the massive star content in the Orion region, focusing on the four principal OB associations (Orion A, B, C, D). The authors first compile a census of the most massive O‑type and early‑B stars using recent optical and infrared surveys, and then determine individual stellar masses and ages by fitting modern stellar isochrones (e.g., PARSEC, Geneva) that incorporate appropriate metallicity and rotation effects. To assess methodological robustness, three distinct representations of the massive‑star population are explored: (1) a “star‑by‑star” approach that uses the measured mass and age of each object, (2) a “cluster” approach that treats each association as a single stellar population characterized by an initial mass function (IMF) and a mean age, and (3) a “global” approach that models the entire Orion complex as a continuous IMF with an age distribution. All three schemes are processed through a newly developed population‑synthesis code that tracks, as a function of time, the emission of Lyman‑continuum photons, the kinetic energy injected by stellar winds and supernovae, and the ejection of the long‑lived radioactive isotopes ^26Al and ^60Fe. The code adopts up‑to‑date mass‑loss prescriptions (including Vink et al. wind recipes and rotational enhancements) and the latest nuclear reaction network yields, allowing a realistic estimate of the isotopic production.

The model outputs are directly compared with two independent observational constraints. First, the 1.809 MeV gamma‑ray line measured by INTEGRAL/SPI indicates a total ^26Al mass of roughly 5.8 × 10⁻⁴ M⊙ in Orion. The synthesis predicts 5.8(+2.7 – 2.5) × 10⁻⁴ M⊙, in excellent agreement. Second, the kinetic energy required to inflate the Eridanus superbubble, inferred from its size and expansion velocity, is estimated at 1–3 × 10⁵² erg. The integrated kinetic energy released by the massive stars over the past ∼10 Myr is 1.8(+1.5 – 0.4) × 10⁵² erg, again matching the observational requirement. These concordances suggest that the standard mass‑loss rates used in the models are sufficient; no drastic reduction (as would be implied by strong wind clumping) is needed to reproduce the data.

Nevertheless, the authors acknowledge that if future studies confirm very strong clumping, the effective mass‑loss rates could be lower, potentially leading to a shortfall in UV luminosity, kinetic feedback, and ^26Al production. In that scenario, additional processes—such as binary mass transfer, enhanced pre‑supernova eruptions, or contributions from lower‑mass stars—would have to be invoked to compensate.

Overall, the study demonstrates the power of a joint kinematic and radioactive‑tracer approach to constrain massive‑star feedback. By simultaneously reproducing three distinct observables—ISM excavation (kinetic energy), ionizing radiation (free‑free emission), and nucleosynthetic ejecta (gamma‑ray lines)—the authors provide a self‑consistent picture of how massive stars shape their environment in Orion. The methodology is readily extensible to other star‑forming complexes, offering a valuable framework for future investigations of stellar feedback, chemical enrichment, and the dynamics of the interstellar medium.