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.     

Cosmological origin of the lowest metallicity halo stars

We explore the predictions of the standard hierarchical clustering scenario of galaxy formation, regarding the numbers and metallicities of PopIII stars that are likely to be found within our Galaxy today. By PopIII we refer to stars formed at large redshift ( z >4), with low metallicities ([ Z /Z_⊙]<−2.5) and in small systems (total mass ≲ 2×10⁸ M_⊙) that are extremely sensitive to stellar feedback, and which through a prescribed merging history end up becoming part of the Milky Way today. An analytic, extended Press–Schechter formalism is used to obtain the mass functions of haloes which will host PopIII stars at a given redshift, and which will end up in Milky Way sized systems today. Each of these is modelled as a mini-galaxy, with a detailed treatment of the dark halo structure, angular momentum distribution, final gas temperature and disc instabilities, all of which determine the fraction of the baryons that are subject to star formation. The use of new primordial metallicity stellar evolutionary models allows us to trace the history of the stars formed, and give accurate estimates of their expected numbers today and their location in L /L_⊙ versus T /K Hertzsprung–Russell (HR) diagrams. A first comparison with observational data suggests that the initial mass function (IMF) of the first stars was increasingly high-mass weighted towards high redshifts, levelling off at z ≳9 at a characteristic stellar mass scale m _s=10–15 M_⊙.     

Non-adiabatic linear pulsation models for low-mass helium stars

Non-adiabatic linear pulsation models have been calculated for low-mass stars with effective temperatures between 16 000 and 35 000 K, and with surface gravities in the range 3⊙, X =0.00, Z =0.02. It is shown that the Z -bump instability persists to low masses ( M ∼0.4 M_⊙) but is suppressed either by a reduction in metallicity Z or by a selective enhancement of the carbon abundance. An unexpected result is the discovery that Z -bump instability persists at hydrogen abundances X >0.3, although the position of the red edge is sensitive to X . We have found that non-radial pulsations are also excited in the same instability region as radial pulsations.
The implications of these results for individual low-mass helium stars are discussed. It is concluded that Z -bump driven pulsations (radial and/or non-radial) may be excited in some helium-rich subdwarf B stars, representing a possible major extension to the class of variable stars represented by the prototype V652 Her.     

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.     

Predictions of the extent of self-enrichment in oxygen of giant metal-poor H ii regions

In general, H  ii regions do not show clear signs of self-enrichment in products from massive stars  ( M ≥ 8 M_⊙)  . In order to explore why I modelled the contamination with Wolf–Rayet star ejecta of metal-poor  ( Z = 0.001)  H  ii regions, ionized either by a  10⁶ M_⊙  cluster of coeval stars (cluster 1) or by a cluster resulting from continuous star formation at a rate of  1 M_⊙ yr⁻¹  (cluster 2). The clusters have   Z = 0.001  and a Salpeter initial mass function from 0.1 to  120 M_⊙  . Independent one-dimensional constant density simulations of the emission-line spectra of unenriched H  ii regions were computed at the discrete ages 1, 2, 3, 4 and 5 Myr, with the photoionization code cloudy , using as input, radiative and mechanical stellar feedbacks predicted by the evolutionary synthesis code starburst99 . Each H  ii region was placed at the outer radius of the adiabatically expanding superbubble of Mac Low & McCray. For models with thermal and ionization balance time-scales of less than 1 Myr, and with oxygen emission-line ratios in agreement with observations, the volume of the superbubble and the H  ii region was uniformly and instantaneously polluted with stellar ejecta predicted by starburst99 . I obtained a maximum oxygen abundance enhancement of 0.025 dex, with cluster 1, at 4 Myr. It would be unobservable.     

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 depletion of carbon by extra mixing in metal-poor giants

There is an apparent dichotomy between the metal-poor  ([Fe/H]≤−2)  yet carbon-normal giants and their carbon-rich counterparts. The former undergo significant depletion of carbon on the red giant branch after they have undergone first dredge-up, whereas the latter do not appear to experience significant depletion. We investigate this in the context that the extra mixing occurs via the thermohaline instability that arises due to the burning of  ³He  . We present the evolution of [C/Fe], [N/Fe] and  ¹²C/¹³C  for three models: a carbon-normal metal-poor star, and two stars that have accreted material from a  1.5 M_⊙  AGB companion, one having received  0.01 M_⊙  of material and the other having received  0.1 M_⊙  . We find the behaviour of the carbon-normal metal-poor stars is well reproduced by this mechanism. In addition, our models also show that the efficiency of carbon-depletion is significantly reduced in carbon-rich stars. This extra-mixing mechanism is able to reproduce the observed properties of both carbon-normal and carbon-rich stars.     

Binary formation with different metallicities: dependence on initial conditions

The fragmentation process in collapsing clouds with various metallicities is studied using three-dimensional nested-grid hydrodynamics. Initial clouds are specified by three parameters: cloud metallicity, initial rotation energy and initial cloud shape. For different combinations of these parameters, we calculate 480 models in total and study cloud evolution, fragmentation conditions, orbital separation and binary frequency. For the cloud to fragment during collapse, the initial angular momentum must be higher than a threshold value, which decreases with decreasing metallicity. Although the exact fragmentation conditions depend also on the initial cloud shape, this dependence is only modest. Our results indicate a higher binary frequency in lower metallicity gas. In particular, with the same median rotation parameter as in the solar neighbourhood, a majority of stars are born as members of binary/multiple systems for  <10⁻⁴ Z_⊙  . With initial mass  <0.1 M_⊙  , if fragments are ejected in embryo from the host clouds by multibody interaction, they evolve to substellar-mass objects. This provides a formation channel for low-mass stars in zero- or low-metallicity environments.     

On the first generation of stars

We argue that the first stars may have spanned the conventional mass range rather than be identified with the very massive objects  (∼100–10³ M_⊙)  favoured by numerical simulations. Specifically, we find that magnetic field generation processes acting in the first protostellar systems suffice to produce fields that exceed the threshold for magneto-rotational instability (MRI) to operate, and thereby allow the MRI dynamo to generate equipartition-amplitude magnetic fields on protostellar mass scales below  ∼50 M_⊙  . Such fields allow primordial star formation to occur at essentially any metallicity by regulating angular momentum transfer, fragmentation, accretion and feedback in much the same way as occurs in conventional molecular clouds.     

Carbon abundances and 12C/13C from globular cluster giants

The behaviour of the  Δν= 2 CO  bands around 2.3 μm was examined by comparing observed and synthetic spectra in stars in globular clusters of different metallicity. Changes in the ¹²C/¹³C isotopic ratio and the carbon abundances were investigated in stars from 3500–4900 K in the galactic globular clusters M71, M5, M3 and M13, covering the metallicity range from −0.7 to −1.6. We found relatively low carbon abundances that are not affected by the value of oxygen abundance. For most giants, the ¹²C/¹³C ratios determined are consistent with the equilibrium value for the CN cycle. This suggests complete mixing on the ascent of the red giant branch, in contrast to the substantially higher values predicted across this range of parameters by the current generation of models. We found some evidence for a larger dispersion of ¹²C/¹³C in giants of M71 of metallicity  [μ]=[M/H]=−0.7  in comparison with the giants of M3, M5 and M13, which are more metal deficient. Finally, we show evidence for lower ¹²C/¹³C in giants of globular clusters with lower metallicities, as predicted by theory.     

IZ photometry of L dwarfs and the implications for brown dwarf surveys

The I − Z colour has recently been shown to be a good temperature indicator for M dwarfs. We present the first IZ photometry of a small sample of L dwarfs ranging in spectral type from L0.5V to L6.0V. We find that the I − Z colour is not a good temperature indicator for objects between L1V and L5V, such objects having colours that overlap with mid M dwarfs. We attribute this to the reduction in the strength of the TiO and VO bands in the L dwarfs, which are the dominant opacity source in the I band for late M dwarfs. Beyond L5V, I − Z appears to be a reasonable indicator. This has important implications for the planning of optical surveys for cool objects in clusters and the field. For example, I − Z will cease to be a good method of identifying brown dwarfs in the Pleiades below around 0.04 M_⊙, and at around 0.075 M_⊙ in the Hyades and Praesepe.     

The colours of the Sun

We compile a sample of Sun-like stars with accurate effective temperatures, metallicities and colours (from the ultraviolet to the near-infrared). A crucial improvement is that the effective temperature scale of the stars has recently been established as both accurate and precise through direct measurement of angular diameters obtained with stellar interferometers. We fit the colours as a function of effective temperature and metallicity, and derive colour estimates for the Sun in the Johnson–Cousins, Tycho, Strömgren, 2MASS and SDSS photometric systems. For  ( B − V )_⊙  , we favour the 'red' colour 0.64 versus the 'blue' colour 0.62 of other recent papers, but both values are consistent within the errors; we ascribe the difference to the selection of Sun-like stars versus interpolation of wider colour– T _eff–metallicity relations.     

The abundance of brown dwarfs

The amount of mass contained in low-mass objects is investigated anew. Instead of using a mass–luminosity relation to convert a luminosity function to a mass function, I predict the mass–luminosity relation from assumed mass functions and the luminosity functions of Jahreiss & Wielen and Gould, Bahcall & Flynn. Comparison of the resulting mass–luminosity relations with data for binary stars constrains the permissible mass functions. If the mass function is assumed to be a power law, the best-fitting slope lies either side of the critical slope, α =−2, below which the mass in low-mass objects is divergent, depending on the luminosity function adopted. If these power-law mass functions are truncated at 0.001 M_⊙, the contribution to the local density from stars lies between 0.013 and 0.10 M_⊙ pc⁻³ depending on the mass at which the mass function is normalized and the adopted value of α . Recent dynamical estimates of the local mass density rule out stellar mass densities above ∼0.05 M_⊙ pc⁻³. Hence, power laws steeper than α =−2 that extend down to 0.001 M_⊙ are allowed only if one adopts an implausible normalization of the mass function. If the mass function is generalized from a power law to a low-order polynomial in log( M ), the mass in stars with M <0.1 M_⊙ is either negligible or strongly divergent, depending on the order of the polynomial adopted.     

Dust formation in primordial Type II supernovae

We have investigated the formation of dust in the ejecta of Type II supernovae (SNe), mostly of primordial composition, to answer the question of where the first solid particles are formed in the Universe. However, we have also considered non-zero progenitor metallicity values up to Z = Z_⊙ . The calculations are based on standard nucleation theory, and the scheme has been tested for the first time on the well-studied case of SN1987A, yielding results that are in agreement with the available data. We find that: (i) the first dust grains are predominantly made of silicates, amorphous carbon (AC), magnetite and corundum; and (ii) the largest grains are the AC ones, with sizes around 300 Å, whereas the other grain types have smaller radii, around 10–20 Å . The grain size distribution depends somewhat on the thermodynamics of the ejecta expansion, and variations in the results by a factor ≈2 might occur within reasonable estimates of the relevant parameters. Also, and for the same reason, the grain size distribution is essentially unaffected by metallicity changes. The predictions on the amount of dust formed are very robust: for Z =0 , we find that SNe with masses in the range (12–35) M_⊙ produce about 0.08 M_⊙≲ M _d≲0.3 M_⊙ of dust per supernova. The above range increases by roughly three times as the metallicity is increased to solar values. We discuss the implications and the cosmological consequences of the results.     

The initial mass function of the rich young cluster NGC 1818 in the Large Magellanic Cloud

We use deep Hubble Space Telescope photometry of the rich, young (∼20- to 45-Myr old) star cluster NGC 1818 in the Large Magellanic Cloud to derive its stellar mass function (MF) down to  ∼0.15 M_⊙  . This represents the deepest robust MF thus far obtained for a stellar system in an extragalactic, low-metallicity  ([Fe/H]≃−0.4 dex)  environment. Combining our results with the published MF for masses above  1.0 M_⊙  , we obtain a complete present-day MF. This is a good representation of the cluster's initial MF (IMF), particularly at low masses, because our observations are centred on the cluster's uncrowded half-mass radius. Therefore, stellar and dynamical evolution of the cluster will not have affected the low-mass stars significantly. The NGC 1818 IMF is well described by both a lognormal and a broken power-law distribution with slopes of  Γ= 0.46 ± 0.10  and  Γ≃−1.35  (Salpeter-like) for masses in the range from 0.15 to  0.8 M_⊙  and greater than  0.8 M_⊙  , respectively. Within the uncertainties, the NGC 1818 IMF is fully consistent with both the Kroupa solar neighbourhood and the Chabrier lognormal mass distributions.     

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.     

Modelling the components of binaries in the Hyades: the dependence of the mixing-length parameter on stellar mass

We present our findings based on a detailed analysis of the binaries of the Hyades, in which the masses of the components are well known. We fit the models of the components of a binary system to observations so as to give the observed total V and B − V of that system and the observed slope of the main sequence in the corresponding parts. According to our findings, there is a very definite relationship between the mixing-length parameter and the stellar mass. The fitting formula for this relationship can be given as  α= 9.19( M /M_⊙− 0.74)^0.053− 6.65  , which is valid for stellar masses greater than  0.77 M_⊙  . While no strict information is gathered for the chemical composition of the cluster, as a result of degeneracy in the colour–magnitude diagram, by adopting   Z = 0.033  and using models for the components of 70 Tau and θ² Tau we find the hydrogen abundance to be   X = 0.676  and the age to be 670 Myr. If we assume that   Z = 0.024  , then   X = 0.718  and the age is 720 Myr. Our findings concerning the mixing-length parameter are valid for both sets of the solution. For both components of the active binary system V818 Tau, the differences between radii of the models with   Z = 0.024  and the observed radii are only about 4 per cent. More generally, the effective temperatures of the models of low-mass stars in the binary systems studied are in good agreement with those determined by spectroscopic methods.