Δ Y/Δ Z from fine structure in the main sequence based on Hipparcos parallaxes

The slope Δ Y /Δ Z is a quantity of interest in relation to stellar evolution, the initial mass function and the determination of the primordial helium abundance. In this paper we estimate Δ Y /Δ Z from fine structure in the main sequence of nearby stars from Hipparcos data for stars with Z  ≤ Z⊙. We find a value of about 3, which is consistent with what has been found in extragalactic H  II regions and with stellar models for suitable upper limits to the initial masses of supernovae according to the initial mass function slope adopted.     

Stellar evolution models for Z = 0.0001 to 0.03

We have calculated a grid of empirically well tested evolutionary tracks with masses M between 0.5 and 50 M⊙, spaced by approximately 0.1 in log M , and with metallicities Z  = 0.0001, 0.0003, 0.001, 0.004, 0.01, 0.02 and 0.03. We use a robust and fast evolution code with a self-adaptive non-Lagrangian mesh, which employs the mixing-length theory but treats convective mixing as a diffusion process, solving simultaneously for the structure and the chemical composition. The hydrogen and helium abundances are chosen as functions of the metallicity: X  = 0.76 − 3.0 Z Y  = 0.24 + 2.0 Z .   Two sets of models were computed, one without and one with a certain amount of enhanced mixing or 'overshooting'. This amount has been empirically chosen by means of various sensitive tests for overshooting: (1) the luminosity of core helium burning (blue loop) giants of well-known mass, (2) the width of the main sequence as defined by double-lined eclipsing binaries with well-measured masses and radii, and (3) the shape and implied stellar distribution of isochrones of various open clusters. The first two tests have been the subject of previous papers, the third test is discussed in this paper. On the basis of these tests, we recommend the use of the overshooting models for masses above about 1.5M ⊙.   We describe here the characteristics of the models, the procedure for constructing isochrones for arbitrary age and metallicity from the models, and the performance of these isochrones for several intermediate-age and old open clusters. All original models are available in electronic form and we describe the means by which they may be obtained.     

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 ages and colours of cool helium-core white dwarf stars

The purpose of this work is to explore the evolution of helium-core white dwarf stars in a self-consistent way with the predictions of detailed non-grey model atmospheres and element diffusion. To this end, we consider helium-core white dwarf models with stellar masses of 0.406, 0.360, 0.327, 0.292, 0.242, 0.196 and 0.169 M_⊙ and follow their evolution from the end of mass-loss episodes, during their pre-white dwarf evolution, down to very low surface luminosities.
We find that when the effective temperature decreases below 4000 K, the emergent spectrum of these stars becomes bluer within time-scales of astrophysical interest. In particular, we analyse the evolution of our models in the colour–colour and in the colour–magnitude diagrams and find that helium-core white dwarfs with masses ranging from ∼0.18 to 0.3 M_⊙ can reach the turn-off in their colours and become blue again within cooling times much less than 15 Gyr and then remain brighter than M _V ≈16.5 . In view of these results, many low-mass helium white dwarfs could have had enough time to evolve to the domain of collision-induced absorption from molecular hydrogen, showing blue colours.     

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.     

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.     

Evolutionary stellar population synthesis at high spectral resolution: optical wavelengths

We present the single stellar population (SSP) synthesis results of our new synthetic stellar atmosphere models library with a spectral sampling of 0.3 Å, covering the wavelength range from 3000 to 7000 Å for a wide range of metallicities (twice solar, solar, half solar and 1/10 solar). The stellar library is composed of 1650 spectra computed with the latest improvements in stellar atmospheres. In particular, it incorporates non-local thermodynamic equilibrium (LTE) line-blanketed models for hot  ( T _eff≥ 27 500 K)  , and LTE line-blanketed models (Phoenix) for cool  (3000 ≤ T _eff≤ 4500 K)  stars. Because of the high spectral resolution of this library, evolutionary synthesis models can be used to predict the strength of numerous weak absorption lines and the evolution of the profiles of the strongest lines over a wide range of ages. The SSP results have been calculated for ages from 1 Myr to 17 Gyr using the stellar evolutionary tracks provided by the Geneva and Padova groups. For young stellar populations, our results have a very detailed coverage of high-temperature stars with similar results for the Padova and Geneva isochrones. For intermediate and old stellar populations, our results, once degraded to a lower resolution, are similar to the ones obtained by other groups (limitations imposed by the stellar evolutionary physics notwidthstanding). The limitations and advantages of our models for the analysis of integrated populations are described. The full set of the stellar library and the evolutionary models are available for retrieval at the websites http://www.iaa.csic.es/~rosa and http://www.iaa.csic.es/~mcs/sed@ , or on request from the first two authors.     

Full evolution of low-mass white dwarfs with helium and oxygen cores

We study the full evolution of low-mass white dwarfs with helium and oxygen cores. We revisit the age dichotomy observed in many white dwarf companions to millisecond pulsar on the basis of white dwarf configurations derived from binary evolution computations. We evolve 11 dwarf sequences for helium cores with final masses of 0.1604, 0.1869, 0.2026, 0.2495, 0.3056, 0.3333, 0.3515, 0.3844, 0.3986, 0.4160 and  0.4481 M_⊙  . In addition, we compute the evolution of five sequences for oxygen cores with final masses of 0.3515, 0.3844, 0.3986, 0.4160 and  0.4481 M_⊙  . A metallicity of   Z = 0.02  is assumed. Gravitational settling, chemical and thermal diffusion are accounted for during the white dwarf regime. Our study reinforces the result that diffusion processes are a key ingredient in explaining the observed age and envelope dichotomy in low-mass helium-core white dwarfs, a conclusion we arrived at earlier on the basis of a simplified treatment for the binary evolution of progenitor stars. We determine the mass threshold where the age dichotomy occurs. For the oxygen white dwarf sequences, we report the occurrence of diffusion-induced, hydrogen-shell flashes, which, as in the case of their helium counterparts, strongly influence the late stages of white dwarf cooling. Finally, we present our results as a set of white dwarf mass–radius relations for helium and oxygen cores.     

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.     

On the metallicity of the Milky Way thin disc and photometric abundance scales

The mean metallicity of the Milky Way thin disc in the solar neighbourhood is still a matter of debate, and we recently proposed an upward revision. Our star sample was drawn from a set of solar neighbourhood dwarfs with photometric metallicities. In a very recent study, it has been suggested that our metallicity calibration, based on Geneva photometry, is biased. We show here that the effect detected is not a consequence of our adopted metallicity scale, and we confirm that our findings are robust. On the contrary, the application to Strömgren photometry of the Schuster & Nissen metallicity scale is problematic. Systematic discrepancies of  ∼0.1–0.3 dex  affect the photometric metallicity determination of metal-rich stars, on the colour interval  0.22 < b − y < 0.59  , i.e. including F and G stars. For F stars, it is shown that this is a consequence of a mismatch between the standard sequence   m ₁( b − y )  of the Hyades used by Schuster & Nissen to calibrate their metallicity scale, and the system of Olsen. It means that although the calibration of Schuster & Nissen and Olsen's photometry are intrinsically correct, they are mutually incompatible for metal-rich F-type stars. For G stars, the discrepancy is most probably the continuation of the same problem, albeit worsened by the lack of spectroscopic calibrating stars. A corrected calibration is proposed that renders the calibration of Schuster & Nissen applicable to the catalogues of Olsen. We also give a simpler calibration referenced to the Hyades sequence, valid over the same colour and metallicity ranges.     

Hipparcos and the age of the Galactic disc

We use the Hipparcos colour–magnitude diagram of field stars with Tycho colours to make a new minimum age estimate for the Galactic disc. The method is based on fits to the red envelope of subgiants in the Hipparcos colour–magnitude diagram with synthetic isochrones covering the range of disc metal abundance. The colours and luminosities of the isochrones as a function of abundance are checked using new techniques involving 'red-clump' stars in the giant branch region and on the main sequence using G and K dwarfs. We derive a minimum disc age of 8 Gyr, in good agreement with other methods.     

The role of chemical composition in RR Lyrae pulsationalproperties – I. Periods

We analyse the periods of theoretical radial pulsators, covering the range of total masses, luminosities, effective temperatures and chemical compositions expected for RR Lyrae variables in both galactic fields and globular clusters.
We show that for fixed values of the structural parameters (mass, luminosity and effective temperature), the period of fundamental and first-overtone pulsators is independent of the helium content ( Y ), whereas it slightly increases as the amount of metals ( Z ) increases. Furthermore, we find that the period along the blue edge for first-overtone pulsation is a function of mass, luminosity and helium content, with a marginal dependence on Z .
On these grounds, new linear relations connecting the periods to stellar parameters are derived. Such new relations should allow a more accurate interpretation of the RR Lyrae observed periods and, in particular, they should help in ascertaining the calibration of the mean absolute magnitude of RR Lyrae stars in terms of metal content.     

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.     

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.     

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_⊙.     

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.     

Intermediate-metallicity, high-velocity stars and Galactic chemical evolution

High signal-to-noise ratio spectra were obtained of 10 high-proper-motion stars having  −1 ≲[Fe/H] < 0  , and a comparable number of disc stars. All but two of the high-proper-motion stars were confirmed to have  [Fe/H] > −1.0  , some approaching solar metallicity, but, even so, earlier measurements overestimated the metallicities and velocities of some of these stars. Models of stellar populations were used to assign membership probabilities to the Galactic components to which the high-velocity stars might belong. Many were found to be more probably thick-disc than halo objects, despite their large space motions, and two might be associated with the inner Galaxy. It may be necessary to reassess contamination of previous halo samples, such as those used to define the metallicity distribution, to account for contamination by high-velocity thick-disc stars, and to consider possible subcomponents of the halo.
The change in [α/Fe] ratios at  [Fe/H]≃−1.0  is often used to constrain the degree and timing of Type Ia supernova nucleosynthesis in Galactic chemical-evolution models. [Ti/Fe] values were measured for eight of the high-velocity stars. Both high- and low-[Ti/Fe] halo stars exist; likewise high- and low-[Ti/Fe] thick-disc stars exist. We conclude that the [Ti/Fe]'break' is not well defined for a given population; nor is there a simple, continuous evolutionary sequence through the break. Implications for the interpretation of the [α/Fe] break in terms of SN Ia time-scales and progenitors are discussed. The range of [Ti/Fe] found for high -velocity (low rotation) thick-disc stars contrasts with that for the low -velocity (high rotation) thick-disc sample studied by Prochaska et al.     

Metallicity and photospheric abundances in field K and M dwarfs

Elemental abundances in late-type stars are of interest in several ways: they determine the location of the stars in the HR diagram and therefore their ages, as well as the atmospheric structure in their middle and upper photospheres. Especially in the case of chromospherically active late-type stars the question arises to what degree the upper photosphere is influenced by the nearby chromosphere. Analysing S/N ∼ 200 and Δλ/λ ∼ 20 000 data, we found a mean metallicity index [M/H] = −0.2 for programme K and M field stars based on an analysis of spectra in the region 5500–9000 Å. We also found that the Ca  I 6162-Å transition is a potential surface gravity indicator for K-type stars. For the chromospheric activity interval 4.4 < log  F _{Mg II}  < 6.6 we did not find any chromospheric activity impact on photospheric and upper photospheric transitions. With the derived metallicity, we confirmed the Li abundance from our previous paper and thus its dependence on the Mg  II chromospheric activity index. The nature of the spectrum for the active M-type star Gl 896A is explained by pure rotation of 14 km s⁻¹. As far as the lithium–rotation relation is concerned, the spectrum of Gl 517 is rotationally broadened as well, by 12 km s⁻¹, and the Li abundance is the second highest in our sample of stars. However, there is no link between very high Li abundance, 2.2 dex, in the K dwarf star Gl 5 and stellar rotation.     

A rigorous statistical test for the existence of supermetallicity in F, G and K dwarfs

In this paper, the first statistically rigorous test is performed to determine whether there are dwarfs in the solar neighbourhood with [Fe/H] > +0.2 dex. The possibility that no such stars exist is rejected with a confidence limit exceeding 99.994 per cent. If the threshold value is raised to +0.38 dex, the corresponding hypothesis is rejected with a confidence limit of 98.5 per cent. It is tentatively found that this upper limit does not depend on the presence or absence of planets orbiting metal-rich stars. In a result that is unchanged from a previous paper, the upper limit for the known metallicities of giants is found to be about 0.2 dex lower than the limit for dwarfs. Stars that might be observed to improve the upper metallicity limit for dwarfs are listed.     

ROSAT All-Sky Survey X-ray and EUV observations of YY Gem and AU Mic

Simultaneous X-ray and extreme ultraviolet (EUV) ( ROSAT XRT and WFC All-Sky Survey) observations of the highly active dMe flare stars YY Gem and AU Mic show that the two stars displayed an unusual type of flaring behaviour. We detect several X-ray and EUV flares superimposed on an enhanced and smoothly varying quiescent background. The two large impulsive-type X-ray flares on YY Gem reach peak X-ray luminosities of     and we estimate that they had similar integrated luminosities (∼6–8×10³³ erg). AU Mic also produced several X-ray and EUV flares, with one very impulsive flare producing a 10-fold increase in XRT count rate. This flare was even larger than the YY Gem flares (peak L _X of     and integrated L _X of    
The     ratio for both stars is at the 'saturation' limit found in rapidly rotating dwarfs and the most active RS CVn stars. We suggest that the gradually varying components are the result of a period of continuous, unresolved flaring activity. Alternatively, they may be the result of the emergence and subsequent decay of a new magnetic active region on the stellar surface of these stars.