OCCAM X. Neutron Capture Abundances with Keck/HIRES & Magellan/MIKE

Abstract: The chemistry of stars provides powerful insight into the history of the Milky Way. With multiple large-sky spectroscopic surveys that are currently available, using chemistry as a means to study the evolution and history of the Milky Way has flourished. Open clusters have long been used as landmarks to calibrate different age dating methods (e.g., gyrochronology and asteroseismology). In this work, we utilize the SDSS-IV/APOGEE-based Open Cluster Chemical Abundances and Mapping (OCCAM) survey as our foundation for new optical observations; enabling us to characterize neutron-capture abundances for known cluster members. For 56 stars in 18 open clusters, we collected high-resolution (R > 50,000), high-S/N (>75 at 5500A), spectra from Keck I and Magellan Baade telescopes. With these data, we derive abundances for 23 elements using BACCHUS, including 7 neutron capture abundances not measurable by APOGEE. Finally, we characterize the radial distribution of these neutron-capture elements in the Milky Way. We find that the second-peak s-process and r-process abundances exhibit relatively flat gradients in the Milky Way. Although not as distinct, the first-peak s-process abundances also have slopes which are shallower than the alpha and iron-peak elements. The differences in the neutron-capture gradients from the lighter elements not just indicates the sources producing these elements are fundamentally different, but that the timescales on which they are produced also differ (especially for the r-process). Moreover, a metallicity dependence of the AGB stars responsible for producing the heaviest s-process abundances may be necessary to consider in Galactic evolution models.

• The paper quantifies precise radial gradients of s- and r-process elements, showing flatter distributions compared to lighter nucleosynthetic families.
• Using high S/N spectra from Keck/HIRES and Magellan/MIKE along with the BACCHUS code, it establishes homogeneous measurements across 56 stars in 18 clusters.
• The paper highlights the influence of metallicity-dependent AGB yields and delayed r-process enrichment, providing robust constraints for Galactic chemical evolution models.

Neutron-Capture Element Abundance Gradients in the Galactic Disk from High-Resolution Spectroscopy

Introduction and Motivation

This study presents an analysis of neutron-capture element abundances in Milky Way open clusters based on high-S/N, high-resolution spectra obtained with Keck/HIRES and Magellan/MIKE, within the OCCAM (Open Cluster Chemical Abundances and Mapping) framework. The primary objective is to precisely characterize the concentrations and radial gradients of elements generated by slow and rapid neutron-capture (s- and r-) processes, disentangling their spatial distributions from those of lighter nucleosynthetic families (α, odd-Z, and Fe-peak elements), and thus refine chemical evolution constraints for the Galactic disk.

Previous studies leveraging large-scale spectroscopic surveys have provided robust constraints on metallicity and α-element gradients in the disk, but those works have been limited in the precision and coverage of heavy element abundances, particularly for n-capture species. With the deep, optical, high-resolution spectra in this study, it is possible to homogeneously measure up to 23 elements, including crucial n-capture elements (Sr, Y, Zr, La, Ba, Ce, Eu, Nd, Sm) that are not accessible to surveys such as APOGEE. The goal is to measure the spatial abundance gradients for these elements and compare them directly to lighter element behaviors, thereby probing the dependence of their production mechanisms on Galactic environment and stellar population properties.

Data, Targets, and Methodology

The sample comprises 56 stars across 18 open clusters with well-defined memberships based on the OCCAM catalog and Gaia EDR3 astrometry. Clusters span a Galactocentric guide radius ($R_\mathrm{guide}$) range of $6$ to $17$ kpc, encompassing both inner and outer disk environments. Stellar atmospheric parameters and detailed element abundances are derived using the BACCHUS code, exploiting spectral synthesis matched to the MARCS LTE model grid. A robust error analysis is undertaken via perturbation of model parameters and reproducibility checks on stars observed from both telescopes.

Clusters are chosen for reliable membership and include both well-sampled bright systems used as abundance calibrators, as well as fainter, more remote clusters granted robust membership probabilities and sufficient S/N. The open cluster mean abundances are computed as the average across validated members, and cluster properties such as age and distance are taken from recent neural-network based catalog measurements, ensuring homogeneous scaling.

Key Results: Radial Abundance Gradients

The principal focus is the radial distribution of n-capture element abundances, compared with metallicity, α, Fe-peak, and odd-Z elements:

• [Fe/H] exhibits a steep negative radial gradient of $-0.065 \pm 0.017$ dex kpc$^{-1}$, fully consistent with both high-resolution and large-sample cluster studies in the literature.
• First-peak s-process elements ([Sr/H], [Y/H], [Zr/H]) show negative gradients that are significantly shallower than $-0.033 \pm 0.021$ to $-0.043 \pm 0.020$ dex kpc$^{-1}$" title="" rel="nofollow" data-turbo="false" class="assistant-link">Fe/H, with Y systematically lower at all radii.
• Second-peak s-process elements ([La/H], [Ce/H]) display near-flat slopes ($-0.026 \pm 0.020$, $-0.015 \pm 0.020$ dex kpc$6$0), while [Ba/H] curiously follows a steeper gradient akin to [Fe/H]; the latter is attributed to strong line saturation and NLTE effects in more metal-rich stars (Ba NLTE corrections not applied here, cf. (2607.00291)).
• r-process elements ([Eu/H], [Sm/H]) manifest slopes that are either only marginally negative or statistically indistinguishable from zero within the uncertainties, while [Nd/H] is consistent with a perfectly flat gradient.
• Fe-peak and α-elements track the metallicity gradient closely, with only minor family-dependent deviations; exceptions include Cu and Zn, whose gradients mimic those of the first-peak s-process elements, reflecting their partial s-process production channels.

These results unambiguously demonstrate that heavy n-capture element gradients are systematically flatter than those of lighter nucleosynthetic families, in qualitative agreement with theoretical predictions for delayed production channels and metallicity-dependent yields in AGB stars for s-process species, and distinct delay-timescales and production rates for r-process elements.

Interpretation and Theoretical Implications

The study finds a clear separation in gradient steepness between nucleosynthetic families, directly tying the observed spatial trends to the dominant enrichment mechanisms:

• Type Ia and II SNe products (α and Fe-peak) build up efficiently in regions of rapid, continuing star formation, leading to steep negative abundance gradients as a function of Galactocentric radius.
• s-process elements, especially from the second peak, are produced in AGB stars. The AGB yields for heavy s-process elements have a strong inverse dependence on initial metallicity: stars with lower [Fe/H] possess greater neutron/seed ratios, favoring production of heavier nuclei (Ce, La, Ba). This metallicity sensitivity acts to flatten abundance gradients for these elements compared to Fe-peak nuclei.
• r-process elements, arising primarily from neutron star mergers and possibly rare core-collapse SNe, are characterized by long and stochastic DTDs, leading to more effective spatial mixing and thus minimal radial gradients. No strong evidence emerges for significant metallicity sensitivity in r-process production, in agreement with recent interpretations ([Kobayashi et al. 2023], [Thielemann et al. 2017]).
• The anomalous behavior of Ba is highlighted as likely due to line formation complexity and NLTE corrections not included here; correcting for these is likely to produce gradients more consistent with the global s-process element pattern, as observed in recent NLTE studies ([Liu et al. 2020]).

The observed flattening of n-capture element gradients is not plausibly attributed to radial migration, since most clusters in this sample are $6$1 Gyr old and thus have not experienced significant blurring effects. The effect is instead intrinsic to the metallicity and site dependence of the enrichment mechanisms.

Broader Impact and Future Prospects

The fine-grained empirical determination of n-capture gradients provides essential calibration data for theoretical models of Galactic chemical evolution. The measured flattening, and its correlation with the s- or r-process contribution fraction, constitutes a stringent test for AGB yield predictions and the assumed DTDs for neutron star mergers. The inability to explain r-process patterns via a single enrichment site or channel makes clear that further constraints on the merger rate and its metallicity dependence remain a priority ([kobayashi23-can-NSMs-explain-rprocess]).

These results also underscore the necessity to incorporate metallicity-dependent AGB star yields for heavy elements in any model of disk evolution. The cluster-based methodology, controlling for age and birth radius, remains uniquely powerful for this purpose, as demonstrated in recent works utilizing open cluster samples from SDSS, GALAH, and Gaia-ESO ([Magrini et al. 2023], [Otto et al. 2026], [Spina et al. 2021]).

As large-scale high-resolution spectroscopic surveys (SDSS-V, 4MOST, WEAVE) continue to expand the sample of faint and remote clusters with sufficient S/N, and as NLTE corrections become routinely incorporated, further improvements in the precision and interpretive power of heavy element abundance gradients are expected. The continued integration of chemical clocks anchored on s-process elements (e.g., Y/Mg, [C/N]) with abundance mapping will also resolve degeneracies in cluster age-labeling and Galactic disk assembly history.

Conclusion

This investigation has clearly demonstrated that n-capture element abundances in open clusters show systematically flatter or statistically insignificant gradients across the Galactic disk, in contrast to α- and Fe-peak elements, and with gradient steepness correlating to s-/r-process origin. The outcome highlights the combined effects of metallicity-dependent AGB yields and long-delay enrichment mechanisms, providing robust validation for chemical evolution models with such features. Further progress will require joint advances in spectroscopic sample coverage, the treatment of atomic physics (especially NLTE), and increasingly sophisticated comparisons to cosmological chemo-dynamical simulations.