Cooling curves and initial models for low-mass white dwarfs (<0.25 M⊙) with helium cores

We present a detailed calculation of the evolution of low-mass (<0.25 M_⊙) helium white dwarfs. These white dwarfs (the optical companions to binary millisecond pulsars) are formed via long-term, low-mass binary evolution. After detachment from the Roche lobe, the hot helium cores have a rather thick hydrogen layer with mass between 0.01 and 0.06 M_⊙. As a result of mixing between the core and outer envelope, the surface hydrogen content ( X _surf) is 0.5–0.35 , depending on the initial value of the heavy element Z and the initial secondary mass. We found that the majority of our computed models experience one or two hydrogen shell flashes. We found that the mass of the helium dwarf in which the hydrogen shell flash occurs depends on the chemical composition. The minimum helium white dwarf mass in which a hydrogen flash takes place is 0.213 M_⊙ ( Z =0.003), 0.198 M_⊙ ( Z =0.01), 0.192 M_⊙ ( Z =0.02) or 0.183 M_⊙ ( Z =0.03). The duration of the flashes (independent of chemical composition) is between a few ×10⁶ and a few ×10⁷ yr. In several flashes the white dwarf radius will increase so much that it forces the model to fill its Roche lobe again. Our calculations show that the cooling history of the helium white dwarf depends dramatically on the thickness of the hydrogen layer. We show that the transition from a cooling white dwarf with a temporarily stable hydrogen-burning shell to a cooling white dwarf in which almost all residual hydrogen is lost in a few thermal flashes (via Roche lobe overflow) occurs between 0.183 and 0.213 M_⊙ (depending on the heavy element value).     

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.     

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.     

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.     

The origin of the magnetic fields in white dwarfs

Magnetic white dwarfs with fields in excess of ∼10⁶ G (the high field magnetic white dwarfs; HFMWDs) constitute about ∼10 per cent of all white dwarfs and show a mass distribution with a mean mass of  ∼0.93 M_⊙  compared to  ∼0.56 M_⊙  for all white dwarfs. We investigate two possible explanations for these observations. First, that the initial–final mass relationship (IFMR) is influenced by the presence of a magnetic field and that the observed HFMWDs originate from stars on the main sequence that are recognized as magnetic (the chemically peculiar A and B stars). Secondly, that the IFMR is essentially unaffected by the presence of a magnetic field, and that the observed HFMWDs have progenitors that are not restricted to these groups of stars. Our calculations argue against the former hypothesis and support the latter. The HFMWDs have a higher than average mass because on the average they have more massive progenitors and not because the IFMR is significantly affected by the magnetic field. A requirement of our model is that ∼40 per cent of main-sequence stars more massive than  ∼4.5 M_⊙  must either have magnetic fields in the range of ∼10–100 G, which is below the current level of detection, or generate fields during subsequent stellar evolution towards the white dwarf phase. In the former case, the magnetic fields of the HFMWDs could be fossil remnants from the main-sequence phase consistent with the approximate magnetic flux conservation.     

The impact of element diffusion on the formation and evolution of helium white dwarf stars

The aim of this work is to investigate the effect of element diffusion on the evolution of helium white dwarfs. To this end, we couple the multicomponent flow equations that describe gravitational settling, chemical and thermal diffusion to an evolutionary code. We compute the evolution of a set of helium white dwarf models with masses ranging from 0.169 to 0.406 M_⊙. In particular, several low-mass white dwarfs have been found in binary systems as companion to millisecond pulsars. In these systems, pulsar emission is activated by mass transfer episodes so that, if we place the zero-age point at the end of such mass transfer, then the pulsar and the white dwarf ages should be equal. Interestingly enough, available models of helium white dwarfs neglect element diffusion. Using such models, good agreement has been found between the ages of the components of the PSR J1012+5307 system. However, recent observations of the PSR B1855+09 system cast doubts on the correctness of such models, which predict a white dwarf age twice as long as the spin-down age of the pulsar. In this work, we find that element diffusion induces thermonuclear hydrogen shell flashes for models in the mass interval 0.18≲ M /M_⊙ ≲ 0.41 . We show, in particular, that the occurrence of these diffusion-induced flashes eventually leads to white dwarf models with hydrogen envelope masses too small to support any further nuclear burning, thus implying much shorter cooling ages than in the case when diffusion is neglected. In particular, excellent agreement is found between the ages of PSR B1855+09 system components, solving the age discrepancy from first principles.     

The triple degenerate star WD 1704+481

WD 1704+481 is a visual binary in which both components are white dwarfs. We present spectra of the H α line of both stars which show that one component (WD 1704+481.2=Sanduleak B=GR 577) is a close binary with two white dwarf components. Thus, WD 1704+481 is the first known triple degenerate star. From radial velocity measurements of the close binary we find an orbital period of 0.1448 d, a mass ratio, q M _bright M _faint, of 0.70±0.03 and a difference in the gravitational redshifts of 11.5±2.3 km s⁻¹. The masses of the close pair of white dwarfs predicted by the mass ratio and gravitational redshift difference combined with theoretical cooling curves are 0.39±0.05 and 0.56±0.07 M_⊙. WD 1704+481 is therefore also likely to be the first example of a double degenerate in which the less massive white dwarf is composed of helium and the other white dwarf is composed of carbon and oxygen.     

Sakurai's object: the ionized nebula at radio wavelengths

Sakurai's object (V4334 Sgr) is a planetary nebula nucleus which is undergoing its final helium shell flash. This is the first of these rare and important events to be observable with non-optical instruments. We report the first radio detection, using a short (2-h) observation with the Very Large Array (VLA) at 4.86 GHz. The radio emission structure is coincident with the 34-arcsec diameter planetary nebula seen in optical emission lines. We find a statistical distance ∼ 3.8 ± 0.6 kpc, with a range of 1.9 <  D  < 5.3 kpc, depending on the planetary nebula (PN) mass. While we have no direct evidence for a new (post-flash) stellar wind, we estimate an upper limit to the mass-loss rate due to any such wind of 1.7 × 10⁻⁷ M⊙ yr⁻¹. The number of emitting knots in the radio-visible nebula indicates an electron density of ∼ 2 × 10⁸ m⁻³ in those knots, and a total emitting ionized mass of ∼ 0.15 M⊙, at an assumed distance of 3.8 kpc. The radio flux density indicates an Hβ flux of ∼ 6 × 10⁻¹⁶ W m⁻², suggesting an extinction E ( B  −  V ) ∼ 1.15, comparable with reddening estimates in the direction of V4334 Sgr.     

An extremely slow nova?

The binary companion to the peculiar F supergiant HD 172481 is shown to be a Mira variable with a pulsation period of 312 d. Its characteristics are within the normal range found for solitary Miras of that period, although its pulsation amplitude and mass-loss rate ̇ ∼3×10⁻⁶ M_⊙ yr⁻¹ are higher than average. Reasons are given for suspecting that the F supergiant, which has L ∼10⁴ L_⊙, is a white dwarf burning hydrogen accreted from its companion.     

The case and fate of HD 75767 – neutron star or supernova?

We report the discovery of the nearby  ( d = 24 pc)  HD 75767 as an eight billion year old quadruple system consisting of a distant M dwarf pair, HD 75767 C–D, in orbit around the known short-period   P = 10.25 d  single-lined binary HD 75767 A–B, the primary of which is a solar-like G star. On the reasonable assumption of synchronous orbital rotation as well as rotational and orbital coplanarity for the inner pair, we get   M _B= 0.96 M_⊙  for the unseen HD 75767 B, that is, the case of a massive white dwarf. Upon future evolution, mass transfer towards HD 75767 B will render the   M _A= 0.96 M_⊙  G-type primary, now a turnoff star, to become a helium white dwarf of   M _A∼ 0.33 M_⊙  . Depending on the mass accretion rate, accretion efficiency and composition of the massive white dwarf, this in turn may result in a collapse of HD 75767 B with the formation of a millisecond pulsar, i.e. the creation of a low-mass binary pulsar (LMBP), or, instead, a Type Ia supernova explosion and the complete disruption of HD 75767 B. Irrespective of which scenario applies, we point to the importance of the distant M dwarfs as the likely agents for the formation of the inner, short-period HD 75767 A–B pair, and hence a path that particularly avoids preceding phases of common envelope evolution.     

A top-heavy stellar initial mass function in starbursts as an explanation for the high mass-to-light ratios of ultra-compact dwarf galaxies

It has been recently shown that the dynamical V -band mass-to-light ratios of compact stellar systems with masses from 10⁶ to  10⁸ M_⊙  are not consistent with the predictions from simple stellar population models. Top-heavy stellar initial mass functions (IMFs) in these so-called ultra-compact dwarf galaxies (UCDs) offer an attractive explanation for this finding, the stellar remnants and retained stellar envelopes providing the unseen mass. We therefore construct a model which quantifies by how much the IMFs of UCDs would have to deviate in the intermediate- and high-mass range from the canonical IMF in order to account for the enhanced   M / L_V   ratio of the UCDs. The deduced high-mass IMF in the UCDs depends on the age of the UCDs and the number of faint products of stellar evolution retained by them. Assuming that the IMF in the UCDs is a three-part power law equal to the canonical IMF in the low-mass range and taking 20 per cent as a plausible choice for the fraction of the remnants of high-mass stars retained by UCDs, the model suggests the exponent of the high-mass IMF to be ≈1.6 if the UCDs are  13 Gyr  old (i.e. almost as old as the Universe) or ≈1.0 if the UCDs are  7 Gyr  old, in contrast to 2.3 for the Salpeter–Massey IMF. If the IMF was as top heavy as suggested here, the stability of the UCDs might have been threatened by heavy mass loss induced by the radiation and evolution of massive stars. The central densities of UCDs must have been in the range  10⁶ to 10⁷ M_⊙ pc⁻³  when they formed with star formation rates of  10 to 100 M_⊙ yr⁻¹  .     

Sakurai's Object: characterizing the near-infrared CO ejecta between 2003 and 2007

We present observations of Sakurai's Object obtained at 1–5 μm between 2003 and 2007. By fitting a radiative transfer model to an echelle spectrum of CO fundamental absorption features around  4.7 μm  , we determine the excitation conditions in the line-forming region. We find  ¹²C/¹³C = 3.5^+2.0_−1.5  , consistent with CO originating in ejecta processed by the very late thermal pulse, rather than in the pre-existing planetary nebula. We demonstrate the existence of  2.2 × 10⁻⁶≤ M _CO≤ 2.7 × 10⁻⁶ M_⊙  of CO ejecta outside the dust, forming a high-velocity wind of  500 ± 80 km s⁻¹  . We find evidence for significant weakening of the CO band and cooling of the dust around the central star between 2003 and 2005. The gas and dust temperatures are implausibly high for stellar radiation to be the sole contributor.     

Low- and intermediate-mass close binary evolution and the initial–final mass relation

Using Eggleton's stellar evolution code, we carry out 150 runs of Population I binary evolution calculations with the initial primary mass between 1 and 8 M_⊙, the initial mass ratio     between 1.1 and 4, and the onset of Roche lobe overflow (RLOF) at an early, middle or late Hertzsprung-gap stage. We assume that RLOF is conservative in the calculations, and find that the remnant mass of the primary may change by more than 40 per cent over the range of initial mass ratio or orbital period, for a given primary mass. This is contrary to the often-held belief that the remnant mass depends only on the progenitor mass if mass transfer begins in the Hertzsprung gap. We fit a formula, with an error less than 3.6 per cent, for the remnant (white dwarf) mass as a function of the initial mass M _1i of the primary, the initial mass ratio q _i and the radius of the primary at the onset of RLOF. We also find that a carbon–oxygen white dwarf with mass as low as 0.33 M_⊙ may be formed if the initial mass of the primary is around 2.5 M_⊙.     

The evolution of a rapidly accreting helium white dwarf to become a low-luminosity helium star

We have examined the evolution of merged low-mass double white dwarfs which become low-luminosity (or high-gravity) extreme helium stars. We have approximated the merging process by the rapid accretion of matter, consisting mostly of helium, on to a helium 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. As the helium-burning shell moves inwards with repeating shell flashes, the effective temperature gradually increases as the star evolves towards the helium main sequence. When the mass interior to the helium‐burning shell is approximately 0.25 M_⊙, the star enters a regime where it is pulsationally unstable. We have obtained radial pulsation periods for these models.
These models have properties very similar to those of the pulsating helium star V652 Her. We have compared the rate of period change of the theoretical models with that observed in V652 Her, as well as with its position on the Hertzsprung–Russell diagram. We conclude that the merger between two helium white dwarfs can produce a star with properties remarkably similar to those observed in at least one extreme helium star, and is a viable model for their evolutionary origin. Such helium stars will evolve to become hot subdwarfs close to the helium main sequence. We also discuss the number of low-luminosity helium stars in the Galaxy expected for our evolution scenario.     

The polarimetric nature of HD 108

V -band polarimetric observations of HD 108 were obtained during 1994 August and September. This is the first time that the temporal polarimetric variability of this star has been investigated. Its percentage polarization and position angle vary randomly, and no evidence of its proposed binary nature is detected. The random variability is consistent with that seen for previously observed Wolf–Rayet stars and OB supergiants. Nightly variations do show some systematic behaviour consistent with the blob ejection model of Underhill & Fehey (1984). From the data it is estimated that the mass-loss rate resulting from blobs is ∼1-10⁻⁷ M⊙ yr⁻¹. The total mass-loss rate of the star is estimated to be ∼1-10⁻⁵ M⊙ yr⁻¹. The stellar rotation rate is estimated to be V _rot∼400 km s⁻¹ with an inclination of i <20c. The data are also interpreted in terms of perturbations occurring in an equatorial disc plus bipolar jets viewed equatorially (Underhill 1994). It is found that the perturbation mass-loss rate is about ∼3-10⁻⁷ M⊙ yr⁻¹, giving a total mass-loss rate for the star of 13-10⁻⁵ M⊙ yr⁻¹. Both models are consistent with HD 108 being an OB supergiant or of the class Ofpe/WR9.     

Impact of cosmic rays on Population III star formation

We explore the implications of a possible cosmic-ray (CR) background generated during the first supernova explosions that end the brief lives of massive Population III stars. We show that such a CR background could have significantly influenced the cooling and collapse of primordial gas clouds in minihaloes around redshifts of   z ∼ 15–20  , provided the CR flux was sufficient to yield an ionization rate greater than about 10⁻¹⁹ s⁻¹ near the centre of the minihalo. The presence of CRs with energies  ≲10⁷  eV would indirectly enhance the molecular cooling in these regions, and we estimate that the resulting lower temperatures in these minihaloes would yield a characteristic stellar mass as low as  ∼10 M_⊙  . CRs have a less-pronounced effect on the cooling and collapse of primordial gas clouds inside more massive dark matter haloes with virial masses  ≳10⁸ M_⊙  at the later stages of cosmological structure formation around   z ∼ 10–15  . In these clouds, even without CR flux the molecular abundance is already sufficient to allow cooling to the floor set by the temperature of the cosmic microwave background.     

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.     

Boundary layer on the surface of a neutron star

In an attempt to model the accretion on to a neutron star in low-mass X-ray binaries, we present 2D hydrodynamical models of the gas flow in close vicinity of the stellar surface. First, we consider a gas pressure-dominated case, assuming that the star is non-rotating. For the stellar mass we take   M _star= 1.4 × 10⁻² M_⊙  and for the gas temperature   T = 5 × 10⁶ K  . Our results are qualitatively different in the case of a realistic neutron star mass and a realistic gas temperature of T ≃ 10⁸ K, when the radiation pressure dominates. We show that to get the stationary solution in a latter case, the star most probably has to rotate with the considerable velocity.     

Active galactic nuclei jet mass loading and truncation by stellar winds

Active galactic nuclei can produce extremely powerful jets. While tightly collimated, the scale of these jets and the stellar density at galactic centres implies that there will be many jet/star interactions, which can mass load the jet through stellar winds. Previous work employed modest wind mass outflow rates, but this does not apply when mass loading is provided by a small number of high mass-loss stars. We construct a framework for jet mass loading by stellar winds for a broader spectrum of wind mass-loss rates than has previously been considered. Given the observed stellar mass distributions in galactic centres, we find that even highly efficient (0.1 Eddington luminosity) jets from supermassive black holes of masses M _BH≲ 10⁴ M_⊙ are rapidly mass loaded and quenched by stellar winds. For  10⁴ M_⊙ < M _BH < 10⁸ M_⊙  , the quenching length of highly efficient jets is independent of the jet's mechanical luminosity. Stellar wind mass loading is unable to quench efficient jets from more massive engines, but can account for the observed truncation of the inefficient M87 jet, and implies a baryon-dominated composition on scales ≳2 kpc therein even if the jet is initially pair plasma dominated.     

PSR J1012+5307: younger than it looks?

Lorimer et al. have recently reported that the spin-down age (∼7 × 10⁹ yr) of the low-mass binary pulsar PSR J1012+5307 is much higher than the cooling age (3 × 10⁸ yr) of its white dwarf companion. The proposed solutions for this discrepancy are outlined and discussed. In particular, the revised cooling age estimate proposed by Alberts et al. agrees with data from other low-mass binary pulsar systems if a transition to the 'classical' cooling regime occurs between ∼0.14 and ∼0.28 M_⊙. If this transition is excluded, PSR J1012+5307 seems to have finished its accretion phase far from the spin-up line.