OB Stars Research Articles

Hydrodynamic simulation of Cygnus OB2: the absence of a cluster wind termination shock

ABSTRACT We perform a large-scale hydrodynamic simulation of a massive star cluster whose stellar population mimics that of the Cygnus OB2 association. The main-sequence stars are first simulated during 1.6 Myr, until a quasi-stationary state is reached. At this time, the three Wolf–Rayet stars observed in Cygnus OB2 are added to the simulation, which continues to 2 Myr. Using a high-resolution grid in the centre of the domain, we can resolve the most massive stars individually, which allows us to probe the kinetic structures at small (parsec) scales. We find that, although the cluster excavates a spherical ‘superbubble’ cavity, the stellar population is too loosely distributed to blow a large-scale cluster wind termination shock, and that collective effects from wind–wind interactions are much less efficient than usually assumed. This challenges our understanding of the ultra-high energy emission observed from the region.

The Long-period Spectroscopic Orbit and Dust Creation in the Wolf–Rayet Binary System WR 125

Several long-period binaries with a carbon-rich Wolf–Rayet star and an OB star produce dust in their wind collisions. In eccentric binaries, this is seen most strongly near periastron passage. The exact conditions leading to dust creation require the orbital properties to be determined, which is difficult owing to their long periods. Recently, the binary system WR 125 (WC7 + O9III) began a dust-creation episode seen through an infrared outburst first detected by NEOWISE-R, which was the first outburst detected since 1991. We present new near- and mid-infrared photometry, which we use to show consistency between the two outbursts and derive an orbital period of 28.12−0.05+0.10 yr. We use a long time series of optical spectra to place the first constraints on its orbital elements, on the assumption that this system will produce dust near periastron. The orbit has a mild eccentricity of 0.29 ± 0.12 and is only derived for the Wolf–Rayet component, as the O star’s radial velocities have noise that is likely larger than the expected semiamplitude of the orbit. We also present SOFIA/FORCAST grism spectroscopy to examine the infrared spectral energy distribution of the dust during this outburst, comparing its properties to other WCd binaries, and deriving a dust temperature of 580 K in 2021. This collection of observations will allow us to plan future observations of this system and place the system in the context of dust-creating Wolf–Rayet binaries.

Expanded Sample of Small Magellanic Cloud Ultraviolet Dust Extinction Curves: Correlations between the 2175 Å Bump, q PAH, Ultraviolet Extinction Shape, and N(H i)/A(V)

The Small Magellanic Cloud (SMC) shows a large variation in ultraviolet (UV) dust extinction curves, ranging from Milky Way (MW) like to significantly steeper curves with no detectable 2175 Å bump. This result is based on a sample of only nine sight lines. From Hubble Space Telescope Space Telescope Imaging Spectrograph and IUE spectra of OB stars, we have measured UV extinction curves along 32 SMC sight lines where eight of these curves were published previously. We find 16 sight lines with steep extinction with no detectable 2175 Å bump, four sight lines with MW-like extinction with a detectable 2175 Å bump, two sight lines with fairly flat UV extinction and weak/absent 2175 Å bumps, and 10 sight lines with unreliable curves due to low SMC dust columns. Our expanded sample shows that the sight lines with and without the 2175 Å bump are located throughout the SMC and not limited to specific regions. The average extinction curve of the 16 bump-less sight lines is very similar to the previous average based on four sight lines. We find no correlation between dust column and the strength of the 2175 Å bump. We test the hypothesis that the 2175 Å bump is due to the same dust grains that are responsible for the mid-infrared carbonaceous (polycyclic aromatic hydrocarbon) emission features and find they are correlated, confirming recent work in the MW. Overall, the slope of the UV extinction increases as the amplitudes of the 2175 Å bump and far-ultraviolet curvature decrease. Finally, the UV slope is correlated with N(H i)/A(V) and the 2175 Å bump and nonlinear far-ultraviolet rise amplitudes are anticorrelated with N(H i)/A(V).

Stellar Population near NGC 2021: Procession of Star Formation in the South Rim of Supergiant Shell LMC 4

Supergiant shells (SGSs) are the largest interstellar structures where heated and enriched gas flows into the host galaxy’s halo. The SGSs in the Large Magellanic Cloud (LMC) are so close that their stars can be resolved with ground-based telescopes to allow studies of star-formation history. Aiming to study the star formation history and energy budget of LMC 4, we have conducted a pilot study of the cluster NGC 2021 and the OB associations in its vicinity near the south rim of LMC 4. We use the Magellanic Cloud Photometric Survey data of the LMC to establish a methodology to examine the stellar population and assess the massive star formation history. We find a radial procession of massive star formation from the northwest part of the OB association LH79 through NGC 2021 to the OB association LH78 in the south. Using the stellar content of NGC 2021 and the assumption of Salpeter’s initial mass function, we estimate that ∼4 supernovae have occurred in NGC 2021, injecting at least 4 × 1051 erg of kinetic energy into the interior of LMC 4.

Runaway BN supergiant star HD 93840: Progenitor of an imminent core-collapse supernova above the Galactic plane

We present a quantitative spectral analysis of the extreme nitrogen-enhanced supergiant HD 93840 (BN1\,Ib) at an intermediate galactic latitude. Based on an optical high-resolution spectrum and complementary ultraviolet and infrared (spectro-)photometry, in addition to Gaia data, we carried out a full characterisation of the star's properties. We used both hydrostatic and unified (photosphere+wind) model atmospheres that account for deviations from local thermodynamic equilibrium. A highly unusual surface CNO-mixing signature and a marked stellar overluminosity compared to the mass imply a binary channel for the star's past evolution. The kinematics shows that it has reached its current position above the Galactic plane as a runaway star, likely ejected by the supernova explosion of its former companion star. Its current bulk composition, with a notably increased mean molecular weight due to core He- and progressed shell H-burning, suggests an advanced evolutionary stage. It is poised to yield a rare core-collapse supernova of a blue supergiant about ten OB star population scale heights above the Galactic disk relatively soon, contributing to the metal enrichment of the circumgalactic medium.

What Is the Nature of Little Red Dots and what Is Not, MIRI SMILES Edition

We study 31 little red dots (LRD) detected by JADES/NIRCam and covered by the SMILES/MIRI survey, of which ∼70% are detected in the two bluest MIRI bands and 40% in redder MIRI filters. The median/quartiles redshifts are z=6.95.97.7 (55% spectroscopic). The spectral slopes flatten in the rest-frame near-infrared, consistent with a 1.6 μm stellar bump but bluer than direct pure emission from active galactic nuclei (AGN) tori. The apparent dominance of stellar emission at these wavelengths for many LRDs expedites stellar mass estimation: the median/quartiles are logM⋆/M⊙=9.49.19.7 . The number density of LRDs is 10−4.0±0.1 Mpc−3, accounting for 14% ± 3% of the global population of galaxies with similar redshifts and masses. The rest-frame near-/mid-infrared (2–4 μm) spectral slope reveals significant amounts of warm dust (bolometric attenuation ∼3–4 mag). Our spectral energy distribution modeling implies the presence of <0.4 kpc diameter knots, heated by either dust-enshrouded OB stars or an AGN producing a similar radiation field, obscured by A(V) > 10 mag. We find a wide variety in the nature of LRDs. However, the best-fitting models for many of them correspond to extremely intense and compact starburst galaxies with mass-weighted ages 5–10 Myr, very efficient in producing dust, with their global energy output dominated by the direct (in the flat rest-frame ultraviolet and optical spectral range) and dust-recycled emission from OB stars with some contribution from an obscured AGN (in the infrared).

Multiwavelength study of the HII region LHA 120-N11 in the Large Magellanic Cloud with eROSITA

Aims. We studied the diffuse X-ray emission around the HII region LHA 120-N11, which is one of the most active star-forming regions in the Large Magellanic Cloud. We want to determine the nature of the diffuse X-ray emission and improve our understanding of its origin including related interactions with the cold interstellar medium. Methods. We analyzed the diffuse X-ray emission observed with the extended Roentgen Survey with an Imaging Telescope Array (eROSITA) on the Spectrum-Roentgen-Gamma mission to determine the physical properties of the hot diffuse X-ray emission. Four spectral extraction regions were defined based on the morphology of the X-ray emission. We also studied HI and CO data, as well as Hα line emission in the optical, and compared them with the properties of the diffuse X-ray emission. Results. The X-ray emission in the four regions is well fitted with an absorbed model consisting of thermal plasma models (vapec) yielding temperatures of kT = ~0.2 keV and kT = 0.8–1.0 keV. The comparison of the X-ray absorption column density and the hydrogen column density shows that the X-ray dark lane located north of N11 is apparently caused by the absorption by HI and CO clouds. By estimating the energy budget of the thermal plasma, we also investigated the heating mechanism of the X-ray emitting plasma. The energy of the diffuse X-ray emission in the superbubble which is a star-forming bubble with a radius of ~120 pc including OB associations LH9, LH10, LH11, and LH13 can be explained by heating from high-mass stars. In the surrounding regions we find that the energy implied by the X-ray emission suggests that additional heating might have been caused by shocks generated by cloud–cloud collisions.

Peering into the Wolf-Rayet phenomenon through [WO] and [WC] stars

Abstract Spectroscopic observations have shown for decades that the Wolf-Rayet (WR) phenomenon is ubiquitous among stars with different initial masses. Although much effort to understand the winds from massive WR stars has been presented in the literature, not much has been done for such type of stars in the low-mass range. Here we present an attempt to understand the winds from [WR]-type stars using results from spectral analyses with the full non-LTE stellar atmosphere code PoWR. These results are put into context with the properties of massive WR stars. We found that WC+[WC] stars and WO+[WO] stars create independent sequences in the mass-loss rate ($\dot{M}$) and modified wind momentum (Dmom) versus luminosity (L) diagrams. Our analysis indicates that even when the winds of WR and [WR] stars become optically thin, there is no breakdown of the general mass-loss trend, contrary to the observed “weak wind phenomenon” in OB stars. We report that all WR-type stars studied here broadly define single sequences in the wind efficiency (η) versus transformed mass-loss rate ($\dot{M}_\mathrm{t}$), the $\dot{M}_\mathrm{t}$–Teff diagram, and the $(L, T_\mathrm{eff}, \dot{M})$ space, which suggest these to be fundamental properties of the WR phenomenon (regardless of the mass range), at least for WR-type stars of the O and C sequences. Our analytical estimations could drive computations of future stellar evolution models for WR-type stars.

Searching for magnetar binaries disrupted by core-collapse supernovae

ABSTRACT Core-collapse supernovae (CCSNe) are considered the primary magnetar formation channel, with 15 magnetars associated with supernova remnants (SNRs). A large fraction of these should occur in massive stellar binaries that are disrupted by the explosion, meaning that $\sim 45~{{\ \rm per\ cent}}$ of magnetars should be nearby high-velocity stars. Here, we conduct a multiwavelength search for unbound stars, magnetar binaries, and SNR shells using public optical (uvgrizy bands), infrared (J, H, K, and Ks bands), and radio (888 MHz, 1.4 GHz, and 3 GHz) catalogues. We use Monte Carlo analyses of candidates to estimate the probability of association with a given magnetar based on their proximity, distance, proper motion, and magnitude. In addition to recovering a proposed magnetar binary, a proposed unbound binary, and 13 of 15 magnetar SNRs, we identify two new candidate unbound systems: an OB star from the Gaia catalogue we associate with SGR J1822.3−1606, and an X-ray pulsar we associate with 3XMM J185246.6 + 003317. Using a Markov Chain Monte Carlo simulation that assumes all magnetars descend from CCSNe, we constrain the fraction of magnetars with unbound companions to $5\lesssim f_u \lesssim 24~{{\ \rm per\ cent}}$, which disagrees with neutron star population synthesis results. Alternate formation channels are unlikely to wholly account for the lack of unbound binaries as this would require $31\lesssim f_{nc} \lesssim 66~{{\ \rm per\ cent}}$ of magnetars to descend from such channels. Our results support a high fraction ($48\lesssim f_m \lesssim 86~{{\ \rm per\ cent}}$) of pre-CCSN mergers, which can amplify fossil magnetic fields to preferentially form magnetars.

Runaway OB Stars in the Small Magellanic Cloud. III. Updated Kinematics and Insights into Dynamical versus Supernova Ejections

We use the kinematics of field OB stars to estimate the frequencies of runaway stars generated by the dynamical ejection scenario (DES), the binary supernova scenario (BSS), and the combined two-step mechanism. We update the proper motions for field OB and OBe stars in the Small Magellanic Cloud (SMC) using Gaia DR3. Our sample now contains 336 stars from the Runaways and Isolated O-Type Star Spectroscopic Survey of the SMC, and we update our algorithm to calculate more accurate velocities compared to those obtained previously from DR2. We find a decrease in median velocity from 39 to 29 km s−1, implying that the proper motions from our previous work were systematically overestimated. We present the velocity distribution for OBe stars and quantitatively compare it to those of noncompact binaries and high-mass X-ray binaries. We confirm that OBe stars appear to be dominated by the BSS and are likely post-SN binary systems, further supporting the mass-transfer model to explain the origin of their emission-line disks. In contrast, normal OB stars may show a bimodal velocity distribution, which may be expected from different processes that occur with dynamical ejections. The kinematics of fast-rotating OB stars are similar to those of normal OB stars rather than OBe stars, suggesting that the origin of their high vrsini is different from that of OBe stars. We update our model parameters describing the kinematic origins of the SMC field population, still confirming that for runaway stars, the DES mechanism dominates, and two-step ejections seem comparable in frequency to pure BSS ejections.

Variability of Blue Supergiants in the LMC with TESS

The blue supergiant (BSG) problem, namely, the overabundance of BSGs inconsistent with classical stellar evolution theory, remains an open question in stellar astrophysics. Several theoretical explanations have been proposed, which may be tested by their predictions for the characteristic time variability. In this work, we analyze the light curves of a sample of 20 BSGs obtained from the Transiting Exoplanet Survey Satellite (TESS) mission. We report a characteristic signal in the low-frequency (f ≲ 2 day−1) range for all our targets. The amplitude spectrum has a peak frequency of ∼0.2 day−1, and we are able to fit it by a modified Lorentzian profile. The signal itself shows strong stochasticity across different TESS sectors, suggesting its driving mechanism happens on short (≲months) timescales. Our signals resemble those obtained for a limited sample of hotter OB stars and yellow supergiants, suggesting their possible common origins. We discuss three possible physical explanations: stellar winds launched by rotation, convection motions that reach the stellar surface, and waves from the deep stellar interior. The peak frequency of the signal favors processes related to the convective zone caused by the iron opacity peak, and the shape of the spectra might be explained by the propagation of high-order, damped gravity waves excited from that zone. We discuss the uncertainties and limitations of all these scenarios.

Implications on star formation rate indicators from H ii regions and diffuse ionized gas in the M101 Group

ABSTRACT We examine the connection between diffuse ionized gas (DIG), H ii regions, and field O and B stars in the nearby spiral M101 and its dwarf companion NGC 5474 using ultra-deep H α narrow-band imaging and archival GALEX UV imaging. We find a strong correlation between DIG H α surface brightness and the incident ionizing flux leaked from the nearby H ii regions, which we reproduce well using simple cloudy simulations. While we also find a strong correlation between H α and co-spatial far-ultraviolet (FUV) surface brightness in DIG, the extinction-corrected integrated UV colours in these regions imply stellar populations too old to produce the necessary ionizing photon flux. Combined, this suggests that H ii region leakage, not field OB stars, is the primary source of DIG in the M101 Group. Corroborating this interpretation, we find systematic disagreement between the H α- and FUV-derived star formation rates (SFRs) in the DIG, with SFRH α &amp;lt;SFRFUV everywhere. Within H ii regions, we find a constant SFR ratio of 0.44 to a limit of ∼10−5 M⊙ yr−1. This result is in tension with other studies of star formation in spiral galaxies, which typically show a declining SFRH α/SFRFUV ratio at low SFR. We reproduce such trends only when considering spatially averaged photometry that mixes H ii regions, DIG, and regions lacking H α entirely, suggesting that the declining trends found in other galaxies may result purely from the relative fraction of diffuse flux, leaky compact H ii regions, and non-ionizing FUV-emitting stellar populations in different regions within the galaxy.

The persistence of high altitude non-equilibrium diffuse ionized gas in simulations of star-forming galaxies

ABSTRACT Widespread, high altitude, diffuse ionized gas with scale heights of around a kiloparsec is observed in the Milky Way and other star-forming galaxies. Numerical radiation-magnetohydrodynamic simulations of a supernova-driven turbulent interstellar medium show that gas can be driven to high altitudes above the galactic mid-plane, but the degree of ionization is often less than inferred from observations. For computational expediency, ionizing radiation from massive stars is often included as a post-processing step assuming ionization equilibrium. We extend our simulations of an Milky Way-like interstellar medium to include the combined effect of supernovae and photoionization feedback from mid-plane OB stars and a population of hot evolved low mass stars. The diffuse ionized gas has densities below 0.1 ${\rm \,cm^{-3}}$, so recombination time-scales can exceed millions of years. Our simulations now follow the time-dependent ionization and recombination of low density gas. The long recombination time-scales result in diffuse ionized gas that persists at large altitudes long after the deaths of massive stars that produce the vast majority of the ionized gas. The diffuse ionized gas does not exhibit the large variability inherent in simulations that adopt ionization equilibrium. The vertical distribution of neutral and ionized gas is close to what is observed in the Milky Way. The volume filling factor of ionized gas increases with altitude resulting in the scale height of free electrons being larger than that inferred from H $\alpha$ emission, thus reconciling the observations of ionized gas made in H $\alpha$ and from pulsar dispersion measurements.

Thermal X-ray emission in the western half of the LMC superbubble 30 Dor C

ABSTRACT While 30 Dor C is a unique superbubble in the Large Magellanic Cloud for its luminous non-thermal X-ray emission, the thermal X-ray emission it emanates has not yet been thoroughly investigated and well constrained. Based on the separate ∼1 Ms deep XMM–Newton and Chandra observations, we report the discovery of the thermally emitting plasma in some portions of the western half of 30 Dor C. The thermal emission can be reproduced by a collisional-ionization-equilibrium plasma model with an average electron temperature of ∼0.4 keV. We find a significant overabundance of the intermediate-mass elements such as O, Ne, Mg, and Si, which may be indicative of a recent supernova explosion in 30 Dor C. Dynamical properties in combination with the information of the OB association LH 90 suggest that the internal post-main-sequence stars dominate the power of the superbubble and blow it out in the past ∼1 Myr.

Stellar wind bubbles of OB stars as Galactic cosmic-ray re-accelerators

Abstract Cosmic rays are highly energetic messengers propagating in magnetized plasma, which are, possibly but not exclusively, accelerated at astrophysical shocks. Amongst the variety of astrophysical objects presenting shocks, the huge circumstellar stellar wind bubbles forming around very massive stars, are potential non-thermal emitters. We present the 1D magneto-hydrodynamical simulation of the evolving magnetized surroundings of a single, OB-type main-sequence $60\, \rm M_{\odot }$ star, which is post-processed to calculate the re-acceleration of pre-existing non-thermal particles of the Galactic cosmic ray background. It is found that the forward shock of such circumstellar bubble can, during the early phase ($1\, \rm Myr$) of its expansion, act as a substantial re-accelerator of pre-existing interstellar cosmic rays. This results in an increasing excess emission flux by a factor of 5, the hadronic component producing γ-rays by π0 decay being more important than those by synchrotron and inverse Compton radiation mechanisms. We propose that this effect is at work in the circumstellar environments of massive stars in general and we conjecture that other nebulae such as the stellar wind bow shocks of runaway massive stars also act as Galactic cosmic-ray re-accelerators. Particularly, this study supports the interpretation of the enhanced hadronic emission flux measured from the surroundings of κ Ori as originating from the acceleration of pre-existing particles at the forward shock of its wind bubble.

A Self-consistent Data-driven Model for Determining Stellar Parameters from Optical and Near-infrared Spectra

Data-driven models, which apply machine learning to infer physical properties from large quantities of data, have become increasingly important for extracting stellar properties from spectra. In general, these methods have been applied to data in one wavelength regime or another. For example, APOGEE Net has been applied to near-IR spectra from the Sloan Digital Sky Survey (SDSS)–V APOGEE survey to predict stellar parameters (T eff, log g, and [Fe/H]) for all stars with T eff from 3000 to 50,000 K, including pre-main-sequence stars, OB stars, main-sequence dwarfs, and red giants. The increasing number of large surveys across multiple wavelength regimes provides the opportunity to improve data-driven models through learning from multiple data sets at once. In SDSS-V, a number of spectra of stars will be observed not just with APOGEE in the near-IR, but also with BOSS in the optical regime. Here, we aim to develop a complementary model, BOSS Net, that will replicate the performance of APOGEE Net in these optical data through label transfer. We further improve the model by extending it to brown dwarfs, as well as white dwarfs, resulting in a comprehensive coverage between 1700 < T eff < 100,000 K and 0 < logg < 10, to ensure BOSS Net can reliably measure parameters of most of the commonly observed objects within this parameter space. We also update APOGEE Net to achieve a comparable performance in the near-IR regime. The resulting models provide a robust tool for measuring stellar evolutionary states, and, in turn, enable characterization of the star-forming history of the Galaxy.

Outward migration in Nascent stellar groups

ABSTRACT As a stellar group forms within its parent molecular cloud, new members first appear in the deep interior. These overcrowded stars continually diffuse outward to the cloud boundary, and even beyond. Observations have so far documented only the interior drift. Those stars that actually leave the cloud form an expanding envelope that I call the ‘stellar mantle.’ Simple fluid models for the cloud and mantle illustrate their basic structure. The mantle’s expansion speed is subsonic with respect to the cloud’s dynamical temperature. I describe, in qualitative terms, how the expanding mantle and Galactic tidal radius might together shape the evolution of specific types of stellar groups. The massive stars in OB associations form in clouds that contract before extruding a substantial mantle. In contrast, the more slowly evolving clouds forming open clusters and T associations have extended mantles that encounter a shrinking tidal radius. These clouds are dispersed by internal stellar outflows. If the remaining group of stars is gravitationally bound, it appears as a long-lived open cluster, truncated by the tidal radius. If the group is unbound, it is a late-stage T association that will soon be torn apart by the tidal force. The ‘distributed’ populations of pre-main sequence stars observed in the outskirts of several star-forming regions are too distant to be stellar mantles. Rather, they could be the remnants of especially low-mass T associations.

Massive OB Star Clusters in Cassiopeia

We present a new investigation of several sparse OB associations in Cassiopeia (Cas OB2, OB4, OB5, OB7, and OB8) to revisit their membership and kinematics in the Gaia era. We compare the spatial distribution of the OB associations’ members, the known OB-type stars from the Reed catalog, and candidate hot stars from the Gaia Data Release 3. We also investigate the reddening distributions across the region to further understand the kinematics of these groups relative to the region’s molecular clouds.

Impacts of the Local arm on the local circular velocity inferred from the Gaia DR3 young stars in the Milky Way

ABSTRACT A simple one-dimensional axisymmetric disc model is applied to the kinematics of O type and B type stars (OB stars) near the Sun obtained from Gaia Data Release 3 catalogue. The model determines the ‘local centrifugal speed’ Vc(R0) – defined as the circular velocity in the Galactocentric rest frame, where the star would move in a near-circular orbit if the potential is axisymmetric with the local potential of the Galaxy. We find that the Vc(R0) values and their gradient vary across the selected region of stars within the solar neighbourhood. By comparing with an N-body/hydrodynamic simulation of a Milky Way-like galaxy, we find that the kinematics of the young stars in the solar neighbourhood is affected by the Local arm, which makes it difficult to measure Vc(R0). However, from the resemblance between the observational data and the simulation, we suggest that the known rotational velocity gap between the Coma Bernices and Hyades-Pleiades moving groups could be driven by the co-rotation resonance of the Local arm, which can be used to infer the azimuthally averaged circular velocity. We find that Vc(R) obtained from the D &amp;lt; 2 kpc sample is well matched with this gap at the position of the Local arm. Hence, we argue that our results from the D &amp;lt; 2 kpc sample, Vc(R0) = 234 ± 2 km s−1, are close to the azimuthally averaged circular velocity rather than the local centrifugal speed, which is influenced by the presence of the Local arm.

Origin of 60Fe nuclei in cosmic rays: the contribution of local OB associations

ABSTRACT The presence of live 60Fe nuclei (lifetime of 3.8 Myr) in cosmic rays detected by the Advanced Composition Explorer/Cosmic Ray Isotope Spectrometer instrument suggests a nearby nucleosynthesis source. 60Fe is primarily produced in core-collapse supernovae, and we aim to clarify whether the detected 60Fe nuclei can be associated with a particular local supernova. We consider 25 OB associations and subgroups located within 1 kpc of the Solar system based on recent Gaia census. A model is developed that combines stellar population synthesis within these OB associations, cosmic ray acceleration within associated superbubbles, and cosmic ray transport to the Solar system. The most critical model parameter impacting 60Fe cosmic ray production is the explodability criterion, which determines if a massive star ends its life as a supernova. Our study points to the Scorpius–Centaurus (Sco–Cen) OB association as the most probable origin of the observed 60Fe nuclei, particularly suggesting they were accelerated in the Sco–Cen superbubble by a young supernova aged ≤500 kyr with a progenitor mass of approximately 13–20 M⊙. A less likely source is the supernova at the origin of the Geminga pulsar 342 kyr ago, if the progenitor originated in the Orion OB1 association. The contribution of local OB associations to the cosmic ray density of stable 56Fe is estimated to be around 20 per cent, with some sensitivity to cosmic ray acceleration efficiency and diffusion coefficient. These findings shed light on the origins of cosmic ray nuclei, connecting them to nucleosynthesis events within our local cosmic neighbourhood.