Light element abundances in carbon-enhanced metal-poor stars

We model the evolution of the abundances of light elements in carbon-enhanced metal-poor (CEMP) stars, under the assumption that such stars are formed by mass transfer in a binary system. We have modelled the accretion of material ejected by an asymptotic giant branch star on to the surface of a companion star. We then examine three different scenarios: one in which the material is mixed only by convective processes, one in which thermohaline mixing is present and a third in which both thermohaline mixing and gravitational settling are taken in to account. The results of these runs are compared to light element abundance measurements in CEMP stars (primarily CEMP- s stars, which are rich in s -processes elements and likely to have formed by mass transfer from an AGB star), focusing on the elements Li, F, Na and Mg. None of the elements is able to provide a conclusive picture of the extent of mixing of accreted material. We confirm that lithium can only be preserved if little mixing takes place. The bulk of the sodium observations suggest that accreted material is effectively mixed but there are also several highly Na and Mg-rich objects that can only be explained if the accreted material is unmixed. We suggest that the available sodium data may hint that extra mixing is taking place on the giant branch, though we caution that the data are sparse.     

Influence of the r-mode instability on hypercritically accreting neutron stars

We have investigated the influence of the r-mode instability on hypercritically accreting neutron stars in close binary systems during their common envelope phases, based on the scenario proposed by Brown et al. On the one hand, neutron stars are heated by the accreted matter at the stellar surface, but on the other hand they are also cooled down by the neutrino radiation. At the same time, the accreted matter transports its angular momentum and mass to the star. We have studied the evolution of the stellar mass, temperature and rotational frequency.
The gravitational-wave-driven instability of the r-mode oscillation strongly suppresses spinning up of the star, the final rotational frequency of which is well below the mass-shedding limit, in fact typically as low as 10 per cent of that of the mass-shedding state. On a very short time-scale the rotational frequency tends to approach a certain constant value and saturates there, as long as the amount of accreted mass does not exceed a certain limit to collapse to a black hole. This implies that a similar mechanism of gravitational radiation to that in the so-called 'Wagoner star' may work in this process. The star is spun up by accretion until the angular momentum loss by gravitational radiation balances the accretion torque. The time-integrated dimensionless strain of the radiated gravitational wave may be large enough to be detectable by gravitational wave detectors such as LIGO II.     

Chemical enrichment by Wolf–Rayet stars: non-solar metallicities

Following on from our recent Paper I, we present theoretical models of Wolf–Rayet (WR) stars for non-solar metallicities from   Z = 0.03  to 0.0001 by mass fraction with different mass-loss rate assumptions. We find that some single WR stars may still form even at the lowest metallicities, but whether this occurs or not depends critically on the upper cut-off point of the initial mass function used. As at solar metallicity, a population of binaries is required to fully reproduce WR star observations. For most scenarios, these binaries dominate the low-metallicity WR population but probably not the enrichment. We find comparable carbon enrichment from single WR stars to that from asymptotic giant branch stars at all metallicities for which data are available, but which of them is the dominant source of carbon depends strongly on the set of asymptotic giant branch yields adopted and the assumed initial mass function. We find an increase in carbon enrichment with increasing metallicity but a decrease in oxygen enrichment, as confirmed by observation.     

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.     

Dark stars at the Galactic Centre – the main sequence

In regions of very high dark matter density such as the Galactic Centre, the capture and annihilation of WIMP dark matter by stars has the potential to significantly alter their evolution. We describe the dark stellar evolution code D ark S tars , and present a series of detailed grids of WIMP-influenced stellar models for main-sequence stars. We describe the changes in stellar structure and main-sequence evolution which occur as a function of the rate of energy injection by WIMPs, for masses of  0.3–2.0 M_⊙  and metallicities   Z = 0.0003–0.02  . We show what rates of energy injection can be obtained using realistic orbital parameters for stars at the Galactic Centre, including detailed consideration of the velocity and density profiles of dark matter. Capture and annihilation rates are strongly boosted when stars follow elliptical rather than circular orbits. If there is a spike of dark matter induced by the supermassive black hole at the Galactic Centre, single solar mass stars following orbits with periods as long as 50 yr and eccentricities as low as 0.9 could be significantly affected. Binary systems with similar periods about the Galactic Centre could be affected on even less eccentric orbits. The most striking observational effect of this scenario would be the existence of a binary consisting of a low-mass protostar and a higher mass evolved star. The observation of low-mass stars and/or binaries on such orbits would either provide a detection of WIMP dark matter, or place stringent limits on the combination of the WIMP mass, spin-dependent nuclear-scattering cross-section, halo density and velocity distribution near the Galactic Centre. In some cases, the derived limits on the WIMP mass and spin-dependent nuclear-scattering cross-section would be of comparable sensitivity to current direct-detection experiments.     

After the Supernova: Runaway Stars and Massive X-ray Binary Populations with Metallicity

One of the main pathways by which massive runaways are thought to be produced is by the disruption of a binary system after the supernova (SN) of one of its components. Under such a scenario, the populations of runaway stars in different phases should reflect the input binary population and its evolution. Conversely, if the system stays together after the SN, a High Mass X-Ray Binary (HMXB) may result. We present simulations exploring the behaviour of such runaway and HMXB populations with metallicity, and compare them to observations. As many as two-thirds of massive runaway stars may be produced by supernovae in binaries. Decreasing metallicity lowers the fraction of O stars which are runaway, but increases the Wolf-Rayet runaway fractions and the number of potential HMXBs.     

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.     

The fragmentation of pre-enriched primordial objects

Recent theoretical investigations have suggested that the formation of the very first stars, forming out of metal-free gas, was fundamentally different from the present-day case. The question then arises which effect was responsible for this transition in the star formation properties. In this paper, we study the effect of metallicity on the evolution of the gas in a collapsing dark matter mini-halo. We model such a system as an isolated 3 σ peak of mass     that collapses at     , using smoothed particle hydrodynamics. The gas has a supposed level of pre-enrichment of either     or 10⁻³ Z_⊙. We assume that H₂ has been radiatively destroyed by the presence of a soft UV background. Metals therefore provide the only viable cooling at temperatures below 10⁴ K. We find that the evolution proceeds very differently for the two cases. The gas in the lower metallicity simulation fails to undergo continued collapse and fragmentation, whereas the gas in the higher metallicity case dissipatively settles into the centre of the dark matter halo. The central gas, characterized by densities     , and a temperature,     , that closely follows that of the cosmic microwave background, is gravitationally unstable and undergoes vigorous fragmentation. We discuss the physical reason for the existence of a critical metallicity,     , and its possible dependence on redshift. Compared with the pure H/He case, the fragmentation of the     gas leads to a larger relative number of low-mass clumps.     

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.     

A revision of the solar neighbourhood metallicity distribution

We present a revised metallicity distribution of dwarfs in the solar neighbourhood. This distribution is centred on solar metallicity. We show that previous metallicity distributions, selected on the basis of spectral type, are biased against stars with solar metallicity or higher. A selection of G-dwarf stars is inherently biased against metal-rich stars and is not representative of the solar neighbourhood metallicity distribution. Using a sample selected on colour, we obtain a distribution where approximately half the stars in the solar neighbourhood have metallicities higher than [Fe/H]=0 . The percentage of mid-metal-poor stars ([Fe/H]<−0.5) is approximately 4 per cent, in agreement with present estimates of the thick disc.
In order to have a metallicity distribution comparable to chemical evolution model predictions, we convert the star fraction to mass fraction, and show that another bias against metal-rich stars affects dwarf metallicity distributions, due to the colour (or spectral type) limits of the samples. Reconsidering the corrections resulting from the increasing thickness of the stellar disc with age, we show that the simple closed-box model with no instantaneous recycling approximation gives a reasonable fit to the observed distribution. Comparisons with the age–metallicity relation and abundance ratios suggest that the simple closed-box model may be a viable model of the chemical evolution of the Galaxy at solar radius.     

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.     

F stars, metallicity and the ages of red galaxies at z > 1

We explore whether the rest-frame near-ultraviolet spectral region, observable in high-redshift galaxies via optical spectroscopy, contains sufficient information to allow the degeneracy between age and metallicity to be lifted. We do this by first testing the ability of evolutionary synthesis models to reclaim the correct metallicity when fitted to the near-ultraviolet spectra of F stars of known (subsolar and supersolar) metallicity. F stars are of particular interest because the rest-frame near-ultraviolet spectra of the oldest known elliptical galaxies at   z > 1  appear to be dominated by F stars near to the main-sequence turn-off.
We find that, in the case of the F stars, where the Hubble Space Telescope ultraviolet spectra have a high signal-to-noise ratio, fitting models in which the metallicity is allowed to vary as a free parameter is rather successful at deriving the correct metallicity. As a result, the estimated turn-off ages of these stars yielded by model-fitting are well constrained. Encouraged by this we have fitted these same variable-metallicity models to the deep, optical spectra of the   z ≃ 1.5 mJy  radio galaxies 53W091 and 53W069 obtained with the Keck telescope. While the age and metallicity are not so easily constrained for these galaxies, we find that even when metallicity is allowed as a free parameter, the best estimates of their ages are still ≥3 Gyr, with ages younger than 2 Gyr now strongly excluded. Furthermore, we find that a search of the entire parameter space of metallicity and star formation history using MOPED leads to the same conclusion. Our results therefore continue to argue strongly against an Einstein–de Sitter universe, and favour a Λ-dominated universe in which star formation in at least these particular elliptical galaxies was completed somewhere in the redshift range   z = 3–5  .     

The circumstellar environment of Wolf–Rayet stars and gamma-ray burst afterglows

We study the evolution of the circumstellar medium of massive stars. We pay particular attention to Wolf-Rayet stars that are thought to be the progenitors of some long gamma-ray bursts (GRBs). We detail the mass-loss rates we use in our stellar evolution models and how we estimate the stellar wind speeds during different phases. With these details we simulate the interactions between the wind and the interstellar medium to predict the circumstellar environment around the stars at the time of core-collapse. We then investigate how the structure of the environment might affect the GRB afterglow. We find that when the afterglow jet encounters the free-wind/stalled-wind interface, rebrightening occurs and a bump is seen in the afterglow light curve. However, our predicted positions of this interface are too distant from the site of the GRB to reach while the afterglow remains observable. The values of the final wind density,   A _*  , from our stellar models are of the same order (≲1) as some of the values inferred from observed afterglow light curves. We do not reproduce the lowest   A _*  values below 0.5 inferred from afterglow observations. For these cases, we suggest that the progenitors could have been a WO-type Wolf–Rayet (WR) star or a very low-metallicity star. Finally, we turn our attention to the matter of stellar wind material producing absorption lines in the afterglow spectra. We discuss the observational signatures of two WR stellar types, WC and WO, in the afterglow light curve and spectra. We also indicate how it may be possible to constrain the initial mass and metallicity of a GRB progenitor by using the inferred wind density and wind velocity.     

Evolution of a 3-M⊙ star from the main sequence to the ZZ Ceti stage: the role played by element diffusion

The purpose of this paper is to present new full evolutionary calculations for DA white dwarf stars with the major aim of providing a physically sound reference frame for exploring the pulsation properties of the resulting models in future communications. Here, white dwarf evolution is followed in a self-consistent way with the predictions of time-dependent element diffusion and nuclear burning. In addition, full account is taken of the evolutionary stages prior to white dwarf formation. In particular, we follow the evolution of a 3-M_⊙ model from the zero-age main sequence (the adopted metallicity is   Z =0.02)  , all the way from the stages of hydrogen and helium burning in the core up to the thermally pulsing phase. After experiencing 11 thermal pulses, the model is forced to evolve towards its white dwarf configuration by invoking strong mass loss episodes. Further evolution is followed down to the domain of the ZZ Ceti stars on the white dwarf cooling branch.
Emphasis is placed on the evolution of the chemical abundance distribution caused by diffusion processes and the role played by hydrogen burning during the white dwarf evolution. We find that discontinuities in the abundance distribution at the start of the cooling branch are considerably smoothed out by diffusion processes by the time the ZZ Ceti domain is reached. Nuclear burning during the white dwarf stage does not represent a major source of energy, as expected for a progenitor star of initially high metallicity. We also find that thermal diffusion lessens even further the importance of nuclear burning.
Furthermore, the implications of our evolutionary models for the main quantities relevant for adiabatic pulsation analysis are discussed. Interestingly, the shape of the Ledoux term is markedly smoother compared with previous detailed studies of white dwarfs. This is translated into a different behaviour of the Brunt–Väisälä frequency.     

X-rays from the open cluster NGC 6633

NGC 6633 is a young, open cluster with a similar age to the Hyades and Praesepe, but probably a lower metallicity. We present the results of ROSAT High Resolution Imager observations of an optically selected catalogue of likely members of NGC 6633. 8 out of 51 NGC 6633 members have been detected, with main-sequence spectral types A to G, above a threshold X-ray luminosity of ≈6–12×10²⁸ erg s⁻¹. We find that NGC 6633 does not contain cool stars that are as X-ray luminous as the most active objects in the Hyades and that the median X-ray luminosity of F-G stars in NGC 6633 is less than that in the Hyades, but probably greater than in Praesepe. However, when X-ray activity is expressed as the X-ray to bolometric flux ratio we find that NGC 6633 and the Hyades are very similar and display similar peak levels of coronal activity. We attribute this discrepancy to a number of possible wide binary systems with higher X-ray (and bolometric) luminosities in the Hyades sample and either a low metallicity in NGC 6633, which makes its cool stars both X-ray and bolometrically less luminous at the same colour, or a distance to NGC 6633 that has been underestimated, which would decrease stellar X-ray luminosities without changing X-ray to bolometric flux ratios.     

Modelling the binary progenitor of Supernova 1993J

We have developed a detailed stellar evolution code capable of following the simultaneous evolution of both stars in a binary system, together with their orbital properties. To demonstrate the capabilities of the code, we investigate potential progenitors for the Type IIb Supernova 1993J, which is believed to have been an interacting binary system prior to its primary exploding. We use our detailed binary stellar evolution code to model this system to determine the possible range of primary and secondary masses that could have produced the observed characteristics of this system, with particular reference to the secondary. Using the luminosities and temperatures for both stars (as determined by Maund et al.) and the remaining mass of the hydrogen envelope of the primary at the time of explosion, we find that if mass transfer is 100 per cent efficient, the observations can be reproduced by a system consisting of a  15 M_⊙  primary and a  14 M_⊙  secondary in an orbit with an initial period of 2100 days. With a mass transfer efficiency of 50 per cent, a more massive system consisting of a  17 M_⊙  primary and a  16 M_⊙  secondary in an initial orbit of 2360 days is needed. We also investigate some of the uncertainties in the evolution, including the effects of tidal interaction, convective overshooting and thermohaline mixing.     

贫金属富碳天琴RR变星的形成

贫金属富碳恒星(Carbon-Enhanced Metal-Poor, CEMP)是研究宇宙早期恒星性质和化学演化的极佳样本,通常认为来自双星.目前发现的贫金属富碳星中有9颗天琴RR变星(RR Lyrae star, RRL),其中至少7颗未表现出任何双星特征.传统双星物质转移模型不足以充分解释贫金属富碳天琴RR变星(CEMP-RR Lyrae)单星的形成.之前研究表明氦白矮星和赫氏空隙星(HG)的并合模型可以解释部分富碳红巨星单星的碳增丰现象,因此贫金属富碳星单星也可能来自氦白矮星和赫氏空隙星的并合模型渠道.通过详细计算的氦白矮星和赫氏空隙星并合模型来检验这一演化渠道,结果表明:该并合模型在后续的演化过程中,其重力加速度、温度、表面碳丰度均能与观测符合较好.由此,氦白矮星和赫氏空隙星并合模型极有可能是贫金属富碳天琴RR变星的形成渠道之一.     

Metallicity effects on the chromospheric activity–age relation for late-type dwarfs

We show that there is a relationship between the age excess, defined as the difference between the stellar isochrone and chromospheric ages, and the metallicity as measured by the index [Fe/H] for late-type dwarfs. The chromospheric age tends to be lower than the isochrone age for metal-poor stars, and the opposite occurs for metal-rich objects. We suggest that this could be an effect of neglecting the metallicity dependence of the calibrated chromospheric emission–age relation. We propose a correction to account for this dependence. We also investigate the metallicity distributions of these stars, and show that there are distinct trends according to the chromospheric activity level. Inactive stars have a metallicity distribution which resembles the metallicity distribution of solar neighbourhood stars, while active stars appear to be concentrated in an activity strip on the log  R '_HK × [Fe/H] diagram. We provide some explanations for these trends, and show that the chromospheric emission–age relation probably has different slopes on the two sides of the Vaughan–Preston gap.     

大质量双星系统的非守恒演化

由于大质量双星系统有强大的星风物质损失,因而在研究其结构和演化时必须考虑星风物质损失,动量损失,物质交换以及由以上原因引起的轨道参量的变化,此外,天文观测又证实,一些大质量双星系统中存在星风冲击波,有Ｘ射线辐射以及有致密天体（白矮星,中子星）的存在,因此在研究大质量双星的演化时,又会遇到在星风冲击波理论及其对演化的影响,双星系统何时会演化成为公共外壳的系统,以及双星系统中如果发生超新星爆发,是否会