The evolution of white dwarfs is a simple gravothermal process of cooling. Since the shape of their luminosity function is sensitive to the characteristic cooling time, it is possible to use its slope to test the existence of additional sources or sinks of energy, such as those predicted by alternative physical theories. The aim of this paper is to study if the changes in the slope of the white dwarf luminosity function around bolometric magnitudes ranging from 8 to 10 and previously attributed to axion emission are, effectively, a consequence of the existence of axions and not an artifact introduced by the star formatio rate. We compute theoretical luminosity functions of the thin and thick disk, and of the stellar halo including axion emission and we compare them with the existing observed luminosity functions. Since these stellar populations have different star formation histories, the slope change should be present in all of them at the same place if it is due to axions or any other intrinsic cooling mechanism. The signature of an unexpected cooling seems to be present in the luminosity functions of the thin and thick disks, as well as in the halo luminosity function. This additional cooling is compatible with axion emission, thus supporting the idea that DFSZ axions, with a mass in the range of 4 to 10 meV, could exist. If this were the case, these axions could be detected by the future solar axioscope IAXO.

The evolution of the remnant of the merger of two white dwarfs is still an open problem. Furthermore, few studies have addressed the case in which the remnant is a magnetic white dwarf with a mass larger than the Chandrasekhar limiting mass. Angular momentum losses might bring the remnant of the merger to the physical conditions suitable for developing a thermonuclear explosion. Alternatively, the remnant may be prone to gravitational or rotational instabilities, depending on the initial conditions reached after the coalescence. Dipole magnetic braking is one of the mechanisms that can drive such losses of angular momentum. However, the timescale on which these losses occur depend on several parameters, like the strength of the magnetic field. In addition, the coalescence leaves a surrounding Keplerian disk that can be accreted by the newly formed white dwarf. Here we compute the post-merger evolution of a super-Chandrasekhar magnetized white dwarf taking into account all the relevant physical processes. These include magnetic torques acting on the star, accretion from the Keplerian disk, the threading of the magnetic field lines through the disk, as well as the thermal evolution of the white dwarf core. We find that the central remnant can reach the conditions suitable to develop a thermonuclear explosion before other instabilities (such as the inverse beta-decay instability or the secular axisymmetric instability) are reached, which would instead lead to gravitational collapse of the magnetized remnant.

Context. Classical novae are thermonuclear explosions hosted by accreting white dwarfs in stellar binary systems. Material piles up on top of the white dwarf star under mildly degenerate conditions, driving a thermonuclear runaway. The energy released by the suite of nuclear processes operating at the envelope (mostly proton-capture reactions and beta-decays) heats the material up to peak temperatures ranging from 100 to 400 MK. In these events, about 10-3 - 10-7 Msun, enriched in CNO and, sometimes, other intermediate-mass elements (e.g., Ne, Na, Mg, Al), are ejected into the interstellar medium. Aims. To date, most of the efforts undertaken in the modeling of classical nova outbursts have focused on the early stages of the explosion and ejection, ignoring the interaction of the ejecta, first with the accretion disk orbiting the white dwarf, and ultimately with the secondary star. Methods. A suite of three-dimensional, SPH simulations of the interaction between the nova ejecta, the accretion disk, and the stellar companion have been performed to fill this gap, aimed at testing the influence of the different parameters (i.e., mass and velocity of the ejecta, mass and geometry of the accretion disk) on the dynamical and chemical properties of the system. Results. We discuss the conditions that lead to the disruption of the accretion disk and to mass loss from the binary system. In addition, we discuss the likelihood of chemical contamination of the stellar secondary induced by the impact with the nova ejecta and its potential effect on the next nova cycle.

We present the second paper of a series of publications aiming at obtaining a better understanding regarding the nature of type Ia supernovae (SNIa) progenitors by studying a large sample of detached F, G and K main sequence stars in close orbits with white dwarf companions (i.e. WD+FGK binaries). We employ the LAMOST (Large Sky Area Multi-Object Fibre Spectroscopic Telescope) data release 4 spectroscopic data base together with GALEX (Galaxy Evolution Explorer) ultraviolet fluxes to identify 1,549 WD+FGK binary candidates (1,057 of which are new), thus doubling the number of known sources. We measure the radial velocities of 1,453 of these binaries from the available LAMOST spectra and/or from spectra obtained by us at a wide variety of different telescopes around the globe. The analysis of the radial velocity data allows us to identify 24 systems displaying more than 3sigma radial velocity variation that we classify as close binaries. We also discuss the fraction of close binaries among WD+FGK systems, which we find to be ~10 per cent, and demonstrate that high-resolution spectroscopy is required to efficiently identify double-degenerate SNIa progenitor candidates.

We present a Monte Carlo population synthesis study of white dwarf-main sequence (WD+MS) binaries in the Galactic disk aimed at reproducing the ensemble properties of the entire population observed by the Sloan Digital Sky Survey (SDSS) Data Release 12. Our simulations take into account all known observational biases and use the most up-to-date stellar evolutionary models. This allows us to perform a sound comparison between the simulations and the observational data. We find that the properties of the simulated and observed parameter distributions agree best when assuming low values of the common envelope efficiency (0.2-0.3), a result that is in agreement with previous findings obtained by observational and population synthesis studies of close SDSS WD+MS binaries. We also show that all synthetic populations that result from adopting an initial mass ratio distribution with a positive slope are excluded by observations. Finally, we confirm that the properties of the simulated WD+MS binary populations are nearly independent of the age adopted for the thin disk, on the contribution of WD+MS binaries from the thick disk (0-17 per cent of the total population) and on the assumed fraction of the internal energy that is used to eject the envelope during the common envelope phase (0.1-0.5).

We use the Sloan Digital Sky Survey Data Release 12, which is the largest available white dwarf catalog to date, to study the evolution of the kinematical properties of the population of white dwarfs in the Galactic disc. We derive masses, ages, photometric distances and radial velocities for all white dwarfs with hydrogen-rich atmospheres. For those stars for which proper motions from the USNO-B1 catalog are available the true three-dimensional components of the stellar space velocity are obtained. This subset of the original sample comprises 20,247 objects, making it the largest sample of white dwarfs with measured three-dimensional velocities. Furthermore, the volume probed by our sample is large, allowing us to obtain relevant kinematical information. In particular, our sample extends from a Galactocentric radial distance $R_{\rm G}=7.8$~kpc to 9.3~kpc, and vertical distances from the Galactic plane ranging from $Z=-0.5$~kpc to 0.5~kpc. We examine the mean components of the stellar three-dimensional velocities, as well as their dispersions with respect to the Galactocentric and vertical distances. We confirm the existence of a mean Galactocentric radial velocity gradient, $\partial\langle V_{\rm R}\rangle/\partial R_{\rm G}=-3\pm5$~km~s$^{-1}$~kpc$^{-1}$. We also confirm North-South differences in $\langle V_{\rm z}\rangle$. Specifically, we find that white dwarfs with $Z>0$ (in the North Galactic hemisphere) have $\langle V_{\rm z}\rangle<0$, while the reverse is true for white dwarfs with $Z<0$. The age-velocity dispersion relation derived from the present sample indicates that the Galactic population of white dwarfs may have experienced an additional source of heating, which adds to the secular evolution of the Galactic disc.

The merger of close double white dwarfs (CDWDs) is one of the favourite evolutionary channels for producing Type Ia supernovae (SN Ia). Unfortunately, current theories of the evolution and formation of CDWDs are still poorly constrained and have several serious uncertainties, which affect the predicted SN Ia rates. Moreover, current observational constraints on this evolutionary pathway for SN Ia mainly rely on only 18 double-lined and/or eclipsing CDWDs with measured orbital and stellar parameters for both white dwarfs. In this paper we present the orbital periods and the individual masses of three new double-lined CDWDs, derived using a new method. This method employs mass ratios, the Halpha core ratios and spectral model-fitting to constrain the masses of the components of the pair. The three CDWDs are WD0028-474 (Porb=9.350 +- 0.007 hours, M1=0.60 +- 0.06 Msun, M2=0.45 +- 0.04 Msun), HE0410-1137 (Porb = 12.208 +- 0.008 hours, M1= 0.51 +- 0.04 Msun, M2= 0.39 +- 0.03 Msun) and SDSSJ031813.25-010711.7 (Porb = 45.908 +- 0.006 hours, among the longest period systems, M1= 0.40 +- 0.05 Msun, M2= 0.49 +- 0.05 Msun). While the three systems studied here will merge in timescales longer than the Hubble time and are expected to become single massive (>~0.9 Msun) white dwarfs rather than exploding as SN Ia, increasing the small sample of CDWDs with determined stellar parameters is crucial for a better overall understanding of their evolution.

Double-degenerate (DD) mergers of carbon-oxygen white dwarfs have recently emerged as a leading candidate for normal Type Ia supernovae (SNe Ia). However, many outstanding questions surround DD mergers, including the characteristics of their light curves and spectra. We have recently identified a spiral instability in the post-merger phase of DD mergers and demonstrated that this instability self-consistently leads to detonation in some cases. We call this the spiral merger SN Ia model. Here, we utilize the SuperNu radiative transfer software to calculate three-dimensional synthetic light curves and spectra of the spiral merger simulation with a system mass of 2.1 $M_\odot$ from Kashyap et al. Because of their large system masses, both violent and spiral merger light curves are slowly declining. The spiral merger resembles very slowly declining SNe Ia, including SN 2001ay, and provides a more natural explanation for its observed properties than other SN Ia explosion models. Previous synthetic light curves and spectra of violent DD mergers demonstrate a strong dependence on viewing angle, which is in conflict with observations. Here, we demonstrate that the light curves and spectra of the spiral merger are less sensitive to the viewing angle than violent mergers, in closer agreement with observation. We find that the spatial distribution of 56Ni and IMEs follows a characteristic hourglass shape. We discuss the implications of the asymmetric distribution of 56Ni for the early-time gamma-ray observations of 56Ni from SN 2014J. We suggest that DD mergers that agree with the light curves and spectra of normal SNe Ia will likely require a lower system mass.

Classical novae are thermonuclear explosions that take place in the envelopes of accreting white dwarfs in binary systems. The material piles up under degenerate conditions, driving a thermonuclear runaway. The energy released by the suite of nuclear processes operating at the envelope heats the material up to peak temperatures about 100 - 400 MK. During these events, about 10-3 - 10-7 Msun, enriched in CNO and, sometimes, other intermediate-mass elements (e.g., Ne, Na, Mg, Al) are ejected into the interstellar medium. To account for the gross observational properties of classical novae (in particular, the large concentrations of metals spectroscopically inferred in the ejecta), models require mixing between the (solar-like) material transferred from the secondary and the outermost layers (CO- or ONe-rich) of the underlying white dwarf. Recent multidimensional simulations have demonstrated that Kelvin-Helmholtz instabilities can naturally produce self-enrichment of the accreted envelope with material from the underlying white dwarf at levels that agree with observations. However, the feasibility of this mechanism has been explored in the framework of CO white dwarfs, while mixing with different substrates still needs to be properly addressed. Three-dimensional simulations of mixing at the core-envelope interface during nova outbursts have been performed with the multidimensional code FLASH, for two types of substrates: CO- and ONe-rich. We show that the presence of an ONe-rich substrate, as in "neon novae", yields larger metallicity enhancements in the ejecta, compared to CO,rich substrates (i.e., non-neon novae). A number of requirements and constraints for such 3-D simulations (e.g., minimum resolution, size of the computational domain) are also outlined.

We present an updated version of the spectroscopic catalogue of white dwarf-main sequence (WDMS) binaries from the Sloan Digital Sky Survey (SDSS). We identify 939 WDMS binaries within the data releases (DR) 9-12 of SDSS plus 40 objects from DR 1-8 that we missed in our previous works, 646 of which are new. The total number of spectroscopic SDSS WDMS binaries increases to 3294. This is by far the largest and most homogeneous sample of compact binaries currently available. We use a decomposition/fitting routine to derive the stellar parameters of all systems identified here (white dwarf effective temperatures, surface gravities and masses, and secondary star spectral types). The analysis of the corresponding stellar parameter distributions shows that the SDSS WDMS binary population is seriously affected by selection effects. We also measure the NaI 8183.27, 8194.81 absorption doublet and Halpha emission radial velocities (RV) from all SDSS WDMS binary spectra identified in this work. 98 objects are found to display RV variations, 62 of which are new. The RV data are sufficient enough to estimate the orbital periods of three close binaries.

We critically discuss the recent observations of the binary system at the center of the bipolar planetary nebula Henize 2-428. We find that the proposed explanation of two equal-mass degenerate objects with a total mass larger than the Chandrasekhar limiting mass that supposedly will merge in less than a Hubble time, possibly leading to a SN~Ia, is controversial. This hypothesis relies on the assumption that the variability of the He II 5412A spectral line is due to two absorption components. Instead, we propose that it can be accounted for by a broad absorption line from the central system on top of which there is a narrow emission line from the nebula. This prompted us to study if the binary system can be made of a degenerate star and a low-mass main sequence companion, or of two degenerate objects of smaller mass. We find that although both scenarios can account for the existence of two symmetric broad minima in the light curve, the second one agrees better with observations. We thus argue that the claim that Henize 2-428 provides observational evidence supporting the double-degenerate scenario for SN~Ia is premature.

The mass function of hydrogen-rich atmosphere white dwarfs has been frequently found to reveal a distinctive high-mass excess near 1Msun. However, a significant excess of massive white dwarfs has not been detected in the mass function of the largest white dwarf catalogue to date from the Sloan Digital Sky Survey. Hence, whether a high-mass excess exists or not has remained an open question. In this work we build the mass function of the latest catalogue of data release 10 SDSS hydrogen-rich white dwarfs, including the cool and faint population (i.e. effective temperatures 6,000 <~ Teff <~ 12,000 K, equivalent to 12 mag <~ Mbol <~ 13 mag). We show that the high-mass excess is clearly present in our mass function, and that it disappears only if the hottest (brightest) white dwarfs (those with Teff >~ 12,000 K, Mbol <~ 12 mag) are considered. This naturally explains why previous SDSS mass functions failed at detecting a significant excess of high-mass white dwarfs. Thus, our results provide additional and robust observational evidence for the existence of a distinctive high-mass excess near 1Msun. We investigate possible origins of this feature and argue that the most plausible scenario that may lead to an observed excess of massive white dwarfs is the merger of the degenerate core of a giant star with a main sequence or a white dwarf companion during or shortly after a common envelope event.

Gaia will identify several 1e5 white dwarfs, most of which will be in the solar neighborhood at distances of a few hundred parsecs. Ground-based optical follow-up spectroscopy of this sample of stellar remnants is essential to unlock the enormous scientific potential it holds for our understanding of stellar evolution, and the Galactic formation history of both stars and planets.

Modern large-scale surveys have allowed the identification of large numbers of white dwarfs. However, these surveys are subject to complicated target selection algorithms, which make it almost impossible to quantify to what extent the observational biases affect the observed populations. The LAMOST (Large Sky Area Multi-Object Fiber Spectroscopic Telescope) Spectroscopic Survey of the Galactic anti-center (LSS-GAC) follows a well-defined set of criteria for selecting targets for observations. This advantage over previous surveys has been fully exploited here to identify a small yet well-characterised magnitude-limited sample of hydrogen-rich (DA) white dwarfs. We derive preliminary LSS-GAC DA white dwarf luminosity and mass functions. The space density and average formation rate of DA white dwarfs we derive are 0.83+/-0.16 x 10^{-3} pc^{-3} and 5.42 +/- 0.08 x 10^{-13} pc^{-3} yr^{-1}, respectively. Additionally, using an existing Monte Carlo population synthesis code we simulate the population of single DA white dwarfs in the Galactic anti-center, under various assumptions. The synthetic populations are passed through the LSS-GAC selection criteria, taking into account all possible observational biases. This allows us to perform a meaningful comparison of the observed and simulated distributions. We find that the LSS-GAC set of criteria is highly efficient in selecting white dwarfs for spectroscopic observations (80-85 per cent) and that, overall, our simulations reproduce well the observed luminosity function. However, they fail at reproducing an excess of massive white dwarfs present in the observed mass function. A plausible explanation for this is that a sizable fraction of massive white dwarfs in the Galaxy are the product of white dwarf-white dwarf mergers.

The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3 Cosmic Vision framework and participated in the final down-selection for a launch slot in 2022-2024. Thanks to the unprecedented combination of effective area and spectral resolution of its main instrument, LOFT will study the behaviour of matter under extreme conditions, such as the strong gravitational field in the innermost regions of accretion flows close to black holes and neutron stars, and the supra-nuclear densities in the interior of neutron stars. The science payload is based on a Large Area Detector (LAD, 10 m 2 effective area, 2-30 keV, 240 eV spectral resolution, 1 deg collimated field of view) and a WideField Monitor (WFM, 2-50 keV, 4 steradian field of view, 1 arcmin source location accuracy, 300 eV spectral resolution). The WFM is equipped with an on-board system for bright events (e.g. GRB) localization. The trigger time and position of these events are broadcast to the ground within 30 s from discovery. In this paper we present the status of the mission at the end of its Phase A study.

This is the product of a Gaia GREAT workshop on "Gaia and the end states of stellar evolution" held at the University of Leicester in April 2011.

We argue that the properties of the Type Ia supernova (SN Ia) SN 2011fe can be best explained within the frame of the core-degenerate (CD) scenario. In the CD scenario a white dwarf (WD) merges with the core of an asymptotic giant branch (AGB) star and forms a rapidly rotating WD, with a mass close to and above the critical mass for explosion. Rapid rotation prevents immediate collapse and/or explosion. Spinning down over a time of 0-10 Gyr brings the WD to explosion. A very long delayed explosion to post-crystallization phase, which lasts for about 2 Gyr leads to the formation of a highly carbon-enriched outer layer. This can account for the carbon-rich composition of the fastest-moving ejecta of SN 2011fe. In reaching the conclusion that the CD scenario best explains the observed properties of SN 2011fe we consider both its specific properties, like a very compact exploding object and carbon rich composition of the fastest-moving ejecta, and the general properties of SNe Ia.

Type Ia supernovae (SNe Ia) play a crucial role as standardizable cosmological candles, though the nature of their progenitors is a subject of active investigation. Recent observational and theoretical work has pointed to merging white dwarf binaries, referred to as the double-degenerate channel, as the possible progenitor systems for some SNe Ia. Additionally, recent theoretical work suggests that mergers which fail to detonate may produce magnetized, rapidly-rotating white dwarfs. In this paper, we present the first multidimensional simulations of the post-merger evolution of white dwarf binaries to include the effect of the magnetic field. In these systems, the two white dwarfs complete a final merger on a dynamical timescale, and are tidally disrupted, producing a rapidly-rotating white dwarf merger surrounded by a hot corona and a thick, differentially-rotating disk. The disk is strongly susceptible to the magnetorotational instability (MRI), and we demonstrate that this leads to the rapid growth of an initially dynamically weak magnetic field in the disk, the spin-down of the white dwarf merger, and to the subsequent central ignition of the white dwarf merger. Additionally, these magnetized models exhibit new features not present in prior hydrodynamic studies of white dwarf mergers, including the development of MRI turbulence in the hot disk, magnetized outflows carrying a significant fraction of the disk mass, and the magnetization of the white dwarf merger to field strengths $\sim 2 \times 10^8$ G. We discuss the impact of our findings on the origins, circumstellar media, and observed properties of SNe Ia and magnetized white dwarfs.

We discuss the impact of residual nuclear burning in the cooling sequences of hydrogen-rich DA white dwarfs with very low metallicity progenitors ($Z=0.0001$). These cooling sequences are appropriate for the study of very old stellar populations. The results presented here are the product of self-consistent, fully evolutionary calculations. Specifically, we follow the evolution of white dwarf progenitors from the zero-age main sequence through all the evolutionary phases, namely the core hydrogen-burning phase, the helium-burning phase, and the thermally pulsing asymptotic giant branch phase to the white dwarf stage. This is done for the most relevant range of main sequence masses, covering the most usual interval of white dwarf masses --- from $0.53\, M_{\sun}$ to $0.83\, M_{\sun}$. Due to the low metallicity of the progenitor stars, white dwarfs are born with thicker hydrogen envelopes, leading to more intense hydrogen burning shells as compared with their solar metallicity counterparts. We study the phase in which nuclear reactions are still important and find that nuclear energy sources play a key role during long periods of time, considerably increasing the cooling times from those predicted by standard white dwarf models. In particular, we find that for this metallicity and for white dwarf masses smaller than about $0.6\, M_{\sun}$, nuclear reactions are the main contributor to the stellar luminosity for luminosities as low as $\log(L/L_{\sun})\simeq -3.2$. This, in turn, should have a noticeable impact in the white dwarf luminosity function of low-metallicity stellar populations.

We argue that the multiple shells of circumstellar material (CSM) and the supernovae (SN) ejecta interaction with the CSM starting 59 days after the explosion of the Type Ia SN (SN Ia) PTF 11kx, are best described by a violent prompt merger. In this prompt merger scenario the common envelope (CE) phase is terminated by a merger of a WD companion with the hot core of a massive asymptotic giant (AGB) star. In most cases the WD is disrupted and accreted onto the more massive core. However, in the rare cases where the merger takes place when the WD is denser than the core, the core will be disrupted and accreted onto the cooler WD. In such cases the explosion might occur with no appreciable delay, i.e., months to years after the termination of the CE phase. This, we propose, might be the evolutionary route that could lead to the explosion of PTF 11kx. This scenario can account for the very massive CSM within ~1000 AU of the exploding PTF 11kx star, for the presence of hydrogen, and for the presence of shells in the CSM.

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

Recent surveys have allowed to derive the white dwarf mass distribution with reasonable accuracy. This distribution shows a noticeable degree of structure that it is often attributed to the evolution of close binaries in general, and to mergers in particular. To analyze if the origin of this structure can be attributed to the merger of double white dwarfs, we have used a simplified population synthesis model that retains the essential processes of formation of double degenerate binaries. Special care has been taken to avoid artifacts introduced by discontinuities in the distribution functions. Our result is that these structures are not probably due to mergers, but they can provide a deep insight on the evolution of close binary systems.

The LOFT mission concept is one of four candidates selected by ESA for the M3 launch opportunity as Medium Size missions of the Cosmic Vision programme. The launch window is currently planned for between 2022 and 2024. LOFT is designed to exploit the diagnostics of rapid X-ray flux and spectral variability that directly probe the motion of matter down to distances very close to black holes and neutron stars, as well as the physical state of ultra-dense matter. These primary science goals will be addressed by a payload composed of a Large Area Detector (LAD) and a Wide Field Monitor (WFM). The LAD is a collimated (<1 degree field of view) experiment operating in the energy range 2-50 keV, with a 10 m^2 peak effective area and an energy resolution of 260 eV at 6 keV. The WFM will operate in the same energy range as the LAD, enabling simultaneous monitoring of a few-steradian wide field of view, with an angular resolution of <5 arcmin. The LAD and WFM experiments will allow us to investigate variability from submillisecond QPO's to year-long transient outbursts. In this paper we report the current status of the project.

White dwarfs are almost completely degenerate objects that cannot obtain energy from thermonuclear sources, so their evolution is just a gravothermal cooling process. Recent improvements in the accuracy and precision of the luminosity function and in pulsational data of variable white dwarfs suggest that they are cooling faster than expected from conventional theory. In this contribution we show that the inclusion of an additional cooling term due to axions able to interact with electrons with a coupling constant g_ae ~(2-7)x10^{-13} allows to fit better the observations.

Within the theoretical framework of some modern unification theories the constants of nature are functions of cosmological time. White dwarfs offer the possibility of testing a possible variation of G and, thus, to place constraints to these theories. We present full white dwarf evolutionary calculations in the case that G decreases with time. White dwarf evolution is computed in a self-consistent way, including the most up-to-date physical inputs, non-gray model atmospheres and a detailed core chemical composition that results from the calculation of the full evolution of progenitor stars. We find that the mechanical structure and the energy balance of white dwarfs are strongly modified by the presence of a varying G. In particular, for certain values of the rate of change of G, the evolution of cool white dwarfs is markedly affected. The impact of a varying G is more notorious in the case of more massive white dwarfs. In view of the recent results reporting that a very accurate white dwarf cooling age can be derived for the old and metal-rich open cluster NGC 6791, our study suggests that this cluster could be a potential target to constrain or detect a ypothetical secular variation of G.

Classical novae are explosive phenomena that take place in stellar binary systems. They are powered by mass transfer from a low-mass, main sequence star onto a white dwarf. The material piles up under degenerate conditions and a thermonuclear runaway ensues. The energy released by the suite of nuclear processes operating at the envelope heats the material up to peak temperatures of ~ (1 - 4) \times 108 K. During these events, about 10-4 - 10-5 M\odot, enriched in CNO and other intermediate-mass elements, are ejected into the interstellar medium. To account for the gross observational properties of classical novae (in particular, a metallicity enhancement in the ejecta above solar values), numerical models assume mixing between the (solar-like) material transferred from the companion and the outermost layers (CO- or ONe-rich) of the underlying white dwarf. The nature of the mixing mechanism that operates at the core-envelope interface has puzzled stellar modelers for about 40 years. Here we investigate the role of Kelvin-Helmholtz instabilities as a natural mechanism for self-enrichment of the accreted envelope with core material. The feasibility of this mechanism is studied by means of the multidimensional code FLASH. Here, we present a series of 9 numerical simulations perfomed in two dimensions aimed at testing the possible influence of the initial perturbation (duration, strength, location, and size), the resolution adopted, or the size of the computational domain on the results. We show that results do not depend substantially on the specific choice of these parameters, demonstrating that Kelvin- Helmholtz instabilities can naturally lead to self-enrichment of the accreted envelope with core material, at levels that agree with observations.

White dwarfs are almost completely degenerate objects that cannot obtain energy from the thermonuclear sources and their evolution is just a gravothermal process of cooling. The simplicity of these objects, the fact that the physical inputs necessary to understand them are well identified, although not always well understood, and the impressive observational background about white dwarfs make them the most well studied Galactic population. These characteristics allow to use them as laboratories to test new ideas of physics. In this contribution we discuss the robustness of the method and its application to the axion case.

The LISA PathFinder DMU (Data Management Unit) flight model was formally accepted by ESA and ASD on 11 February 2010, after all hardware and software tests had been successfully completed. The diagnostics items are scheduled to be delivered by the end of 2010. In this paper we review the requirements and performance of this instrumentation, specially focusing on the Radiation Monitor and the DMU, as well as the status of their programmed use during mission operations, on which work is ongoing at the time of writing.

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 stellar masses 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 compute self-consistently 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, like 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 Lyman alpha quasi-molecular opacity. The calculations are extended down to an effective temperature of 2,500 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.

Classical novae are powered by thermonuclear runaways that occur on the white dwarf component of close binary systems. During these violent stellar events, whose energy release is only exceeded by gamma-ray bursts and supernova explosions, about 10-4 10-5 Msun of material is ejected into the interstellar medium. Because of the high peak temperatures attained during the explosion, Tpeak ~ (1-4)x10+8 K, the ejecta are enriched in nuclear-processed material relative to solar abundances, containing significant amounts of 13C, 15N, and 17O and traces of other isotopes. The origin of these metal enhancements observed in the ejecta is not wellknown and has puzzled theoreticians for about 40 years. In this paper, we present new 2-D simulations of mixing at the core-envelope interface. We show that Kelvin-Helmholtz instabilities can naturally lead to self-enrichment of the solar-like accreted envelopes with material from the outermost layers of the underlying white dwarf core, at levels that agree with observations.

The nature of the several microlensing events observed by the MACHO team towards the Large Magellanic Cloud (LMC) is still a subject of debate. Low-mass substellar objects and stars with masses larger than ~M_{sun} have been ruled out as major components of a Massive Astrophysical Compact Halo Object (MACHO) Galactic halo, while stars of half a solar mass seem to be viable candidates. Main sequence stars have been already discarded, and there are tight restrictions on the role played by white dwarfs with hydrogen-dominated atmospheres. In this paper we evaluate the contribution to the dark matter content of the Galaxy of white dwarfs with hydrogen-deficient atmospheres. For this purpose we use a Monte Carlo simulator which incorporates up-to-date evolutionary sequences of white dwarfs with hydrogen-rich and hydrogen-deficient atmospheres. We also take into account detailed descriptions of the thick disk and the halo of our Galaxy as well as of a reliable model of the LMC. We find that the contribution of white dwarfs with hydrogen-deficient atmospheres moderately increases the theoretical estimate of the optical depth with respect to the value obtained when only hydrogen-rich white dwarfs are considered. We also find that the contribuiton of the thick disk population of white dwarfs is comparable to the halo contribution. However, the contributions of both the halo and the thick disk white-dwarf populations are still insufficient to explain the number of events observed by the MACHO team. Finally, we find that the contribution to the halo dark matter of the entire population under study is less than 10% at the 95% conficence level.

White dwarfs are the end-product of the lifes of intermediate- and low-mass stars and their evolution is described as a simple cooling process. Recently, it has been possible to determine with an unprecedented precision their luminosity function, that is, the number of stars per unit volume and luminosity interval. We show here that the shape of the bright branch of this function is only sensitive to the averaged cooling rate of white dwarfs and we propose to use this property to check the possible existence of axions, a proposed but not yet detected weakly interacting particle. Our results indicate that the inclusion of the emission of axions in the evolutionary models of white dwarfs noticeably improves the agreement between the theoretical calculations and the observational white dwarf luminosity function. The best fit is obtained for m_a cos^2 \beta ~ 5 meV, where m_a is the mass of the axion and cos^2 \beta is a free parameter. We also show that values larger than 10 meV are clearly excluded. The existing theoretical and observational uncertainties do not allow yet to confirm the existence of axions, but our results clearly show that if their mass is of the order of few meV, the white dwarf luminosity function is sensitive enough to detect their existence.

The evolution of white dwarfs can be described as a simple cooling process. Recently, it has been possible to determine with an unprecedented precision their luminosity function, that is, the number of stars per unit volume and luminosity interval. Since the shape of the bright branch of this function is only sensitive to the average cooling rate, we use this property to check the possible existence of axions, a proposed but not yet detected weakly interacting particle. We show here that the inclusion of the axion emissivity in the evolutionary models of white dwarfs noticeably improves the agreement between the theoretical calculations and the observational white dwarf luminosity function, thus providing the first positive indication that axions could exist. Our results indicate that the best fit is obtained for m_a cos^2beta ~ 2-6 meV, where m_a is the mass of the axion and cos^2beta is a free parameter, and that values larger than 10 meV are clearly excluded.

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 semi-empirical initial-final mass relationship arising from our study covers the range of initial masses from 1.0 to 6.5 M_sun, 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 nature of the several microlensing events observed by the MACHO team towards the LMC still remains controversial. Low-mass substellar objects and stars with masses larger than ~1 M_{sun} have been ruled out as major components of a MACHO galactic halo, while stars of half solar masses are the most probable candidates. In this paper we assess jointly the relative contributions of both red dwarfs and white dwarfs to the mass budget of the galactic halo. In doing so we use a Monte Carlo simulator which incorporates up-to-date evolutionary sequences of both red dwarfs and white dwarfs as well as detailed descriptions of our Galaxy and of the LMC. We explore the complete mass range between 0.08 and 1 M_{sun} as possible microlensing candidates and we compare the synthetic populations obtained with our simulator with the results obtained by the MACHO and EROS experiments. Our results indicate that, despite that the contribution of the red dwarf population increases by a factor of 2 the value of the optical depth obtained when taking into account the white dwarf population alone, it is still insufficient to explain the number of events claimed by the MACHO team. Finally, we find that the contribution to the halo dark matter of the whole population under study is smaller than 10% at the 95% conficence level.

G39-27/289 is a common proper motion pair formed by a white dwarf (WD0433+270) and a main-sequence star (BD+26 730) that apparently has been classified as a member of the Hyades open cluster. Previous studies of the white dwarf component yielded a cooling time of ~4 Gyr. Although it has not been pointed out before explicitly, this result is 6 times larger than the age of the Hyades cluster, giving rise to an apparent conflict between the physics of white dwarfs and cluster main-sequence fitting. We investigate whether this system belongs to the Hyades cluster and, accordingly, give a plausible explanation to the nature of the white dwarf member. From kinematic and chemical composition considerations we believe that the system was a former member of the Hyades cluster and therefore has an evolutionary link with it. However, the evidence is not conclusive. With regards to the nature of the white dwarf component, we find that two core compositions (C/O and Fe) are compatible with the observed effective temperature and radius. These compositions yield very different cooling times of ~4 Gyr and ~1 Gyr, respectively. We distinguish two posssible scenarios. If the pair does not belong to the Hyades cluster but only to the Hyades stream, this would indicate that such stream contains rather old objects and definitely not coeval with the cluster. This has interesting consequences for Galactic dynamics. On the contrary, if the white dwarf has an Fe core, it will have a cooling time compatible with the Hyades age. This is a tantalizing result that would have implications for the thermonuclear explosion of white dwarfs and explosion theories of degenerate nuclei.

The initial-final mass relationship of white dwarfs, which is poorly constrained, is of paramount importance for different aspects in modern astrophysics. From an observational perspective, most of the studies up to now have been done using white dwarfs in open clusters. In order to improve the initial-final mass relationship we explore the possibility of deriving a semi-empirical relation studying white dwarfs in common proper motion pairs. We have acquired long-slit spectra of the white dwarf members of the selected common proper motion pairs, as well as high resolution spectra of their companions. From these observations, a full analysis of the two members of each common proper motion pair has lead to the initial and final masses of the white dwarfs. These observations have allowed us to provide updated information for the white dwarfs, since some of them were misclassified. This work is the first one in using common proper motion pairs to improve the initial-final mass relationship, and has also allowed to cover the poorly explored low-mass domain. As in the case of studies based on white dwarfs in open clusters, the distribution of the semi-empirical data presents a large scatter, which is higher than the expected uncertainties in the derived values. This suggests that the initial-final mass relationship may not be a single-valued function.

The interpretation of microlensing results towards the Large Magellanic Cloud (LMC) still remains controversial. White dwarfs have been proposed to explain these results and, hence, to contribute significantly to the mass budget of our Galaxy. However, several constraints on the role played by regular carbon-oxygen white dwarfs exist. Massivewhite dwarfs are thought to be made of a mixture of oxygen and neon. Correspondingly, their cooling rate is larger than those of typical carbon-oxygen white dwarfs and they fade to invisibility in short timescales. Consequently, they constitute a good candidate for explaining the microlensing results. Here, we examine in detail this hypothesis by using the most recent and up-to-date cooling tracks for massive white dwarfs and a Monte Carlo simulator which takes into account the most relevant Galactic inputs. We find that oxygen-neon white dwarfs cannot account for a substantial fraction of the microlensing depth towards the LMC, independently of the adopted initial mass function, although some microlensing events could be due to oxygen--neon white dwarfs. The white dwarf population contributes at most a 5% to the mass of the Galactic halo.

Only the 1/V_max method has been employed so far for observationally determining the white dwarf luminosity function, whereas for other kind of luminosity functions several other methods have been frequently used. Moreover, the procedures to determine the white dwarf luminosity function are not free of biases. These biases have two different origins: they can either be of statistical nature or a consequence of the measurement errors. In a previous paper we carried out an in-depth study of the first category of biases for several luminosity function estimators. In this paper we focus on the biases introduced by the measurement errors and on the effects of the degree of contamination of the input sample used to build the disc white dwarf luminosity function by different kinematical populations. To assess the extent of these biases we use a Monte Carlo simulator to generate a controlled synthetic population and analyse the behaviour of the disc white dwarf luminosity function for several assumptions about the magnitude of the measurement errors and for several degrees of contamination, comparing the performances of the most robust luminosity function estimators under such conditions.

GD 362 is a massive white dwarf with a spectrum suggesting a H-rich atmosphere which also shows very high abundances of Ca, Mg, Fe and other metals. However, for pure H-atmospheres the diffusion timescales are so short that very extreme assumptions have to be made to account for the observed abundances of metals. The most favored hypothesis is that the metals are accreted from either a dusty disk or from an asteroid belt. Here we propose that the envelope of GD 362 is dominated by He, which at these effective temperatures is almost completely invisible in the spectrum. This assumption strongly alleviates the problem, since the diffusion timescales are much larger for He-dominated atmospheres. We also propose that the He-dominated atmosphere of GD 362 is likely to be the result of the merger of a binary white dwarf.

We compute and analyze the evolution of primordial stars of masses at the ZAMS between 5 M_sun and 10 M_sun, with and without overshooting. Our main goals are to determine the nature of the remnants of massive intermediate-mass primordial stars and to check the influence of overshooting in their evolution. Our calculations cover stellar evolution from the main sequence phase until the formation of the degenerate cores and the thermally pulsing phase. We have obtained the values for the limiting masses of Population III progenitor stars leading to carbon-oxygen and oxygen-neon compact cores. Moreover, we have also obtained the limiting mass for which isolated primordial stars would lead to core-collapse supernovae after the end of the main central burning phases. Considering a moderate amount of overshooting the mass thresholds at the ZAMS for the formation of carbon-oxygen and oxygen-neon degenerate cores shifts to smaller values by about 2 M_sun. As a by-product of our calculations, we have also obtained the structure and composition profiles of the resulting compact remnants. Opposite to what happens with solar metallicity objects, the final fate of primordial stars is not straightforward determined from the mass of the compact cores at the end of carbon burning. Instead, the small mass-loss rates typically associated to stellar winds of low metallicity stars might allow the growth of the resulting degenerate cores up to the Chandrasekhar mass, on time scales one or two orders of magnitude shorter than the time required to loose the envelope. This would lead to the formation of supernovae for initial masses as small as about 5 M_sun.

Over the years, several methods have been proposed to compute galaxy luminosity functions, from the most simple ones -counting sample objects inside a given volume- to very sophisticated ones -like the C- method, the STY method or the Choloniewski method, among others. However, only the V/Vmax method is usually employed in computing the white dwarf luminosity function and other methods have not been applied so far to the observational sample of spectroscopically identified white dwarfs. Moreover, the statistical significance of the white dwarf luminosity function has also received little attention and a thorough study still remains to be done. In this paper we study, using a controlled synthetic sample of white dwarfs generated using a Monte Carlo simulator, which is the statistical significance of the white dwarf luminosity function and which are the expected biases. We also present a comparison between different estimators for computing the white dwarf luminosity function. We find that for sample sizes large enough the V/Vmax method provides a reliable characterization of the white dwarf luminosity function, provided that the input sample is selected carefully. Particularly, the V/Vmax method recovers well the position of the cut-off of the white dwarf luminosity function. However, this method turns out to be less robust than the Choloniewski method when the possible incompletenesses of the sample are taken into account. We also find that the Choloniewski method performs better than the V/Vmax method in estimating the overall density of white dwarfs, but misses the exact location of the cut-off of the white dwarf luminosity function.

We compute the emission of gravitational radiation from pulsating white dwarfs. This is done by using an up-to-date stellar evolutionary code coupled with a state-of-the-art pulsational code. The emission of gravitational waves is computed for a standard 0.6 solar masses white dwarf with a liquid carbon-oxygen core and a hydrogen-rich envelope, for a massive DA white dwarf with a partially crystallized core for which various l=2 modes have been observed (BPM 37093) and for PG 1159-035, the prototype of the GW Vir class of variable stars, for which several quadrupole modes have been observed as well. We find that these stars do not radiate sizeable amounts of gravitational waves through their observed pulsation g-modes, in line with previous studies. We also explore the possibility of detecting gravitational waves radiated by the f-mode and the p-modes. We find that in this case the gravitational wave signal is very large and, hence, the modes decay very rapidly. We also discuss the possible implications of our calculations for the detection of gravitational waves from pulsating white dwarfs within the framework of future space-borne interferometers like LISA.

One of the most promising space missions of ESA is the astrometric satellite Gaia, which will provide very precise astrometry and multicolour photometry, for all 1.3 billion objects to V~20, and radial velocities with accuracies of a few km/s for most stars brighter than V ~ 17. Consequently, full homogeneous six-dimensional phase-space information for a huge number of stars will become available. Our Monte Carlo simulator has been used to estimate the number of white dwarfs potentially observable by Gaia. From this we assess which would be the white dwarf luminosity functions which Gaia will obtain and discuss in depth the scientific returns of Gaia in the specific field of white dwarf populations. Scientific attainable goals include, among others, a reliable determination of the age of the Galactic disk, a better knowledge of the halo of the Milky Way and the reconstruction of the star formation history of the Galactic disk. Our results also demonstrate the potential impact of a mission like Gaia in the current understanding of the white dwarf cooling theory.

Gaia is an ambitious space observatory devoted to obtain the largest and most precise astrometric catalogue of astronomical objects from our Galaxy and beyond. On-board processing and transmission of the huge amount of data generated by the instruments is one of its several technological challenges. The measurement time tags are critical for the scientific results of the mission, so they must be measured and transmitted with the highest precision - leading to an important telemetry channel occupation. In this paper we present the optimisation of time data, which has resulted in a useful software tool. We also present how time data is adapted to the Packet Telemetry standard. The several communication layers are illustrated and a method for coding and transmitting the relevant data is described as well. Although our work is focused on Gaia, the timing scheme and the corresponding tools can be applied to any other instrument or mission with similar operational principles.

Gaia will perform astrometric measurements with an unprecedented resolution. Consequently, the electronics of the Astro instrument must time tag every measurement with a precision of a few nanoseconds. Hence, it requires a high stability clock signal, for which a Rb-type spacecraft master clock has been baselined. The distribution of its signal and the generation of clock subproducts must maintain these high accuracy requirements. We have developed a software application to simulate generic clock frameworks. The most critical clock structures for Gaia have also been identified, and its master clock has been parameterised.

The Payload Data Handling System (PDHS) of Gaia is a technological challenge, since it will have to process a huge amount of data with limited resources. Its main tasks include the optimal codification of science data, its packetisation and its compression, before being stored on-board ready to be transmitted. Here we describe a set of proposals for its design, as well as some simulators developed to optimise and test these proposals.

We present the results of a detailed set of one-zone models that account for the coupling between pulsation and convection following the original prescriptions of Stellingwerf (1986). Motivated by the arbitrary nature of the input parameters adopted in this theoretical framework, we computed several sequences of models that cover a substantial fraction of the parameter space and a longer integration time. We also included a turbulent pressure term and found that this physical mechanism plays a crucial role in the pulsation characteristics of the models by removing the sharp discontinuities along the light and the velocity curves showed by models that do not account for turbulent pressure. Finally, we investigated the vibrational and the pulsational stability of completely convective models. We consider as the most important finding of the present work the identification of a well-defined region in the parameter space where they approach limit-cycle stability. Several numerical experiments performed by adopting different values of the adiabatic exponent and of the shell thickness indicate that the coupling between pulsation and convection is the key driving mechanism for LPVs, a finding supported by recent theoretical predictions.

The final stages of the evolution of electron--degenerate ONe cores, resulting from carbon burning in ``heavy weight'' intermediate--mass stars ($8 M_{\sun}\la M \la 11 M_{\sun}$) and growing in mass, either from carbon burning in a shell or from accretion of matter in a close binary system, are examined in the light of their detailed chemical composition. In particular, we have modelled the evolution taking into account the abundances of the following minor nuclear species, which result from the previous evolutionary history: $^{12}$C, $^{23}$Na, $^{24}$Mg, and $^{25}$Mg. Both $^{23}$Na and $^{25}$Mg give rise to Urca processes, which are found to be unimportant for the final outcome of the evolution. $^{24}$Mg was formerly considered a major component of ONe cores (hence called ONeMg cores), but updated evolutionary calculations in this mass range have severely reduced its abundance. Nevertheless, we have parameterized it and we have found that the minimum amount of $^{24}$Mg required to produce NeO burning at moderate densities is $\sim 23%$, a value exceedingly high in the light of recent evolutionary models. Finally, we have determined that models with relatively small abundances of unburnt carbon ($X(^{12}$C)$\sim 0.015$) could be a channel to explosion at low to moderate density ($\sim 1\times 10^9$ g cm$^{-3}$). This is clearly below the current estimate for the explosion/collapse threshold and would have interesting consequences.

This paper is designed to 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 M_sun star 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 overshooting, salt fingers and convection is considered. We employed the double-diffusive mixing-length theory of convection for fluids with composition gradients (Grossman & Taam 1996). 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 PG1159 stars and DB white dwarf evolution are studied in the frame 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 PG1159-DB-DQ. In this regard, the implications for the double-layered chemical structure in pulsating DB white dwarfs is investigated. Another aspect of the investigation concerned the consequences of mass-loss episodes during the PG1159 stage for the chemical stratification of the outer layer of DB and DQ white dwarfs.

The evolution of a zero metallicity 9 M_s star is computed, analyzed and compared with that of a solar metallicity star of identical ZAMS mass. Our computations range from the main sequence until the formation of a massive oxygen-neon white dwarf. Special attention has been payed to carbon burning in conditions of partial degeneracy as well as to the subsequent thermally pulsing Super-AGB phase. The latter develops in a fashion very similar to that of a solar metallicity 9 M_s star, as a consequence of the significant enrichment in metals of the stellar envelope that ensues due to the so-called third dredge-up episode. The abundances in mass of the main isotopes in the final ONe core resulting from the evolution are X(^{16}O) approx 0.59, X(^{20}Ne) approx 0.28 and X(^{24}Mg) approx 0.05. This core is surrounded by a 0.05 M_s buffer mainly composed of carbon and oxygen, and on top of it a He envelope of mass 10^{-4} M_s

We compute the emission of gravitational radiation from the merging of a close white dwarf binary system. This is done for a wide range of masses and compositions of the white dwarfs, ranging from mergers involving two He white dwarfs, mergers in which two CO white dwarfs coalesce to mergers in which a massive ONe white dwarf is involved. In doing so we follow the evolution of binary system using a Smoothed Particle Hydrodynamics code. Even though the coalescence process of the white dwarfs involves considerable masses, moving at relatively high velocities with a high degree of asymmetry we find that the signature of the merger is not very strong. In fact, the most prominent feature of the coalescence is that in a relatively small time scale (of the order of the period of the last stable orbit, typically a few minutes) the sources stop emitting gravitational waves. We also discuss the possible implications of our calculations for the detection of the coalescence within the framework of future space-borne interferometers like LISA.

We explore the adiabatic pulsational properties of massive white dwarf stars with hydrogen-rich envelopes and oxygen/neon and carbon/oxygen cores. To this end, we compute the cooling of massive white dwarf models for both core compositions taking into account the evolutionary history of the progenitor stars and the chemical evolution caused by time-dependent element diffusion. In particular, for the oxygen/neon models, we adopt the chemical profile resulting from repeated carbon-burning shell flashes expected in very massive white dwarf progenitors. For carbon/oxygen white dwarfs we consider the chemical profiles resulting from phase separation upon crystallization. For both compositions we also take into account the effects of crystallization on the oscillation eigenmodes. We find that the pulsational properties of oxygen/neon white dwarfs are notably different from those made of carbon/oxygen, thus making asteroseismological techniques a promising way to distinguish between both types of stars and, hence, to obtain valuable information about their progenitors.

We present in this work new pulsational calculations for improved carbon-oxygen DA white dwarf models suitable for the study of massive ZZ Ceti stars. The background models employed in this study, presented in detail in a recent paper by Althaus et al. (2003), are the result of the complete evolution of massive white dwarf progenitors from the zero-age main sequence through the Asymptotic Giant Branch (AGB) and mass loss phases to the white dwarf regime. Abundance changes are accounted for by means of a full coupling between nuclear evolution and time-dependent mixing due to convection, salt fingers, and diffusive overshoot. In addition, time-dependent element diffusion for multicomponent gases has been considered during the white dwarf evolution. Crystallization and chemical rehomogenization due to phase separation upon crystallization in the core of our models have been fully considered. The effects of crystallization on the period spectrum of these massive white dwarf models are assessed by means of a detailed pulsational analysis. We find that the theoretical pulsation spectrum is strongly modified when crystallization is considered, in particular concerning the mode trapping properties of the equilibrium models. We also discuss at some length the implications of our study for BPM 37093, the most massive ZZ Ceti star presently known. We find that if BPM 37093 has a stellar mass of $\approx 1.00$ \msun its observed spectrum could bear the signature of overshoot episodes during the helium core burning.

The interpretation of microlensing results towards the Large Magellanic Cloud (LMC) still remains controversial. Whereas white dwarfs have been proposed to explain these results and, hence, to contribute significantly to the mass budget of our Galaxy, there are as well several constraints on the role played by white dwarfs. In this paper we analyze self-consistently and simultaneously four different results, namely, the local halo white dwarf luminosity function, the microlensing results reported by the MACHO team towards the LMC, the results of Hubble Deep Field (HDF) and the results of the EROS experiment, for several initial mass functions and halo ages. We find that the proposed log-normal initial mass functions do not contribute to solve the problem posed by the observed microlensing events and, moreover, they overproduce white dwarfs when compared to the results of the HDF and of the EROS survey. We also find that the contribution of hydrogen-rich white dwarfs to the dynamical mass of the halo of the Galaxy cannot be more than $\sim 4%$.

We derive the formulae for the time variation of the gravitational "constant" G and of the fine structure "constant" \alpha in various models with extra dimensions and analyze their consistency with the available observational data for distant supernovae. We find that the reported variation of \alpha translates into a small variation of G that makes distant supernovae to appear brighter, in contradiction with recent observations of high z supernovae. The significance of these results within the framework of some cosmological scenarios is also discussed. We find, however, that the magnitude of the effect is not large enough to safely discard the models with extra dimensions studied here.

We analyse the impact of the initial abundances of the underlying white dwarf core on the nucleosynthesis accompanying classical nova outbursts, in the framework of hydrodynamic models of the explosion. Specifically, we take into account the chemical stratification of the white dwarf. It turns out that the presence of a thick CO-buffer on top of the ONe-rich core, as results from detailed calculations of previous evolution of the progenitor star, may lead to significant amounts of both 7Li and 26Al, after an outburst that, due to the lack of neon isotopes in the ejecta, would be misclassified as a non-neon nova (i.e., CO nova).

In this paper, we revisit the problem of the determination of the frequency of occurrence of galactic nova outbursts which involve an oxygen-neon (ONe) white dwarf. The improvement with respect to previous work on the subject derives from the fact that we use the results that our evolutionary calculations provide for the final mass and for the chemical profiles of intermediate-to-massive primary components of close binary systems. In particular, the final evolutionary stages, such as the carbon burning phase, have been carefully followed for the whole range of masses of interest. The chemical profiles obtained with our evolutionary code are of interest in determining the chemical composition of the ejecta after being processed through the thermonuclear runaway, although such other factors as the efficiency of the mixing between the accreted material and that of the underlying white dwarf must also be considered. In our calculations of the frequency of occurrence of nova outbursts involving an ONe white dwarf, we also take into account the observational selection effects introduced by the different recurrence times of the outbursts and by the spatial distribution of novae. In spite of the very different evolutionary sequences, we find that approximately 1/3 of the novae observed in outburst should involve an oxygen-neon white dwarf, in agreement with previous theoretical estimates.

There is growing observational evidence that the irregular changes in the light curves of certain variable stars might be due to deterministic chaos. Supporting these conclusions, several simple models of non-linear oscillators have been shown to be capable of reproducing the observed complex behaviour. In this work, we introduce a non-linear, non-adiabatic one-zone model intended to reveal the factors leading to irregular luminosity variations in some pulsating stars. We have studied and characterized the dynamical behaviour of the oscillator as the input parameters are varied. The parametric study implied values corresponding to stellar models in the family of Long Period Variables and in particular of Mira-type stars. We draw the attention on certain solutions that reproduce with reasonable accuracy the observed behaviour of some peculiar Mira variables.

The evolution of a star of initial mass 9 M_s, and Z = 0.02 in a Close Binary System is followed in the presence of different mass companions in order to study their influence on the final evolutionary stages and, in particular, on the structure and composition of the remnant components. We study two extreme cases. In the first one the mass of the secondary is 8 M_s, whereas in the second one the mass was assumed to be 1 M_s. For the first of those cases we have also explored the possible outcomes of both conservative and non-conservative mass-loss episodes. During the first mass transfer episode, several differences arise between the models. The system with the more extreme mass ratio is not able to survive the 1st. Roche lobe overflow, and spiral-in of the secondary onto the envelope of the primary is most likely. The system formed by two stars of comparable mass undergoes two mass transfer episodes in which the primary is the donor. We have performed two sets of calculations corresponding to this case in order to account for conservative and non-conservative mass transfer during the first mass loss episode. One of our main results is that for the non-conservative case the secondary becomes a Super-AGB. Such a star undergoes a final dredge-up episode, similar to that of a single star of comparable mass. The primary components do not undergo a Super-AGB phase, but instead a carbon-oxygen white dwarf is formed in both cases, before reversal mass transfer occurs. However, given the extreme mass ratios at this stage between the components of the binary system, the possibility of merger episodes remains likely. We also discuss the presumable final outcomes of the system and possible observational counterparts.

The interpretation of the old, cool white dwarfs recently found by Oppenheimer et al. (2001) is still controversial. Whereas these authors claim that they have finally found the elusive ancient halo white dwarf population that contributes significantly to the mass budget of the galactic halo, there have been several other contributions that argue that these white dwarfs are not genuine halo members but, instead, thick disk stars. We show here that the interpretation of this sample is based on the adopted distances, which are obtained from a color--magnitude calibration, and we demonstrate that when the correct distances are used a sizeable fraction of these putative halo white dwarfs belong indeed to the disk population. We also perform a maximum likelihood analysis of the remaining set of white dwarfs and we find that they most likely belong to the thick disk population. However, another possible explanation is that this sample of white dwarfs has been drawn from a 1:1 mixture of the halo and disk white dwarf populations.

When intermediate mass stars reach their last stages of evolution they show pronounced oscillations. This phenomenon happens when these stars reach the so-called Asymptotic Giant Branch (AGB), which is a region of the Hertzsprung-Russell diagram located at about the same region of effective temperatures but at larger luminosities than those of regular giant stars. The period of these oscillations depends on the mass of the star. There is growing evidence that these oscillations are highly correlated with mass loss and that, as the mass loss increases, the pulsations become more chaotic. In this paper we study a simple oscillator which accounts for the observed properties of this kind of stars. This oscillator was first proposed and studied by Icke et al. [Astron.Astrophys. 258, 341 (1992)] and we extend their study to the region of more massive and luminous stars - the region of Super-AGB stars. The oscillator consists of a periodic nonlinear perturbation of a linear Hamiltonian system. The formalism of dynamical systems theory has been used to explore the associated Poincare map for the range of parameters typical of those stars. We have studied and characterized the dynamical behaviour of the oscillator as the parameters of the model are varied, leading us to explore a sequence of local and global bifurcations. Among these, a tripling bifurcation is remarkable, which allows us to show that the Poincare map is a nontwist area preserving map. Meandering curves, hierarchical-islands traps and sticky orbits also show up. We discuss the implications of the stickiness phenomenon in the evolution and stability of the Super-AGB stars.

Recent high-redshift Type Ia supernovae results can be used to set new bounds on a possible variation of the gravitational constant $G$. If the local value of $G$ at the space-time location of distant supernovae is different, it would change both the kinetic energy release and the amount of $^{56}$Ni synthesized in the supernova outburst. Both effects are related to a change in the Chandrasekhar mass $M_{Ch} \propto G^{-3/2}$. In addition, the integrated variation of $G$ with time would also affect the cosmic evolution and therefore the luminosity distance relation. We show that the later effect in the magnitudes of Type Ia supernovae is typically several times smaller than the change produced by the corresponding variation of the Chandrasekhar mass. We investigate in a consistent way how a varying $G$ could modify the Hubble diagram of Type Ia supernovae and how these results can be used to set upper bounds to a hypothetical variation of $G$. We find $G/G_0 \la 1.1$ and $G'/G \la 10^{-11} yr^{-1}$ at redshifts $z\simeq 0.5$. These new bounds extend the currently available constrains on the evolution of $G$ all the way from solar and stellar distances to typical scales of Gpc/Gyr, i.e. by more than 15 orders of magnitudes in time and distance.

The evolution of a star of initial mass 10 $M_{\odot}$, and metallicity $Z = 0.02$ in a Close Binary System (CBS) is followed from its main sequence until an ONe degenerate remnant forms. Restrictions have been made on the characteristics of the companion as well as on the initial orbital parameters in order to avoid the occurrence of reversal mass transfer before carbon is ignited in the core. The system undergoes three mass loss episodes. The first and second ones are a consequence of a case B Roche lobe overflow. During the third mass loss episode stellar winds may play a role comparable to, or even more important than Roche lobe overflow. In this paper, we extend the previously existing calculations of stars of intermediate mass belonging to close binary systems by following carefully the carbon burning phase of the primary component. We also propose different possible outcomes for our scenario and discuss the relevance of our findings. In particular, our main result is that the resulting white dwarf component of mass $1.1 M_\odot$ more likely has a core composed of oxygen and neon, surrounded by a mantle of carbon-oxygen rich material. The average abundances of the oxygen-neon rich core are $X({\rm O}^{16})=0.55$, $X({\rm Ne}^{20})=0.28$, $X({\rm Na}^{23})=0.06$ and $X({\rm Mg}^{24})=0.05$. This result has important consequences for the Accretion Induced Collapse scenario. The average abundances of the carbon-oxygen rich mantle are $X({\rm O}^{16})=0.55$, and $X({\rm C}^{12})=0.43$. The existence of this mantle could also play a significant role in our understanding of cataclysmic variables.

White dwarfs are well studied objects. The relative simplicity of their physics allows to obtain very detailed models which can be ultimately compared with their observed properties. Among white dwarfs there is a specific class of stars, known as ZZ-Ceti objects, which have a hydrogen-rich envelope and show periodic variations in their light curves. G117-B15A belongs to this particular set of stars. The luminosity variations have been successfully explained as due to g-mode pulsations. G117-B15A has been recently claimed to be the most stable optical clock ever found, being the rate of change of its 215.2 s period very small: \dot{P}= (2.3 +- 1.4)x10^{-15} s s^-1, with a stability comparable to that of the most stable millisecond pulsars. The rate of change of the period is closely related to its cooling timescale, which can be accurately computed. In this paper we study the pulsational properties of G117-B15A and we use the observed rate of change of the period to impose constraints on the axion emissivity and, thus, to obtain a preliminary upper bound to the mass of the axion. This upper bound turns out to be 4cos^{2}{\beta} meV at the 95% confidence level. Although there are still several observational and theoretical uncertainties, we conclude that G117-B15A is a very promising stellar object to set up constraints on particle physics.

We analyze in detail various possible sources of systematic errors on the distances of globular clusters derived by fitting a local template DA white dwarf sequence to the cluster counterpart (the so-called WD-fitting technique). We find that the unknown thickness of the hydrogen layer of white dwarfs in clusters plays a non negligible role. For reasonable assumptions - supported by the few sparse available observational constraints - about the unknown mass and thickness of the hydrogen layer for the cluster white dwarfs, a realistic estimate of the systematic error on the distance is within +-0.10 mag. However, particular combinations of white dwarf masses and envelope thicknesses - which at present cannot be excluded a priori - could produce larger errors. Contamination of the cluster DA sequence by non-DA white dwarfs introduces a very small systematic error of about -0.03 mag in the Mv/(V-I) plane, but in the Mv/(B-V) plane the systematic error amounts to ~ +0.20 mag. Contamination by white dwarfs with helium cores should not influence appreciably the WD-fitting distances. Finally, we obtain a derivative D((m-M)v)/D(E(B-V))~ -5.5 for the WD-fitting distances, which is very similar to the dependence found when using the Main Sequence fitting technique.

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/Lo)=-4.5 for an otherwise typical 0.6 Mo 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.

There are at least three ways in which a varying gravitational constant $G$ could affect the interpretation of the recent high-redhisft Type Ia supernovae results. If the local value of $G$ at the space-time location of distant supernovae is different, it would change both the thermonuclear energy release and the time scale of the supernova outburst. In both cases the effect is related to a change in the Chandrasekhar mass $M_{\rm Ch}\propto G^{-3/2}$. Moreover the integrated variation of $G$ with time would also affect cosmic evolution and therefore the luminosity distance relation. Here we investigate in a consistent way how these different effects of a varying $G$ could change the current interpretation of the Hubble diagram of Type Ia supernovae. We parametrize the variation of $G$ using scalar-tensor theories of gravity, such as the Jordan-Brans-Dicke theory or its extensions. It is remarkable that Dirac's hypothesis that $G$ should decrease with time can qualitatively explain the observed $\Delta m \simeq 0.2$ mag decrease at $z\simeq 0.5$ (with respect to a decelerating universe) and, at the same time, reduce the duration of the risetimes of distant Type Ia supernovae as recently reported.

The evolution of white dwarfs is a cooling process that depends on the energy stored in the core and on the way in which it is transferred through the envelope. In this paper we show that despite some (erroneous) claims, the redistribution of chemical elements ensuing the crystallization of C/O white dwarfs provides between the 10% and the 20% of the total energy released during the crystallization process, depending on the internal chemical composition, which is not negligible at all, given the present state of the art of the white dwarf cooling theory.

The white dwarf luminosity function has proven to be an excellent tool to study some properties of the galactic disk such as its age and the past history of the local star formation rate. The existence of an observational luminosity function for halo white dwarfs could provide valuable information about its age, the time that the star formation rate lasted, and could also constrain the shape of the allowed Initial Mass Functions (IMF). However, the main problem is the scarce number of white dwarfs already identified as halo stars. In this Letter we show how an artificial intelligence algorithm can be succesfully used to classify the population of spectroscopically identified white dwarfs allowing us to identify several potential halo white dwarfs and to improve the significance of its luminosity function.

In order to understand the dynamical and chemical evolution of our Galaxy it is of fundamental importance to study the local neighborhood. White dwarf stars are ideal candidates to probe the history of the solar neighborhood, since these ``fossil'' stars have very long evolutionary time-scales and, at the same time, their evolution is relatively well understood. In fact, the white dwarf luminosity function has been used for this purpose by several authors. However, a long standing problem arises from the relatively poor statistics of the samples, especially at low luminosities. In this paper we assess the statistical reliability of the white dwarf luminosity function by using a Monte Carlo approach.

Halo white dwarfs can provide important information about the properties and evolution of the galactic halo. In this paper we compute, assuming a standard IMF and updated models of white dwarf cooling, the expected luminosity function, both in luminosity and in visual magnitude, for different star formation rates. We show that a deep enough survey (limiting magnitude > 20) could provide important information about the halo age and the duration of the formation stage. We also show that the number of white dwarfs produced using the recently proposed biased IMFs cannot represent a large fraction of the halo dark matter if they are constrained by the presently observed luminosity function. Furthermore, we show that a robust determination of the bright portion of the luminosity function can provide strong constraints on the allowable IMF shapes.

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_{\sun}$ 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. (1988) and Oswalt et al. (1996) gives and 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 Gyrs when it is neglected.