We analyze the effect of the sedimentation of $^{22}$Ne on the local white dwarf luminosity function by studying scenarios under different Galactic metallicity models. We make use of an up-to-date population synthesis code based on Monte Carlo techniques to derive the synthetic luminosity function. Constant solar metallicity models are not able to simultaneously reproduce the peak and cut-off of the white dwarf luminosity function. The extra release of energy due to $^{22}$Ne sedimentation piles up more objects in brighter bins of the faint end of the luminosity function. The contribution of a single burst thick disk population increases the number of stars in the magnitude interval centered around $M_{\rm bol}=15.75$. Among the metallicity models studied, the one following a Twarog's profile is disposable. Our best fit model was obtained when a dispersion in metallicities around the solar metallicity value is considered along with a $^{22}$Ne sedimentation model, a thick disk contribution and an age of the thin disk of $8.8\pm0.2$ Gyr. Our population synthesis model is able to reproduce the local white dwarf luminosity function with a high degree of precision when a dispersion in metallicities around the solar value model is adopted. Although the effects of $^{22}$Ne sedimentation are only marginal and the contribution of a thick disk population is minor, both of them help in better fitting the peak and the cut-off regions of the white dwarf luminosity function.

Type Ia supernovae (SNe Ia) originate from the thermonuclear explosion of carbon-oxygen white dwarfs (CO WDs), giving rise to luminous optical transients. A relatively common variety of subluminous SNe Ia events, referred to as SNe Iax, are believed to arise from the failed detonation of a CO WD. In this paper, we explore failed detonation SNe Ia in the context of the double-degenerate channel of merging white dwarfs. In particular, we have carried out the first fully three-dimensional simulation of the merger of a ONe WD with a CO WD. While the hot, tidally-disrupted carbon-rich disk material originating from the CO WD secondary is readily susceptible to detonation in such a merger, the ONe WD primary core is not. This merger yields a failed detonation, resulting in the ejection of a small amount of mass, and leaving behind a kicked, super-Chandrasekhar ONe WD remnant enriched by the fallback of the products of nuclear burning. The resulting outburst is a rapidly-fading optical transient with a small amount of radioactive $^{56}$Ni powering the light curve. Consequently, the ONe-CO WD merger naturally produces a very faint and rapidly-fading transient, fainter even than the faintest Type Iax events observed to date, such as SN $2008$ha and SN 2010ae. More massive ONe primaries than considered here may produce brighter and longer-duration transients.

Since Dirac stated his Large Number Hypothesis the space-time variation of fundamental constants has been an active subject of research. Here we analyze the possible spatial variation of two fundamental constants: the fine structure constant {\alpha} and the speed of light c. We study the effects of such variations on the luminosity distance and on the peak luminosity of Type Ia supernovae (SNe Ia). For this, we consider the change of each fundamental constant separately and discuss a dipole model for its variation. Elaborating upon our previous work, we take into account the variation of the peak luminosity of Type Ia supernovae resulting from the variation of each of these fundamental constants. Furthermore, we also include the change of the energy release during the explosion, which was not studied before in the literature. We perform a statistical analysis to compare the predictions of the dipole model for {\alpha} and c variation with the Union2.1 and JLA compilations of SNe Ia. For this, we also allow the nuisance parameters of the distance estimator $\mu_0$ and the cosmological density matter $\Omega_m$ to vary. As a result of our analysis we obtain a first estimate of the possible spatial variation of the speed of light c. On the other hand, we find that there is no significant difference between the several phenomenological models studied here and the standard cosmological model, in which fundamental constants do not vary at all. Thus, we conclude that the actual set of data of Type Ia supernovae does not allow to verify the hypothetical spatial variation of fundamental constants.

Recent Hubble Space Telescope observations have unveiled the white dwarf cooling sequence of the Galactic bulge. Although the degenerate sequence can be well fitted employing the most up-to-date theoretical cooling sequences, observations show a systematic excess of red objects that cannot be explained by the theoretical models of single carbon-oxygen white dwarfs of the appropriate masses. Here we present a population synthesis study of the white dwarf cooling sequence of the Galactic bulge that takes into account the populations of both single white dwarfs and binary systems containing at least one white dwarf. These calculations incorporate state-of-the-art cooling sequences for white dwarfs with hydrogen-rich and hydrogen-deficient atmospheres, for both white dwarfs with carbon-oxygen and helium cores, and also take into account detailed prescriptions of the evolutionary history of binary systems. Our Monte Carlo simulator also incorporates all the known observational biases. This allows us to model with a high degree of realism the white dwarf population of the Galactic bulge. We find that the observed excess of red stars can be partially attributed to white dwarf plus main sequence binaries, and to cataclysmic variables or dwarf novae. Our best fit is obtained with a higher binary fraction and an initial mass function slope steeper than standard values, as well as with the inclusion of differential reddening and blending. Our results also show that the possible contribution of double degenerate systems or young and thick-disk bulge stars is negligible.

Increasing observational and theoretical evidence points to binary white dwarf mergers as the origin of some if not most normal Type Ia supernovae (SNe Ia). In this paper, we discuss the post-merger evolution of binary white dwarf (WD) mergers, and their relevance to the double-degenerate channel of SNe Ia. We present 3D simulations of carbon-oxygen (C/O) WD binary systems undergoing unstable mass transfer, varying both the total mass and the mass ratio. We demonstrate that these systems generally give rise to a one-armed gravitational spiral instability. The spiral density modes transport mass and angular momentum in the disk even in the absence of a magnetic field, and are most pronounced for secondary-to-primary mass ratios larger than $0.6$. We further analyze carbon burning in these systems to assess the possibility of detonation. Unlike the case of a $1.1 + 1.0 M_{\odot}$ C/O WD binary, we find that WD binary systems with lower mass and smaller mass ratios do not detonate as SNe Ia up to $\sim8-22$ outer dynamical times. Two additional models do however undergo net heating, and their secular increase in temperature could possibly result in a detonation on timescales longer than those considered here.

Recent advances in time-domain astronomy have uncovered a new class of optical transients with timescales shorter than typical supernovae and a wide range of peak luminosities. Several subtypes have been identified within this broad class, including Ca-rich transients, .Ia supernovae, and fast/bright transients. We examine the predic- tions from a state-of-the-art grid of three-dimensional simulations of dynamical white dwarf interactions in the context of these fast optical transients. We find that for colli- sions involving carbon-oxygen or oxygen-neon white dwarfs the peak luminosities and durations of the light curves in our models are in good agreement with the properties of fast/bright transients. When one of the colliding white dwarfs is made of helium the properties of the light curves are similar to those of Ca-rich gap transients. The model lightcurves from our white dwarf collisions are too slow to reproduce those of .Ia SNe, and too fast to match any normal or peculiar Type Ia supernova.

We present a set of full evolutionary sequences for white dwarfs with hydrogen-deficient atmospheres. We take into account the evolutionary history of the progenitor stars, all the relevant energy sources involved in the cooling, element diffusion in the very outer layers, and outer boundary conditions provided by new and detailed non-gray white dwarf model atmospheres for pure helium composition. These model atmospheres are based on the most up-to-date physical inputs. Our calculations extend down to very low effective temperatures, of $\sim 2\,500$~K, provide a homogeneous set of evolutionary cooling tracks that are appropriate for mass and age determinations of old hydrogen-deficient white dwarfs, and represent a clear improvement over previous efforts, which were computed using gray atmospheres.

Isolated magnetic white dwarfs have field strengths ranging from kilogauss to gigagauss. However, the origin of the magnetic field has not been hitherto elucidated. Whether these fields are fossil, hence the remnants of original weak magnetic fields amplified during the course of the evolution of their progenitor stars, or are the result of binary interactions or, finally, they are produced by other internal physical mechanisms during the cooling of the white dwarf itself, remains a mystery. At sufficiently low temperatures white dwarfs crystallize. Upon solidification, phase separation of its main constituents, 12C and 16O, and of the impurities left by previous evolution occurs. This process leads to the formation of a Rayleigh-Taylor unstable liquid mantle on top of a solid core. This convective region, as it occurs in Solar System planets like the Earth and Jupiter, can produce a dynamo able to yield magnetic fields of strengths of up to 0.1 MG, thus providing a mechanism that could explain magnetism in single white dwarfs.

The age-metallicity relation (AMR) is a fundamental observational constraint for understanding how the Galactic disc formed and evolved chemically in time. However, there is not yet an agreement on the observational properties of the AMR for the solar neighbourhood, primarily due to the difficulty in obtaining accurate stellar ages for individual field stars. We have started an observational campaign for providing the much needed observational input by using wide white dwarf-main sequence (WDMS) binaries. White dwarfs are natural clocks and can be used to derive accurate ages. Metallicities can be obtained from the main sequence companions. Since the progenitors of white dwarfs and the main sequence stars were born at the same time, WDMS binaries provide a unique opportunity to observationally constrain in a robust way the properties of the AMR. In this work we present the AMR derived from analysing a pilot sample of 23 WDMS binaries and provide clear observational evidence for the lack of correlation between age and metallicity at young and intermediate ages (0-7 Gyrs).

White dwarfs are the final remnants of low- and intermediate-mass stars. Their evolution is essentially a cooling process that lasts for $\sim 10$ Gyr. Their observed properties provide information about the history of the Galaxy, its dark matter content and a host of other interesting astrophysical problems. Examples of these include an independent determination of the past history of the local star formation rate, identification of the objects responsible for the reported microlensing events, constraints on the rate of change of the gravitational constant, and upper limits to the mass of weakly interacting massive particles. To carry on these tasks the essential observational tools are the luminosity and mass functions of white dwarfs, whereas the theoretical tools are the evolutionary sequences of white dwarf progenitors, and the corresponding white dwarf cooling sequences. In particular, the observed white dwarf luminosity function is the key manifestation of the white dwarf cooling theory, although other relevant ingredients are needed to compare theory and observations. In this review we summarize the recent attempts to empirically determine the white dwarf luminosity function for the different Galactic populations. We also discuss the biases that may affect its interpretation. Finally, we elaborate on the theoretical ingredients needed to model the white dwarf luminosity function, paying special attention to the remaining uncertainties, and we comment on some applications of the white dwarf cooling theory. Astrophysical problems for which white dwarf stars may provide useful leverage in the near future are also discussed.

We employ an asteroseismic model of L19-2, a relatively massive ($M_{\star} \sim 0.75 M_{\sun}$) and hot ($T_{\rm eff} \sim 12\,100$ K) pulsating DA (H-rich atmosphere) white dwarf star (DAV or ZZ Ceti variable), and use the observed values of the temporal rates of period change of its dominant pulsation modes ($\Pi \sim 113$ s and $\Pi \sim 192$ s), to derive a new constraint on the mass of the axion, the hypothetical non-barionic particle considered as a possible component of the dark matter of the Universe. If the asteroseismic model employed is an accurate representation of L19-2, then our results indicate hints of extra cooling in this star, compatible with emission of axions of mass $m_{\rm a} \cos^2 \beta \lesssim 25$ meV or an axion-electron coupling constant of $g_{\rm ae} \lesssim 7 \times 10^{-13}$.

A feasible design for a magnetic diagnostics subsystem for eLISA will be based on that of its precursor mission, LISA Pathfinder. Previous experience indicates that magnetic field estimation at the positions of the test masses has certain complications. This is due to two reasons. The first one is that magnetometers usually back-act due to their measurement principles (i.e., they also create their own magnetic fields), while the second is that the sensors selected for LISA Pathfinder have a large size, which conflicts with space resolution and with the possibility of having a sufficient number of them to properly map the magnetic field around the test masses. However, high-sensitivity and small-size sensors that significantly mitigate the two aforementioned limitations exist, and have been proposed to overcome these problems. Thus, these sensors will be likely selected for the magnetic diagnostics subsystem of eLISA. Here we perform a quantitative analysis of the new magnetic subsystem, as it is currently conceived, and assess the feasibility of selecting these sensors in the final configuration of the magnetic diagnostic subsystem.

The white dwarf luminosity function is an important tool to understand the properties of the Solar neighborhood, like its star formation history, and its age. Here we present a population synthesis study of the white dwarf population within 40~pc from the Sun, and compare the results of this study with the properties of the observed sample. We use a state-of-the-art population synthesis code based on Monte Carlo techniques, that incorporates the most recent and reliable white dwarf cooling sequences, an accurate description of the Galactic neighborhood, and a realistic treatment of all the known observational biases and selection procedures. We find a good agreement between our theoretical models and the observed data. In particular, our simulations reproduce a previously unexplained feature of the bright branch of the white dwarf luminosity function, which we argue is due to a recent episode of star formation. We also derive the age of the Solar neighborhood employing the position of the observed cut-off of the white dwarf luminosity function, obtaining ~8.9+-0.2 Gyr. We conclude that a detailed description of the ensemble properties of the population of white dwarfs within 40pc of the Sun allows us to obtain interesting constraints on the history of the Solar neighborhood.

Isolated magnetic white dwarfs have field strengths ranging from kilogauss to gigagauss, and constitute an interesting class of objects. The origin of the magnetic field is still the subject of a hot debate. Whether these fields are fossil, hence the remnants of original weak magnetic fields amplified during the course of the evolution of the progenitor of white dwarfs, or on the contrary, are the result of binary interactions or, finally, other physical mechanisms that could produce such large magnetic fields during the evolution of the white dwarf itself, remains to be elucidated. In this work we review the current status and paradigms of magnetic fields in white dwarfs, from both the theoretical and observational points of view.

NGC 6397 is one of the most interesting, well observed and theoretically studied globular clusters. The existing wealth of observations allows us to study the reliability of the theoretical white dwarf cooling sequences of low metallicity progenitors,to determine its age and the percentage of unresolved binaries, and to assess other important characteristics of the cluster, like the slope of the initial mass function, or the fraction of white dwarfs with hydrogen deficient atmospheres. We present a population synthesis study of the white dwarf population of NGC 6397. In particular, we study the shape of the color-magnitude diagram, and the corresponding magnitude and color distributions. We do this using an up-to-date Monte Carlo code that incorporates the most recent and reliable cooling sequences and an accurate modeling of the observational biases. We find a good agreement between our theoretical models and the observed data. In particular, we find that this agreement is best for those cooling sequences that take into account residual hydrogen burning. This result has important consequences for the evolution of progenitor stars during the thermally-pulsing asymptotic giant branch phase, since it implies that appreciable third dredge-up in low-mass, low-metallicity progenitors is not expected to occur. Using a standard burst duration of 1.0 Gyr, we obtain that the age of the cluster is 12.8+0.50-0.75 Gyr. Larger ages are also compatible with the observed data, but then realistic longer durations of the initial burst of star formation are needed to fit the luminosity function. We conclude that a correct modeling of the white dwarf opulation of globular clusters, used in combination with the number counts of main sequence stars provides an unique tool to model the properties of globular clusters.

White dwarfs are the fossils left by the evolution of low-and intermediate-mass stars, and have very long evolutionary timescales. This allows us to use them to explore the properties of old populations, like the Galactic halo. We present a population synthesis study of the luminosity function of halo white dwarfs, aimed at investigating which information can be derived from the currently available observed data. We employ an up-to-date population synthesis code based on Monte Carlo techniques, that incorporates the most recent and reliable cooling sequences for metal poor progenitors as well as an accurate modeling of the observational biases. We find that because the observed sample of halo white dwarfs is restricted to the brightest stars only the hot branch of the white dwarf luminosity function can be used for such purposes, and that its shape function is almost insensitive to the most relevant inputs, like the adopted cooling sequences, the initial mass function, the density profile of the stellar spheroid, or the adopted fraction of unresolved binaries. Moreover, since the cut-off of the observed luminosity has not been yet determined only lower limits to the age of the halo population can be placed. We conclude that the current observed sample of the halo white dwarf population is still too small to obtain definite conclusions about the properties of the stellar halo, and the recently computed white dwarf cooling sequences which incorporate residual hydrogen burning should be assessed using metal-poor globular clusters.

The core-degenerate (CD) scenario for type Ia supernovae (SN Ia) involves the merger of the hot core of an asymptotic giant branch (AGB) star and a white dwarf, and might contribute a non-negligible fraction of all thermonuclear supernovae. Despite its potential interest, very few studies, and based on only crude simplifications, have been devoted to investigate this possible scenario, compared with the large efforts invested to study some other scenarios. Here we perform the first three-dimensional simulations of the merger phase, and find that this process can lead to the formation of a massive white dwarf, as required by this scenario. We consider two situations, according to the mass of the circumbinary disk formed around the system during the final stages of the common envelope phase. If the disk is massive enough, the stars merge on a highly eccentric orbit. Otherwise, the merger occurs after the circumbinary disk has been ejected and gravitational wave radiation has brought the stars close to the Roche lobe radius on a nearly circular orbit. Not surprisingly, the overall characteristics of the merger remnants are similar to those found for the double-degenerate (DD) scenario, independently of the very different core temperature and of the orbits of the merging stars. They consist of a central massive white dwarf, surrounded by a hot, rapidly rotating corona and a thick debris region.

We conduct numerical simulations of the interacting ejecta from an exploding CO white dwarf (WD) with a He WD donor in the double-detonation scenario for Type Ia supernovae (SNe Ia), and study the possibility of exploding the companion WD. We also study the long time imprint of the collision on the supernova remnant. When the donor He WD has a low mass, M_WD = 0.2Msun, it is at a distance of ~ 0.08Rsun from the explosion, and helium is not ignited. The low mass He WD casts an 'ejecta shadow' behind it. By evolving the ejecta for longer times, we find that the outer parts of the shadowed side are fainter and its boundary with the ambient gas is somewhat flat. More massive He WD donors, M_WD ~ 0.4Msun, must be closer to the CO WD to transfer mass. At a distance of a < 0.045Rsun helium is detonated and the He WD explodes, leading to a triple detonation scenario. In the explosion of the donor WD approximately 0.15Msun of unburned helium is ejected. This might be observed as a peculiar type Ib supernova.

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.

Thermonuclear, or Type Ia supernovae (SNe Ia), originate from the explosion of carbon--oxygen white dwarfs, and serve as standardizable cosmological candles. However, despite their importance, the nature of the progenitor systems that give rise to SNe Ia has not been hitherto elucidated. Observational evidence favors the double-degenerate channel in which merging white dwarf binaries lead to SNe Ia. Furthermore, significant discrepancies exist between observations and theory, and to date, there has been no self-consistent merger model that yields a SNe Ia. Here we show that a spiral mode instability in the accretion disk formed during a binary white dwarf merger leads to a detonation on a dynamical timescale. This mechanism sheds light on how white dwarf mergers may frequently yield SNe Ia.

The large-number hypothesis conjectures that fundamental constants may vary. Accordingly, the spacetime variation of fundamental constants has been an active subject of research for decades. Recently, using data obtained with large telescopes a phenomenological model in which the fine structure constant might vary spatially has been proposed. We test whether this hypothetical spatial variation of {\alpha}, which follows a dipole law, is compatible with the data of distant thermonuclear supernovae. Unlike previous works, in our calculations we consider not only the variation of the luminosity distance when a varying {\alpha} is adopted, but we also take into account the variation of the peak luminosity of Type Ia supernovae resulting from a variation of {\alpha}. This is done using an empirical relation for the peak bolometric magnitude of thermonuclear supernovae that correctly reproduces the results of detailed numerical simulations. We find that there is no significant difference between the several phenomenological models studied here and the standard one, in which {\alpha} does not vary spatially. We conclude that the present set of data of Type Ia supernovae is not able to distinguish the standard model from the dipole models, and thus cannot be used to discard nor to confirm the proposed spatial variation of {\alpha}.

We follow the mass expelled during the WD-WD merger process in a particular case of the Double-Degenerate (DD) scenario for Type Ia supernovae (SNe Ia), and find that the interaction of the SN ejecta with the resulting wind affects the early (first day) light curve in a way that may be in conflict with some SN Ia observations, if the detonation occurs shortly after the merger (i.e., $10^3~{\rm sec} \lesssim t_{\rm exp} \lesssim 1~{\rm day}$). The main source of the expelled mass is a disk-wind, or jets that are launched by the accretion disk around the more massive WD during the viscous phase of the merger. This disk-originated matter (DOM) will be shocked and heated by the SN ejecta from an explosion, leading to additional radiation in the early lightcurve. This enhanced early radiation could then be interpreted as an explosion originating from a progenitor having an inferred radius of one solar radius or more, in conflict with observations of SN 2011fe.

47 Tucanae is one of the most interesting and well observed and theoretically studied globular clusters. This allows us to study the reliability of our understanding of white dwarf cooling sequences, to confront different methods to determine its age, and to assess other important characteristics, like its star formation history. Here we present a population synthesis study of the cooling sequence of the globular cluster 47 Tucanae. In particular, we study the distribution of effective temperatures, the shape of the color-magnitude diagram, and the corresponding magnitude and color distributions. We do so using an up-to-date population synthesis code based on Monte Carlo techniques, that incorporates the most recent and reliable cooling sequences and an accurate modeling of the observational biases. We find a good agreement between our theoretical models and the observed data. Thus, our study, rules out previous claims that there are still missing physics in the white dwarf cooling models at moderately high effective temperatures. We also derive the age of the cluster using the termination of the cooling sequence, obtaining a good agreement with the age determinations using the main-sequence turn-off. Finally, we find that the star formation history of the cluster is compatible with that btained using main sequence stars, which predict the existence of two distinct populations. We conclude that a correct modeling of the white dwarf population of globular clusters, used in combination with the number counts of main sequence stars provides an unique tool to model the properties of globular clusters.

Pulsating white dwarf stars can be used as astrophysical laboratories to constrain the properties of weakly interacting particles. Comparing the cooling rates of these stars with the expected values from theoretical models allows us to search for additional sources of cooling due to the emission of axions, neutralinos, or neutrinos with magnetic dipole moment. In this work, we derive an upper bound to the neutrino magnetic dipole moment using an estimate of the rate of period change of the pulsating DB white dwarf star PG 1351+489. By comparing the theoretical rate of change of period expected for this star with the rate of change of period with time of PG 1351+489, we assess the possible existence of additional cooling by neutrinos with magnetic dipole moment. Our models suggest the existence of some additional cooling in this pulsating DB white dwarf, consistent with a non-zero magnetic dipole moment. Our upper limit for the neutrino magnetic dipole moment is somewhat less restrictive than, but still compatible with, other limits inferred from the white dwarf luminosity function or from the color-magnitude diagram of the Globular cluster M5. Further improvements of the measurement of the rate of period change of the dominant pulsation mode of PG 1351+489 will be necessary to confirm our bound.

Detached WD+MS PCEBs are perhaps the most suitable objects for testing predictions of close-compact binary-star evolution theories, in particular, CE evolution. The population of WD+MS PCEBs has been simulated by several authors in the past and compared with observations. However, most of those predictions did not take the possible contributions to the envelope ejection from additional sources of energy (mostly recombination energy) into account. Here we update existing binary population models of WD+MS PCEBs by assuming that a fraction of the recombination energy available within the envelope contributes to ejecting the envelope. We performed Monte Carlo simulations of 10^7 MS+MS binaries for 9 different models using standard assumptions for the initial primary mass function, binary separations, and initial-mass-ratio distribution and evolved these systems using the publicly available BSE code. Including a fraction of recombination energy leads to a clear prediction of a large number of long orbital period (>~10 days) systems mostly containing high-mass WDs. The fraction of systems with He-core WD primaries increases with the CE efficiency and the existence of very low-mass He WDs is only predicted for high values of the CE efficiency (>~0.5). All models predict on average longer orbital periods for PCEBs containing C/O-core WDs than for PCEBs containing He WDs. This effect increases with increasing values of both efficiencies. Longer periods after the CE phase are also predicted for systems containing more massive secondary stars. The initial-mass-ratio distribution affects the distribution of orbital periods, especially the distribution of secondary star masses. Our simulations, in combination with a large and homogeneous observational sample, can provide constraints on the values of the CE efficiencies, as well as on the initial-mass-ratio distribution for MS+MS binary stars.

We compute the possible observational signatures of white dwarf dynamical interactions in dense stellar environments. Specifically, we compute the emission of gravitational waves, and we compare it with the sensitivity curves of planned space-borne gravitational wave detectors. We also compute the light curves for those interactions in which a detonation occurs, and one of the stars is destroyed, as well as the corresponding neutrino luminosities. We find that for the three possible outcomes of these interactions - which are the formation of an eccentric binary system, a lateral collision in which several mass transfer episodes occur, and a direct one in which just a single mass transfer episode takes place - only those in which an eccentric binary are formed are likely to be detected by the planned gravitational wave mission eLISA, while more sensitive detectors would be able to detect the signals emitted in lateral collisions. On the other hand, the light curves (and the thermal neutrino emission) of these interactions are considerably different, producing both very powerful outbursts and low luminosity events. Finally, we also calculate the X-ray signature produced in the aftermath of those interactions for which a merger occurs. We find that the temporal evolution follows a power law with the same exponent found in the case of the mergers of two neutron stars, although the total energy released is smaller.

It has been known for a long time that stellar metallicity plays a significant role in the determination of the ages of the different Galactic stellar populations, when main sequence evolutionary tracks are employed. Here we analyze the role that metallicity plays on the white dwarf luminosity function of the Galactic disk, which is often used to determine its age. We employ a Monte Carlo population synthesis code that accounts for the properties of the population of Galactic disk white dwarfs. Our code incorporates the most up-to-date evolutionary cooling sequences for white dwarfs with hydrogen-rich and hydrogen-deficient atmospheres for both carbon-oxygen and oxygen-neon cores. We use two different models to assess the evolution of the metallicity, one in which the adopted metallicity is constant with time, but with a moderate dispersion, and a second one in which the metallicity increases with time. We found that our theoretical results are in a very satisfactory agreement with the observational luminosity functions obtained from the Sloan Digital Sky Survey (SDSS) and from the SuperCOSMOS Sky Survey (SSS), independently of the adopted age-metallicity law. In particular, we found that the age-metallicity law has no noticeable impact in shaping the bright branch of the white dwarf luminosity function, and that the position of its cut-off is almost insensitive to the adopoted age-metallicity relationship. Because the shape of the bright branch of the white dwarf luminosity function is insensitive to the age-metallicity law, it can be safely employed to test the theoretical evolutionary sequences, while due to the limited sensitivity of the position of the drop-off to the distribution of metallicities, its location provides a robust indicator of the age of the Galactic disk.

Detached white dwarf + main sequence (WD+MS) systems represent the simplest population of post-common envelope binaries (PCEBs). Since the ensemble properties of this population carries important information about the characteristics of the common-envelope (CE) phase, it deserves close scrutiny. However, most population synthesis studies do not fully take into account the effects of the observational selection biases of the samples used to compare with the theoretical simulations. Here we present the results of a set of detailed Monte Carlo simulations of the population of WD+MS binaries in the Sloan Digital Sky Survey (SDSS) Data Release 7. We used up-to-date stellar evolutionary models, a complete treatment of the Roche lobe overflow episode, and a full implementation of the orbital evolution of the binary systems. Moreover, in our treatment we took into account the selection criteria and all the known observational biases. Our population synthesis study allowed us to make a meaningful comparison with the available observational data. In particular, we examined the CE efficiency, the possible contribution of internal energy, and the initial mass ratio distribution (IMRD) of the binary systems. We found that our simulations correctly reproduce the properties of the observed distribution of WD+MS PCEBs. In particular, we found that once the observational biases are carefully taken into account, the distribution of orbital periods and of masses of the WD and MS stars can be correctly reproduced for several choices of the free parameters and different IMRDs, although models in which a moderate fraction (<=10%) of the internal energy is used to eject the CE and in which a low value of CE efficiency is used (<=0.3) seem to fit better the observational data. We also found that systems with He-core WDs are over-represented in the observed sample, due to selection effects.

We present a coherent and detailed Monte Carlo simulation of the population of hot white dwarfs. We assess the statistical significance of the hot end of the white dwarf luminosity function and the role played by the bolometric corrections of hydrogen-rich white dwarfs at high effective temperatures. We use the most up-to-date stellar evolutionary models and implement a full description of the observational selection biases to obtain realistic simulations of the observed white dwarf population. Our theoretical results are compared with the luminosity function of hot white dwarfs obtained from the Sloan Digital Sky Survey (SDSS), for both DA and non-DA white dwarfs. We find that the theoretical results are in excellent agreement with the observational data for the population of white dwarfs with hydrogen deficient atmospheres (non-DA white dwarfs). For the population of white dwarfs with hydrogen-rich atmospheres (white dwarfs of the DA class), our simulations show some discrepancies with the observations for the brightest luminosity bins. These discrepancies can be attributed to the way in which the masses of the white dwarfs contributing to this luminosity bin have been computed, as most of them have masses smaller than the theoretical lower limit for carbon-oxygen white dwarfs. We conclude that the way in which the observational luminosity function of hot white dwarfs is obtained is very sensitive to the particular implementation of the method used to derive the masses of the sample. We also provide a revised luminosity function for hot white dwarfs with hydrogen-rich atmospheres.

A secular variation of $G$ modifies the structure and evolutionary time scales of white dwarfs. Using an state-of-the-art stellar evolutionary code, an up-to-date pulsational code, and a detailed population synthesis code we demonstrate that the effects of a running $G$ are obvious both in the properties of individual white dwarfs, and in those of the white dwarf populations in clusters. Specifically, we show that the white dwarf evolutionary sequences depend on both the value of $\dot G/G$, and on the value of $G$ when the white dwarf was born. We show as well that the pulsational properties of variable white dwarfs can be used to constrain $\dot G/G$. Finally, we also show that the ensemble properties of of white dwarfs in clusters can also be used to set upper bounds to $\dot G/G$. Precisely, the tightest bound --- $\dot G/G \sim -1.8 10^{-12}$ yr$^{-1}$ --- is obtained studying the population of the old, metal-rich, well populated, open cluster NGC 6791. Less stringent upper limits can be obtained comparing the theoretical results obtained taking into account the effects of a running $G$ with the measured rates of change of the periods of two well studied pulsating white dwarfs, G117--B15A and R548. Using these white dwarfs we obtain $\dot G/G\sim -1.8\times 10^{-10}$ yr$^{-1}$, and $\dot G/G\sim -1.3\times 10^{-10}$ yr$^{-1}$, respectively, which although less restrictive than the previous bound, can be improved measuring the rate of change of the period of massive white dwarfs.

In old, dense stellar systems collisions of white dwarfs are a rather frequent phenomenon. Here we present the results of a comprehensive set of Smoothed Particle Hydrodynamics simulations of close encounters of white dwarfs aimed to explore the outcome of the interaction and the nature of the final remnants for different initial conditions. Depending on the initial conditions and the white dwarf masses, three different outcomes are possible. Specifically, the outcome of the interaction can be either a direct or a lateral collision or the interaction can result in the formation of an eccentric binary system. In those cases in which a collision occurs, the infalling material is compressed and heated such that the physical conditions for a detonation may be reached during the most violent phases of the merger. While we find that detonations occur in a significant number of our simulations, in some of them the temperature increase in the shocked region rapidly lifts degeneracy, leading to the quenching of the burning. We thus characterize under which circumstances a detonation is likely to occur as a result of the impact of the disrupted star on the surface of the more massive white dwarf. Finally, we also study which interactions result in bound systems, and in which ones the more massive white dwarf is also disrupted as a consequence of the dynamical interaction. The sizable number of simulations performed in this work allows to find how the outcome of the interaction depends on the distance at closest approach, and on the masses of the colliding white dwarfs, and which is the chemical pattern of the nuclearly processed material. Finally, we also discuss the influence of the masses and core chemical compositions of the interacting white dwarfs and the different kinds of impact in the properties of the remnants.

An accurate assessment of white dwarf cooling times is paramount to place white dwarf cosmochronology of Galactic populations on more solid grounds. This issue is particularly relevant in view of the enhanced observational capabilities provided by the next generation of Extremely Large Telescopes, that will offer more avenues to employ white dwarfs as probes of Galactic evolution and test-beds of fundamental physics. We estimate for the first time the consistency of results obtained from two independent and widely used evolutionary codes (BaSTI and LPCODE) for white dwarf models with fixed mass and chemical stratification, when the same input physics is employed in both codes. We considered 0.55Msun white dwarf models with both pure carbon and uniform carbon-oxygen (50/50 mass fractions) core. We have assessed for the first time the maximum possible accuracy in the current estimates of white dwarf cooling times, resulting only from the different implementations of the stellar evolution equations and homogeneous input physics in two independent stellar evolution codes. This accuracy amounts to 2% at luminosities lower than log(L/L_{\sun})= -1.5$. This agreement is true for both pure carbon and uniform carbon-oxygen stratification core models, and also when the release of latent heat and carbon-oxygen phase separation are considered. This difference is smaller than the uncertainties in cooling times due to the present uncertainties in the white dwarf chemical stratification. We have also determined quantitatively and explained the effect of varying equation of state, low-temperature radiative opacities and electron conduction opacities in our calculation. Finally, we extend the scope of our work by providing tabulations of our cooling sequences and the required input physics, that can be used as a comparison test of cooling times obtained from other white dwarf evolutionary codes.

A secular variation of the gravitational constant modifies the structure and evolutionary time scales of white dwarfs. Using an state-of-the-art stellar evolutionary code and an up-to-date pulsational code we compute the effects of a secularly varying $G$ on the pulsational properties of variable white dwarfs. Comparing the the theoretical results obtained taking into account the effects of a running $G$ with the observed periods and measured rates of change of the periods of two well studied pulsating white dwarfs, G117--B15A and R548, we place constraints on the rate of variation of Newton's constant. We derive an upper bound $\dot G/G\sim -1.8\times 10^{-10}$ yr$^{-1}$ using the variable white dwarf G117--B15A, and $\dot G/G\sim -1.3\times 10^{-10}$ yr$^{-1}$ using R548. Although these upper limits are currently less restrictive than those obtained using other techniques, they can be improved in a future measuring the rate of change of the period of massive white dwarfs.

It has long been accepted that a possible mechanism for explaining the existence of magnetic white dwarfs is the merger of a binary white dwarf system, as there are viable mechanisms for producing sustainable magnetic fields within the merger product. However, the lack of rapid rotators in the magnetic white dwarf population has been always considered a problematic issue of this scenario. Smoothed Particle Hydrodynamics simulations show that in mergers in which the two white dwarfs have different masses a disc around the central compact object is formed. If the central object is magnetized it can interact with the disc through its magnetosphere. The torque applied by the disc changes the spin of the star, whereas the transferred angular momentum from the star to the disc determines the properties of the disc. In this work we build a model for the disc evolution under the effect of magnetic accretion, and for the angular momentum evolution of the star, which can be compared with the observations. Our model predicts that the magnetospheric interaction of magnetic white dwarfs with their discs results in a significant spin down, and we show that for magnetic white dwarfs with relatively strong fields (larger than 10 MG) the observed rotation periods of the magnetic white dwarf population can be reproduced. We also investigate whether turbulence can be sustained during the late phases of the evolution of the system. When a critical temperature below which turbulence is not sustained is introduced into the model, the periods of the three fast rotating, strongly magnetic, massive white dwarfs in the solar neighborhood are recovered.

It has long been accepted that a possible mechanism for explaining the existence of magnetic white dwarfs is the merger of a binary white dwarf system, as there are viable mechanisms for producing sustainable magnetism within the merger product. However, the lack of rapid rotators in the magnetic white dwarf population has been always considered a problematic issue of this scenario. In order to explain this discrepancy we build a model in which the interaction between the magnetosphere of the star and the disk induces angular momentum transfer. Our model predicts that the magnetospheric interaction of magnetic white dwarfs with their disks results in a significant spin down, and we show that the observed rotation period of REJ 0317-853, which is suggested to be a product of a double degenerate merger, can be reproduced.

We present a detailed Monte Carlo simulation of the population of the hot branch of the white dwarf luminosity function. We used the most up-to-date stellar evolutionary models and we implemented a full description of the observational selection biases. Our theoretical results are compared with the luminosity function of hot white dwarfs obtained from the Sloan Digital Sky Survey (SDSS), for both DA and non-DA white dwarfs. For non-DA white dwarfs we find an excellent agreement with the observational data, while for DA white dwarfs our simulations show some discrepancies with the observations for the brightest luminosity bins, those corresponding to L>= 10 L_sun.

High-field magnetic white dwarfs have been long suspected to be the result of stellar mergers. However, the nature of the coalescing stars and the precise mechanism that produces the magnetic field are still unknown. Here we show that the hot, convective, differentially rotating corona present in the outer layers of the remnant of the merger of two degenerate cores is able to produce magnetic fields of the required strength that do not decay for long timescales. We also show, using an state-of-the-art Monte Carlo simulator, that the expected number of high-field magnetic white dwarfs produced in this way is consistent with that found in the solar neighborhood.

White dwarf evolution is essentially a gravothermal cooling process, which,for cool white dwarfs, sensitively depends on the treatment of the outer boundary conditions. We provide detailed outer boundary conditions appropriate for computing the evolution of cool white dwarfs employing detailed non-gray model atmospheres for pure H composition. We also explore the impact on the white dwarf cooling times of different assumptions for energy transfer in the atmosphere of cool white dwarfs. Detailed non-gray model atmospheres are computed taken into account non-ideal effects in the gas equation of state and chemical equilibrium, collision-induced absorption from molecules, and the Lyman alpha quasi-molecular opacity. Our results show that the use of detailed outer boundary conditions becomes relevant for effective temperatures lower than 5800 and 6100K for sequences with 0.60 and 0.90 M_sun, respectively. Detailed model atmospheres predict ages that are up to approx 10% shorter at log L/L_sun=-4 when compared with the ages derived using Eddington-like approximations at tau_Ross=2/3. We also analyze the effects of various assumptions and physical processes of relevance in the calculation of outer boundary conditions. In particular, we find that the Ly_alpha red wing absorption does not affect substantially the evolution of white dwarfs. White dwarf cooling timescales are sensitive to the surface boundary conditions for T_eff < 6000K. Interestingly enough, non-gray effects have little consequences on these cooling times at observable luminosities. In fact, collision-induced absorption processes, which significantly affect the spectra and colors of old white dwarfs with hydrogen-rich atmospheres, have not noticeable effects in their cooling rates, except throughout the Rosseland mean opacity.

The Large Observatory for X-ray Timing (LOFT) is one of the four candidate ESA M3 missions considered for launch in the time-frame of 2022. It is specifically designed to perform fast X-ray timing and probe the status of the matter near black holes and neutron stars. The LOFT scientific payload consists of a Large Area Detector and a Wide Field Monitor. The LAD is a 10 m^2-class pointed instrument with high spectral (200 eV @ 6 keV) and timing (< 10 {\mu}s) resolution over the 2-80 keV range. It is designed to observe persistent and transient X-ray sources with a very large dynamic range from a few mCrab up to an intensity of 15 Crab. An unprecedented large throughput (~280.000 cts/s from the Crab) is achieved with a segmented detector, making pile-up and dead-time, often worrying or limiting focused experiments, secondary issues. We present the on-board data handling concept that follows the highly segmented and hierarchical structure of the instrument from the front-end electronics to the on-board software. The system features customizable observation modes ranging from event-by-event data for sources below 0.5 Crab to individually adjustable time resolved spectra for the brighter sources. On-board lossless data compression will be applied before transmitting the data to ground.

We employ a state-of-the-art asteroseismological model of G117-B15A, the archetype of the H-rich atmosphere (DA) white dwarf pulsators (also known as DAV or ZZ Ceti variables), and use the most recently measured value of the rate of period change for the dominant mode of this pulsating star to derive a new constraint on the mass of axion, the still conjectural non-barionic particle considered as candidate for dark matter of the Universe. Assuming that G117-B15A is truly represented by our asteroseismological model, and in particular, that the period of the dominant mode is associated to a pulsation g-mode trapped in the H envelope, we find strong indications of the existence of extra cooling in this star, compatible with emission of axions of mass m_a \cos^2 \beta = 17.4^{+2.3}_{-2.7} meV.

The origin of R Coronae Borealis stars has been elusive for over 200 years. Currently, two theories for their formation have been presented. These are the Final Flash scenario, in which a dying asymptotic giant branch (AGB) star throws off its atmosphere to reveal the hydrogen poor, heavily processed material underneath, and the double degenerate scenario, in which two white dwarfs merge to produce a new star with renewed vigour. Some theories predict that the temperatures reached during the latter scenario would destroy any lithium originally present in the white dwarfs. The observed lithium content of some R Coronae Borealis stars, therefore, is often interpreted as an indication that the Final Flash scenario best describes their formation. In this paper, it is shown that lithium production can, indeed, occur in the merging of a helium white dwarf with a carbon-oxygen white dwarf if their chemical composition, particularly that of 3He, is fully considered. The production mechanism is described in detail, and the sensitivity of lithium production to the merger environment is investigated. Nucleosynthesis post-processing calculations of smoothed-particle hydrodynamics (SPH) tracer particles are performed to show that any lithium produced in these environments will be concentrated towards the cloud of material surrounding the R CrB star. Measurements of the lithium content of these clouds would, therefore, provide a valuable insight into the formation mechanism of R CrB stars.

High-field magnetic white dwarfs have been long suspected to be the result of stellar mergers. However, the nature of the coalescing stars and the precise mechanism that produces the magnetic field are still unknown. Here we show that the hot, convective, differentially rotating corona present in the outer layers of the remnant of the merger of two degenerate cores is able to produce magnetic fields of the required strength that do not decay for long timescales. We also show, using an state-of-the-art Monte Carlo simulator, that the expected number of high-field magnetic white dwarfs produced in this way is consistent with that found in the Solar neighborhood.

Type Ia supernovae are thought to be the outcome of the thermonuclear explosion of a white dwarf in a close binary system. Two possible scenarios, not necessarily incompatible, have been advanced. One assumes a white dwarf that accretes matter from a nondegenerate companion (the single degenerate scenario), the other assumes two white dwarfs that merge as a consequence of the emission of gravitational waves (the double degenerate scenario). The delay time distribution of star formation bursts strongly suggests that the DD scenario should be responsible of the late time explosions, but this contradicts the common wisdom that the outcome of the merging of two white dwarfs is an accretion induced collapse to a neutron star. In this contribution we review some of the most controversial issues of this problem.

Cool white dwarfs are reliable and independent stellar chronometers. The most common white dwarfs have carbon-oxygen dense cores. Consequently, the cooling ages of very cool white dwarfs sensitively depend on the adopted phase diagram of the carbon-oxygen binary mixture. A new phase diagram of dense carbon-oxygen mixtures appropriate for white dwarf interiors has been recently obtained using direct molecular dynamics simulations. In this paper, we explore the consequences of this phase diagram in the evolution of cool white dwarfs. To do this we employ a detailed stellar evolutionary code and accurate initial white dwarf configurations, derived from the full evolution of progenitor stars. We use two different phase diagrams, that of Horowitz et al. (2010), which presents an azeotrope, and the phase diagram of Segretain & Chabrier (1993), which is of the spindle form. We computed the evolution of 0.593 and 0.878M_sun white dwarf models during the crystallization phase, and we found that the energy released by carbon-oxygen phase separation is smaller when the new phase diagram of Horowitz et al. (2010) is used. This translates into time delays that are on average a factor about 2 smaller than those obtained when the phase diagram of Segretain & Chabrier (1993) is employed. Our results have important implications for white dwarf cosmochronology, because the cooling ages of very old white dwarfs are different for the two phase diagrams. This may have a noticeable impact on the age determinations of very old globular clusters, for which the white dwarf color-magnitude diagram provides an independent way of estimating their age.

We present a nucleosynthesis study of the merger of a 0.4 solar masses helium white dwarf with a 0.8 solar masses carbon-oxygen white dwarf, coupling the thermodynamic history of Smoothed Particle Hydrodynamics particles with a post-processing code. The resulting chemical abundance pattern, particularly for oxygen and fluorine, is in qualitative agreement with the observed abundances in R Coronae Borealis stars.

We perform an asteroseismological study on the DAV star G117-B15A on the basis of a modern set of fully evolutionary DA white dwarf models that have consistent chemical profiles at the core and the envelope. We found an asteroseismological model for G117-B15A that closely reproduces its observed pulsation periods. Then, we use the most recently measured value of the rate of period change for the dominant mode of this pulsating star to impose a preliminary upper limit to the mass of the axion.

We investigate the internal rotation profile of the GW Vir (PG1159-type) star PG 0122+200 by employing an asteroseismological model that closely reproduces the observed periods of this star. We adopt a forward approach and two inversion methods based on the rotational splitting of the pulsation frequencies to explore the properties of the rotation of PG 0122+200. We found evidence for differential rotation in this star.

NGC 6791 is a well-studied, metal-rich open cluster that is so close to us that can be imaged down to luminosities fainter than that of the termination of its white dwarf cooling sequence, thus allowing for an in-depth study of its white dwarf population. We use a Monte Carlo simulator that employs up-to-date evolutionary cooling sequences for white dwarfs with hydrogen-rich and hydrogen-deficient atmospheres, with carbon-oxygen and helium cores. The cooling sequences for carbon-oxygen cores account for the delays introduced by both Ne^22 sedimentation in the liquid phase and by carbon-oxygen phase separation upon crystallization. We do not find evidence for a substantial fraction of helium-core white dwarfs, and hence our results support the suggestion that the origin of the bright peak of the white dwarf luminosity function can only be attributed to a population of unresolved binary white dwarfs. Moreover, our results indicate that the number distribution of secondary masses of the population of unresolved binaries has to increase with increasing mass ratio between the secondary and primary components of the progenitor system. We also find that the observed cooling sequence appears to be able to constrain the presence of progenitor sub-populations with different chemical compositions and the fraction of non-DA white dwarfs. Our simulations place interesting constraints on important characteristics of the stellar populations of NGC 6791. In particular, we find that the fraction of single helium-core white dwarfs must be smaller than 5%, that a sub-population of stars with zero metallicity must be <12%, while if the adopted metallicity of the sub-population is solar the upper limit is ~8%. Finally, we also find that the fraction of non-DA white dwarfs in this particular cluster is surprinsingly small <6%.

Many hydrogen deficient stars are characterised by surface abundance patterns that are hard to reconcile with conventional stellar evolution. Instead, it has been suggested that they may represent the result of a merger episode between a helium and a carbon-oxygen white dwarf. In this Letter, we present a nucleosynthesis study of the merger of a 0.4 M_sol helium white dwarf with a 0.8 M_sol carbon-oxygen white dwarf, by coupling the thermodynamic history of Smoothed Particle Hydrodynamics particles with a post-processing code. The resulting chemical abundance pattern, particularly for oxygen and fluorine, is in qualitative agreement with the observed abundances in R Coronae Borealis stars.

A variation of the gravitational constant over cosmological ages modifies the main sequence lifetimes and white dwarf cooling ages. Using an state-of-the-art stellar evolutionary code we compute the effects of a secularly varying G on the main sequence ages and, employing white dwarf cooling ages computed taking into account the effects of a running G, we place constraints on the rate of variation of Newton's constant. This is done using the white dwarf luminosity function and the distance of the well studied open Galactic cluster NGC 6791. We derive an upper bound G'/G ~ -1.8 10^{-12} 1/yr. This upper limit for the secular variation of the gravitational constant compares favorably with those obtained using other stellar evolutionary properties, and can be easily improved if deep images of the cluster allow to obtain an improved white dwarf luminosity function.

LISA Pathfinder is a science and technology demonstrator of the European Space Agency within the framework of its LISA mission, which aims to be the first space-borne gravitational wave observatory. The payload of LISA Pathfinder is the so-called LISA Technology Package, which is designed to measure relative accelerations between two test masses in nominal free fall. Its disturbances are monitored and dealt by the diagnostics subsystem. This subsystem consists of several modules, and one of these is the magnetic diagnostics system, which includes a set of four tri-axial fluxgate magnetometers, intended to measure with high precision the magnetic field at the positions of the test masses. However, since the magnetometers are located far from the positions of the test masses, the magnetic field at their positions must be interpolated. It has been recently shown that because there are not enough magnetic channels, classical interpolation methods fail to derive reliable measurements at the positions of the test masses, while neural network interpolation can provide the required measurements at the desired accuracy. In this paper we expand these studies and we assess the reliability and robustness of the neural network interpolation scheme for variations of the locations and possible offsets of the magnetometers, as well as for changes in environmental conditions. We find that neural networks are robust enough to derive accurate measurements of the magnetic field at the positions of the test masses in most circumstances.

Gravitational settling of 22Ne in cooling white dwarfs can affect the outcome of thermonuclear supernovae. We investigate how the supernova energetics and nucleosynthesis are in turn influenced by this process. We use realistic chemical profiles derived from state-of-the-art white dwarf cooling sequences. The cooling sequences provide a link between the white dwarf chemical structure and the age of the supernova progenitor system. The cooling sequence of a 1 M_sun white dwarf was computed until freezing using an up-to-date stellar evolutionary code. We computed explosions of both Chandrasekhar mass and sub-Chandrasekhar mass white dwarfs, assuming spherical symmetry and neglecting convective mixing during the pre-supernova carbon simmering phase to maximize the effects of chemical separation. Neither gravitational settling of 22Ne nor chemical differentiation of 12C and 16O have an appreciable impact on the properties of Type Ia supernovae, unless there is a direct dependence of the flame properties (density of transition from deflagration to detonation) on the chemical composition. At a fixed transition density, the maximum variation in the supernova magnitude obtained from progenitors of different ages is ~0.06 magnitudes, and even assuming an unrealistically large diffusion coefficient of 22Ne it would be less than ~0.09 mag. However, if the transition density depends on the chemical composition (all other things being equal) the oldest SNIa can be as much as 0.4 magnitudes brighter than the youngest ones (in our models the age difference is 7.4 Gyr). In addition, our results show that 22Ne sedimentation cannot be invoked to account for the formation of a central core of stable neutron-rich Fe-group nuclei in the ejecta of sub-Chandrasekhar models, as required by observations of Type Ia supernovae.

Abridged. White dwarf stars are the final evolutionary stage of the vast majority of stars, including our Sun. The study of white dwarfs has potential applications to different fields of astrophysics. In particular, they 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, like our Galaxy and open and globular clusters. In addition, the high densities and temperatures characterizing white dwarfs allow to use these stars as cosmic laboratories for studying physical processes under extreme conditions that cannot be achieved in terrestrial laboratories. They 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. In this work, we review the essentials of the physics of white dwarf stars. 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 an unique opportunity to peer into their internal structure that would otherwise remain hidden from observers. We will show that this allows to measure with unprecedented precision the stellar masses 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.

Motivated by the strong discrepancy between the main sequence turn-off age and the white dwarf cooling age in the metal-rich open cluster NGC 6791, we compute a grid of white dwarf evolutionary sequences that incorporates for the first time the energy released by the processes of 22Ne sedimentation and of carbon/oxygen phase separation upon crystallization. The grid covers the mass range from 0.52 to 1.0 Msun, and it is appropriate for the study of white dwarfs in metal-rich clusters. The evolutionary calculations are based on a detailed and self-consistent treatment of the energy released from these two processes, as well as on the employment of realistic carbon/oxygen profiles, of relevance for an accurate evaluation of the energy released by carbon/oxygen phase separation. We find that 22Ne sedimentation strongly delays the cooling rate of white dwarfs stemming from progenitors with high metallicities at moderate luminosities, whilst carbon/oxygen phase separation adds considerable delays at low luminosities. Cooling times are sensitive to possible uncertainties in the actual value of the diffusion coefficient of 22Ne. Changing the diffusion coefficient by a factor of 2, leads to maximum age differences of approx. 8-20% depending on the stellar mass. We find that the magnitude of the delays resulting from chemical changes in the core is consistent with the slow down in the white dwarf cooling rate that is required to solve the age discrepancy in NGC 6791.

We compute new chemical profiles for the core and envelope of white dwarfs appropriate for pulsational studies of ZZ Ceti stars. These profiles are extracted from the complete evolution of progenitor stars, evolved through the main sequence and the thermally-pulsing asymptotic giant branch (AGB) stages, and from time-dependent element diffusion during white dwarf evolution. We discuss the importance of the initial-final mass relationship for the white dwarf carbon-oxygen composition. In particular, we find that the central oxygen abundance may be underestimated by about 15% if the white dwarf mass is assumed to be the hydrogen-free core mass before the first thermal pulse. We also discuss the importance for the chemical profiles expected in the outermost layers of ZZ Ceti stars of the computation of the thermally-pulsing AGB phase and of the phase in which element diffusion is relevant. We find a strong dependence of the outer layer chemical stratification on the stellar mass. In particular, in the less massive models, the double-layered structure in the helium layer built up during the thermally-pulsing AGB phase is not removed by diffusion by the time the ZZ Ceti stage is reached. Finally, we perform adiabatic pulsation calculations and discuss the implications of our new chemical profiles for the pulsational properties of ZZ Ceti stars. We find that the whole $g-$mode period spectrum and the mode-trapping properties of these pulsating white dwarfs as derived from our new chemical profiles are substantially different from those based on chemical profiles widely used in existing asteroseismological studies. Thus, we expect the asteroseismological models derived from our chemical profiles to be significantly different from those found thus far.

NGC 6791 is a well studied open cluster1 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 significantly 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 crystallise 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 significant 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, solving the age discrepancy for NGC 6791.

The collision of two white dwarfs is a quite frequent event in dense stellar systems, like globular clusters and galactic nuclei. In this paper we present the results of a set of simulations of the close encounters and collisions of two white dwarfs. 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 \times 10^5) and covering a wide range of velocities and initial distances of the colliding white dwarfs. We discuss in detail when the initial eccentric binary white dwarf survives the closest approach, when a lateral collision in which several mass transfer episodes occur is the outcome of the newly formed binary system, and which range of input parameters leads to a direct collision, in which only one mass transfer episode occurs. We also discuss the characteristics of the final configuration and we assess the possible observational signatures of the merger, such as the associated gravitational waveforms and the fallback luminosities. We find that the overall evolution of the system and the main characteristics of the final object agree with those found in previous studies. We also find that the fallback luminosities are close to 10^48 erg/s. Finally, we find as well that in the case of lateral and direct collisions the gravitational waveforms are characterized by large-amplitude peaks which are followed by a ring-down phase, while in the case in which the binary white dwarf survives the closest approach, the gravitational pattern shows a distinctive behavior, typical of eccentric systems.

Axions are the natural consequence of the introduction of the Peccei-Quinn symmetry to solve the strong CP problem. All the efforts to detect such elusive particles have failed up to now. Nevertheless, it has been recently shown that the luminosity function of white dwarfs is best fitted if axions with a mass of a few meV are included in the evolutionary calculations. Our aim is to show that variable white dwarfs can provide additional and independent evidence about the existence of axions. The evolution of a white dwarf is a slow cooling process that translates into a secular increase of the pulsation periods of some variable white dwarfs, the so-called DAV and DBV types. Since axions can freely escape from such stars, their existence would increase the cooling rate and, consequently, the rate of change of the periods as compared with the standard ones. The present values of the rate of change of the pulsation period of G117-B15A are compatible with the existence of axions with the masses suggested by the luminosity function of white dwarfs, in contrast with previous estimations. Furthermore, it is shown that if such axions indeed exist, the drift of the periods of pulsation of DBV stars would be noticeably perturbed.

We present full evolutionary calculations appropriate for the study of hot hydrogen-deficent DO white dwarfs, PG 1159 stars, and DB white dwarfs. White dwarf sequences are computed for a wide range of stellar masses and helium envelopes on the basis of a complete treatment of the evolutionary history of progenitors stars, including the core hydrogen and helium burning phases, the thermally-pulsing AGB phase, and the born-again episode that is responsible for the hydrogen deficiency. We also provide colors and magnitudes for the new sequences for $T_{\rm eff} < 40 000$ K, where the NLTE effects are not dominant. These new calculations provide an homogeneous set of evolutionary tracks appropriate for mass and age determinations for both PG 1159 stars and DO white dwarfs. The calculations are extended down to an effective temperature of 7 000 K. We applied these new tracks to redetermine stellar masses and ages of all known DO white dwarfs with spectroscopically-determined effective temperatures and gravities, and compare them with previous results. We also compare for the first time consistent mass determinations for both DO and PG 1159 stars, and find a considerably higher mean mass for the DO white dwarfs. We discuss as well the chemical profile expected in the envelope of variable DB white dwarfs from the consideration of the evolutionary history of progenitor stars. Finally, we present tentative evidence for a different evolutionary channel, other than that involving the PG 1159 stars, for the formation of hot, hydrogen-deficient white dwarfs.

Hot DQ white dwarfs constitute a new class of white dwarf stars, uncovered recently within the framework of the SDSS project. Recently, three hot DQ white dwarfs have been reported to exhibit photometric variability with periods compatible with pulsation g-modes. Here, we present a nonadiabatic pulsation analysis of the recently discovered carbon-rich hot DQ white dwarf stars. One of our main aims is to test the convective-mixing picture for the origin of hot DQs through the pulsational properties. Our study relies on the full evolutionary models of hot DQ white dwarfs recently developed by Althaus et al. (2009), that consistently cover the whole evolution from the born-again stage to the white dwarf cooling track. Specifically, we present a stability analysis of white dwarf models from stages before the blue edge of the DBV instability strip until the domain of the hot DQ white dwarfs, including the transition DB --> hot DQ white dwarf. We found that at evolutionary phases in which the models are characterized by He-dominated atmospheres, they exhibit unstable $g$-mode pulsations typical of DBV stars, and when the models become DQ white dwarfs with carbon-dominated atmospheres, they continue being pulsationally unstable with similar characteristics than DB models, and in agreement with the periods detected in variable hot DQ white dwarfs. Our calculations provide strong support to the convective-mixing picture for the formation of hot DQs. In particular, our results point to the existence of pulsating DB white dwarfs with very thin He-rich envelopes, which after passing the DBV instability strip become variable hot DQ stars. The existence of these DB stars with very thin envelopes could be investigated through asteroseismology.

Motivated by the recent detection of single and binary He-core white dwarfs in metal-rich clusters, we present a full set of evolutionary calculations and colors appropriate for the study of such white dwarfs. The paper is also aimed at investigating whether stable hydrogen burning may constitute a main source of energy for massive He-core white dwarfs resulting from high-metallicity progenitors. White dwarf sequences are derived by taking into account the evolutionary history of progenitor stars with supersolar metallicities. We also incorporate a self-consistent, time-dependent treatment of gravitational settling and chemical diffusion, as well as of the residual nuclear burning. We find that the influence of residual nuclear burning during the late stages of white dwarf evolution is strongly dependent on the occurrence of chemical diffusion at the base of the hydrogen-rich envelope. When no diffusion is considered, residual hydrogen burning strongly influences the advanced stages of white dwarf cooling, introducing evolutionary delays of several Gyr. By contrast, when diffusion is taken into account the role of residual nuclear burning is strongly mitigated, and the evolution is dictated only by the thermal content stored in the ions. In addition, for all of our sequences, we provide accurate color and magnitudes on the basis of new and improved non gray model atmospheres which explicitly include Ly$\alpha$ quasi-molecular opacity.

We present the results of a set of high-resolution simulations of the merging process of two white dwarfs. In order to do so, we use an up-to-date Smoothed Particle Hydrodynamics code which incorporates very detailed input physics and an improved treatment of the artificial viscosity. Our simulations have been done using a large number of particles (4x10^5) and cover 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 that 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, like the associated gravitational waveforms and the fallback X-ray flares, and we study the long-term evolution of the coalescence.

We present the first full evolutionary calculations aimed at exploring the origin of hot DQ white dwarfs. These calculations consistently cover the whole evolution from the born-again stage to the white dwarf cooling track. Our calculations provide strong support to the diffusive/convective-mixing picture for the formation of hot DQs. We find that the hot DQ stage is a short-lived stage and that the range of effective temperatures where hot DQ stars are found can be accounted for by different masses of residual helium and/or different initial stellar masses. In the frame of this scenario, a correlation between the effective temperature and the surface carbon abundance in DQs should be expected, with the largest carbon abundances expected in the hottest DQs. From our calculations, we suggest that most of the hot DQs could be the cooler descendants of some PG1159 stars characterized by He-rich envelopes markedly smaller than those predicted by the standard theory of stellar evolution. At least for one hot DQ, the high-gravity white dwarf SDSS J142625.70+575218.4, an evolutionary link between this star and the massive PG1159 star H1504+65 is plausible.

We explore different evolutionary scenarios to explain the helium deficiency observed in H1504+65, the most massive known PG1159 star. We concentrate mainly on the possibility that this star could be the result of mass loss shortly after the born-again and during the subsequent evolution through the [WCL] stage. This possibility is sustained by recent observational evidence of extensive mass-loss events in Sakurai's object and is in line with the recent finding that such mass losses give rise to PG1159 models with thin helium-rich envelopes and large rates of period change, as demanded by the pulsating star PG1159-035. We compute the post born again evolution of massive sequences by taking into account different mass-loss rate histories. Our results show that stationary winds during the post-born-again evolution fail to remove completely the helium-rich envelope so as to explain the helium deficiency observed in H1504+65. Stationary winds during the Sakurai and [WCL] stages only remove at most half of the envelope surviving the violent hydrogen burning during the born-again phase. In view of our results, the recently suggested evolutionary connection born-again stars --> H1504+65 --> white dwarfs with carbon-rich atmospheres is difficult to sustain unless the whole helium-rich envelope could be ejected by non-stationary mass-loss episodes during the Sakurai stage.

We reexamine the theoretical instability domain of pulsating DB white dwarfs (DBV or V777 Her variables). We performed an extensive $g$-mode nonadiabatic pulsation analysis of DB evolutionary models considering a wide range of stellar masses, for which the complete evolutionary stages of their progenitors from the ZAMS, through the thermally pulsing AGB and born-again phases, the domain of the PG1159 stars, the hot phase of DO white dwarfs, and then the DB white dwarf stage have been considered. We explicitly account for the evolution of the chemical abundance distribution due to time-dependent chemical diffusion processes. We examine the impact of the different prescriptions of the MLT theory of convection and the effects of small amounts of H in the almost He-pure atmospheres of DB stars on the precise location of the theoretical blue edge of the DBV instability strip.

We present evolutionary calculations and colors for massive white dwarfs with oxygen-neon cores for masses between 1.06 and 1.28 Mo. The evolutionary stages computed cover the luminosity range from log(L/Lo) approx. 0.5 down to -5.2. Our cooling sequences are based on evolutionary calculations that take into account the chemical composition expected from massive white dwarf progenitors that burned carbon in partially degenerate conditions. The use of detailed non-gray model atmospheres provides us with accurate outer boundary conditions for our evolving models at low effective temperatures. We examine the cooling age, colors and magnitudes of our sequences. We find that massive white dwarfs are characterized by very short ages to such an extent that they reach the turn-off in their colors and become blue at ages well below 10 Gyr. Extensive tabulations for massive white dwarfs, accessible from our web site, are also presented.

We study the formation and evolution of DA white dwarfs, the progenitors of which have experienced a late thermal pulse (LTP) shortly after the departure from the thermally pulsing AGB. To this end, we compute the complete evolution of an initially 2.7 Mo star all the way from the zero-age main sequence to the white dwarf stage. We find that most of the original H-rich material of the post-AGB remnant is burnt during the post-LTP evolution, with the result that, at entering its white dwarf cooling track, the remaining H envelope becomes 10^-6 Mo in agreement with asteroseismological inferences for some ZZ Ceti stars.