Search

Recent observations of volume-limited samples of magnetic white dwarfs (WD) have revealed a higher incidence of magnetism in older WDs. Specifically, these studies indicate that magnetism is more prevalent in WDs with fully or partially crystallized cores compared to those with entirely liquid cores. This has led to the recognition of a crystallization-driven dynamo as an important mechanism for explaining magnetism in isolated WDs. However, recent simulations challenged the capability of this mechanism to match both the incidence of magnetism and the field strengths detected in WDs. In this letter, we explore an alternative hypothesis for the surface emergence of magnetic fields in isolated WDs. WDs with masses $\gtrsim 0.55 M_\odot$ are the descendants of main-sequence stars with convective cores capable of generating strong dynamo magnetic fields. This idea is supported by asteroseismic evidence of strong magnetic fields buried within the interiors of red giant branch stars. Assuming that these fields are disrupted by subsequent convective zones, we have estimated magnetic breakout times for WDs. Due to the significant uncertainties in breakout times stemming from the treatment of convective boundaries and mass loss rates, we cannot provide a precise prediction for the emergence time of the main-sequence dynamo field. However, we can predict that this emergence should occur during the WD phase for WDs with masses $\gtrsim 0.65 M_\odot$. We also find that the magnetic breakout is expected to occur earlier in more massive WDs, consistently with observations from volume-limited samples and the well-established fact that magnetic WDs tend to be more massive than non-magnetic ones. Moreover, within the uncertainties of stellar evolutionary models, we find that the emergence of main-sequence dynamo magnetic fields can account for a significant portion of the magnetic WDs., Comment: accepted for publication in A&A letters

We used the Javalambre Photometric Local Universe Survey (J-PLUS) DR2 photometry in twelve optical bands over 2176 deg2 to estimate the fraction of white dwarfs with presence of CaII H+K absorption along the cooling sequence. We compared the J-PLUS photometry against metal-free theoretical models to estimate the equivalent width in the J0395 passband of 10 nm centered at 395 nm (EW_J0395), a proxy to detect calcium absorption. A total of 4399 white dwarfs within 30000 > Teff > 5500 K and mass M > 0.45 Msun were analyzed. Their EW_J0395 distribution was modeled using two populations, corresponding to polluted and non-polluted systems, to estimate the fraction of calcium white dwarfs (f_Ca) as a function of Teff. The probability for each individual white dwarf of presenting calcium absorption, pca, was also computed. The comparison with both the measured Ca/He abundance and the metal pollution from spectroscopy shows that EW_J0395 correlates with the presence of calcium. The fraction of calcium white dwarfs increases from f_Ca = 0 at Teff = 13500 K to f_Ca = 0.15 at Teff = 5500 K. We compare our results with the fractions derived from the 40 pc spectroscopic sample and from SDSS spectra. The trend found in J-PLUS observations is also present in the 40 pc sample, however SDSS shows a deficit of metal-polluted objects at Teff < 12000 K. Finally, we found 39 white dwarfs with pca > 0.99. Twenty of them have spectra presented in previous studies, whereas we observed six additional targets. These 26 objects were all confirmed as metal-polluted systems. The J-PLUS optical data provide a robust statistical measurement for the presence of CaII H+K absorption in white dwarfs. We find a 15 +- 3 % increase in the fraction of calcium white dwarfs from Teff = 13500 K to 5500 K, which reflects their selection function in the optical from the total population of metal-polluted systems., Comment: Submitted to Astronomy and Astrophysics. Comments are welcome

We present new results on PHL 5038AB, a widely separated binary system composed of a white dwarf and a brown dwarf, refining the white and brown dwarf parameters and determining the binary separation to be $66^{+12}_{-24}$~AU. New spectra of the white dwarf show calcium absorption lines suggesting the hydrogen-rich atmosphere is weakly polluted, inferring the presence of planetesimals in the system, which we determine are in an S-type orbit around the white dwarf in orbits closer than 17-32 AU. We do not detect any infrared excess that would indicate the presence of a disc, suggesting all dust present has either been totally accreted or is optically thin. In this system, we suggest the metal pollution in the white dwarf atmosphere can be directly attributed to the presence of the brown dwarf companion disrupting the orbits of planetesimals within the system., Comment: 8 pages, 5 figures, 3 tables. Accepted for publication in MNRAS

White dwarfs with a F, G or K type companion represent the last common ancestor for a plethora of exotic systems throughout the galaxy, though to this point very few of them have been fully characterised in terms of orbital period and component masses, despite the fact several thousand have been identified. Gaia data release 3 has examined many hundreds of thousands of systems, and as such we can use this, in conjunction with our previous UV excess catalogues, to perform spectral energy distribution fitting in order to obtain a sample of 206 binaries likely to contain a white dwarf, complete with orbital periods, and either a direct measurement of the component masses for astrometric systems, or a lower limit on the component masses for spectroscopic systems. Of this sample of 206, four have previously been observed with Hubble Space Telescope spectroscopically in the ultraviolet, which has confirmed the presence of a white dwarf, and we find excellent agreement between the dynamical and spectroscopic masses of the white dwarfs in these systems. We find that white dwarf plus F, G or K binaries can have a wide range of orbital periods, from less than a day to many hundreds of days. A large number of our systems are likely post-stable mass transfer systems based on their mass/period relationships, while others are difficult to explain either via stable mass transfer or standard common envelope evolution., Comment: 17 pages, 10 figures, Accepted for publication in MNRAS

The Wide-field Spectroscopic Telescope (WST) is proposed as a new facility dedicated to the efficient delivery of spectroscopic surveys. This white paper summarises the initial concept as well as the corresponding science cases. WST will feature simultaneous operation of a large field-of-view (3 sq. degree), a high multiplex (20,000) multi-object spectrograph (MOS) and a giant 3x3 sq. arcmin integral field spectrograph (IFS). In scientific capability these requirements place WST far ahead of existing and planned facilities. Given the current investment in deep imaging surveys and noting the diagnostic power of spectroscopy, WST will fill a crucial gap in astronomical capability and work synergistically with future ground and space-based facilities. This white paper shows that WST can address outstanding scientific questions in the areas of cosmology; galaxy assembly, evolution, and enrichment, including our own Milky Way; origin of stars and planets; time domain and multi-messenger astrophysics. WST's uniquely rich dataset will deliver unforeseen discoveries in many of these areas. The WST Science Team (already including more than 500 scientists worldwide) is open to the all astronomical community. To register in the WST Science Team please visit https://www.wstelescope.com/for-scientists/participate, Comment: 194 pages, 66 figures. Comments are welcome (wstelescope@gmail.com)

Ultra-short period compact binaries are important sources of gravitational waves, which include e.g. the progenitors of type Ia supernovae or the progenitors of merger episodes that may lead to massive and magnetic single white dwarfs. J0526+5934 is an unresolved compact binary star with an orbital period of 20.5 minutes that belongs to this category. The visible component of J0526+5934 has been recently claimed to be a hot sub-dwarf star with a CO white dwarf companion. Our aim is to provide strong observational plus theoretical evidence that the primary star is rather an extremely-low mass white dwarf, although the hot subdwarf nature cannot be completely ruled out. We analyse optical spectra together with time-series photometry of the visible component of J0526+5934 to constrain its orbital and stellar parameters. We also employ evolutionary sequences for low-mass white dwarfs to derive independent values of the primary mass. From the analysis of our observational data, we find a stellar mass for the primary star in J0526+5934 of 0.26+-0.05 Msun, which perfectly matches the 0.237+-0.035 Msun independent measurement we derived from the theoretical evolutionary models. This value is considerably lower than the theoretically expected and generally observed mass range of hot subdwarf stars, but falls well within the mass limit values of extremely low-mass white dwarfs. We conclude J0526+5934 is the fifth ultra-short period detached double white dwarf currently known., Comment: Submitted to A&A

White dwarf (WD) stars are considered cosmic laboratories to study the physics of dense plasma. Furthermore, the use of WD stars as cosmic clocks to date stellar populations and main sequence companions demands an appropriate understanding of the WD physics in order to provide precise ages for these stars. We aim at studying exsolution in the interior of WD stars, a process in which a crystallized ionic binary mixture separates into two solid solutions with different fractions of the constituents. Depending on the parent solid mixture composition, this process can release or absorb heat, thus leading to a delay or a speed-up of WD cooling. Relying on accurate phase diagrams for exsolution, we have modeled this process in hydrogen-rich WDs with both carbon-oxygen and oxygen-neon core composition, with masses ranging from 0.53 to 1.29Msun and from 1.10 to 1.29Msun, respectively. Exsolution is a slow process that takes place at low luminosities (log(L/Lsun)$\lesssim$-2.75) and effective temperatures (Teff$\lesssim$18 000K) in WDs. We find that exsolution begins at brighter luminosities in CO than in ONe WDs of the same mass. Massive WDs undergo exsolution at brighter luminosities than their lower-mass counterparts. The net effect of exsolution on the WD cooling times depends on the stellar mass and the exact chemical profile. For standard core chemical profiles and preferred assumptions regarding miscibility gap microphysics, the cooling delay can be as large as ~0.35 Gyrs at L/Lsun ~ -5. We have neglected a chemical redistribution possibly associated with this process, which could lead to a further cooling delay. Exsolution has a marginal effect on the WD cooling times and, accordingly, we find no WD branches on the Gaia color magnitude diagram associated with it. However, exsolution in massive WDs can alter the faint end of the WD luminosity function, thus impacting WD cosmochronology., Comment: Accepted for publication in Astronomy & Astrophysics

Magnetic activity and rotation are related to the age of low-mass main sequence stars. To further constrain these relations, we study a sample of 574 main sequence stars members of common proper motion pairs with white dwarfs, identified thanks to Gaia astrometry. We use the white dwarfs as age indicators, while the activity indexes and rotational velocities are obtained from the main sequence companions using standard procedures. We find that stars older than 5 Gyr do not display Halpha nor Caii H&K emission unless they are fast rotators due to tidal locking from the presence of unseen companions and that the rotational velocities tend to decrease over time, thus supporting the so-called gyrochronology. However, we also find moderately old stars (~2-6 Gyr) that are active presumably because they rotate faster than they should for their given ages. This indicates that they may be suffering from weakened magnetic braking or that they possibly evolved through wind accretion processes in the past. The activity fractions that we measure for all stars younger than 5 Gyr range between ~10-40 per cent. This is line with the expectations, since our sample is composed of F, G, K and early M stars, which are thought to have short (<2 Gyr) activity lifetimes. Finally, we observe that the Halpha fractional luminosities and the R'HK indexes for our sample of (slowly rotating) stars show a spread (-4>log(LHalpha/Lbol); log(R'HK)> -5) typically found in inactive M stars or weakly active/inactive F, G, K stars., Comment: Acceted for publication by MNRAS

Planetary nebulae in Galactic open star clusters are rare objects; only three are known to date. They are of particular interest because their distance can be determined with high accuracy, allowing one to characterize the physical properties of the planetary nebula and its ionizing central star with high confidence. Here we present the first quantitative spectroscopic analysis of a central star in an open cluster, namely the faint nucleus of IPHASX J055226.2$+$323724 in M37. This cluster contains 14 confirmed white dwarf members, which were previously used to study the initial-to-final-mass relation of white dwarfs, and six additional white dwarf candidates. We performed an atmosphere modeling of spectra taken with the 10m Gran Telescopio Canarias. The central star is a hot hydrogen-deficient white dwarf with an effective temperature of 90,000 K and spectral type PG1159 (helium- and carbon-rich). We know it is about to transform into a helium-rich DO white dwarf because the relatively low atmospheric carbon abundance indicates ongoing gravitational settling of heavy elements. The star belongs to a group of hot white dwarfs that exhibit ultrahigh-excitation spectral lines possibly emerging from shock-heated material in a magnetosphere. We find a relatively high stellar mass of $M= 0.85^{+0.13}_{-0.14}$ M$_\odot$. This young white dwarf is important for the semi-empirical initial-final mass relation because any uncertainty related to white-dwarf cooling theory is insignificant with respect to the pre-white-dwarf timescale. Its post-asymptotic-giant-branch age of $170,000-480,000$ yr suggests that the extended planetary nebula is extraordinarily old. We also performed a spectroscopic analysis of the six other white dwarf candidates of M37, confirming one as a cluster member., Comment: Accepted for publication in A&A

The third data release of Gaia has provided approximately 220 million low resolution spectra. Among these, about 100,000 correspond to white dwarfs. The magnitude of this quantity of data precludes the possibility of performing spectral analysis and type determination by human inspection. In order to tackle this issue, we explore the possibility of utilising a machine learning approach, based on a Random Forest algorithm. We aim to analyze the viability of the Random Forest algorithm for the spectral classification of the white dwarf population within 100 pc from the Sun, based on the Hermite coefficients of Gaia spectra. We utilized the assigned spectral type from the Montreal White Dwarf Database for training and testing our Random Forest algorithm. Once validated, our algorithm model is applied to the rest of unclassified white dwarfs within 100 pc. First, we started by classifying the two major spectral type groups of white dwarfs: hydrogen-rich (DA) and hydrogen-deficient (non-DA). Next, we explored the possibility of classifying the various spectral subtypes, including in some cases the secondary spectral types. Our Random Forest classification presented a very high recall (>80%) for DA and DB white dwarfs, and a very high precision (>90%) for DB, DQ and DZ white dwarfs. As a result we have assigned a spectral type to 9,446 previously unclassified white dwarfs: 4,739 DAs, 76 DBs (60 of them DBAs), 4,437 DCs, 132 DZs and 62 DQs (9 of them DQpec). Despite the low resolution of Gaia spectra, the Random Forest algorithm applied to the Gaia spectral coefficients proves to be a highly valuable tool for spectral classification., Comment: 15 pages, 17 figures, 3 tables. Submitted to Astronomy & Astrophysics

The characterization of white dwarf atmospheres is crucial for accurately deriving stellar parameters such as effective temperature, mass, and age. We aim to classify the population of white dwarfs up to 500 pc into hydrogen-rich or hydrogen-deficient atmospheres based on Gaia spectra and to derive an accurate spectral type-temperature distribution of white dwarfs as a function of the effective temperature for the largest observed unbiased sample of these objects. We took advantage of the recent Gaia low-resolution spectra available for 76,657 white dwarfs up to 500 pc. We calculated synthetic J-PAS narrow-band photometry and fitted the spectral energy distribution of each object with up-to-date models for hydrogen-rich and helium-rich white dwarf atmospheres. We estimated the probability for a white dwarf to have a hydrogen-rich atmosphere and validated the results using the Montreal White Dwarf Database. Finally, precise effective temperature values were derived for each object using La Plata evolutionary models. We have successfully classified a total of 65,310 white into DAs and non-DAs with an accuracy of 94%. An unbiased subsample of nearly 34,000 objects was built, from which we computed a precise spectral distribution spanning an effective temperature range from 5,500 to 40,000 K, while accounting for potential selection effects. Some characteristic features of the spectral evolution, such as the deficit of helium-rich stars at T_eff $\approx$35,000-40,000 K and in the range 22,000 < T_eff < 25,000 K, as well as a gradual increase from 18,000K to T_eff $\approx$7,000K, where the non-DA stars percentage reaches its maximum of 41%, followed by a decrease for cooler temperatures, are statistically significant. These findings will provide precise constraints for the proposed models of spectral evolution., Comment: 7 pages, 11 pages. Accepted for publication in A&A

Planets orbiting close to hot stars experience intense extreme-ultraviolet radiation, potentially leading to atmosphere evaporation and to thermal dissociation of molecules. However, this extreme regime remains mainly unexplored due to observational challenges. Only a single known ultra-hot giant planet, KELT-9b, receives enough ultraviolet radiation for molecular dissociation, with a day-side temperature of ~4,600K. An alternative approach uses irradiated brown dwarfs as hot-Jupiter analogues. With atmospheres and radii similar to those of giant planets, brown dwarfs orbiting close to hot Earth-sized white-dwarf stars can be directly detected above the glare of the star. Here we report observations revealing an extremely irradiated low-mass companion to the hot white dwarf WD0032-317. Our analysis indicates a day-side temperature of ~8,000K, and a day-to-night temperature difference of ~6,000K. The amount of extreme-ultraviolet radiation (with wavelengths 100-912\r{A}) received by WD0032-317B is equivalent to that received by planets orbiting close to stars as hot as a late B-type stars, and about 5,600 times higher than that of KELT-9b. With a mass of ~75-88 Jupiter masses, this near-hydrogen-burning-limit object is potentially one of the most massive brown dwarfs known., Comment: Authors' version of the article published in Nature Astronomy (DOI https://doi.org/10.1038/s41550-023-02048-z)

We employ the La Plata stellar evolution code, LPCODE, to compute the first set of constant rest-mass carbon-oxygen ultra-massive white dwarf evolutionary sequences for masses higher than 1.29 Msun that fully take into account the effects of general relativity on their structural and evolutionary properties. In addition, we employ the LP-PUL pulsation code to compute adiabatic g-mode Newtonian pulsations on our fully relativistic equilibrium white dwarf models. We find that carbon-oxygen white dwarfs more massive than 1.382 Msun become gravitationally unstable with respect to general relativity effects, being this limit higher than the 1.369 Msun we found for oxygen-neon white dwarfs. As the stellar mass approaches the limiting mass value, the stellar radius becomes substantially smaller compared with the Newtonian models. Also, the thermo-mechanical and evolutionary properties of the most massive white dwarfs are strongly affected by general relativity effects. We also provide magnitudes for our cooling sequences in different passbands. Finally, we explore for the first time the pulsational properties of relativistic ultra-massive white dwarfs and find that the period spacings and oscillation kinetic energies are strongly affected in the case of most massive white dwarfs. We conclude that the general relativity effects should be taken into account for an accurate assessment of the structural, evolutionary, and pulsational properties of white dwarfs with masses above 1.30 Msun., Comment: 12 pages, 12 figures, accepted for publication in MNRAS. arXiv admin note: text overlap with arXiv:2208.14144

The ESA Gaia space mission has revealed a bifurcation of the white dwarf (WD) sequence on the color magnitude diagram in two branches: A and B. While the A branch consists mostly of WDs with H-rich atmospheres, the B branch is not completely understood. Although invoked to be populated mainly by He-rich WDs, the B branch overlaps a $\sim 0.8M_\odot$ evolutionary track with a pure He envelope, fact that would imply an unexpected peak in the WD mass distribution. In cold He-rich WDs, it is expected that the outer convective zone penetrates into deep C-rich layers, thus leading to a slight C contamination in their surfaces at $\sim 10,000$K. Here we aim at studying the Gaia bifurcation as the natural consequence of C dredge-up by convection in cold He-dominated WDs. Relying on accurate atmosphere models, we provide a new set of evolutionary models for He-rich WDs employing different prescriptions for the C enrichment. On the basis of these models, we made a population synthesis study of the Gaia 100pc WD sample to constrain the models that best fit the bifurcation. Our study shows that He-rich WD models with a slight C contamination below the optical detection limit can accurately reproduce the Gaia bifurcation. We refer to these stars as stealth DQ WDs because they do not exhibit detectable C signatures in their optical spectra, but the presence of C in their atmosphere produces a continuum absorption favouring the emission in bluer wavelengths, thereby creating the B branch of the bifurcation. Also, we show that the mass distribution for He-rich WDs obtained when a stealth C contamination is considered is consistent with the mass distribution for H-rich WDs and with the standard evolutionary channels for their formation. We conclude that stealth DQ WDs can account for the lower branch in the Gaia bifurcation. The C signatures of these stars could be detectable in Ultra-Violet spectra., Comment: 9 pages, 8 figures, accepted for publication in Astronomy & Astrophysics

Wide-field time-domain photometric sky surveys are now finding hundreds of eclipsing white dwarf plus M dwarf binaries, a population encompassing a wealth of information and potential insight into white dwarf and close binary astrophysics. Precise follow-up observations are essential in order to fully constrain these systems and capitalise on the power of this sample. We present the first results from our program of high-speed, multi-band photometric follow-up. We develop a method to measure temperatures, (model-dependent) masses, and radii for both components from the eclipse photometry alone and characterize 34 white dwarf binaries, finding general agreement with independent estimates using an alternative approach while achieving around a factor of two increase in parameter precision. In addition to these parameter estimates, we discover a number of interesting systems -- finding four with sub-stellar secondaries, doubling the number of eclipsing examples, and at least six where we find the white dwarf to be strongly magnetic, making these the first eclipsing examples of such systems and key to investigating the mechanism of magnetic field generation in white dwarfs. We also discover the first two pulsating white dwarfs in detached and eclipsing post-common-envelope binaries -- one with a low-mass, likely helium core, and one with a relatively high mass, towards the upper end of the known sample of ZZ Cetis. Our results demonstrate the power of eclipse photometry, not only as a method of characterising the population, but as a way of discovering important systems that would have otherwise been missed by spectroscopic follow-up., Comment: 12 pages with a 5 page appendix and 14 figures. Accepted for publication in MNRAS

The third data release of Gaia has provided low resolution spectra for ~100,000 white dwarfs (WDs) that, together with the excellent photometry and astrometry, represent an unrivalled benchmark for the study of this population. In this work, we first built a highly-complete volume-limited sample consisting in 12,718 WDs within 100 pc from the Sun. The use of VOSA tool allowed us to perform an automated fitting of their spectral energy distributions to different atmospheric models. In particular, the use of spectrally derived J-PAS photometry from Gaia spectra led to the classification of DA and non-DA WDs with an accuracy >90%, tested in already spectroscopically labelled objects. The excellent performance achieved was extended to practically the whole population of WDs with effective temperatures above 5500 K. Our results show that, while the A branch of the Gaia WD Hertzsprung-Russell diagram is practically populated by DA WDs, the B branch is largely formed by non-DAs (65%). The remaining 35% of DAs within the B branch implies a second peak at ~0.8 Mo in the DA-mass distribution. Additionally, the Q branch and its extension to lower temperatures can be observed for both DA and non-DA objects due to core crystallisation. Finally, we derived a detailed spectral evolution function, which confirms a slow increase of the fraction of non-DAs as the effective temperature decreases down to 10,500 K, where it reaches a maximum of 36% and then decreases for lower temperatures down to ~31%., Comment: Accepted for publication in MNRAS

Virtually all binaries consisting of a white dwarf with a non-degenerate companion can be classified as either close post-interaction systems (with orbital periods of a few days or less), or wide systems (with periods longer than decades), in which both components have effectively evolved as single stars. Binaries with periods between these two extremes can help constrain common envelope efficiency, or highlight alternative pathways towards the creation of compact binaries. To date such binaries have remained mostly elusive. Here we present three white dwarfs in binaries with evolved subgiant stars with orbital periods of 41, 52 and 461 d. Using Hubble Space Telescope spectroscopy we find that all three systems contain low mass white dwarfs ($\leq$0.4 M$_{\odot}$). One system, TYC 8394$-$1331$-$1, is the inner binary of a hierarchical triple, where the white dwarf plus subgiant binary is orbited by a more distant companion star. These binaries were likely formed from a phase of stable but non-conservative mass transfer, as opposed to common envelope evolution. All three systems will undergo a common envelope phase in the future, but the two shorter period systems are expected to merge during this event, while the longest period system is likely to survive and create a close binary with two low mass white dwarfs., Comment: 13 pages, 9 figures. Accepted for publication in MNRAS

We present a spectroscopic survey of 248 white dwarf candidates within 40 pc of the Sun; of these 244 are in the southern hemisphere. Observations were performed mostly with the Very Large Telescope (X-Shooter) and Southern Astrophysical Research Telescope. Almost all candidates were selected from $\textit{Gaia}$ Data Release 3 (DR3). We find a total of 246 confirmed white dwarfs, 209 of which had no previously published spectra, and two main-sequence star contaminants. Of these, 100 white dwarfs display hydrogen Balmer lines, 69 have featureless spectra, and two show only neutral helium lines. Additionally, 14 white dwarfs display traces of carbon, while 37 have traces of other elements that are heavier than helium. We observe 36 magnetic white dwarfs through the detection of Zeeman splitting of their hydrogen Balmer or metal spectral lines. High spectroscopic completeness (> 97 per cent) has now been reached, such that we have 1058 confirmed $\textit{Gaia}$ DR3 white dwarfs out of 1083 candidates within 40 pc of the Sun at all declinations., Comment: 49 pages, 19 figures. Accepted by MNRAS on 8 November, 2022

The white dwarf binary pathways survey is dedicated to studying the origin and evolution of binaries containing a white dwarf and an intermediate-mass secondary star of the spectral type A, F, G, or K (WD+AFGK). Here we present CPD-65\,264, a new post common envelope binary with an orbital period of 1.37\,days that contains a massive white dwarf ($ 0.86\pm 0.06\,\mathrm{M}_{\odot}$) and an intermediate-mass ($1.00\pm0.05\,\mathrm{M}_{\odot}$) main-sequence secondary star. We characterized the secondary star and measured the orbital period using high-resolution optical spectroscopy. The white dwarf parameters are determined from HST spectroscopy. In addition, TESS observations revealed that up to 19 per cent of the surface of the secondary is covered with starspots. Small period changes found in the light curve indicate that the secondary is the second example of a G-type secondary star in a post-common envelope binary with latitudinal differential rotation. Given the relatively large mass of the white dwarf and the short orbital period, future mass transfer will be dynamically and thermally stable and the system will evolve into a cataclysmic variable. The formation of the system can be understood assuming common envelope evolution without contributions from energy sources besides orbital energy. CPD-65\,264 is the seventh post-common envelope binaries with intermediate-mass secondaries that can be understood assuming a small efficiency in the common envelope energy equation, in agreement with findings for post-common envelope binaries with M-dwarf or sub-stellar companions., Comment: 10 pages, 10 figures, 3 tables. Accepted for publication in MNRAS

We present the first set of constant rest-mass ultra-massive oxygen/neon white dwarf cooling tracks with masses larger than 1.29 Msun which fully take into account the effects of general relativity on their structural and evolutionary properties. We have computed the full evolution sequences of 1.29, 1.31, 1.33, 1.35, and 1.369 Msun white dwarfs with the La Plata stellar evolution code, LPCODE. For this work, the standard equations of stellar structure and evolution have been modified to include the full effects of general relativity. For comparison purposes, the same sequences have been computed but for the Newtonian case. According to our calculations, the evolutionary properties of the most massive white dwarfs are strongly modified by general relativity effects. In particular, the resulting stellar radius is markedly smaller in the general relativistic case, being up to 25% smaller than predicted by the Newtonian treatment for the more massive ones. We find that oxygen/neon white dwarfs more massive than 1.369 Msun become gravitationally unstable with respect to general relativity effects. When core chemical distribution due to phase separation on crystallization is considered, such instability occurs at somewhat lower stellar masses, greater than 1.36 Msun. In addition, cooling times for the most massive white dwarf sequences result in about a factor of two smaller than in the Newtonian case at advanced stages of evolution. Finally, a sample of white dwarfs has been identified as ideal candidates to test these general relativistic effects. We conclude that the general relativity effects should be taken into account for an accurate assessment of the structural and evolutionary properties of the most massive white dwarfs., Comment: 12 pages, accepted for publication in Astronomy and Astrophysics

Ultra-massive white dwarfs ($ 1.05 \rm M_\odot \lesssim M_{WD}$) are particularly interesting objects that allow us to study extreme astrophysical phenomena such as type Ia supernovae explosions and merger events. Traditionally, ultra-massive white dwarfs are thought to harbour oxygen-neon (ONe) cores. However, recent theoretical studies and new observations suggest that some ultra-massive white dwarfs could harbour carbon-oxygen (CO) cores. Although several studies have attempted to elucidate the core composition of ultra-massive white dwarfs, to date, it has not been possible to distinguish them through their observed properties. Here, we present a new method for revealing the core-chemical composition in ultra-massive white dwarfs that is based on the study of magnetic fields generated by convective mixing induced by the crystallization process. ONe white dwarfs crystallize at higher luminosities than their CO counterparts. Therefore, the study of magnetic ultra-massive white dwarfs in the particular domain where ONe cores have reached the crystallization conditions but CO cores have not, may provide valuable support to their ONe core-chemical composition, since ONe white dwarfs would display signs of magnetic fields and CO would not. We apply our method to eight white dwarfs with magnetic field measurements and we suggest that these stars are candidate ONe white dwarfs., Comment: Accepted for publication in MNRAS Letters

The third Gaia data release provides low-resolution spectra for around 200 million sources. It is expected that a sizeable fraction of them contain a white dwarf (WD), either isolated, or in a binary system with a main-sequence (MS) companion, i.e. a white dwarf-main sequence (WDMS) binary. Taking advantage of a consolidated Random Forest algorithm used in the classification of WDs, we extend it to study the feasibility of classifying Gaia WDMS binary spectra. The Random Forest algorithm is first trained with a set of synthetic spectra generated by combining individual WD and MS spectra for the full range of effective temperatures and surface gravities. Moreover, with the aid of a detailed population synthesis code, we simulate the Gaia spectra for the above mentioned populations. For evaluating the performance of the models, a set of metrics are applied to our classifications. Our results show that for resolving powers above ~300 the accuracy of the classification depends exclusively on the SNR of the spectra, while below that value the SNR should be increased as the resolving power is reduced to maintain a certain accuracy. The algorithm is then applied to the already classified SDSS WDMS catalogue, revealing that the automated classification exhibits an accuracy comparable (or even higher) to previous classification methods. Finally, we simulate the {\it Gaia} spectra, showing that our algorithm is able to correctly classify nearly 80% the synthetic WDMS spectra. Our algorithm represents a useful tool in the analysis and classification of real Gaia WDMS spectra. Even for those spectra dominated by the flux of the MS stars, the algorithm reaches a high degree of accuracy (60%)., Comment: 15 pages, 11 figures, 2 tables, submitted to A&A

We identify Gaia0007-1605AC as the first inner brown dwarf-white dwarf binary of a hierarchical triple system in which the outer component is another white dwarf (Gaia0007-1605B). From optical/near-infrared spectroscopy obtained at the Very Large Telescope with the X-Shooter instrument and/or from Gaia photometry plus SED fitting, we determine the effective temperatures and masses of the two white dwarfs (12018+-68 K, 0.54+-0.01 Msun for Gaia0007-1605A and 4445+-116 K, 0.56+-0.05 Msun for Gaia0007-1605B) and the effective temperature of the brown dwarf (1850+-50 K; corresponding to a spectral type L3+-1). By analysing the available TESS light curves of Gaia0007-1605AC we detect a signal at 1.0446+-0.0015 days with an amplitude of 6.25 ppt, which we interpret as the orbital period modulated from irradiation effects of the white dwarf on the brown dwarf's surface. This drives us to speculate that the inner binary evolved through a common envelope phase in the past. Using the outer white dwarf as a cosmochronometer and analysing the kinematic properties of the system, we conclude that the triple system is about 10 Gyr old., Comment: ApJL in press

Establishing a large sample of post common envelope binaries (PCEBs) that consist of a white dwarf plus an intermediate mass companion star of spectral type AFGK, offers the potential to provide new constraints on theoretical models of white dwarf binary formation and evolution. Here we present a detailed analysis of two new systems, TYC 110-755-1 and TYC 3858-1215-1. Based on radial velocity measurements we find the orbital periods of the two systems to be $\sim$ 0.85 and $\sim$ 1.64 days, respectively. In addition, HST spectroscopy of TYC 110-755-1 allowed us to measure the mass of the white dwarf in this system (0.78 M$_\odot$). We furthermore analysed TESS high time resolution photometry and find both secondary stars to be magnetically extremely active. Differences in the photometric and spectroscopic periods of TYC 110-755-1 indicate that the secondary in this system is differentially rotating. Finally, studying the past and future evolution of both systems, we conclude that the common envelope efficiency is likely similar in close white dwarf plus AFGK binaries and PCEBs with M-dwarf companions and find a wide range of possible evolutionary histories for both systems. While TYC 3858-1215-1 will run into dynamically unstable mass transfer that will cause the two stars to merge and evolve into a single white dwarf, TYC 110-755-1 is a progenitor of a cataclysmic variable system with an evolved donor star., Comment: Accepted for publication in MNRAS, 14 pages, 11 figures

The $Gaia$ mission has provided an unprecedented wealth of information about the white dwarf population of our Galaxy. In particular, our studies show that the sample up to 100\,pc from the Sun can be considered as practically complete. This fact allows us to estimate a precise fraction of double-degenerate ($1.18\pm 0.10$\%) and white dwarf plus main-sequence stars ($6.31\pm0.23$\%) among all white dwarfs through comoving pairs identification. With the aid of a detailed population synthesis code we are able to reproduce synthetic white dwarf populations with nearly identical fractions as the observed ones, thus obtaining valuable information about the binary fraction, $f_{\rm b}$, initial mass ratio distribution, $n(q)$, and initial separation distribution, $f(a)$, among other parameters. Our best-fit model is achieved within a $1\sigma$ confidence level for $f(a)\propto a^{-1}$, $n(q)\propto q^{n_q}$, with $n_q=-1.13^{+0.12}_{-0.10}$ and $f_{\rm b}=0.32\pm 0.02$. The fraction of white dwarf mergers generated by this model is $9\sim16\%$, depending on the common-envelope treatment. As sub-products of our modelling we find that around $1\sim3\%$ of the white dwarf population are unresolved double-degenerates and that only $\sim1\%$ of all white dwarfs contain a He-core. Finally, only a mild kick during white dwarf formation seems to be necessary for fitting the observed sky separation of double-degenerate systems., Comment: 14 pages, 8 figures, 5 tables, accepted for publication in MNRAS

We present the first detailed 3D kinematic analysis of a sample of 3133 white dwarfs that use Gaia astrometry plus radial velocities, which were measured either by Gaia or by ground based spectroscopic observations. The sample includes either isolated white dwarfs that have direct radial velocity measurements, or white dwarfs that belong to common proper motion pairs that contain nondegenerate companions with available radial velocities. A subset of common proper motion pairs also have metal abundances that have been measured by large scale spectroscopic surveys or by our own follow-up observations. We used the white dwarfs as astrophysical clocks, by determining their masses and total ages via interpolation with dedicated evolutionary models. We also used the nondegenerate companions in common proper motions to chemically tag of the population. Combining accurate radial velocities with Gaia astrometry and proper motions, we derived the velocity components of our sample in the Galactic rest frame and their Galactic orbital parameters. The sample is mostly located within ~300 pc from the Sun. It contains predominantly (90-95%) thin-disk stars with close-to-circular Galactic orbits, while the remaining 5-10% of stars have more eccentric trajectories and belong to the thick disk. We identified seven isolated white dwarfs and two common proper motion pairs as halo members. We determined the age - velocity dispersion relation for the thin-disk members, which shows signatures of dynamical heating and saturation after 4-6 Gyr. We observed a mild anti-correlation between [Fe/H] and the radial component of the average - velocity dispersion, showing that dynamical mixing of populations takes place in the Galactic disk, as was detected via the analysis of other samples of FGK stars. [Abridged], Comment: Accepted for publication (Raddi et al. 2022, A&A, 658, A22). 16 pages, 14 figures, 9 tables. Tables 2, 4, 5, 6, and 7 are available as pdf files and in electronic form via http://cdsarc.u-strasbg.fr/viz-bin/cat/J/A+A/658/A22

We use the data provided by the Gaia Early Data Release 3 to search for a highly-complete volume-limited sample of unresolved binaries consisting of a white dwarf and a main sequence companion (i.e. WDMS binaries) within 100pc. We select 112 objects based on their location within the Hertzsprung-Russell diagram, of which 97 are new identifications. We fit their spectral energy distributions (SED) with a two-body fitting algorithm implemented in VOSA (Virtual Observatory SED Analyser) to derive the effective temperatures, luminosities and radii (hence surface gravities and masses) of both components. The stellar parameters are compared to those from the currently largest catalogue of close WDMS binaries, from the Sloan Digital Sky Survey (SDSS). We find important differences between the properties of the Gaia and SDSS samples. In particular, the Gaia sample contains WDMS binaries with considerably cooler white dwarfs and main sequence companions (some expected to be brown dwarfs). The Gaia sample also shows an important population of systems consisting of cool and extremely low-mass white dwarfs, not present in the SDSS sample. Finally, using a Monte Carlo population synthesis code, we find that the volume-limited sample of systems identified here seems to be highly complete (~80+-9 per cent), however it only represents ~9 per cent of the total underlying population. The missing ~91 per cent includes systems in which the main sequence companions entirely dominate the SEDs. We also estimate an upper limit to the total space density of close WDMS binaries of ~(3.7+-1.9)x10^{-4} pc{-3}., Comment: Accepted for publication by MNRAS

Infrared excesses around white dwarf stars indicate the presence of various astrophysical objects of interest, including companions and debris disks. In this second paper of a series, we present follow-up observations of infrared excess candidates from Gaia and unWISE discussed in the first paper, Paper I. We report space-based infrared photometry at 3.6 and 4.5 micron for 174 white dwarfs from the Spitzer Space Telescope and ground-based near-infrared J, H, and K photometry of 235 white dwarfs from Gemini Observatory with significant overlap between Spitzer and Gemini observations. This data is used to confirm or rule-out the observed unWISE infrared excess. From the unWISE-selected candidate sample, the most promising infrared excess sample comes from both colour and flux excess, which has a Spitzer confirmation rate of 95%. We also discuss a method to distinguish infrared excess caused by stellar or sub-stellar companions from potential dust disks. In total, we confirm the infrared excess around 62 white dwarfs, 10 of which are likely to be stellar companions. The remaining 52 bright white dwarf with infrared excess beyond two microns has the potential to double the known sample of white dwarfs with dusty exoplanetary debris disks. Follow-up high-resolution spectroscopic studies of a fraction of confirmed excess white dwarfs in this sample have discovered emission from gaseous dust disks. Additional investigations will be able to expand the parameter space from which dust disks around white dwarfs are found., Comment: 23 pages, 8 figures, 8 tables, 3 of the tables are available digitally, and 329 SEDs are available as digital content

The age-metallicity relation is a fundamental tool for constraining the chemical evolution of the Galactic disc. In this work we analyse the observational properties of this relation using binary stars that have not interacted consisting of a white dwarf - from which we can derive the total age of the system - and a main sequence star - from which we can derive the metallicity as traced by the [Fe/H] abundances. Our sample consists of 46 widely separated, but unresolved spectroscopic binaries identified within the Sloan Digital Sky Survey, and 189 white dwarf plus main sequence common proper motion pairs identified within the second data release of Gaia. This is currently the largest white dwarf sample for which the metallicity of their progenitors have been determined. We find a flat age-metallicity relation displaying a scatter of [Fe/H] abundances of approximately 0.5 dex around the solar metallicity at all ages. This independently confirms the lack of correlation between age and metallicity in the solar neighbourhood that is found in previous studies focused on analysing single main sequence stars and open clusters., Comment: Accepted for publication by MNRAS

We analyze the volume-limited nearly complete 100 pc sample of 95 halo white dwarf candidates identified by the second data release of Gaia. Based on a detailed population synthesis model, we apply a method that relies on Gaia astrometry and photometry to accurately derive the individual white dwarf parameters (mass, radius, effective temperature, bolometric luminosity and age). This method is tested with 25 white dwarfs of our sample for which we took optical spectra and performed spectroscopic analysis. We build and analyse the halo white dwarf luminosity function, for which we find for the first time possible evidences of the cut-off at its faintest end, leading to an age estimate of $\simeq12\pm0.5 $Gyr. The mass distribution of the sample peaks at $0.589\,M_{\odot}$, with $71\%$ of the white dwarf masses below $0.6\,M_{\odot}$ and just two massive white dwarfs of more than $0.8\,M_{\odot}$. From the age distribution we find three white dwarfs with total ages above 12 Gyr, of which J1312-4728 is the oldest white dwarf known with an age of $12.41\pm0.22 $Gyr. We prove that the star formation history is mainly characterised by a burst of star formation that occurred from 10 to 12 Gyr in the past, but extended up to 8 Gyr. We also find that the peak of the star formation history is centered at around 11 Gyr, which is compatible with the current age of the Gaia-Enceladus encounter. Finally, $13\%$ of our halo sample is contaminated by high-speed young objects (total age<7 Gyr). The origin of these white dwarfs is unclear but their age distribution may be compatible with the encounter with the Sagittarius galaxy., Comment: 15 pages, 9 figures, 2 tables; accepted for publication in MNRAS

Constraints from surveys of post common envelope binaries (PCEBs) consisting of a white dwarf plus an M-dwarf companion have led to significant progress in our understanding of the formation of close white dwarf binary stars with low-mass companions. The white dwarf binary pathways project aims at extending these previous surveys to larger secondary masses, i.e. secondary stars of spectral type AFGK. Here we present the discovery and observational characterization of three PCEBs with G-type secondary stars and orbital periods between 1.2 and 2.5 days. Using our own tools as well as MESA we estimate the evolutionary history of the binary stars and predict their future. We find a large range of possible evolutionary histories for all three systems and identify no indications for differences in common envelope evolution compared to PCEBs with lower mass secondary stars. Despite their similarities in orbital period and secondary spectral type, we estimate that the future of the three systems are very different: TYC 4962-1205-1 is a progenitor of a cataclysmic variable system with an evolved donor star, TYC 4700-815-1 will run into dynamically unstable mass transfer that will cause the two stars to merge, and TYC 1380-957-1 may appear as super soft source before becoming a rather typical cataclysmic variable star., Comment: Accepted for publication in MNRAS, 13 pages, 9 figures, 5 tables

(Abridged abstract) We explore the formation of ultra-massive (M_{\rm WD} \gtrsim 1.05 M_\sun$), carbon-oxygen core white dwarfs resulting from single stellar evolution. We also study their evolutionary and pulsational properties and compare them with those of the ultra-massive white dwarfs with oxygen-neon cores resulting from carbon burning in single progenitor stars, and with binary merger predictions. We consider two single-star evolution scenarios for the formation of ultra-massive carbon-oxygen core white dwarfs that involve rotation of the degenerate core after core helium burning and reduced mass-loss rates in massive asymptotic giant-branch stars. We compare our findings with the predictions from ultra-massive white dwarfs resulting from the merger of two equal-mass carbon-oxygen core white dwarfs, by assuming complete mixing between them and a carbon-oxygen core for the merged remnant. The resulting ultra-massive carbon-oxygen core white dwarfs evolve markedly slower than their oxygen-neon counterparts. Our study strongly suggests the formation of ultra-massive white dwarfs with carbon-oxygen core from single stellar evolution. We find that both the evolutionary and pulsation properties of these white dwarfs are markedly different from those of their oxygen-neon core counterparts and from those white dwarfs with carbon-oxygen core that might result from double degenerate mergers. This can eventually be used to discern the core composition of ultra-massive white dwarfs and their formation scenario., Comment: 15 pages, 10 figures, 2 tables. To be published in Astronomy and Astrophysics

Close white dwarf binaries consisting of a white dwarf and an A, F, G or K type main sequence star, henceforth close WD+AFGK binaries, are ideal systems to understand the nature of type Ia supernovae progenitors and to test binary evolution models. In this work we identify 775 WD+AFGK candidates from TGAS (The Tycho-Gaia Astrometric Solution) and Gaia Data Release 2 (DR2), a well-defined sample of stars with available parallaxes, and we measure radial velocities (RVs) for 275 of them with the aim of identifying close binaries. The RVs have been measured from high resolution spectra obtained at the Xinglong 2.16m Telescope and the San Pedro M\'artir 2.12m Telescope and/or from available LAMOST DR6 (low-resolution) and RAVE DR5 (medium-resolution) spectra. We identify 23 WD+AFGK systems displaying more than 3$\sigma$ RV variation among 151 systems for which the measured values are obtained from different nights. Our WD+AFGK binary sample contains both AFGK dwarfs and giants, with a giant fraction $\sim$43%. The close binary fractions we determine for the WD+AFGK dwarf and giant samples are $\simeq$24% and $\simeq$15%, respectively. We also determine the stellar parameters (i.e. effective temperature, surface gravity, metallicity, mass and radius) of the AFGK companions with available high resolution spectra. The stellar parameter distributions of the AFGK companions that are members of close and wide binary candidates do not show statistically significant differences., Comment: 22 pages, 17 figures, 9 tables. Accepted for publication in The Astrophysical Journal

UCAC2 46706450 is a late-type star with an ultraviolet (UV) excess. It was considered a candidate to establish a sample of FGK stars with white dwarf (WD) companions that can be used to test binary evolution models. To verify the WD nature of the companion, UV spectroscopy was performed by Parsons et al. (2016). By a detailed model-atmosphere analysis we show that the UV source is an extremely hot WD with effective temperature $T_\mathrm{eff}$ = $105\,000\pm5000$ K, mass $M/M_\odot = 0.54\pm0.02$, radius $R/R_\odot = 0.040^{+0.005}_{-0.004}$, and luminosity $L/L_\odot= 176^{+55}_{-49}$, i.e., the compact object is just about to enter the WD cooling sequence. Investigating spectra of the cool star ($T_\mathrm{eff}$ = $4945\pm250$ K) we found that it is a K-type subgiant with $M/M_\odot = 0.8-2.4$, $R/R_\odot = 5.9^{+0.7}_{-0.5}$, and $L/L_\odot= 19^{+5}_{-5}$, that is rapidly rotating with $v \sin(i)=81$ km s$^{-1}$. Optical light curves reveal a period of two days and an o-band peak-to-peak amplitude of 0.06 mag. We suggest, that it is caused by stellar rotation in connection with star spots. With the radius we infer an extremely high rotational velocity of $v_{\mathrm{rot}}=151^{+18}_{-13}$ km s$^{-1}$, thus marking the star as one of the most rapidly rotating subgiants known. This explains chromospheric activity observed by H$\alpha$ emission and emission-line cores in CaII H and K as well as NUV flux excess. From equal and constant radial velocities of the WD and the K subgiant as well as from a fit to the spectral energy distribution we infer that they form a physical, wide though unresolved binary system. Both components exhibit similar metal abundances and show iron-group elements with slightly oversolar (up to 0.6 dex) abundance, meaning that atomic diffusion in the WD atmosphere is not yet active due to a residual, weak radiation-driven wind. (abridged), Comment: 13 pages, accepted for publication in A&A

White dwarf stars are the most common end point of stellar evolution. Of special interest are the ultramassive white dwarfs, as they are related to type Ia Supernovae explosions, merger events, and Fast Radio Bursts. Ultramassive white dwarfs are expected to harbour oxygen-neon (ONe) cores as a result of single standard stellar evolution. However, a fraction of them could have carbon-oxygen (CO) cores. Recent studies, based on the new observations provided by the {\it Gaia} space mission, indicate that a small fraction of the ultramassive white dwarfs experience a strong delay in their cooling, which cannot be attributed only to the occurrence of crystallization, thus requiring an unknown energy source able to prolong their life for long periods of time. In this study we find that the energy released by $^{22}$Ne sedimentation in the deep interior of ultramassive white dwarfs with CO cores and high $^{22}$Ne content is consistent with the long cooling delay of these stellar remnants. On the basis of a synthesis study of the white dwarf population, based on Monte Carlo techniques, we find that the observations revealed by {\it Gaia} can be explained by the existence of these prolonged youth ultramassive white dwarfs. Although such a high $^{22}$Ne abundance is not consistent with the standard evolutionary channels, our results provide sustain to the existence of CO-core ultramassive white dwarfs and to the occurrence of $^{22}$Ne sedimentation., Comment: 5 pages, 4 figures. Accepted for publication in A&A Letters

We present a spectroscopic survey of 230 white dwarf candidates within 40 pc of the Sun from the William Herschel Telescope and Gran Telescopio Canarias. All candidates were selected from $\textit{Gaia}$ Data Release 2 (DR2) and in almost all cases had no prior spectroscopic classifications. We find a total of 191 confirmed white dwarfs and 39 main-sequence star contaminants. The majority of stellar remnants in the sample are relatively cool ($\langle T_{\rm eff} \rangle$ = 6200 K), showing either hydrogen Balmer lines or a featureless spectrum, corresponding to 89 DA and 76 DC white dwarfs, respectively. We also recover two DBA white dwarfs and 9--10 magnetic remnants. We find two carbon-bearing DQ stars and 14 new metal-rich white dwarfs. This includes the possible detection of the first ultra-cool white dwarf with metal lines. We describe three DZ stars for which we find at least four different metal species, including one which is strongly Fe- and Ni-rich, indicative of the accretion of a planetesimal with core-Earth composition. We find one extremely massive (1.31 $\pm$ 0.01 M$_{\odot}$) DA white dwarf showing weak Balmer lines, possibly indicating stellar magnetism. Another white dwarf shows strong Balmer line emission but no infrared excess, suggesting a low-mass sub-stellar companion. High spectroscopic completeness ($>$99%) has now been reached for $\textit{Gaia}$ DR2 sources within 40 pc sample, in the northern hemisphere ($\delta >$ 0 deg) and located on the white dwarf cooling track in the Hertzsprung-Russell diagram. A statistical study of the full northern sample is presented in a companion paper., Comment: 37 pages, 25 figures. Accepted for publication in MNRAS

We present an overview of the sample of northern hemisphere white dwarfs within 40 pc of the Sun detected from $Gaia$ Data Release 2 (DR2). We find that 521 sources are spectroscopically confirmed degenerate stars, 111 of which were first identified as white dwarf candidates from $Gaia$ DR2 and followed-up recently with the William Herschel Telescope and Gran Telescopio Canarias. Three additional white dwarf candidates remain spectroscopically unobserved and six unresolved binaries are known to include a white dwarf but were not in our initial selection of white dwarfs in the $Gaia$ DR2 Hertzsprung-Russell diagram (HRD). Atmospheric parameters are calculated from $Gaia$ and Pan-STARRS photometry for all objects in the sample, confirming most of the trends previously observed in the much smaller 20 pc sample. Local white dwarfs are overwhelmingly consistent with Galactic disc kinematics, with only four halo candidates. We find that DAZ white dwarfs are significantly less massive than the overall DA population ($\overline{M}_\mathrm{DAZ} = 0.59\,\mathrm{M}_\odot$, $\overline{M}_\mathrm{DA} = 0.66\,\mathrm{M}_\odot$). It may suggest that planet formation is less efficient at higher mass stars, producing more massive white dwarfs. We detect a sequence of crystallised white dwarfs in the mass range from $0.6\ \lesssim M/\mathrm{M}_\odot \lesssim\ 1.0\,$ and find that the vast majority of objects on the sequence have standard kinematic properties that correspond to the average of the sample, suggesting that their nature can be explained by crystallisation alone. We also detect 56 wide binaries including a white dwarf and 26 double degenerates., Comment: 20 pages, 14 figures. Accepted for publication in MNRAS. Online tables available at http://deneb.astro.warwick.ac.uk/phrgwr/40pcTables/index.html

We perform an analysis of the single white dwarf and the double degenerate binary populations in the solar neighbourhood following a population synthesis approach to investigate the effects of unresolved double degenerates in the white dwarf luminosity function. We consider all unresolved synthetic binaries to be associated with fictitious effective temperatures and surface gravities that are obtained in the same way as if these objects were observed as single point sources. We evaluate the effects of unresolved double white dwarfs assuming that the synthetic samples are observed both by the magnitude-limited SDSS and the volume-limited Gaia surveys, the latter limited to a distance of no more than 100pc. We find that, for our standard model, the impact of unresolved double degenerates in the white dwarf luminosity function derived from the Gaia sample is nearly negligible. Unresolved double degenerates are hence expected to have no effect on the age of the Galactic disc, nor on the star formation history from this population. However, for the SDSS sample, the effect of unresolved double degenerates is significant at the brighter bins (Mbol<6.5 mag), with the fraction of such systems reaching ~40% of the total white dwarf population at Mbol=6 mag. This indicates unresolved double degenerates may influence the constraints on the star formation history derived from the SDSS white dwarf sample., Comment: accepted for publication by MNRAS

We analyse the 100pc Gaia white dwarf volume-limited sample by means of VOSA (Virtual Observatory SED Analyser) with the aim of identifying candidates for displaying infrared excesses. Our search focuses on the study of the spectral energy distribution (SED) of 3,733 white dwarfs with reliable infrared photometry and GBP-GRP colours below 0.8 mag, a sample which seems to be nearly representative of the overall white dwarf population. Our search results in 77 selected candidates, 52 of which are new identifications. For each target we apply a two-component SED fitting implemented in VOSA to derive the effective temperatures of both the white dwarf and the object causing the excess. We calculate a fraction of infrared-excess white dwarfs due to the presence of a circumstellar disk of 1.6+-0.2%, a value which increases to 2.6+-0.3% if we take into account incompleteness issues. Our results are in agreement with the drop in the percentage of infrared excess detections for cool (<8,000K) and hot (>20,000K) white dwarfs obtained in previous analyses. The fraction of white dwarfs with brown dwarf companions we derive is ~0.1-0.2%., Comment: Accepted for publication by MNRAS

We analyzed the velocity space of the thin and thick-disk Gaia white dwarf population within 100 pc looking for signatures of the Hercules stellar stream. We aimed to identify those objects belonging to the Hercules stream and, by taking advantage of white dwarf stars as reliable cosmochronometers, to derive a first age distribution. We applied a kernel density estimation to the $UV$ velocity space of white dwarfs. For the region where a clear overdensity of stars was found, we created a 5-D space of dynamic variables. We applied a hierarchichal clustering method, HDBSCAN, to this 5-D space, identifying those white dwarfs that share similar kinematic characteristics. Finally, under general assumptions and from their photometric properties, we derived an age estimate for each object. The Hercules stream was firstly revealed as an overdensity in the $UV$ velocity space of the thick-disk white dwarf population. Three substreams were then found: Hercules $a$ and Hercules $b$, formed by thick-disk stars with an age distribution peaked $4\,$Gyr in the past and extended to very old ages; and Hercules $c$, with a ratio of 65:35 thin:thick stars and a more uniform age distribution younger than 10 Gyr., Comment: 5 pages, 3 figures, accepted for publication in A&A Letters

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., Comment: 11 pages, 9 figures and 2 tables; accepted for publication in A&A

Gaia-DR2 has provided an unprecedented number of white dwarf candidates of our Galaxy. In particular, it is estimated that Gaia-DR2 has observed nearly 400,000 of these objects and close to 18,000 up to 100 pc from the Sun. This large quantity of data requires a thorough analysis in order to uncover their main Galactic population properties, in particular the thin and thick disk and halo components. Taking advantage of recent developments in artificial intelligence techniques, we make use of a detailed Random Forest algorithm to analyse an 8-dimensional space (equatorial coordinates, parallax, proper motion components and photometric magnitudes) of accurate data provided by Gaia-DR2 within 100 pc from the Sun. With the aid of a thorough and robust population synthesis code we simulated the different components of the Galactic white dwarf population to optimize the information extracted from the algorithm for disentangling the different population components. The algorithm is first tested in a known simulated sample achieving an accuracy of 85.3%. Our methodology is thoroughly compared to standard methods based on kinematic criteria demonstrating that our algorithm substantially improves previous approaches. Once trained, the algorithm is then applied to the Gaia-DR2 100 pc white dwarf sample, identifying 12,227 thin disk, 1,410 thick disk and 95 halo white dwarf candidates, which represent a proportion of 74:25:1, respectively. Hence, the numerical spatial densities are $(3.6\pm0.4)\times10^{-3}\,{\rm pc^{-3}}$, $(1.2\pm0.4)\times10^{-3}\,{\rm pc^{-3}}$ and $(4.8\pm0.4)\times10^{-5}\,{\rm pc^{-3}}$ for the thin disk, thick disk and halo components, respectively. The populations thus obtained represent the most complete and volume-limited samples to date of the different components of the Galactic white dwarf population., Comment: 18 pages, 11 figures and 3 tables. Accepted for publication in MNRAS

Interacting binaries containing white dwarfs can lead to a variety of outcomes that range from powerful thermonuclear explosions, which are important in the chemical evolution of galaxies and as cosmological distance estimators, to strong sources of low frequency gravitational wave radiation, which makes them ideal calibrators for the gravitational low-frequency wave detector LISA mission. However, current theoretical evolution models still fail to explain the observed properties of the known populations of white dwarfs in both interacting and detached binaries. Major limitations are that the existing population models have generally been developed to explain the properties of sub-samples of these systems, occupying small volumes of the vast parameter space, and that the observed samples are severely biased. The overarching goal for the next decade is to assemble a large and homogeneous sample of white dwarf binaries that spans the entire range of evolutionary states, to obtain precise measurements of their physical properties, and to further develop the theory to satisfactorily reproduce the properties of the entire population. While ongoing and future all-sky high- and low-resolution optical spectroscopic surveys allow us to enlarge the sample of these systems, high-resolution ultraviolet spectroscopy is absolutely essential for the characterization of the white dwarfs in these binaries. The Hubble Space Telescope is currently the only facility that provides ultraviolet spectroscopy, and with its foreseeable demise, planning the next ultraviolet mission is of utmost urgency., Comment: 5 pages, 2 figures Science White Paper submitted to Astro2020 Decadal Survey

Cool subdwarfs are metal-poor low-mass stars that formed during the early stages of the evolution of our Galaxy. Because they are relatively rare in the vicinity of the Sun, we know of few cool subdwarfs in the solar neighbourhood, and none with both the mass and the radius accurately determined. This hampers our understanding of stars at the low-mass end of the main-sequence. Here we report the discovery of SDSSJ235524.29+044855.7 as an eclipsing binary containing a cool subdwarf star, with a white dwarf companion. From the light-curve and the radial-velocity curve of the binary we determine the mass and the radius of the cool subdwarf and we derive its effective temperature and luminosity by analysing its spectral energy distribution. Our results validate the theoretical mass-radius-effective temperature-luminosity relations for low-mass low-metallicity stars., Comment: Accepted

Gaia16aye was a binary microlensing event discovered in the direction towards the northern Galactic disc and was one of the first microlensing events detected and alerted to by the Gaia space mission. Its light curve exhibited five distinct brightening episodes, reaching up to I=12 mag, and it was covered in great detail with almost 25,000 data points gathered by a network of telescopes. We present the photometric and spectroscopic follow-up covering 500 days of the event evolution. We employed a full Keplerian binary orbit microlensing model combined with the motion of Earth and Gaia around the Sun to reproduce the complex light curve. The photometric data allowed us to solve the microlensing event entirely and to derive the complete and unique set of orbital parameters of the binary lensing system. We also report on the detection of the first-ever microlensing space-parallax between the Earth and Gaia located at L2. The properties of the binary system were derived from microlensing parameters, and we found that the system is composed of two main-sequence stars with masses 0.57$\pm$0.05 $M_\odot$ and 0.36$\pm$0.03 $M_\odot$ at 780 pc, with an orbital period of 2.88 years and an eccentricity of 0.30. We also predict the astrometric microlensing signal for this binary lens as it will be seen by Gaia as well as the radial velocity curve for the binary system. Events such as Gaia16aye indicate the potential for the microlensing method of probing the mass function of dark objects, including black holes, in directions other than that of the Galactic bulge. This case also emphasises the importance of long-term time-domain coordinated observations that can be made with a network of heterogeneous telescopes., Comment: accepted for publication in A&A, 24 pages, 10 figures, tables with the data will be available electronically

We report the discovery of J1953-1019, the first resolved triple white dwarf system. The triplet consists of an inner white dwarf binary and a wider companion. Using Gaia DR2 photometry and astrometry combined with our follow-up spectroscopy, we derive effective temperatures, surface gravities, masses and cooling ages of the three components. All three white dwarfs have pure-hydrogen (DA) atmospheres, masses of 0.60-0.63 Msun and cooling ages of 40-290 Myr. We adopt eight initial-to-final mass relations to estimate the main sequence progenitor masses (which we find to be similar for the three components, 1.6-2.6 Msun) and lifetimes. The differences between the derived cooling times and main sequence lifetimes agree for most of the adopted initial-to-final mass relations, hence the three white dwarfs in J1953-1019 are consistent with coeval evolution. Furthermore, we calculate the projected orbital separations of the inner white dwarf binary (303.25 +- 0.01 au) and of the centre of mass of the inner binary and the outer companion (6398.97 +- 0.09 au). From these values, and taking into account a wide range of possible configurations for the triplet to be currently dynamically stable, we analyse the future evolution of the system. We find that a collision between the two inner white dwarfs due to Lidov-Kozai oscillations is unlikely, though if it occurs it could result in a sub-Chandrasekhar Type Ia supernova explosion., Comment: accepted for publication by MNRAS

Double white dwarf binaries with merger timescales smaller than the Hubble time and with a total mass near the Chandrasekhar limit (i.e. classical Chandrasekhar population) or with high-mass primaries (i.e. sub-Chandrasekhar population) are potential supernova type Ia (SNIa) progenitors. However, we have not yet unambiguously confirmed the existence of these objects observationally, a fact that has been often used to criticise the relevance of double white dwarfs for producing SNIa. We analyse whether this lack of detections is due to observational effects. To that end we simulate the double white dwarf binary population in the Galaxy and obtain synthetic spectra for the SNIa progenitors. We demonstrate that their identification, based on the detection of Halpha double-lined profiles arising from the two white dwarfs in the synthetic spectra, is extremely challenging due to their intrinsic faintness. This translates into an observational probability of finding double white dwarf SNIa progenitors in the Galaxy of (2.1+-1.0)x10^{-5} and (0.8+-0.4)x10^{-5} for the classical Chandrasekhar and the sub-Chandrasekhar progenitor populations, respectively. Eclipsing double white dwarf SNIa progenitors are found to suffer from the same observational effect. The next generation of large-aperture telescopes are expected to help in increasing the probability for detection by ~1 order of magnitude. However, it is only with forthcoming observations such as those provided by LISA that we expect to unambiguously confirm or disprove the existence of double white dwarf SNIa progenitors and to test their importance for producing SNIa., Comment: Re-submitted to MNRAS after applying all referee's suggestions

M dwarfs are prime targets in the hunt for habitable worlds around other stars. This is due to their abundance as well as their small radii and low masses and temperatures, which facilitate the detection of temperate, rocky planets in orbit around them. However, the fundamental properties of M dwarfs are difficult to constrain, often limiting our ability to characterise the planets they host. Here we test several theoretical relationships for M dwarfs by measuring 23 high precision, model-independent masses and radii for M dwarfs in binaries with white dwarfs. We find a large scatter in the radii of these low-mass stars, with 25 per cent having radii consistent with theoretical models while the rest are up to 12 per cent over-inflated. This scatter is seen in both partially- and fully-convective M dwarfs. No clear trend is seen between the over-inflation and age or metallicity, but there are indications that the radii of slowly rotating M dwarfs are more consistent with predictions, albeit with a similar amount of scatter in the measurements compared to more rapidly rotating M dwarfs. The sample of M dwarfs in close binaries with white dwarfs appears indistinguishable from other M dwarf samples, implying that common envelope evolution has a negligible impact on their structure. We conclude that theoretical and empirical mass-radius relationships lack the precision and accuracy required to measure the fundamental parameters of M dwarfs well enough to determine the internal structure and bulk composition of the planets they host., Comment: 12 pages, 10 figures, accepted for publication in MNRAS

We present a catalogue of 73,221 white dwarf candidates extracted from the astrometric and photometric data of the recently published Gaia DR2 catalogue. White dwarfs were selected from the Gaia Hertzsprung-Russell diagram with the aid of the most updated population synthesis simulator. Our analysis shows that Gaia has virtually identified all white dwarfs within 100 pc from the Sun. Hence, our sub-population of 8,555 white dwarfs within this distance limit and the colour range considered, $-\,0.52<(G_{\rm BP}-G_{\rm RP})<0.80$, is the largest and most complete volume-limited sample of such objects to date. From this sub-sample we identified 8,343 CO-core and 212 ONe-core white dwarf candidates and derived a white dwarf space density of $4.9\pm0.4\times10^{-3}\,{\rm pc^{-3}}$. A bifurcation in the Hertzsprung-Russell diagram for these sources, which our models do not predict, is clearly visible. We used the Virtual Observatory tool VOSA to derive effective temperatures and luminosities for our sources by fitting their spectral energy distributions, that we built from the UV to the NIR using publicly available photometry through the Virtual Observatory. From these parameters, we derived the white dwarf radii. Interpolating the radii and effective temperatures in hydrogen-rich white dwarf cooling sequences, we derived the surface gravities and masses. The Gaia 100 pc white dwarf population is clearly dominated by cool ($\sim$ 8,000 K) objects and reveals a significant population of massive ($M \sim 0.8 M_{\odot}$) white dwarfs, of which no more than $\sim$ $30-40 \%$ can be attributed to hydrogen-deficient atmospheres, and whose origin remains uncertain., Comment: Accepted in MNRAS (2018 July 16)

We present the data release (DR) 5 catalogue of white dwarf-main sequence (WDMS) binaries from the Large Area Multi-Object fiber Spectroscopic Telescope (LAMOST). The catalogue contains 876 WDMS binaries, of which 757 are additions to our previous LAMOST DR1 sample and 357 are systems that have not been published before. We also describe a LAMOST-dedicated survey that aims at obtaining spectra of photometrically-selected WDMS binaries from the Sloan Digital Sky Survey (SDSS) that are expected to contain cool white dwarfs and/or early type M dwarf companions. This is a population under-represented in previous SDSS WDMS binary catalogues. We determine the stellar parameters (white dwarf effective temperatures, surface gravities and masses, and M dwarf spectral types) of the LAMOST DR5 WDMS binaries and make use of the parameter distributions to analyse the properties of the sample. We find that, despite our efforts, systems containing cool white dwarfs remain under-represented. Moreover, we make use of LAMOST DR5 and SDSS DR14 (when available) spectra to measure the Na I {\lambda}{\lambda} 8183.27, 8194.81 absorption doublet and/or H{\alpha} emission radial velocities of our systems. This allows identifying 128 binaries displaying significant radial velocity variations, 76 of which are new. Finally, we cross-match our catalogue with the Catalina Surveys and identify 57 systems displaying light curve variations. These include 16 eclipsing systems, two of which are new, and nine binaries that are new eclipsing candidates. We calculate periodograms from the photometric data and measure (estimate) the orbital periods of 30 (15) WDMS binaries., Comment: 14 pages, 10 figures, accepted for publication on MNRAS