Free Access
Issue
A&A
Volume 657, January 2022
Article Number A50
Number of page(s) 10
Section Stellar structure and evolution
DOI https://doi.org/10.1051/0004-6361/202141775
Published online 05 January 2022

© ESO 2022

1. Introduction

Carbon-enhanced metal-poor (CEMP) stars form an important class of metal-poor giants, sub-giants, and dwarfs, with a large fraction of them showing enhanced abundances of heavy elements (see Beers & Christlieb 2005; Frebel 2018 for a general review). Among the different types of CEMP stars, the CEMP-s stars are enriched with products of s-process nucleosynthesis, the CEMP-r stars are enriched with the products of r-process nucleosynthesis, and CEMP-r/s stars are enriched with products of i-process nucleosynthesis. Understanding the diverse abundance patterns exhibited by different groups of CEMP stars that are believed to be associated with different formation mechanisms has been a challenge. In Goswami et al. (2021), we present a detailed analysis and discussion on the classification criteria of CEMP stars, as well as the formation scenarios of CEMP stars put forward by different authors (Cowan & Rose 1977; Hill et al. 2000; Qian & Wasserburg 2003; Cohen et al. 2003; Jonsell et al. 2006; Campbell & Lattanzio 2008; Campbell et al. 2010; Stancliffe et al. 2011; Herwig et al. 2011; Doherty et al. 2015; Abate et al. 2016; Jones et al. 2016; Banerjee et al. 2018; Clarkson et al. 2018; Denissenkov et al. 2017; Côté et al. 2018). In this paper, we report an extremely metal-poor carbon-enhanced star, HE 1005–1439, whose surface chemical composition is found to be enriched with both s-process and i-process nucleosynthesis that forms a new class of object with a distinct abundance pattern. The peculiar abundance pattern, observed for the first time in a CEMP star, was investigated based on a parametric-model-based analysis that revealed almost equal contributions from both the s-process and the i-process to its surface chemical composition. We examined various production mechanisms and formation scenarios for this object. A formation scenario involving effective proton ingestion episodes (PIEs) triggering i-process nucleosynthesis followed by s-process asymptotic giant branch (AGB) nucleosynthesis with limited third-dredge-up (TDU) episodes seems to be most promising for this type of object.

Literature surveys show that this object had been studied earlier by different groups (Aoki et al. 2007; Schuler et al. 2008; Yong et al. 2013; Caffau et al. 2017). However, these studies were limited by the number of elements for which abundances were estimated. As the abundances of the neutron-capture elements except Ba are not reported in the literature, we re-visited the object and estimated the atmospheric parameters as well as abundances of ten light elements and twelve heavy elements. The derived abundances were then carefully investigated with an aim to understand the origin and formation mechanism(s) of the object.

2. Source of spectra

The wavelength-calibrated high-resolution (R ∼ 50 000) spectra of HE 1005–1439 used in this study are retrieved from the SUBARU archive1. Spectra obtained at four different epochs using the high-dispersion spectrograph (HDS) (Noguchi et al. 2002) attached to the 8.2m Subaru Telescope cover the 3515–6780 Å wavelength range. The spectra acquired on October 26, 2002 (single exposure), and December 8, 2003 (five exposures) cover 4020–6780 Å with a gap of ∼70 Å from 5370–5440 Å, and the spectra acquired on October 28, 2002 (single exposure) and May 26, 2003 (three exposures) cover the 3515–5270 Å wavelength region with a gap of 15 Å from 4380–4395 Å. For our studies, spectra obtained on the same dates are combined to increase the signal-to-noise ratio (S/N). The sample spectra of the programme star at two different wavelength regions are shown in Fig. 1.

thumbnail Fig. 1.

Sample spectra of the programme star in the 4055–4065 Å (upper panel) and 5160–5190 Å (bottom panel) wavelength regions.

3. Radial velocities and stellar atmospheric parameters

We measured the radial velocity of HE 1005–1439 from spectra acquired at four different epochs using several clean unblended lines. The estimated radial velocities with 98.54 ± 0.75 km s−1 (October 26, 2002), 99.17 ± 0.90 km s−1 (October 28, 2002), 48.95 ± 0.46 km s−1 (May 26, 2003), and 103.53 ± 0.46 km s−1 (December 8, 2003) show that the object is a radial velocity variable, and that it is likely to be in a binary system.

We determined the photometric temperatures of HE 1005–1439 using broad-band colours with colour-temperature calibrations available for giants (Alonso et al. 1999) based on the infrared flux method (IRFM) following the procedure described in Goswami et al. (2006, 2016). The 2MASS photometric magnitudes for J, H, and K are taken from Cutri et al. (2003). The photometric temperatures corresponding to (BV), (VK), (JH), and (JK) at [Fe/H] ∼ −3.0 are found to be Teff(BV) = 4650 K, Teff(VK) = 5323 K, Teff(JH) = 4817 K, and Teff(JK) = 5017 K, respectively. As the photometric temperature Teff(JK) is independent of metallicity (Alonso et al. 1996, 1999), we used Teff(JK) as an initial guess for estimating the spectroscopic temperature of HE 1005–1439 through an iterative process of selecting the appropriate model atmosphere.

Following the detailed procedure described in Goswami et al. (2021), the stellar atmospheric parameters Teff, log g, micro-turbulent velocity (ζ), and metallicity [Fe/H] of the programme star are derived using 45 clean lines of Fe I and five lines of Fe II. The excitation potentials of these lines are in the 0.0–5.0 eV range. The list of lines with the measured equivalent widths and atomic line information is presented in Table A.1.

An updated version of MOOG software by Sneden (Sneden 1973) that assumes local thermodynamic equilibrium (LTE) conditions was used for our analysis. The model atmospheres are used from the Kurucz grid of model atmospheres with no convective overshooting2. The solar abundances are adopted from (Asplund et al. 2009). The effective temperature and the micro-turbulent velocity are fixed by the conventional methods of excitation potential balance and equivalent width balance. The ionisation equilibrium method of equating the abundance of Fe derived from neutral and ionised Fe lines determines the surface gravity. The metallicity is given by the abundance of Fe derived from Fe I and Fe II lines. The estimates of stellar parameters along with the literature values are listed in Table 1.

Table 1.

Derived atmospheric parameters of our programme star, and literature values.

Our estimate of effective temperature is 170 K higher than that of Aoki et al. (2007). Aoki et al. (2007) estimated the effective temperature from (VK), (VR), (VI), and (RI) photometric colour indices, using the empirical temperature calibration scale of Alonso et al. (1999). In general, due to the presence of C2 and CN molecular bands in carbon enhanced stars, (VK) gives lower temperature than the other colour indices. However, in the case of HE 1005–1439 temperature estimate from Teff–(VK) calibration relation is found to be much higher than those obtained from other colour indices. This was also noticed by Aoki et al. (2007) and due to this discrepancy, adopted a lower value for temperature (∼5000 K) closer to those derived using different colour indices. Yong et al. (2013) estimated the effective temperature of HE 1005–1439 using IRFM adopting the colour-temperature relations given in Ramírez & Meléndez (2005). Our estimate of effective temperature is close to that of Yong et al. (2013).

Our estimate of micro-turbulent velocity (ζ) is lower than that of Aoki et al. (2007) and Yong et al. (2013). For their analysis, Yong et al. (2013) used the same measured equivalent widths and lines reported by Aoki et al. (2007). While Aoki et al. (2007) used 19 lines of Fe I covering a range of equivalent widths 16.7 mÅ–88.0 Å, we have used 45 lines of Fe I covering a range of equivalent widths 15.6 Å–127.1 Å. With more clean lines covering a good range in line strength, we are confident about our estimated micro-turbulent velocity.

Our measured surface gravity log g is similar to that of Aoki et al. (2007). The log g derived by Aoki et al. (2007) is based on the ionisation equilibrium method similar to the one that we have used. Yong et al. (2013) used Y2 isochrones (Demarque et al. 2004) to determine log g, assuming an age of 10 Gyr, and [α/Fe] = +0.3 and reported a higher log g (∼2.55) for this object.

We estimated the mass of HE 1005–1439 from its position in the Hertzsprung-Russell (HR) diagram (log(L/L) versus log(Teff) plot). The value of the parallax (=0.2733 mas) is taken from Gaia Collaboration (2018) and the V (=13.52) magnitude is taken from SIMBAD. The bolometric correction is determined based on the empirical calibration equations of Alonso et al. (1999). Interstellar extinction for HE 1005–1439 is calculated using the formulae by Chen et al. (1998). The value of log(L/L) is found to be ∼1.83. We have used the updated BaSTI-IAC evolutionary tracks3 (Hidalgo et al. 2018) generated for [Fe/H] = −3.2 and [α/Fe] = 0.4, including overshooting and diffusion, to estimate the mass of the star. The mass of the object is found to be 0.8 M. The log g value is calculated using the following relation:

log ( g / g ) = log ( M / M ) + 4 log ( T eff / T eff ) + 0.4 ( M bol M bol ) . $$ \begin{aligned}&\log ({g}/{g}_{\odot }) = \log (M/M_{\odot }) + 4\log (T_{\rm eff}/T_{\rm eff\odot }) \nonumber \\&\qquad \qquad \quad + 0.4(M_{\rm bol} - M_{\rm bol\odot }) .\end{aligned} $$(1)

The adopted solar values are log g⊙ = 4.44, Teff⊙ = 5770 K and Mbol⊙ = 4.74 mag (Yang et al. 2016). Using this method, we found a log g (∼2.32) comparable to that of Yong et al. (2013). This value is much larger than our spectroscopic log g value ∼1.8. We note that the evolutionary tracks and isochrones highly depend on the opacity in the stellar atmosphere. As BaSTI-IAC evolutionary tracks are generated using normal carbon and without considering the influence of high carbon, log g values determined using such evolutionary tracks and isochrones may lead to erroneous estimates. This may also explain the discrepancy between the log g values obtained by us from spectroscopy and that reported by Yong et al. (2013). We used the spectroscopic log g value for our analysis.

4. Results and discussions

4.1. Abundance analysis

We determined the elemental abundances by measuring the equivalent widths of the absorption lines due to neutral and ionised species of several elements and/or by applying a spectrum synthesis technique using the radiative transfer code MOOG (Sneden 1973) that assumes LTE and model atmospheres from the Kurucz grid of model atmospheres with no convective overshooting4. Elemental abundances of C, Na, Mg, Ca, Sc, Ti, Cr, and Mn, iron-peak elements Co and Ni, and neutron-capture elements Sr, Y, Ba, La, Ce, Pr, Nd, Eu, Dy, Er, Hf, and Pb were estimated. We used the method of spectrum synthesis calculations for elements showing hyperfine splitting (e.g. Sc, Mn, Ba, La, and Eu). The lines used in the analysis with the measured equivalent widths and atomic line information are presented in Table A.1. Atomic line information, such as the excitation potential and log gf values, were taken from the Kurucz database of atomic line list. The abundance results along with the literature values are presented in Table 2.

Table 2.

Elemental abundances in HE 1005–1439.

4.1.1. Light elements

We used spectrum synthesis calculations for estimating the abundances of carbon from molecular bands of carbon. The abundance of carbon was derived using the CH band near 4310 Å and C2 bands near 5160 Å and 5635 Å (Fig. 2). The slight difference in A(C) derived from C2 and CH molecular bands might have appeared due to the difference in S/N in the regions of the CH and C2 bands. The S/N of the spectra near 5200 Å (i.e. near the C2 bands) is ∼130, while the S/N is ∼50 near 4320 Å (i.e. near the CH band). Carbon is found to be enhanced in HE 1005–1439 with [C/Fe] ∼ 2.37. The estimated carbon isotopic ratio 12C/13C, obtained from spectrum synthesis calculations of the C2 band near 4740 Å (Fig. 3) is ∼5.0. We have taken the line lists of C2 bands at 5165 Å, 5635 Å, and 4740 Å, and of the CH band at 4310 Å from the ‘linemake’5 atomic and molecular line database.

thumbnail Fig. 2.

Spectral synthesis fits of C2 bands around 5165 Å (top panel) and 5635 Å (bottom panel). The dotted lines indicate the synthesised spectra and the solid lines indicate the observed spectra. Two alternative synthetic spectra are shown corresponding to Δ[C/Fe] = +0.3 (long-dashed line) and Δ[C/Fe] = −0.3 (short-dashed line).

thumbnail Fig. 3.

Spectral synthesis fits (dotted curves) of the C2 features around 4740 Å obtained with the adopted C abundance and 12C/13C value (dotted curve). The observed spectrum is shown by a solid curve. Two alternative fits with 12C/13C ∼ 1 (short-dashed line) and 12 (long-dashed line) are shown to illustrate the sensitivity of the line strengths to the isotopic carbon abundance ratios.

While HE 1005–1439 is enhanced in Na ([Na/Fe] = 1.37), Mg ([Mg/Fe] = 0.53) is moderately enhanced. Gehren et al. (2004) found that the systematic NLTE corrections required on the abundances of Na and Mg in case of metal-poor halo stars are −0.4 and +0.1, respectively.

The other light elements, Ca, Sc, Ti, Co, and Ni, are moderately enhanced in HE 1005–1439. The abundances of Cr and Mn are found to be sub-solar. However, Bergemann & Gehren (2008) and Bergemann et al. (2010) have found that in case of metal-poor giants, abundances of Co and Mn are underestimated in LTE calculations. Bergemann & Gehren (2008) performed NLTE calculations for Mn on a sample of fourteen stars and found that the NLTE abundances of Mn in all the sample stars are higher than the LTE abundances, in fact, the NLTE correction may go up to 0.5–0.7 dex at low metallicities. NLTE calculations done by Bergemann et al. (2010) for Co I lines revealed that the corrections may vary from 0.1–0.8 dex depending on the effective temperature and metallicity. They found that deviation from LTE is larger in the case of giants than it is for dwarfs. Our estimates of the abundances of light elements are found to match closely with the literature values (calculated with LTE assumption) of the programme star within error bars. Hyperfine splitting contributions of the lines (Table A.1) used for spectrum synthesis calculations of Sc and Mn are taken from linemake.

4.1.2. Neutron-capture elements

While HE 1005–1439 is found to be enhanced ([X/Fe] > 1.0) in heavy elements Ba, La, Ce, Pr, Nd, Er, Hf, and Pb, the elements Sr, Y, Eu, and Dy are moderately enhanced ([X/Fe] in the range 0.26 to 0.72). We could not detect lines due to Tc. Hyperfine splitting contributions for the spectrum synthesis calculations of the lines Ba II 5853.67 Å & Ba II 6141.71 Å are taken from McWilliam (1998) and for La II 4921.78 Å line the hyperfine splitting contributions are taken from Jonsell et al. (2006). For the spectrum synthesis calculations of the lines Ba II 6496.90 Å, La II 4086.71 Å, La II 4123.22 Å, Eu II 4129.73 Å, and Eu II 4205.04 Å, the hyperfine splitting contributions are taken from linemake. The abundance of Pb is estimated using the spectrum synthesis calculation of Pb I 4057.81 Å.

4.2. Abundance uncertainties

We estimated the total uncertainties on the elemental abundances as discussed in Goswami et al. (2021). Two components, namely random error (σran = σ s N $ \dfrac{\sigma_{s}}{\sqrt N} $, where σs represents the standard deviation of the abundance of a particular species derived using N number of lines of that species) and systematic error (σsys), contribute to the total uncertainties. While the random error arises due to the uncertainties on the factors like oscillator strength, equivalent width measurement, and line blending, the systematic error arises due to the uncertainties in estimating the stellar atmospheric parameters. Finally, the uncertainties on [X/Fe] are derived as follows:

σ [ X / Fe ] 2 = σ log ϵ 2 + σ [ Fe / H ] 2 , $$ \begin{aligned} \sigma ^{2}_{\rm [X/Fe]} = \sigma ^{2}_{\log \epsilon } + \sigma ^{2}_{\rm [Fe/H]} ,\end{aligned} $$(2)

σ log ϵ 2 = σ ran 2 + σ sys 2 , $$ \begin{aligned} \sigma ^{2}_{\log \epsilon } = \sigma ^{2}_{\rm ran} + \sigma ^{2}_{\rm sys} ,\end{aligned} $$(3)

σ sys 2 = ( δ log ϵ δ T ) 2 σ T eff 2 + ( δ log ϵ δ log g ) 2 σ log g 2 + ( δ log ϵ δ ζ ) 2 σ ζ 2 + ( δ log ϵ δ [ Fe / H ] ) 2 σ [ Fe / H ] 2 , $$ \begin{aligned}&\sigma ^{2}_{\rm sys} = \left( \dfrac{\delta \log \epsilon }{\delta T}\right)^{2}\sigma ^{2}_{T_{\rm eff}} + \left( \dfrac{\delta \log \epsilon }{\delta \log {g}}\right)^{2}\sigma ^{2}_{\log {g}}\nonumber \\&\qquad + \left( \dfrac{\delta \log \epsilon }{\delta \zeta }\right)^{2}\sigma ^{2}_{\zeta } + \left( \dfrac{\delta \log \epsilon }{\delta \mathrm{[Fe/H]}}\right)^{2}\sigma ^{2}_{\rm [Fe/H]}, \end{aligned} $$(4)

where, σTeff = 100 K, σlog g = 0.2 dex, σζ = 0.2 km s−1, and σ[Fe/H] = 0.15 dex represent the typical uncertainties on the stellar atmospheric parameters Teff, log g, ζ, and [Fe/H], respectively. We evaluated the partial derivatives appearing in Eq. (4) for the programme star, varying the stellar parameters Teff, log g, ζ, and [Fe/H] by ±100 K, ±0.2 dex, ±0.2 km s−1, and ±0.2 dex, respectively. We note that the uncorrelated nature of the uncertainties arising from the different stellar parameters in Eq. (4) may lead to the overestimation of the calculated uncertainties on log ϵ and [X/Fe]. The estimated uncertainties σlog ϵ and σ[X/Fe] are listed in Cols. 5 and 8, respectively, of Table 2.

4.3. Interpretation of results

Our analysis shows the object HE 1005–1439 to be an extremely metal-poor star with [Fe/H] = −3.03 in accordance with Aoki et al. (2007). Carbon and neutron-capture elements are found to be enhanced in the programme star. Following the classification criteria of CEMP stars (Goswami et al. 2021), the object is found to belong to the CEMP-s sub-group (Fig. 4a). However, as shown in Fig. 4b, the ratio of heavy-s process (hs) elements (Ba, La, Ce, and Nd) to the light-s process (ls) elements (Sr and Y), [hs/ls] (∼0.92) is closer to the value at which the CEMP-r/s stars peak (1.06) (Goswami et al. 2021).

thumbnail Fig. 4.

Filled red circles and filled black squares respectively represent literature CEMP-r/s and CEMP-s stars compiled by Goswami et al. (2021), and the blue star represents the programme star. Panel a: grid formed by the dotted red lines bound by 0.0 < [La/Eu] < 0.6 and 0.0 < [Ba/Eu] < 1.0 indicates the region defined for CEMP-r/s stars by Goswami et al. (2021). The grid formed by the black dashed lines bound by 0.5 < [La/Eu] < 0.7 represents the region where [Eu/Fe] > 1.0 classifies the stars as CEMP-r/s and [Eu/Fe] < 1.0 classifies the stars as CEMP-s. Panel b: red dashed line at [hs/ls] = 1.06 and the black dashed line at [hs/ls] = 0.65 represent the peaks of [hs/ls] for CEMP-r/s and CEMP-s, respectively, as shown by Goswami et al. (2021) in Fig. 13a.

Using the s-process yields calculated using the FRUITY6 model (Straniero et al. 2006; Cristallo et al. 2009b, 2011, 2015) at the same metallicity (z = 0.00002) of the programme star, and considering different masses, we were not able to reproduce the observed abundance pattern of heavy elements. In Fig. 5a, we show a comparison of the observed elemental abundances with the AGB model yields (normalised to the La abundance of HE 1005–1439), calculated for M = 1.3 M and M = 2.0 M. As can be seen in the top panel of the residual plot of Fig. 5d, the s-process AGB models over-produce the light s-process elements Sr and Y, under-produce the elements Pr, Er, and Hf, and over-produce the third s-process peak element Pb.

thumbnail Fig. 5.

Examples of theoretical model fits with the observed abundances of the star. In panels a, b, and c, the points with error bars indicate the observed abundances. (a) s-process model fits for M = 1.3 M & 2.0 M. (b) i-process model fits for n = 1012 cm−3 and 1014 cm−3. (c) Best-fit from parametric-model. (d) Residual plot.

In Fig. 5b, we show a comparison of the i-process model yields of Hampel et al. (2016) (normalised to the La abundance of HE 1005–1439), calculated for n = 1012 cm−3 and n = 1014 cm−3, with the observed elemental abundance pattern of the programme star. The i-process models (Hampel et al. 2016) alone at neutron-densities n = 1012–1015 cm−3 cannot satisfactorily reproduce the observed abundance pattern of the programme star either. We can see from the middle panel of residual plot (Fig. 5d) that the i-process model with n = 1012 cm−3 fits the light s-process elements Sr and Y satisfactorily but under-produces Ce and Pr. Again, i-process model with n = 1014 cm−3 under produces the light s-process elements Sr and Y but over-produces Ba, Eu, and Er. The diverse abundance pattern observed in HE 1005–1439, which could not be explained either by s-process AGB nucleosynthesis or by the i-process alone, prompted us to explore alternate production mechanisms that might have influenced its surface chemical composition.

4.3.1. Parametric-model-based study

We performed a parametric-model-based study to delineate the contributions of s-, i- and r-processes to the observed abundances of heavy elements of the programme star. We used s-process model yields of the FRUITY model at different masses (M = 1.3 − 2.0 M), the Solar System r-process residual pattern (stellar model) given in Arlandini et al. (1999), and i-process model yields of Hampel et al. (2016) at different neutron-densities (n = 1012–1015 cm−3). We excluded the element Pb from the parametric-model based study as the i-process model yields of Pb are not reported by Hampel et al. (2016). We normalised the elemental abundances of the models to the La abundance of HE 1005–1439. The observed elemental abundances of HE 1005–1439 are then fitted with the parametric-model function log ϵj = Cs Nsj + CiNij + CrNrj, where Nsj, Nij, and Nrj indicate the normalised abundance from the s-process, i-process, and r-process, respectively. Here, Cs, Ci, and Cr indicate the component coefficients corresponding to contributions from the s-process, i-process, and r-process, respectively. In order to find the best fit, we calculated χ2 for all the possible combinations of models of s-process (M = 1.3 − 2.0 M) and i-process (n = 1012–1015 cm−3) along with solar r-process residues. The minimum χ2 is achieved for a combination of an s-process model with M = 2.0 M and an i-process model with n = 1014 cm−3 with no contribution from the r-process. The parametric-model function gives an excellent fit to the light s-process elements Sr and Y and a satisfactory fit to the heavier neutron-capture elements. The best parametric-model fit, where the contribution from both s- and i-process are similar (Cs = 0.56, Ci = 0.44), is shown in Fig. 5c and the bottom panel of the residual plot (Fig. 5d).

4.3.2. Origin of the programme star: Possible formation scenarios

The parametric-model-based study clearly established that the surface chemical composition of HE 1005−1439 is influenced by almost equal contributions of s-process AGB nucleosynthesis and i-process nucleosynthesis. We attempted to capture a formation scenario for this object involving a binary picture. Our assumption on binarity is profoundly based on the radial velocity variations observed on a few epochs. The low value of 12C/13C (∼5.0) measured for HE 1005–1439 also points towards the extrinsic nature of carbon, and hence the heavy elements in the star. In the intrinsic carbon stars (as they are in AGB phase), third-dredge-up (TDU) episodes bring 12C and s-process material to the surface, and the 12C/13C ratio can increase up to > 100 depending on the initial mass of the object (Karakas & Lattanzio 2014). However, in extrinsic carbon giants, such as HE 1005–1439, the 12C/13C ratio decreases due to the first dredge-up (FDU), which brings 13C produced in the internal CNO cycle to the stellar atmosphere (Smith et al. 1993). Attributing the s-process contribution to the AGB mass-transfer scenario in the binary system, we carefully examined if any of the proposed formation scenarios of i-process nucleosynthesis available in literature could explain the i-process contribution in the observed abundance pattern of HE 1005−1439.

We investigated if the very late thermal pulse (VLTP) scenario proposed to explain the peculiar abundance pattern of the Sakurai object (V4334 Sagittarii) by Herwig et al. (2011) could explain the abundance peculiarities of HE 1005−1439. However, a characteristic property of VLTP in pre-white dwarf of producing light s-process elements 2 dex more than that of heavy s-process elements is not observed in HE 1005−1439.

Low-metallicity or zero-metallicity massive (20–30 M) stars (Banerjee et al. 2018) and super-AGB stars (9–11 M) (Jones et al. 2016) can pollute the ISM with i-process yields. However, an AGB mass-transfer scenario in a binary system formed from the i-process-enriched ISM also cannot explain the peculiar abundance pattern of the programme star due to its inadequacy to describe the low observed Pb abundance in HE 1005–1439. As shown in Fig. 5a and the top panel of Fig. 5d, the observed Pb abundance of HE 1005–1439 is about 0.30–0.60 dex lower than the prediction of the FRUITY model (Cristallo et al. 2009b, 2016) at [Fe/H] ∼ −3.0. Pb is the main product of the s-process in AGB stars at low metallicities (Cristallo et al. 2009b). The low Pb abundance of HE 1005–1439 is difficult to explain by the scenario of pre-enrichment of ISM with i-process material or a scenario involving any external source of i-process material such as rapidly accreting white dwarfs (Denissenkov et al. 2019). This is because in these cases Pb would be more than (if the i-process neutron exposure (τ) of the progenitor is sufficient to produce enough Pb) or at a similar level to the s-process prediction, but it would not be less, as Pb is not destroyed in β-decay.

We propose that the surface chemical composition of HE 1005–1439 may be attributed to mass-transfer from a now extinct AGB companion with both s- and i-process nucleosynthesis occurring under suitable conditions during its evolution at different thermal pulses. Such a formation scenario may not be unlikely as many studies on the evolution of AGB stars have shown that neutron densities required for the s-process and the i-process can be achieved in the intershell region with the help of the partial mixing of protons in the radiative conditions (Herwig 2000; Denissenkov & Tout 2003; Herwig et al. 2003; Cristallo et al. 2009b, 2011; Piersanti et al. 2013; Karakas & Lattanzio 2014) and efficient PIEs in the convective conditions (Hollowell et al. 1990; Fujimoto et al. 2000; Iwamoto et al. 2004; Campbell & Lattanzio 2008; Lau et al. 2009; Cristallo et al. 2009a, 2016; Choplin et al. 2021), respectively.

In a certain study (Cristallo et al. 2009a) on the AGB evolution of a 1.5 M model with [Fe/H] = −2.45 without α enhancement, it was found that a strong PIE is followed by a deep TDU. Due to the PIE, H-burning occurring in high-temperature convective conditions creates a huge amount of 13C, which leads to an efficient neutron production with neutron densities of the order of 1015 cm−3. The study also noted that the convective He-shell splits into two sub-shells when the energy released by proton capture reactions slightly exceeds the energy production at the base of the convective shell. In the lower sub-shell, the 13C (α, n) 16O reaction produces i-process neutron density, and the nucleosynthesis path goes away from the valley of β stability producing a very high local [hs/ls] ratio. The upper sub-shell, where the CNO cycle is the main energy source, is later engulfed by the envelope. After that, a standard TP-AGB phase follows with repetitive TPs and TDUs, and s-process nucleosynthesis occurs. The number of expected stars experiencing PIEs is significantly reduced with the introduction of α-element enhancements (Cristallo et al. 2016).

Cristallo et al. (2009a) noted that the minimum mass of the models to experience TDU is significantly lowered by PIEs. Their model, after the PIE, encounters 25 additional TDU episodes, each of which is followed by radiative burning of the 13C pocket. Although the final abundance pattern is a combination of i- and s-processes, excessive s-process nucleosynthesis after the PIE would remove the trace of the i-process and thus reduce the [hs/ls] ratio.

In a recent study on the evolution of a 1 M object at low metallicity ([Fe/H] = −2.5), Choplin et al. (2021) noticed three convective instabilities occurring at the beginning of the TP-AGB phase. The main PIE occurs during the third instability, which produces high neutron densities of about 4.3 × 1014 cm−3, and rich i-process nucleosynthesis occurs. We note that although they are of similar metallicities, while in the model of Choplin et al. (2021) no further TPs were possible after the main PIE, in that of Cristallo et al. (2009a) 25 TPs after the PIE masked the contribution of the i-process and resulted in an s-process surface abundance pattern.

We propose that a model that undergoes PIEs during the beginning of the TP-AGB phase (producing an i-process abundance pattern) followed by limited TPs (producing an s-process abundance pattern) might explain the abundance peculiarity of HE 1005–1439. This scenario is also likely to explain the low-Pb abundance in the programme star. The abundance of Pb depends on the time-integrated neutron-exposure (τ). From Fig. 1 of Hampel et al. (2019), it is clear that as neutron density increases, production of Pb starts at higher τ. Therefore, a low Pb abundance may indicate that due to the operation of both i- and s-processes in succession, neither of the processes received sufficient exposure time to produce enough Pb expected at [Fe/H] ∼ −3.0.

5. Conclusions

A detailed, high-resolution spectroscopic study of the extremely metal-poor (EMP) star HE 1005–1439 revealed a peculiar abundance pattern different from those typically exhibited by CEMP stars. While the CEMP stars’ classification criteria place the object in the CEMP-s group, the value of [hs/ls] ∼ 0.9 is closer to the value (1.06) at which CEMP-r/s stars peak. The abundance pattern could not be reproduced either by the s-process or the i-process model predictions alone. However, a parametric-model based analysis clearly indicated that similar contributions from both the s- and i-process might have resulted in the observed abundance pattern of HE 1005–1439. We propose that the origin of the observed peculiar abundance pattern may be attributed to mass transfer from a now extinct AGB companion where both i- and s-process nucleosynthesis took place during various stages of the AGB evolution with PIEs triggering i-process followed by s-process AGB nucleosynthesis with a few TDU episodes.

Several uncertainties such as initial mass, metallicity, treatment of opacities, nuclear rates, and mixing mechanisms affect the theoretical understanding of AGB stars. The observational constraints derived from the programme star and the proposed scenario might provide important clues for a better understanding of the interplay between PIEs and partial mixing of protons in the intershell region and also the conditions resulting in a pure s- or i-process surface abundance pattern in low-mass, low-metallicity AGB stars. We believe that the proposed scenario will be helpful in explaining the overlap of [hs/ls] ratio in CEMP-s and CEMP-r/s stars (see Fig. 13 of Goswami et al. 2021) and the smooth transition of elemental abundances from the CEMP-s to CEMP-r/s regime (Goswami et al. 2021).


5

linemake contains laboratory atomic data (transition probabilities, hyperfine and isotopic substructures) published by the Wisconsin Atomic Physics and the Old Dominion Molecular Physics groups. These lists and accompanying line list assembly software have been developed by C. Sneden and are curated by V. Placco at https://github.com/vmplacco/linemake

Acknowledgments

Funding from DST SERB project EMR/2016/005283 is gratefully acknowledged. We thank Melanie Hampel for the i-process yields. We thank the referee Elisabetta Caffau for many constructive suggestions and useful comments. This work made use of the SIMBAD astronomical database, operated at CDS, Strasbourg, France, the NASA ADS, USA and data from the European Space Agency (ESA) mission Gaia (https://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www.cosmos.esa.int/web/gaia/dpac/consortium).

References

  1. Abate, C., Stancliffe, R. J., & Liu, Z.-W. 2016, A&A, 587, A50 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. Alonso, A., Arribas, S., & Martinez-Roger, C. 1996, A&A, 313, 873 [NASA ADS] [Google Scholar]
  3. Alonso, A., Arribas, S., & Martínez-Roger, C. 1999, A&AS, 140, 261 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Aoki, W., Beers, T. C., Christlieb, N., et al. 2007, ApJ, 655, 492 [NASA ADS] [CrossRef] [Google Scholar]
  5. Arlandini, C., Käppeler, F., Wisshak, K., et al. 1999, ApJ, 525, 886 [Google Scholar]
  6. Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 [NASA ADS] [CrossRef] [Google Scholar]
  7. Banerjee, P., Qian, Y.-Z., & Heger, A. 2018, ApJ, 865, 120 [Google Scholar]
  8. Beers, T. C., & Christlieb, N. 2005, ARA&A, 43, 531 [NASA ADS] [CrossRef] [Google Scholar]
  9. Bergemann, M., & Gehren, T. 2008, A&A, 492, 823 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Bergemann, M., Pickering, J. C., & Gehren, T. 2010, MNRAS, 401, 1334 [CrossRef] [Google Scholar]
  11. Caffau, E., Bonifacio, P., Starkenburg, E., et al. 2017, Astron. Nachr., 338, 686 [NASA ADS] [CrossRef] [Google Scholar]
  12. Campbell, S. W., & Lattanzio, J. C. 2008, A&A, 490, 769 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Campbell, S. W., Lugaro, M., & Karakas, A. I. 2010, A&A, 522, L6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Chen, B., Vergely, J. L., Valette, B., & Carraro, G. 1998, A&A, 336, 137 [NASA ADS] [Google Scholar]
  15. Choplin, A., Siess, L., & Goriely, S. 2021, A&A, 648, A119 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Clarkson, O., Herwig, F., & Pignatari, M. 2018, MNRAS, 474, L37 [Google Scholar]
  17. Cohen, J. G., Christlieb, N., Qian, Y.-Z., & Wasserburg, G. J. 2003, ApJ, 588, 1082 [Google Scholar]
  18. Côté, B., Denissenkov, P., Herwig, F., et al. 2018, ApJ, 854, 105 [Google Scholar]
  19. Cowan, J. J., & Rose, W. K. 1977, ApJ, 212, 149 [Google Scholar]
  20. Cristallo, S., Straniero, O., Gallino, R., et al. 2009a, ApJ, 696, 797 [NASA ADS] [CrossRef] [Google Scholar]
  21. Cristallo, S., Piersanti, L., Straniero, O., et al. 2009b, PASA, 26, 139 [Google Scholar]
  22. Cristallo, S., Piersanti, L., Straniero, O., et al. 2011, ApJS, 197, 17 [NASA ADS] [CrossRef] [Google Scholar]
  23. Cristallo, S., Straniero, O., Piersanti, L., & Gobrecht, D. 2015, ApJS, 219, 40 [Google Scholar]
  24. Cristallo, S., Karinkuzhi, D., Goswami, A., Piersanti, L., & Gobrecht, D. 2016, ApJ, 833, 181 [Google Scholar]
  25. Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2003, VizieR Online Data Catalog: II/246 [Google Scholar]
  26. Demarque, P., Woo, J.-H., Kim, Y.-C., & Yi, S. K. 2004, ApJS, 155, 667 [NASA ADS] [CrossRef] [Google Scholar]
  27. Denissenkov, P. A., & Tout, C. A. 2003, MNRAS, 340, 722 [NASA ADS] [CrossRef] [Google Scholar]
  28. Denissenkov, P. A., Herwig, F., Battino, U., et al. 2017, ApJ, 834, L10 [Google Scholar]
  29. Denissenkov, P. A., Herwig, F., Woodward, P., et al. 2019, MNRAS, 488, 4258 [Google Scholar]
  30. Doherty, C. L., Gil-Pons, P., Siess, L., Lattanzio, J. C., & Lau, H. H. B. 2015, MNRAS, 446, 2599 [Google Scholar]
  31. Frebel, A. 2018, Ann. Rev. Nucl. Part. Sci., 68, 237 [Google Scholar]
  32. Fujimoto, M. Y., Ikeda, Y., & Iben, I. Jr. 2000, ApJ, 529, L25 [NASA ADS] [CrossRef] [Google Scholar]
  33. Gaia Collaboration (Katz, D., et al.) 2018, A&A, 616, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Gehren, T., Liang, Y. C., Shi, J. R., Zhang, H. W., & Zhao, G. 2004, A&A, 413, 1045 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Goswami, A., Aoki, W., Beers, T. C., et al. 2006, MNRAS, 372, 343 [Google Scholar]
  36. Goswami, A., Aoki, W., & Karinkuzhi, D. 2016, MNRAS, 455, 402 [NASA ADS] [CrossRef] [Google Scholar]
  37. Goswami, P. P., Rathour, R. S., & Goswami, A. 2021, A&A, 649, A49 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Hampel, M., Stancliffe, R. J., Lugaro, M., & Meyer, B. S. 2016, ApJ, 831, 171 [Google Scholar]
  39. Hampel, M., Karakas, A. I., Stancliffe, R. J., Meyer, B. S., & Lugaro, M. 2019, ApJ, 887, 11 [Google Scholar]
  40. Herwig, F. 2000, A&A, 360, 952 [NASA ADS] [Google Scholar]
  41. Herwig, F., Langer, N., & Lugaro, M. 2003, ApJ, 593, 1056 [NASA ADS] [CrossRef] [Google Scholar]
  42. Herwig, F., Pignatari, M., Woodward, P. R., et al. 2011, ApJ, 727, 89 [Google Scholar]
  43. Hidalgo, S. L., Pietrinferni, A., Cassisi, S., et al. 2018, ApJ, 856, 125 [Google Scholar]
  44. Hill, V., Barbuy, B., Spite, M., et al. 2000, A&A, 353, 557 [NASA ADS] [Google Scholar]
  45. Hollowell, D., Iben, I. Jr., & Fujimoto, M. Y. 1990, ApJ, 351, 245 [NASA ADS] [CrossRef] [Google Scholar]
  46. Iwamoto, N., Kajino, T., Mathews, G. J., Fujimoto, M. Y., & Aoki, W. 2004, ApJ, 602, 377 [Google Scholar]
  47. Jones, S., Ritter, C., Herwig, F., et al. 2016, MNRAS, 455, 3848 [Google Scholar]
  48. Jonsell, K., Barklem, P. S., Gustafsson, B., et al. 2006, A&A, 451, 651 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  49. Karakas, A. I., & Lattanzio, J. C. 2014, PASA, 31 [CrossRef] [Google Scholar]
  50. Lau, H. H. B., Stancliffe, R. J., & Tout, C. A. 2009, MNRAS, 396, 1046 [CrossRef] [Google Scholar]
  51. McWilliam, A. 1998, AJ, 115, 1640 [NASA ADS] [CrossRef] [Google Scholar]
  52. Noguchi, K., Aoki, W., Kawanomoto, S., et al. 2002, PASJ, 54, 855 [NASA ADS] [CrossRef] [Google Scholar]
  53. Piersanti, L., Cristallo, S., & Straniero, O. 2013, ApJ, 774, 98 [NASA ADS] [CrossRef] [Google Scholar]
  54. Qian, Y. Z., & Wasserburg, G. J. 2003, ApJ, 588, 1099 [Google Scholar]
  55. Ramírez, I., & Meléndez, J. 2005, ApJ, 626, 465 [Google Scholar]
  56. Schuler, S. C., Margheim, S. J., Sivarani, T., et al. 2008, AJ, 136, 2244 [NASA ADS] [CrossRef] [Google Scholar]
  57. Smith, V. V., Coleman, H., & Lambert, D. L. 1993, ApJ, 417, 287 [CrossRef] [Google Scholar]
  58. Sneden, C. A. 1973, PhD Thesis, The University of Texas at Austin [Google Scholar]
  59. Stancliffe, R. J., Dearborn, D. S. P., Lattanzio, J. C., Heap, S. A., & Campbell, S. W. 2011, ApJ, 742, 121 [Google Scholar]
  60. Straniero, O., Gallino, R., & Cristallo, S. 2006, Nucl. Phys. A, 777, 311 [NASA ADS] [CrossRef] [Google Scholar]
  61. Yang, G.-C., Liang, Y.-C., Spite, M., et al. 2016, Res. Astron. Astrophys., 16, 19 [Google Scholar]
  62. Yong, D., Norris, J. E., Bessell, M. S., et al. 2013, ApJ, 762, 26 [Google Scholar]

Appendix A: Line list

Table A.1.

Equivalent widths (in mÅ) of lines used for the calculation of elemental abundances.

All Tables

Table 1.

Derived atmospheric parameters of our programme star, and literature values.

Table 2.

Elemental abundances in HE 1005–1439.

Table A.1.

Equivalent widths (in mÅ) of lines used for the calculation of elemental abundances.

All Figures

thumbnail Fig. 1.

Sample spectra of the programme star in the 4055–4065 Å (upper panel) and 5160–5190 Å (bottom panel) wavelength regions.

In the text
thumbnail Fig. 2.

Spectral synthesis fits of C2 bands around 5165 Å (top panel) and 5635 Å (bottom panel). The dotted lines indicate the synthesised spectra and the solid lines indicate the observed spectra. Two alternative synthetic spectra are shown corresponding to Δ[C/Fe] = +0.3 (long-dashed line) and Δ[C/Fe] = −0.3 (short-dashed line).

In the text
thumbnail Fig. 3.

Spectral synthesis fits (dotted curves) of the C2 features around 4740 Å obtained with the adopted C abundance and 12C/13C value (dotted curve). The observed spectrum is shown by a solid curve. Two alternative fits with 12C/13C ∼ 1 (short-dashed line) and 12 (long-dashed line) are shown to illustrate the sensitivity of the line strengths to the isotopic carbon abundance ratios.

In the text
thumbnail Fig. 4.

Filled red circles and filled black squares respectively represent literature CEMP-r/s and CEMP-s stars compiled by Goswami et al. (2021), and the blue star represents the programme star. Panel a: grid formed by the dotted red lines bound by 0.0 < [La/Eu] < 0.6 and 0.0 < [Ba/Eu] < 1.0 indicates the region defined for CEMP-r/s stars by Goswami et al. (2021). The grid formed by the black dashed lines bound by 0.5 < [La/Eu] < 0.7 represents the region where [Eu/Fe] > 1.0 classifies the stars as CEMP-r/s and [Eu/Fe] < 1.0 classifies the stars as CEMP-s. Panel b: red dashed line at [hs/ls] = 1.06 and the black dashed line at [hs/ls] = 0.65 represent the peaks of [hs/ls] for CEMP-r/s and CEMP-s, respectively, as shown by Goswami et al. (2021) in Fig. 13a.

In the text
thumbnail Fig. 5.

Examples of theoretical model fits with the observed abundances of the star. In panels a, b, and c, the points with error bars indicate the observed abundances. (a) s-process model fits for M = 1.3 M & 2.0 M. (b) i-process model fits for n = 1012 cm−3 and 1014 cm−3. (c) Best-fit from parametric-model. (d) Residual plot.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.