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.     

Diffusion in helium white dwarf stars

This paper is aimed at exploring the effects of diffusion on the structure and evolution of low-mass helium white dwarfs. To this end, we solve the multicomponent flow equations describing gravitational settling and chemical and thermal diffusion. The diffusion calculations are coupled to an evolutionary code in order to follow the cooling of low-mass, helium core white dwarf models having envelopes made up of a mixture of hydrogen and helium, as recently suggested by detailed evolutionary calculations for white dwarf progenitors in binary systems. We find that diffusion causes hydrogen to float and the other elements to sink over time-scales shorter than evolutionary time-scales. This produces a noticeable change in the structure of the outer layers, making the star inflate. Thus, in order to compute accurately the mass–radius relation for low-mass helium white dwarfs we need to account for the diffusion processes during (at least) the white dwarf stages of the evolution of these objects. This should be particularly important when studying the general characteristics of binary systems containing a helium white dwarf and a pulsar.
In addition, we present an analytic, approximate model for the outer layers of the white dwarf aimed at interpreting the physical reasons for the change in the surface gravity for low-mass white dwarfs induced by diffusion.     

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.     

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.     

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

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

Non-adiabatic linear pulsation models for low-mass helium stars

Non-adiabatic linear pulsation models have been calculated for low-mass stars with effective temperatures between 16 000 and 35 000 K, and with surface gravities in the range 3⊙, X =0.00, Z =0.02. It is shown that the Z -bump instability persists to low masses ( M ∼0.4 M_⊙) but is suppressed either by a reduction in metallicity Z or by a selective enhancement of the carbon abundance. An unexpected result is the discovery that Z -bump instability persists at hydrogen abundances X >0.3, although the position of the red edge is sensitive to X . We have found that non-radial pulsations are also excited in the same instability region as radial pulsations.
The implications of these results for individual low-mass helium stars are discussed. It is concluded that Z -bump driven pulsations (radial and/or non-radial) may be excited in some helium-rich subdwarf B stars, representing a possible major extension to the class of variable stars represented by the prototype V652 Her.     

Models of rapidly rotating neutron stars: remnants of accretion-induced collapse

Equilibrium models of differentially rotating nascent neutron stars are constructed, which represent the result of the accretion-induced collapse of rapidly rotating white dwarfs. The models are built in a two-step procedure: (1) a rapidly rotating pre-collapse white dwarf model is constructed; (2) a stationary axisymmetric neutron star having the same total mass and angular momentum distribution as the white dwarf is constructed. The resulting collapsed objects consist of a high-density central core of size roughly 20 km, surrounded by a massive accretion torus extending over 1000 km from the rotation axis. The ratio of the rotational kinetic energy to the gravitational potential energy of these neutron stars ranges from 0.13 to 0.26, suggesting that some of these objects may have a non-axisymmetric dynamical instability that could emit a significant amount of gravitational radiation.     

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

Searching for proto‐brown dwarfs: Extending near IR spectroscopy of protostars below the hydrogen burning limit

Recent observations of nearby star forming regions have offered evidence that young brown dwarfs undergo a period of mass accretion analogous to the T Tauri phase observed in young stars. Brown dwarf analogs to stellar protostars, however, have yet to be definitively observed. These young, accreting objects would shed light on the nature of the dominant brown dwarf formation process, as well as provide ideal laboratories to investigate the dependence of the accretion mechanism on protostellar mass. Recent near infrared surveys have identified candidate proto‐brown dwarfs and characterized low mass protostars in nearby star forming regions. These techniques allow near infrared spectra to diagnose the effective temperature, accretion luminosity, magnetic field strength and rotation velocity of young low mass stars across the stellar/substellar boundary. The lowest mass proto‐brown dwarfs (M < 40 M_Jup), however, will prove challenging to observe given current near IR observational capabilities. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

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 helium abundance and ΔY/ΔZ in lower main-sequence stars

We use nearby K dwarf stars to measure the helium-to-metal enrichment ratio  Δ Y /Δ Z   , a diagnostic of the chemical history of the solar neighbourhood. Our sample of K dwarfs has homogeneously determined effective temperatures, bolometric luminosities and metallicities, allowing us to fit each star to the appropriate stellar isochrone and determine its helium content indirectly. We use a newly computed set of Padova isochrones which cover a wide range of helium and metal content.
Our theoretical isochrones have been checked against a congruous set of main-sequence binaries with accurately measured masses, to discuss and validate their range of applicability. We find that the stellar masses deduced from the isochrones are usually in excellent agreement with empirical measurements. Good agreement is also found with empirical mass-luminosity relations.
Despite fitting the masses of the stars very well, we find that anomalously low helium content (lower than primordial helium) is required to fit the luminosities and temperatures of the metal-poor K dwarfs, while more conventional values of the helium content are derived for the stars around solar metallicity.
We have investigated the effect of diffusion in stellar models and the assumption of local thermodynamic equilibrium (LTE) in deriving metallicities. Neither of these is able to resolve the low-helium problem alone and only marginally if the cumulated effects are included, unless we assume a mixing-length which is strongly decreasing with metallicity. Further work in stellar models is urgently needed.
The helium-to-metal enrichment ratio is found to be  Δ Y /Δ Z = 2.1 ± 0.9  around and above solar metallicity, consistent with previous studies, whereas open problems still remain at the lowest metallicities. Finally, we determine the helium content for a set of planetary host stars.     

The brown dwarf desert as a consequence of orbital migration

We show that the dearth of brown dwarfs in short-period orbits around Solar-mass stars – the brown dwarf desert – can be understood as a consequence of inward migration within an evolving protoplanetary disc. Brown dwarf secondaries forming at the same time as the primary star have masses which are comparable to the initial mass of the protoplanetary disc. Subsequent disc evolution leads to inward migration, and destruction of the brown dwarf, via merger with the star. This is in contrast with massive planets, which avoid this fate by forming at a later epoch when the disc is close to being dispersed. Within this model, a brown dwarf desert arises because the mass at the hydrogen-burning limit is coincidentally comparable to the initial disc mass for a Solar mass star. Brown dwarfs should be found in close binaries around very low mass stars, around other brown dwarfs, and around Solar-type stars during the earliest phases of star formation.     

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.     

On rejuvenation in massive binary systems

We introduce a set of stellar models for massive stars whose evolution has been affected by mass transfer in a binary system, at a range of metallicities. As noted by other authors, the effect of such mass transfer is frequently more than just rejuvenation. We find that, whilst stars with convective cores which have accreted only H-rich matter rejuvenate as expected, those stars which have accreted He-rich matter (e.g. at the end stages of conservative mass transfer) evolve in a way that is qualitatively similar to rejuvenated stars of much higher metallicity. Thus, the effects of non-conservative evolution depend strongly on whether He-rich matter is amongst the portion accreted or ejected. This may lead to a significant divergence in binary evolution paths with only a small difference in initial assumptions. We compare our models to observed systems and find approximate formulae for the effect of mass accretion on the effective age and metallicity of the resulting star.     

The mass ratio distribution of short-period double degenerate stars

Short-period double degenerates (DDs) are close white dwarf–white dwarf binary stars which are the result of the evolution of interacting binary stars. We present the first definitive measurements of the mass ratio for two DDs, WD 0136+768 and WD 1204+450, and an improved measurement of the mass ratio for WD 0957−666. We compare the properties of the six known DDs with measured mass ratios to the predictions of various theoretical models. We confirm the result that standard models for the formation of DDs do not predict sufficient DDs with mass ratios close to 1. We also show that the observed difference in cooling ages between white dwarfs in DDs is a useful constraint on the initial mass ratio of the binary. A more careful analysis of the properties of the white dwarf pair WD 1704+481.2 leads us to conclude that the brighter white dwarf is older than its fainter companion. This is the opposite of the usual case for DDs and is caused by the more massive white dwarf being smaller and cooling faster. The mass ratio in the sense (mass of younger star)/(mass of older star) is then  1.43±0.06  rather than the value of  0.70±0.03  given previously.     

Detection of superhumps in the VY Scl-type nova-like variable KR Aur

We report on the detection of negative superhumps in KR Aur, a typical VY Scl star. The observations were obtained with a multi-channel photometer over 107 h. The analysis of the data clearly revealed brightness variations with a period of 3.771 (±0.005) h. This is 3.5 per cent shorter than P _orb, suggesting that the observed oscillation is a negative superhump. Negative superhumps in VY Scl stars are widespread. The discovery of powerful soft X-rays from V751 Cyg suggests that VY Scl stars may contain white dwarfs, on to which nuclear burning of the accreted material occurs. If this suspicion is correct, it is possible that the powerful radiation emerging from the white dwarf may cause a tilt of the accretion disc to the orbital plane, and its retrograde precession may produce the negative superhumps seen in VY Scl stars.     

High surface brightness dwarf galaxies in the Fornax cluster

We report the discovery, in an Extreme Ultraviolet Explorer ( EUVE ) short-wavelength spectrum, of an unresolved hot white dwarf companion to the 5th magnitude B5Vp star HR 2875. This is the first time that a non-interacting white dwarf+B star binary has been discovered: previously, the earliest type of star known with a white dwarf companion was Sirius (A1V). As the white dwarf must have evolved from a main-sequence progenitor with a mass greater than that of a B5V star (≯6.0 M⊙), this places a lower limit on the maximum mass for white dwarf progenitors, with important implications for our knowledge of the initial–final mass relation. Assuming a pure-hydrogen atmospheric composition, we constrain the temperature of the white dwarf to be between 39 000 and 49 000 K. We also argue that this degenerate star is likely to have a mass significantly greater than the mean mass for white dwarf stars (≈0.55 M⊙). Finally, we suggest that other bright B stars (e.g. θ Hya) detected in the extreme ultraviolet surveys of the ROSAT Wide Field Camera and EUVE may also be hiding hot white dwarf companions.     

Kinematics of the ultracompact helium accretor AM Canum Venaticorum

We report on the results from a five-night campaign of high-speed spectroscopy of the 17-min binary AM Canum Venaticorum (AM CVn), obtained with the 4.2-m William Herschel Telescope on La Palma.
We detect a kinematic feature that appears to be entirely analogous to the 'central spike' known from the long-period, emission-line AM CVn stars GP Com, V396 Hya and SDSS J124058.03−015919.2, which has been attributed to the accreting white dwarf. Assuming that the feature indeed represents the projected velocity amplitude and phase of the accreting white dwarf, we derive a mass ratio   q = 0.18 ± 0.01  for AM CVn. This is significantly higher than the value found in previous, less direct measurements. We discuss the implications for AM CVn's evolutionary history and show that a helium star progenitor scenario is strongly favoured. We further discuss the implications for the interpretation of AM CVn's superhump behaviour, and for the detectability of its gravitational-wave signal with the Laser Interferometer Space Antenna ( LISA ).
In addition, we demonstrate a method for measuring the circularity or eccentricity of AM CVn's accretion disc, using stroboscopic Doppler tomography. We test the predictions of an eccentric, precessing disc that are based on AM CVn's observed superhump behaviour. We limit the effective eccentricity in the outermost part of the disc, where the resonances that drive the eccentricity are thought to occur, to   e = 0.04 ± 0.01  , which is smaller than previous models indicated.     

贫金属富碳天琴RR变星的形成

贫金属富碳恒星(Carbon-Enhanced Metal-Poor, CEMP)是研究宇宙早期恒星性质和化学演化的极佳样本,通常认为来自双星.目前发现的贫金属富碳星中有9颗天琴RR变星(RR Lyrae star, RRL),其中至少7颗未表现出任何双星特征.传统双星物质转移模型不足以充分解释贫金属富碳天琴RR变星(CEMP-RR Lyrae)单星的形成.之前研究表明氦白矮星和赫氏空隙星(HG)的并合模型可以解释部分富碳红巨星单星的碳增丰现象,因此贫金属富碳星单星也可能来自氦白矮星和赫氏空隙星的并合模型渠道.通过详细计算的氦白矮星和赫氏空隙星并合模型来检验这一演化渠道,结果表明:该并合模型在后续的演化过程中,其重力加速度、温度、表面碳丰度均能与观测符合较好.由此,氦白矮星和赫氏空隙星并合模型极有可能是贫金属富碳天琴RR变星的形成渠道之一.     

Very low mass stars in binaries: a theoretical look

We have undertaken a series of hydrodynamic + N ‐body simulations in order to explore the binary properties of young stars. We find that multiple stars are a natural outcome of collapsing turbulent flows, with a high incidence of N > 2 multiples, specially among the higher mass objects. We find a positive correlation of multiplicity with primary mass and a companion frequency that decreases with age, during the first few Myr after formation. Binary brown dwarfs are rarely formed, in conflict with observations. Brown dwarfs as companions are predominantly found orbiting binaries or triples at large separations. The paucity of ultra low mass and low mass ratio binaries has been investigated further, and we tentatively conclude that their formation is intricately related to an appropriate selection of initial conditions and an accurate modelling of disc accretion and evolution. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)