The evolution of low-metallicity asymptotic giant branch stars and the formation of carbon-enhanced metal-poor stars

We investigate the behaviour of asymptotic giant branch (AGB) stars between metallicities   Z = 10⁻⁴  and 10⁻⁸. We determine which stars undergo an episode of flash-driven mixing, where protons are ingested into the intershell convection zone, as they enter the thermally pulsing AGB phase and which undergo third dredge-up. We find that flash-driven mixing does not occur above a metallicity of   Z = 10⁻⁵  for any mass of star and that stars above  2 M_⊙  do not experience this phenomenon at any metallicity. We find carbon ingestion (CI), the mixing of carbon into the tail of hydrogen-burning region, occurs in the mass range  2 M_⊙  to around  4 M_⊙  . We suggest that CI may be a weak version of the flash-driven mechanism. We also investigate the effects of convective overshooting on the behaviour of these objects. Our models struggle to explain the frequency of Carbon-Enhanced Metal-Poor (CEMP) stars that have both significant carbon and nitrogen enhancement. Carbon can be enhanced through flash-driven mixing, CI or just third dredge-up. Nitrogen can be enhanced through hot bottom burning and the occurrence of hot dredge-up also converts carbon into nitrogen. The C/N ratio may be a good indicator of the mass of the primary AGB stars.     

Self-enrichment by asymptotic giant branch stars in globular clusters: comparison between intermediate and high metallicities

We present theoretical evolutionary sequences of intermediate-mass stars  ( M = 3 − 6.5 M_⊙)  with metallicity   Z = 0.004  . Our goal is to test whether the self-enrichment scenario by massive asymptotic giant branch stars may work for the high-metallicity globular clusters, after previous works by the same group showed that the theoretical yields by this class of objects can reproduce the observed trends among the abundances of some elements, namely the O–Al and O–Na anticorrelations, at intermediate metallicities, i.e.  [Fe/H]=−1.3  . We find that the increase in the metallicity favours only a modest decrease of the luminosity and the temperature at the bottom of the envelope for the same core mass, and also the efficiency of the third dredge-up is scarcely altered. On the contrary, differences are found in the yields, due to the different impact that processes with the same efficiency have on the overall abundance of envelopes with different metallicities. We expect the same qualitative patterns as in the intermediate-metallicity case, but the slopes of some of the relationships among the abundances of some elements are different. We compare the sodium–oxygen anticorrelation for clusters of intermediate metallicity ( Z ≈ 10⁻³) and clusters of metallicity large as in these new models. Although the observational data are still too scarce, the models are consistent with the observed trends, provided that only stars of   M ≳ 5 M_⊙  contribute to self-enrichment.     

An explosive end to intermediate-mass zero-metallicity stars and early Universe nucleosynthesis

We use the Cambridge stellar evolution code stars to model the evolution of 5 and  7 M_⊙  zero-metallicity stars. With enhanced resolution at the hydrogen- and helium-burning shell in the asymptotic giant branch (AGB) phases, we are able to model the entire thermally pulsing AGB (TP-AGB) phase. The helium luminosities of the thermal pulses are significantly lower than in higher metallicity stars so there is no third dredge-up. The envelope is enriched in nitrogen by hot-bottom burning of carbon that was previously mixed in during second dredge-up. There is no s -process enrichment owing to the lack of third dredge-up. The thermal pulses grow weaker as the core mass increases and they eventually cease. From then on the star enters a quiescent burning phase which lasts until carbon ignites at the centre of the star when the CO core mass is  1.36 M_⊙  . With such a high degeneracy and a core mass so close to the Chandrasekhar mass, we expect these stars to explode as type 1.5 supernovae, very similar to type Ia supernovae but inside a hydrogen-rich envelope.     

Metal-rich carbon stars in the Sagittarius dwarf spheroidal galaxy

We present spectroscopic observations from the Spitzer Space Telescope of six carbon-rich asymptotic giant branch (AGB) stars in the Sagittarius dwarf spheroidal galaxy (Sgr dSph) and two foreground Galactic carbon stars. The band strengths of the observed C₂H₂ and SiC features are very similar to those observed in Galactic AGB stars. The metallicities are estimated from an empirical relation between the acetylene optical depth and the strength of the SiC feature. The metallicities are higher than those of the Large Magellanic Cloud, and close to Galactic values. While the high metallicity could imply an age of around 1 Gyr, for the dusty AGB stars, the pulsation periods suggest ages in excess of 2 or 3 Gyr. We fit the spectra of the observed stars using the dusty radiative transfer model and determine their dust mass-loss rates to be in the range  1.0–3.3 × 10⁻⁸ M_⊙ yr⁻¹  . The two Galactic foreground carbon-rich AGB stars are located at the far side of the solar circle, beyond the Galactic Centre. One of these two stars shows the strongest SiC feature in our present Local Group sample.     

The helium abundance and ΔY/ΔZ in lower main-sequence stars

We use nearby K dwarf stars to measure the helium-to-metal enrichment ratio  Δ Y /Δ Z   , a diagnostic of the chemical history of the solar neighbourhood. Our sample of K dwarfs has homogeneously determined effective temperatures, bolometric luminosities and metallicities, allowing us to fit each star to the appropriate stellar isochrone and determine its helium content indirectly. We use a newly computed set of Padova isochrones which cover a wide range of helium and metal content.
Our theoretical isochrones have been checked against a congruous set of main-sequence binaries with accurately measured masses, to discuss and validate their range of applicability. We find that the stellar masses deduced from the isochrones are usually in excellent agreement with empirical measurements. Good agreement is also found with empirical mass-luminosity relations.
Despite fitting the masses of the stars very well, we find that anomalously low helium content (lower than primordial helium) is required to fit the luminosities and temperatures of the metal-poor K dwarfs, while more conventional values of the helium content are derived for the stars around solar metallicity.
We have investigated the effect of diffusion in stellar models and the assumption of local thermodynamic equilibrium (LTE) in deriving metallicities. Neither of these is able to resolve the low-helium problem alone and only marginally if the cumulated effects are included, unless we assume a mixing-length which is strongly decreasing with metallicity. Further work in stellar models is urgently needed.
The helium-to-metal enrichment ratio is found to be  Δ Y /Δ Z = 2.1 ± 0.9  around and above solar metallicity, consistent with previous studies, whereas open problems still remain at the lowest metallicities. Finally, we determine the helium content for a set of planetary host stars.     

The origin of the lead-rich stars in the Galactic halo: investigation of model parameters for the s-process

Several stars at the low-metallicity extreme of the Galactic halo show large spreads of lead and associated 'heavy' s-process elements ([Pb/hs]). Theoretically, an s-process pattern should be obtained from an AGB star with a fixed metallicity and initial mass. For the third dredge-up and the s-process model, several important properties depend primarily on the core mass of AGB stars. Zijlstra reported that the initial-to-final mass relation steepens at low metallicity, due to low mass-loss efficiency. This might affect the model parameters of the AGB stars, e.g. the overlap factor and the neutron irradiation time, in particular at low metallicity. The calculated results do indeed show that the overlap factor and the neutron irradiation time are significantly small at low metallicities, especially for  3.0 M_⊙ AGB  stars. The scatter of [Pb/hs] found in low metallicities can therefore be explained naturally when varying the initial mass of the low-mass AGB stars.     

The stellar populations of starburst galaxies through near-infrared spectroscopy

We study the central (inner few hundred parsecs) stellar populations of four starburst galaxies (NGC 34, 1614, 3310 and 7714) in the near-infrared (NIR), from 0.8 to 2.4 μm, by fitting combinations of stellar population models of various ages and metallicities. The NIR spectra of these galaxies feature many absorption lines. For the first time, we fit simultaneously as much as 15 absorption features in the NIR. The observed spectra are best explained by stellar populations containing a sizable amount (20–56 per cent by mass) of ∼1-Gyr-old stellar population with thermally pulsing asymptotic giant branch stars. We found that the metallicity of the stars which dominates the light is solar. Metallicities substantially different from solar give a worse fit. Though the ages and metallicities we estimate using the NIR spectroscopy are in agreement with values from the literature based on the ultraviolet/optical, we find older ages and a larger age spread. This may be due to the fact that the optical is mostly sensitive to the last episode of star formation, while the NIR better maintains the record of previous stellar generations. Another interesting result is that the reddening estimated from the whole NIR spectrum is considerably lower than that based on emission lines. Finally, we find a good agreement of the free emission-line spectrum with photoionization models, using as input spectral energy distribution the synthetic composite template we derived as best fit.     

On rejuvenation in massive binary systems

We introduce a set of stellar models for massive stars whose evolution has been affected by mass transfer in a binary system, at a range of metallicities. As noted by other authors, the effect of such mass transfer is frequently more than just rejuvenation. We find that, whilst stars with convective cores which have accreted only H-rich matter rejuvenate as expected, those stars which have accreted He-rich matter (e.g. at the end stages of conservative mass transfer) evolve in a way that is qualitatively similar to rejuvenated stars of much higher metallicity. Thus, the effects of non-conservative evolution depend strongly on whether He-rich matter is amongst the portion accreted or ejected. This may lead to a significant divergence in binary evolution paths with only a small difference in initial assumptions. We compare our models to observed systems and find approximate formulae for the effect of mass accretion on the effective age and metallicity of the resulting star.     

The red giant branch tip and bump of the Leo II dwarf spheroidal galaxy

We present V and I photometry of a  9.4 × 9.4 arcmin²  field centred on the dwarf spheroidal galaxy Leo II. The tip of the red giant branch (TRGB) is identified at   I ^TRGB= 17.83 ± 0.03  and adopting  〈[M/H]〉=−1.53 ± 0.2  from the comparison of RGB stars with Galactic templates, we obtain a distance modulus  ( m − M )₀= 21.84 ± 0.13  , corresponding to a distance   D = 233 ± 15 kpc  . Two significant bumps have been detected in the luminosity function of the RGB. The fainter bump (B1, at   V = 21.76 ± 0.05  ) is the RGB bump of the dominant stellar population while the actual nature of the brightest one (B2, at   V = 21.35 ± 0.05  ) cannot be firmly assessed on the basis of the available data; it may be due to the asymptotic giant branch clump of the main population or it may be a secondary RGB bump. The luminosity of the main RGB bump (B1) suggests that the majority of RGB stars in Leo II belong to a population that is ≳4 Gyr younger than the classical Galactic globular clusters. The stars belonging to the He-burning red clump are shown to be significantly more centrally concentrated than RR Lyrae and blue horizontal branch stars, probing the existence of an age/metallicity radial gradient in this remote dwarf spheroidal.     

Massive stars

We describe the present state of massive star research seen from the viewpoint of stellar evolution, with special emphasis on close binaries. Statistics of massive close binaries are reasonably complete for the Solar neighbourhood. We defend the thesis that within our knowledge, many scientific results where the effects of binaries are not included, have an academic value, but may be far from reality. In chapter I, we summarize general observations of massive stars where we focus on the HR diagram, stellar wind mass loss rates, the stellar surface chemistry, rotation, circumstellar environments, supernovae. Close binaries can not be studied separately from single stars and vice versa. First, the evolution of single stars is discussed (chapter I). We refer to new calculations with updated stellar wind mass loss rate formalisms and conclusions are proposed resulting from a comparison with representative observations. Massive binaries are considered in chapter II. Basic processes are briefly described, i.e. the Roche lobe overflow and mass transfer, the common envelope process, the spiral-in process in binaries with extreme mass ratio, the effects of mass accretion and the merging process, the implications of the (asymmetric) supernova explosion of one of the components on the orbital parameters of the binary. Evolutionary computations of interacting close binaries are discussed and general conclusions are drawn. The enormous amount of observational data of massive binaries is summarized. We separately consider the non-evolved and evolved systems. The latter class includes the semi-detached and contact binaries, the WR binaries, the X-ray binaries, the runaways, the single and binary pulsars. A general comparison between theoretical evolution and observations is combined with a discussion of specially interesting binaries: the evolved binaries HD 163181, HD 12323, HD 14633, HD 193516, HD 25638, HD 209481, Per and Sgr; the WR+OB binary V444 Cyg; the high mass X-ray binaries Vela X-1, Wray 977, Cyg X-1; the low mass X-ray binaries Her X-1 and those with a black hole candidate; the runaway Pup, the WR+compact companion candidates Cyg X-3, HD 50896 and HD 197406. We finally propose an overall evolutionary model of massive close binaries as a function of primary mass, mass ratio and orbital period. Chapter III deals with massive star population synthesis with a realistic population of binaries. We discuss the massive close binary frequency, mass ratio and period distribution, the observations that allow to constrain possible asymmetries during the supernova explosion of a massive star. We focuss on the comparison between observed star numbers (as a function of metallicity) and theoretically predicted numbers of stellar populations in regions of continuous star formation and in starburst regions. Special attention is given to the O-type star/WR star/red supergiant star population, the pulsar and binary pulsar population, the supernova rates. Received 17 July 1998     

A near-infrared and optical photometric study of the Sculptor dwarf spheroidal galaxy: implications for the metallicity spread

The Sculptor dwarf spheroidal galaxy has a giant branch with a significant spread in colour, symptomatic of an intrinsic age–metallicity spread. We present here a detailed study of the Sculptor giant branch and horizontal branch (HB) morphology, combining new near-infrared photometry from the Cambridge Infrared Survey Instrument (CIRSI), with optical data from the European Southern Observatory Wide Field Imager. For a Sculptor-like old and generally metal-poor system, the position of red giant branch (RGB) and asymptotic giant branch (AGB) stars on the colour–magnitude diagram (CMD) is mainly metallicity dependent. The advantage of using optical–near-infrared colours is that the position of the RGB locus is much more sensitive to metallicity than with optical colours alone. In contrast the HB morphology is strongly dependent on both metallicity and age. Therefore a detailed study of both the RGB in optical–near-infrared colours and the HB can help break the age–metallicity degeneracy. Our measured photometric width of the Sculptor giant branch corresponds to a range in metallicity of 0.75 dex. We detect the RGB and AGB bumps in both the near-infrared and the optical luminosity functions, and derive from them a mean metallicity of  [M/H]=−1.3 ± 0.1  . From isochrone fitting we derive a mean metallicity of  [Fe/H]=−1.42  with a dispersion of 0.2 dex. These photometric estimators are for the first time consistent with individual metallicity measurements derived from spectroscopic observations. No spatial gradient is detected in the RGB morphology within a radius of 13 arcmin, twice the core radius. On the other hand, a significant gradient is observed in the HB morphology index, confirming the 'second parameter problem' present in this galaxy. These observations are consistent with an early extended period of star formation continuing in time for a few Gyr.     

Low-mass binary-induced outflows from asymptotic giant branch stars

A significant fraction of planetary nebulae (PNe) and protoplanetary nebulae (PPNe) exhibit aspherical, axisymmetric structures, many of which are highly collimated. The origin of these structures is not entirely understood, however, recent evidence suggests that many observed PNe harbour binary systems, which may play a role in their shaping. In an effort to understand how binaries may produce such asymmetries, we study the effect of low-mass  (<0.3 M_⊙)  companions (planets, brown dwarfs and low-mass main-sequence stars) embedded into the envelope of a  3.0-M_⊙  star during three epochs of its evolution [red giant branch, asymptotic giant branch (AGB), interpulse AGB]. We find that common envelope evolution can lead to three qualitatively different consequences: (i) direct ejection of envelope material resulting in a predominately equatorial outflow, (ii) spin-up of the envelope resulting in the possibility of powering an explosive dynamo-driven jet and (iii) tidal shredding of the companion into a disc which facilitates a disc-driven jet. We study how these features depend on the secondary's mass and discuss observational consequences.     

The accretion of brown dwarfs and planets by giant stars – II. Solar-mass stars on the red giant branch

This paper extends our previous study of planet/brown dwarf accretion by giant stars to solar-mass stars located on the red giant branch. The model assumes that the planet is dissipated at the bottom of the convective envelope of the giant star. The evolution of the giant is then followed in detail. We analyse the effects of different accretion rates and different initial conditions. The computations indicate that the accretion process is accompanied by a substantial expansion of the star, and, in the case of high accretion rates, hot bottom burning can be activated. The possible observational signatures that accompany the engulfing of a planet are also extensively investigated. They include the ejection of a shell and a subsequent phase of IR emission, an increase in the ⁷Li surface abundance and a potential stellar metallicity enrichment, spin-up of the star because of the deposition of orbital angular momentum, the possible generation of magnetic fields and the related X-ray activity caused by the development of shear at the base of the convective envelope, and the effects on the morphology of the horizontal branch in globular clusters. We propose that the IR excess and high Li abundance observed in 4–8 per cent of the G and K giants originate from the accretion of a giant planet, a brown dwarf or a very low-mass star.     

Neutron star retention and millisecond pulsar production in globular clusters

We investigate the conditions by which neutron star retention in globular clusters is favoured. We find that neutron stars formed in massive binaries are far more likely to be retained. Such binaries are likely to then evolve into contact before encountering other stars, possibly producing a single neutron star after a common envelope phase. A large fraction of the single neutron stars in globular clusters are then likely to exchange into binaries containing moderate-mass main-sequence stars, replacing the lower-mass components of the original systems. These binaries will become intermediate-mass X-ray binaries (IMXBs), once the moderate-mass star evolves off the main sequence, as mass is transferred on to the neutron star, possibly spinning it up in the process. Such systems may be responsible for the population of millisecond pulsars (MSPs) that has been observed in globular clusters. Additionally, the period of mass-transfer (and thus X-ray visibility) in the vast majority of such systems will have occurred 5–10 Gyr ago, thus explaining the observed relative paucity of X-ray binaries today, given the MSP population.     

Red giant stellar collisions in the Galactic Centre

We show that collisions with stellar-mass black holes can partially explain the absence of bright giant stars in the Galactic Centre, first noted by Genzel et al. We show that the missing objects are low-mass giants and asymptotic giant branch stars in the range  1–3 M_⊙  . Using detailed stellar evolution calculations, we find that to prevent these objects from evolving to become visible in the depleted K bands, we require that they suffer collisions on the red giant branch, and we calculate the fractional envelope mass losses required. Using a combination of smoothed particle hydrodynamic calculations, restricted three-body analysis and Monte Carlo simulations, we compute the expected collision rates between giants and black holes, and between giants and main-sequence stars in the Galactic Centre. We show that collisions can plausibly explain the missing giants in the  10.5 < K < 12  band. However, depleting the brighter  ( K < 10.5)  objects out to the required radius would require a large population of black hole impactors which would in turn deplete the  10.5 < K < 12  giants in a region much larger than is observed. We conclude that collisions with stellar-mass black holes cannot account for the depletion of the very brightest giants, and we use our results to place limits on the population of stellar-mass black holes in the Galactic Centre.     

Rotation and lithium in the late-type stars of NGC 2516

We present CCD photometry and high-resolution spectroscopy of low-mass stars in the open cluster NGC 2516, which has an age of about 150 Myr and may have a much lower metallicity than the Pleiades. 24 probable F to early K type, single cluster members have been identified from their photometry and radial velocities, along with three possible spectroscopic binaries. The projected equatorial velocities are measured and compared with younger and older clusters. Several fast rotating late G /early K stars are seen, but all hotter stars have v _e sin  i  < 20 km s⁻¹. The data are consistent with angular momentum loss models with spin-down time-scales that increase from tens of Myr for G stars to hundreds of Myr for K stars. The observed X-ray activity is consistent with the currently accepted rotation–activity paradigm. Lithium abundances are derived from the Li  i 6708-Å line. The pattern of Li depletion is indistinguishable from that in the Pleiades, including a spread in the K0 stars, where the most rapid rotators suffer the least Li depletion. The observations argue in favour of either a metallicity in the range −0.1 < [Fe/H]< 0.0 for NGC 2516, or a lower metallicity and extra Li depletion through non-standard mixing modes which occurs on time-scales of only ∼ 50 Myr. Neither our low signal-to-noise ratio spectroscopy nor our photometry can constrain [Fe/H] sufficiently to decide between these possibilities. A detailed spectroscopic chemical abundance analysis is urgently required.     

100-yr mass-loss modulations on the asymptotic giant branch

We analyse the differences in infrared circumstellar dust emission between oxygen-rich Mira and non-Mira stars, and find that they are statistically significant. In particular, we find that these stars segregate in the K–[12] versus [12]–[25] colour–colour diagram, and have distinct properties of the IRAS LRS spectra, including the peak position of the silicate emission feature. We show that the infrared emission from the majority of non-Mira stars cannot be explained within the context of standard steady-state outflow models.
The models can be altered to fit the data for non-Mira stars by postulating non-standard optical properties for silicate grains, or by assuming that the dust temperature at the inner envelope radius is significantly lower (300–400 K) than typical silicate grain condensation temperatures (800–1000 K) . We argue that the latter is more probable and provide detailed model fits to the IRAS LRS spectra for 342 stars. These fits imply that two-thirds of non-Mira stars and one-third of Mira stars do not have hot dust (>500 K) in their envelopes.
The absence of hot dust can be interpreted as a recent (∼100 yr) decrease in the mass-loss rate. The distribution of best-fitting model parameters agrees with this interpretation and strongly suggests that the mass loss resumes on similar time-scales. Such a possibility appears to be supported by a number of spatially resolved observations (e.g. recent Hubble Space Telescope images of the multiple shells in the Egg Nebula) and is consistent with new dynamical models for mass loss on the asymptotic giant branch.     

Merged binary white dwarf evolution: rapidly accreting carbon–oxygen white dwarfs and the progeny of extreme helium stars

We have examined the evolution of merged low-mass double white dwarfs that become luminous helium stars. We have approximated the merging process by the rapid accretion of matter, consisting mostly of helium, on to a carbon–oxygen (CO) white dwarf. After a certain mass is accumulated, a helium shell flash occurs, the radius and luminosity increase and the star becomes a yellow giant. Mass accretion is stopped artificially when the total mass reaches a pre-determined value. When the mass above the helium-burning shell becomes small enough, the star evolves blueward almost horizontally in the Hertzsprung–Russell diagram. The theoretical models for the merger of a 0.6-M_⊙ CO white dwarf with a 0.3-M_⊙ He white dwarf agree very well with the observed locations of extreme helium stars in the  log  T _eff–log  g   diagram, with their observed rates of blueward evolution, and with luminosities and masses obtained from their pulsations. Together with predicted merger rates for  CO+He  white dwarf pairs, the evolutionary time-scales are roughly consistent with the observed numbers of extreme helium stars. Predicted surface carbon and oxygen abundances can be consistent with the observed values if carbon and oxygen produced in the helium shell during a previous asymptotic giant branch phase are assumed to exist in the helium zone of the initial CO white dwarfs. These results establish the  CO+He  white dwarf merger as the best, if not only, viable model for the creation of extreme helium stars and, by association, the majority of R Coronae Borealis stars.     

Asymptotic giant branch stars in the Phoenix dwarf galaxy

JHK _s near-infrared photometry of stars in the Phoenix dwarf galaxy is presented and discussed. Combining these data with the optical photometry of Massey et al. allows a rather clean separation of field stars from Phoenix members. The discovery of a Mira variable ( P = 425 d), which is almost certainly a carbon star, leads to an estimate of the distance modulus of 23.10 ± 0.18 that is consistent with other estimates and indicates the existence of a significant population of age ∼2 Gyr. The two carbon stars of Da Costa have   M _bol=−3.8  and are consistent with belonging to a population of similar age; some other possible members of such a population are identified. A Da Costa non-carbon star is  Δ K _s∼ 0.3  mag brighter than these two carbon stars. It may be an asymptotic giant branch star of the dominant old population. The nature of other stars lying close to it in the   K _s, ( J − K _s)  diagram needs studying.