Enrique García-Berro

Evolutionary and pulsational properties of white dwarf stars

White dwarf stars are the final evolutionary stage of the vast majority of stars, including our Sun. Since the coolest white dwarfs are very old objects, the present population of white dwarfs contains a wealth of information on the evolution of stars from birth to death, and on the star formation rate throughout the history of our Galaxy. Thus, the study of white dwarfs has potential applications in different fields of astrophysics. In particular, white dwarfs can be used as independent reliable cosmic clocks, and can also provide valuable information about the fundamental parameters of a wide variety of stellar populations, such as our Galaxy and open and globular clusters. In addition, the high densities and temperatures characterizing white dwarfs allow these stars to be used as cosmic laboratories for studying physical processes under extreme conditions that cannot be achieved in terrestrial laboratories. Last but not least, since many white dwarf stars undergo pulsational instabilities, the study of their properties constitutes a powerful tool for applications beyond stellar astrophysics. In particular, white dwarfs can be used to constrain fundamental properties of elementary particles such as axions and neutrinos and to study problems related to the variation of fundamental constants. These potential applications of white dwarfs have led to renewed interest in the calculation of very detailed evolutionary and pulsational models for these stars. In this work, we review the essentials of the physics of white dwarf stars. We enumerate the reasons that make these stars excellent chronometers, and we describe why white dwarfs provide tools for a wide variety of applications. Special emphasis is placed on the physical processes that lead to the formation of white dwarfs as well as on the different energy sources and processes responsible for chemical abundance changes that occur along their evolution. Moreover, in the course of their lives, white dwarfs cross different pulsational instability strips. The existence of these instability strips provides astronomers with a unique opportunity to peer into their internal structure that would otherwise remain hidden from observers. We will show that this allows one to measure stellar masses with unprecedented precision and to infer their envelope thicknesses, to probe the core chemical stratification, and to detect rotation rates and magnetic fields. Consequently, in this work, we also review the pulsational properties of white dwarfs and the most recent applications of white dwarf asteroseismology.

The initial-final mass relationship of white dwarfs revisited: effect on the luminosity function and mass distribution

The initial-final mass relationship connects the mass of a white dwarf with the mass of its progenitor in the main-sequence. Although this function is of fundamental importance to several fields in modern astrophysics, it is not well constrained either from the theoretical or the observational points of view. In this work we revise the present semi-empirical initial-final mass relationship by re-evaluating the available data. The distribution obtained from grouping all our results presents a considerable dispersion, which is larger than the uncertainties. We have carried out a weighted least-squares linear fit of these data and a careful analysis to give some clues on the dependence of this relationship on some parameters such as metallicity or rotation. The semiempirical initial-final mass relationship arising from our study covers the range of initial masses from 1.0 to 6.5 M⊙, including in this way the low-mass domain, poorly studied until recently. Finally, we have also performed a test of the initial-final mass relationship by studying its effect on the luminosity function and on the mass distribution of white dwarfs. This was done by using different initial-final mass relationships from the literature, including the expression derived in this work, and comparing the results obtained with the observational data from the Palomar Green Survey and the Sloan Digital Sky Survey (SDSS). We find that the semi-empirical initial-final mass relationship derived here gives results in good agreement with the observational data, especially in the case of the white dwarf mass distribution.

The cooling of co white dwarfs: influence of the internal chemical distribution

White dwarfs are the remnants of stars of low and intermediate masses on the main sequence. Since they have exhausted all of their nuclear fuel, their evolution is just a gravothermal process. The release of energy only depends on the detailed internal structure and chemical composition and on the properties of the envelope equation of state and opacity; its consequences on the cooling curve (i.e., the luminosity vs. time relationship) depend on the luminosity at which this energy is released. The internal chemical profile depends on the rate of the 12C(α, γ)16O reaction as well as on the treatment of convection. High reaction rates produce white dwarfs with oxygen-rich cores surrounded by carbon-rich mantles. This reduces the available gravothermal energy and decreases the lifetime of white dwarfs. In this paper we compute detailed evolutionary models providing chemical profiles for white dwarfs having progenitors in the mass range from 1.0 to 7 M☉, and we examine the influence of such profiles in the cooling process. The influence of the process of separation of carbon and oxygen during crystallization is decreased as a consequence of the initial stratification, but it is still important and cannot be neglected. As an example, the best fit to the luminosity functions of Liebert et al. and Oswalt et al. gives an age of the disk of 9.3 and 11.0 Gyr, respectively, when this effect is taken into account, and only 8.3 and 10.0 Gyr when it is neglected.

New cooling sequences for old white dwarfs

We present full evolutionary calculations appropriate for the study of hydrogen-rich DA white dwarfs. This is done by evolving white dwarf progenitors from the zero-age main sequence, through the core hydrogen-burning phase, the helium-burning phase, and the thermally pulsing asymptotic giant branch phase to the white dwarf stage. Complete evolutionary sequences are computed for a wide range of stellarmasses and for two different metallicities, Z = 0.01, which is representative of the solar neighborhood, and Z = 0.001, which is appropriate for the study of old stellar systems, like globular clusters. During the white dwarf cooling stage, we self-consistently compute the phase in which nuclear reactions are still important, the diffusive evolution of the elements in the outer layers and, finally, we also take into account all the relevant energy sources in the deep interior of the white dwarf, such as the release of latent heat and the release of gravitational energy due to carbon–oxygen phase separation upon crystallization. We also provide colors and magnitudes for these sequences, based on a new set of improved non-gray white dwarf model atmospheres, which include the most up-to-date physical inputs like the Lyα quasi-molecular opacity. The calculations are extended down to an effective temperature of 2500 K. Our calculations provide a homogeneous set of evolutionary cooling tracks appropriate for mass and age determinations of old DA white dwarfs and for white dwarf cosmochronology of the different Galactic populations.

The formation and evolution of hydrogen-deficient post-AGB white dwarfs: The emerging chemical profile and the expectations for the PG 1159-DB-DQ evolutionary connection

We explore the formation and evolution of hydrogen-deficient post-AGB white dwarfs. To this end, we compute the complete evolution of an initially 2.7 Mstar from the zero-age main sequence through the thermally pulsing and mass-loss phases to the white dwarf stage. Particular attention is given to the chemical abundance changes during the whole evolution. A time-dependent scheme for the simultaneous treatment of abundance changes caused by nuclear reactions, diffusive overshoot- ing, salt fingers and convection is considered. We employed the double-diffusive mixing-length theory of convection for fluids with composition gradients. The study can therefore be considered as a test of its performance in low-mass stars. Also, time- dependent element diffusion for multicomponent gases is taken into account during the white dwarf evolution. The evolutionary stages corresponding to the last helium thermal pulse on the early white-dwarf cooling branch and the following born-again episode are carefully explored. Relevant aspects for PG 1159 stars and DB white dwarf evolution are studied in the framework of these new evolutionary models that take into account the history of the white dwarf progenitor. The scope of the calculations is extended to the domain of the helium-rich, carbon-contaminated DQ white dwarfs with the aim of exploring the plausibility of the evolutionary connection PG 1159-DB-DQ. In this regard, the implications for the double-layered chemical structure in pulsating DB white dwarfs is investigated. We examine the consequences of mass-loss episodes during the PG 1159 stage for the chemical stratification of the outer layer of DB and DQ white dwarfs.

A white dwarf cooling age of 8 Gyr for NGC 6791 from physical separation processes

NGC 6791 is a well studied open cluster that it is so close to us that can be imaged down to very faint luminosities. The main-sequence turn-off age (∼8 Gyr) and the age derived from the termination of the white dwarf cooling sequence (∼6 Gyr) are very different. One possible explanation is that as white dwarfs cool, one of the ashes of helium burning, 22Ne, sinks in the deep interior of these stars. At lower temperatures, white dwarfs are expected to crystallize and phase separation of the main constituents of the core of a typical white dwarf (12C and 16O) is expected to occur. This sequence of events is expected to introduce long delays in the cooling times, but has not hitherto been proven. Here we report that, as theoretically anticipated, physical separation processes occur in the cores of white dwarfs, resolving the age discrepancy for NGC 6791.

On the Evolution of Stars that Form Electron-degenerate Cores Processed by Carbon Burning. V. Shell Convection Sustained by Helium Burning, Transient Neon Burning, Dredge-out, Urca Cooling, and Other Properties of an 11 M? Population I Model Star

A stellar model of mass 11 M☉ and Population I composition is evolved from the hydrogen-burning main sequence through the core carbon-burning phase. In contrast with 9, 10, and 10.5 M☉ models studied in earlier papers of this series, carbon burning is ignited at the center of the 11 M☉ model. Like the 10.5 M☉ model, the 11 M☉ model experiences a dredge-out episode at the end of the carbon-burning phase. At the beginning of this episode, a semiconvective zone forms at the base of the hydrogen-rich envelope and carries hydrogen inward in mass toward the outer edge of a fully convective zone that is sustained by helium burning at its base. Hydrogen diffuses into the helium-rich convective zone untila hydrogen shell flash occurs. Helium burning dies out and the outer edge of the convective layer,sustained by fluxes due to hydrogen burning, extends outward in mass through hydrogen-rich material, mixing freshly synthesized nuclei outward. Then, hydrogen burning dies out and the outer edge of the convective shell, now sustained primarily by fluxes due to the release of gravothermal energy, moves outward until it reaches the inner edge of the convective envelope. Freshly synthesized material is then convected to the surface. Mixing during the final phase of homogenization in the convective envelope is maintained by fluxes due to the release of gravothermal energy. At the end of the dredge-out phase, the surface nitrogen abundance has decreased and the C/N ratio has changed from less than unity to larger than unity, showing that mixing has extended into regions where helium burning has manufactured substantial quantities of 12C and destroyed 14N. Prior to the dredge-out phase, neon burning is narrowly averted, and, after the dredge-out phase, neutrino losses due to electron capture and decay reactions between A=25 and A=23 isotopes in and above convective Urca shells cool the inner portions of the electron-degenerate oxygen-neon (ONe) core. Ultimately, the model becomes a thermally pulsing super-asymptotic giant branch (TPSAGB) star with an ONe core of mass ~1.368 M☉. Hydrogen and helium burning over a period of ~1.4×104 yr of TPSAGB evolution add a carbon-oxygen layer of mass ~0.014 M☉ to the electron-degenerate core. Then, electron captures on products of carbon burning lead to the collapse of the core into a neutron star and expulsion of the envelope in a weak Type II supernova explosion. The ratio of helium to hydrogen in the ejecta is approximately twice solar.

The Ages of Very Cool Hydrogen-rich White Dwarfs

The evolution of white dwarfs is essentially a cooling process that depends primarily on the energy stored in their degenerate cores and on the transparency of their envelopes. In this paper we compute accurate cooling sequences for carbon-oxygen white dwarfs with hydrogen dominated atmospheres for the full range of masses of interest. For this purpose we use the most accurate available physical inputs for both the equation of state and opacities of the envelope and for the thermodynamic quantities of the degenerate core. We also investigate the role of the latent heat in the computed cooling sequences. We present separately cooling sequences in which the effects of phase separation of the carbon-oxygen binary mixture upon crystallization have been neglected, and the delay introduced in the cooling times when this mechanism is properly taken into account, in order to compare our results with other published cooling sequences which do not include a treatment of this phenomenon. We find that the cooling ages of very cool white dwarfs with pure hydrogen atmospheres have been systematically underestimated by roughly 1.5 Gyr at log(L/L☉) = -4.5 for an otherwise typical ~0.6 M☉ white dwarf, when phase separation is neglected. If phase separation of the binary mixture is included, then the cooling ages are further increased by roughly 10%. Cooling tracks and cooling isochrones in several color-magnitude diagrams are presented as well.

Axions and the cooling of white dwarf stars

W hite dwarfsare the end productofthe lifes ofinterm ediate-and low-m assstarsand their evolution is described as a sim ple cooling process. Recently,it has been possible to determ ine with an unprecedented precision their lum inosity function,that is,the num ber ofstars per unit volum e and lum inosity interval. W e show here thatthe shape ofthe brightbranch ofthis function isonly sensitive to the averaged cooling rate ofwhite dwarfsand we propose to use thisproperty to check the possible existence ofaxions,a proposed but not yet detected weakly interacting particle. Our resultsindicatethattheinclusion oftheem ission ofaxionsin theevolutionarym odelsofwhitedwarfs noticeably im provesthe agreem entbetween the theoreticalcalculationsand the observationalwhite

High-resolution smoothed particle hydrodynamics simulations of the merger of binary white dwarfs

Context. The coalescence of two white dwarfs is the final outcome of a sizeable fraction of binary stellar systems. Moreover, this process has been proposed to explain several interesting astrophysical phenomena. Aims. We present the results of a set of high-resolution simulations of the merging process of two white dwarfs. Methods. We use an up-to-date smoothed particle hydrodynamics code that incorporates very detailed input physics and an improved treatment of the artificial viscosity. Our simulations have been done using a large number of particles (∼4 × 10 5 ) and covering the full range of masses and chemical compositions of the coalescing white dwarfs. We also compare the time evolution of the system during the first phases of the coalescence with what is obtained using a simplified treatment of mass transfer; we discuss in detail the characteristics of the final configuration; we assess the possible observational signatures of the merger, such as the associated gravitational waveforms and the fallback X-ray flares; and we study the long-term evolution of the coalescence. Results. The mass transfer rates obtained during the first phases of the merger episode agree with the theoretical expectations. In all the cases studied, the merged configuration is a central compact object surrounded by a self-gravitating Keplerian disk, except in the case where two equal-mass white dwarfs coalesce. Conclusions. We find that the overall evolution the system and the main characteristics of the of the final object agree with other previous studies in which lower resolutions were used. We also find that the fallback X-ray luminosities are close to 10 47 erg/s. The gravitational waveforms are characterized by the sudden disappearance of the signal in a few orbital periods.