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.

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.

Ultracompact accreting binary systems each consist of a stellar remnant accreting helium-enriched material from a compact donor star. Such binaries include two related sub-classes, AM CVn-type binaries and helium cataclysmic variables, in both of which the central star is a white dwarf. We present a spectroscopic and photometric study of six accreting binaries with orbital periods in the range of 40--70 min, including phase-resolved VLT spectroscopy and high-speed ULTRACAM photometry. Four of these are AM CVn systems and two are helium cataclysmic variables. For four of these binaries we are able to identify orbital periods (of which three are spectroscopic). SDSS J1505+0659 has an orbital period of 67.8 min, significantly longer than previously believed, and longer than any other known AM CVn binary. We identify a WISE infrared excess in SDSS J1505+0659 that we believe to be the first direct detection of an AM CVn donor star in a non-direct impacting binary. The mass ratio of SDSS J1505+0659 is consistent with a white dwarf donor. CRTS J1028-0819 has an orbital period of 52.1 min, the shortest period of any helium cataclysmic variable. MOA 2010-BLG-087 is co-aligned with a K-class star that dominates its spectrum. ASASSN-14ei and ASASSN-14mv both show a remarkable number of echo outbursts following superoutbursts (13 and 10 echo outbursts respectively). ASASSN-14ei shows an increased outburst rate over the years following its superoutburst, perhaps resulting from an increased accretion rate.

We present a study of the hypothesis that white dwarfs undergo a spectral change from hydrogen- to helium-dominated atmospheres using a volume-limited photometric sample drawn from the $Gaia$ DR2 catalogue, the Sloan Digital Sky Survey (SDSS) and the $Galaxy~Evolution~Explorer$ (GALEX). We exploit the strength of the Balmer jump in hydrogen-atmosphere DA white dwarfs to separate them from helium-dominated objects in SDSS colour space. Across the effective temperature range from 20000K to 9000K we find that 22% of white dwarfs will undergo a spectral change, with no spectral evolution being ruled out at 5$\sigma$. The most likely explanation is that the increase in He-rich objects is caused by the convective mixing of DA stars with thin hydrogen layers, in which helium is dredged up from deeper layers by a surface hydrogen convection zone. The rate of change in the fraction of He-rich objects as a function of temperature, coupled with a recent grid of 3D radiation-hydrodynamic simulations of convective DA white dwarfs - which include the full overshoot region - lead to a discussion on the distribution of total hydrogen mass in white dwarfs. We find that 60% of white dwarfs must have a hydrogen mass larger than $M_{\rm H}/M_{\rm WD} = 10^{-10}$, another 25% have masses in the range $M_{\rm H}/M_{\rm WD} = 10^{-14}-10^{-10}$, and 15% have less hydrogen than $M_{\rm H}/M_{\rm WD} = 10^{-14}$. These results have implications for white dwarf asteroseismology, stellar evolution through the asymptotic giant branch (AGB) and accretion of planetesimals onto white dwarfs.

White dwarfs are stellar embers depleted of nuclear energy sources that cool over billions of years. These stars, which are supported by electron degeneracy pressure, reach densities of 1e7 grams per cubic centimetre in their cores. It has been predicted that a first-order phase transition occurs during white-dwarf cooling, leading to the crystallization of the non-degenerate carbon and oxygen ions in the core, which releases a considerable amount of latent heat and delays the cooling process by about one billion years. However, no direct observational evidence of this effect has been reported so far. Here we report the presence of a pile-up in the cooling sequence of evolving white dwarfs within 100 parsecs of the Sun, determined using photometry and parallax data from the Gaia satellite. Using modelling, we infer that this pile-up arises from the release of latent heat as the cores of the white dwarfs crystallize. In addition to the release of latent heat, we find strong evidence that cooling is further slowed by the liberation of gravitational energy from element sedimentation in the crystallizing cores. Our results describe the energy released by crystallization in strongly coupled Coulomb plasmas, and the measured cooling delays could help to improve the accuracy of methods used to determine the age of stellar populations from white dwarfs.

We present a theoretical description of macroscopic diffusion caused by convective overshoot in pure-hydrogen DA white dwarfs using three-dimensional (3D), closed-bottom, radiation hydrodynamics CO$^5$BOLD simulations. We rely on a new grid of deep 3D white dwarf models in the temperature range 11400 K $\leq T_{\mathrm{eff}} \leq$ 18000 K where tracer particles and a tracer density are used to derive macroscopic diffusion coefficients driven by convective overshoot. These diffusion coefficients are compared to microscopic diffusion coefficients from one-dimensional structures. We find that the mass of the fully mixed region is likely to increase by up to 2.5 orders of magnitude while inferred accretion rates increase by a more moderate order of magnitude. We present evidence that an increase in settling time of up to 2 orders of magnitude is to be expected which is of significance for time-variability studies of polluted white dwarfs. Our grid also provides the most robust constraint on the onset of convective instabilities in DA white dwarfs to be in the effective temperature range from 18000 to 18250 K.

Practically all known planet hosts will evolve into white dwarfs, and large parts of their planetary systems will survive this transition - the same is true for the solar system beyond the orbit of Mars. Spectroscopy of white dwarfs accreting planetary debris provides the most accurate insight into the bulk composition of exo-planets. Ground-based spectroscopic surveys of ~260, 000 white dwarfs detected with Gaia will identify >1000 evolved planetary systems, and high-throughput high-resolution space-based ultraviolet spectroscopy is essential to measure in detail their abundances. So far, evidence for two planetesimals orbiting closely around white dwarfs has been obtained, and their study provides important constraints on the composition and internal structure of these bodies. Major photometric and spectroscopic efforts will be necessary to assemble a sample of such close-in planetesimals that is sufficiently large to establish their properties as a population, and to deduce the architectures of the outer planetary systems from where they originated. Mid-infrared spectroscopy of the dusty disks will provide detailed mineralogical information of the debris, which, in combination with the elemental abundances measured from the white dwarf spectroscopy, will enable detailed physical modelling of the chemical, thermodynamic, and physical history of the accreted material. Flexible multi-epoch infrared observations are essential to determine the physical nature, and origin of the variability observed in many of the dusty disks. Finally, the direct detection of the outer reservoirs feeding material to the white dwarfs will require sensitive mid- and far-infrared capabilities.

We report on a comparison of spectroscopic analyses for hydrogen (DA) and helium atmosphere (DB) white dwarfs with Gaia Data Release 2 (DR2) parallaxes and photometry. We assume a reddening law and a mass-radius relation to connect the effective temperatures (Teff) and surface gravities (log g) to masses and radii. This allows the comparison of two largely independent sets of fundamental parameters for 7039 DA and 521 DB stars with high-quality observations. This subset of the Gaia white dwarf sample is large enough to detect systematic trends in the derived parameters. We find that spectroscopic and photometric parameters generally agree within uncertainties when the expectation of a single star is verified. Gaia allows the identification of a small systematic offset in the temperature scale between the two techniques, as well as confirming a small residual high-mass bump in the DA mass distribution around 11,000-13,000 K. This assessment of the accuracy of white dwarf fundamental parameters derived from Gaia is a first step in understanding systematic effects in related astrophysical applications such as the derivation of the local stellar formation history, initial-to-final mass relation, and statistics of evolved planetary systems.

We present a catalogue of white dwarf candidates selected from the second data release of $Gaia$ (DR2). We used a sample of spectroscopically confirmed white dwarfs from the Sloan Digital Sky Survey (SDSS) to map the entire space spanned by these objects in the $Gaia$ Hertzsprung-Russell diagram. We then defined a set of cuts in absolute magnitude, colour, and a number of $Gaia$ quality flags to remove the majority of contaminating objects. Finally, we adopt a method analogous to the one presented in our earlier SDSS photometric catalogues to calculate a probability of being a white dwarf ($P_{\mathrm{WD}}$) for all $Gaia$ sources which passed the initial selection. The final catalogue is composed of $486\,641$ stars with calculated $P_{\mathrm{WD}}$ from which it is possible to select a sample of $\simeq 260\,000$ high-confidence white dwarf candidates in the magnitude range $8<G<21$. By comparing this catalogue with a sample of SDSS white dwarf candidates we estimate an upper limit in completeness of $85$ per cent for white dwarfs with $G \leq 20$ mag and $T_{\mathrm{eff}} > 7000$ K, at high Galactic latitudes ($|b|>20^{\circ}$). However, the completeness drops at low Galactic latitudes, and the magnitude limit of the catalogue varies significantly across the sky as a function of $Gaia$'s scanning law. We also provide the list of objects within our sample with available SDSS spectroscopy. We use this spectroscopic sample to characterise the observed structure of the white dwarf distribution in the H-R diagram.

The expanding ejecta from a classical nova remains hot enough ($\sim10^{4}\, {\rm K}$) to be detected in thermal radio emission for up to years after the cessation of mass loss triggered by a thermonuclear instability on the underlying white dwarf (WD). Nebular spectroscopy of nova remnants confirms the hot temperatures observed in radio observations. During this same period, the unstable thermonuclear burning transitions to a prolonged period of stable burning of the remnant hydrogen-rich envelope, causing the WD to become, temporarily, a super-soft X-ray source. We show that photoionization heating of the expanding ejecta by the hot WD maintains the observed nearly constant temperature of $(1-4)\times10^4\mathrm{~K}$ for up to a year before an eventual decline in temperature due to either the cessation of the supersoft phase or the onset of a predominantly adiabatic expansion. We simulate the expanding ejecta using a one-zone model as well as the Cloudy spectral synthesis code, both incorporating the time-dependent WD effective temperatures for a range of masses from $0.60\ M_{\odot}$ to $1.10\ M_{\odot}$. We show that the duration of the nearly isothermal phase depends most strongly on the velocity and mass of the ejecta and that the ejecta temperature depends on the WD's effective temperature, and hence its mass.