Evolution Of White Dwarfs Research Articles

WDEC: A Code for Modeling White Dwarf Structure and Pulsations

Abstract The White Dwarf Evolution Code (WDEC), written in Fortran, makes models of white dwarf stars. It is fast, versatile, and includes the latest physics. The code evolves hot (∼100,000 K) input models down to a chosen effective temperature by relaxing the models to be solutions of the equations of stellar structure. The code can also be used to obtain g-mode oscillation modes for the models. WDEC has a long history going back to the late 1960s. Over the years, it has been updated and re-packaged for modern computer architectures and has specifically been used in computationally intensive asteroseismic fitting. Generations of white dwarf astronomers and dozens of publications have made use of the WDEC, although the last true instrument paper is the original one, published in 1975. This paper discusses the history of the code, necessary to understand why it works the way it does, details the physics and features in the code today, and points the reader to where to find the code and a user guide.

Search for trans-iron elements in hot, helium-rich white dwarfs with the HST Cosmic Origins Spectrograph

The metal abundances in the atmospheres of hot white dwarfs (WDs) entering the cooling sequence are determined by the preceding Asymptotic Giant Branch (AGB) evolutionary phase and, subsequently, by the onset of gravitational settling and radiative levitation. In this paper, we investigate three hot He-rich WDs, which are believed to result from a late He-shell flash. During such a flash, the He-rich intershell matter is dredged up and dominates the surface chemistry. Hence, in contrast to the usual H-rich WDs, their spectra allow direct access to s-process element abundances in the intershell that were synthesized during the AGB stage. In order to look for trans-iron group elements (atomic number Z > 29), we performed a non-local thermodynamic equilibrium model atmosphere analysis of new ultraviolet spectra taken with the Cosmic Origins Spectrograph aboard the Hubble Space Telescope. One of our program stars is of PG 1159 spectral type; this star, PG 1707+427, has effective temperature Teff = 85 000 K, and surface gravity logg = 7.5. The two other stars are DO white dwarfs: WD 0111+002 has Teff = 58 000 K and log g = 7.7, and PG 0109+111 has Teff = 70 000 K and log g = 8.0. These stars trace the onset of element diffusion during early WD evolution. While zinc is the only trans-iron element we could detect in the PG 1159 star, both DOs exhibit lines from Zn, Ga, Ge, Se; one additionally exhibits lines from Sr, Sn, Te, and I and the other from As. Generally, the trans-iron elements are very abundant in the DOs, meaning that radiative levitation must be acting. Most extreme is the almost six orders of magnitude oversolar abundance of tellurium in PG 0109+111. In terms of mass fraction, it is the most abundant metal in the atmosphere. The two DOs join the hitherto unique hot DO RE 0503−289, in which 14 trans-iron elements had even been identified.

Thermal evolution of old white dwarfs

This work is devoted to a description of thermodynamic properties of Coulomb crystals which are expected to form in white dwarf interiors. Effects of magnetic field, isotopic impurities, polarization of the electron background and crystal lattice type on the thermal evolution of white dwarfs are discussed. It is shown that the electron polarization could play a noticeable role in the cooling process. While other parameters in concern do not make a significant impact.

The importance of Urca-process cooling in accreting ONe white dwarfs

We study the evolution of accreting oxygen-neon (ONe) white dwarfs (WDs), with a particular emphasis on the effects of the presence of the carbon-burning products $\mathrm{^{23}Na}$ and $\mathrm{^{25}Mg}$. These isotopes lead to substantial cooling of the WD via the $\mathrm{^{25}Mg}$-$\mathrm{^{25}Na}$, $\mathrm{^{23}Na}$-$\mathrm{^{23}Ne}$, and $\mathrm{^{25}Na}$-$\mathrm{^{25}Ne}$ Urca pairs. We derive an analytic formula for the peak Urca-process cooling rate and use it to obtain a simple expression for the temperature to which the Urca process cools the WD. Our estimates are equally applicable to accreting carbon-oxygen WDs. We use the Modules for Experiments in Stellar Astrophysics (MESA) stellar evolution code to evolve a suite of models that confirm these analytic results and demonstrate that Urca-process cooling substantially modifies the thermal evolution of accreting ONe WDs. Most importantly, we show that MESA models with lower temperatures at the onset of the $\mathrm{^{24}Mg}$ and $\mathrm{^{24}Na}$ electron captures develop convectively unstable regions, even when using the Ledoux criterion. We discuss the difficulties that we encounter in modeling these convective regions and outline the potential effects of this convection on the subsequent WD evolution. For models in which we do not allow convection to operate, we find that oxygen ignites around a density of $\log(\rho_{\rm c}/\rm g\,cm^{-3}) \approx 9.95$, very similar to the value without Urca cooling. Nonetheless, the inclusion of the effects of Urca-process cooling is an important step in producing progenitor models with more realistic temperature and composition profiles which are needed for the evolution of the subsequent oxygen deflagration and hence for studies of the signature of accretion-induced collapse.

Helium ignition in rotating magnetized CO white dwarfs leading to fast and faint rather than classical Type Ia supernovae

Based mostly on stellar models which do not include rotation, CO white dwarfs which accrete helium at rates of about $\sim 10^{-8}~\mathrm{M}_\odot/\mathrm{yr}$ have been put forward as candidate progenitors for a number of transient astrophysical phenomena, including supernovae of Type Ia, and the peculiar and fainter Type Iax supernovae. Here we study the impact of accretion-induced spin-up including the subsequent magnetic field generation, angular momentum transport, and viscous heating on the white dwarf evolution up to the point of helium ignition. We resolve the structure of the helium accreting white dwarf models with a one dimensional Langrangian hydrodynamic code, modified to include rotational and magnetic effects. We find magnetic angular momentum transport, which leads to quasi solid-body rotation, profoundly impacts the evolution of the white dwarf models. Our rotating lower mass ($0.54$ and $0.82~\mathrm{M}_\odot$) models accrete up to 50\% more mass up to ignition compared to the non-rotating case and the opposite for our more massive models. Furthermore, we find that rotation leads to up to 10-times smaller helium ignition densities, except for the lowest adopted initial white dwarf mass. Ignition densities of order $10^6\,\mathrm{g/cm}^3$ are only found for the lowest accretion rates and for large amounts of accreted helium $\gtrsim 0.4\,\mathrm{M}_\odot$. However, correspondingly massive donor stars would transfer mass at much higher rates. We expect explosive He-shell burning to mostly occur as deflagrations and at $\dot{M}>2\cdot10^{-8}~\mathrm{M}_\odot/\mathrm{yr}$, regardless of white dwarf mass. Our results imply that helium accretion onto CO white dwarfs at the considered rates is unlikely to lead to explosions of the CO core or to classical Type,Ia supernovae, but may instead produce events like faint and fast hydrogen-free supernovae.

Tidal dissipation and evolution of white dwarfs around massive black holes: an eccentric path to tidal disruption

Awhite dwarf (WD)captured into a high-eccentricity orbit around amassive black hole (MBH) may undergo many pericentre passages before tidal disruption. During these passages, the tidal potential of the MBH excites internal oscillations or waves in the WD, and the dissipation of these oscillations can significantly influence the physical properties of the WD prior to its disruption. We calculate the amplitude of the tidally excited gravity (buoyancy) waves in the WD as a function of the pericentre distance and eccentricity for realistic WD models, under the assumption that these outgoing gravity waves are efficiently dissipated in the outer layers of the WD by non-linear effects or radiative damping. We obtain fitting formulae for the tidal energy and angular momentum transfer rates as well as the tidal heating rate.We find that these dynamical tides are much weaker than gravitational radiation in driving the orbital decay of theWD–MBH binary, and they are also inefficient in changing theWD spin during the orbital evolution. Incorporating our computed tidal dissipation rate into a MESA-based WD evolution code, we find that tidal heating can lead to appreciable brightening of theWD and may induce runaway fusion in the hydrogen envelope well before the WD undergoes tidal disruption.

Asteroseismology of ZZ Ceti stars with fully evolutionary white dwarf models

ZZ Ceti stars are pulsating white dwarfs with a carbon-oxygen core build up during the core helium burning and thermally pulsing Asymptotic Giant Branch phases. Through the interpretation of their pulsation periods by means of asteroseismology, details about their origin and evolution can be inferred. The whole pulsation spectrum exhibited by ZZ Ceti stars strongly depends on the inner chemical structure. At present, there are several processes affecting the chemical profiles that are still not accurately determined. We present a study of the impact of the current uncertainties of the white dwarf formation and evolution on the expected pulsation properties of ZZ Ceti stars. Our analysis is based on a set of carbon-oxygen core white dwarf models with masses $0.548$ and $0.837 M_{\sun}$ derived from full evolutionary computations from the ZAMS to the ZZ Ceti domain. We have considered models in which we varied the number of thermal pulses, the amount of overshooting, and the $^{12}$C$(\alpha,\gamma)^{16}$O reaction rate within their uncertainties. We explore the impact of these major uncertainties in prior evolution on the chemical structure and the expected pulsation spectrum. We find that these uncertainties yield significant changes in the $g$-mode pulsation periods. We conclude that the uncertainties in the white dwarf progenitor evolution should be be taken into account in detailed asterseismological analysis of these pulsating stars.

The evolution of white dwarfs resulting from helium-enhanced, low-metallicity progenitor stars

Some globular clusters host multiple stellar populations with different chemical abundance patterns. This is particularly true for $\omega$ Centauri, which shows clear evidence of a helium- enriched sub-population characterized by a helium abundance as high as $Y= 0.4$. We present a whole and consistent set of evolutionary tracks from the ZAMS to the white dwarf stage appropriate for the study of the formation and evolution of white dwarfs resulting from the evolution of helium-rich progenitors. Different issues of the white dwarf evolution and their helium-rich progenitors have been explored. In particular, the final mass of the remnants, the role of overshooting during the thermally-pulsing phase, and the cooling of the resulting white dwarfs differ markedly from the evolutionary predictions of progenitor stars with standard initial helium abundance. Finally, the pulsational properties of the resulting white dwarfs are also explored. We find that, for the range of initial masses explored in this paper, the final mass of the helium-rich progenitors is markedly larger than the final mass expected from progenitors with the usual helium abundance. We also find that progenitors with initial mass smaller than $M_\star \simeq 0.65\,M_\odot$ evolve directly into helium-core white dwarfs in less than 14~Gyr, and that for larger progenitor masses the evolution of the resulting low-mass carbon-oxygen white dwarfs is dominated by residual nuclear burning. For helium-core white dwarfs, we find that they evolve markedly faster than their counterparts coming from standard progenitors. Also, in contrast with what occurs for white dwarfs resulting from progenitors with the standard helium abundance, the impact of residual burning on the cooling time of white dwarfs is not affected by the occurrence of overshooting during the thermally-pulsing phase of progenitor stars.

Asteroseismic test of rotational mixing in low-mass white dwarfs

We exploit the recent discovery of pulsations in mixed-atmosphere (He/H), extremely low-mass white dwarf precursors (ELM proto-WDs) to test the proposition that rotational mixing is a fundamental process in the formation and evolution of low-mass helium core white dwarfs. Rotational mixing has been shown to be a mechanism able to compete efficiently against gravitational settling, thus accounting naturally for the presence of He, as well as traces of metals such as Mg and Ca, typically found in the atmospheres of ELM proto-WDs. Here we investigate whether rotational mixing can maintain a sufficient amount of He in the deeper driving region of the star, such that it can fuel, through HeII-HeIII ionization, the observed pulsations in this type of stars. Using state-of-the-art evolutionary models computed with MESA, we show that rotational mixing can indeed explain qualitatively the very existence and general properties of the known pulsating, mixed-atmosphere ELM proto-WDs. Moreover, such objects are very likely to pulsate again during their final WD cooling phase.

First axion bounds from a pulsating helium-rich white dwarf star

The Peccei-Quinn mechanism proposed to solve the CP problem of Quantum Chromodynamics has as consequence the existence of axions, hypothetical weakly interacting particles whose mass is constrained to be on the sub-eV range. If these particles exist and interact with electrons, they would be emitted from the dense interior of white dwarfs, becoming an important energy sink for the star. Due to their well known physics, white dwarfs are good laboratories to study the properties of fundamental particles such as the axions. We study the general effect of axion emission on the evolution of helium-rich white dwarfs and on their pulsational properties. To this aim, we calculate evolutionary helium-rich white dwarf models with axion emission, and assess the pulsational properties of this models. Our results indicate that the rates of change of pulsation periods are significantly affected by the existence of axions. We are able for the first time to independently constrain the mass of the axion from the study of pulsating helium-rich white dwarfs. To do this, we use an estimation of the rate of change of period of the pulsating white dwarf PG 1351+489 corresponding to the dominant pulsation period. From an asteroseismological model of PG 1351+489 we obtain gae < 3.3 × 10−13 for the axion-electron coupling constant, or macos2β ≲ 11.5 meV for the axion mass. This constraint is relaxed to gae < 5.5 × 10−13 (macos2β ≲ 19.5 meV), when no detailed asteroseismological model is adopted for the comparison with observations.

White dwarf evolutionary sequences for low-metallicity progenitors: The impact of third dredge-up

We present new white dwarf evolutionary sequences for low-metallicity progenitors. White dwarf sequences have been derived from full evolutionary calculations that take into account the entire history of progenitor stars, including the thermally-pulsing and the post-asymptotic giant branch phases. We show that for progenitor metallicities in the range 0.00003--0.001, and in the absence of carbon enrichment due to the occurrence of a third dredge-up episode, the resulting H envelope of the low-mass white dwarfs is thick enough to make stable H burning the most important energy source even at low luminosities. This has a significant impact on white dwarf cooling times. This result is independent of the adopted mass-loss rate during the thermally-pulsing and post-AGB phases, and the planetary nebulae stage. We conclude that in the absence of third dredge-up episodes, a significant part of the evolution of low-mass white dwarfs resulting from low-metallicity progenitors is dominated by stable H burning. Our study opens the possibility of using the observed white dwarf luminosity function of low-metallicity globular clusters to constrain the efficiency of third dredge up episodes during the thermally-pulsing AGB phase of low-metallicity progenitors.

Discovery of ZZ Cetis in detached white dwarf plus main-sequence binaries

We present the first results of a dedicated search for pulsating white dwarfs (WDs) in detached white dwarf plus main-sequence binaries. Candidate systems were selected from a catalogue of WD+MS binaries, based on the surface gravities and effective temperatures of the WDs. We observed a total of 26 systems using ULTRACAM mounted on ESO's 3.5m New Technology Telescope (NTT) at La Silla. Our photometric observations reveal pulsations in seven WDs of our sample, including the first pulsating white dwarf with a main-sequence companion in a post common envelope binary, SDSSJ1136+0409. Asteroseismology of these new pulsating systems will provide crucial insight into how binary interactions, particularly the common envelope phase, affect the internal structure and evolution of WDs. In addition, our observations have revealed the partially eclipsing nature of one of our targets, SDSSJ1223-0056.

Globular cluster interstellar media: ionized and ejected by white dwarfs

UV radiation from white dwarfs can efficiently clear Galactic globular clusters (GCs) of their intra-cluster medium (ICM). This solves the problem of the missing ICM in clusters, which is otherwise expected to build up to easily observable quantities. To show this, we recreate the ionizing flux in 47 Tuc, following randomly generated stars through their AGB, post-AGB and white dwarf evolution. Each white dwarf can ionize all the material injected into the cluster by stellar winds for ~3 Myr of its evolution: ~40 such white dwarfs exist at any point. Every GC's ICM should be ionized. The neutral cloud in M15 should be caused by a temporary overdensity. A pressure-supported ICM will expand over the cluster's tidal radius, where it will be truncated, allowing Jeans escape. The modelled Jeans mass-loss rate approximates the total stellar mass-loss rate, allowing efficient clearing of ICM. Any cluster's ICM mass should equal the mass injected by its stars over the sound-travel time between the cluster core and tidal radius. We predict ~11.3 solar masses of ICM within 47 Tuc, cleared over ~4 Myr, compared to a dynamical timescale of 4.3 Myr. We present a new mass hierarchy, discussing the transition between globular clusters dwarf galaxies.

The timescale of low-mass proto-helium white dwarf evolution

A large number of low-mass (< 0.20 M_sun) helium white dwarfs (He WDs) have recently been discovered. The majority of these are orbiting another WD or a millisecond pulsar (MSP) in a close binary system; a few examples are found to show pulsations or to have a main-sequence star companion. There appear to be discrepancies between the current theoretical modelling of such low-mass He WDs and a number of key observed cases, indicating that their formation scenario remains to be fully understood. Here we investigate the formation of detached proto-He WDs in close-orbit low-mass X-ray binaries (LMXBs). Our prime focus is to examine the thermal evolution and the contraction phase towards the WD cooling track and investigate how this evolution depends on the WD mass. Our calculations are then compared to the most recent observational data. Numerical calculations with a detailed stellar evolution code were used to trace the mass-transfer phase in a large number of close-orbit LMXBs. Subsequently, we followed the evolution of the detached low-mass proto-He WDs, including stages with residual shell hydrogen burning and vigorous flashes caused by unstable CNO burning. We find that the time between Roche-lobe detachment until the low-mass proto-He WD reaches the WD cooling track is typically Delta_t_proto = 0.5 - 2 Gyr, depending systematically on the WD mass and therefore on its luminosity. The lowest WD mass for developing shell flashes is ~0.21 M_sun for progenitor stars of mass M2 <= 1.5 M_sun (and ~0.18 M_sun for M2 = 1.6 M_sun). The long timescale of low-mass proto-He WD evolution can explain a number of recent observations, including some MSP systems hosting He WD companions with very low surface gravities and high effective temperatures. We find no evidence for Delta_t_proto to depend on the occurrence of flashes and thus question the suggested dichotomy in thermal evolution of proto-WDs.

MID-INFRARED HIGH-CONTRAST IMAGING OF HD 114174 B: AN APPARENT AGE DISCREPANCY IN A “SIRIUS-LIKE” BINARY SYSTEM

We present new observations of the faint "Sirius-like" companion discovered to orbit HD 114174. Previous attempts to image HD 114174 B at mid-infrared wavelengths using NIRC2 at Keck have resulted in a non-detection. Our new L'-band observations taken with the Large Binocular Telescope and LMIRCam recover the companion ($\Delta L$ = 10.15 $\pm$ 0.15 mag, $\rho$ = 0.675'' $\pm$ 0.016'') with a high signal-to-noise ratio (10 $\sigma$). This measurement represents the deepest L' high-contrast imaging detection at sub-arcsecond separations to date, including extrasolar planets. We confirm that HD 114174 B has near-infrared colors consistent with the interpretation of a cool white dwarf ($J-L'$ = 0.76 $\pm$ 0.19 mag, $K-L'$ = 0.64 $\pm$ 0.20). New model fits to the object's spectral energy distribution indicate a temperature $T_{\rm eff}$ = 4260 $\pm$ 360 K, surface gravity log g = 7.94 $\pm$ 0.03, a cooling age t$_{c} \approx$ 7.8 Gyr, and mass $M$ = 0.54 $\pm$ 0.01 $M_{\odot}$. We find that the cooling age given by theoretical atmospheric models do not agree with the age of HD 114174 A derived from both isochronological and gyrochronological analyses. We speculate on possible scenarios to explain the apparent age discrepancy between the primary and secondary. HD 114174 B is a nearby benchmark white dwarf that will ultimately enable a dynamical mass estimate through continued Doppler and astrometric monitoring. Efforts to characterize its physical properties in detail will test theoretical atmospheric models and improve our understanding of white dwarf evolution, cooling, and progenitor masses.

Limits on the neutrino magnetic dipole moment from the luminosity function of hot white dwarfs

Recent determinations of the white dwarf luminosity function (WDLF) from very large surveys have extended our knowledge of the WDLF to very high luminosities. This, together with the availability of new full evolutionary white dwarf models that are reliable at high luminosities, have opened the possibility of testing particle emission in the core of very hot white dwarfs, where neutrino processes are dominant. We use the available WDLFs from the Sloan Digital Sky Survey and the SuperCOSMOS Sky Survey to constrain the value of the neutrino magnetic dipole moment ($\mu_\nu$). We constructed theoretical WDLFs for different values of $\mu_\nu$ and performed a $\chi^2$-test to derive constraints on the value of $\mu_\nu$. We also constructed a unified WDLF by averaging the SDSS and SSS and estimated the uncertainties by taking into account the differences between the WDLF at each magnitude bin. Then we compared all WDLFs with theoretical WDLFs.Comparison between theoretical WDLFs and both the SDSS and the averaged WDLF indicates that $\mu_\nu$ should be $\mu_\nu<5\times 10^{-12}\, e\hbar/(2m_e c)$. In particular, a $\chi^2$-test on the averaged WDLF suggests that observations of the disk WDLF exclude values of $\mu_\nu>5\times 10^{-12}e\hbar/(2m_e c)$ at more than a 95\% confidence level, even when conservative estimates of the uncertainties are adopted. Our study shows that modern WDLFs, which extend to the high-luminosity regime, are an excellent tool for constraining the emission of particles in the core of hot white dwarfs. However, discrepancies between different WDLFs suggest there might be some relevant unaccounted systematic errors. A larger set of completely independent WDLFs, as well as more detailed studies of the theoretical WDLFs and their own uncertainties, is desirable to explore the systematic uncertainties behind this constraint.

Asteroseismology of the DAV star HS 0507+0434B, including core-composition profiles

The DAV star, HS 0507+0434B, was observed by Fu et al. in 2007 and from 2009 December to 2010 January. There were a total of six triplets with nearly equal split, which were identified as l = 1 modes caused by rotation. In order to fit the six l = 1 modes, grids of white dwarf models are generated using the white dwarf evolution code (wdec). For the core-composition profiles, we choose the linear fittings to carbon profile of white dwarf models from modules for experiments in stellar astrophysics (mesa), which can be considered as the results of the real nuclear-burning process. Coupled with diffusion equilibrium H/He and He/C mixtures, we make grids of models in wdec and perform asteroseismology of HS 0507+0434B. There is a total error of 9.50 s for our best-fitting model, which is smaller than the result of Fu et al. (i.e. a total error of 22.1 s), when fitting the six l = 1 modes. The two other previously identified pulsation mode frequencies (286.1 and 743.40 s) can also be well fitted by our best-fitting model. The model parameters are Teff = 11 450 K, log g = 8.088, M* = 0.640 M⊙, log (MH/M*) = −6 and log (MHe/M*) = −3.

Diffusion and Coulomb Separation of Ions in Dense Matter

We analyze diffusion equations in strongly coupled Coulomb mixtures of ions in dense stellar matter. Strong coupling of ions in the presence of gravitational forces and electric fields (induced by plasma polarization in the presence of gravity) produces a specific diffusion current which can separate ions with the same A/Z (mass to charge number) ratios but different Z. This Coulomb separation of ions can be important for the evolution of white dwarfs and neutron stars.

White dwarf cooling sequences and cosmochronology

The evolution of white dwarfs is a simple gravothermal process. This means that their luminosity function, i.e. the number of white dwarfs per unit bolometric magnitude and unit volume as a function of bolometric magnitude, is a monotonically increasing function that decreases abruptly as a consequence of the finite age of the Galaxy. The precision and the accuracy of the white dwarf luminosity functions obtained with the recent large surveys together with the improved quality of the theoretical models of evolution of white dwarfs allow to feed the hope that in a near future it will be possible to reconstruct the history of the different Galactic populations.

Constraints on white dwarfs structure and evolution from asteroseismology

In the last two decades, considerable efforts have been devoted to the study of white dwarfs pulsators. Owing to ground-based multi-sites observational campaigns, and more recently to the long time-based high-precision photometric observations with the Kepler satellite, a large number of pulsating white dwarfs have been studied through the asteroseismology method. I emphasize here the accuracy that this method achieves on the determination of the total mass, one of the fundamental stellar parameter that can be derived from asteroseismology, among others. Then I discuss the issue of the age determination, emphasizing the contribution of the neutrinos cooling and the importance of the outer layer mass fraction. I show that some observations throw some doubts on the assumption that the cooling is the only cause of the observed frequency variations. Finally I summarize the results demonstrating that the entire stellar angular momentum is lost during the evolutionary phases prior to the white dwarf stage. While the asteroseismology of white dwarf stars has given a rich harvest of results on their internal structure and evolution, there are still some challenges to face. Determining accurate mass of the outer layers, specifically for the DA white dwarfs, is one of them. An increasing number of white dwarf pulsators are now known to exhibit frequency and amplitude temporal variations. The physical mechanism producing such variabilities has to be identified.