# Unlocking the mystery of strontium synthesis in the early Galaxy through analysis of barium isotopes in very metal-poor stars

T. Sitnova, L. Lombardo, L. Mashonkina, F. Rizzuti, G. Cescutti, C. Hansen, P. Bonifacio, E. Caffau, A. Koch-Hansen, G. Meynet, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

T. Sitnova, L. Lombardo, L. Mashonkina, F. Rizzuti, G. Cescutti, et al.. Unlocking the mystery of strontium synthesis in the early Galaxy through analysis of barium isotopes in very metal-poor stars. Astronomy & Astrophysics - A&A, 2025, 699, pp.A262. 10.1051/0004-6361/202555073. obspm-05365083

HAL Id: obspm-05365083

https://hal-obspm.ccsd.cnrs.fr/obspm-05365083v1

Submitted on 16 Nov 2025

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Unlocking the mystery of strontium synthesis in the early Galaxy through analysis of barium isotopes in very metal-poor stars*

T. M. Sitnova<fnref n="1" />,*, L. Lombardo<fnref n="2" />, L. I. Mashonkina<fnref n="1" />, F. Rizzuti<fnref n="3" />,4,5, G. Cescutti<fnref n="6" />,4,5, C. J. Hansen<fnref n="2" />, P. Bonifacio<fnref n="7" />, E. Caffau<fnref n="7" />, A. Koch-Hansen<fnref n="8" />, G. Meynet<fnref n="9" />, and R. Fernandes de Melo<fnref n="2" />

Received 8 April 2025 / Accepted 5 June 2025

## ABSTRACT

**Aims.** We determine the contributions of the rapid (r) and slow (s) neutron capture processes to the Ba isotope mixture, along with Ba, Eu, and Sr NLTE abundances, in a sample of very metal-poor stars. The selected stars formed before the contribution from the main s-process in low- and intermediate-mass stars became significant. Some of our sample stars are enhanced in Sr, with [Sr/Ba] reaching up to 0.7. These stars gained their high Sr abundance from a poorly understood process, sometimes referred to in the literature as a light element primary process, which may appear to be a weak s-process or a weak r-process. Our aim is to uncover the nature of this additional Sr source.

**Methods.** The abundances derived from the resonance Ba II 4554 and 4934 Å lines are influenced by the adopted Ba isotope mixture. We computed Ba isotope mixtures corresponding to different r- to s-process contributions (pure r-process, 80%/20%, 50%/50% and 12%/88%, i.e. solar ratio) and determined the corresponding abundances from the Ba II resonance lines in each sample star. Additionally, we determined Ba abundances from weak subordinate Ba II lines, which are unaffected by the adopted Ba isotope mixture. We then compared the Ba abundances derived from the subordinate lines with those from the Ba II resonance lines.

**Results.** We find a higher s-process contribution to Ba isotopes in stars with greater [Sr/Eu] and [Sr/Ba] overabundances, suggesting that the additional Sr synthesis was due to the early s-process occurring in massive stars. Using Sr-enhanced stars, we estimate the [Sr/Ba] ratio produced by the early s-process and obtain $[\text{Sr/Ba}]_{\text{earlyS}} = 1.1 \pm 0.2$. The derived value should be regarded as an upper limit, as we cannot definitively exclude the possibility of a contribution to Sr from the weak r-process, which produces Sr but not Ba. Regarding the potential synthesis of Sr and Ba in the i-process in massive stars, our results for Ba isotopes and element abundances argue that there was no detectable contribution from this process within the error bars in our sample stars.

**Conclusions.** In the early Galaxy, before significant main s-process enrichment, barium and strontium were produced primarily by the main r-process and the early s-process, which occurred in rapidly rotating massive stars.

**Key words.** stars: abundances – stars: atmospheres – Galaxy: evolution – Galaxy: halo

## 1. Introduction

The synthesis of light neutron-capture (n-capture) elements (Sr, Y, and Zr) has been puzzling astrophysicists for several decades. Detailed observational constraints on the origin of heavy elements and the early Galactic chemical evolution can be derived from spectroscopic studies of very metal-poor (VMP, [Fe/H]<fnref n="1" /><−<fnref n="2" />) stars, which, as the oldest objects, provide a unique opportunity to study the earliest epoch of the Galaxy formation.

The literature data report a 2 dex spread in their [Sr/Ba] ratio, ranging from [Sr/Ba] = −0.4 in stars where n-capture elements originate purely from the main rapid (r) n-capture process (Sneden et al. 2003; Barklem et al. 2005; Hayek et al. 2009; Mashonkina et al. 2010; Roederer et al. 2018; Shah et al. 2024), up to a high [Sr/Ba] of 1.5 dex (see, for example, pioneering studies and some more recent work by Spite & Spite 1978; Gratton & Sneden 1994; McWilliam et al. 1995; Mashonkina & Gehren 2001; Aoki et al. 2005; François et al. 2007; Sneden et al. 2008; Roederer 2013; Hansen et al. 2012, 2014, and many others studies). This suggests the presence of an additional source of light n-capture elements beyond the main r-process. The main slow (s) n-capture process does not fit this role, as it occurs in asymptotic giant branch (AGB) stars with masses from 1 to 6 M⊙, and its contribution to the Galactic chemical evolution becomes

<footnote n="1">¹ Institute of Astronomy, Russian Academy of Sciences, Pyatnitskaya 48, 119017 Moscow, Russia</footnote>

<footnote n="2">² Goethe University Frankfurt, Institute for Applied Physics (IAP), Max-von-Laue-Str. 12, 60438 Frankfurt am Main, Germany</footnote>

<footnote n="3">³ Heidelberger Institut für Theoretische Studien, Schloss-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany</footnote>

<footnote n="4">⁴ INAF, Osservatorio Astronomico di Trieste, via Tiepolo 11, 34143 Trieste, Italy</footnote>

<footnote n="5">⁵ INFN, Sezione di Trieste, via Valerio 2, 34134 Trieste, Italy</footnote>

<footnote n="6">⁶ Dipartimento di Fisica, Sezione di Astronomia, Università di Trieste, Via G. B. Tiepolo 11, 34143 Trieste, Italy</footnote>

<footnote n="7">⁷ LIRA, Observatoire de Paris, Université PSL, Sorbonne Université, Université Paris Cité, CY Cergy Paris Université, CNRS, 92190 Meudon, France</footnote>

<footnote n="8">⁸ Zentrum für Astronomie der Universität Heidelberg, Astronomisches Rechen-Institut, Mönchhofstr. 12, 69120 Heidelberg, Germany</footnote>

<footnote n="9">⁹ Department of Astronomy, University of Geneva, Chemin Pegasi 51, 1290 Versoix, Switzerland</footnote>

<footnote>* Based on observations collected at the European Southern Observatory Science Archive and the Keck Observatory Archive.</footnote>

<footnote>** Corresponding author: sitamih@gmail.com</footnote>

<footnote n="1">¹ We use a standard designation, $[\text{X/Y}] = \log(\text{N}_\text{X}/\text{N}_\text{Y})_* - \log(\text{N}_\text{X}/\text{N}_\text{Y})_\odot$, where N_X and N_Y are total number densities of elements X and Y, respectively.</footnote>
significant after a delay of at least 1 Gyr, detectable in stars
with [Fe/H] >-2 (Travaglio et al. 2004; Simmerer et al. 2004;
Käppeler et al. 2011; Kobayashi et al. 2020).

The need for an additional source of light n-capture elements operating in the early epoch of Galactic chemical evolution was first postulated by Travaglio et al. (2004), and it was named the light element primary process (LEPP). However, the LEPP was introduced as a phenomenological process without any specific hypotheses regarding its astrophysical site or nucleosynthesis mechanism. The nature and nucleosynthesis mechanism of this extra source of the light trans-iron elements remain debated. However, it is generally believed to be linked to massive stars, as they evolved rapidly and began to enrich the interstellar gas at the earliest phases of Galactic chemical evolution. Several hypotheses have been proposed for the source associated with the s-process. This includes the weak s-process in massive stars during the hydrostatic core He-burning phase (Raiteri et al. 1991b) and the shell carbon-burning phase (Raiteri et al. 1991a), as well as the s-process in fast-rotating, extremely metal-poor (EMP, [Fe/H] ≤ −3) massive stars (Pignatari et al. 2008; Chiappini et al. 2011; Cescutti et al. 2013; Cescutti & Chiappini 2014; Rizzuti et al. 2021) occurring in the same phases. For observers, a key difference between the above predictions is that models incorporating stellar rotation predict the production of the second-peak s-process elements (Ba, Ce, La, etc.), whereas standard weak s-process models stop at the first-peak s-process elements (Sr, Y, Zr). The s-process models in massive stars have been explored across a wide range of masses, metallicities, and rotational velocities (Frischknecht et al. 2012, 2016; Choplin et al. 2018; Limongi & Chieffi 2018). A key signature of the s-process is the production of even Ba isotopes ¹³⁴Ba and ¹³⁶Ba, which are s-process-only isotopes. Consequently, the chemical evolution model of Cescutti et al. (2013) and Cescutti & Chiappini (2014) predicts that VMP stars with high [Sr/Ba] should exhibit an enhanced value of the even isotopes. This results from their formation from s-process material produced by rapidly rotating massive stars.

Other proposed scenarios include the following: charged-particle reactions in the high-entropy wind of young neutron stars (Woosley & Hoffman 1992; Qian & Wasserburg 2007); the weak r-process in core collapse supernovae (Izutani et al. 2009; Arcones & Bliss 2014); the vp-process in proton-rich supernovae ejecta (Fröhlich et al. 2006; Eichler et al. 2018; Ghosh et al. 2022); the intermediate n-capture process (i-process, Cowan & Rose 1977) in low-metallicity AGB stars with masses below 4 M⊙ (Choplin et al. 2021, 2024); and the i-process in massive VMP stars (Banerjee et al. 2018). These processes, except for the i-process, produce Sr but not Ba, while the i-process yields both Sr and Ba, with a predominance of odd barium isotopes ¹³⁵Ba and ¹³⁷Ba (Choplin et al. 2024). However, the i-process struggles to explain VMP stars with high [Sr/Ba]: the first scenario (Choplin et al. 2024) is attributed to low- and intermediate-mass stars, which enrich the interstellar gas with a delay, while the second one (Banerjee et al. 2018) does not efficiently produce Sr, predicting close to or below solar [Sr/Ba] ratios.

To better understand the source of additional Sr, we aim to
answer (i) whether the source of extra Sr is the s-process or r-
process, (ii) whether this source produces a significant amount of
Ba, and (iii) which Ba isotopes are formed. Our key diagnostics
for identifying the n-capture element sources are: (i) [Ba/Eu] as
a proxy for r-process contribution; (ii) [Sr/Ba] as a tracer of the
extra Sr production source in the earliest epoch; and (iii) the Ba

isotope ratio as a key indicator distinguishing between the s- and
r-process.

The paper is structured as follows. In Sect. 2, we describe our sample stars and observations. The abundance determination method is presented in Sect. 3. A discussion of our findings is provided in Sect. 4, and our conclusions are summarised in Sect. 5.

## 2. Stellar sample and observations

Our sample consists of 16 VMP giant stars from the Chemical Evolution of R-process Elements in Stars (CERES) survey (Lombardo et al. 2022, hereafter LB22). The selected stars have normal carbon abundances (Fernandes de Melo et al. 2024), and their chemical composition is typical to the Milky Way (MW) halo stars with the corresponding metallicity. The majority of the sample stars fall within a metallicity range of -3.10 < [Fe/H] <-2.45, with one star having [Fe/H] = -2.15. Regarding the origin of the n-capture elements in the sample stars, they either have [Ba/H] <-2.4 (Lombardo et al. 2025, hereafter LH25) and formed before the contribution from the main s-process in low- and intermediate-mass stars became significant (Simmerer et al. 2004, LH25), or they are r-process enhanced stars with -2.4 < [Ba/H] <-2.2. Three of the sample stars are strongly r-process-enhanced r-II type stars with [Eu/Fe]<fnref n="2" /> = 1.3, while seven stars are r-I type, with 0.3 < [Eu/Fe] < 0.8. Thus, the selected stars trace the earliest epoch of Galactic chemical evolution and n-capture element synthesis.

For Ba isotope ratio analysis, we selected stars with high-resolution observed spectra covering both Ba II resonance lines at 4554 Å and 4934 Å. Additionally, the equivalent widths (EWs) of these lines should not exceed 140 mÅ. Sixteen of the fifty-two CERES stars meet these criteria, with their Ba II resonance line EWs ranging from 42 mÅ to 135 mÅ. The selected stars, along with their stellar atmosphere parameters, are listed in Table 1. Our selection criteria provide a suitable stellar sample, where the Ba II resonance lines are strong enough to trace Ba isotope ratios while remaining unsaturated enough to allow for accurate abundance measurements.

For our sample stars, high-resolution spectra were obtained by LB22 with the UV-Visual Echelle Spectrograph (UVES) at the UT2 Kueyen Telescope, with high signal-to-noise (S/N) ratios >150 per pixel at 3900 Å. Details of the individual LB22 spectra, including S/N ratio, spectral resolution, exposure time, observation date, and settings, are provided in Table A.1 of LB22. The LB22 spectra cover a wide wavelength range and allow us to analyse the lines of Ba II at 5853 Å, 6141 Å, 6494 Å, and 4934 Å, but not the Ba II 4554 Å line, due to a wavelength gap between 4520 Å and 4780 Å. Therefore, in this study, we complement the observations from LB22 with additional spectra covering the Ba II 4554 Å line, obtained either with UVES or the High-Resolution Echelle Spectrometer (HIRES) at the Keck I Telescope. The adopted spectra are of high quality, with a spectral resolution of $\lambda/\Delta\lambda > 40\,000$ and a signal-to-noise ratio of S/N > 100. These spectra are available from the European Southern Observatory Science Archive and the Keck Observatory Archive, respectively. The proposal IDs and PIs for the spectra used in the Ba II 4554 Å analysis are listed in Table 1.

<footnote n="2">² Here, we employ the notation of Christlieb et al. (2004) for the r-process enhanced star classification.</footnote>
<table><caption>**Table 1.** Stellar sample, atmospheric parameters, and characteristics of the observed spectra used for the Ba II 4554 Å line analysis.</caption><thead><tr><th>Name<sub>SIMBAD</sub></th><th>Name<sub>CERES</sub></th><th>T<sub>eff</sub><sup>1</sup>, K</th><th>log g cm s<sup>-2</sup></th><th>[Fe/H]</th><th>ξ<sub>t</sub>, km s<sup>-1</sup></th><th>S/T<sup>2</sup></th><th>Program ID; PI</th></tr></thead><tbody><tr><td>HD 4306</td><td>CES0045-0932</td><td>5020</td><td>2.29</td><td>-2.95</td><td>1.60</td><td>H</td><td>C36H; S. Castro</td></tr><tr><td>HD 13979</td><td>CES0215-2554</td><td>5080</td><td>2.00</td><td>-2.73</td><td>1.90</td><td>H</td><td>U09H; R. Kraft</td></tr><tr><td>HD 27928</td><td>CES0422-3715</td><td>5100</td><td>2.46</td><td>-2.45</td><td>1.60</td><td>H</td><td>C23H; I. Ivans</td></tr><tr><td>HD 107752</td><td>CES1222+1136</td><td>4830</td><td>1.72</td><td>-2.91</td><td>1.75</td><td>H</td><td>C04H; J. Cohen</td></tr><tr><td>HD 108317</td><td>CES1226+0518</td><td>5340</td><td>2.84</td><td>-2.38</td><td>1.50</td><td>U</td><td>165.N-0276(A); R. Cayrel</td></tr><tr><td>HD 122563</td><td>CES1402+0941</td><td>4680</td><td>1.35</td><td>-2.79</td><td>2.00</td><td>U</td><td>266.D-5655(A); S. Bagnulo</td></tr><tr><td>HD 126587</td><td>CES1427-2214</td><td>4910</td><td>1.99</td><td>-3.05</td><td>1.60</td><td>H</td><td>U65H; M. Bolte</td></tr><tr><td>HD 128279</td><td>CES1436-2906</td><td>5280</td><td>3.15</td><td>-2.15</td><td>1.30</td><td>U</td><td>71.B-0529(A); D. Silva</td></tr><tr><td>HE0336-2412</td><td>CES0338-2402</td><td>5240</td><td>2.78</td><td>-2.81</td><td>1.50</td><td>H</td><td>C2230; M. Brown</td></tr><tr><td>HE1320-1339</td><td>CES1322-1355</td><td>4960</td><td>1.81</td><td>-2.93</td><td>1.90</td><td>U</td><td>170.D-0010(G); N. Christlieb</td></tr><tr><td>HE2229-4153</td><td>CES2232-4138</td><td>5190</td><td>2.76</td><td>-2.58</td><td>1.50</td><td>U</td><td>170.D-0010(G); N. Christlieb</td></tr><tr><td>HE2327-5642</td><td>CES2330-5626</td><td>5030</td><td>2.31</td><td>-3.10</td><td>1.55</td><td>U</td><td>170.D-0010(G); N. Christlieb</td></tr><tr><td>BD+23 3130</td><td>CES1732+2344</td><td>5370</td><td>2.82</td><td>-2.57</td><td>1.50</td><td>H</td><td>U080Hb; J. Prochaska</td></tr><tr><td>BPS CS31078-018</td><td>CES0301+0616</td><td>5220</td><td>3.01</td><td>-2.93</td><td>1.30</td><td>H</td><td>U10H; M. Bolte</td></tr><tr><td>BPS CS29491-069</td><td>CES2231-3238</td><td>5220</td><td>2.67</td><td>-2.77</td><td>1.50</td><td>U</td><td>170.D-0010(G); N. Christlieb</td></tr><tr><td>BPS CS22186-023</td><td>CES0419-3651</td><td>5090</td><td>2.29</td><td>-2.81</td><td>1.60</td><td>U</td><td>68.B-0320(A); N. Christlieb</td></tr></tbody></table>

Notes. ¹The effective temperatures from LB22 have been rounded to the nearest 10 K. ²Spectrograph/Telescope: U – UVES/UT2, H – HIRES/Keck I.

# 3. Abundance analysis

## 3.1. Codes and model atmospheres

In this study, we used classical 1D model atmospheres from the MARCS model grid (Gustafsson et al. 2008), interpolated for the given $T_{eff}$, log g, and [Fe/H] of the stars, while ATLAS12 model atmospheres (Kurucz 2005) were adopted in the CERES project (LB22 and LH25). We verified that using either MARCS or ATLAS12 models yields consistent results.

We solved the coupled radiative transfer and statistical equilibrium equations using the DETAIL code (Giddings 1981; Butler 1984), incorporating the updated opacity package presented by Mashonkina et al. (2011). For synthetic spectra calculations, we used the SYNTHV_NLTE code (Tsymbal et al. 2019), attached to the IDL BINMAG code (Kochukhov 2018). This method allowed us to obtain the best fit to the observed line profiles while accounting for the non-local thermodynamic equilibrium (NLTE) effects via pre-calculated departure coefficients (the ratio between NLTE and LTE atomic level populations) for a given model atmosphere. When fitting the line profiles, the abundance of the element of interest was varied alongside the macroturbulent velocity ($v_{mac}$) and the radial velocity ($v_r$).

The spectral synthesis line list was extracted from a recent version of the Vienna Atomic Line Database (VALD, Pakhomov et al. 2019; Ryabchikova et al. 2015), which provides isotopic and hyperfine structure components for a number of studied chemical elements. We adopted the oscillator strengths recommended by VALD.

## 3.2. Atmospheric parameters and impact of their uncertainties on abundance determination

The selected sample stars are VMP giants with parallaxes ranging from 0.17 to 7.63 mas (Gaia Collaboration 2021). Photometric effective temperatures ($T_{eff}$) and distance-based surface gravities (log g) are adopted from LB22. The uncertainty in $T_{eff}$ is 100 K, and the uncertainty in log g amounts to 0.08 dex. In this study, we focused on the analysis of Ba II lines and the

comparison of barium abundances derived from the resonance (strong) and subordinate (weak) lines.

The uncertainties in $T_{eff}$ and log g affect abundances from different Ba II lines in the same way and do not significantly impact the abundance difference between the resonance and subordinate lines of Ba II. For example, in a model atmosphere with $T_{eff}/\log g$-[Fe/H] = 5100/2.46/-2.5, a shift in $T_{eff}$ of 100 K results in an abundance shift of 0.10 and 0.07 dex for the Ba II 4554 Å and 5853 Å lines, respectively, with EWs of 120 mÅ and 25 mÅ. Overall, this shift in $T_{eff}$ produces a minor effect of 0.03 dex on the abundance difference between the two lines. In the same model atmosphere, a change in log g of 0.08 dex results in a similar abundance shift of 0.03 dex for both lines and does not affect the abundance difference.

Regarding variations in the microturbulent velocity ($\xi_t$), an uncertainty of 0.2 km s$^{-1}$ has a significant effect on the abundance derived from the resonance line, causing a shift of 0.10 dex, while the corresponding shift for the Ba II 5853 Å line is negligible. In our sample stars, the EWs of the Ba II resonance lines vary from 60 to 130 mÅ. For the weaker Ba II resonance line, with EW$_{4554}$ of 60 mÅ, an uncertainty in $\xi_t$ of 0.2 km s$^{-1}$ results in an abundance shift of 0.04 dex.

## 3.3. Determination of $\xi_t$

LB22 relies on $\xi_t$ calculated using an empirical formula provided in Mashonkina et al. (2017a, hereafter MJ17). The adopted formula was derived from $\xi_t$ determinations based on Fe I and Fe II lines in a sample of VMP giants and provides an accuracy of 0.2 km s$^{-1}$. In our study, the uncertainty in $\xi_t$ is crucial, since it mostly contributes to an uncertainty in barium abundance from the resonance lines. We refine $\xi_t$ using Ti II lines and adopt $\xi_t$ values that yield consistent titanium abundances from Ti II lines with different EWs.

Unlike the conventional use of iron lines for $\xi_t$ determinations, we selected Ti II lines because they have accurate oscillator strengths based on the laboratory measurements of Wood et al. (2013), exhibit minimal NLTE effects, and are convenient to use.
<table><caption>**Table 2.** NLTE abundance corrections and EWs in mÅ for Ti II lines as a function of $T_{\text{eff}}$, log g, and [Fe/H].</caption><thead><tr><th colspan="8">λ, Å species E<sub>exc</sub>, eV log gf</th></tr><tr><th>T<sub>eff</sub>, K<br>[Fe/H]<sub>1</sub><br>[Fe/H]<sub>1</sub></th><th>log g<sub>1</sub><br>EW<sub>1</sub><br>Δ<sub>1</sub></th><th>log g<sub>2</sub><br>EW<sub>2</sub><br>Δ<sub>2</sub></th><th>log g<sub>3</sub><br>EW<sub>3</sub><br>Δ<sub>3</sub></th><th>log g<sub>4</sub><br>EW<sub>4</sub><br>Δ<sub>4</sub></th><th>log g<sub>5</sub><br>EW<sub>5</sub><br>Δ<sub>5</sub></th><th>log g<sub>6</sub><br>EW<sub>6</sub><br>Δ<sub>6</sub></th></tr></thead><tbody><tr><td>4417.710, Ti II, E<sub>exc</sub> = 1.16, log gf = -1.19</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>T<sub>eff</sub> = 5000</td><td>1.0</td><td>1.5</td><td>2.0</td><td>2.5</td><td>3.0</td><td>3.5</td></tr><tr><td>-2.0</td><td>-1</td><td>124</td><td>106</td><td>92</td><td>79</td><td>-1</td></tr><tr><td>-2.0</td><td>-1</td><td>-0.094</td><td>-0.085</td><td>-0.065</td><td>-0.049</td><td>-1</td></tr></tbody></table>

**Notes.** The table is accessible in a machine-readable format at the CDS. A portion is shown to illustrate its format and content. If EW = -1 and $\Delta_{\text{NLTE}} = -1$, this means that EW is either <3 mÅ or not computed in a model atmosphere with given parameters.

The strongest lines of Ti II are of comparable strength to the Ba II resonance lines in the corresponding model atmospheres. For example, in the case of the Ba II 4554 Å line, the equivalent widths are EW<sub>4554</sub> = 60 mÅ and 124 mÅ, for HE0336–2412 and HD 27928, respectively. For the Ti II 4468 Å line, the equivalent widths are EW<sub>4417</sub> = 85 mÅ and 109 mÅ in the same two stars, respectively.

Although the NLTE effects for Ti II are not large (Bergemann 2011; Sitnova et al. 2016; Mallinson et al. 2022), they are not negligible. We account for deviations from LTE using the Ti I–II model atom from Sitnova et al. (2016), updated with data for inelastic collisions with hydrogen atoms (Ti II + H) based on quantum-mechanical calculations presented in Sitnova et al. (2020). A grid of NLTE abundance corrections (the difference between NLTE and LTE abundance) for Ti II lines covering our stellar parameter range is provided in Table 2. For each sample star, we select the microturbulent velocity that provides consistent abundances from weak and strong Ti II lines, with abundance slopes as a function of EW of $\le 10^{-4}$ dex per mÅ.

When comparing the derived $\xi_t$ values with those calculated from the empirical formula of MJ17, we find slightly lower values, with the difference never exceeding $0.25 \text{ km s}^{-1}$. For CES1427–2214, the star with the largest revision in $\xi_t$ (0.25 km s$^{-1}$), we show the NLTE abundances from Ti II as a function of EW for $\xi_t$ calculated using the MJ17 formula and those derived from Ti II lines (Fig. 1). Using $\xi_t = 1.85 \text{ km s}^{-1}$ results in a significant trend in both NLTE and LTE, with slopes of $-2.0 \cdot 10^{-3}$ and $-1.3 \cdot 10^{-3}$ dex per mÅ, respectively. In contrast, a $\xi_t$ value lower by $0.25 \text{ km s}^{-1}$ results in consistent abundances from Ti II lines of different strengths, with slopes of $-5 \cdot 10^{-5}$ and $6 \cdot 10^{-4}$ dex per mÅ in NLTE and LTE, respectively. For CES1427–2214, we ensured that the $\xi_t$ derived from Ti II lines yielded consistent NLTE abundances from Fe I and Fe II lines of different strength (Fig. 1, bottom panel). See MJ17 for the line list and the NLTE method of iron abundance determination.

Observations show that stars with higher log g have lower $\xi_t$ (see, for example, the empirical formula for $\xi_t$ derived for dwarfs and subgiants in Sitnova et al. 2015). The majority of our sample stars are giants at the base of the RGB, with log g $\approx 2.5$, whereas the stellar sample of MJ17 consists mostly of cool giants, with 45 out of 59 stars having log g less than 2.0. The difference in stellar parameters between our sample and the stars used to derive the empirical formula may explain the slightly lower $\xi_t$ values obtained in this study. We find systematically lower $\xi_t$ values in

<figure><img src="image_3.png" /><figcaption>**Fig. 1.** Top panel: NLTE abundances from Ti II lines in CES1427–2214 as a function of the EW computed for $\xi_t = 1.60 \text{ km s}^{-1}$ (squares) and $1.85 \text{ km s}^{-1}$ (circles). Bottom panel: NLTE abundances from Fe I (triangles) and Fe II (asterisks) lines in CES1427–2214 computed for $\xi_t = 1.60 \text{ km s}^{-1}$.</figcaption></figure>

stars with higher log g, whereas consistent $\xi_t$ values are obtained for the four stars with log g <2.0.

Using our new measurements, we derived an empirical formula suitable for $\xi_t$ determination in giants at the base of the RGB: $\xi_t = -1.12 + 8.37 \cdot 10^{-4} \cdot T_{\text{eff}} - 6.48 \cdot 10^{-1} \cdot \log g$.

Here, we assume that $\xi_t$ does not depend on [Fe/H], as the corresponding coefficients in the MJ17 and Sitnova et al. (2015) formulae are small and have opposite signs, suggesting no significant correlation between $\xi_t$ and [Fe/H]. When comparing $\xi_t$ measurements for individual stars with predictions from the above formula, we find that the difference does not exceed $0.10 \text{ km s}^{-1}$ in absolute value, with a standard deviation of $0.05 \text{ km s}^{-1}$.

### 3.4. Analysis of the Ba II lines

We analysed five lines of Ba II: the strong resonance lines at $\lambda = 4554$ Å and 4934 Å, and the subordinate lines at $\lambda = 5853$ Å, 6141 Å, and 6496 Å. There is a Fe I 4934.0 Å line with a lower level energy of $E_{\text{exc}} = 4.15 \text{ eV}$ and an oscillator strength of log gf = -0.58 in the blue wing of the Ba II 4934.1 Å line. Fifteen of our sample stars have [Fe/H] $\le -2.4$, and blending with the Fe I 4934.0 Å line is negligible. For the most metal-rich sample star, HD 128279 (CES1436–2906), with [Fe/H] = -2.15, we show the observed spectrum around the Ba II 4934.1 Å line, along with synthetic spectra calculated with and without the iron line (Fig. 2).

The subordinate line Ba II 6141.71 Å is blended with the Fe I 6141.73 Å line, which has $E_{\text{exc}} = 3.60 \text{ eV}$ and log gf = -1.46. We accounted for this line in our spectrum synthesis using the iron abundance derived in LB22. For the most metal-rich star, varying the iron abundance within an uncertainty of 0.13 dex (LB22) results in negligible changes in the synthetic spectrum and does not affect barium abundance determination.

#### 3.4.1. NLTE effects

To account for departures from LTE, we adopted the Ba II model atom from Mashonkina & Belyaev (2019). The NLTE effects for
<figure><img src="image_4.png" /><figcaption>**Fig. 2.** Ba II 4934 Å line profile in the observed spectrum (circles) of the most metal-rich sample star CES1436–2906 with [Fe/H] = −2.15 along with a synthetic best-fit NLTE spectrum (solid line) and a spectrum, calculated neglecting the Fe I 4934 Å line (dashed line).</figcaption></figure>

Ba II are moderate, leading to either weaker or stronger lines compared to LTE, depending on the line strength and stellar atmosphere parameters. In our sample stars, NLTE abundance corrections for individual Ba II lines do not exceed 0.13 dex in absolute value.

We draw attention to an intriguing line formation phenomenon observed for the Ba II resonance line in some stars. Based on the ratio of the oscillator strengths of the Ba II 4934 Å and Ba II 4554 Å lines ($f_{4934}/f_{4554} = 0.60$), one might expect a similar ratio for their corresponding EWs, assuming the spectral lines lie within the linear portion of the curve of growth. Although our lines are stronger and extend beyond the linear portion, in the three stars with the weakest Ba II lines, the ratio EW$_{4934}$/EW$_{4554}$ is approximately 0.7. Notably, this ratio increases for stronger Ba II lines, though it never reaches unity in LTE. In seven of our sample stars with strong resonance lines (EW > 110 mÅ), the Ba II 4554 Å and Ba II 4934 Å lines have nearly the same strength, and in four of them, the EW of the Ba II 4934 Å line slightly exceeds that of the Ba II 4554 Å line.

When Ba II resonance lines are stronger than ≈100 mÅ, NLTE leads to stronger lines compared to LTE, and the NLTE abundance corrections are negative. For the Ba II 4934 Å line, NLTE corrections are larger in absolute value than for the Ba II 4554 Å line. This paradox is explained by the fact that the NLTE effects for the Ba II 4934 Å line result in more efficient line strengthening compared to the Ba II 4554 Å line. Figure 3 shows the departure coefficients for the ground state of Ba II (6s) and the upper levels of the resonance transitions in a model atmosphere of a VMP giant. The Ba II 4554 Å and 4934 Å lines form in the 6s–6p (J = 3/2) and 6s–6p (J = 1/2) transitions, respectively.

Strong resonance lines form across a broad atmospheric depth range, up to $log\tau_{5000} = -2.1$. In the high atmospheric layers, the upper levels are underpopulated relative to LTE and to the ground state, and the ratio $S_v/B_v \sim b_{up}/b_{low} < 1$, where $S_v$ and $B_v$ are the line source function and the Planck function, respectively. A decrease in the line source function results in a lower flux and a stronger line core compared to LTE. In deeper atmospheric layers with $log\tau_{5000} > -1$, the situation is reversed: radiative pumping in the resonance transition results in a relative overpopulation of the upper levels relative to the ground state, and weakening the line compared to LTE. This effect is smaller for the Ba II 4934 Å line, which has a lower oscillator strength than the Ba II 4554 Å line. Due to the reduced line weakening

<figure><img src="image_5.png" /><figcaption>**Fig. 3.** Departure coefficients of the ground state of Ba II (6s) and the upper levels of the resonance transitions in a model atmosphere of a VMP giant. The vertical dashes indicate the line core formation depths of the resonance lines of different strength: 130 mÅ and 60 mÅ at $log\tau_{5000} = -2.1$ and $-1.2$, respectively.</figcaption></figure>

in the deep layers, the Ba II 4934 Å line may ultimately appear stronger than the Ba II 4554 Å line.

### 3.4.2. Impact of Ba isotope ratio on abundance determination

Barium abundance in stars can be represented by up to five isotopes: $^{134}$Ba, $^{135}$Ba, $^{136}$Ba, $^{137}$Ba, and the most abundant isotope, $^{138}$Ba. The isotopes $^{134}$Ba and $^{136}$Ba are s-only isotopes and cannot be produced in the r-process due to the presence of the stable r-only isotopes $^{134}$Xe and $^{136}$Xe, which block the $\beta$-decay path and prevent the formation of $^{134}$Ba and $^{136}$Ba (Rolfs et al. 1987).

Odd isotopes are subject to hyperfine splitting (HFS) of the energy levels, and a higher fraction of odd isotopes results in a broader line profile and greater total absorbed energy. The HFS primarily affects the ground state, meaning that Ba II resonance lines can serve as a diagnostic of the Ba isotope ratio. In contrast, subordinate lines are unaffected by the adopted isotope ratio and can be used as reliable indicators of barium abundance. By determining the abundance from subordinate lines, one can infer the Ba isotope ratio, which can then be translated into the relative contributions of the r- and s-processes. The idea of estimating the r/s mixture of a star was first proposed by Cowley & Frey (1989), who pointed out the importance of accounting for HFS in Ba abundance determinations. This method has been applied to MP stars in the MW (Mashonkina & Zhao 2006; Mashonkina et al. 2008; Mashonkina & Belyaev 2019) and the Sculptor dwarf spheroidal (dSph) galaxy (Jablonka et al. 2015).

Magain & Zhao (1993) and Magain (1995) suggested a different method based on measuring the broadening of the Ba II resonance lines and applied it to the VMP star HD 140283. Depending on the quality of the observed spectrum and the adopted line formation calculation method, different authors report different results for Ba isotope origin in HD 140283: a pure r-process in 1D LTE (Magain & Zhao 1993; Magain 1995), a solar mixture in 1D LTE (Lambert & Allende Prieto 2002) and in 3D LTE (Collet et al. 2009), a pure s-process in 1D LTE (Gallagher et al. 2010), and a mostly r-process in 3D LTE (Gallagher et al. 2015).

For barium, we determined abundances using different isotope ratios corresponding to the following r- to s-process fractions: 100%/0% (pure r-process), 80%/20%, 50%/50%, 12%/88% (solar), and 0%/100% (pure s-process). We adopted s-process and r-process residual Ba isotope yields from Prantzos et al. (2020). A larger r-process fraction results in a higher proportion of odd isotopes and a lower abundance derived from the resonance lines compared to that obtained with a lower r-process fraction. For each star, we provided abundances from individual spectral lines along with their measured EWs (Table 3).
<table><caption>**Table 3.** NLTE and LTE abundances and EWs of Ba II lines in the sample stars.</caption><thead><tr><th rowspan="2">Ba II lines:</th><th colspan="5">4554 Å</th><th colspan="5">4934 Å</th><th rowspan="2">5853 Å</th><th rowspan="2">6141 Å</th><th rowspan="2">6496 Å</th></tr><tr><th>100/0</th><th>80/20</th><th>50/50</th><th>15/85</th><th>0/100</th><th>100/0</th><th>80/20</th><th>50/50</th><th>12/88</th><th>0/100</th></tr></thead><tbody><tr><td>CES1322–1355</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>log ε<sub>NLTE</sub></td><td>-1.06</td><td>-1.02</td><td>-0.94</td><td>-0.76</td><td>-0.64</td><td>-1.11</td><td>-1.07</td><td>-1.00</td><td>-0.83</td><td>-0.72</td><td>-0.95</td><td>-0.96</td><td>-0.95</td></tr><tr><td>log ε<sub>LTE</sub></td><td>-1.07</td><td>-1.03</td><td>-0.95</td><td>-0.77</td><td>-0.65</td><td>-1.13</td><td>-1.09</td><td>-1.02</td><td>-0.85</td><td>-0.74</td><td>-1.02</td><td>-1.04</td><td>-0.99</td></tr><tr><td>EW, mÅ</td><td>111.0</td><td>111.0</td><td>111.0</td><td>111.0</td><td>111.0</td><td>102.1</td><td>102.1</td><td>102.1</td><td>102.1</td><td>102.1</td><td>11.2</td><td>47.0</td><td>40.0</td></tr></tbody></table>

**Notes.** The table is accessible in a machine-readable format at the CDS. A portion is shown to illustrate its format and content.

<figure><img src="image_6.png" /><figcaption>**Fig. 4.** NLTE abundances from the Ba II 4554 Å and 4934 Å lines as a function of r-process contribution to Ba isotopes in the two sample stars. The EWs of the resonance lines are indicated.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Fig. 5.** Top panel: Ba II 4554 Å line profile in the observed spectrum (circles) of HE0336–2412 (CES0338–2402). Synthetic best-fit NLTE spectra, derived using pure r-process (dashed line) and pure s-process (dotted line) Ba isotope mixtures are almost indistinguishable, demonstrating that a reasonable fit can be achieved with any adopted isotope ratio. Vertical dashes show the relative contribution of different Ba isotopes to the r-process and s-process; see legend for designations. Bottom panel: Same as the top panel for HD 27928 (CES0422–3715).</figcaption></figure>

Figure 4 illustrates the impact of the adopted Ba isotope mixture on barium abundance derived from the resonance lines. We plot the abundances from Ba II 4554 Å and 4934 Å lines in two sample stars: one with the weakest lines (HE0336–2412/CES0338–2402) and another with the strongest lines (r-II type star HE2327–5642/CES2330–5626). The stronger the Ba II resonance lines, the more sensitive their abundances are to the adopted Ba isotope mixture. For example, for the Ba II 4554 Å line with EW = 59 mÅ, the abundance difference between pure r- and pure s-process Ba isotope mixtures is 0.30 dex. The corresponding difference in HE2327–5642, where EW<sub>4554</sub> = 133 mÅ, amounts to 0.34 dex. Notably, in HE2327–5642, the Ba II 4934 Å line is slightly more sensitive to the Ba isotope mixture, as it is slightly stronger than the Ba II 4554 Å line due to the effect discussed earlier.

<table><caption>**Table 4.** NLTE and LTE abundances and EWs from individual Sr II and Eu II lines in the sample stars.</caption><thead><tr><th rowspan="2">λ, Å:</th><th colspan="2">Sr II</th><th colspan="3">Eu II</th></tr><tr><th>4077</th><th>4215</th><th>3819</th><th>4129</th><th>4205</th></tr></thead><tbody><tr><td>CES1732+2344</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>log ε<sub>NLTE</sub></td><td>0.14</td><td>0.14</td><td>-1.73</td><td>-1.55</td><td>-1.59</td></tr><tr><td>log ε<sub>LTE</sub></td><td>0.24</td><td>0.26</td><td>-1.86</td><td>-1.67</td><td>-1.70</td></tr><tr><td>EW, mÅ</td><td>107.9</td><td>97.9</td><td>12.6</td><td>9.4</td><td>9.3</td></tr></tbody></table>

**Notes.** The table is accessible in a machine-readable format at the CDS. A portion is shown to illustrate its format and content.

The adopted Ba isotope mixture has a significant impact on the abundance derived from resonance lines, while the observed line profiles can be fitted when adopting any Ba isotope mixture. To demonstrate the quality of spectral fitting, we selected the UVES spectrum of HE0336–2412 (CES0338–2402), which has a spectral resolution of R = 107 200, higher than that of most sample stars, whose spectra were observed at R ≈45 000. Since HE0336–2412 has one of the weakest Ba II 4554 Å lines (EW = 59 mÅ), we also include the spectrum of HD 27928 (CES0422–3715) with EW<sub>4554</sub> = 124 mÅ, obtained using HIRES with R = 47 700 (Fig. 5). For illustration, we present the best-fit synthetic NLTE spectra of the Ba II 4554 Å line computed with pure r-process and pure s-process Ba isotope ratios. Regardless

of whether an r- or s-process Ba isotope mixture is used, we achieve a satisfactory fit to the observed spectra.

For CES0338–2402, the best-fit spectrum parameters are log ε(Ba)/v<sub>mac</sub> = -1.24/3.6 km s<sup>-1</sup> (pure r-process) and
<table><caption>**Table 5.** NLTE abundance ratios and Ba odd isotope fractions of sample stars.</caption><thead><tr><th>Name</th><th>[Sr/H]</th><th>[Ba/H]</th><th>[Eu/H]</th><th>[Sr/Ba]</th><th>[Ba/Eu]</th><th>[Sr/Eu]</th><th>[Eu/Fe]</th><th>F<sub>odd</sub></th></tr></thead><tbody><tr><td>CES0419-3651</td><td>-2.97 (0.13)</td><td>-3.45 (0.09)</td><td>&lt;-3.09 (0.08)</td><td>0.48 (0.11)</td><td>&gt;-0.36 (0.06)</td><td>&gt;0.12 (0.09)</td><td>&lt;-0.30</td><td>0.15<sup>+0.15</sup><sub>-0.05</sub></td></tr><tr><td>CES1402+0941</td><td>-2.83 (0.14)</td><td>-3.58 (0.07)</td><td>-3.02 (0.08)</td><td>0.74 (0.11)</td><td>-0.55 (0.02)</td><td>0.19 (0.11)</td><td>-0.25</td><td>0.14<sup>+0.09</sup><sub>-0.04</sub></td></tr><tr><td>CES0045-0932</td><td>-3.05 (0.13)</td><td>-3.56 (0.07)</td><td>&lt;-2.83 (0.08)</td><td>0.51 (0.09)</td><td>&gt;-0.73 (0.04)</td><td>&gt;-0.21 (0.10)</td><td>&lt;0.10</td><td>0.11<sup>+0.12</sup><sub>-0.01</sub></td></tr><tr><td>CES1436-2906</td><td>-2.74 (0.13)</td><td>-2.58 (0.08)</td><td>-2.02 (0.08)</td><td>-0.15 (0.11)</td><td>-0.56 (0.06)</td><td>-0.71 (0.10)</td><td>0.11</td><td>0.23<sup>+0.13</sup><sub>-0.08</sub></td></tr><tr><td>CES0338-2402</td><td>-3.03 (0.14)</td><td>-3.25 (0.09)</td><td>&lt;-2.65 (0.08)</td><td>0.22 (0.12)</td><td>&gt;-0.60 (0.06)</td><td>&gt;-0.38 (0.10)</td><td>&lt;0.14</td><td>0.20<sup>+0.18</sup><sub>-0.08</sub></td></tr><tr><td>CES0215-2554</td><td>-3.03 (0.14)</td><td>-3.06 (0.07)</td><td>-2.40 (0.08)</td><td>0.03 (0.10)</td><td>-0.67 (0.02)</td><td>-0.64 (0.10)</td><td>0.31</td><td>0.23<sup>+0.13</sup><sub>-0.07</sub></td></tr><tr><td>CES1322-1355</td><td>-2.62 (0.14)</td><td>-3.12 (0.07)</td><td>-2.45 (0.08)</td><td>0.51 (0.11)</td><td>-0.67 (0.01)</td><td>-0.16 (0.11)</td><td>0.46</td><td>0.39<sup>+0.20</sup><sub>-0.13</sub></td></tr><tr><td>CES1732+2344</td><td>-2.74 (0.13)</td><td>-2.87 (0.08)</td><td>-2.08 (0.08)</td><td>0.13 (0.10)</td><td>-0.79 (0.05)</td><td>-0.66 (0.09)</td><td>0.47</td><td>0.26<sup>+0.27</sup><sub>-0.13</sub></td></tr><tr><td>CES2232-4138</td><td>-2.47 (0.13)</td><td>-2.65 (0.07)</td><td>-1.92 (0.09)</td><td>0.18 (0.10)</td><td>-0.73 (0.04)</td><td>-0.55 (0.10)</td><td>0.64</td><td>0.32<sup>+0.16</sup><sub>-0.10</sub></td></tr><tr><td>CES1222+1136</td><td>-2.81 (0.14)</td><td>-2.89 (0.07)</td><td>-2.23 (0.08)</td><td>0.07 (0.10)</td><td>-0.66 (0.02)</td><td>-0.58 (0.10)</td><td>0.66</td><td>0.69<sup>+0.06</sup><sub>-0.23</sub></td></tr><tr><td>CES1427-2214</td><td>-2.92 (0.13)</td><td>-2.93 (0.07)</td><td>-2.29 (0.08)</td><td>0.01 (0.10)</td><td>-0.63 (0.03)</td><td>-0.62 (0.10)</td><td>0.74</td><td>0.58<sup>+0.17</sup><sub>-0.17</sub></td></tr><tr><td>CES0422-3715</td><td>-2.59 (0.13)</td><td>-2.42 (0.08)</td><td>-1.67 (0.08)</td><td>-0.18 (0.10)</td><td>-0.75 (0.04)</td><td>-0.92 (0.10)</td><td>0.76</td><td>0.57<sup>+0.18</sup><sub>-0.18</sub></td></tr><tr><td>CES1226+0518</td><td>-2.41 (0.14)</td><td>-2.24 (0.07)</td><td>-1.57 (0.08)</td><td>-0.18 (0.10)</td><td>-0.67 (0.02)</td><td>-0.84 (0.10)</td><td>0.79</td><td>0.46<sup>+0.24</sup><sub>-0.15</sub></td></tr><tr><td>CES2330-5626</td><td>-2.78 (0.13)</td><td>-2.49 (0.07)</td><td>-1.70 (0.08)</td><td>-0.30 (0.09)</td><td>-0.79 (0.01)</td><td>-1.08 (0.09)</td><td>1.38</td><td>0.75</td></tr><tr><td>CES0301+0616</td><td>-2.63 (0.13)</td><td>-2.33 (0.07)</td><td>-1.57 (0.08)</td><td>-0.30 (0.10)</td><td>-0.77 (0.03)</td><td>-1.06 (0.10)</td><td>1.34</td><td>0.75<sup>-0.11</sup></td></tr><tr><td>CES2231-3238</td><td>-2.65 (0.15)</td><td>-2.29 (0.08)</td><td>-1.45 (0.08)</td><td>-0.35 (0.12)</td><td>-0.84 (0.05)</td><td>-1.19 (0.11)</td><td>1.30</td><td>0.75<sup>-0.06</sup></td></tr></tbody></table>

Notes. Solar abundances are taken from Lodders (2021). Total uncertainties are given in parenthesis.

-0.94/4.2 km s⁻¹ (pure s-process). The corresponding values for CES0422-3715 are log ε(Ba)/v<sub>mac</sub> = -0.30/3.0 km s⁻¹ and 0.07/3.3 km s⁻¹. In our abundance determination method, v<sub>mac</sub> is treated as a free parameter when fitting individual spectral lines. From analysis of different spectral lines in a given star, we estimate a typical uncertainty in v<sub>mac</sub> as 2 km s⁻¹. This uncertainty prevents us from determining the Ba isotope mixture solely from the Ba II resonance line profiles.

### 3.5. Sr and Eu abundances

We determined Sr abundances from the Sr II 4077 Å and 4215 Å lines (Table 4.) To account for the NLTE effects, we relied on the Sr II model atom presented in Mashonkina et al. (2022) and the NLTE abundance corrections from the INASAN database<fnref n="3" /> (Mashonkina et al. 2023). There are two reasons why we determined Sr abundances in this study instead of taking them from LB22. First, for six of the selected sample stars, Sr abundances are missing in LB22. Secondly, a revision in ξᵢ affects the Sr II lines, as they are strong in our sample stars, and Sr abundances might be up to 0.2 dex higher than those computed with the ξᵢ from the MJ17 formula.

We determined Eu abundances in our sample stars using the Eu II model atom of Mashonkina & Gehren (2000). LH25 relies on the Eu II 3819 Å and 4129 Å lines. The first of these is difficult to analyse, as this line is located in the blue wing of a strong Fe I 3819.62 Å line and, in addition, its blue wing is blended with the Cr I 3819.56 Å line. In stars with weak Eu lines, the above uncertainty, along with the uncertainty in the continuum placement, results in a systematically lower abundance from the Eu II 3819 Å compared to that from the Eu 4129 Å line. Instead of this line, we included in our analysis the Eu II 4205 Å line, which is free of blending in the metallicity range we deal with in this study. The average Eu abundances are calculated from Eu II 4129 Å and 4205 Å lines. Three of the sample stars have very weak lines of Eu II with EW < 6 mÅ, and we used their Eu abundances with caution. These measurements are marked with a < symbol in Table 5.

### 3.6. Element abundances and uncertainties

Table 5 presents the NLTE abundance ratios and fractions of odd Ba isotopes F<sub>odd</sub> = (N(<sup>135</sup>Ba)+N(<sup>137</sup>Ba))/N(Ba). We adopted the following solar abundances from Lodders (2021): log ε(Sr)⊙ = 2.88, log ε(Ba)⊙ = 2.17, and log ε(Eu)⊙ = 0.51.

Uncertainties in [X/H] ratios were computed, including uncertainties in stellar atmosphere parameters along with the dispersion of the single line measurements around the mean σ<sub>st</sub> = √(Σ(log ε - log εᵢ)<sup>2</sup>/(N - 1)), where N is the total number of lines. For our Sr, Ba, and Eu abundance determination, we assume that the dispersion and the systematic uncertainties caused by variations in stellar parameters are uncorrelated. This assumption is valid since, for each of the studied elements, their spectral lines belong to the same species, have nearly the same intensity, and originate from atomic levels with similar or identical E<sub>exc</sub>. Therefore, changes in T<sub>eff</sub>, log g, and ξᵢ do not contribute to line-to-line scatter, allowing the total abundance ratio [X/H] uncertainty, to be calculated as follows:

$$ \sigma^2_{[X/H]} = \sigma^2(X) + \sigma^2(T_{\text{eff}}) + \sigma^2(\log g) + \sigma^2(\xi_i). $$

When computing abundance ratios [X/Y], where X and Y are Sr, Ba, and Eu, we neglected the uncertainty in log g, since a shift in log g produces the same abundance shift for the above elements. The corresponding uncertainties caused by changes in T<sub>eff</sub> and ξᵢ are reduced compared to those for absolute abundances, since changes in these parameters affect the Sr II, Ba II, and Eu II spectral lines in the same way. For abundance ratios X and Y, the uncertainty in their ratio is calculated as follows:

$$ \sigma^2_{[X/Y]} = \sigma^2(X) + \sigma^2(Y) + \sigma^2(\Delta T_{\text{eff}}) + \sigma^2(\Delta \xi_i). $$

## 4. Discussion

### 4.1. Abundance trends

Figure 6 shows the [Sr/Ba] ratio as a function of [Ba/Eu] in our sample stars. The [Ba/Eu] ratio ranges from -0.87 in an r-II type star to -0.36 in a star with low Eu abundance. Three r-II sample stars exhibit consistently low [Ba/Eu] values, with an average ratio of [Ba/Eu]<sub>r</sub> = -0.81 ± 0.05, which aligns with the pure r-process ratio predictions [Ba/Eu]<sub>r</sub> = -0.8 (Bisterzo et al. 2014), -0.9 (Prantzos et al. 2020), -0.7 (Arlandini et al. 1999), and an

<footnote n="3">³ https://spectrum.inasan.ru/nLTE2/</footnote>
<figure><img src="image_8.png" /></figure>

**Fig. 6.** Abundance ratios in r-II (circles), r-I (triangles) stars, and sample stars with normal [Eu/Fe] < −0.3 (squares). Smaller symbols represent stars with [Fe/H] > −2.5. An x symbol indicates a pure r-process ratio derived from average abundances in r-II sample stars. Different lines show abundance ratios calculated for mixtures of r-process and s-process with varying [Sr/Ba]$_{earlyS}$. See the legend for designations.

empirical ratio [Ba/Eu]$_{r}$ = −0.87 ± 0.06 based on NLTE analysis of r-II stars (Mashonkina & Christlieb 2014). The [Sr/Ba] ratio spans a wide range, from [Sr/Ba]$_{r}$ = −0.31 ± 0.02 in r-II stars to a maximum of 0.74.

Figure 7 presents the abundance difference between the subordinate and resonance lines, ΔBa (sub.-res.), computed with different r- to s-process fractions (100%/0%, 80%/20%, 50%/50%, and 12%/88%) as a function of abundance ratios. We observe an increasing trend: the fraction of Ba even isotopes rises with [Ba/Eu], suggesting an increasing s-process contribution. It is unsurprising that the s-process produces Ba even isotopes and results in an increasing [Ba/Eu]. The key point is that an increasing trend is also observed when plotting ΔBa (sub.-res.) against [Sr/Eu] and [Sr/Ba] (Fig. 7, middle and bottom panels). This supports the idea of a common synthesis of additional Sr (beyond that produced via the main r-process) and Ba even isotopes, implying that the extra Sr originates from the s-process. This, in turn, clearly indicates that Sr and Ba are co-produced in this source. We do not claim that this was the only source, however, its contribution to Sr synthesis in the early Galaxy was dominant. Hereafter, we refer to this process as the early s-process.

## 4.2. Literature data on Ba isotope measurements

### 4.2.1. Data on individual sample stars

Applying the abundance comparison method in 1D NLTE for HD 122563, Mashonkina et al. (2008) found $F_{\text{odd}} = 0.22 \pm 0.15$, which agrees with our measurement of $F_{\text{odd}} = 0.14_{-0.04}^{+0.09}$. A slightly higher value of $F_{\text{odd}} = 0.39 \pm 0.06$ was obtained by Mashonkina & Belyaev (2019) when using modified stellar atmosphere parameters of HD 122563 and an updated Ba II model atom. Using the profile fitting method in 1D LTE, Gallagher et al. (2012) found a non-physical result, yielding a negative $F_{\text{odd}} = -0.12 \pm 0.07$ in HD 122563.

For HD 4306, Cescutti et al. (2021) found a good agreement between the observed spectra and 1D LTE synthetic spectra of the Ba II resonance lines when assuming a pure s-process Ba isotope ratio.

For the r-process enhanced r-I star HD 108317, Mashonkina & Belyaev (2019) found consistent within 0.05 dex NLTE abundances from the resonance and subordinate lines of Ba II when using $F_{\text{odd}} = 0.46$. Our measurement of $F_{\text{odd}} = 0.46_{-0.15}^{+0.24}$ is in agreement with Mashonkina & Belyaev (2019).

<figure><img src="image_9.png" /></figure>

**Fig. 7.** Abundance difference between the subordinate and the resonance lines for different r- and s-process fractions (see the legend for designations). The arrows indicate the stars with weak lines of Eu II with EW < 6 mÅ.

For the strongly r-process enhanced r-II star HE2327−5642, Mashonkina et al. (2010) found a pure r-process Ba isotope ratio in line with our findings.

### 4.2.2. Data on MP stars derived from the abundance comparison method

Mashonkina & Zhao (2006, hereafter MZ06) determined NLTE abundances of Eu and Ba, using the r-process and solar Ba isotope mixture for a sample of MW stars. We combined these measurements with Sr NLTE abundance determinations for the same stars presented by Mashonkina et al. (2007). In total, we found 24 MP thick disc and halo dwarfs with Sr, Ba, and Eu abundances along with Ba isotope ratio determinations. The original publications of MZ06 and Mashonkina et al. (2007) do not include abundances from individual Ba II lines. However, we retrieved them from surviving handwritten notes and digitised the data, as presented in Table A.1 in the Appendix.

The majority of the selected MZ06 stars have metallicities in the range −1.4 < [Fe/H] < −0.3, with the two stars at [Fe/H] = −2.1. The MZ06 sample traces more advanced galactic chemical evolution stages compared to those in the present study and is affected by the main s-process in low- to intermediate-mass AGB stars producing Sr and Ba. However, regardless of the source of Sr and Ba in the MZ06 subsample, the data align well with our findings in the ΔBa − [Ba/Eu] and ΔBa−[Sr/Eu] diagrams (Fig. 8).
<figure><img src="image_10.png" /><figcaption>Fig. 8. Abundance difference between the subordinate and the resonance lines for pure r-process (circles) and solar Ba isotope mixture (squares) in our sample stars (filled symbols) and data from MZ06 (open symbols). For comparison with MZ06, we plot our $\Delta$Ba (sub.-res.) computed with the r-process Ba isotope mixture from Arlandini et al. (1999).</figcaption></figure>

Jablonka et al. (2015) reports a Ba isotope ratio close to the s-process in two VMP stars in Sculptor dSph. These stars exhibit [Sr/Ba] ratios of 0.35 (ET0381) and -0.02 (Scl03059), both higher than the r-process value [Sr/Ba]$_r = -0.31$. This finding supports our hypothesis regarding the synthesis of additional Sr via the early s-process.

### 4.3. [Sr/Ba] ratio in the early s-process

We aim to determine the chemical properties of the early s-process by subtracting the contribution of the main r-process from the Sr, Ba, and Eu abundances of the sample stars. Here, we assume that the main r-process and the early s-process are the two major sources of n-capture elements in the early Galaxy. Potential additional sources of n-capture elements are discussed in the following section.

The key point is that Sr and Ba originate from both the main r-process and the early s-process, while Eu is exclusively produced by the main r-process. This fact enables using Eu as a reference element. To estimate the [Sr/Ba]$_{earlys}$ ratio produced by the early s-process, we use the following data: the [Sr/Eu] and [Ba/Eu] ratios in the sample stars; the pure r-process abundance ratios from the r-II sample stars: $[Sr/Eu]_r = -1.1$, $[Ba/Eu]_r = -0.8$, and $[Sr/Ba]_r = -0.3$; solar abundances, $\log \varepsilon(Sr)_\odot = 2.88$, $\log \varepsilon(Ba)_\odot = 2.17$, $\log \varepsilon(Eu)_\odot = 0.51$ from Lodders (2021).

In a given star, we can compute the number of Sr particles relative to the number of Eu particles. For example, in the Sun, this value amounts to $N(Sr)/N(Eu) = 224$. In r-II stars, the corresponding ratio is $N(Sr)/N(Eu) = 18$. Now let us compute the ratio in a given Sr-enhanced VMP star with $[Sr/Eu] = 0.1$, $[Sr/Ba] = 0.7$, and $[Ba/Eu] = -0.6$. In this star, $[Sr/Eu] = 0.1$ corresponds to $\log N(Sr)/N(Eu) = [Sr/Eu] + \log N(Sr)/N(Eu)_\odot = 2.45$, or $N(Sr)/N(Eu) = 282$. In other words, 18 Sr particles originate from the main r-process, while the remaining 264 Sr particles originate from the early s-process. Applying the same reasoning to the number of Ba particles, we find $N(Ba)/N(Eu) = 7$ for the pure r-process material and $N(Ba)/N(Eu) = 11$ in the example star. In this star, 7 Ba particles are produced by the main

r-process, while the remaining 4 Ba particles originate from the early s-process. Combining our results for Sr and Ba in the star, we find that, in the early s-process, $N(Sr)/N(Ba) = 264/4$, or $[Sr/Ba] = \log(264/4) - \log(N(Sr)/N(Ba))_\odot = 1.1$.

Thus, we find that the example star with $[Sr/Ba] = 0.7$ can be represented as a mixture of main r-process material with $[Sr/Ba]_r = -0.3$ and early s-process material with $[Sr/Ba]_{earlys} = 1.1$. From these values, we determine that in the example star, the main r-process and the early s-process are mixed in proportions of 75% and 25%, respectively.

Now, let us formalise this reasoning and define the procedure for estimating $[Sr/Ba]_{earlys}$. In a given star, its Sr, Ba, and Eu abundances can be represented as a mixture of r-process and early s-process material combined in a proportion $k_{mix} = N_r/N_{earlys}$, where $N_r$ and $N_{earlys}$ denote the total number of particles produced by the r-process and the early s-process, respectively. Let us define $N(Sr)_{earlys}/N(Ba)_{earlys}$ as $X_{earlys}$ and express the number of Sr, Ba, and Eu particles in a given star relative to $N(Ba)_{earlys}$ with

$$N(Sr)_* = X_{earlys} + k_{mix} \cdot N(Sr)_r/N(Ba)_r$$

$$N(Ba)_* = 1 + k_{mix}$$

$$N(Eu)_* = k_{mix} \cdot N(Eu)_r/N(Ba)_r$$

Using the measured abundance ratios in a given star, we deduce:

$$k_{mix} = \left( \frac{N(Ba)_*/N(Eu)_*}{N(Ba)_r/N(Eu)_r} - 1 \right)^{-1}$$

$$X_{earlys} = k_{mix} \cdot \left( \frac{N(Sr)_*/N(Eu)_*}{N(Ba)_r/N(Eu)_r} - N(Sr)_r/N(Ba)_r \right)$$

and finally, $[Sr/Ba]_{earlys} = \log X_{earlys} - (\log \varepsilon(Sr)_\odot - \log \varepsilon(Ba)_\odot)$.

To evaluate the effectiveness of $[Sr/Ba]_{earlys}$ synthesis in the early s-process, we selected sample stars with a significant contribution from this process – those with $[Sr/Ba] > 0.1$. From seven stars with $-0.79 < [Ba/Eu] < -0.36$ and $0.13 < [Sr/Ba] < 0.74$, we found an average $[Sr/Ba]_{earlys} = 1.02 \pm 0.24$. When excluding two stars with weak Eu lines, we obtained $[Sr/Ba]_{earlys} = 1.08 \pm 0.17$.

The selected seven stars exhibit different n-capture element abundance ratios; however, they consistently yield high $[Sr/Ba]_{earlys}$ values. For example, CES1732+2344, with $[Sr/Ba] = 0.13$ and $[Ba/Eu] = -0.79$, and CES1402+0941, with $[Sr/Ba] = 0.74$ and $[Ba/Eu] = -0.55$, provide close $[Sr/Ba]_{earlys}$ values of 1.3 and 1.1, respectively. We attribute the variations in n-capture element abundance ratios among these stars to differing relative contributions of the r-process and the early s-process to their chemical composition. For the two stars mentioned above, we estimate that r-process and early s-process material are mixed in the following proportions: 95%/5% in CES1732+2344 and 55%/45% in CES1402+0941.

The remaining sample stars are either r-I or r-II type stars, dominated by the r-process. An exception is the most metal-rich sample star, CES1436-2906, with $[Fe/H] = -2.15$. Its chemical composition can be explained by a mixture of r-process and the main s-process, which produces $[Sr/Ba] = 0$ (Fig. 6).

For stars with moderate $[Sr/Ba]$ enhancement, we found that the additional source of Sr synthesis is characterised by $[Sr/Ba]_{earlys} = 0.5$ (Fig. 6). This can be explained either by the main s-process contribution in addition to the early s-process or by variations in the effectiveness of Sr synthesis in the early s-process. The latter explanation seems more plausible, as the above stars have metallicities in the range $-3 < [Fe/H] < -2.8$, where the first intermediate-mass stars did not evolve yet to an AGB stage and the contribution of the main s-process is unlikely. Variations in the effectiveness of Sr synthesis in the
<figure><img src="image_11.png" /><figcaption>**Fig. 9.** R-process contribution to the 13 sample stars (excluding the three r-II stars) derived from the abundance ratios method (filled symbols) and Ba isotopes method (open symbols) as a function of [Sr/Eu]. Designations are the same as in Fig. 6.</figcaption></figure>

early s-process are supported by different [Sr/Ba] enhancement observed in different systems: stars in ultra-faint dwarf galaxies exhibit systematically lower [Sr/Ba] ratios compared to those found in classical dSphs (Mashonkina et al. 2017b; Ji et al. 2019; Reichert et al. 2020; Sitnova et al. 2021).

## 4.4. Potential sources of Sr and Ba synthesis beyond the main r-process and early s-process

When estimating $[\text{Sr/Ba}]_{\text{earlys}}$, we assume that the main r-process and the early s-process are the two primary n-capture element sources in the early Galaxy. However, we must consider the potential contribution of the i-process from massive stars (Banerjee et al. 2018) and weak r-process (Izutani et al. 2009) to the chemical composition of the sample stars. We do not consider a potential contribution from the vp process (Ghosh et al. 2022), as its nucleosynthesis yields for Sr and Ba are not available.

Similar to the r-process, the i-process follows a nucleosynthesis path through neutron-rich nuclei, followed by beta decay. Consequently, similar to the r-process, the i-process encounters the stable r-only isotopes $^{134}$Xe and $^{136}$Xe, which block the β-decay path and prevent the formation of $^{134}$Ba and $^{136}$Ba. As a result, the i-process is expected to effectively produce odd Ba isotopes $^{135}$Ba and $^{137}$Ba. For low-mass VMP AGB stars where the i-process occurs in their interiors, Choplin et al. (2021) found that $^{137}$Ba is surprisingly more abundant than $^{138}$Ba, which dominates in the solar system, as well as in the main s-process and r-process (Arlandini et al. 1999; Bisterzo et al. 2014). In other words, in terms of Ba isotopes, the contribution from the i-process appears similar to that of the r-process, increasing the fraction of odd Ba isotopes.

The r-process contribution to the chemical composition of a given star can be determined through element abundance ratios, as described above, along with another method based on Ba isotope analysis and achieving consistent abundances from the resonance and subordinate lines. For 13 sample stars, excluding the three r-II type stars, we computed the corresponding r-process fractions that result in $\Delta$Ba (sub.-res.) = 0. Figure 9 presents a comparison of r-process contributions derived from the two methods. Both methods yield results consistent within the error bars, though the uncertainties are sometimes uncomfortably large. If there were a significant contribution from the i-process, it would appear as a higher r-process contribution derived from the Ba isotope method compared to that from the element abundance ratio method, which is not observed in our data. The consistency of the two methods supports our hypothesis that, within our detection limits, the r-process and the early s-process are the two major Ba sources in the early Galaxy.

Regarding Sr synthesis, we demonstrated that, in addition to the main r-process, it is synthesised in the early s-process alongside the Ba even isotopes. However, the weak r-process may also produce a certain amount of Sr (Izutani et al. 2009). How can we identify a possible contribution to Sr abundance from the weak r-process in our sample stars? Let us take another look at Fig. 6. When excluding stars with [Fe/H] > −2.5, which may be contaminated by contribution from the main s-process, and r-II type stars, one can note that the remaining stars can be split into two groups according to their predicted $[\text{Sr/Ba}]_{\text{earlys}}$: five stars (CES0419–3651, CES0338–2402, CES0215–2554, CES1222+1136, and CES1427–2214) predict $[\text{Sr/Ba}]_{\text{earlys}} \approx 0.5$, while other five stars (CES1322–1355, CES1402+0941, CES1732+2344, CES0045–0932, CES2232–4138) predict $[\text{Sr/Ba}]_{\text{earlys}} \approx 1.0$ (see Fig. 6, the dashed and solid line, respectively.) Let us assume that the higher predicted $[\text{Sr/Ba}]_{\text{earlys}}$ can be interpreted as significant contamination by Sr produced in the weak r-process in addition to that Sr produced in the early s-process. Therefore, stars from the first group could be considered as having lower (or negligible) contribution from the weak r-process compared to stars from the second group.

Now let us examine the weak r-process yields available in Izutani et al. (2009). To trace the contribution of the weak r-process, one would need an element produced exclusively by that process. Sr, Y, and Zr are not suitable for distinguishing between contributions from the weak r-process and the early s-process. Figure 12 in Izutani et al. (2009) demonstrates that the weak r-process produces Zn as efficiently as Sr, Y, and Zr. In contrast, the s-process in VMP massive rotating stars produces a negligible amount of Zn (see Fig. 29 in Limongi & Chieffi 2018). Therefore, we select Zn as a tracer of the weak r-process. Following our hypothesis regarding the different contributions of the weak r-process to stars in the first and second groups, we might expect them to exhibit distinct Zn abundances, with higher values in the second group stars. Using the Zn abundances determined in LB22, we computed the average $[\text{Zn/Fe}] = 0.27 \pm 0.08$ for the first group and $0.30 \pm 0.12$ for the second group. Therefore, we found no significant difference in $[\text{Zn/Fe}]$ ratios between the two groups. In both groups, the stars span the same metallicity range, with $[\text{Fe/H}]$ from −2.9 to −2.5, and we found no systematic difference in the other chemical element abundances between the two groups. Based on our results, we found no strong evidence supporting a contribution from the weak r-process to the chemical composition of our sample stars.

## 4.5. Comparison with chemical evolution model predictions

We employed a stochastic chemical evolution model of the Galactic halo, making use of the GEMS code described in Rizzuti et al. (2025) and based on the stochastic models of Cescutti (2008) and Rizzuti et al. (2021). The model traces the early evolutionary stages within the first 1 Gyr of the halo formation, with a timestep of 1 Myr. Stochasticity is introduced by dividing the halo into 10 000 cubic regions, each considered isolated and containing the typical mass of gas swept by a SN-II explosion. Within each region, star formation, gas inflow and outflow are modelled, and gas recycling occurs at each step due to enrichment from evolving stars of different masses.

In the model, we included three sources of n-capture element synthesis: the main r-process, the s-process in rotating massive stars, and the s-process in AGB stars; however, the latter only contributes to the chemical evolution at $[\text{Fe/H}] > -2$. We assumed that the main r-process occurs in merging neutron stars with progenitor masses of 9–50 $M_{\odot}$. A binary fraction of
<figure><img src="image_12.png" /></figure>

**Fig. 10.** Comparison of the derived abundance ratios with the simulated data for long-living stars in the chemical evolution model. Designations are the same as in Fig. 6. Colour bar indicates the predicted fraction of odd Ba isotopes.

0.018 and a fixed coalescence time-scale of 1 Myr are adopted according to Matteucci et al. (2014) and Cescutti et al. (2015). We employed the empirical r-process ratios $[\text{Sr/Ba}]_r = -0.3$ and $[\text{Ba/Eu}]_r = -0.8$ for consistency with this study. The outcome of our chemical evolution results would not change significantly if magneto-rotational SNe are adopted as producers of r-process, instead of merging neutron stars with a very short delay. For the s-process in rotating massive stars (13–120 $M_{\odot}$), we used yields from the Limongi & Chieffi (2018) grid, ejected through stellar winds and explosions. For the main s-process in AGB stars, we used the yields from the FRUITY database (Cristallo et al. 2011). Iron yields were taken from Limongi & Chieffi (2018). For more details on the code, see Rizzuti et al. (2025).

Figure 10 compares the $[\text{Sr/Ba}]$ and $[\text{Ba/Eu}]$ ratios derived for our sample stars with the simulated data for long-living stars in the chemical evolution model. For consistency, we extracted model data for stars with metallicities in the range $-3.1 < [\text{Fe/H}] < -2.1$, matching our sample. The derived abundance ratios generally agree with model predictions. However, reproducing the highest observed $[\text{Sr/Ba}]$ values requires a slightly higher predicted $[\text{Sr/Ba}]$, which could be achieved by assuming slower rotation velocities in massive stars with $-3.1 < [\text{Fe/H}] < -2.1$.

Figure 11 compares the derived fraction of odd Ba isotopes as a function of $[\text{Sr/Ba}]$ with chemical evolution model predictions. The observed values span a wide range, with $F_{\text{odd}}$ varying from pure s-process to pure r-process values. Notably, the model predicts a significant number of stars with $[\text{Sr/Ba}] < -0.3$, originating from material produced by massive VMP stars with the highest rotational velocities, as faster rotation leads to lower $[\text{Sr/Ba}]$ production. Although our sample does not include stars enriched by such objects, we note that VMP stars with $[\text{Sr/Ba}] < -0.3$ should exhibit a low $F_{\text{odd}} < 0.2$. Determining Ba isotopic fractions in such stars would provide key constraints on the rotational velocities of massive VMP stars.

### 4.6. Other light n-capture elements: Y and Zr

We employ Sr as a tracer of light n-capture elements sources, although Y and Zr are also recognised as light n-capture elements and exhibit similar behaviour in VMP stars. Our conclusion regarding the joint synthesis of Sr and the even Ba isotopes can also be extended to Y and Zr. Figure 12 shows $F_{\text{odd}}$ as a function of $[\text{Sr/Ba}]$, $[\text{Y/Ba}]$, and $[\text{Zr/Ba}]$, with Y and Zr abundances

<figure><img src="image_13.png" /></figure>

**Fig. 11.** Comparison of the derived fraction of odd Ba isotopes with the simulated data for long-living stars in the chemical evolution model. Designations are the same as in Fig. 6. Colour bar indicates the predicted number of stars in each area of the diagram.

<figure><img src="image_14.png" /></figure>

**Fig. 12.** $F_{\text{odd}}$ as a function of $[\text{X/Ba}]$ (see legend) in the sample stars.

taken from LB22. All three elements display a similar trend: $F_{\text{odd}}$ decreases with increasing $[\text{X/Ba}]$, supporting the idea of co-production of Sr, Y, and Zr alongside the even Ba isotopes.

## 5. Conclusions

We presented a spectroscopic analysis of 16 VMP stars and determined their Sr, Ba, and Eu NLTE abundances, along with the contributions of the r- and s-processes to their Ba isotope mixture. For Ba isotope ratio analysis, we employed a method based on abundance comparisons between the resonance and subordinate lines of Ba II. This method is most effective within a specific stellar atmosphere and barium abundance range, where the EWs of the Ba II resonance lines range from 60 to 130 mÅ, making them sensitive to the Ba isotope ratio while avoiding excessive saturation.

Our results provide observational constraints on the source of Sr operating at the early epoch of Galactic chemical evolution. Our findings are summarised as follows:

* We find a higher s-process contribution to Ba isotopes in stars with greater $[\text{Sr/Eu}]$ ratio, suggesting that the additional Sr synthesis was due to the early s-process occurring in massive stars. Our finding matches the Galactic chemical evolution model predictions of Cescutti et al. (2013), Cescutti & Chiappini (2014), and Rizzuti et al. (2021).

* The early s-process produces not only Sr but also Ba; therefore it is not a weak s-process. The most likely candidate for the early s-process is a non-standard s-process occurring in massive stars (Frischknecht et al. 2012, 2016; Choplin et al. 2018; Limongi & Chieffi 2018).

* Using Sr-enhanced stars, we estimate that the $[\text{Sr/Ba}]$ ratio produced by the early s-process could be as high as
[Sr/Ba]$_{earlyS}$ = 1.1 ± 0.2. Regarding the potential synthesis of Sr and Ba in the i-process in massive stars, our results for Ba isotopes and element abundances argue that there was no detectable contribution from this process within the error bars in our sample stars.

## Data availability

The full Tables 2, 3, and 4 are available at the CDS https://cdsarc.cds.unistra.fr/viz-bin/cat/J/A+A/699/A262

**Acknowledgements.** The authors thank the referee for carefully reading the manuscript and providing valuable feedback. We acknowledge A. Alencastro Puls for reading the manuscript and providing feedback. TS is grateful to Yu. V. Pakhomov for providing his code for echelle orders merging. LL and CJH acknowledge the support by the State of Hesse within the Research Cluster ELEMENTS (Project ID 500/10.006). FR is a fellow of the Alexander von Humboldt Foundation. FR acknowledges support by the Klaus Tschira Foundation. FR and GC acknowledge the National Recovery and Resilience Plan (NRRP), Mission 4, Component 2, Investment 1.1, Call for tender No. 104 published on 2.2.2022 by the Italian Ministry of University and Research (MUR), funded by the European Union - NextGenerationEU - Project 'Cosmic POT' (PI: L. Magrini) Grant No. 2022X4TM3H by the MUR. GC thanks for the support INAF (Large Grant 2023, EPOCH and the MiniGrant 2022 Checs). GC and CJH acknowledge the European Union (ChETEC-INFRA, project no. 101008324). CJH acknowledges HFHF. PB acknowledges support from the ERC advanced grant N. 835087 - SPIAKID. This study has made use of the Keck Observatory Archive, which is operated by the W. M. Keck Observatory and the NASA Exoplanet Science Institute, under contract with the National Aeronautics and Space Administration.

## References

Aoki, W., Honda, S., Beers, T. C., et al. 2005, ApJ, 632, 611

Arcones, A., & Bliss, J. 2014, J. Phys. G Nuclear Phys., 41, 044005

Arlandini, C., Käppeler, F., Wisshak, K., et al. 1999, ApJ, 525, 886

Banerjee, P., Qian, Y.-Z., & Heger, A. 2018, ApJ, 865, 120

Barklem, P. S., Christlieb, N., Beers, T. C., et al. 2005, A&A, 439, 129

Bergemann, M. 2011, MNRAS, 413, 2184

Bisterzo, S., Travaglio, C., Gallino, R., Wiescher, M., & Käppeler, F. 2014, ApJ, 787, 10

Butler, K. 1984, Ph.D. Thesis, University of London, UK

Cescutti, G. 2008, A&A, 481, 691

Cescutti, G., & Chiappini, C. 2014, A&A, 565, A51

Cescutti, G., Chiappini, C., Hirschi, R., Meynet, G., & Frischknecht, U. 2013, A&A, 553, A51

Cescutti, G., Romano, D., Matteucci, F., Chiappini, C., & Hirschi, R. 2015, A&A, 577, A139

Cescutti, G., Morossi, C., Franchini, M., et al. 2021, A&A, 654, A164

Chiappini, C., Frischknecht, U., Meynet, G., et al. 2011, Nature, 472, 454

Choplin, A., Hirschi, R., Meynet, G., et al. 2018, A&A, 618, A133

Choplin, A., Siess, L., & Goriely, S. 2021, A&A, 648, A119

Choplin, A., Siess, L., Goriely, S., & Martinet, S. 2024, A&A, 684, A206

Christlieb, N., Beers, T. C., Barklem, P. S., et al. 2004, A&A, 428, 1027

Collet, R., Asplund, M., & Nissen, P. E. 2009, PASA, 26, 330

Cowan, J. J., & Rose, W. K. 1977, ApJ, 212, 149

Cowley, C. R., & Frey, M. 1989, ApJ, 346, 1030

Cristallo, S., Piersanti, L., Straniero, O., et al. 2011, ApJS, 197, 17

Eichler, M., Nakamura, K., Takiwaki, T., et al. 2018, J. Phys. G Nucl. Phys., 45, 014001

Fernandes de Melo, R., Lombardo, L., Alencastro Puls, A., et al. 2024, A&A, 691, A220

François, P., Depange, E., Hill, V., et al. 2007, A&A, 476, 935

Frischknecht, U., Hirschi, R., & Thielemann, F. K. 2012, A&A, 538, L2

Frischknecht, U., Hirschi, R., Pignatari, M., et al. 2016, MNRAS, 456, 1803

Fröhlich, C., Martínez-Pinedo, G., Liebendörfer, M., et al. 2006, Phys. Rev. Lett., 96, 142502

Gaia Collaboration (Brown, A. G. A., et al.) 2021, A&A, 649, A1

Gallagher, A. J., Ryan, S. G., García Pérez, A. E., & Aoki, W. 2010, A&A, 523, A24

Gallagher, A. J., Ryan, S. G., Hosford, A., et al. 2012, A&A, 538, A118

Gallagher, A. J., Ludwig, H. G., Ryan, S. G., & Aoki, W. 2015, A&A, 579, A94

Ghosh, S., Wolfe, N., & Fröhlich, C. 2022, ApJ, 929, 43

Giddings, J. 1981, Ph.D. Thesis, University of London, UK

Gratton, R. G., & Sneden, C. 1994, A&A, 287, 927

Gustafsson, B., Edvardsson, B., Eriksson, K., et al. 2008, A&A, 486, 951

Hansen, C. J., Primas, F., Hartman, H., et al. 2012, A&A, 545, A31

Hansen, C. J., Montes, F., & Arcones, A. 2014, ApJ, 797, 123

Hayek, W., Wiesendahl, U., Christlieb, N., et al. 2009, A&A, 504, 511

Izutani, N., Umeda, H., & Tominaga, N. 2009, ApJ, 692, 1517

Jablonka, P., North, P., Mashonkina, L., et al. 2015, A&A, 583, A67

Ji, A. P., Simon, J. D., Frebel, A., Venn, K. A., & Hansen, T. T. 2019, ApJ, 870, 83

Käppeler, F., Gallino, R., Bisterzo, S., & Aoki, W. 2011, Rev. Mod. Phys., 83, 157

Kobayashi, C., Karakas, A. I., & Lugaro, M. 2020, ApJ, 900, 179

Kochukhov, O. 2018, Astrophysics Source Code Library [record ascl:1805.015]

Kurucz, R. L. 2005, Mem. Soc. Astron. Ital. Suppl., 8, 14

Lambert, D. L., & Allende Prieto, C. 2002, MNRAS, 335, 325

Limongi, M., & Chieffi, A. 2018, ApJS, 237, 13

Lodders, K. 2021, Space Sci. Rev., 217, 44

Lombardo, L., Bonifacio, P., François, P., et al. 2022, A&A, 665, A10

Lombardo, L., Hansen, C. J., Rizzuti, F., et al. 2025, A&A, 693, A293

Magain, P. 1995, A&A, 297, 686

Magain, P., & Zhao, G. 1993, A&A, 268, L27

Mallinson, J. W. E., Lind, K., Amarsi, A. M., et al. 2022, A&A, 668, A103

Mashonkina, L. I., & Belyaev, A. K. 2019, Astron. Lett., 45, 341

Mashonkina, L., & Christlieb, N. 2014, A&A, 565, A123

Mashonkina, L., & Gehren, T. 2000, A&A, 364, 249

Mashonkina, L., & Gehren, T. 2001, A&A, 376, 232

Mashonkina, L., & Zhao, G. 2006, A&A, 456, 313

Mashonkina, L. I., Vinogradova, A. B., Ptitsyn, D. A., Khokhlova, V. S., & Chernetoova, T. A. 2007, Astron. Rep., 51, 903

Mashonkina, L., Zhao, G., Gehren, T., et al. 2008, A&A, 478, 529

Mashonkina, L., Christlieb, N., Barklem, P. S., et al. 2010, A&A, 516, A46

Mashonkina, L., Gehren, T., Shi, J.-R., Korn, A. J., & Grupp, F. 2011, A&A, 528, A87

Mashonkina, L., Jablonka, P., Pakhomov, Y., Sitnova, T., & North, P. 2017a, A&A, 604, A129

Mashonkina, L., Jablonka, P., Sitnova, T., Pakhomov, Y., & North, P. 2017b, A&A, 608, A89

Mashonkina, L., Pakhomov, Y. V., Sitnova, T., et al. 2022, MNRAS, 509, 3626

Mashonkina, L., Pakhomov, Y., Sitnova, T., et al. 2023, MNRAS, 524, 3526

Matteucci, F., Romano, D., Arcones, A., Korobkin, O., & Rosswog, S. 2014, MNRAS, 438, 2177

McWilliam, A., Preston, G. W., Sneden, C., & Searle, L. 1995, AJ, 109, 2757

Pakhomov, Y. V., Ryabchikova, T. A., & Piskunov, N. E. 2019, Astron. Rep., 63, 1010

Pignatari, M., Gallino, R., Meynet, G., et al. 2008, ApJ, 687, L95

Prantzos, N., Abia, C., Cristallo, S., Limongi, M., & Chieffi, A. 2020, MNRAS, 491, 1832

Qian, Y. Z., & Wasserburg, G. J. 2007, Phys. Rep., 442, 237

Raiteri, C. M., Busso, M., Gallino, R., & Picchio, G. 1991a, ApJ, 371, 665

Raiteri, C. M., Busso, M., Gallino, R., Picchio, G., & Pulone, L. 1991b, ApJ, 367, 228

Reichert, M., Hansen, C. J., Hanke, M., et al. 2020, A&A, 641, A127

Rizzuti, F., Cescutti, G., Matteucci, F., et al. 2021, MNRAS, 502, 2495

Rizzuti, F., Cescutti, G., Molaro, P., et al. 2025, A&A, 698, A118

Roederer, I. U. 2013, AJ, 145, 26

Roederer, I. U., Sakari, C. M., Placco, V. M., et al. 2018, ApJ, 865, 129

Rolfs, C., Trautvetter, H. P., & Rodney, W. S. 1987, Rep. Progr. Phys., 50, 233

Ryabchikova, T., Piskunov, N., Kurucz, R. L., et al. 2015, Phys. Scr, 90, 054005

Shah, S. P., Ezzeddine, R., Roederer, I. U., et al. 2024, MNRAS, 529, 1917

Simmerer, J., Sneden, C., Cowan, J. J., et al. 2004, ApJ, 617, 1091

Sitnova, T., Zhao, G., Mashonkina, L., et al. 2015, ApJ, 808, 148

Sitnova, T. M., Mashonkina, L. I., & Ryabchikova, T. A. 2016, MNRAS, 461, 1000

Sitnova, T. M., Yakovleva, S. A., Belyaev, A. K., & Mashonkina, L. I. 2020, Astron. Lett., 46, 120

Sitnova, T. M., Mashonkina, L. I., Tатарников, А. М., et al. 2021, MNRAS, 504, 1183

Sneden, C., Cowan, J. J., Lawler, J. E., et al. 2003, ApJ, 591, 936

Sneden, C., Cowan, J. J., & Gallino, R. 2008, ARA&A, 46, 241

Spite, M., & Spite, F. 1978, A&A, 67, 23

Travaglio, C., Gallino, R., Arnone, E., et al. 2004, ApJ, 601, 864

Tsymbal, V., Ryabchikova, T., & Sitnova, T. 2019, in Astronomical Society of the Pacific Conference Series, 518, Physics of magnetic stars, eds. I. I. Romanyuk, I. A. Yakunin, & D. O. Kudryavtsev, 247

Wood, M. P., Lawler, J. E., Sneden, C., & Cowan, J. J. 2013, ApJS, 208, 27

Woosley, S. E., & Hoffman, R. D. 1992, ApJ, 395, 202
## Appendix A: Additional table

<table><caption>**Table A.1.** Stellar parameters and abundances from individual lines of Sr II, Ba II, and Eu II from MZ06 and Mashonkina et al. (2007).</caption><thead><tr><th rowspan="2">Name</th><th rowspan="2">Pop.*</th><th rowspan="2">T<sub>eff</sub>, K</th><th rowspan="2">log g</th><th rowspan="2">[Fe/H]</th><th rowspan="2">ξ<sub>t</sub>, km s<sup>-1</sup></th><th rowspan="2">[Mg/Fe]</th><th colspan="2">Ba II 4554</th><th colspan="2">Ba II 5853</th><th colspan="2">log ε</th><th colspan="2">Sr II 4215</th><th colspan="2">Eu II 4129</th></tr><tr><th>r, NLTE</th><th>sol., NLTE</th><th>NLTE</th><th>LTE</th><th>NLTE</th><th>LTE</th><th>NLTE</th><th>LTE</th><th>NLTE</th><th>LTE</th><th>NLTE</th><th>LTE</th></tr></thead><tbody><tr><td>HD 65583</td><td>t</td><td>5320</td><td>4.55</td><td>-0.73</td><td>0.80</td><td>0.39</td><td>1.41</td><td>1.49</td><td>1.43</td><td>1.46</td><td>1.44</td><td>1.54</td><td>2.20</td><td>2.21</td><td>0.27</td><td>0.23</td></tr><tr><td>HD 62301</td><td>t</td><td>5940</td><td>4.06</td><td>-0.69</td><td>1.30</td><td>0.30</td><td>1.38</td><td>1.56</td><td>1.48</td><td>1.54</td><td>1.43</td><td>1.65</td><td>2.18</td><td>2.25</td><td>0.21</td><td>0.15</td></tr><tr><td>HD 30649</td><td>t</td><td>5820</td><td>4.28</td><td>-0.47</td><td>1.20</td><td>0.35</td><td>1.60</td><td>1.68</td><td>1.65</td><td>1.72</td><td>1.62</td><td>1.82</td><td>2.34</td><td>2.36</td><td>0.39</td><td>0.35</td></tr><tr><td>72 Her</td><td>t</td><td>5735</td><td>4.24</td><td>-0.34</td><td>1.00</td><td>0.38</td><td>1.71</td><td>1.83</td><td>1.74</td><td>1.81</td><td>1.74</td><td>1.94</td><td>2.52</td><td>2.54</td><td>0.54</td><td>0.51</td></tr><tr><td>HD 69611</td><td>t</td><td>5820</td><td>4.18</td><td>-0.60</td><td>1.20</td><td>0.43</td><td>1.39</td><td>1.55</td><td>1.52</td><td>1.57</td><td>1.48</td><td>1.69</td><td>2.35</td><td>2.37</td><td>0.30</td><td>0.26</td></tr><tr><td>HD 102158</td><td>t</td><td>5760</td><td>4.24</td><td>-0.46</td><td>1.10</td><td>0.40</td><td>1.63</td><td>1.70</td><td>1.63</td><td>1.70</td><td>1.61</td><td>1.83</td><td>2.44</td><td>2.46</td><td>0.42</td><td>0.38</td></tr><tr><td>HD 22879</td><td>t</td><td>5870</td><td>4.27</td><td>-0.86</td><td>1.20</td><td>0.44</td><td>1.24</td><td>1.39</td><td>1.34</td><td>1.37</td><td>1.31</td><td>1.48</td><td>2.13</td><td>2.17</td><td>0.10</td><td>0.04</td></tr><tr><td>HD 201891</td><td>t</td><td>5940</td><td>4.24</td><td>-1.05</td><td>1.20</td><td>0.41</td><td>1.05</td><td>1.20</td><td>1.10</td><td>1.10</td><td>1.11</td><td>1.21</td><td>1.85</td><td>1.89</td><td>-0.09</td><td>-0.16</td></tr><tr><td>HD 84937</td><td>h</td><td>6350</td><td>4.03</td><td>-2.07</td><td>1.70</td><td>0.36</td><td>0.12</td><td>0.19</td><td></td><td></td><td>0.14</td><td>-0.04</td><td>0.73</td><td>0.69</td><td>-0.96</td><td>-1.12</td></tr><tr><td>HD 68017</td><td>t</td><td>5630</td><td>4.45</td><td>-0.40</td><td>0.90</td><td>0.34</td><td>1.59</td><td>1.73</td><td>1.69</td><td>1.74</td><td>1.69</td><td>1.84</td><td>2.43</td><td>2.45</td><td>0.42</td><td>0.38</td></tr><tr><td>HD 45282</td><td>h</td><td>5280</td><td>3.12</td><td>-1.52</td><td>1.40</td><td>0.37</td><td>0.62</td><td>0.81</td><td>0.64</td><td>0.66</td><td>0.59</td><td>0.79</td><td>1.32</td><td>1.36</td><td>-0.38</td><td>-0.44</td></tr><tr><td>HD 18757</td><td>t</td><td>5710</td><td>4.34</td><td>-0.28</td><td>1.00</td><td>0.32</td><td>1.69</td><td>1.85</td><td>1.82</td><td>1.88</td><td>1.81</td><td>1.98</td><td>2.56</td><td>2.58</td><td>0.51</td><td>0.47</td></tr><tr><td>HD 3795</td><td>t</td><td>5370</td><td>3.82</td><td>-0.64</td><td>1.00</td><td>0.39</td><td>1.59</td><td>1.68</td><td>1.59</td><td>1.68</td><td>1.58</td><td>1.79</td><td>2.32</td><td>2.34</td><td>0.46</td><td>0.43</td></tr><tr><td>HD 103095</td><td>h</td><td>5110</td><td>4.66</td><td>-1.35</td><td>0.85</td><td>0.28</td><td>0.80</td><td>0.90</td><td>0.85</td><td>0.84</td><td>0.87</td><td>0.91</td><td>1.46</td><td>1.47</td><td>-0.26</td><td>-0.29</td></tr><tr><td>HD 194598</td><td>h</td><td>6060</td><td>4.27</td><td>-1.12</td><td>1.45</td><td>0.29</td><td>0.95</td><td>1.12</td><td>1.11</td><td>1.09</td><td>1.01</td><td>1.08</td><td>1.69</td><td>1.75</td><td>-0.00</td><td>-0.07</td></tr><tr><td>HD 37124</td><td>t</td><td>5610</td><td>4.44</td><td>-0.44</td><td>0.90</td><td>0.32</td><td>1.56</td><td>1.66</td><td>1.66</td><td>1.71</td><td>1.63</td><td>1.78</td><td>2.39</td><td>2.41</td><td>0.40</td><td>0.36</td></tr><tr><td>HD 10519</td><td>t</td><td>5710</td><td>4.00</td><td>-0.64</td><td>1.10</td><td>0.43</td><td>1.50</td><td>1.63</td><td>1.53</td><td>1.59</td><td>1.51</td><td>1.71</td><td>2.33</td><td>2.34</td><td>0.28</td><td>0.24</td></tr><tr><td>HD 222794</td><td>t</td><td>5620</td><td>3.94</td><td>-0.69</td><td>1.20</td><td>0.40</td><td>1.34</td><td>1.47</td><td>1.41</td><td>1.48</td><td>1.44</td><td>1.65</td><td>2.25</td><td>2.26</td><td>0.23</td><td>0.18</td></tr><tr><td>HD 64606</td><td>t</td><td>5320</td><td>4.54</td><td>-0.89</td><td>0.95</td><td>0.37</td><td>1.15</td><td>1.28</td><td>1.25</td><td>1.27</td><td>1.20</td><td>1.29</td><td>1.99</td><td>2.01</td><td>0.14</td><td>0.10</td></tr><tr><td>HD 132142</td><td>t</td><td>5240</td><td>4.58</td><td>-0.39</td><td>0.70</td><td>0.31</td><td>1.67</td><td>1.75</td><td>1.73</td><td>1.76</td><td>1.73</td><td>1.83</td><td>2.41</td><td>2.42</td><td>0.41</td><td>0.38</td></tr><tr><td>HD 112758</td><td>t</td><td>5240</td><td>4.62</td><td>-0.43</td><td>0.70</td><td>0.36</td><td>1.60</td><td>1.68</td><td>1.64</td><td>1.67</td><td>1.66</td><td>1.76</td><td>2.34</td><td>2.35</td><td>0.39</td><td>0.36</td></tr><tr><td>HD 144579</td><td>t</td><td>5330</td><td>4.59</td><td>-0.69</td><td>0.80</td><td>0.38</td><td>1.45</td><td>1.48</td><td>1.44</td><td>1.48</td><td></td><td></td><td>2.16</td><td>2.17</td><td>0.31</td><td>0.28</td></tr><tr><td>HD 148816</td><td>t</td><td>5880</td><td>4.07</td><td>-0.78</td><td>1.20</td><td>0.41</td><td>1.21</td><td>1.37</td><td>1.30</td><td>1.35</td><td>1.29</td><td>1.49</td><td>2.07</td><td>2.09</td><td>0.06</td><td>0.00</td></tr><tr><td>HD 193901</td><td>h</td><td>5780</td><td>4.46</td><td>-1.08</td><td>0.90</td><td>0.18</td><td>1.07</td><td>1.19</td><td>1.12</td><td>1.13</td><td>1.06</td><td>1.16</td><td>1.70</td><td>1.73</td><td>-0.10</td><td>-0.16</td></tr></tbody></table>

1 – 't' and 'h' stellar population types correspond to the thick disc and halo, respectively.

Here we present the digitalised hand-written line-by-line abundances used for average abundance calculations in MZ06 and Mashonkina et al. (2007). We only include MP stars where abundances of the three elements: Sr, Ba, and Eu, are available.