Showing 1–50 of 141 results for author: García-Berro, E

3D Hydrodynamical Simulations of Helium-Ignited Double-degenerate White Dwarf Mergers

Authors: Niranjan Roy, Vishal Tiwari, Alexey Bobrick, Daniel Kosakowski, Robert Fisher, Hagai B. Perets, Rahul Kashyap, Pablo Lorén-Aguilar, Enrique García-Berro

Abstract: The origins of type Ia supernovae (SNe Ia) are still debated. Some of the leading scenarios involve a double detonation in double white dwarf (WD) systems. In these scenarios, helium shell detonation occurs on top of a carbon-oxygen (CO) WD, which then drives the detonation of the CO-core, producing a SN Ia. Extensive studies have been done on the possibility of a double helium detonation, followi… ▽ More The origins of type Ia supernovae (SNe Ia) are still debated. Some of the leading scenarios involve a double detonation in double white dwarf (WD) systems. In these scenarios, helium shell detonation occurs on top of a carbon-oxygen (CO) WD, which then drives the detonation of the CO-core, producing a SN Ia. Extensive studies have been done on the possibility of a double helium detonation, following a dynamical helium mass-transfer phase onto a CO-WD. However, 3D self-consistent modeling of the double-WD system, the mass transfer, and the helium shell detonation have been little studied. Here we use 3D hydrodynamical simulations to explore this case in which a helium detonation occurs near the point of Roche lobe overflow of the donor WD and may lead to an SN Ia through the dynamically driven double-degenerate double-detonation (D6) mechanism. We find that the helium layer of the accreting primary WD does undergo a detonation, while the underlying carbon-oxygen core does not, leading to an extremely rapid and faint nova-like transient instead of a luminous SN Ia event. This failed core detonation suggests that D6 SNe Ia may be restricted to the most massive carbon-oxygen primary WDs. We highlight the nucleosynthesis of the long-lived radioisotope $^{44}$Ti during explosive helium burning, which may serve as a hallmark both of successful as well as failed D6 events which subsequently detonate as classical double-degenerate mergers. △ Less

Submitted 25 July, 2022; v1 submitted 20 April, 2022; originally announced April 2022.

Comments: 11 pages, 4 figures, 1 table. Published in the Astrophysical Journal Letters

Journal ref: 2022 ApJL 932 L24

The effects of $^{22}$Ne sedimentation and metallicity on the local 40 pc white dwarf luminosity function

Authors: Jordi Tononi, Santiago Torres, Enrique García-Berro, María E. Camisassa, Leandro G. Althaus, Alberto Rebassa-Mansergas

Abstract: 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 t… ▽ More 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. △ Less

Submitted 19 June, 2019; originally announced June 2019.

Comments: 11 pages, 9 figures and 2 tables; accepted for publication in A&A

Journal ref: A&A 628, A52 (2019)

Double-Degenerate Carbon-Oxygen and Oxygen-Neon White Dwarf Mergers: A New Mechanism for Faint and Rapid Type Ia Supernovae

Authors: Rahul Kashyap, Tazkera Haque, Pablo Lorén-Aguilar, Enrique García-Berro, Robert T. Fisher

Abstract: 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… ▽ More 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. △ Less

Submitted 31 October, 2018; originally announced November 2018.

Comments: 13 pages, 8 figures, accepted to The Astrophysical Journal. Dedicated to the memory of our late colleague Enrique García-Berro

Axions and the luminosity function of white dwarfs. The thin and thick disks, and the halo

Authors: Jordi Isern, Enrique Garcia-Berro, Santiago Torres, Roxana Cojocaru, Silvia Catalan

Abstract: 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 dw… ▽ More 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. △ Less

Submitted 30 April, 2018; originally announced May 2018.

Comments: 8 pages, 6 figures, accepted for publication in MNRAS

Spatial variation of fundamental constants: testing models with thermonuclear supernovae

Authors: Carolina Negrelli, Lucila Kraiselburd, Susana J. Landau, Enrique García-Berro

Abstract: 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 α 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… ▽ More 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 α 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 α 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 $μ_0$ and the cosmological density matter $Ω_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. △ Less

Submitted 4 April, 2018; originally announced April 2018.

Comments: 20 pages, 3 figures. Accepted for publication in IJMPD

Journal ref: International Journal of Modern Physics D, Vol. 27, No. 09, 1850099 (2018)

The Spin Evolution of Fast-Rotating, Magnetized Super-Chandrasekhar White Dwarfs in the Aftermath of White Dwarf Mergers

Authors: L. Becerra, J. A. Rueda, P. Loren-Aguilar, E. Garcia-Berro

Abstract: 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… ▽ More 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. △ Less

Submitted 4 April, 2018; originally announced April 2018.

Comments: Accepted for publication in The Astrophysical Journal

The population of single and binary white dwarfs of the Galactic bulge

Authors: S. Torres, E. García-Berro, R. Cojocaru, A. Calamida

Abstract: 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 w… ▽ More 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. △ Less

Submitted 1 February, 2018; originally announced February 2018.

Comments: 10 pages, 9 figures, accepted for publication in MNRAS

Three-dimensional simulations of the interaction between the nova ejecta, the accretion disk, and the companion star

Authors: Joana Figueira, Jordi Jose, Enrique Garcia-Berro, Simon W. Campbell, Domingo Garcia-Senz, Shazrene Mohamed

Abstract: 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 temperat… ▽ More 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. △ Less

Submitted 22 December, 2017; originally announced December 2017.

Comments: Accepted for Publication in Astronomy and Astrophysics

Journal ref: A&A 613, A8 (2018)

The white dwarf binary pathways survey - II. Radial velocities of 1453 FGK stars with white dwarf companions from LAMOST DR4

Authors: A. Rebassa-Mansergas, J. J. Ren, P. Irawati, E. Garcia-Berro, S. G. Parsons, M. R. Schreiber, B. T. Gaensicke, P. Rodriguez-Gil, X. Liu, C. Manser, S. Palomo-Nevado, F. Jimenez-Ibarra, R. Costero, J. Echevarria, R. Michel, M. Zorotovic, M. Hollands, Z. Han, A. Luo, E. Villaver, X. Kong

Abstract: 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 spect… ▽ More 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. △ Less

Submitted 30 August, 2017; originally announced August 2017.

Comments: accepted for publication by MNRAS

The population of white dwarf-main sequence binaries in the SDSS DR 12

Authors: R. Cojocaru, A. Rebassa-Mansergas, S. Torres, E. Garcia-Berro

Abstract: 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 perfor… ▽ More 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). △ Less

Submitted 16 May, 2017; originally announced May 2017.

Comments: accepted for publication by MNRAS

One-Armed Spiral Instability in Double-Degenerate Post-Merger Accretion Disks

Authors: Rahul Kashyap, Robert Fisher, Enrique García-Berro, Gabriela Aznar-Siguán, Suoqing Ji, Pablo Lorén-Aguilar

Abstract: 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… ▽ More 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. △ Less

Submitted 12 April, 2017; v1 submitted 5 April, 2017; originally announced April 2017.

Comments: Accepted to Astrophysical Journal

White dwarf dynamical interactions and fast optical transients

Authors: Enrique García-Berro, Carles Badenes, Gabriela Aznar-Siguán, Pablo Lorén-Aguilar

Abstract: 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 simulat… ▽ More 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. △ Less

Submitted 29 March, 2017; originally announced March 2017.

Comments: Accepted for publication in MNRAS

The kinematics of the white dwarf population from the SDSS DR12

Authors: Borja Anguiano, Alberto Rebassa-Mansergas, Enrique Garcia-Berro, Santiago Torres, Ken Freeman, Tomaz Zwitter

Abstract: 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… ▽ More 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. △ Less

Submitted 27 March, 2017; originally announced March 2017.

Comments: Accepted for publication by MNRAS

Evolutionary sequences for hydrogen-deficient white dwarfs

Authors: María E. Camisassa, Leandro G. Althaus, René D. Rohrmann, Enrique García-Berro, Santiago Torres, Alejandro H. Córsico, Felipe C. Wachlin

Abstract: 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.… ▽ More 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. △ Less

Submitted 15 March, 2017; originally announced March 2017.

Comments: 40 pages, 13 figures. To be published in ApJ

A common origin of magnetism from planets to white dwarfs

Authors: Jordi Isern, Enrique García-Berro, Baybars Külebi, Pablo Lorén-Aguilar

Abstract: 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 oth… ▽ More 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. △ Less

Submitted 21 February, 2017; v1 submitted 6 February, 2017; originally announced February 2017.

Comments: The Astrophysical Journal Letters, 836:L28 (5pp), 2017 February 20

Orbital periods and component masses of three double white dwarfs

Authors: A. Rebassa-Mansergas, S. G. Parsons, E. Garcia-Berro, B. T. Gaensicke, M. R. Schreiber, M. Rybicka, D. Koester

Abstract: 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… ▽ More 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. △ Less

Submitted 1 December, 2016; originally announced December 2016.

Comments: Accepted for publication by MNRAS

The age-metallicity relation in the solar neighbourhood from a pilot sample of white dwarf-main sequence binaries

Authors: A. Rebassa-Mansergas, B. Anguiano, E. García-Berro, K. C. Freeman, R. Cojocaru, C. J. Manser, A. F. Pala, B. T. Gänsicke, X. -W. Liu

Abstract: 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 campa… ▽ More 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). △ Less

Submitted 10 August, 2016; originally announced August 2016.

Comments: Accepted for publication by MNRAS

The white dwarf luminosity function

Authors: Enrique García-Berro, Terry D. Oswalt

Abstract: 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 s… ▽ More 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. △ Less

Submitted 8 August, 2016; originally announced August 2016.

Comments: Accepted for publication in New Astronomy Reviews, 3 August 2016

Three-dimensional simulations of turbulent convective mixing in ONe and CO classical nova explosions

Authors: Jordi Casanova, Jordi Jose, Enrique Garcia-Berro, Steven N. Shore

Abstract: 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, enric… ▽ More 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. △ Less

Submitted 28 June, 2016; originally announced June 2016.

Comments: Accepted for publication in Astronomy & Astrophysics

An asteroseismic constraint on the mass of the axion from the period drift of the pulsating DA white dwarf star L19-2

Authors: Alejandro H. Córsico, Alejandra D. Romero, Leandro G. Althaus, Enrique García-Berro, Jordi Isern, S. O. Kepler, Marcelo M. Miller Bertolami, Denis J. Sullivan, Paul Chote

Abstract: 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 ($Π\sim 113$ s and $Π\sim 192$ s), to derive a new constraint on the mass of the axion,… ▽ More 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 ($Π\sim 113$ s and $Π\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 β\lesssim 25$ meV or an axion-electron coupling constant of $g_{\rm ae} \lesssim 7 \times 10^{-13}$. △ Less

Submitted 12 July, 2016; v1 submitted 20 May, 2016; originally announced May 2016.

Comments: 20 pages, 6 figures, 3 tables. Revised version. Accepted for publication in the Journal of Cosmology and Astroparticle Physics. arXiv admin note: text overlap with arXiv:1205.6180

Interpolation of the magnetic field at the test masses in eLISA

Authors: I. Mateos, M. Díaz-Aguiló, J. Ramos-Castro, E. García-Berro, A. Lobo

Abstract: 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… ▽ More 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. △ Less

Submitted 28 April, 2016; originally announced April 2016.

Comments: 14 pages, 8 figures

Journal ref: Classical and Quantum Gravity, 32, 165003, 2015

The SDSS spectroscopic catalogue of white dwarf-main sequence binaries: new identifications from DR9-12

Authors: A. Rebassa-Mansergas, J. J. Ren, S. G. Parsons, B. T. Gaensicke, M. R. Schreiber, E. Garcia-Berro, X. -W. Liu, D. Koester

Abstract: 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… ▽ More 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. △ Less

Submitted 3 March, 2016; originally announced March 2016.

Comments: 14 pages, 11 figures, accepted for publication by MNRAS

The white dwarf population within 40 pc of the Sun

Authors: S. Torres, E. García-Berro

Abstract: 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… ▽ More 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. △ Less

Submitted 8 February, 2016; originally announced February 2016.

Comments: 8 pages, 7 figures, accepted for publication in A&A

Journal ref: A&A 588, A35 (2016)

Magnetic White Dwarfs: Observations, Theory, and Future Prospects

Authors: Enrique García-Berro, Mukremin Kilic, S. O. Kepler

Abstract: 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… ▽ More 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. △ Less

Submitted 29 October, 2015; originally announced October 2015.

Comments: 19 pages, 1 figure

Light Curves and Spectra from a Thermonuclear Explosion of a White Dwarf Merger

Authors: Daniel R. van Rossum, Rahul Kashyap, Robert Fisher, Ryan T. Wollaeger, Enrique Garcia-Berro, Gabriela Aznar-Siguan, Suoqing Ji, Pablo Loren-Aguilar

Abstract: 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-… ▽ More 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. △ Less

Submitted 16 August, 2016; v1 submitted 14 October, 2015; originally announced October 2015.

Comments: in ApJ 2016

Report number: LANL Report Number: LA-UR-15-27535

The white dwarf population of NGC 6397

Authors: S. Torres, E. García-Berro, L. G. Althaus, M. E. Camisassa

Abstract: 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 i… ▽ More 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. △ Less

Submitted 31 July, 2015; originally announced July 2015.

Comments: 9 pages, 6 figures, accepted for publication in A&A

Revisiting the luminosity function of single halo white dwarfs

Authors: R. Cojocaru, S. Torres, L. G. Althaus, J. Isern, E. García-Berro

Abstract: 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 availabl… ▽ More 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. △ Less

Submitted 22 July, 2015; originally announced July 2015.

Comments: 9 pages, 9 figures, accepted for publication in A&A

The mass function of hydrogen-rich white dwarfs: robust observational evidence for a distinctive high-mass excess near 1Msun

Authors: A. Rebassa-Mansergas, M. Rybicka, X. -W. Liu, Z. Han, E. Garcia-Berro

Abstract: 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 thi… ▽ More 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. △ Less

Submitted 22 June, 2015; originally announced June 2015.

Comments: Accepted for publication by MNRAS

Stellar archaeology with Gaia: the Galactic white dwarf population

Authors: Boris Gaensicke, Pier-Emmanuel Tremblay, Martin Barstow, Giuseppe Bono, Matt Burleigh, Sarah Casewell, Vik Dhillon, Jay Farihi, Enrique Garcia-Berro, Stephan Geier, Nicola Gentile-Fusillo, JJ Hermes, Mark Hollands, Alina Istrate, Stefan Jordan, Christian Knigge, Christopher Manser, Tom Marsh, Gijs Nelemans, Anna Pala, Roberto Raddi, Thomas Tauris, Odette Toloza, Dimitri Veras, Klaus Werner , et al. (1 additional authors not shown)

Abstract: 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. 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. △ Less

Submitted 8 June, 2015; originally announced June 2015.

Comments: Summary of a talk at the 'Multi-Object Spectroscopy in the Next Decade' conference in La Palma, March 2015, to be published in ASP Conference Series (editors Ian Skillen & Scott Trager)

DA white dwarfs from the LSS-GAC survey DR1: the preliminary luminosity and mass functions and formation rate

Authors: A. Rebassa-Mansergas, X. -W. Liu, R. Cojocaru, H. -B. Yuan, S. Torres, E. Garcia-Berro, M. -X. Xiang, Y. Huang, D. Koester, Y. Hou, G. Li, Y. Zhang

Abstract: 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 ant… ▽ More 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. △ Less

Submitted 18 March, 2015; originally announced March 2015.

Comments: 23 pages, 14 figures and 5 tables. Accepted for publication by MNRAS

Smoothed Particle Hydrodynamics simulations of the core-degenerate scenario for Type Ia supernovae

Authors: G. Aznar-Siguán, E. García-Berro, P. Lorén-Aguilar, N. Soker, A. Kashi

Abstract: 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 wi… ▽ More 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. △ Less

Submitted 15 April, 2015; v1 submitted 9 March, 2015; originally announced March 2015.

Comments: 17 pages, 10 figures. Accepted for publication in MNRAS

Is the central binary system of the planetary nebula Henize 2-428 a Type Ia supernova progenitor?

Authors: Enrique Garcia-Berro, Noam Soker, Leandro G. Althaus, Ignasi Ribas, Juan C. Morales

Abstract: 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 assu… ▽ More 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. △ Less

Submitted 25 October, 2015; v1 submitted 5 March, 2015; originally announced March 2015.

Comments: Accepted for publication in New Astronomy

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

Authors: Leandro G. Althaus, María E. Camisassa, Marcelo M. Miller Bertolami, Alejandro H. Córsico, Enrique García-Berro

Abstract: 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 e… ▽ More 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. △ Less

Submitted 12 February, 2015; originally announced February 2015.

Comments: To be published in Astronomy and Astrophysics. 12 pages, 11 figures

Spiral Disk Instability Can Drive Thermonuclear Explosions in Binary White Dwarf Mergers

Authors: Rahul Kashyap, Robert Fisher, Enrique García-Berro, Gabriela Aznar-Siguán, Suoqing Ji, Pablo Lorén-Aguilar

Abstract: 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 I… ▽ More 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. △ Less

Submitted 30 January, 2015; v1 submitted 22 January, 2015; originally announced January 2015.

Comments: Final version (as in ApJL) with minor edits

The variation of the fine structure constant: testing the dipole model with thermonuclear supernovae

Authors: Lucila Kraiselburd, Susana J. Landau, Carolina Negrelli, Enrique García-Berro

Abstract: 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… ▽ More 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 α, 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 α is adopted, but we also take into account the variation of the peak luminosity of Type Ia supernovae resulting from a variation of α. 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 α 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 α. △ Less

Submitted 11 December, 2014; v1 submitted 10 December, 2014; originally announced December 2014.

Comments: 8 pages, 4 figures. Accepted for publication in Ap&SS

Journal ref: Astrophysics and Space Science 357 (2015) 4

The response of a helium white dwarf to an exploding type Ia supernova

Authors: Oded Papish, Noam Soker, Enrique García-Berro, Gabriela Aznar-Siguán

Abstract: 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.0… ▽ More 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. △ Less

Submitted 16 February, 2015; v1 submitted 5 October, 2014; originally announced October 2014.

Comments: Final version. Accepted for publication in mnras

The white dwarf cooling sequence of 47 Tucanae

Authors: Enrique García-Berro, Santiago Torres, Leandro G. Althaus, Marcelo M. Miller Bertolami

Abstract: 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… ▽ More 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. △ Less

Submitted 2 October, 2014; originally announced October 2014.

Comments: 6 pages and 4 figures accepted for publication in A & A

The Large Observatory For x-ray Timing

Authors: M. Feroci, J. W. den Herder, E. Bozzo, D. Barret, S. Brandt, M. Hernanz, M. van der Klis, M. Pohl, A. Santangelo, L. Stella, A. Watts, J. Wilms, S. Zane, M. Ahangarianabhari, C. Albertus, M. Alford, A. Alpar, D. Altamirano, L. Alvarez, L. Amati, C. Amoros, N. Andersson, A. Antonelli, A. Argan, R. Artigue , et al. (320 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 29 August, 2014; v1 submitted 27 August, 2014; originally announced August 2014.

Comments: Proc. SPIE 9144, Space Telescopes and Instrumentation 2014: Ultraviolet to Gamma Ray, 91442T

Constraining the double-degenerate scenario for Type Ia supernovae from merger ejected matter

Authors: Naveh Levanon, Noam Soker, Enrique García-Berro

Abstract: 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.,… ▽ More 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. △ Less

Submitted 7 December, 2014; v1 submitted 6 August, 2014; originally announced August 2014.

Comments: 7 pages. Accepted to MNRAS

White paper: Gaia and the end states of stellar evolution

Authors: M. A Barstow, S. L. Casewell, S. Catalan, C. Copperwheat, B. Gaensicke, E. Garcia-Berro, N. Hambly, U. Heber, J. Holberg, J. Isern, S. Jeffery, S. Jordan, K. A. Lawrie, A. E. Lynas-Gray, T. Maccarone, T. Marsh, S. Parsons, R. Silvotti, J. Subasavage, S. Torres, P. Wheatley

Abstract: 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. 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. △ Less

Submitted 15 August, 2014; v1 submitted 23 July, 2014; originally announced July 2014.

Monte Carlo simulations of post-common-envelope white dwarf + main sequence binaries: The effects of including recombination energy

Authors: M. Zorotovic, M. R. Schreiber, E. García-Berro, J. Camacho, S. Torres, A. Rebassa-Mansergas, B. T. Gänsicke

Abstract: 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… ▽ More 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. △ Less

Submitted 11 July, 2014; originally announced July 2014.

Comments: 11 pages, 10 figures, accepted for publication in A&A

Journal ref: A&A 568, A68 (2014)

On the possible observational signatures of white dwarf dynamical interactions

Authors: G. Aznar-Siguán, E. García-Berro, M. Magnien, P. Lorén-Aguilar

Abstract: 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,… ▽ More 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. △ Less

Submitted 1 July, 2014; originally announced July 2014.

Comments: 14 pages, 11 figures. Accepted for publication in MNRAS

Constraining the neutrino magnetic dipole moment from white dwarf pulsations

Authors: Alejandro H. Córsico, Leandro G. Althaus, Marcelo M. Miller Bertolami, S. O. Kepler, Enrique García-Berro

Abstract: 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 boun… ▽ More 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. △ Less

Submitted 25 July, 2014; v1 submitted 23 June, 2014; originally announced June 2014.

Comments: 18 pages, 10 figures, 3 tables. Accepted for publication in Journal of Cosmology and Astroparticle Physics

The effects of metallicity on the Galactic disk population of white dwarfs

Authors: Ruxandra Cojocaru, Santiago Torres, Jordi Isern, Enrique García-Berro

Abstract: 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 populati… ▽ More 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. △ Less

Submitted 29 April, 2014; originally announced April 2014.

Comments: 7 pages, 5 figures, accepted for publication in A&A

Monte Carlo simulations of post-common-envelope white dwarf + main sequence binaries: comparison with the SDSS DR7 observed sample

Authors: J. Camacho, S. Torres, E. García-Berro, M. Zorotovic, M. R. Schreiber, A. Rebassa-Mansergas, A. Nebot Gómez-Morán, B. T. Gänsicke

Abstract: 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 observatio… ▽ More 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. △ Less

Submitted 22 April, 2014; originally announced April 2014.

Comments: 15 pages, 7 figures, accepted for publication in A&A

A population synthesis study of the luminosity function of hot white dwarfs

Authors: S. Torres, E. García-Berro, J. Krzesinski, S. J. Kleinman

Abstract: 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 obs… ▽ More 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. △ Less

Submitted 10 February, 2014; originally announced February 2014.

Comments: 6 pages, 5 figures, accepted for publication in A&A

The explosion of supernova 2011fe in the frame of the core-degenerate scenario

Authors: Noam Soker, Enrique Garcia-Berro, Leandro G. Althaus

Abstract: 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 expl… ▽ More 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. △ Less

Submitted 10 October, 2013; v1 submitted 2 September, 2013; originally announced September 2013.

Comments: Accepted by MNRAS Letters

White dwarf constraints on a varying $G$

Authors: Enrique García-Berro, Santiago Torres, Leandro G. Althaus, Alejandro H. Córsico, Pablo Lorén-Aguilar, Alejandra D. Romero, Jordi Isern

Abstract: 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. Specifica… ▽ More 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. △ Less

Submitted 25 August, 2013; originally announced August 2013.

Comments: 6 pages, 3 figures. To be published in the proceedings of the conference "Varying fundamental constants and dynamical dark energy" (8 - 13 July 2013, Sexten Center for Astrophysics)

Quiescent nuclear burning in low-metallicity white dwarfs

Authors: Marcelo M. Miller Bertolami, Leandro G. Althaus, Enrique Garcia-Berro

Abstract: 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 prog… ▽ More 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. △ Less

Submitted 9 August, 2013; originally announced August 2013.

Comments: 4 pages, 3 figures. Accepted for publication in ApJ Letters

Detonations in white dwarf dynamical interactions

Authors: Gabriela Aznar-Siguán, Enrique García-Berro, Pablo Lorén-Aguilar, Jordi José, Jordi Isern

Abstract: 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… ▽ More 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. △ Less

Submitted 27 June, 2013; originally announced June 2013.

Comments: 18 pages, 6 figures. Accepted for publication in MNRAS