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 possible companions of young radio pulsars

We discuss the formation of pulsars with massive companions in eccentric orbits. We demonstrate that the probability for a non-recycled radio pulsar to have a white dwarf as a companion is comparable to that of having an old neutron star as a companion. Special emphasis is given to PSR B1820−11 and PSR B2303+46. Based on population synthesis calculations we argue that PSR B1820−11 and PSR B2303+46 could very well be accompanied by white dwarfs with mass ≳1.1 M_⊙. For PSR B1820−11, however, we cannot exclude the possibility that its companion is a main-sequence star with a mass between ∼0.7 M_⊙ and ∼5 M_⊙.     

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.     

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.     

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.     

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.     

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.     

Two-Micron All-Sky Survey J01542930+0053266: a new eclipsing M dwarf binary system

We report on Two-Micron All-Sky Survey (2MASS) J01542930+0053266, a faint eclipsing system composed of two M dwarfs. The variability of this system was originally discovered during a pilot study of the 2MASS Calibration Point Source Working Data base. Additional photometry from the Sloan Digital Sky Survey yields an eight-passband light curve from which we derive an orbital period of  2.639 0157 ± 0.000 0016  d. Spectroscopic followup confirms our photometric classification of the system, which is likely composed of M0 and M1 dwarfs. Radial velocity measurements allow us to derive the masses  (M₁= 0.66 ± 0.03 M_⊙; M₂= 0.62 ± 0.03 M_⊙)  and radii  (R₁= 0.64 ± 0.08 R_⊙; R₂= 0.61 ± 0.09 R_⊙)  of the components, which are consistent with empirical mass–radius relationships for low-mass stars in binary systems. We perform Monte Carlo simulations of the light curves which allow us to uncover complicated degeneracies between the system parameters. Both stars show evidence of Hα emission, something not common in early-type M dwarfs. This suggests that binarity may influence the magnetic activity properties of low-mass stars; activity in the binary may persist long after the dynamos in their isolated counterparts have decayed, yielding a new potential foreground of flaring activity for next generation variability surveys.     

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 mass of the neutron star in the binary millisecond pulsar PSR J 1012 + 5307

We have measured the radial velocity variation of the white dwarf secondary in the binary system containing the millisecond pulsar PSR J 1012 + 5307. Combined with the orbital parameters of the radio pulsar, we infer a mass ratio q (≡ M ₁/ M ₂) = 10.5 ± 0.5. Our optical spectroscopy has also allowed us to determine the mass of the white dwarf companion by fitting the spectrum to a grid of DA model atmospheres: we estimate M ₂ = 0.16 ± 0.02 M⊙, and hence the mass of the neutron star is 1.64 ± 0.22 M⊙, where the error is dominated by that of M ₂. The orbital inclination is 52 ± 4°. For an initial neutron star mass of ∼ 1.4 M⊙, only a few tenths of a solar mass at most has been successfully accreted over the lifetime of the progenitor low-mass X-ray binary. If the initial mass of the secondary was ∼ 1 M⊙, our result suggests that the mass transfer may have been non-conservative.     

The mass and radius of the M dwarf companion to GD 448

We present spectroscopy and photometry of GD 448, a detached white dwarf – M dwarf binary with a period of 2.47 h. We find that the Na  I  8200-Å feature is composed of narrow emission lines, owing to irradiation of the M dwarf by the white dwarf, within broad absorption lines that are essentially unaffected by heating. Combined with an improved spectroscopic orbit and gravitational redshift measurement from spectra of the Hα line, we are able to derive masses for the white dwarf and M dwarf directly (0.41 ± 0.01 and 0.096 ± 0.004 M_⊙, respectively). We use a simple model of the Ca II emission lines to establish the radius of the M dwarf assuming the emission from its surface to be proportional to the incident flux per unit area from the white dwarf. The radius derived is 0.125 ± 0.020 R_⊙. The M dwarf appears to be a normal main-sequence star in terms of its mass and radius, and is less than half the size of its Roche lobe. The thermal time-scale of the M dwarf is much longer than the cooling age of the white dwarf, so we conclude that the M dwarf was unaffected by the common-envelope phase. The anomalous width of the Hα emission from the M dwarf remains to be explained, but the strength of the line may be due to X-ray heating of the M dwarf owing to accretion on to the white dwarf from the M dwarf wind.     

Discovery of a widely separated ultracool dwarf–white dwarf binary

We present the discovery of the widest known ultracool dwarf–white dwarf binary. This binary is the first spectroscopically confirmed widely separated system from our target sample. We have used the Two-Micron All-Sky Survey (2MASS) and SuperCOSMOS archives in the southern hemisphere, searching for very widely separated ultracool dwarf–white dwarf binaries, and find one common proper motion system, with a separation of 3650–5250 au at an estimated distance of 41–59 pc, making it the widest known system of this type. Spectroscopy reveals 2MASS J0030−3740 is a DA white dwarf with   T _eff= 7600 ± 100 K, log( g ) = 7.79–8.09  and   M _WD= 0.48–0.65 M_⊙  . We spectroscopically type the ultracool dwarf companion (2MASS J0030−3739) as M9 ± 1 and estimate a mass of  0.07–0.08 M_⊙,  T _eff= 2000–2400 K  and  log( g ) = 5.30–5.35  , placing it near the mass limit for brown dwarfs. We estimate the age of the system to be >1.94 Gyr (from the white dwarf cooling age and the likely length of the main-sequence lifetime of the progenitor) and suggest that this system and other such wide binaries can be used as benchmark ultracool dwarfs.     

KPD 0422 + 5421: a new short-period subdwarf B/white dwarf binary

The subdwarf B (sdB) star KPD 0422 + 5421 was discovered to be a single-lined spectroscopic binary with a period of P  = 0.090 1795 ± (3 × 10⁻⁷) d (2 h 10 min). The U B light curves display an ellipsoidal modulation with amplitudes of ≈ 0.02 mag. The sdB star contributes nearly all of the observed flux. This and the absence of any reflection effect suggest that the unseen companion star is small (i.e. R _comp ≈ 0.01 R_⊙) and therefore degenerate. We modelled the U B light curves and derived i  = 78.05° ± 0.50° and a mass ratio of q  =  M _comp/ M _sdB = 0.87 ± 0.15. The sdB star fills 69 per cent of its Roche lobe. These quantities may be combined with the mass function of the companion [ f ( M ) = 0.126 ± 0.028 M_⊙] to derive M _sdB = 0.72 ± 0.26 M_⊙ and M _comp = 0.62 ±  0.18 M_⊙. We used model spectra to derive the effective temperature, surface gravity and helium abundance of the sdB star. We found T _eff = 25 000 ± 1500 K, log g  = 5.4 ± 0.1 and [He/H] = −1.0. With a period of 2 h 10 min, KPD 0422 + 5421 has one of the shortest known orbital periods of a detached binary. This system is also one of only a few known binaries that contain a subdwarf B star and a white dwarf. Thus KPD 0422 + 5421 represents a relatively unobserved, and short-lived, stage of binary star evolution.     

Discovery of 10 pulsars in an Arecibo drift-scan survey

We present the results of a 430-MHz survey for pulsars conducted during the upgrade to the 305-m Arecibo radio telescope. Our survey covered a total of 1147 deg² of sky using a drift-scan technique. We detected 33 pulsars, 10 of which were not known prior to the survey observations. The highlight of the new discoveries is PSR J0407+1607, which has a spin period of 25.7 ms, a characteristic age of 1.5 Gyr and is in a 1.8-yr orbit about a low-mass  (>0.2 M_⊙)  companion. The long orbital period and small eccentricity  ( e = 0.0009)  make the binary system an important new addition to the ensemble of binary pulsars suitable to test for violations of the strong equivalence principle. We also report on our initially unsuccessful attempts to detect optically the companion to J0407+1607, which imply that its absolute visual magnitude is >12.1. If, as expected on evolutionary grounds, the companion is an He white dwarf, our non-detection implies a cooling age of least 1 Gyr.     

New Praesepe white dwarfs and the initial mass–final mass relation

We report the spectroscopic confirmation of four further white dwarf members of Praesepe. This brings the total number of confirmed white dwarf members to 11, making this the second largest collection of these objects in an open cluster identified to date. This number is consistent with the high-mass end of the initial mass function of Praesepe being Salpeter in form. Furthermore, it suggests that the bulk of Praesepe white dwarfs did not gain a substantial recoil kick velocity from possible asymmetries in their loss of mass during the asymptotic giant branch phase of evolution. By comparing our estimates of the effective temperatures and the surface gravities of WD0833+194, WD0840+190, WD0840+205 and WD0843+184 to modern theoretical evolutionary tracks, we have derived their masses to be in the range  0.72–0.76 M_⊙  and their cooling ages ∼300 Myr. For an assumed cluster age of 625 ± 50 Myr, the inferred progenitor masses are between 3.3 and  3.5 M_⊙  . Examining these new data in the context of the initial mass–final mass relation, we find that it can be adequately represented by a linear function  ( a ₀= 0.289 ± 0.051,  a ₁= 0.133 ± 0.015)  over the initial mass range 2.7–6  M_⊙  . Assuming an extrapolation of this relation to larger initial masses is valid and adopting a maximum white dwarf mass of  1.3 M_⊙  , our results support a minimum mass for core-collapse supernovae progenitors in the range  ∼6.8–8.6 M_⊙  .     

Mass transfer in eccentric binaries: the new oil-on-water smoothed particle hydrodynamics technique

To measure the onset of mass transfer in eccentric binaries, we have developed a two-phase smoothed particle hydrodynamics (SPH) technique. Mass transfer is important in the evolution of close binaries, and a key issue is to determine the separation at which mass transfer begins. The circular case is well understood and can be treated through the use of the Roche formalism. To treat the eccentric case, we use a newly developed two-phase system. The body of the donor star is made up from high-mass water particles, whilst the atmosphere is modelled with low-mass oil particles. Both sets of particles take part fully in SPH interactions. To test the technique, we model circular mass-transfer binaries containing a  0.6 M_⊙  donor star and a  1 M_⊙  white dwarf; such binaries are thought to form cataclysmic variable (CV) systems. We find that we can reproduce a reasonable CV mass-transfer rate, and that our extended atmosphere gives a separation that is too large by approximately 16 per cent, although its pressure scale height is considerably exaggerated. We use the technique to measure the semimajor axis required for the onset of mass transfer in binaries with a mass ratio of   q = 0.6  and a range of eccentricities. Comparing to the value obtained by considering the instantaneous Roche lobe at pericentre, we find that the radius of the star required for mass transfer to begin decreases systematically with increasing eccentricity.     

Mass transfer in eccentric binary stars

The concept of Roche lobe overflow is fundamental to the theory of interacting binaries. Based on potential theory, it is dependent on all the relevant material corotating in a single frame of reference. Therefore if the mass losing star is asynchronous with the orbital motion or the orbit is eccentric, the simple theory no longer applies and no exact analytical treatment has been found. We use an analytic approximation whose predictions are largely justified by smoothed particle hydrodynamic simulations (SPH). We present SPH simulations of binary systems with the same semi-major axis   a = 5.55 R_⊙  , masses   M ₁= 1 M_⊙, M ₂= 2 M_⊙  and radius   R ₁= 0.89 R_⊙  for the primary star but with different eccentricities   e = 0.4, 0.5, 0.6  and 0.7. In each case the secondary star is treated as a point mass. When   e = 0.4  no mass is lost from the primary while at   e = 0.7  catastrophic mass transfer, partly through the L₂ point, takes place near periastron. This would probably lead to common-envelope evolution if star 1 were a giant or to coalescence for a main-sequence star. In between, at   e ≥ 0.5  , some mass is lost through the L₁ point from the primary close to periastron. However, rather than being all accreted by the secondary, some of the stream appears to leave the system. Our results indicate that the radius of the Roche lobe is similar to circular binaries when calculated for the separation and angular velocity at periastron. Part of the mass loss occurs through the L₂ point.     

Superhumps in the helium dwarf nova KL Draconis

We report on the discovery of a 25.5-min superhump period for the suspected helium dwarf nova system KL Draconis in a high state. The presence of superhumps combined with the previously observed helium spectrum and large-amplitude photometric variations confirm that KL Dra is an AM CVn system similar to CR Bootis, V803 Cen and CP Eridani. We also find a low-state photometric period at 25.0 min that we suggest may be the orbital period. With this assumption, we estimate   q =0.075  ,   M ₁=0.76 M_⊙  and   M ₂=0.057 M_⊙  .