EDP Sciences
Home All issues Volume 592 (August 2016) A&A, 592 (2016) A29 Full HTML
Free Access
Issue
A&A
Volume 592, August 2016
Article Number A29
Number of page(s) 16
Section Stellar structure and evolution
DOI https://doi.org/10.1051/0004-6361/201628540
Published online 12 July 2016

© ESO, 2016

1. Introduction

Currently observable metal-poor stars ([Fe/H] 1 < −2) are relatively unevolved low-mass objects that formed within the first couple of billion years after the Big Bang. Because these stars have not modified their surface composition by internal nucleosynthesis, they are expected to carry the signature of the chemical evolution of these early epochs, providing us with the means to study this long-gone era. The number of known metal-poor stars in our Galaxy has exploded thanks to large-scale photometric and spectroscopic surveys such as the HK survey (Beers et al. 1985, 1992), and more recently the Hamburg/ESO survey (Christlieb et al. 2001, 2008), the Sloan Digital Sky Survey (SDSS; e.g. York et al. 2000; Ahn et al. 2012, 2014), and its sub-survey, the Sloan Extension for Galactic Understanding and Exploration (SEGUE; Yanny et al. 2009). A common finding of these surveys is that a substantial fraction of all metal-poor stars are relatively carbon-rich with [C/Fe] ≳ 1.0. This fraction is around 10% at [Fe/H] ≈ −2 and increasing towards lower metallicities (e.g. Lucatello et al. 2006; Carollo et al. 2012; Lee et al. 2013; Placco et al. 2014).

These so-called carbon-enhanced metal-poor (CEMP) stars are further classified into CEMP -no, -s, -r, and -r/s sub-classes depending on whether they show enhancements of elements produced by slow (s) or rapid (r) neutron-capture nucleosynthesis. Most CEMP stars with [Fe/H] > −3 display significant s-process element enrichment and are classified either as CEMP-s ([Ba/Fe] > 1 and [Ba/Eu] > 0.5) or CEMP-r/s ([Ba/Fe] > 1 and 0 < [Ba/Eu] < 0.5) stars (Beers & Christlieb 2005)1. While the abundance patterns of these stars have been linked to nucleosynthesis occurring in asymptotic giant branch (AGB) stars (Herwig 2005; Sneden et al. 2008; Bisterzo et al. 2010, 2011; Lugaro et al. 2012), most s-process-rich CEMP stars are not luminous enough to be AGB stars. However, results from radial velocity monitoring are consistent with all of them being in binaries which is not the case for the other sub-classes (Lucatello et al. 2005; Starkenburg et al. 2014; Hansen et al. 2015, 2016a,b). Therefore, these stars are generally thought to be a product of mass transfer from an AGB companion that later became a white dwarf 2. As noted by Lucatello et al. (2005), this would make CEMP-s stars the low-metallicity analogues of CH and Ba stars (McClure & Woodsworth 1990; Jorissen et al. 2016). If this is the case, observations of the secondary can, in principle, be used to infer the nucleosynthesis of low-metallicity AGB stars with initial masses around solar and above, which, although important for galactic and globular cluster chemical evolution, no longer exist in the Local Universe (e.g. Travaglio et al. 1999, 2001; Valiante et al. 2009; Kobayashi et al. 2011; Bisterzo et al. 2014; D’Ercole et al. 2010; Ventura et al. 2014).

A comparison of AGB nucleosynthesis models with abundances of CEMP-s stars is straightforward only if the accreted material remains on the surface of the star. This will certainly not be the case once it evolves off the main sequence and develops a deep convective envelope. But as demonstrated by Stancliffe et al. (2007), this is also unlikely on the main sequence because the higher mean molecular weight of the accreted material should trigger thermohaline mixing (Ulrich 1972; Kippenhahn et al. 1980). Furthermore, gravitational settling of heavier elements could both modify the extent of this mixing and the subsequent evolution of the secondary (Stancliffe & Glebbeek 2008; Thompson et al. 2008; Stancliffe 2009). If the overall effect is to dilute the accreted material by uniformly mixing it throughout some portion of the secondary, a comparison between AGB nucleosynthesis models and abundances of CEMP-s stars is still possible, provided this amount of dilution can be estimated (e.g. as attempted by Bisterzo et al. 2011, 2012). However, the rather impartial (leading to similar dilution of all elements) process of settling will be counteracted by the highly selective process of radiative levitation (also known as radiative accelerations), a process in which the ions in the stellar plasma gain a net outward momentum from absorption of the diffusing photons. Metal-poor stars with masses around 0.8 M⊙ have very shallow convective envelopes, and, in absence of any counteracting processes, large abundance anomalies can result (e.g. Richard et al. 2002a,b). If radiative levitation is important during the post-mass-transfer evolution of CEMP-s stars, the interpretation of their abundances in the context of AGB nucleosynthesis gets considerably more complicated. In this paper we model the main sequence evolution of CEMP-s stars, including the effect of radiative levitation, to investigate whether this is the case.

We focus primarily on the evolution of carbon and iron surface abundances. Abundances of s-process elements are not modelled. However, levitation is expected to have a much greater impact on iron than on carbon (Gonzalez et al. 1995; Seaton 1997, 2007). We therefore attempt to constrain its overall importance for CEMP-s stars by investigating these two elements.

2. Methods

We use the stellar evolution code STARS originally written by Eggleton (1971, 1972), Eggleton et al. (1973) and since improved by many authors (e.g. Pols et al. 1995; Stancliffe & Eldridge 2009). The version used in this work tracks the abundances of the nuclear species 1H, 3He, 4He, , 14N, 16O, 20Ne, 24Mg, 28Si, and 56Fe, the last three of which were previously not tracked in detail. The mass fraction Xi of each species i is governed by an advection-diffusion equation: (1)where the first term on the right-hand side accounts for convective mixing, thermohaline mixing, and concentration diffusion (Dmix is the sum of the individual diffusion coefficients Dconv, Dμ, and Di, respectively), the second term describes the net effect from atomic diffusion, and the last term, Ri, accounts for nuclear processing3. Inside convective regions Dconv is obtained from the mixing length theory (MLT; Böhm-Vitense 1958) using a solar-calibrated value of αMLT = 2.0 that is fully consistent with the models presented in this paper (see Sect. 5 for details). Near the convective boundaries the convective mixing coefficient takes the form from Eggleton (1972) for numerical stability reasons. The diffusion coefficient for thermohaline mixing is taken from Denissenkov (2010) assuming a finger length-to-diameter ratio of 0.5 as constrained by their numerical simulations. This assumption of more blob-like than finger-like structures is in accord with Kippenhahn et al. (1980) and results in relatively inefficient thermohaline mixing.

Following Stancliffe & Glebbeek (2008) we treat atomic diffusion in a trace approximation, which allows the diffusion velocity of elements other than hydrogen to be written as (Talon 2008) (2)where g and gr are the gravitational and radiative acceleration, respectively; k is the Boltzmann constant; T is temperature; μ is the mean molecular weight of the stellar plasma; is the molecular weight of element i with an atomic mass mi and a mean charge ; and αT,i is the thermal diffusion coefficient. The velocity of hydrogen follows from mass conservation: XHwH = − ∑ i ≠ HXiwi. The diffusion coefficients Di and αT,i are taken from Paquette et al. (1986).

We simulate the accretion of AGB ejecta by adding mass of a given composition to our models. We fix the accretion rate to 10-6 M⊙yr-1 following Stancliffe et al. (2007), and the accreted composition to the average composition of the ejecta from the models of Lugaro et al. (2012). These yields, together with the zero-age main sequence (ZAMS) abundances (Asplund et al. 2009), are given in Table 1. Mass loss is not included in our models until Sect. 5.2.

Table 1

Chemical composition of the secondaries on the zero-age main sequence (ZAMS; abundance distribution from Asplund et al. 2009, scaled to Z = 10-4) and of the ejecta from the AGB models of Lugaro et al. (2012).

2.1. Opacity and radiative accelerations

We compute the radiative acceleration of each element using the monochromatic data from version 3.3 of the Opacity Project (OP) database (Badnell et al. 2005; Seaton 2007). The OP data consists of cross-sections σi(u ≡ hν/kT) and electron scattering corrections ai(u) for 17 chemical elements between H and Ni in a temperature range between log 10T = 3.5 and log 10T = 8.0. The monochromatic data are used to compute the Rosseland mean opacity κR and, for each element i, a factor γi that is proportional to the radiative acceleration of the respective element:

Here Ni is the number fraction of element i, and v(u) is the OP frequency variable (5)With these quantities the radiative accelerations are given by (6)where lr is the radiative luminosity, and c is the speed of light.

The OP team have created the OPserver module (Mendoza et al. 2007) which is intended to facilitate the computation of accelerations in stellar evolution calculations. We have made some changes to this module in coupling it to the STARS code. First, we have made it possible to store multiple opacity (and acceleration) tables in memory at the same time. This requires computing the opacity corresponding to a given chemical composition only once. When this composition is encountered during evolution, one only needs to interpolate to the required temperature and density in the corresponding table. The second modification is the same as made by Hu et al. (2011) in their incorporation of the OPserver module in a version of the STARS code – instead of calculating the acceleration for multiple relative abundances of a given element, we calculate the accelerations of all elements in a given mixture. Finally, we have added a routine that computes the mean charge of each element from the OP data. These charges are computed on the same temperature and density grid as the opacity and are used to calculate the molecular weights of the elements.

Calculating both the opacity and the accelerations from the monochromatic OP data makes the models self-consistent in that changes in relative abundances modify the structure of the star through the opacity, which, in turn, changes the accelerations. Unfortunately, the OP opacities do not include any contribution from conduction, which becomes important after the main sequence when the central regions of the star become increasingly degenerate. Since in many cases we follow the evolution all the way up the giant branch, we use the opacity tables introduced in the code by Eldridge & Tout (2004) for regions hotter than log 10T = 7.3. These tables are based on the OPAL opacities (Iglesias & Rogers 1996), supplemented by the low-temperature opacities of Alexander & Ferguson (1994) and the conductive opacities of Hubbard & Lampe (1969) and Canuto (1970). While switching to the tabulated OPAL opacities means that changes in the relative abundances no longer modify the structure at high temperatures, by then the effects of atomic diffusion have already started to disappear because of the first dredge-up (FDU), and none of our results depend on this choice (see Sect. 5).

2.2. Grid selection

Our simulations cover a range of primary masses M1, accreted masses ΔM, and initial secondary masses M2,i (or, equivalently, final masses M2,f). In this work we consider those systems that are the most probable in the synthetic populations computed by Abate et al. (2015). According to their work, typical masses are M1 ≃ 0.9–1.25 M⊙, M2,f ≃ 0.8–0.9 M⊙, and ΔM ≃ 0.05–0.2 M⊙. These accreted masses and final masses of the secondaries are larger than considered in a related earlier study by Stancliffe & Glebbeek (2008). Therefore, we also consider some systems with smaller ΔM values, namely 0.001 and 0.01 M⊙.

In summary, we evolve stellar models with initial masses of 0.60, 0.65, 0.70, 0.75, and 0.80 M⊙ and metallicity Z = 10-4 ([Fe/H] = −2.14) starting from the pre-main-sequence. At the ages listed in Table 1, somewhere between 0.001 and 0.2 M⊙ of material of the corresponding composition is added to the models at a rate of 10-6 M⊙yr-1 yielding CEMP-s stellar models with masses between 0.8 and 0.95 M⊙. These models are evolved up to the core helium flash or an age of 16 Gyr, whichever comes first.

3. Results

Two sets of models were initially evolved: in one set only thermohaline mixing, gravitational settling, and thermal diffusion were active; in the other, radiative levitation was also included. Table A.1 lists some properties of these systems, including the [C/Fe] ratio at the surface at key points of the evolution: after thermohaline mixing, at the point where the convective envelope is smallest in mass (near the turn-off), and after first dredge-up.

3.1. An illustrative model sequence

To understand how the evolution of surface abundances is influenced by the different physical processes included in our simulations, let us consider a particular model sequence in detail. Figure 1 illustrates the case of a secondary with an initial mass of 0.75 M⊙ that accretes 0.05 M⊙ of material from a 1.25 M⊙ primary. Multiple stages of evolution can be distinguished.

thumbnail Fig. 1

Evolution and abundances of a M2,i = 0.75 M⊙ secondary accreting ΔM = 0.05 M⊙ of material from a M1 = 1.25 M⊙ primary. The labels highlight specific parts of the evolution discussed in the text. The model sequences with and without radiative levitation overlap at this scale of the HRD. Interior abundance profiles are shown at a few of the stages indicated in the upper panels: before mass transfer (“1”, solid); before thermohaline mixing (“4”, long-dashed); after thermohaline mixing (“5”, dotted); post-mass-transfer main sequence (“6”, dot-dashed); minimum of convective envelope mass (“7”, dot-dot-dashed); during first dredge-up (“9”, short-dashed). The vertical lines in the lower panels indicate the position of the base of the convective envelope at the respective time. The interior profiles of carbon with and without levitation nearly coincide, and only the case with levitation is shown.
Open with DEXTER

Prior to mass transfer the secondary slowly evolves as a 0.75 M⊙ main sequence star (the part of the evolution labelled “1” in Figs. 1a and b). During this stage gravitational settling dominates, and the abundance of every element other than hydrogen decreases at the surface.

At t = 3.06 Gyr mass transfer begins, and the surface composition quickly becomes equal to that of the accreted material (“2”). During the accretion the star becomes hotter and more luminous. This is common for many system configurations in which the Kelvin-Helmholtz timescale of the secondary becomes comparable to the accretion timescale. In some models the effective temperature and luminosity can reach values as high as 9500 K and 30 L⊙, respectively. But once accretion stops (“3”), both luminosity and temperature rapidly drop, resulting in loops in the Hertzsprung-Russell diagram (HRD) as the star settles back on the main sequence. Generally, these loops are more characteristic of secondaries with larger initial masses.

Shortly after accretion stops, the accreted material starts to mix with the original material of the secondary as a result of the thermohaline instability (“4”). As shown by Fig. 1c, some of the interior is already mixed by the time the surface abundances change. The mixing takes only a few hundred million years (about 150 Myr in this case) and is over before the star has settled back on the main sequence (“5”). Ultimately, the surface carbon abundance is reduced by about 0.8 dex (regardless of radiative levitation), whereas the iron abundance is barely affected because it is virtually the same in the original and accreted compositions.

Over the rest of the post-mass-transfer main sequence lifetime the abundances are again modified by atomic diffusion (“6”). At first gravitational settling prevails over radiative levitation, and the surface becomes increasingly hydrogen-rich as all heavier elements settle out of the surface convection zone. As the star nears the turn-off, this convection zone becomes ever more superficial, and radiative effects become increasingly important (Fig. 1d). Once an element’s radiative velocity at the base of the convection zone exceeds its settling velocity, the surface abundance of this element increases. This is typically the case with iron. In contrast, if an element’s settling velocity is always greater than its radiative velocity, the surface abundance of this element continues to decrease (although less so than in the case when radiative effects are ignored). This is always the case with helium and carbon. The behaviour of other elements is not readily predicted because of the non-monotonic shape of the radiative accelerations (as a function of temperature) and the outward movement of the base of the envelope (decreasing temperature at the base). Therefore, the abundance of most elements alternates between increasing at those times when the radiative velocity exceeds the settling velocity and decreasing at others (e.g. see Fig. 2 of Richard et al. 2002b).

Eventually, the abundance anomalies, i.e. their values relative to those after thermohaline mixing, reach their maxima (“7”). At this stage the difference between the two sets of models is greatest – compared to the abundances after thermohaline mixing, in models without levitation the abundances of all elements are reduced (dotted lines in Fig. 1b), whereas in models with levitation this is not always the case (solid lines in Fig. 1b; in both cases only carbon and iron are shown for clarity), and readily levitated elements can be over-abundant. The anomalies are maximal shortly after the turn-off when the convective envelope is smallest (Figs. 1c and d)4. This occurs at the same time in models with and without levitation.

Next, as the convective envelope grows in mass, the material in it is mixed with that of the immediately adjacent, previously radiative layers. In models without diffusion no change in surface abundances would occur until the envelope reached depths where CN cycling had occurred (i.e. at first dredge-up). With diffusion, however, the composition of the envelope is different from that of the radiative layers below, and therefore the effect of the deepening of the envelope is to first undo all the work done by diffusion (“8”). When the envelope mass has reached a few thousandths of a solar mass, all surface evidence of atomic diffusion has been erased, and the abundances are similar to those after thermohaline mixing (“9”).

First dredge-up homogenizes the composition in the layers above a mass coordinate of 0.3–0.35 M⊙ (0.34 M⊙ in this case). What effect this has on the surface abundances depends on how this depth compares to the depth of thermohaline mixing (mthm = 0.37 M⊙ in this case). If thermohaline mixing is not as deep as the maximum depth reached by the envelope at FDU, the accreted material is further diluted with the original material of the secondary. Otherwise, most abundances do not change. However, some of the accreted carbon will then have been converted into nitrogen. As shown by Stancliffe et al. (2007), during late FDU (which ends at around log L ≈ 1.5; “10”) this nitrogen is dredged up to the surface.

Finally, after the luminosity bump (“11”) 3He-burning reduces the mean molecular weight above the hydrogen-burning shell. Thus, a μ-inversion, which is magnified by the settling of 4He (Michaud et al. 2010), develops between the shell and the receding convective envelope – a situation again unstable to thermohaline mixing. This alters the surface abundance of nitrogen by 0.1 dex at most. The much greater carbon abundance remains essentially unchanged. The abundance change after the bump is much smaller than found by Stancliffe et al. (2009), because the thermohaline mixing coefficient in this work is about 103 times smaller.

This model sequence illustrates the role each physical process plays in all models with atomic diffusion. We see that diffusion modifies the surface composition on the main sequence both before and after mass transfer. This modification is greatest around the turn-off, when the convective envelope is shallowest (point “7” in Fig. 1). We now turn to discussing the expected abundance changes for all CEMP-s stars in this evolutionary stage.

3.2. Abundance anomalies near the turn-off

During the main sequence the mass of the convective envelope, Menv, of a low-mass star decreases. Therefore, the timescale for atomic diffusion, which is proportional to roughly the square root of Menv (Michaud 1977), also decreases. In nearly all of our CEMP-s models the envelope mass reaches a minimum of less than 10-4 M⊙ around the turn-off. The corresponding timescales are short enough, compared to the nuclear timescale, for atomic diffusion to notably modify the surface composition. Figure 2 summarizes the extent of the abundance variations in our models with diffusion (Table A.1). Specifically, the figure shows the [Fe/H], [C/H], and [C/Fe] abundances at the time when the convective envelope mass reaches the minimum in each of the CEMP-s models. In models with envelope masses always above about 2 × 10-5 M⊙, gravitational settling prevails; however, the abundances are decreased only by up to a factor of two from their values after thermohaline mixing. But in models with even smaller Menv at the turn-off, abundances are modified by a factor of ten or more, and radiative levitation becomes important (Fig. 2a). The model discussed in Sect. 3.1 is close to this limit – its envelope mass has a minimum of about 1.2 × 10-5 M⊙. At this minimum its [C/Fe] is 1.71 when levitation is ignored versus 1.12 when it is included (in both cases down from 1.78 after thermohaline mixing), and [Fe/H] is − 2.74 and − 2.11, respectively.

thumbnail Fig. 2

a) [Fe/H], b) [C/H], and c) [C/Fe] in each of the CEMP-s models at the point where the mass of the convective envelope is smallest, i.e. just before first dredge-up. Models with and without radiative levitation are plotted with black and grey symbols, respectively. All values for envelope masses below 10-7 M⊙ denote upper or lower limits. The [Fe/H] plot also shows the metallicity evolution in a 0.8 M⊙ model sequence with no accretion (solid lines). Similarly, the [C/Fe] plot shows the evolution in a model corresponding to M1 = 1.25 M⊙, M2,i = 0.8 M⊙, and ΔM = 0.05 M⊙. Both sequences without levitation stop during the first dredge-up. The criterion for being classified as a carbon-enhanced metal-poor star ([C/Fe] ≥ 1.0) is from Beers & Christlieb (2005).
Open with DEXTER

The results from many simulations plotted in Figs. 2a and 2c form two sequences corresponding to model sets with and without levitation. As shown by the solid lines, to a decent approximation we can interpret these sequences as describing the abundance evolution in a single simulation as the envelope mass changes. For example, as the envelope mass decreases from 10-5 to 10-6 M⊙, [C/Fe] decreases by about two orders of magnitude in models with levitation because, while carbon continues to settle, iron is levitated. On the other hand, in models without levitation [C/Fe] does not change, because both elements settle at similar rates (Figs. 2a and b). At still smaller envelope masses, [C/Fe] increases because of the low degree of ionization of iron at the base of the envelope. The small mean charge of iron gives a large diffusion coefficient because of the approximately dependence (Paquette et al. 1986) and, consequently, a large settling velocity.

Figure 2 shows that for given post-thermohaline-mixing abundances the abundance evolution prior to FDU can be parametrized as a function of only Menv. But once FDU starts and the envelope deepens, a kind of hysteresis is seen in that the abundances at a given Menv are not the same as they were prior to FDU. This is because diffusion has modified the radiative layers below the envelope in the meantime.

Carbon-enhanced metal-poor stars are distinguished from other metal-poor stars based on their [C/Fe] value. Assuming that atomic diffusion is correctly predicted and no additional mixing processes operate in the radiative regions below the envelope, Fig. 2c implies that they must have envelope masses larger than 10-5 M⊙. Otherwise, they would not be classified as carbon-enhanced. While models without levitation do not have such a limit, the metallicity rapidly decreases below this value because of settling – at Menv ≈ 10-6 M⊙ the surface [Fe/H] ≈ −4, which is much lower than typical of CEMP-s stars.

For a given combination of AGB and CEMP-s star masses (M1 and M2,f, respectively), the convective envelope is deeper in models with larger accreted mass. For example, a 0.8 M⊙ CEMP-s star with an initial mass of 0.6 M⊙ retains a more massive envelope than one with an initial mass of 0.7 M⊙. This is due to the higher average opacity of these stars (more metal-rich stars maintain thicker convective envelopes for the same reason). The difference in Menv can be a factor of 2–10 (depending on M1, M2,f, and the range of M2,i), which can lead to substantially different abundances when the envelopes are small (Fig. 2).

As can be seen from Fig. 2 and Table A.1, diffusion is extremely efficient in many of our models, leading to unrealistic abundance anomalies (e.g. [C/Fe] < −1 or [C/Fe] > 4). In most of our more massive CEMP-s models (M2,f ≥ 0.85 M⊙) diffusion is so efficient that our code is incapable of resolving the steep abundance gradients developing at the base of the envelope, and we are forced to stop the computations before the main-sequence turn-off. Such massive CEMP-s stars are nevertheless probable according to population synthesis calculations (Abate et al. 2015) and would help explain the properties of some CEMP-s RR Lyrae stars (Stancliffe et al. 2013), so we would like to explore their connection to observations. Therefore, we proceed by assuming that diffusion, and possibly thermohaline mixing, is inhibited throughout these stars for one reason or another (leaving open the nature and cause of the underlying mechanism) and evolve two sets of model sequences without atomic diffusion: one set with thermohaline mixing and one without. The results from these simulations are summarized in Table A.2.

These models have some key differences in global properties and surface abundances from the model sequences with atomic diffusion. First, the surface abundances do not change prior to mass transfer. More importantly, after thermohaline mixing has reached equilibrium, no further abundance changes occur until FDU (i.e. between the stages labelled “5” and “9” in Figs. 1a and b). The importance of FDU still depends on the depth of thermohaline mixing, as in models with diffusion. In models without thermohaline mixing the surface abundances do not change until FDU, during which the accreted material is invariably diluted by mixing throughout most of the star (down to a mass coordinate of 0.3–0.35 M⊙).

In agreement with previous studies, models with diffusion are younger (by a few percent) at a given evolutionary stage and have lower effective temperatures at the turn-off than models without diffusion, primarily because of the gravitational settling of helium throughout the star (e.g. Castellani & degl’Innocenti 1999; VandenBerg et al. 2002; Bressan et al. 2012). Our non-accreting models with diffusion are about 150 K cooler than models without diffusion (compare Table A.1 and A.2). Among our CEMP-s models, those with thermohaline mixing but without diffusion are generally between 100 to 300 K hotter than models without both. The latter are cooler because of the high opacity of the outer layers owing to the metal-richness of the accreted material. CEMP-s models with diffusion likely fall somewhere in between, but we cannot do a proper comparison because our massive models with diffusion do not make it to turn-off for numerical reasons.

3.3. Thermohaline mixing

The fraction of the star mixed by thermohaline convection, which we shall call the thermohaline mixing efficiency, essentially depends on the mean molecular weight of the accreted material and its total mass compared to the final mass of the star. The more helium- and metal-rich the accreted material, the greater its molecular weight compared to the initial composition, and the greater the portion of the star that gets mixed. From the last column of Table 1 we therefore expect that, for a given amount of accreted material, thermohaline mixing should be most efficient when that material comes from a primary of 1.5 M⊙, and least efficient when it comes from a primary of 1.0 M⊙, which is indeed the case (Fig. 3).

thumbnail Fig. 3

Thermohaline mixing efficiency (fraction of the star that is mixed) as a function of the ratio between accreted mass and final mass for different primaries. Black and grey symbols correspond to models with and without diffusion, respectively. Each symbol represents a unique combination of M1, M2,f, and ΔM.
Open with DEXTER

Furthermore, the greater the amount of the high-μ material, the deeper the mixing must be for the μ-gradient to be removed. If an amount ΔM of AGB ejecta with a mean molecular weight μa is mixed with M2,f−ΔM−mthm of the unpolluted material with an average molecular weight μi (<μa) before the μ-gradient is removed, a mixed region M2,f−mthm with molecular weight μf results (Fig. 4). Equating the states before and after mixing, one gets that the removal of the μ-gradient implies a linear relationship between the accreted-to-final mass ratio and mixing efficiency. Indeed, a linear relationship is a reasonable approximation in the range 0.05 ≲ ΔM/M2,f ≲ 0.2 of accreted-to-final mass ratios (Fig. 3). However, higher and lower ratios of ΔM/M2,f require special consideration.

thumbnail Fig. 4

Schematic illustration of the structure of the star a) before and b) after thermohaline mixing. The shaded area around mthm and below indicates a region in which the mean molecular weight has been raised as a result of nuclear processing.
Open with DEXTER

First, a high ΔM/M2,f corresponds to a case where a lot of mass is transferred to a low-mass star, which implies that the mass transfer takes place when the secondary is still nearly on the ZAMS. Two outcomes are possible. If the accreted material has a sufficiently high molecular weight compared to the molecular weight throughout the star, thermohaline mixing affects the whole star, and the composition is nearly homogenized. For example, this occurs when a 0.65 M⊙ secondary accretes 0.2 M⊙ of material from a 1.5 M⊙ primary. On the other hand, if the accreted material has a lower molecular weight than some central region of the star, increasing the accreted mass will not lead to a much deeper mixing because of the steepness of the μ-gradient near the center. The mixing efficiency can therefore decrease for higher ΔM/M2,f (although at this point almost all of the star will be mixed, anyway).

Second, in models with diffusion low accreted-to-final mass ratios lead to relatively inefficient thermohaline mixing. Contrary to models without diffusion, where the molecular weight is constant outside of nuclear burning regions, in these models there is a stabilizing μ-gradient throughout the star owing to gravitational settling. This presents a “μ-barrier” to thermohaline mixing that must be overcome for mixing to happen (Thompson et al. 2008). If only a tiny amount of material is accreted (ΔM ≲ 0.001 M⊙), thermohaline mixing can be almost completely inhibited (left bottom corner of Fig. 3). Nevertheless, even in these cases the surface carbon content is depleted by a factor of two or more because the mass of the mixed region, M2,f−mthm, is always greater than ΔM.

For ΔM> 0.01 M⊙ the μ-barrier is largely overwhelmed. However, the mixing is nevertheless slightly less efficient than in models without diffusion, even for large amounts of accretion. Similar conclusions were reached by Stancliffe & Glebbeek (2008). Overall, the surface carbon abundance is typically reduced by somewhere between 0.3 to 1 dex, depending on the relative amount of the accreted material and its molecular weight. Since radiative accelerations have almost no influence on the molecular weight profile deep in the star, they have almost no influence on the efficiency of thermohaline mixing.

4. Comparison with observations

We have presented four sets of models of CEMP-s stars. Two of the sets comprise models with thermohaline mixing and atomic diffusion (one set with, one without radiative levitation). The other two sets comprise models without diffusion (one set with, one without thermohaline mixing). We now compare the four sets of models with observations of CEMP stars.

The largest data set of Galactic metal-poor stellar spectra currently available is that from the SDSS and SEGUE surveys. Lee et al. (2013) used the spectra collected by SDSS/SEGUE to derive the stellar parameters and carbon abundances ([C/Fe]) in close to 250 000 stars, around 10 000 of which have [Fe/H] < −2. We now use this homogeneous metal-poor sample (Lee et al., priv. comm.) to compare the observed carbon abundances with our models.

Comparing the observed [C/Fe] abundances with models must be done with care. Use of a fixed metallicity ([Fe/H]) range might not be adequate because of the diffusion of iron (Figs. 5a and b). Only the 0.8 M⊙ models with levitation generally predict [Fe/H] to remain within a factor of about two from the initial value ([Fe/H] = −2.14) throughout the evolution, whereas the models without levitation have [Fe/H] ≲ −2.5 near the turn-off. Meanwhile, most of the 0.85 M⊙ models have [Fe/H] < −3 (without levitation) or [Fe/H] ≳ −1 (with levitation).

thumbnail Fig. 5

Evolution of [Fe/H] (upper panels), [C/H] (middle panels) and [C/Fe] (lower panels) in CEMP-s models of 0.8 M⊙ (left panels) and 0.85 M⊙ (right panels) with diffusion. Thick lines are models with radiative levitation, whereas thin lines are those without. Empty circles mark an age of 10 Gyr, and filled circles mark an age of 13.8 Gyr, i.e. the part of the track between the circles covers the expected age range of CEMP-s stars. The small, grey circles show the metal-poor stars observed in the Sloan Digital Sky Survey in the metallicity range − 2.5 ≤ [Fe/H] ≤ −2.0 (Lee et al. 2013).
Open with DEXTER

It is safer to first consider [C/H] because the largest [C/H] values should be close to 0.5 dex, independent of metallicity. In the metallicity range typical of CEMP-s stars (−2.5 ≤ [Fe/H] ≤ −2.0), this is indeed the case for subgiants (log g ≲ 3.7), but very few turn-off stars (log g ≈ 4) seem to have [C/H] > 0 (Figs. 5c and d). Is this difference in the maximum [C/H] values between the two groups in the observational data (Δ [C/H] ≈ 0.5) evidence of gravitational settling of carbon in CEMP stars? And does the similar 0.5 dex difference seen in the [C/Fe] data (Figs. 5e and f) then imply that iron is levitated just enough that its abundance stays roughly constant throughout the evolution? This seems rather unlikely because then the carbon-normal population should plausibly also have lower [C/H] (and [C/Fe]) values near the turn-off, which is not the case. On the contrary, the observations suggest that the carbon abundance in carbon-normal metal-poor stars is increasing on the main sequence until they reach the turn-off. This is exactly the opposite behaviour that atomic diffusion predicts! Moreover, there is no obvious candidate for a physical process that could cause the surface carbon abundance to increase on the main sequence.

The carbon-normal metal-poor stars listed in the Stellar Abundances for Galactic Archeology database (SAGA; Suda et al. 2008, 2011) do not seem to show a similar behaviour of increasing carbon abundance on the main sequence (C. Abate, priv. comm.), although the relatively small number statistics (23 stars with −2.5 ≤ [Fe/H] ≤ −2.0, [C/Fe] < 1 and log g > 4.0) and the heterogeneity of the data could perhaps hide such a trend. Unfortunately, whether there is a real difference in the upper [C/H] and [C/Fe] values between turn-off stars and subgiants remains unclear.

Most of the models are at odds with the Lee et al. (2013) data, regardless of whether the abundance differences between turn-off stars and subgiants are caused by atomic diffusion. For example, while some of the 0.8 M⊙ models predict abundances that are consistent with the observations, they only do so at very late times, i.e. at ages exceeding the age of the Universe (13.8 Gyr; Hinshaw et al. 2013). At earlier times stars with low-mass AGB companions (which thus accreted mass later) and/or with low initial masses (which were thus less evolved at the point of mass transfer) are still relatively unevolved. They should be observable as carbon-rich low-luminosity (log g ≳ 4.1) objects. But such objects are conspicuous by their absence in the Lee et al. (2013) results (at all metallicities; their Fig. 6). Since there are plenty of carbon-normal low-luminosity stars in the data, it is difficult to imagine how this could be a selection effect, particularly since a few carbon-rich dwarfs have been found in the SDSS data in detailed abundance studies (Aoki et al. 2008; Behara et al. 2010)5. By contrast, the 0.85 M⊙ CEMP-s models are sufficiently evolved, but diffusion is so efficient that unrealistic abundances (e.g. [C/Fe] < −0.5 with levitation or [C/Fe] > 3.5 without levitation) are predicted in nearly all such stars around the turn-off (Fig. 5f). Clearly, in these more massive CEMP-s stars some physical process must be countering atomic diffusion, at least near their surface.

The disagreement between observations and models concerning the existence of low-luminosity carbon-rich stars has little to do with atomic diffusion – even if we identify a process that neatly counteracts diffusion near the surface, models will still predict many carbon-rich unevolved objects. In fact, if this process were to counteract diffusion throughout the star, the tension with observations would increase, because diffusion starves the core of fuel and accelerates the evolution, making the star spend less time on the main sequence. Perhaps, the SDSS sample indicates that the mass ratio (q ≡ M2,i/M1) in CEMP-s progenitor systems is biased towards unity, contrary to the common assumption of a flat distribution. If q is close to one, the mass transfer occurs relatively late, giving the secondary more time to evolve before it becomes a CEMP-s star.

Overall, models without diffusion are better able to envelop the observational data (Fig. 6). For example, CEMP-s models that have accreted less carbon-rich material (here coming from 1 M⊙ primaries) have [C/Fe] < 1.5, with [C/Fe] increasing with accreted mass. Models that have accreted more carbon-rich material (from primaries with masses 1.25 M⊙, 1.5 M⊙) can have [C/Fe] ≳ 2 throughout the evolution. Models with thermohaline mixing seem to be preferred, because no sharp change in [C/Fe] at FDU is evident from the data. Lower-mass models (M2,f ≈ 0.8 M⊙) without diffusion would predict similar abundance evolution (Stancliffe et al. 2007; Stancliffe & Glebbeek 2008). However, most of them, coming from systems with relatively low-mass AGB companions compared to earlier works, would be consistent with the data only for t > 13.8 Gyr.

thumbnail Fig. 6

Similar to Figs. 5e and f, but here the lines correspond to CEMP-s models of 0.9 M⊙ without diffusion. Thick lines are models with thermohaline mixing, whereas thin lines are those without.

Open with DEXTER

Last, there are some objects with [C/Fe] ≳ 2.5 whose surface gravities (log g ≲ 3) imply that they are close to the end of FDU if not past it. How might we explain such objects, assuming that they were polluted by an AGB companion? Since carbon is not produced in the star, it is hard to imagine how the surface carbon abundance could be above the one in the accreted material. This limits the possible primaries to those that are able to produce at least this much carbon. From the models of Lugaro et al. (2012) these are AGB stars with 1.25 ≲ M1/M⊙ ≲ 3. Lower-mass AGB stars do not produce enough carbon, and higher-mass stars convert the carbon into nitrogen in the lower part of their convective envelope (a process known as hot-bottom-burning). Moreover, a large amount of mass must be transferred, because the combined carbon reduction from thermohaline mixing and FDU can only be about 0.5 dex at most (the maximum [C/Fe] given by the AGB models of Lugaro et al. 2012, is about + 3.2). This implies an accreted-to-initial mass ratio, ΔM/M2,i, between about 0.25, if all the mixing occurs during FDU (where the envelope mass grows to about 0.5 M⊙), and 0.5, if thermohaline mixing is efficient (as it will be in the absence of some inhibitory process because of how unevolved the progenitor of the CEMP-s star must be). The respective accreted-to-final mass ratios, ΔM/M2,f, are between 0.2 and 0.35. Hence, the progenitor systems of the most carbon-rich evolved stars must be small mass-ratio binaries in which a lot of mass has been transferred to the secondary. It may be difficult to account for such stars without also predicting too many low-luminosity carbon-rich stars from cases where less mass has been transferred.

5. Discussion

The results presented in this paper may lead one to wonder whether the abundance anomalies predicted by our diffusion models are overestimated. We have run multiple tests to address this concern. First, we have tried to reproduce the results of Richard et al. (2002b). In particular, we have compared the abundance evolution in a M = 0.8 M⊙ model with an initial composition taken from their Table 1 (Z = 1.7 × 10-4; [Fe/H] = −2.31). We obtain good agreement in terms of the temperature and luminosity at the turn-off, although our model is longer-lived by about 0.5 Gyr. The abundance anomalies from settling and/or levitation of He, C, N, and Fe agree within 0.1 dex. For other elements the abundances differ by 0.3–0.4 dex. Given that in our model the convective envelope mass is smaller by about a factor of two throughout most of the evolution (and the minimum size of the convective envelope in our model is only about 3.7 × 10-6 M⊙, whereas Richard et al. (2002b) get about 2.5 × 10-5 M⊙), these differences are plausible. Judging from their Fig. 2, a smaller envelope mass in their model would lead to greater over-abundances of O, Ne, and Mg, and a smaller over-abundance of Si. All of these changes would reduce the discrepancies between their model and ours.

Second, we have tested whether the large abundance anomalies predicted by our models stem specifically from our simplified treatment of diffusion. For this purpose we have ported into our code the relevant parts from the code used by Hu et al. (priv. comm.)6 and run a M = 0.8 M⊙, Z = 10-4 diffusive model without radiative levitation (with the ZAMS chemical composition from Table 1). In this model diffusion reduces the helium and metal abundances in the envelope on much shorter timescales. The same M = 0.8 M⊙ model run with the MESA code (Paxton et al. 2011, 2015) yields similar results, which is reassuring given that the treatment of diffusion is based on the work of Burgers (1969) and Thoul et al. (1994) in both MESA and Hu et al. (2010). The abundances in the MESA model are depleted to a much greater degree: 5 dex for helium (compared to our 2.5 dex) and 5–6 dex (3–4 dex) for metals, even though the envelope masses are in a good agreement (the minimum envelope mass is 4.5 × 10-7 M⊙ in STARS and 6.0 × 10-7 M⊙ in MESA). The conclusion from all three of these tests is the same – if anything, the diffusion models presented here underestimate the amount of diffusion that we would get from a more rigorous treatment.

We have performed spatial resolution tests in three systems (denoted by Λ in Table A.1) by varying the default number of meshpoints (999) by a factor of two. All models give consistent results (within a couple of percent) in terms of the global properties, depth of thermohaline mixing, abundance anomalies after turn-off, and post-FDU abundances. The size of the convective envelope at minimum is consistent within ten percent. We thus conclude that the models are sufficiently resolved.

Our approach of interpolating the opacities and accelerations from tables computed during the run time necessitates the introduction of some numerical parameters. These parameters control mainly the amount by which some species has to change to warrant the computation of a new table. We have done extensive tests to make sure that our results do not depend on the choice of these parameters, i.e. the tables are computed often enough. As stated earlier, we set the temperature above which we use the old opacity tables that include conduction (Eldridge & Tout 2004) to log T = 7.3. We have since included the conductive opacities from Cassisi et al. (2007) in our code and made sure that use of OP opacities above log T = 7.3 would have virtually no effect on any of our results.

The size of the convective envelope throughout the evolution depends somewhat on the choice of the mixing-length parameter, with larger values resulting in more massive envelopes. Our value, αMLT = 2.0, is based on a calibration between the radius, effective temperature, and luminosity of a Z = 0.0142, M = 1 M⊙ diffusive model with OP opacities at an age of 4.56 Gyr and the Sun (our αMLT value is slightly smaller than the value of 2.025 presented by Stancliffe et al. 2016, because of the different opacities). Stars of masses, metallicities, and evolutionary stages different from the Sun should have other values of αMLT (e.g. Trampedach et al. 2014), but meaningful quantitative predictions are virtually impossible. In a 0.8 M⊙, Z = 10-4 model increasing or decreasing αMLT by 5% accordingly changes the envelope mass by about 50%, which, given that Menv< 10-6 M⊙, translates into substantial changes in the surface abundances (Fig. 2). Since theoretical models suggest that at lower metallicities one should use lower αMLT values (at least for − 0.6 < [Fe/H] < + 0.3; Bonaca et al. 2012), it is unlikely that we have overestimated the importance of diffusion by underestimating the value of αMLT.

5.1. Missing mixing processes

The strong abundance anomalies predicted by diffusive models are not observed in CEMP stars. This suggests that in real stars atomic diffusion is inhibited by some physical process that we have not included in our models. While we leave a more in-depth investigation of possible culprits to future work, we examine a simple test case here. We add a “turbulent” diffusion term to Dmix in Eq. (1)as proposed by Richer et al. (2000): (7)This type of parametrization extends the surface mixing region down to where the local temperature is somewhat larger than T0. Richard et al. (2005) find that observations of lithium abundances in population II stars require T0 ≈ 106 K. Similar values have been found to reproduce the small systematic abundance differences between turn-off and giant stars in old globular clusters (Korn et al. 2006; Nordlander et al. 2012; Gruyters et al. 2014). Nevertheless, the ad-hoc nature of this prescription should be kept in mind. With this prescription we can primarily constrain the depth to which some form of mixing must occur to reconcile the models with observations, but not the physical processes responsible for this mixing.

We test the effect of turbulent diffusion on the evolution of a 0.75 M⊙ star accreting 0.1 M⊙ of material from a 1 M⊙ primary. While in the absence of turbulence the resulting 0.85 M⊙ model shows extremely large abundance anomalies ([Fe/H] > −0.5 with levitation and [Fe/H] < −9.2 without levitation), turbulence with D0 = 400 (as used by Richard et al. 2005) and log T0 = 6.0 completely negates them (Fig. 7). Even much smaller turbulent diffusion coefficients (e.g. D0 = 1) suffice to erase the anomalies. Indeed, the key parameter here is T0 – as long as the mixing region remains large enough (log T0 ≳ 5.5 or Menv ≳ 10-4 M⊙), atomic diffusion is strongly suppressed. In terms of global properties, models with more pervasive turbulence are hotter and therefore more closely resemble models without diffusion. This can be seen from comparing the turbulent models with the model with thermohaline mixing only (solid grey line) in Fig. 7c7.

thumbnail Fig. 7

Effect of turbulence on the evolution of a) carbon and b) iron mass fractions, and c) the HRD of a 0.75 M⊙ star accreting 0.1 M⊙ of material from a 1 M⊙ primary. Large abundance variations are expected after accretion in absence of turbulence (black and red dashes). Inclusion of turbulence as in Eq. (7)erases all abundance signatures of atomic diffusion on the post-mass-transfer main sequence (dark blue), even for small turbulent diffusion coefficients (orange). Only when the temperature parameter is reduced to log T0 ≲ 5.5 do the abundance variations start to reappear (light blue and magenta). The thermohaline-mixing-only model (solid grey) is hotter than the models with turbulence, because in the latter diffusion still modifies the layers with T ≳ T0, and thermohaline mixing is not as deep. For clarity, the HRD only shows the post-mass-transfer part of the evolutionary tracks.
Open with DEXTER

Figure 7 also shows that, while the turbulent diffusion prescription of Richer et al. (2000) can inhibit diffusion in the outer layers of a star, it has almost no influence on thermohaline mixing. This is not to say, however, that some form of turbulence (e.g. rotationally driven horizontal turbulence; Denissenkov & Pinsonneault 2008) could not inhibit thermohaline mixing as well. Rather, investigating this requires treating the different processes together instead of considering them as independent and simply adding the individual diffusion coefficients (Maeder et al. 2013). Such a treatment is beyond the scope of this paper.

5.2. Mass loss

So far we have ignored mass loss. Simple estimates imply that it may be too important to neglect. With a mass-loss rate comparable to the current solar value, 2–3 × 10-14 M⊙yr-1 (Wang 1998), a star will lose more than 10-4 M⊙ over the roughly 1010 years it spends on the main sequence. This is a very large amount compared to the envelope masses of our models (Fig. 2) and could greatly interfere with atomic diffusion.

In the absence of mass loss the convective envelope of a star moves outwards in mass until the beginning of FDU. But when mass loss erodes the surface, this outward movement is halted and eventually reversed while the star is still on the main sequence. As the envelope now moves inwards, the surface abundances reflect the composition of the progressively deeper layers that get exposed. Qualitatively, if the mass-loss rate is sufficiently high, the removal of the outer layers is so fast that diffusion has not had enough time to modify the newly exposed layers, and only small abundance anomalies can develop (Swenson 1995). On the other hand, mass-loss rates below some limit must be negligible and have essentially no effect on the surface abundances.

We now estimate what mass-loss rates are necessary to prevent the development of abundance anomalies and what mass-loss rates are negligible. For simplicity, we use Reimers (1975) mass-loss formula with different factors η: (8)Here we only consider metal-poor 0.8 M⊙ and 0.85 M⊙ models without accretion. We find that if η ≳ 0.1, the effects of atomic diffusion are almost entirely erased in both models (Fig. 8). This translates to mass-loss rates of a few times 10-14–10-13 M⊙yr-1 throughout the main sequence. In contrast, when the mass-loss rate falls below about 10-16 M⊙yr-1, the abundance evolution proceeds as in models without mass loss. Intermediate mass-loss rates result in less extreme but non-negligible abundance variations. We note that the lost material is assumed to have the same composition as the surface at that time. Depending on the mass-loss rate and mechanism, some elements may be lost more readily than others, leading to more complicated abundance variations (e.g. Michaud et al. 1987; Vick et al. 2010).

thumbnail Fig. 8

Envelope mass and abundance evolution in metal-poor 0.8 M⊙ (left panels) and 0.85 M⊙ (right panels) models with different Reimers-type mass-loss rates. The values in the parentheses are the mass-loss rates at log L ≈ 0; the average mass-loss rate during the main sequence is about 50% higher as a result of the mass, radius, and luminosity scaling in Reimers’ law.
Open with DEXTER

Is it reasonable to expect somewhat super-solar mass-loss rates from CEMP-s stars on the main sequence? That is difficult to say. The form of mass loss in these stars is presumably the same as in normal low-mass main sequence stars – as magnetized winds originating in a corona that is heated by turbulent dissipation of Alfvén waves (e.g. Suzuki 2007; Cranmer & Saar 2011). While mass-loss rates in these stars are too small to be directly observable, indirect measurements based on the interaction of the wind with the interstellar medium yield values within about an order of magnitude of the solar mass-loss rate (Wood et al. 2002, 2005). Various theoretical models also normally predict mass-loss rates in this range (e.g. Holzwarth & Jardine 2007; Cranmer & Saar 2011; Johnstone et al. 2015). While most of these models concern stars with near solar metallicity, to first order the mass-loss rates are not expected to depend on metallicity. We note that the mass-loss rate in our test models increases over time because of the LR/M scaling. This is because there is no rotational dependence in the Reimers’ mass-loss law, which is reasonable given that it was derived from observations of red giants (i.e. slow rotators). But CEMP-s stars could have high rotation rates after the mass transfer phase, if the transferred material carries with it some angular momentum. Over most of the post-mass-transfer main sequence evolution their mass-loss rate would then be decreasing as the wind carried away the excess angular momentum (e.g. for solar-type stars Johnstone et al. 2015, give Ṁ ~ Ω1.3 ~ t-0.75), which should result in higher average mass-loss rates. This scenario may even have further complications aside from mass loss, because rapid rotation could directly lead to enhanced chemical mixing in the star.

In any case, the mass-loss rates of CEMP-s dwarfs are likely at least a few times 10-15 M⊙ yr-1 throughout their evolution. Such mass-loss rates should at least moderate the effects of atomic diffusion and help explain why turn-off stars with extreme abundance anomalies are not observed. Additionally, or alternatively, some form of turbulence might play a role. As discussed by Vick et al. (2010), two models, one with mass loss and one with turbulent diffusion, that predict the same surface abundances do not necessarily have the same internal abundance profiles. In a model with turbulence the abundance profiles are flat down to some depth (e.g. determined by T0 in Eq. (7)). In a model with mass loss no mixing is enforced outside the convective region, so the abundance profiles will not be flat unless the mass-loss rate is large (Ṁ ≳ 10-13 M⊙ yr-1). Asteroseismic measurements sensitive to the internal structure of a star might, in principle, be able to distinguish between the two types of models (e.g. van Saders & Pinsonneault 2012). However, in practice the difference in the internal structure might be too small at these low metallicities for such measurements to be possible.

6. Conclusions

In this paper we present stellar evolution models of s-process-rich carbon-enhanced metal-poor (CEMP-s) stars under the assumption that they form when a low-mass metal-poor star accretes material from an AGB companion. Motivated by results from binary population synthesis calculations of Abate et al. (2015), our models cover current CEMP-s star masses between 0.8 and 0.95 M⊙, deriving from initial secondary masses between 0.6 and 0.8 M⊙, and initial primary masses between 0.9 and 1.5 M⊙. Our main focus is the post-mass-transfer evolution of the surface abundances of carbon and iron driven by thermohaline mixing and atomic diffusion, including radiative levitation.

Our simulations with atomic diffusion indicate that CEMP-s stars should show large surface abundance variations on the main sequence, particularly as they approach the turn-off. This is because they have very shallow convective envelopes (Menv ≲ 10-4 M⊙ throughout most of the evolution and perhaps as little as 10-8 M⊙ near the turn-off) and, therefore, short diffusion timescales. In stars whose envelope masses fall below about 10-5 M⊙ (which happens in most of our models, including nearly all those with M> 0.8 M⊙) the abundances should vary by a factor of about ten. This factor rapidly increases with decreasing envelope mass, resulting in unrealistic abundances (Fig. 2). But even though our treatment of diffusion is not as detailed as in some other works, we do not find evidence that the surface abundance variations predicted by our models are exaggerated.

Radiative levitation has only a minor influence on carbon but a large one on iron. Whereas in diffusive models without levitation the metallicity ([Fe/H]) of the star decreases until first dredge-up, in models with levitation the metallicity can increase as the star evolves along the main sequence. Consequently, models with levitation predict reduced carbon enhancements ([C/Fe]) around the turn-off. This implies a systematic difference between the [C/Fe] values of stars near the turn-off and those at the beginning of FDU. Unfortunately, whether there is any such difference is difficult to establish even from the largest homogeneous sample of observational data (from SDSS; Lee et al. 2013). Any such difference, however, would clearly be smaller than predicted by most of our models (Fig. 5). And while some of our 0.8 M⊙ models do predict only a small variation in [C/Fe], at ages typical of metal-poor halo stars most of them are still relatively unevolved and should be visible as carbon-rich low-luminosity objects. Very few such objects have been observed.

Although they too would predict many low-luminosity carbon-rich stars, models without atomic diffusion are generally much more successful at covering the range of observations (Fig. 6). We thus conclude that atomic diffusion cannot be acting alone near the surface convection zone of real CEMP-s stars and needs to be largely counteracted by some other physical process(es). For example, a turbulent diffusive process like proposed by Richard et al. (2002b) can suppress surface abundance variations almost entirely by extending the mixing region to depths where the temperature exceeds about 106 K (about 10-4 M⊙ from the surface; Fig. 7). Additionally, at least the most extreme abundance variations (corresponding to stars with the smallest envelopes) should also be moderated by mass loss. In fact, a mass-loss rate of a few times the current solar value sustained throughout the evolution could on its own prevent substantial abundance anomalies from developing (Fig. 8).

While this work has primarily dealt with carbon and iron, given how divergent their abundance evolution is expected to be, these conclusions should extend to other elements, including those produced by neutron capture. The common assumption that the material coming from the AGB companion has simply been diluted by some factor after accretion onto the CEMP-s star is likely not too far from the truth.

0

The relative abundance of element A with respect to element B is where C is the number or mass fraction.

1

Slightly different distinctions between CEMP stars enriched in s-process elements have also been proposed (Jonsell et al. 2006; Masseron et al. 2010; Allen et al. 2012).

2

Hansen et al. (2016b) have argued that four stars in their observed sample of 22 CEMP-s stars appear to be single. Regardless, the mass transfer scenario should apply to most CEMP-s stars.

3

Other symbols in Eq. (1)have their usual meaning, namely: t is time; ρ is the mass density; r and m are the radial and mass coordinate, respectively. The diffusion velocity wi is defined in Eq. (2).

4

In fact, the convective envelope has already slightly grown in mass. For a short time the diffusion timescale is shorter than the evolutionary timescale.

5

Two examples are the CEMP-no (or -r) star SDSS0036-10 ([Fe/H] = −2.4, log g = 4.5, [C/Fe] = 2.3, [Ba/Fe] = 0.3) and the CEMP-s (or -r/s) star SDSS2047+10 ([Fe/H] = −2.1, log g = 4.5, [C/Fe] = 2.0, [Ba/Fe] = 1.5) from Aoki et al. (2008). Behara et al. (2010) present the CEMP-r/s star SDSSJ0912+0216 ([Fe/H] = −2.5, log g = 4.5, [C/Fe] ≈ 1.5, [Ba/Fe] = 1.5, [Eu/Fe] = 1.2), which may, however, be more evolved with log g ≈ 4.0 (Aoki et al. 2013).

6

In their code the full set of Burgers flow equations (Burgers 1969) is solved, and their treatment of diffusion is thus valid for arbitrary compositions.

7

The model with D0 = 400 and log T0 = 6.0 is quite different from the basic model with no diffusion or thermohaline mixing, although their tracks almost coincide. In the basic model the accreted material remains on the surface of the star. In the turbulent model the material is diluted, but not as much as in a model without diffusion because of the stabilizing μ-gradient in layers where T ≳ T0.

Acknowledgments

We thank the anonymous referee for comments that have helped improve the clarity of the paper. We thank Carlo Abate for help with the model grid selection, constructive comments on this manuscript, and many useful discussions. We also thank Young Sun Lee and Timothy Beers for sharing the SDSS CEMP data, Haili Hu for sharing her code, and Evert Glebbeek and Olivier Richard for useful discussions. R.J.S. is the recipient of a Sofja Kovalevskaja Award from the Alexander von Humboldt Foundation.

References

Appendix A: Additional tables

Table A.1

Results from simulations including atomic diffusion.

Table A.2

Results from simulations without atomic diffusion.

All Tables

Table 1

Chemical composition of the secondaries on the zero-age main sequence (ZAMS; abundance distribution from Asplund et al. 2009, scaled to Z = 10-4) and of the ejecta from the AGB models of Lugaro et al. (2012).

In the text
Table A.1

Results from simulations including atomic diffusion.

In the text
Table A.2

Results from simulations without atomic diffusion.

In the text

All Figures

thumbnail Fig. 1

Evolution and abundances of a M2,i = 0.75 M⊙ secondary accreting ΔM = 0.05 M⊙ of material from a M1 = 1.25 M⊙ primary. The labels highlight specific parts of the evolution discussed in the text. The model sequences with and without radiative levitation overlap at this scale of the HRD. Interior abundance profiles are shown at a few of the stages indicated in the upper panels: before mass transfer (“1”, solid); before thermohaline mixing (“4”, long-dashed); after thermohaline mixing (“5”, dotted); post-mass-transfer main sequence (“6”, dot-dashed); minimum of convective envelope mass (“7”, dot-dot-dashed); during first dredge-up (“9”, short-dashed). The vertical lines in the lower panels indicate the position of the base of the convective envelope at the respective time. The interior profiles of carbon with and without levitation nearly coincide, and only the case with levitation is shown.
Open with DEXTER
In the text
thumbnail Fig. 2

a) [Fe/H], b) [C/H], and c) [C/Fe] in each of the CEMP-s models at the point where the mass of the convective envelope is smallest, i.e. just before first dredge-up. Models with and without radiative levitation are plotted with black and grey symbols, respectively. All values for envelope masses below 10-7 M⊙ denote upper or lower limits. The [Fe/H] plot also shows the metallicity evolution in a 0.8 M⊙ model sequence with no accretion (solid lines). Similarly, the [C/Fe] plot shows the evolution in a model corresponding to M1 = 1.25 M⊙, M2,i = 0.8 M⊙, and ΔM = 0.05 M⊙. Both sequences without levitation stop during the first dredge-up. The criterion for being classified as a carbon-enhanced metal-poor star ([C/Fe] ≥ 1.0) is from Beers & Christlieb (2005).
Open with DEXTER
In the text
thumbnail Fig. 3

Thermohaline mixing efficiency (fraction of the star that is mixed) as a function of the ratio between accreted mass and final mass for different primaries. Black and grey symbols correspond to models with and without diffusion, respectively. Each symbol represents a unique combination of M1, M2,f, and ΔM.
Open with DEXTER
In the text
thumbnail Fig. 4

Schematic illustration of the structure of the star a) before and b) after thermohaline mixing. The shaded area around mthm and below indicates a region in which the mean molecular weight has been raised as a result of nuclear processing.
Open with DEXTER
In the text
thumbnail Fig. 5

Evolution of [Fe/H] (upper panels), [C/H] (middle panels) and [C/Fe] (lower panels) in CEMP-s models of 0.8 M⊙ (left panels) and 0.85 M⊙ (right panels) with diffusion. Thick lines are models with radiative levitation, whereas thin lines are those without. Empty circles mark an age of 10 Gyr, and filled circles mark an age of 13.8 Gyr, i.e. the part of the track between the circles covers the expected age range of CEMP-s stars. The small, grey circles show the metal-poor stars observed in the Sloan Digital Sky Survey in the metallicity range − 2.5 ≤ [Fe/H] ≤ −2.0 (Lee et al. 2013).
Open with DEXTER
In the text
thumbnail Fig. 6

Similar to Figs. 5e and f, but here the lines correspond to CEMP-s models of 0.9 M⊙ without diffusion. Thick lines are models with thermohaline mixing, whereas thin lines are those without.

Open with DEXTER
In the text
thumbnail Fig. 7

Effect of turbulence on the evolution of a) carbon and b) iron mass fractions, and c) the HRD of a 0.75 M⊙ star accreting 0.1 M⊙ of material from a 1 M⊙ primary. Large abundance variations are expected after accretion in absence of turbulence (black and red dashes). Inclusion of turbulence as in Eq. (7)erases all abundance signatures of atomic diffusion on the post-mass-transfer main sequence (dark blue), even for small turbulent diffusion coefficients (orange). Only when the temperature parameter is reduced to log T0 ≲ 5.5 do the abundance variations start to reappear (light blue and magenta). The thermohaline-mixing-only model (solid grey) is hotter than the models with turbulence, because in the latter diffusion still modifies the layers with T ≳ T0, and thermohaline mixing is not as deep. For clarity, the HRD only shows the post-mass-transfer part of the evolutionary tracks.
Open with DEXTER
In the text
thumbnail Fig. 8

Envelope mass and abundance evolution in metal-poor 0.8 M⊙ (left panels) and 0.85 M⊙ (right panels) models with different Reimers-type mass-loss rates. The values in the parentheses are the mass-loss rates at log L ≈ 0; the average mass-loss rate during the main sequence is about 50% higher as a result of the mass, radius, and luminosity scaling in Reimers’ law.
Open with DEXTER
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.