pith. machine review for the scientific record.
sign in

arxiv: 2604.15197 · v1 · submitted 2026-04-16 · 🌌 astro-ph.SR · astro-ph.GA

Stellar nucleosynthesis in the era of large surveys: s-process polluted binaries in GALAH DR4

A. Escorza , S. Vitali , D. Godoy-Rivera , S. Shetye , H. Van Winckel , G. Bustos , S. Goriely , L. Siess
show 6 more authors
M. Abdul-Masih T. Masseron D. A. Garc\'ia-Hern\'andez A. Ardern-Arentsen P. Jofr\'e S. Van Eck
This is my paper

Pith reviewed 2026-05-10 09:51 UTC · model grok-4.3

classification 🌌 astro-ph.SR astro-ph.GA
keywords s-process elementsbinary starsstellar abundancesGALAH surveychemically peculiar starsAGB starsmass transfer
0
0 comments X

The pith

GALAH DR4 yields 1073 new candidate stars showing s-process element overabundances from binary mass transfer.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper compiles a large, homogeneous sample of stars with s-process overabundances by applying abundance thresholds calibrated on a high-resolution validation set to the lower-resolution GALAH DR4 catalog. These chemically peculiar stars are expected to form when a companion on the asymptotic giant branch transfers s-process enriched material. A sympathetic reader cares because the resulting list is nearly five times larger than the set of previously confirmed cases and shares the same heavy-to-light s-element ratios, giving a stronger basis for testing nucleosynthesis and binary-evolution calculations. The work also notes that the candidates already display the elevated binary fraction predicted for this formation channel.

Core claim

Using thresholds on [s/Fe], [Y/Fe], [Zr/Fe], [Ba/Fe] and [La/Fe] derived from a validation sample of UVES and HERMES spectra, the authors flag 1073 stars in GALAH DR4 as s-process-polluted candidates. These objects span a wide range of stellar parameters and exhibit [hs/ls] ratios that closely match those of known barium and related stars, while their measured binary fraction exceeds that of the full survey catalog. The sample is almost five times larger than the number of currently confirmed s-process-rich stars and has been processed uniformly, opening systematic comparisons with models of AGB nucleosynthesis and binary mass transfer.

What carries the argument

Abundance thresholds on s-process elements calibrated from high-resolution validation spectra and applied to the automatic GALAH DR4 analysis pipeline.

If this is right

  • The sample size permits statistical studies of s-process pollution across a broad range of stellar parameters.
  • Uniform GALAH processing enables direct, quantitative comparison of observed abundances with AGB nucleosynthesis yields.
  • The elevated binary fraction among candidates supports the mass-transfer formation channel.
  • Limited orbital data for only 7 percent of the sample points to the need for radial-velocity monitoring to measure periods and eccentricities.
  • The catalog supplies a ready target list for confronting binary-evolution models with observed occurrence rates.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • High-resolution follow-up of the full candidate list could tighten the abundance cuts and reduce contamination.
  • The same threshold method could be applied to other large surveys to produce even bigger homogeneous samples.
  • The candidates offer a new route to measure the efficiency of s-process production in low-mass AGB stars through population statistics.

Load-bearing premise

Abundance thresholds derived from a small set of high-resolution spectra can be transferred directly to the lower-resolution, automatically analyzed GALAH DR4 data without creating large numbers of false positives or missing true cases.

What would settle it

High-resolution follow-up spectra of a random subset of the 1073 candidates that fail to show the expected s-process enhancements at the levels predicted by the thresholds.

Figures

Figures reproduced from arXiv: 2604.15197 by A. Ardern-Arentsen, A. Escorza, D. A. Garc\'ia-Hern\'andez, D. Godoy-Rivera, G. Bustos, H. Van Winckel, L. Siess, M. Abdul-Masih, P. Jofr\'e, S. Goriely, S. Shetye, S. Van Eck, S. Vitali, T. Masseron.

Figure 1
Figure 1. Figure 1: The high-resolution validation sample observed with HERMES (green circles) and UVES (purple squares) plotted over a clean GALAH DR3 catalogue (Buder et al. 2021) containing stars for which the spectroscopic analysis was reliable and good-quality abundances of Y, Zr, Ba, and La were obtained. Note that we plot the HR sample against the GALAH DR3 catalogue rather than the GALAH DR4 equivalent because the HR … view at source ↗
Figure 2
Figure 2. Figure 2: Example spectra of three stars from the HERMES (first row) and UVES (second and third rows) samples spanning different stellar parameters, illustrating the overall quality of our line fits. The observed spectrum is shown as a black dashed line, while the best-fitting synthetic spectrum of the single absorption line is overplotted with solid green (for HERMES) or purple (for UVES) lines. The shaded region m… view at source ↗
Figure 3
Figure 3. Figure 3: Difference (this work - GALAH DR4) between our HR results and the values published in GALAH DR4 (Buder et al. 2025). The panels show, from left to right, the differences in effective temperature, surface gravity and metallicity. Each panel illustrates HERMES targets with circles and UVES targets with squares, and the symbols are colour-coded as a function of the derived metallicity. The yellow regions indi… view at source ↗
Figure 4
Figure 4. Figure 4: Same as [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Kiel diagram, Teff -[Fe/H] diagram, and [Fe/H]-[s/Fe] diagram for our s-process rich stars (orange circles). We plot it over the full heavy￾metal catalogue (gray dots) for reference. Each panel includes the typical errorbar that GALAH DR4 reports for these parameters [PITH_FULL_IMAGE:figures/full_fig_p009_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: [Y/Fe], [Zr/Fe], [Ba/Fe], and [La/Fe] for our s-rich catalogue (orange circles). As in [PITH_FULL_IMAGE:figures/full_fig_p009_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Comparison between our s-rich catalogue (orange circles) and the Ba stars presented in de Castro et al. (2016) (black crosses) in the [s/Fe]-metallicity plane. our candidate sample by a factor of ∼2, and the fact that we can reproduce the [hs/ls] distribution with s-process nucleosynthesis models (see Sect. 7.1 below) indicates that this is not necessary. 7. Discussion 7.1. Comparison with confirmed Ba sta… view at source ↗
Figure 9
Figure 9. Figure 9: shows the 77 stars in our s-rich catalogue with periods and eccentricities in the NSS catalogue together with the Ba stars published by Jorissen et al. (2019); Escorza et al. (2019b); North et al. (2000); Escorza & De Rosa (2023), the CH stars from Jorissen et al. (2019); Escorza et al. (2019b) and the CEMP-s from Jorissen et al. (2016); Hansen et al. (2016) for ref￾erence. We crossmatched our s-rich candi… view at source ↗
Figure 10
Figure 10. Figure 10: Radial velocity variability of the stars monitored with HERMES. TYC 179-2074-1 is not included because we only have one spectrum. rently confirmed polluted stars. About 700 of our candidates have metallicities representative of Ba stars ([Fe/H] ≳ −0.5), while most of the rest have −1.25 ≲ [Fe/H] ≲ −0.5, which would traditionally classify them as CH stars. The most impor￾tant difference between our candida… view at source ↗
read the original abstract

Binary interactions during the AGB phase can lead to the formation of chemically peculiar stars with overabundances of s-process elements. Only a few hundreds of these stars have been subject to detailed chemical or dynamical studies. This work aims at compiling a systematic sample of s-process-polluted candidates using GALAH DR4. We also want to compare their properties with those of confirmed s-polluted stars to have stronger evidence of their nature. GALAH DR4 uses neural networks and automatic spectral analysis methods as well as data of a lower spectral resolution than normally used to characterise these objects. Because of this, we built a validation sample, for which we obtained UVES@VLT and HERMES@Mercator high-resolution spectra. We compare our stellar parameters and abundances with those of the survey and use this validation to define the thresholds that a star in GALAH DR4 must pass to be flagged as a good s-process-rich candidate. Based on our comparisons, we define thresholds on [s/Fe], [Y/Fe], [Zr/Fe], [Ba/Fe], and [La/Fe]. We identified 1073 stars in GALAH DR4 that are good candidates to be s-process polluted stars, covering a broad parameter space. They share many similarities with the samples of confirmed s-rich stars, especially their ratios of heavy over light s-elements ([hs/ls]), which strengthen our confidence in the purity of the sample. We find that only 7% of the candidates have measured orbital periods and eccentricities, limiting for now a full comparison with confirmed Ba and related stars. However, their binary fraction is, as expected, higher than the one we found for the full GALAH DR4 catalogue. Our sample of candidates is almost five times larger than the number of currently confirmed polluted stars. This and the fact that it has been homogeneously treated by GALAH open very interesting avenues to confront nucleosynthesis and binary evolution models.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper compiles a sample of 1073 s-process-polluted candidate stars from GALAH DR4 by calibrating abundance thresholds on [s/Fe], [Y/Fe], [Zr/Fe], [Ba/Fe], and [La/Fe] against a validation set of high-resolution UVES and HERMES spectra. The candidates are shown to share key properties with confirmed s-rich stars, especially [hs/ls] ratios, and exhibit a higher binary fraction than the parent GALAH catalogue, providing a homogeneous dataset five times larger than existing confirmed samples for nucleosynthesis and binary-evolution studies.

Significance. If the selection holds, the work supplies a large, uniformly analyzed sample that can statistically constrain AGB s-process yields and binary mass-transfer models, extending beyond the few hundred previously studied objects while leveraging the survey's broad parameter coverage.

major comments (2)
  1. [§4] §4 (threshold definition and application): The abundance thresholds are set by comparing the small validation sample to GALAH DR4 values, yet the manuscript provides no quantitative propagation of the differing error distributions, line-blending effects, and parameter covariances between R~28,000 automated analysis and higher-resolution spectra into an estimated false-positive rate when the cuts are applied across the full GALAH parameter space.
  2. [§5.1–5.2] §5.1–5.2 (post-selection validation): The reported similarity in [hs/ls] ratios is presented as supporting evidence for sample purity, but because this diagnostic is evaluated only on stars already passing the thresholds, it does not independently test the reliability of the threshold transfer or quantify contamination from stars whose abundance errors are correlated with the selection criteria.
minor comments (2)
  1. [§3] The text should explicitly state the size and parameter coverage of the validation sample (e.g., number of stars, range in Teff, log g, [Fe/H]) so readers can judge how representative it is for the full GALAH DR4 domain.
  2. [Figures] Figure captions and axis labels for abundance comparison plots would benefit from inclusion of typical uncertainty values or contours to illustrate the scatter between GALAH and validation measurements.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive and detailed report. The comments highlight important aspects of our threshold calibration and validation strategy. We address each major comment below and have revised the manuscript to incorporate additional discussion and clarifications where feasible.

read point-by-point responses

  1. Referee: [§4] §4 (threshold definition and application): The abundance thresholds are set by comparing the small validation sample to GALAH DR4 values, yet the manuscript provides no quantitative propagation of the differing error distributions, line-blending effects, and parameter covariances between R~28,000 automated analysis and higher-resolution spectra into an estimated false-positive rate when the cuts are applied across the full GALAH parameter space.

    Authors: We agree that a comprehensive propagation of error distributions, line-blending differences, and parameter covariances would strengthen the analysis. In the revised manuscript we have added a dedicated paragraph in §4 that quantifies the typical abundance offsets and scatter between the GALAH DR4 values and our high-resolution validation spectra, and we provide a simple estimate of the implied contamination level based on the observed dispersion in the validation set. A full Monte Carlo simulation that folds in all covariances across the entire GALAH parameter space would require substantial additional computational effort and is beyond the scope of the present work; however, the multi-dimensional nature of the thresholds and the subsequent consistency of the candidate sample with known s-process stars provide supporting evidence that the false-positive rate remains low. revision: partial

  2. Referee: [§5.1–5.2] §5.1–5.2 (post-selection validation): The reported similarity in [hs/ls] ratios is presented as supporting evidence for sample purity, but because this diagnostic is evaluated only on stars already passing the thresholds, it does not independently test the reliability of the threshold transfer or quantify contamination from stars whose abundance errors are correlated with the selection criteria.

    Authors: The referee correctly notes that the [hs/ls] comparison is performed on the selected sample. In the revised version we have expanded §5.1 to include a direct comparison of [hs/ls] distributions for stars lying just below the adopted thresholds, and we discuss the degree to which abundance errors could be correlated with the selection cuts. Because [hs/ls] is a ratio constructed from elements not used in the primary threshold definition and because the thresholds themselves are multi-dimensional, the observed agreement still supplies useful supporting evidence. We acknowledge that this does not constitute a fully independent blind test and have added a brief statement to this effect in the text. revision: partial

Circularity Check

0 steps flagged

No significant circularity: thresholds derived from independent high-resolution validation spectra

full rationale

The paper obtains new UVES and HERMES high-resolution spectra for a validation sample, compares them directly to GALAH DR4 automated abundances, and uses those comparisons to set explicit numerical thresholds on [s/Fe], [Y/Fe], [Zr/Fe], [Ba/Fe], and [La/Fe]. These thresholds are then applied to the full GALAH DR4 catalogue to select the 1073 candidates. The subsequent observation that the selected stars exhibit [hs/ls] ratios similar to confirmed s-rich stars is presented as supporting evidence after selection and is not part of the threshold definition. No equations, fitted parameters, self-citations, or ansatzes are shown to reduce the selection procedure to the paper's own inputs by construction. The derivation chain therefore remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The work rests on standard domain assumptions about s-process nucleosynthesis and binary mass transfer rather than introducing new free parameters or entities; thresholds are empirically set from validation data.

free parameters (1)
  • abundance thresholds on [s/Fe], [Y/Fe], [Zr/Fe], [Ba/Fe], [La/Fe]
    Defined by comparison to the validation sample of high-resolution spectra; specific numerical values are not stated in the abstract.
axioms (1)
  • domain assumption Binary interactions during the AGB phase produce s-process element overabundances that are detectable via specific abundance ratios.
    Invoked throughout the abstract as the physical basis for identifying polluted stars.

pith-pipeline@v0.9.0 · 5751 in / 1291 out tokens · 27026 ms · 2026-05-10T09:51:49.283164+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

2 extracted references · 1 canonical work pages

  1. [1]

    R., Stancliffe, R

    Abate, C., Pols, O. R., Stancliffe, R. J., et al. 2015, A&A, 581, A62 Abdurro’uf, Accetta, K., Aerts, C., et al. 2022, ApJS, 259, 35 Aguado, D. S., Belokurov, V ., Myeong, G. C., et al. 2021, ApJ, 908, L8 Alvarez, R. & Plez, B. 1998, A&A, 330, 1109 Amarsi, A. M., Nordlander, T., Barklem, P. S., et al. 2018, A&A, 615, A139 Bailer-Jones, C. A. L., Rybizki, ...

    work page arXiv 2015

  2. [2]

    2000), the APASS Landolt-Sloan BVgri photometric catalogue (Munari et al

    and the V magnitudes from the Tycho-2 catalogue (Høg et al. 2000), the APASS Landolt-Sloan BVgri photometric catalogue (Munari et al. 2014), or the Fourth US Naval Observatory CCD Astrograph Catalogue (UCAC4; Zacharias et al. 2013). The exposure times, number of spectra (N), and final SNR values, refer to the spectra used (and coadded when N>1) for the sp...

    2000