Towards a unified injection model of short-lived radioisotopes in N-body simulations of star-forming regions

Recent research provides compelling evidence that the decay of short-lived radioisotopes (SLRs), such as 26Al, provided the bulk of energy for heating and desiccation of volatile-rich planetesimals in the early Solar System. However, it remains unclear whether the early Solar System was highly enriched relative to other planetary systems with similar formation characteristics. While the Solar System possesses an elevated level of SLR enrichment compared to the interstellar medium, determining SLR enrichment of individual protoplanetary disks observationally has not been performed and is markedly more difficult. We use N-body simulations to estimate enrichment of SLRs in star-forming regions through two likely important SLR sources, stellar winds from massive stars and supernovae. We vary the number of stars and the radii of the star-forming regions and implement two models of stellar wind SLR propagation for the radioisotopes 26Al and 60Fe. We find that for 26Al enrichment the Solar System is at the upper end of the expected distribution, while for the more supernovae dependent isotope 60Fe we find that the Solar System is comparatively very highly enriched. Furthermore, combined with our previous research, these results suggest that the statistical role of 26Al-driven desiccation on exoplanet bulk composition may be underestimated in typical interpretations of the low-mass exoplanet census, and that 60Fe is even less influential as a source of heating than previously assumed.