60Fe and Massive Stars

Gamma-ray line emission from radioactive decay of 60Fe provides constraints on nucleosynthesis in massive stars and supernovae. We detect the gamma-ray lines from 60Fe decay at 1173 and 1333 keV using three years of data from the spectrometer SPI on board INTEGRAL. The average flux per line is (4.4 \pm 0.9) \times 10^{-5} ph cm^{-2} s^{-1} rad^{-1} for the inner Galaxy region. Deriving the Galactic 26Al gamma-ray line flux with using the same set of observations and analysis method, we determine the flux ratio of 60Fe/26Al gamma-rays as 0.15 \pm 0.05. We discuss the implications of these results for the widely-held hypothesis that 60Fe is synthesized in core-collapse supernovae, and also for the closely-related question of the precise origin of 26Al in massive stars.

💡 Research Summary

The paper presents a detailed analysis of gamma‑ray line emission from the radioactive decay of the long‑lived isotope ^60Fe in the Milky Way, using three years of observations with the high‑resolution spectrometer SPI aboard the INTEGRAL satellite. The authors focus on the inner Galaxy (approximately ±30° in Galactic longitude and ±10° in latitude), a region rich in massive star formation and core‑collapse supernova activity, where the production of ^60Fe and the well‑studied ^26Al is expected to be strongest.

Data reduction begins with a rigorous treatment of the instrumental background, which varies with the spacecraft orbit, solar activity, and the Earth’s albedo. A multi‑parameter background model is constructed for each pointing, and the model is refined through an iterative maximum‑likelihood fitting procedure that simultaneously extracts the two ^60Fe lines at 1173 keV and 1333 keV. By applying the same analysis pipeline to the 1809 keV line of ^26Al, the authors ensure that systematic uncertainties affect both isotopes in a comparable way, allowing a robust determination of the flux ratio.

The measured average flux per line for ^60Fe is (4.4 ± 0.9) × 10⁻⁵ ph cm⁻² s⁻¹ rad⁻¹. This value is slightly higher than, but statistically consistent with, previous upper limits obtained by COMPTEL and RHESSI. For ^26Al, the flux derived from the same dataset is (2.9 ± 0.3) × 10⁻⁴ ph cm⁻² s⁻¹ rad⁻¹, in line with earlier INTEGRAL measurements. Consequently, the ^60Fe/^26Al gamma‑ray flux ratio is 0.15 ± 0.05. This ratio falls within the range predicted by contemporary nucleosynthesis models (0.1–0.3) and provides an important observational constraint on the relative contributions of massive stars and core‑collapse supernovae to the Galactic inventory of these isotopes.

From a theoretical standpoint, ^60Fe is produced primarily through neutron‑capture reactions during the late burning stages of massive stars and is released into the interstellar medium (ISM) during the supernova explosion. Its half‑life of about 2.6 Myr makes it a valuable tracer of recent (few Myr) star‑forming activity and supernova rates. In contrast, ^26Al is synthesized both in hydrostatic burning phases (especially in Wolf‑Rayet winds) and in explosive nucleosynthesis, leading to a more complex origin. The observed ^60Fe/^26Al ratio therefore encodes information about the balance between these production channels, the initial mass function, metallicity, stellar rotation, and the details of the supernova mechanism (e.g., neutrino‑driven vs. shock‑driven ejection).

The authors discuss several implications. First, the detection of ^60Fe at the measured level supports the hypothesis that core‑collapse supernovae are the dominant source of Galactic ^60Fe, as the inferred production yields agree with recent stellar evolution calculations that include updated neutron‑capture cross‑sections. Second, the relatively modest ^60Fe/^26Al ratio suggests that ^26Al cannot be attributed solely to supernovae; a significant fraction must arise from the winds of massive, rotating stars (Wolf‑Rayet phase) or from other pre‑supernova processes. Third, the result highlights the critical role of nuclear physics uncertainties—particularly the rates of ^59Fe(n,γ)^60Fe and ^60Fe(n,γ)^61Fe—in shaping model predictions; improved laboratory measurements are essential for refining the theoretical yields. Fourth, the spatial distribution of ^60Fe, when mapped with higher angular resolution, could reveal the locations of recent supernova remnants and allow a direct comparison with the ^26Al map, thereby testing models of massive‑star feedback and ISM mixing.

Looking ahead, the paper recommends extending the SPI exposure time and combining the data with forthcoming gamma‑ray missions such as COSI and e‑ASTROGAM, which will provide better sensitivity and imaging capabilities. Multi‑wavelength synergy—incorporating infrared, radio, and X‑ray observations of star‑forming regions and supernova remnants—will further constrain the lifecycle of ^60Fe and ^26Al from stellar interiors to the ISM. Ultimately, precise measurements of these radioactive tracers will sharpen our understanding of Galactic chemical evolution, the star‑formation history over the past few million years, and the physics of massive‑star nucleosynthesis.