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

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

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.     

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

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.     

Evolutionary calculations of carbon dredge-up in helium envelope white dwarfs

We investigate the evolution of cooling helium atmosphere white dwarfs using a full evolutionary code, specifically developed to follow the effects of element diffusion and gravitational settling on white dwarf cooling. The major difference between this work and previous work is that we use more recent opacity data from the OPAL project. Since, in general, these opacities are higher than those available 10 years ago, at a given effective temperature, convection zones go deeper than in models with older opacity data. Thus convective dredge-up of observationally detectable carbon in helium atmosphere white dwarfs can occur for thicker helium layers than found by Pelletier et al. We find that the range of observed C to He ratios in different DQ white dwarfs of similar effective temperature is well explained by a range of initial helium layer mass between 10⁻³ and 10⁻² M⊙, in good agreement with stellar evolution theory, assuming a typical white dwarf mass of 0.6 M⊙. We also predict that oxygen will be present in DQ white dwarf atmospheres in detectable amounts if the helium layer mass is near the lower limit compatible with stellar evolution theory. Determination of the oxygen abundance has the potential of providing information on the profile of oxygen in the core and hence on the important ¹²C(α,γ)¹⁶O reaction rate.     

The DODO survey – II. A Gemini direct imaging search for substellar and planetary mass companions around nearby equatorial and Northern hemisphere white dwarfs

The aim of the Degenerate Objects around Degenerate Objects (DODO) survey is to search for very low-mass brown dwarfs and extrasolar planets in wide orbits around white dwarfs via direct imaging. The direct detection of such companions would allow the spectroscopic investigation of objects with temperatures much lower  (<500 K)  than the coolest brown dwarfs currently observed. These ultra-low-mass substellar objects would have spectral types >T8.5, and so could belong to the proposed Y dwarf spectral sequence. The detection of a planet around a white dwarf would prove that such objects can survive the final stages of stellar evolution and place constraints on the frequency of planetary systems around their progenitors (with masses between 1.5 and 8   M_⊙  , i.e. early B to mid-F). This paper presents the results of a multi epoch J band common proper motion survey of 23 nearby equatorial and Northern hemisphere white dwarfs. We rule out the presence of any common proper motion companions, with limiting masses determined from the completeness limit of each observation, to 18 white dwarfs. For the remaining five targets, the motion of the white dwarf is not sufficiently separated from the non-moving background objects in each field. These targets require additional observations to conclusively rule out the presence of any common proper motion companions. From our completeness limits, we tentatively suggest that  ≲5 per cent  of white dwarfs have substellar companions with   T _eff≳ 500 K  between projected physical separations of 60–200 au.     

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

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.     

L dwarfs in the Hyades

We present the results of a proper motion survey of the Hyades to search for brown dwarfs, based on UKIRT Deep Sky Survey (UKIDSS) and Two-Micron All Sky Survey (2MASS) data. This survey covers  ∼275 deg²  to a depth of   K ∼ 15  mag, equivalent to a mass of  ∼0.05 M_⊙  assuming a cluster age of 625 Myr. The discovery of 12 L dwarf Hyades members is reported. These members are also brown dwarfs, with masses between  0.05 < M < 0.075 M_⊙  . A high proportion of these L dwarfs appear to be photometric binaries.     

A SCUBA survey of Orion – the low-mass end of the core mass function

We have re-analysed all of the Submillimetre Common User Bolometer Array (SCUBA) archive data of the Orion star-forming regions. We have put together all of the data taken at different times by different groups. Consequently, we have constructed the deepest submillimetre maps of these regions ever made. There are four regions that have been mapped: Orion A North and South, and Orion B North and South. We find that two of the regions, Orion A North and Orion B North, have deeper sensitivity and completeness limits, and contain a larger number of sources, so we concentrate on these two. We compare the data with archive data from the Spitzer Space Telescope to determine whether or not a core detected in the submillimetre is pre-stellar in nature. We extract all of the pre-stellar cores from the data and make a histogram of the core masses. This can be compared to the stellar initial mass function (IMF). We find the high-mass core mass function (CMF) follows a roughly Salpeter-like slope, just like the IMF, as seen in previous work. Our deeper maps allow us to see that the CMF turns over at,  ∼1.3 M_⊙  about a factor of 4 higher than our completeness limit. This turnover has never previously been observed, and is only visible here due to our much deeper maps. It mimics the turnover seen in the stellar IMF at  ∼0.1 M_⊙  . The low-mass side of the CMF is a power law with an exponent of, 0.35 ± 0.2 which is consistent with the low-mass slope of the young cluster IMF of 0.3 ± 0.1. This shows that the CMF continues to mimic the shape of the IMF all the way down to the lower completeness limit of these data at  ∼0.3 M_⊙  .     

Stellar forensics — I. Cooling curves

The presence of low-mass, degenerate secondaries in millisecond pulsar binaries offers the opportunity to determine an age for the binary system independently of the rotational properties of the pulsar. To this end, we present here a detailed calculation of the evolution of a grid of low-mass (< 0.05 M⊙) helium core white dwarfs. We investigate the effects of different hydrogen layer masses and provide results for well-known optical bandpasses. We supplement the OPAL opacity calculations with our own calculations for low effective temperatures ( T _eff < 6000 K) and also provide fitting formulae for the gravity as a function of mass and effective temperature. In Paper II we apply these results to individual cases.     

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 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 low-mass initial mass function in the young cluster NGC 6611

NGC 6611 is the massive young cluster (2–3 Myr) that ionizes the Eagle Nebula. We present very deep photometric observations of the central region of NGC 6611 obtained with the Hubble Space Telescope and the following filters: ACS/WFC F775W and F850LP and NIC2 F110W and F160W, loosely equivalent to ground-based IZJH filters. This survey reaches down to   I ∼ 26 mag  . We construct the initial mass function (IMF) from  ∼1.5 M_⊙  well into the brown dwarf regime (down to  ∼0.02 M_⊙  ). We have detected 30–35 brown dwarf candidates in this sample. The low-mass IMF is combined with a higher-mass IMF constructed from the ground-based catalogue from Oliveira et al. We compare the final IMF with those of well-studied star-forming regions: we find that the IMF of NGC 6611 more closely resembles that of the low-mass star-forming region in Taurus than that of the more massive Orion Nebula Cluster. We conclude that there seems to be no severe environmental effect in the IMF due to the proximity of the massive stars in NGC 6611.     

A new, efficient stellar evolution code for calculating complete evolutionary tracks

We present a new stellar evolution code and a set of results, demonstrating its capability at calculating full evolutionary tracks for a wide range of masses and metallicities. The code is fast and efficient, and is capable of following through all evolutionary phases, without interruption or human intervention. It is meant to be used also in the context of modelling the evolution of dense stellar systems, for performing live calculations for both normal star models and merger products.
The code is based on a fully implicit, adaptive-grid numerical scheme that solves simultaneously for structure, mesh and chemical composition. Full details are given for the treatment of convection, equation of state, opacity, nuclear reactions and mass loss.
Results of evolutionary calculations are shown for a solar model that matches the characteristics of the present sun to an accuracy of better than 1 per cent; a  1 M_⊙  model for a wide range of metallicities; a series of models of stellar Populations I and II, for the mass range 0.25 to  64 M_⊙  , followed from pre-main-sequence to a cool white dwarf or core collapse. An initial–final mass relationship is derived and compared with previous studies. Finally, we briefly address the evolution of non-canonical configurations, merger products of low-mass main-sequence parents.     

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

Single and binary evolution of Population III stars and their supernova explosions

We present stellar evolution calculations for Population III stars for both single- and binary-star evolutions. Our models include 10- and  16.5-M_⊙  single stars and a  10-M_⊙  model star that undergoes an episode of accretion resulting in a final mass of  16.1 M_⊙  . For comparison, we present the evolution of a solar heavy element abundance model. We use the structure from late-stage evolution models to calculate simulated supernova light curves. Light curve comparisons are made between accretion and non-accretion progenitor models, and models for single-star evolution of comparable masses. Where possible, we make comparisons to previous works. Similar investigations have been carried out, but primarily for solar or near-solar heavy metal abundance stars and not including both the evolution and the supernova explosions in one work.