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).     

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

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

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.     

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.     

Low- and intermediate-mass close binary evolution and the initial–final mass relation – III. Conservative case with convective overshooting and non-conservative case without overshooting

We have completed our series of calculation of case B evolution of Population I close binaries by including a conservative case with convective overshooting and a non-conservative case without convective overshooting in this paper. To apply the results easily to binary population synthesis, we fit two formulae for the remnant mass as a function of the initial mass of the primary, the radius of the primary at the onset of Roche lobe overflow (RLOF) and the initial mass ratio, with errors less than 1.9 per cent and 3.6 per cent for the conservative and non-conservative cases, respectively. We have also made comparisons between the results in this paper and those of previous studies in this series to examine the effects of our assumptions.
We find that the remnant mass depends on when RLOF begins in the Hertzsprung gap and the dependence increases with the mass of the primary. Helium–carbon–oxygen white dwarfs are formed in this series of calculation. Both the non-conservative assumption and convective overshooting make RLOF more stable: some binaries with initial mass ratios   q _i= 4.0  are dynamically stable during RLOF in all cases except for the conservative case without overshooting. From these models, we have discussed potential electron-degenerate oxygen–neon cores. Finally, we show the conditions under which one evolution scenario may be more likely than the others according to the results of previous authors.     

An iterative filter to reconstruct planetary transit signals in the presence of stellar variability

We report the identification, from a photometric, astrometric and spectroscopic study, of a massive white dwarf member of the nearby, approximately solar metallicity, Coma Berenices open star cluster (Melotte 111). We find the optical to near-infrared energy distribution of WD 1216+260 to be entirely consistent with that of an isolated DA and determine the effective temperature and surface gravity of this object to be   T _eff= 15 739⁺¹⁹⁷₋₁₉₆ K  and  log  g = 8.46^+0.03_−0.02  . We set tight limits on the mass of a putative cool companion,   M ≳ 0.036 M_⊙  (spatially unresolved) and   M ≳ 0.034 M_⊙  (spatially resolved and   a ≲ 2500 au  ). Based on the predictions of CO core, thick H layer evolutionary models we determine the mass and cooling time of WD 1216+260 to be   M _WD= 0.90 ± 0.04 M_⊙  and  τ_cool= 363⁺⁴⁶₋₄₁ Myr  , respectively. For an adopted cluster age of  τ= 500 ± 100 Myr  we infer the mass of its progenitor star to be   M _init= 4.77^+5.37_−0.97 M_⊙  . We briefly discuss this result in the context of the form of the stellar initial mass–final mass relation.     

Masses for the Local Group and the Milky Way

We use the very large Millennium Simulation of the concordance Λ cold dark matter cosmogony to calibrate the bias and error distribution of Timing Argument estimators of the masses of the Local Group and of the Milky Way. From a large number of isolated spiral–spiral pairs similar to the Milky Way/Andromeda system, we find the interquartile range of the ratio of timing mass to true mass to be a factor of 1.8, while the 5 and 95 per cent points of the distribution of this ratio are separated by a factor of 5.7. Here, we define true mass as the sum of the 'virial' masses, M ₂₀₀, of the two dominant galaxies. For present best values of the distance and approach velocity of Andromeda, this leads to a median likelihood estimate of the true mass of the Local Group of  5.27 × 10¹² M_⊙  or  log  M _LG/M_⊙= 12.72  , with an interquartile range of [12.58, 12.83] and a 5–95 per cent range of [12.26, 13.01]. Thus, a 95 per cent lower confidence limit on the true mass of the Local Group is  1.81 × 10¹² M_⊙  . A timing estimate of the Milky Way's mass based on the large recession velocity observed for the distant satellite Leo I works equally well, although with larger systematic uncertainties. It gives an estimated virial mass for the Milky Way of  2.43 × 10¹² M_⊙  with a 95 per cent lower confidence limit of  0.80 × 10¹² M_⊙  .     

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.     

Quantifying the coexistence of massive black holes and dense nuclear star clusters

In large spheroidal stellar systems, such as elliptical galaxies, one invariably finds a  10⁶–10⁹ M_⊙  supermassive black hole at their centre. In contrast, within dwarf elliptical galaxies one predominantly observes a  10⁵–10⁷ M_⊙  nuclear star cluster. To date, few galaxies have been found with both types of nuclei coexisting and even less have had the masses determined for both central components. Here, we identify one dozen galaxies housing nuclear star clusters and supermassive black holes whose masses have been measured. This doubles the known number of such hermaphrodite nuclei – which are expected to be fruitful sources of gravitational radiation. Over the host spheroid (stellar) mass range  10⁸–10¹¹ M_⊙  , we find that a galaxy's nucleus-to-spheroid (baryon) mass ratio is not a constant value but decreases from a few per cent to ∼0.3 per cent such that  log[( M _BH+ M _NC)/ M _sph]=−(0.39 ± 0.07) log[ M _sph/10¹⁰ M_⊙]− (2.18 ± 0.07)  . Once dry merging commences and the nuclear star clusters disappear, this ratio is expected to become a constant value.
As a byproduct of our investigation, we have found that the projected flux from resolved nuclear star clusters is well approximated with Sérsic functions having a range of indices from ∼0.5 to ∼3, the latter index describing the Milky Way's nuclear star cluster.     

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.     

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.     

The component masses of the cataclysmic variable V347 Puppis

We present time-resolved spectroscopy and photometry of the double-lined eclipsing cataclysmic variable V347 Pup (=LB 1800). There is evidence of irradiation on the inner hemisphere of the secondary star, which we correct for using a model to give a secondary-star radial velocity of   K _R= 198 ± 5 km s⁻¹  . The rotational velocity of the secondary star in V347 Pup is found to be   v sin  i = 131 ± 5 km s⁻¹  and the system inclination is   i = 840 ± 23  . From these parameters we obtain masses of   M ₁= 0.63 ± 0.04 M_⊙  for the white dwarf primary and   M ₂= 0.52 ± 0.06 M_⊙  for the M0.5V secondary star, giving a mass ratio of   q = 0.83 ± 0.05  . On the basis of the component masses, and the spectral type and radius of the secondary star in V347 Pup, we find tentative evidence for an evolved companion. V347 Pup shows many of the characteristics of the SW Sex stars, exhibiting single-peaked emission lines, high-velocity S-wave components and phase-offsets in the radial velocity curve. We find spiral arms in the accretion disc of V347 Pup and measure the disc radius to be close to the maximum allowed in a pressureless disc.     

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.     

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

Composite spectraPaper 9: HR 5983

HR 5983 is a sixth-magnitude composite-spectrum binary system. It has a circular orbit and a period of 108 d. We separate the component spectra and show that they are of types ∼G6 IIIa and A2.5 IV, with a mass ratio (giant:dwarf) of ∼1.04. We find that the secondary has a mass close to 2.6 M_⊙ and has already evolved somewhat from the main sequence, and that the primary appears to have reached the helium-burning 'blue loop' stage and is in synchronous rotation. The system resembles Capella in several respects, and may represent an earlier phase in the evolution of the latter.     

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.     

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⁻¹  .     

UBV photometry, UV spectroscopy and radio observations of the peculiar binary V Sagittae

We have performed high-speed UBV photometric observations on the peculiar binary V Sagittae. Using three new eclipse timings we update the orbital ephemeris and convert it to a dynamical time-scale (TDB). We also searched for quasi-periodic oscillations but did not detect them. Using the Wilson–Devinney algorithm we have modelled the light curve to find the stellar parameters of V Sge. We find that the system is a detached binary but that the primary star is very close to filling its Roche lobe, while the secondary star fills 90 per cent of its Roche lobe volume. We find temperatures of the primary and the secondary star to be T ₁=41 000 K and T ₂=22 000 K. We find i =72° and masses of 0.8 M_⊙ and 3.3 M_⊙ for the primary and secondary stars respectively. De-archived Hubble Space Telescope ( HST ) spectroscopy of V Sge shows evidence of mass loss via a wind or winds. In addition we report radio observations of V Sge during an optical high state at 2 cm, 3.6 cm and 6 cm wavelengths. The 3.6 cm emission is increased by a factor of more than six compared with an earlier detection in a previous optical high state.