The initial-final mass relation (IFMR) links the birth mass of a star to the mass of the compact remnant left at its death. While the relevance of the IFMR across astrophysics is universally acknowledged, not all of its fine details have yet been resolved. A new analysis of a few carbon-oxygen white dwarfs in old open clusters of the Milky Way led us to identify a kink in the IFMR, located over a range of initial masses, $1.65 \lesssim M_{\rm i}/M_{\odot} \lesssim 2.10$. The kink's peak in WD mass of $\approx 0.70-0.75 \, M_{\odot}$ is produced by stars with $M_{\rm i} \simeq 1.8 - 1.9 \, M_{\odot}$, corresponding to ages of about $1.8 - 1.7 $ Gyr. Interestingly, this peak coincides with the initial mass limit between low-mass stars that develop a degenerate helium core after central hydrogen exhaustion, and intermediate-mass stars that avoid electron degeneracy. We interpret the IFMR kink as the signature of carbon star formation in the Milky Way. This finding is critical to constraining the evolution and chemical enrichment of low-mass stars, and their impact on the spectrophotometric properties of galaxies.

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.

The mass of a star is the most fundamental parameter for its structure, evolution, and final fate. It is particularly important for any kind of stellar archaeology and characterization of exoplanets. There exists a variety of methods in astronomy to estimate or determine it. In this review we present a significant number of such methods, beginning with the most direct and model-independent approach using detached eclipsing binaries. We then move to more indirect and model-dependent methods, such as the quite commonly used isochrone or stellar track fitting. The arrival of quantitative asteroseismology has opened a completely new approach to determine stellar masses and to complement and improve the accuracy of other methods. We include methods for different evolutionary stages, from the pre-main sequence to evolved (super)giants and final remnants. For all methods uncertainties and restrictions will be discussed. We provide lists of altogether more than 200 benchmark stars with relative mass accuracies between $[0.3,2]\%$ for the covered mass range of $M\in [0.1,16]\,M_\odot$, $75\%$ of which are stars burning hydrogen in their core and the other $25\%$ covering all other evolved stages. We close with a recommendation how to combine various methods to arrive at a "mass-ladder" for stars.

We report the discovery of a white dwarf exhibiting deep, irregularly shaped transits, indicative of circumstellar planetary debris. Using Zwicky Transient Facility DR2 photometry of ZTF$\,$J013906.17+524536.89 and follow-up observations from the Las Cumbres Observatory, we identify multiple transit events that recur every ${\approx}\,107.2\,$d, much longer than the $4.5{-}4.9\,$h orbital periods observed in WD$\,$1145+017, the only other white dwarf known with transiting planetary debris. The transits vary in both depth and duration, lasting $15{-}25\,$d and reaching $20{-}45\,\%$ dips in flux. Optical spectra reveal strong Balmer lines, identifying the white dwarf as a DA with $T_{\mathrm{eff}}=10{,}530\pm140\,\mathrm{K}$ and $\log(g)=7.86\pm0.06$. A $\mathrm{Ca\,II\,K}$ absorption feature is present in all spectra both in and out of transit. Spectra obtained during one night at roughly $15\,\%$ transit depth show increased $\mathrm{Ca\,II\,K}$ absorption with a model atmospheric fit suggesting $[\mathrm{Ca/H}]=-4.6\pm0.3$, whereas spectra taken on three nights out of transit have $[\mathrm{Ca/H}]$ of -5.5, -5.3, and -4.9 with similar uncertainties. While the $\mathrm{Ca\,II\,K}$ line strength varies by only 2-sigma, we consider a predominantly interstellar origin for Ca absorption unlikely. We suggest a larger column density of circumstellar metallic gas along the line of site or increased accretion of material onto the white dwarf's surface are responsible for the Ca absorption, but further spectroscopic studies are required. In addition, high-speed time series photometry out of transit reveals variability with periods of 900 and 1030$\,$s, consistent with ZZ Ceti pulsations.

Context. We present our findings on 18 formerly known ZZ Ceti stars observed by the TESS space telescope in 120s cadence mode during the survey observation of the southern ecliptic hemisphere. Aims. We focus on the frequency analysis of the space-based observations, comparing the results with the findings of the previous ground-based measurements. The frequencies detected by the TESS observations can serve as inputs for future asteroseismic analyses. Methods. We performed standard pre-whitening of the data sets to derive the possible pulsation frequencies of the different targets. In some cases, we fitted Lorentzians to the frequency groups that emerged as the results of short-term amplitude/phase variations that occurred during the TESS observations. Results. We detected more than 40 pulsation frequencies in seven ZZ Ceti stars observed in the 120s cadence by TESS, with better than 0.1 microHz precision. We found that HE 0532-5605 may be a new outbursting ZZ Ceti. Ten targets do not show any significant pulsation frequencies in their Fourier transforms, due to a combination of their intrinsic faintness and/or crowding on the large TESS pixels. We also detected possible amplitude/phase variations during the TESS observations in some cases. Such behaviour in these targets was not previously identified from ground-based observations.

White dwarfs are dense, cooling stellar embers consisting mostly of carbon and oxygen, or oxygen and neon (with a few percent carbon) at higher initial stellar masses. These stellar cores are enveloped by a shell of helium which in turn is usually surrounded by a layer of hydrogen, generally prohibiting direct observation of the interior composition. However, carbon is observed at the surface of a sizeable fraction of white dwarfs, sometimes with traces of oxygen, and it is thought to be dredged-up from the core by a deep helium convection zone. In these objects only traces of hydrogen are found as large masses of hydrogen are predicted to inhibit hydrogen/helium convective mixing within the envelope. We report the identification of WDJ055134.612+413531.09, an ultra-massive (1.14 $M_\odot$) white dwarf with a unique hydrogen/carbon mixed atmosphere (C/H=0.15 in number ratio). Our analysis of the envelope and interior indicates that the total hydrogen and helium mass fractions must be several orders of magnitude lower than predictions of single star evolution: less than $10^{-9.5}$ and $10^{-7.0}$, respectively. Due to the fast kinematics ($129\pm5$ km/s relative to the local standard of rest), large mass, and peculiar envelope composition, we argue that WDJ0551+4135 is consistent with formation from the merger of two white dwarfs in a tight binary system.

We present the first discovery from the COol Companions ON Ultrawide orbiTS (COCONUTS) program, a large-scale survey for wide-orbit planetary and substellar companions. We have discovered a co-moving system COCONUTS-1, composed of a hydrogen-dominated white dwarf (PSO J058.9855+45.4184; $d=31.5$ pc) and a T4 companion (PSO J058.9869+45.4296) at a $40.6''$ (1280 au) projected separation. We derive physical properties for COCONUTS-1B from (1) its near-infrared spectrum using cloudless Sonora atmospheric models, and (2) its luminosity and the white dwarf's age ($7.3_{-1.6}^{+2.8}$ Gyr) using Sonora evolutionary models. The two methods give consistent temperatures and radii, but atmospheric models infer a lower surface gravity and therefore an unphysically young age. Assuming evolutionary model parameters ($T_{\rm eff}=1255^{+6}_{-8}$ K, $\log{g}=5.44^{+0.02}_{-0.03}$ dex, $R=0.789^{+0.011}_{-0.005}$ R$_{\rm Jup}$), we find cloudless model atmospheres have brighter Y- and J-band fluxes than the data, suggesting condensate clouds have not fully dispersed around 1300 K. The W2 flux (4.6 $\mu$m) of COCONUTS-1B is fainter than models, suggesting non-equilibrium mixing of CO. To investigate the gravity dependence of the L/T transition, we compile all 60 known L6-T6 benchmarks and derive a homogeneous set of temperatures, surface gravities, and masses. As is well-known, young, low-gravity late-L dwarfs have significantly fainter, redder near-infrared photometry and $\approx200-300$ K cooler temperatures than old, high-gravity objects. Our sample now reveals such gravity dependence becomes weaker for T dwarfs, with young objects having comparable near-infrared photometry and $\approx100$ K cooler temperatures compared to old objects. Finally, we find that young objects have a larger amplitude J-band brightening than old objects, and also brighten at H band as they cross the L/T transition.

Identifying planets around O-type and B-type stars is inherently difficult; the most massive known planet host has a mass of only about $3M_{\odot}$. However, planetary systems which survive the transformation of their host stars into white dwarfs can be detected via photospheric trace metals, circumstellar dusty and gaseous discs, and transits of planetary debris crossing our line-of-sight. These signatures offer the potential to explore the efficiency of planet formation for host stars with masses up to the core-collapse boundary at $\approx 8M_{\odot}$, a mass regime rarely investigated in planet formation theory. Here, we establish limits on where both major and minor planets must reside around $\approx 6M_{\odot}-8M_{\odot}$ stars in order to survive into the white dwarf phase. For this mass range, we find that intact terrestrial or giant planets need to leave the main sequence beyond approximate minimum star-planet separations of respectively about 3 and 6 au. In these systems, rubble pile minor planets of radii 10, 1.0, and 0.1 km would have been shorn apart by giant branch radiative YORP spin-up if they formed and remained within, respectively, tens, hundreds and thousands of au. These boundary values would help distinguish the nature of the progenitor of metal-pollution in white dwarf atmospheres. We find that planet formation around the highest mass white dwarf progenitors may be feasible, and hence encourage both dedicated planet formation investigations for these systems and spectroscopic analyses of the highest mass white dwarfs.

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.

The main objective of the present work is to extend these investigations by computing the gravity and limb-darkening coefficients for white dwarf atmosphere models with hydrogen, helium, or mixed compositions (types DA, DB, and DBA). We computed gravity and limb-darkening coefficients for DA, DB, and DBA white dwarfs atmosphere models, covering the transmission curves of the Sloan, UBVRI, Kepler, TESS, and Gaia photometric systems. Specific calculations for the HiPERCAM instrument were also carried out. For all calculations of the limb-darkening coefficients we used the least-squares method. Concerning the effects of tidal and rotational distortions, we also computed for the first time the gravity-darkening coefficients $y(\lambda)$ for white dwarfs using the same models of stellar atmospheres as in the case of limb-darkening. A more general differential equation was introduced to derive these quantities, including the partial derivative $\left(\partial{\ln I_o(\lambda)}/{\partial{\ln g}}\right)_{T_{\rm eff}}$. Six laws were adopted to describe the specific intensity distribution: linear, quadratic, square root, logarithmic, power-2, and a more general one with four coefficients. The computations are presented for the chemical compositions log[H/He] = $-$10.0 (DB), $-$2.0 (DBA) and He/H = 0 (DA), with log g varying between 5.0 and 9.5 and effective temperatures between 3750 K-100\,000 K. For effective temperatures higher than 40\,000 K, the models were also computed adopting nonlocal thermal equilibirum (DA). The adopted mixing-length parameters are ML2/$\alpha = 0.8$ (DA case) and 1.25 (DB and DBA). The results are presented in the form of 112 tables. Additional calculations, such as for other photometric systems and/or different values of log[H/He], $\log g,$ and T$_{\rm eff}$ can be performed upon request.

We report the discovery of a detached double white dwarf binary with an orbital period of $\approx20.6\,\rm minutes$, PTF J053332.05+020911.6. The visible object in this binary, PTF J0533+0209B, is a $\approx 0.17\,M_\odot$ mass white dwarf with a helium-dominated atmosphere containing traces of hydrogen (DBA). This object exhibits ellipsoidal variations due to tidal deformation, and is the visible component in a single-lined spectroscopic binary with a velocity semi-amplitude of $K_B=618.7\pm6.9 \, \rm km \, s^{-1}$. We have detected significant orbital decay due to the emission of gravitational radiation, and we expect that the Laser Interferometer Space Antenna (LISA) will detect this system with a signal to noise of $8.4^{+4.2}_{-3.0}$ after four years of operation. Because this system already has a well determined orbital period, radial velocity semi-amplitude, temperature, atmospheric composition, surface gravity, and orbital decay rate, a LISA signal will help fully constrain the properties of this system by providing a direct measurement of its inclination. Thus, this binary demonstrates the synergy between electromagnetic and gravitational radiation for constraining the physical properties of an astrophysical object.

In the era of modern digital sky surveys, uncertainties in the flux of stellar standards are commonly the dominant systematic error in photometric calibration and can often affect the results of higher-level experiments. The Hubble Space Telescope (HST) spectrophotometry, which is based on computed model atmospheres for three hot (Teff > 30,000 K) pure-hydrogen (DA) white dwarfs, is currently considered the most reliable and internally consistent flux calibration. However many next generation facilities (e.g. Harmoni on E-ELT, Euclid and JWST) will focus on IR observations, a regime in which white dwarf calibration has not yet been robustly tested. Cool DA white dwarfs have energy distributions that peak close to the optical or near-IR, do not have shortcomings from UV metal line blanketing, and have a reasonably large sky density ($\simeq$ 4 deg$^{-2}$ at G < 20), making them, potentially, excellent calibrators. Here we present a pilot study based on STIS+WFC3 observations of two bright DA white dwarfs to test whether targets cooler than current hot primary standards (Teff < 20,000 K) are consistent with the HST flux scale. We also test the robustness of white dwarf models in the IR regime from an X-shooter analysis of Paschen lines and by cross-matching our previously derived Gaia white dwarf catalogue with observations obtained with 2MASS, UKIDSS, VHS, and WISE.

We perform a calibration of the mixing length parameter at the bottom boundary of the convection zone for helium-dominated atmospheres of white dwarfs. This calibration is based on a grid of 3D DB (pure-helium) and DBA (helium-dominated with traces of hydrogen) model atmospheres computed with the CO5BOLD code, and a grid of 1D DB and DBA envelope structures. The 3D models span a parameter space of hydrogen-to-helium abundances between -10.0 and -2.0, surface gravities between 7.5 and 9.0 and effective temperatures between 12000 K and 34000 K. The 1D envelopes cover a similar atmospheric parameter range, but are also calculated with different values of the mixing length parameter, namely ML2/alpha between 0.4 and 1.4. The calibration is performed based on two definitions of the bottom boundary of the convection zone, the Schwarzschild and the zero convective flux boundaries. Thus, our calibration is relevant for applications involving the bulk properties of the convection zone including its total mass, which excludes the spectroscopic technique. Overall, the calibrated ML2/alpha is smaller than what is commonly used in evolutionary models and theoretical determinations of the blue edge of the instability strip for pulsating DB and DBA stars. With calibrated ML2/alpha we are able to deduce more accurate convection zone sizes needed for studies of planetary debris mixing and dredge-up of carbon from the core. We highlight this by calculating examples of metal-rich 3D DBAZ models and finding their convection zone masses. Mixing length calibration represents the first step of in-depth investigations of convective overshoot in white dwarfs with helium-dominated atmospheres.

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 the photometric calibration of the twelve optical passbands observed by the Javalambre Photometric Local Universe Survey (J-PLUS). The proposed calibration method has four steps: (i) definition of a high-quality set of calibration stars using Gaia information and available 3D dust maps; (ii) anchoring of the J-PLUS gri passbands to the Pan-STARRS photometric solution, accounting for the variation of the calibration with the position of the sources on the CCD; (iii) homogenization of the photometry in the other nine J-PLUS filters using the dust de-reddened instrumental stellar locus in (X - r) versus (g - i) colours, where X is the filter to calibrate. The zero point variation along the CCD in these filters was estimated with the distance to the stellar locus. Finally, (iv) the absolute colour calibration was obtained with the white dwarf locus. We performed a joint Bayesian modelling of eleven J-PLUS colour-colour diagrams using the theoretical white dwarf locus as reference. This provides the needed offsets to transform instrumental magnitudes to calibrated magnitudes outside the atmosphere. The uncertainty of the J-PLUS photometric calibration, estimated from duplicated objects observed in adjacent pointings and accounting for the absolute colour and flux calibration errors, are ~19 mmag in u, J0378 and J0395, ~11 mmag in J0410 and J0430, and ~8 mmag in g, J0515, r, J0660, i, J0861, and z. We present an optimized calibration method for the large area multi-filter J-PLUS project, reaching 1-2% accuracy within an area of 1 022 square degrees without the need for long observing calibration campaigns or constant atmospheric monitoring. The proposed method will be adapted for the photometric calibration of J-PAS, that will observe several thousand square degrees with 56 narrow optical filters.

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.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

Observational evidence of white dwarf planetary systems is dominated by the remains of exo-asteroids through accreted metals, debris discs, and orbiting planetesimals. However, exo-planets in these systems play crucial roles as perturbing agents, and can themselves be perturbed close to the white dwarf Roche radius. Here, we illustrate a procedure for computing the tidal interaction between a white dwarf and a near-spherical solid planet. This method determines the planet's inward and/or outward drift, and whether the planet will reach the Roche radius and be destroyed. We avoid constant tidal lag formulations and instead employ the self-consistent secular Darwin-Kaula expansions from Bou\'{e} & Efroimsky (2019), which feature an arbitrary frequency dependence on the quality functions. We adopt wide ranges of dynamic viscosities and spin rates for the planet in order to straddle many possible outcomes, and provide a foundation for the future study of individual systems with known or assumed rheologies. We find that: (i) massive Super-Earths are destroyed more readily than minor planets (such as the ones orbiting WD 1145+017 and SDSS J1228+1040), (ii) low-viscosity planets are destroyed more easily than high-viscosity planets, and (iii) the boundary between survival and destruction is likely to be fractal and chaotic.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

The mechanisms of the formation and evolution of the Milky Way are encoded in the orbits, chemistry and ages of its stars. With the 4MOST MIlky way Disk And BuLgE Low-Resolution Survey (4MIDABLE-LR) we aim to study kinematic and chemical substructures in the Milky Way disc and bulge region with samples of unprecedented size out to larger distances and greater precision than conceivable with Gaia alone or any other ongoing or planned survey. Gaia gives us the unique opportunity for target selection based almost entirely on parallax and magnitude range, hence increasing the efficiency in sampling larger Milky Way volumes with well-defined and effective selection functions. Our main goal is to provide a detailed chrono-chemo-kinematical extended map of our Galaxy and the largest Gaia follow-up down to $G = 19$ magnitudes (Vega). The complex nature of the disc components (for example, large target densities and highly structured extinction distribution in the Milky Way bulge and disc area), prompted us to develop a survey strategy with five main sub-surveys that are tailored to answer the still open questions about the assembly and evolution of our Galaxy, while taking full advantage of the Gaia data.

The Maunakea Spectroscopic Explorer (MSE) is a planned 11.25-m aperture facility with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. A rebirth of the 3.6m Canada-France-Hawaii Telescope on Maunakea, MSE will use 4332 fibers operating at three different resolving powers (R ~ 2500, 6000, 40000) across a wavelength range of 0.36-1.8mum, with dynamical fiber positioning that allows fibers to match the exposure times of individual objects. MSE will enable spectroscopic surveys with unprecedented scale and sensitivity by collecting millions of spectra per year down to limiting magnitudes of g ~ 20-24 mag, with a nominal velocity precision of ~100 m/s in high-resolution mode. This white paper describes science cases for stellar astrophysics and exoplanet science using MSE, including the discovery and atmospheric characterization of exoplanets and substellar objects, stellar physics with star clusters, asteroseismology of solar-like oscillators and opacity-driven pulsators, studies of stellar rotation, activity, and multiplicity, as well as the chemical characterization of AGB and extremely metal-poor stars.

The semi-empirical initial-final mass relation (IFMR) connects spectroscopically analyzed white dwarfs in star clusters to the initial masses of the stars that formed them. Most current stellar evolution models, however, predict that stars will evolve to white dwarfs $\sim$0.1 M$_\odot$ less massive than that found in the IFMR. We first look at how varying theoretical mass-loss rates, third dredge-up efficiencies, and convective-core overshoot may help explain the differences between models and observations. These parameters play an important role at the lowest masses (M$_{\rm initial}$ $<$ 3 M$_\odot$). At higher masses, only convective-core overshoot meaningfully affects white dwarf mass, but alone it likely cannot explain the observed white dwarf masses nor why the IFMR scatter is larger than observational errors predict. These higher masses, however, are also where rotational mixing in main sequence stars begins to create more massive cores, and hence more massive white dwarfs. This rotational mixing also extends a star's lifetime, making faster rotating progenitors appear like less massive stars in their semi-empirical age analysis. Applying the observed range of young B-dwarf rotations to the MIST or SYCLIST rotational models demonstrates a marked improvement in reproducing both the observed IFMR data and its scatter. The incorporation of both rotation and efficient convective-core overshoot significantly improves the match with observations. This work shows that the IFMR provides a valuable observational constraint on how rotation and convective-core overshoot affect the core evolution of a star.

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.

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.

We present the initial-final mass relation (IFMR) based on the self-consistent analysis of Sirius B and 79 white dwarfs from 13 star clusters. We have also acquired additional signal on eight white dwarfs previously analyzed in the NGC 2099 cluster field, four of which are consistent with membership. These reobserved white dwarfs have masses ranging from 0.72 to 0.97 M$_\odot$, with initial masses from 3.0 to 3.65 M$_\odot$, where the IFMR has an important change in slope that these new data help to observationally confirm. In total, this directly measured IFMR has small scatter ($\sigma$ = 0.06 M$_\odot$) and spans from progenitors of 0.85 to 7.5 M$_\odot$. Applying two different stellar evolutionary models to infer two different sets of white dwarf progenitor masses shows that when the same model is also used to derive the cluster ages, the resulting IFMR has weak sensitivity to the adopted model at all but the highest initial masses ($>$5.5 M$_\odot$). The non-linearity of the IFMR is also clearly observed with moderate slopes at lower masses (0.08 M$_{\rm final}$/M$_{\rm initial}$) and higher masses (0.11 M$_{\rm final}$/M$_{\rm initial}$) that are broken up by a steep slope (0.19 M$_{\rm final}$/M$_{\rm initial}$) between progenitors from 2.85 to 3.6 M$_\odot$. This IFMR shows total stellar mass loss ranges from 33\% of M$_{\rm initial}$ at 0.83 M$_\odot$ to 83\% of M$_{\rm initial}$ at 7.5 M$_\odot$. Testing this total mass loss for dependence on progenitor metallicity, however, finds no detectable sensitivity across the moderate range of -0.15 $<$ [Fe/H] $<$ +0.15.

We present the first grid of 3D simulations for the pure-helium atmospheres of DB white dwarfs. The simulations were computed with the CO5BOLD radiation-hydrodynamics code and cover effective temperatures and surface gravities between 12000 K < Teff < 34000 K and 7.5 < log(g) (cgs units) < 9.0, respectively. In this introductory work, synthetic spectra calculated from the 3D simulations are compared to appropriate 1D model spectra under a differential approach. This results in the derivation of 3D corrections for the spectroscopically derived atmospheric parameters of DB stars with respect to the 1D ML2/alpha = 1.25 mixing-length parameterisation. No significant Teff corrections are found for the V777 Her instability strip region, and therefore no 3D revision is expected for the empirical blue and red edges of the strip. However, large log(g) corrections are found in the range 12000 K < Teff < 23000 K for all log(g) values covered by the 3D grid. These corrections indicate that 1D model atmospheres overpredict log(g), reminiscent of the results found from 3D simulations of pure-hydrogen white dwarfs. The next step will be to compute 3D simulations with mixed helium and hydrogen atmospheres to comprehend the full implications for the stellar parameters of DB and DBA white dwarfs.

WD 1145+017 is currently the only white dwarf known to exhibit periodic transits of planetary debris as well as absorption lines from circumstellar gas. We present the first simultaneous fast optical spectrophotometry and broad-band photometry of the system, obtained with the Gran Telescopio Canarias (GTC) and the Liverpool Telescope (LT), respectively. The observations spanned $5.5$ h, somewhat longer than the $4.5$-h orbital period of the debris. Dividing the GTC spectrophotometry into five wavelength bands reveals no significant colour differences, confirming grey transits in the optical. We argue that absorption by an optically thick structure is a plausible alternative explanation for the achromatic nature of the transits that can allow the presence of small-sized ($\sim\mu$m) particles. The longest ($87$ min) and deepest ($50$ per cent attenuation) transit recorded in our data exhibits a complex structure around minimum light that can be well modelled by multiple overlapping dust clouds. The strongest circumstellar absorption line, Fe II $\lambda$5169, significantly weakens during this transit, with its equivalent width reducing from a mean out-of-transit value of $2$ \AA\ to $1$ \AA\ in-transit, supporting spatial correlation between the circumstellar gas and dust. Finally, we made use of the Gaia Data Release 2 and archival photometry to determine the white dwarf parameters. Adopting a helium-dominated atmosphere containing traces of hydrogen and metals, and a reddening $E(B-V)=0.01$ we find $T_\mathrm{eff}=15\,020 \pm 520$ K, $\log g=8.07\pm0.07$, corresponding to $M_\mathrm{WD}=0.63\pm0.05\ \mbox{$\mathrm{M}_{\odot}$}$ and a cooling age of $224\pm30$ Myr.

Using Gaia DR2 data, we present an up-to-date sample of white dwarfs within 20 pc of the Sun. In total we identified 139 systems in Gaia DR2, nine of which are new detections, with the closest of these located at a distance of 13.05 pc. We estimated atmospheric parameters for all stellar remnants based on the Gaia parallaxes and photometry. The high precision and completeness of the Gaia astrometry allowed us to search for wide binary companions. We re-identified all known binaries where both components have accurate DR2 astrometry, and established the binarity of one of the nine newly identified white dwarfs. No new companions were found to previously known 20 pc white dwarfs. Finally, we estimated the local white dwarf space-density to be $(4.49\pm0.38)\times10^{-3}$ pc$^{-3}$, having given careful consideration to the distance-dependent Gaia completeness, which misses known objects at short distances, but is close to complete for white dwarfs near 20 pc.

With typical periods of order 10 minutes, the pulsation signatures of ZZ Ceti variables (pulsating hydrogen-atmosphere white dwarf stars) are severely undersampled by long-cadence (29.42 minutes per exposure) K2 observations. Nyquist aliasing renders the intrinsic frequencies ambiguous, stifling precision asteroseismology. We report the discovery of two new ZZ Cetis in long-cadence K2 data: EPIC 210377280 and EPIC 220274129. Guided by 3-4 nights of follow-up, high-speed (<=30 s) photometry from McDonald Observatory, we recover accurate pulsation frequencies for K2 signals that reflected 4-5 times off the Nyquist with the full precision of over 70 days of monitoring (~0.01 muHz). In turn, the K2 observations enable us to select the correct peaks from the alias structure of the ground-based signals caused by gaps in the observations. We identify at least seven independent pulsation modes in the light curves of each of these stars. For EPIC 220274129, we detect three complete sets of rotationally split ell=1 (dipole mode) triplets, which we use to asteroseismically infer the stellar rotation period of 12.7+/-1.3 hr. We also detect two sub-Nyquist K2 signals that are likely combination (difference) frequencies. We attribute our inability to match some of the K2 signals to the ground-based data to changes in pulsation amplitudes between epochs of observation. Model fits to SOAR spectroscopy place both EPIC 210377280 and EPIC 220274129 near the middle of the ZZ Ceti instability strip, with Teff = 11590+/-200 K and 11810+/-210 K, and masses 0.57+/-0.03 Msun and 0.62+/-0.03 Msun, respectively.

The intriguing WN4b star WR6 has been known to display epoch-dependent spectroscopic, photometric and polarimetric variability for several decades. In this paper, we set out to verify if a simplified analytical model in which Corotating Interaction Regions (CIRs) threading an otherwise spherical wind is able to reproduce the many broadband continuum light curves from the literature with a reasonable set of parameters. We modified the optically thin model we developed in Ignace, St-Louis & Proulx-Giraldeau (2015) to approximately account for multiple scattering and used it to fit 13 separate datasets of this star. By including two CIRs in the wind, we obtained reasonable fits for all datasets with coherent values for the inclination of the rotation axis ($i_0=166^{\circ}$) and for its orientation in the plane of the sky, although in the latter case we obtained two equally acceptable values ($\psi=63^{\circ}$ and $\psi=152^{\circ}$) from the polarimetry. Additional line profile variation simulations using the Sobolev approximation for the line transfer allowed us to eliminate the $\psi=152^{\circ}$ solution. With the adopted configuration ($i_0=166^{\circ}$ and $\psi=63^{\circ}$), we were able to reproduce all datasets relatively well with two CIRs located near the stellar equator and always separated by $\sim 90^{\circ}$ in longitude. The epoch-dependency comes from the fact that these CIRs migrate along the surface of the star. Density contrasts smaller than a factor of two and large opening angles for the CIR ($\beta \geq 35^{\circ}$) were found to best reproduce the type of spectroscopic variability reported in the literature.

Using images from the Hubble Space Telescope Advanced Camera for Surveys, we measure the rate of cooling of white dwarfs in the globular cluster 47 Tucanae and compare it to modelled cooling curves. We examine the effects of the outer convective envelope reaching the nearly isothermal degenerate core and the release of latent heat during core crystallisation on the white dwarf cooling rates. For white dwarfs typical of 47 Tuc, the onset of these effects occur at similar times. The latent heat released during crystallisation is a small heat source. In contrast, the heat reservoir of the degenerate core is substantially larger. When the convective envelope reaches the nearly isothermal interior of the white dwarf, the star becomes brighter than it would be in the absence of this effect. Our modelled cooling curves that include this convective coupling closely match the observed luminosity function of the white dwarfs in 47 Tuc.

Around 10% of white dwarfs exhibit global magnetic structures with fields ranging from 1 kG to hundreds of MG. Recently, the first radiation magnetohydrodynamics simulations of the atmosphere of white dwarfs showed that convection should be suppressed in their photospheres for magnetic fields with strengths B $\gtrsim$ 50 kG. These predictions are in agreement with our knowledge of stellar physics (e.g. energy transfer in strong magnetic field regions of the solar photosphere), but have yet to be directly confirmed from white dwarf observations. We obtained COS far-UV spectroscopy of the weakly magnetic, hydrogen-atmosphere, white dwarf WD2105-820 and of three additional non-magnetic, convective remnants (all in the $T_{\mathrm{eff}}$ range 9000-11,000 K). We fitted both the COS and the already available optical spectra with convective and radiative atmospheric models. As expected, we find that for two of the non-magnetic comparison stars only convective model fits predict consistent $T_{\mathrm{eff}}$ values from both the optical and the FUV spectra. In contrast, for WD2105-820 only the best fitting radiative model produced consistent results.

We present photometry and spectroscopy for 27 pulsating hydrogen-atmosphere white dwarfs (DAVs, a.k.a. ZZ Ceti stars) observed by the Kepler space telescope up to K2 Campaign 8, an extensive compilation of observations with unprecedented duration (>75 days) and duty cycle (>90%). The space-based photometry reveals pulsation properties previously inaccessible to ground-based observations. We observe a sharp dichotomy in oscillation mode linewidths at roughly 800 s, such that white dwarf pulsations with periods exceeding 800 s have substantially broader mode linewidths, more reminiscent of a damped harmonic oscillator than a heat-driven pulsator. Extended Kepler coverage also permits extensive mode identification: We identify the spherical degree of 61 out of 154 unique radial orders, providing direct constraints of the rotation period for 20 of these 27 DAVs, more than doubling the number of white dwarfs with rotation periods determined via asteroseismology. We also obtain spectroscopy from 4m-class telescopes for all DAVs with Kepler photometry. Using these homogeneously analyzed spectra we estimate the overall mass of all 27 DAVs, which allows us to measure white dwarf rotation as a function of mass, constraining the endpoints of angular momentum in low- and intermediate-mass stars. We find that 0.51-to-0.73-solar-mass white dwarfs, which evolved from 1.7-to-3.0-solar-mass ZAMS progenitors, have a mean rotation period of 35 hr with a standard deviation of 28 hr, with notable exceptions for higher-mass white dwarfs. Finally, we announce an online repository for our Kepler data and follow-up spectroscopy, which we collect at http://www.k2wd.org.

The upcoming NASA Transiting Exoplanet Survey Satellite (TESS) will obtain space-based uninterrupted light curves for a large sample of bright white dwarfs distributed across the entire sky, providing a very rich resource for asteroseismological studies and the search for transits from planetary debris. We have compiled an all-sky catalogue of ultraviolet, optical, and infrared photometry as well as proper motions, which we propose as an essential tool for the preliminary identification and characterisation of potential targets. We present data for 1864 known white dwarfs and 305 high-probability white dwarf candidates brighter than 17 mag. We describe the spectroscopic follow-up of 135 stars, of which 82 are white dwarfs and 25 are hot subdwarfs. The new confirmed stars include six pulsating white dwarf candidates (ZZ Cetis), and nine white dwarf binaries with a cool main-sequence companion. We identify one star with a spectroscopic distance of only 25 pc from the Sun. Around the time TESS is launched, we foresee that all white dwarfs in this sample will have trigonometric parallaxes measured by the ESA Gaia mission next year.

Using data from the extended Kepler mission in K2 Campaign 10 we identify two eclipsing binaries containing white dwarfs with cool companions that have extremely short orbital periods of only 71.2 min (SDSS J1205-0242, a.k.a. EPIC 201283111) and 72.5 min (SDSS J1231+0041, a.k.a. EPIC 248368963). Despite their short periods, both systems are detached with small, low-mass companions, in one case a brown dwarf, and the other case either a brown dwarf or a low-mass star. We present follow-up photometry and spectroscopy of both binaries, as well as phase-resolved spectroscopy of the brighter system, and use these data to place preliminary estimates on the physical and binary parameters. SDSS J1205-0242 is composed of a $0.39\pm0.02$M$_\odot$ helium-core white dwarf which is totally eclipsed by a $0.049\pm0.006$M$_\odot$ ($51\pm6$M$_J$) brown dwarf companion, while SDSS J1231+0041 is composed of a $0.56\pm0.07$M$_\odot$ white dwarf which is partially eclipsed by a companion of mass $\lesssim 0.095$M$_\odot$. In the case of SDSS J1205-0242 we look at the combined constraints from common-envelope evolution and brown dwarf models; the system is compatible with similar constraints from other post common-envelope binaries given the current parameter uncertainties, but has potential for future refinement.

Using patterns in the oscillation frequencies of a white dwarf observed by K2, we have measured the fastest rotation rate, 1.13(02) hr, of any isolated pulsating white dwarf known to date. Balmer-line fits to follow-up spectroscopy from the SOAR telescope show that the star (SDSSJ0837+1856, EPIC 211914185) is a 13,590(340) K, 0.87(03) solar-mass white dwarf. This is the highest mass measured for any pulsating white dwarf with known rotation, suggesting a possible link between high mass and fast rotation. If it is the product of single-star evolution, its progenitor was a roughly 4.0 solar-mass main-sequence B star; we know very little about the angular momentum evolution of such intermediate-mass stars. We explore the possibility that this rapidly rotating white dwarf is the byproduct of a binary merger, which we conclude is unlikely given the pulsation periods observed.

An observational constraint on the contribution of double degenerates to Type Ia supernovae requires multiple radial velocity measurements of ideally thousands of white dwarfs. This is because only a small fraction of the double degenerate population is massive enough, with orbital periods short enough, to be considered viable Type Ia progenitors. We show how the radial velocity information available from public surveys such as the Sloan Digital Sky Survey can be used to pre-select targets for variability, leading to a ten-fold reduction in observing time required compared to an unranked or random survey. We carry out Monte Carlo simulations to quantify the detection probability of various types of binaries in the survey and show that this method, even in the most pessimistic case, doubles the survey size of the largest survey to date (the SPY survey) in less than 15 per cent of the required observing time. Our initial follow-up observations corroborate the method, yielding 15 binaries so far (eight known and seven new), as well as orbital periods for four of the new binaries.

The broadening of the hydrogen lines during flares is thought to result from increased charge (electron, proton) density in the flare chromosphere. However, disagreements between theory and modeling prescriptions have precluded an accurate diagnostic of the degree of ionization and compression resulting from flare heating in the chromosphere. To resolve this issue, we have incorporated the unified theory of electric pressure broadening of the hydrogen lines into the non-LTE radiative transfer code RH. This broadening prescription produces a much more realistic spectrum of the quiescent, A0 star Vega compared to the analytic approximations used as a damping parameter in the Voigt profiles. We test recent radiative-hydrodynamic (RHD) simulations of the atmospheric response to high nonthermal electron beam fluxes with the new broadening prescription and find that the Balmer lines are over-broadened at the densest times in the simulations. Adding many simultaneously heated and cooling model loops as a "multithread" model improves the agreement with the observations. We revisit the three-component phenomenological flare model of the YZ CMi Megaflare using recent and new RHD models. The evolution of the broadening, line flux ratios, and continuum flux ratios are well-reproduced by a multithread model with high-flux nonthermal electron beam heating, an extended decay phase model, and a "hot spot" atmosphere heated by an ultrarelativistic electron beam with reasonable filling factors: 0.1%, 1%, and 0.1% of the visible stellar hemisphere, respectively. The new modeling motivates future work to understand the origin of the extended gradual phase emission.

The Gaia Data Release 1 (DR1) sample of white dwarf parallaxes is presented, including 6 directly observed degenerates and 46 white dwarfs in wide binaries. This data set is combined with spectroscopic atmospheric parameters to study the white dwarf mass-radius relationship (MRR). Gaia parallaxes and G magnitudes are used to derive model atmosphere dependent white dwarf radii, which can then be compared to the predictions of a theoretical MRR. We find a good agreement between Gaia DR1 parallaxes, published effective temperatures (Teff) and surface gravities (log g), and theoretical MRRs. As it was the case for Hipparcos, the precision of the data does not allow for the characterisation of hydrogen envelope masses. The uncertainties on the spectroscopic atmospheric parameters are found to dominate the error budget and current error estimates for well-known and bright white dwarfs may be slightly optimistic. With the much larger Gaia DR2 white dwarf sample it will be possible to explore the MRR over a much wider range of mass, Teff, and spectral types.

The unprecedented extent of coverage provided by Kepler observations recently revealed outbursts in two hydrogen-atmosphere pulsating white dwarfs (DAVs) that cause hours-long increases in the overall mean flux of up to 14%. We have identified two new outbursting pulsating white dwarfs in K2, bringing the total number of known outbursting white dwarfs to four. EPIC 211629697, with T_eff = 10,780 +/- 140 K and log(g) = 7.94 +/- 0.08, shows outbursts recurring on average every 5.0 d, increasing the overall flux by up to 15%. EPIC 229227292, with T_eff = 11,190 +/- 170 K and log(g) = 8.02 +/- 0.05, has outbursts that recur roughly every 2.4 d with amplitudes up to 9%. We establish that only the coolest pulsating white dwarfs within a small temperature range near the cool, red edge of the DAV instability strip exhibit these outbursts.

We revisit the properties and astrophysical implications of the field white dwarf mass distribution in preparation of Gaia applications. Our study is based on the two samples with the best established completeness and most precise atmospheric parameters, the volume-complete survey within 20 pc and the Sloan Digital Sky Survey (SDSS) magnitude-limited sample. We explore the modelling of the observed mass distributions with Monte Carlo simulations, but find that it is difficult to constrain independently the initial mass function (IMF), the initial-to-final-mass relation (IFMR), the stellar formation history (SFH), the variation of the Galactic disk vertical scale height as a function of stellar age, and binary evolution. Each of these input ingredients has a moderate effect on the predicted mass distributions, and we must also take into account biases owing to unidentified faint objects (20 pc sample), as well as unknown masses for magnetic white dwarfs and spectroscopic calibration issues (SDSS sample). Nevertheless, we find that fixed standard assumptions for the above parameters result in predicted mean masses that are in good qualitative agreement with the observed values. It suggests that derived masses for both studied samples are consistent with our current knowledge of stellar and Galactic evolution. Our simulations overpredict by 40-50% the number of massive white dwarfs (M > 0.75 Msun) for both surveys, although we can not exclude a Salpeter IMF when we account for all biases. Furthermore, we find no evidence of a population of double white dwarf mergers in the observed mass distributions.

With the Keck I Low-Resolution Imaging Spectrometer we have observed nine white dwarf candidates in the very rich open cluster NGC 2099 (M37). The spectroscopy shows seven to be DA white dwarfs, one to be a DB white dwarf, and one to be a DZ white dwarf. Three of these DA white dwarfs are consistent with singly evolved cluster membership: an ultramassive (1.28$^{+0.05}_{-0.08}$ M$_\odot$) and two intermediate-mass (0.70 and 0.75 M$_\odot$) white dwarfs. Analysis of their cooling ages allows us to calculate their progenitor masses and establish new constraints on the initial-final mass relation. The intermediate-mass white dwarfs are in strong agreement with previous work over this mass regime. The ultramassive white dwarf has $V$ = 24.5, $\sim$2 mag fainter than the other two remnants. The spectrum of this star has lower quality, so the derived stellar properties (e.g., T$_{\rm eff}$, log g) have uncertainties that are several times higher than the brighter counterparts. We measure these uncertainties and establish the star's final mass as the highest-mass white dwarf discovered thus far in a cluster, but we are unable to calculate its progenitor mass because at this high mass and cooler T$_{\rm eff}$ its inferred cooling age is highly sensitive to its mass. At the highest temperatures, however, this sensitivity of cooling age to an ultramassive white dwarf's mass is only moderate. This demonstrates that future investigations of the upper-mass end of the initial-final mass relation must identify massive, newly formed white dwarfs (i.e., in young clusters with ages 50-150 Myr).

We present a robust statistical analysis of the white dwarf cooling sequence in 47 Tucanae. We combine HST UV and optical data in the core of the cluster, Modules for Experiments in Stellar Evolution (MESA) white dwarf cooling models, white dwarf atmosphere models, artificial star tests, and a Markov Chain Monte Carlo (MCMC) sampling method to fit white dwarf cooling models to our data directly. We use a technique known as the unbinned maximum likelihood to fit these models to our data without binning. We use these data to constrain neutrino production and the thickness of the hydrogen layer in these white dwarfs. The data prefer thicker hydrogen layers $(q_\mathrm{H}=3.2\e{-5})$ and we can strongly rule out thin layers $(q_\mathrm{H}=10^{-6})$. The neutrino rates currently in the models are consistent with the data. This analysis does not provide a constraint on the number of neutrino species.

We have observed a sample of 10 white dwarf candidates in the rich open cluster NGC 2323 (M50) with the Keck Low-Resolution Imaging Spectrometer. The spectroscopy shows eight to be DA white dwarfs, with six of these having high S/N appropriate for our analysis. Two of these white dwarfs are consistent with singly evolved cluster membership, and both are high mass $\sim$1.07 M$_\odot$, and give equivalent progenitor masses of 4.69 M$_\odot$. To supplement these new high-mass white dwarfs and analyze the initial-final mass relation (IFMR), we have also looked at 30 white dwarfs from publicly available data that are mostly all high-mass ($\gtrsim$0.9 M$_\odot$). These original published data exhibited significant scatter, and to test if this scatter is true or simply the result of systematics, we have uniformly analyzed the white dwarf spectra and have adopted thorough photometric techniques to derive uniform cluster parameters for their parent clusters. The resulting IFMR scatter is significantly reduced, arguing that mass-loss rates are not stochastic in nature and that within the ranges of metallicity and mass analyzed in this work mass loss is not highly sensitive to variations in metallicity. Lastly, when adopting cluster ages based on Y$^2$ isochrones, the slope of the high-mass IFMR remains steep and consistent with that found from intermediate-mass white dwarfs, giving a linear IFMR from progenitor masses between 3 to 6.5 M$_\odot$. In contrast, when adopting the slightly younger cluster ages based on PARSEC isochrones, the high-mass IFMR has a moderate turnover near an initial mass of 4 M$_\odot$.

We investigated the prospects for systematic searches of white dwarfs at low Galactic latitudes, using the VLT Survey Telescope (VST) H$\alpha$ Photometric Survey of the Galactic plane and Bulge (VPHAS+). We targeted 17 white dwarf candidates along sightlines of known open clusters, aiming to identify potential cluster members. We confirmed all the 17 white dwarf candidates from blue/optical spectroscopy, and we suggest five of them to be likely cluster members. We estimated progenitor ages and masses for the candidate cluster members, and compared our findings to those for other cluster white dwarfs. A white dwarf in NGC 3532 is the most massive known cluster member (1.13 M$_{\odot}$), likely with an oxygen-neon core, for which we estimate an $8.8_{-4.3}^{+1.2}$ M$_{\odot}$ progenitor, close to the mass-divide between white dwarf and neutron star progenitors. A cluster member in Ruprecht 131 is a magnetic white dwarf, whose progenitor mass exceeded 2-3 M$_{\odot}$. We stress that wider searches, and improved cluster distances and ages derived from data of the ESA Gaia mission, will advance the understanding of the mass-loss processes for low- to intermediate-mass stars.

We report the discovery of 42 white dwarfs in the original Kepler mission field, including nine new confirmed pulsating hydrogen-atmosphere white dwarfs (ZZ Ceti stars). Guided by the Kepler-INT Survey (KIS), we selected white dwarf candidates on the basis of their U-g, g-r, and r-H_alpha photometric colours. We followed up these candidates with high-signal-to-noise optical spectroscopy from the 4.2-m William Herschel Telescope. Using ground-based, time-series photometry, we put our sample of new spectroscopically characterized white dwarfs in the context of the empirical ZZ Ceti instability strip. Prior to our search, only two pulsating white dwarfs had been observed by Kepler. Ultimately, four of our new ZZ Cetis were observed from space. These rich datasets are helping initiate a rapid advancement in the asteroseismic investigation of pulsating white dwarfs, which continues with the extended Kepler mission, K2.

We present the first radiation magnetohydrodynamics simulations of the atmosphere of white dwarf stars. We demonstrate that convective energy transfer is seriously impeded by magnetic fields when the plasma-beta parameter, the thermal to magnetic pressure ratio, becomes smaller than unity. The critical field strength that inhibits convection in the photosphere of white dwarfs is in the range B = 1-50 kG, which is much smaller than the typical 1-1000 MG field strengths observed in magnetic white dwarfs, implying that these objects have radiative atmospheres. We have then employed evolutionary models to study the cooling process of high-field magnetic white dwarfs, where convection is entirely suppressed during the full evolution (B > 10 MG). We find that the inhibition of convection has no effect on cooling rates until the effective temperature (Teff) reaches a value of around 5500 K. In this regime, the standard convective sequences start to deviate from the ones without convection owing to the convective coupling between the outer layers and the degenerate reservoir of thermal energy. Since no magnetic white dwarfs are currently known at the low temperatures where this coupling significantly changes the evolution, effects of magnetism on cooling rates are not expected to be observed. This result contrasts with a recent suggestion that magnetic white dwarfs with Teff < 10,000 K cool significantly slower than non-magnetic degenerates.

We present an extended grid of mean three-dimensional (3D) spectra for low-mass, pure-hydrogen atmosphere DA white dwarfs (WDs). We use CO5BOLD radiation-hydrodynamics 3D simulations covering Teff = 6000-11,500 K and logg = 5-6.5 (cgs units) to derive analytical functions to convert spectroscopically determined 1D temperatures and surface gravities to 3D atmospheric parameters. Along with the previously published 3D models, the 1D to 3D corrections are now available for essentially all known convective DA WDs (i.e., logg = 5-9). For low-mass WDs, the correction in temperature is relatively small (a few per cent at the most), but the surface gravities measured from the 3D models are lower by as much as 0.35 dex. We revisit the spectroscopic analysis of the extremely low-mass (ELM) WDs, and demonstrate that the 3D models largely resolve the discrepancies seen in the radius and mass measurements for relatively cool ELM WDs in eclipsing double WD and WD + milli-second pulsar binary systems. We also use the 3D corrections to revise the boundaries of the ZZ Ceti instability strip, including the recently found ELM pulsators.

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

We have expanded the sample of observed white dwarfs in the rich open cluster NGC 2099 (M37) with the Keck Low-Resolution Imaging Spectrometer. Of 20 white dwarf candidates, the spectroscopy shows 19 to be true white dwarfs with 14 of these having high S/N. We find 11 of these 14 to be consistent with singly evolved cluster members. They span a mass range of $\sim$0.7 to 0.95 M$_\odot$, excluding a low-mass outlier, corresponding to progenitor masses of $\sim$3 to 4 M$_\odot$. This region of the initial final mass relation (IFMR) has large scatter and a slope that remains to be precisely determined. With this large sample of white dwarfs that belong to a single age and metallicity population, we find an initial-final mass relation of (0.171$\pm$0.057)M$_{\rm initial}$+0.219$\pm$0.187 M$_\odot$, significantly steeper than the linear relation adopted over the full observed white dwarf mass range in many previous studies. Comparison of this new relation from the solar metallicity NGC 2099 to 18 white dwarfs in the metal-rich Hyades and Praesepe shows that their IFMR also has a consistently steep slope. This strong consistency also suggests that there is no significant metallicity dependence of the IFMR at this mass and metallicity range. As a result, the IFMR can be more reliably determined with this broad sample of 29 total white dwarfs giving M$_{\rm final}$=(0.163$\pm$0.022)M$_{\rm initial}$+0.238$\pm$0.071 M$_\odot$ from M$_{\rm initial}$ of 3 to 4 M$_\odot$. A steep IFMR in this mass range indicates that the full IFMR is nonlinear.

Characterizing the local space density of double degenerate binary systems is a complementary approach to broad sky surveys of double degenerates to determine the expected rates of white dwarf binary mergers, in particular those that may evolve into other observable phenomena such as extreme helium stars, Am CVn systems, and supernovae Ia. However, there have been few such systems detected in local space. We report here the discovery that WD 1242$-$105, a nearby bright WD, is a double-line spectroscopic binary consisting of two degenerate DA white dwarfs of similar mass and temperature, despite it previously having been spectroscopically characterized as a single degenerate. Follow-up photometry, spectroscopy, and trigonometric parallax have been obtained in an effort to determine the fundamental parameters of each component of this system. The binary has a mass ratio of 0.7 and a trigonometric parallax of 25.5 mas, placing it at a distance of 39 pc. The system's total mass is 0.95 M$_\odot$ and has an orbital period of 2.85 hours, making it the strongest known gravitational wave source ($\log h = -20.78$) in the mHz regime. Because of its orbital period and total mass, WD 1242$-$105 is predicted to merge via gravitational radiation on a timescale of 740 Myr, which will most likely not result in a catastrophic explosion.

Using images from the Hubble Space Telescope Wide-Field Camera 3, we measure the rate of diffusion of stars through the core of the globular cluster 47 Tucanae using a sample of young white dwarfs identified in these observations. This is the first direct measurement of diffusion due to gravitational relaxation. We find that the diffusion rate $\kappa\approx 10-13$ arcsecond$^2$ Myr$^{-1}$ is consistent with theoretical estimates of the relaxation time in the core of 47 Tucanae of about 70 Myr.

A calibration of the mixing-length parameter in the local mixing-length theory (MLT) is presented for the lower part of the convection zone in pure-hydrogen atmosphere white dwarfs. The parameterization is performed from a comparison of 3D CO5BOLD simulations with a grid of 1D envelopes with a varying mixing-length parameter. In many instances, the 3D simulations are restricted to the upper part of the convection zone. The hydrodynamical calculations suggest, in those cases, that the entropy of the upflows does not change significantly from the bottom of the convection zone to regions immediately below the photosphere. We rely on this asymptotic entropy value, characteristic of the deep and adiabatically stratified layers, to calibrate 1D envelopes. The calibration encompasses the convective hydrogen-line (DA) white dwarfs in the effective temperature range 6000 < Teff (K) < 15,000 and the surface gravity range 7.0 < log g < 9.0. It is established that the local MLT is unable to reproduce simultaneously the thermodynamical, flux, and dynamical properties of the 3D simulations. We therefore propose three different parameterizations for these quantities. The resulting calibration can be applied to structure and envelope calculations, in particular for pulsation, chemical diffusion, and convective mixing studies. On the other hand, convection has no effect on the white dwarf cooling rates until there is a convective coupling with the degenerate core below Teff ~ 5000 K. In this regime, the 1D structures are insensitive to the MLT parameterization and converge to the mean 3D results, hence remain fully appropriate for age determinations.

The study of the stellar formation history in the solar neighborhood is a powerful technique to recover information about the early stages and evolution of the Milky Way. We present a new method which consists of directly probing the formation history from the nearby stellar remnants. We rely on the volume complete sample of white dwarfs within 20 pc, where accurate cooling ages and masses have been determined. The well characterized initial-final mass relation is employed in order to recover the initial masses (1 < M/Msun < 8) and total ages for the local degenerate sample. We correct for moderate biases that are necessary to transform our results to a global stellar formation rate, which can be compared to similar studies based on the properties of main-sequence stars in the solar neighborhood. Our method provides precise formation rates for all ages except in very recent times, and the results suggest an enhanced formation rate for the solar neighborhood in the last 5 Gyr compared to the range 5 < Age (Gyr) < 10. Furthermore, the observed total age of ~10 Gyr for the oldest white dwarfs in the local sample is consistent with the early seminal studies that have determined the age of the Galactic disk from stellar remnants. The main shortcoming of our study is the small size of the local white dwarf sample. However, the presented technique can be applied to larger samples in the future.

The measurement of precise absolute fluxes for stellar sources has been pursued with increased vigor since the discovery of the dark energy and the realization that its detailed understanding requires accurate spectral energy distributions (SEDs) of redshifted Ia supernovae in the rest frame. The flux distributions of spectrophotometric standard stars were initially derived from the comparison of stars to laboratory sources of known flux but are now mostly based on calculated model atmospheres. For example, pure hydrogen white dwarf (WD) models provide the basis for the HST CALSPEC archive of flux standards. The basic equations for quantitative spectrophotometry and photometry are explained in detail. Several historical lab based flux calibrations are reviewed; and the SEDs of stars in the major on-line astronomical databases are compared to the CALSPEC reference standard spectrophotometry. There is good evidence that relative fluxes from the visible to the near-IR wavelength of ~2.5 micron are currently accurate to 1% for the primary reference standards; and new comparisons with lab flux standards show promise for improving that precision.

We have identified eight to ten new cool white dwarfs from the Large Area Survey (LAS) Data Release 9 of the United Kingdom InfraRed Telescope (UKIRT) Infrared Deep Sky Survey (UKIDSS). The data set was paired with the Sloan Digital Sky Survey (SDSS) to obtain proper motions and a broad ugrizYJHK wavelength coverage. Optical spectroscopic observations were secured at Gemini Observatory and confirm the degenerate status for eight of our targets. The final sample includes two additional white dwarf candidates with no spectroscopic observations. We rely on improved 1D model atmospheres and new multi-dimensional simulations with CO5BOLD to review the stellar parameters of the published LAS white dwarf sample along with our additional discoveries. Most of the new objects possess very cool atmospheres with effective temperatures below 5000 K, including two pure-hydrogen remnants with a cooling age between 8.5 and 9.0 Gyr, and tangential velocities in the range 40 km/s < vtan < 60 km/s. They are likely thick disk 10-11 Gyr-old objects. In addition we find a resolved double degenerate system with vtan ~ 155 km/s and a cooling age between 3.0 and 5.0 Gyr. These white dwarfs could be disk remnants with a very high velocity or former halo G stars. We also compare the LAS sample with earlier studies of very cool degenerates and observe a similar deficit of helium-dominated atmospheres in the range 5000 < Teff (K) < 6000. We review the possible explanations for the spectral evolution from helium-dominated towards hydrogen-rich atmospheres at low temperatures.

Context. White dwarfs can be used to study the structure and evolution of the Galaxy by analysing their luminosity function and initial mass function. Among them, the very cool white dwarfs provide the information for the early ages of each population. Because white dwarfs are intrinsically faint only the nearby (about 20 pc) sample is reasonably complete. The Gaia space mission will drastically increase the sample of known white dwarfs through its 5-6 years survey of the whole sky up to magnitude V = 20-25. Aims. We provide a characterisation of Gaia photometry for white dwarfs to better prepare for the analysis of the scientific output of the mission. Transformations between some of the most common photometric systems and Gaia passbands are derived. We also give estimates of the number of white dwarfs of the different galactic populations that will be observed. Methods. Using synthetic spectral energy distributions and the most recent Gaia transmission curves, we computed colours of three different types of white dwarfs (pure hydrogen, pure helium, and mixed composition with H/He= 0.1). With these colours we derived transformations to other common photometric systems (Johnson-Cousins, Sloan Digital Sky Survey, and 2MASS). We also present numbers of white dwarfs predicted to be observed by Gaia. Results. We provide relationships and colour-colour diagrams among different photometric systems to allow the prediction and/or study of the Gaia white dwarf colours. We also include estimates of the number of sources expected in every galactic population and with a maximum parallax error. Gaia will increase the sample of known white dwarfs tenfold to about 200 000. Gaia will be able to observe thousands of very cool white dwarfs for the first time, which will greatly improve our understanding of these stars and early phases of star formation in our Galaxy.

We present ULTRACAM photometry and X-Shooter spectroscopy of the eclipsing double white dwarf binary CSS 41177, the only such system that is also a double-lined spectroscopic binary. Combined modelling of the light curves and radial velocities yield masses and radii for both white dwarfs without the need to assume mass-radius relations. We find that the primary white dwarf has a mass of M1 = 0.38(2) Msun and a radius of R1 = 0.0222(4) Rsun. The secondary white dwarf's mass and radius are M2 = 0.32(1) Msun and R2 = 0.0207(4) Rsun, and its temperature and surface gravity (T2 = 11678(313) K, log(g2) = 7.32(2)) put it close to the white dwarf instability strip. However, we find no evidence for pulsations to roughly 0.5% relative amplitude. Both masses and radii are consistent with helium white dwarf models with thin hydrogen envelopes of 0.0001 Mstar. The two stars will merge in 1.14 Gyr due to angular momentum loss via gravitational wave emission.

We establish new constraints on the intermediate-mass range of the initial-final mass relation by studying white dwarfs in four young star clusters, and apply the results to study the evolution of stars on the thermally pulsing asymptotic giant branch (TP-AGB). We show that the stellar core mass on the AGB grows rapidly from 10% to 30% for stars with $M_{\rm initial}$ = 1.6 to 2.0 $M_\odot$. At larger masses, the core-mass growth decreases steadily to $\sim$10% at $M_{\rm initial}$ = 3.4 $M_\odot$. These observations are in excellent agreement with predictions from the latest TP-AGB evolutionary models in Marigo et al. (2013). We also compare to models with varying efficiencies of the third dredge-up and mass loss, and demonstrate that the process governing the growth of the core is largely the stellar wind, while the third dredge-up plays a secondary, but non-negligible role. Based on the new white dwarf measurements, we perform an exploratory calibration of the most popular mass-loss prescriptions in the literature. Finally, we estimate the lifetime and the integrated luminosity of stars on the TP-AGB to peak at $t$ $\sim$ 3 Myr and $E$ = 1.2 $\times$ 10$^{10}$ $L_\odot$ yr for $M_{\rm initial}$ $\sim$ 2 $M_\odot$ ($t$ $\sim$ 2 Myr for luminosities brighter than the RGB tip at $\log(L/L_{\odot})$ $>$ 3.4), decreasing to $t$ = 0.4 Myr and $E$ = 6.1 $\times$ 10$^{9}$ $L_\odot$ yr for stars with $M_{\rm initial}$ $\sim$ 3.5 $M_\odot$. The implications of these results are discussed with respect to general population synthesis studies that require correct modeling of the TP-AGB phase of stellar evolution.

We present the first grid of mean three-dimensional (3D) spectra for pure-hydrogen (DA) white dwarfs based on 3D model atmospheres. We use CO5BOLD radiation-hydrodynamics 3D simulations instead of the mixing-length theory for the treatment of convection. The simulations cover the effective temperature range of 6000 < Teff (K) < 15,000 and the surface gravity range of 7 < log g < 9 where the large majority of DAs with a convective atmosphere are located. We rely on horizontally averaged 3D structures (over constant Rosseland optical depth) to compute <3D> spectra. It is demonstrated that our <3D> spectra can be smoothly connected to their 1D counterparts at higher and lower Teff where the 3D effects are small. Analytical functions are provided in order to convert spectroscopically determined 1D effective temperatures and surface gravities to 3D atmospheric parameters. We apply our improved models to well studied spectroscopic data sets from the Sloan Digital Sky Survey and the White Dwarf Catalog. We confirm that the so-called high-log g problem is not present when employing <3D> spectra and that the issue was caused by inaccuracies in the 1D mixing-length approach. The white dwarfs with a radiative and a convective atmosphere have derived mean masses that are the same within ~0.01 Msun, in much better agreement with our understanding of stellar evolution. Furthermore, the 3D atmospheric parameters are in better agreement with independent Teff and log g values from photometric and parallax measurements.

We present a grid of optical (3200--7200 \AA) synthetic spectra calculated with Tlusty/Synspec. The new NLTE model atmospheres include the most recent hydrogen Stark broadening profiles; were calculated in opacity sampling and limited to pure H/He composition. The grid covers the observed parameter space of (He-)sdB and (He-)sdO stars, therefore it is suitable for the homogeneous spectral analyses of such evolved stars.

3D model atmospheres for giants, dwarfs, and white dwarfs, computed with the CO5BOLD code and part of the CIFIST grid, have been used for spectroscopic and asteroseismic studies. Unlike existing plane-parallel 1D structures, these simulations predict the spatially and temporally resolved emergent intensity so that granulation can be analysed, which provides insights on how convective energy transfer operates in stars. The wide range of atmospheric parameters of the CIFIST 3D simulations (3600 < Teff (K) < 13,000 and 1 < log g < 9) allows the comparison of convective processes in significantly different environments. We show that the relative intensity contrast is correlated with both the Mach and Peclet numbers in the photosphere. The horizontal size of granules varies between 3 and 10 times the local pressure scale height, with a tight correlation between the factor and the Mach number of the flow. Given that convective giants, dwarfs, and white dwarfs cover the same range of Mach and Peclet numbers, we conclude that photospheric convection operates in a very similar way in those objects.

A sequence of pure-hydrogen CO5BOLD 3D model atmospheres of DA white dwarfs is presented for a surface gravity of log g = 8 and effective temperatures from 6000 to 13,000 K. We show that convective properties, such as flow velocities, characteristic granulation size and intensity contrast of the granulation patterns, change significantly over this range. We demonstrate that these 3D simulations are not sensitive to numerical parameters unlike the 1D structures that considerably depend on the mixing-length parameters. We conclude that 3D spectra can be used directly in the spectroscopic analyses of DA white dwarfs. We confirm the result of an earlier preliminary study that 3D model spectra provide a much better characterization of the mass distribution of white dwarfs and that shortcomings of the 1D mixing-length theory are responsible for the spurious high-log g determinations of cool white dwarfs. In particular, the 1D theory is unable to account for the cooling effect of the convective overshoot in the upper atmospheres.

We present an empirical determination of the white dwarf cooling sequence in the globular cluster 47 Tucanae. Using spectral models, we determine temperatures for 887 objects from Wide Field Camera 3 data, as well as 292 objects from data taken with the Advanced Camera for Surveys. We make the assumption that the rate of white dwarf formation in the cluster is constant. Stellar evolution models are then used to determine the rate at which objects are leaving the main sequence, which must be the same as the rate at which objects are arriving on the white dwarf sequence in our field. The result is an empirically derived relation between temperature ($T_{eff}$) and time ($t$) on the white dwarf cooling sequence. Comparing this result to theoretical cooling models, we find general agreement with the expected slopes between 20,000K and 30,000K and between 6,000K and 20,000K, but the transition to the Mestel cooling rate of $T_{eff} \propto t^{-0.4}$ is found to occur at hotter temperatures, and more abruptly than is predicted by any of these models.

The Hyades cluster is known to harbour ten so-called classical white dwarf members. Numerous studies through the years have predicted that more than twice this amount of degenerate stars should be associated with the cluster. Using the PPMXL catalog of proper motions and positions, a recent study proposed 17 new white dwarf candidates. We review the membership of these candidates by using published spectroscopic and photometric observations, as well as by simulating the contamination from field white dwarfs. In addition to the ten classical Hyades white dwarfs, we find six white dwarfs that may be of Hyades origin and three more objects that have an uncertain membership status due to their unknown or imprecise atmospheric parameters. Among those, two to three are expected as field stars contamination. Accurate radial velocity measurements will confirm or reject the candidates. One consequence is that the longstanding problem that no white dwarf older than ~340 Myr appears to be associated with the cluster remains unsolved.

We report the identification of LSR J0745+2627 in the United Kingdom InfraRed Telescope Infrared Deep Sky Survey (UKIDSS) Large Area Survey (LAS) as a cool white dwarf with kinematics and age compatible with the thick-disk/halo population. LSR J0745+2627 has a high proper motion (890 mas/yr) and a high reduced proper motion value in the J band (H_J=21.87). We show how the infrared-reduced proper motion diagram is useful for selecting a sample of cool white dwarfs with low contamination. LSR J0745+2627 is also detected in the Sloan Digital Sky Survey (SDSS) and the Wide-field Infrared Survey Explorer (WISE). We have spectroscopically confirmed this object as a cool white dwarf using X-Shooter on the Very Large Telescope. A detailed analysis of its spectral energy distribution reveals that its atmosphere is compatible with a pure-H composition model with an effective temperature of 3880+-90 K. This object is the brightest pure-H ultracool white dwarf (Teff<4000 K) ever identified. We have constrained the distance (24-45 pc), space velocities and age considering different surface gravities. The results obtained suggest that LSR J0745+2627 belongs to the thick-disk/halo population and is also one of the closest ultracool white dwarfs.

We present a new distance determination to the Galactic globular cluster 47 Tucanae by fitting the spectral energy distributions of its white dwarfs to pure hydrogen atmosphere white dwarf models. Our photometric dataset is obtained from a 121 orbit Hubble Space Telescope program using the Wide Field Camera 3 UVIS/IR channels, capturing F390W, F606W, F110W, and F160W images. These images cover more than 60 square arcmins and extend over a radial range of 5-13.7 arcmin (6.5-17.9 pc) within the globular cluster. Using a likelihood analysis, we obtain a best fitting unreddened distance modulus of (m - M)o=13.36+/-0.02+/-0.06 corresponding to a distance of 4.70+/-0.04+/-0.13 kpc, where the first error is random and the second is systematic. We also search the white dwarf photometry for infrared excess in the F160W filter, indicative of debris disks or low mass companions, and find no convincing cases within our sample.

The surface gravities of cool (Teff < 13,000 K) hydrogen-atmosphere DA white dwarfs, determined from spectroscopic analyses, are found to be significantly higher than the canonical value of log g ~ 8 expected for these stars. It was recently concluded that a problem with the treatment of convective energy transport within the framework of the mixing-length theory was the most plausible explanation for this high-log g problem. We pursue the investigation of this discrepancy by computing model spectra of cool convective white dwarfs from a small sequence (11,300 K < Teff < 12,800 K) of 3D hydrodynamical model atmospheres, which feature a sophisticated treatment of convection and radiative transfer. Our approach is to proceed with a differential analysis between 3D and standard 1D models. We find that the 3D spectra predict significantly lower surface gravities, with corrections of the right amplitude as a function of effective temperature to obtain values of log g ~ 8 on average. We conclude that the surface gravity distribution of cool convective DA white dwarfs is much closer to that of hotter radiative objects when using, for the treatment of the convection, 3D models instead of the mixing-length framework.

We present the results of a search for cool white dwarfs in the United Kingdom InfraRed Telescope (UKIRT) Infrared Deep Sky Survey (UKIDSS) Large Area Survey (LAS). The UKIDSS LAS photometry was paired with the Sloan Digital Sky Survey (SDSS) to identify cool hydrogen-rich white dwarf candidates by their neutral optical colors and blue near-infrared colors, as well as faint Reduced Proper Motion magnitudes. Optical spectroscopy was obtained at Gemini Observatory, and showed the majority of the candidates to be newly identified cool degenerates, with a small number of G- to K-type (sub)dwarf contaminants. Our initial search of 280 deg2 of sky resulted in seven new white dwarfs with effective temperature T_eff ~ 6000 K. The current followup of 1400 deg2 of sky has produced thirteen new white dwarfs. Model fits to the photometry show that seven of the newly identified white dwarfs have 4120 K <= T_eff <= 4480 K, and cooling ages between 7.3 Gyr and 8.7 Gyr; they have 40 km/s <= v_tan <= 85 km/s and are likely to be thick disk 10-11 Gyr-old objects. The other half of the sample has 4610 K <= T_eff <= 5260 K, cooling ages between 4.3 Gyr and 6.9 Gyr, and 60 km/s <= v_tan <= 100 km/s. These are either thin disk remnants with unusually high velocities, or lower-mass remnants of thick disk or halo late-F or G stars.

The absolute flux calibration of the James Webb Space Telescope will be based on a set of stars observed by the Hubble and Spitzer Space Telescopes. In order to cross-calibrate the two facilities, several A, G, and white dwarf (WD) stars are observed with both Spitzer and Hubble and are the prototypes for a set of JWST calibration standards. The flux calibration constants for the four Spitzer IRAC bands 1-4 are derived from these stars and are 2.3, 1.9, 2.0, and 0.5% lower than the official cold-mission IRAC calibration of Reach et al. (2005), i.e. in agreement within their estimated errors of ~2%. The causes of these differences lie primarily in the IRAC data reduction and secondarily in the SEDs of our standard stars. The independent IRAC 8 micron band-4 fluxes of Rieke et al. (2008) are about 1.5 +/- 2% higher than those of Reach et al. and are also in agreement with our 8 micron result.

We present an improved spectroscopic and photometric analysis of hydrogen-line DA white dwarfs from the Sloan Digital Sky Survey Data Release 4 based on model atmospheres that include improved Stark broadening profiles with non-ideal gas effects. We also perform a careful visual inspection of all spectroscopic fits with high signal-to-noise ratios (S/N > 12) and present improved atmospheric parameters (Teff and log g) for each white dwarf. Through a comparison of spectroscopic and photometric temperatures, we report the discovery of 35 DA+DB/DC double degenerate candidates and 2 helium-rich DA stars. We also determine that a cutoff at S/N = 15 optimizes the size and quality of the sample for computing the mean mass of DA white dwarfs, for which we report a value of 0.613 Msun. In the following step, we compare our results to previous analyses of the SDSS DR4 and find a good agreement if we account for the shift produced by the improved Stark profiles. Finally, the properties of DA white dwarfs in the SDSS are weighed against those of the Villanova White Dwarf Catalog sample of Gianninas et al. We find systematically lower masses (by about 3% on average), a difference that we trace back to the data reduction procedure of the SDSS. We conclude that a better understanding of these differences will be important to determine the absolute temperature scale and mean mass of DA white dwarfs.

We present optical spectroscopy and near-infrared photometry of 126 cool white dwarfs in the Sloan Digital Sky Survey (SDSS). Our sample includes high proper motion targets selected using the SDSS and USNO-B astrometry and a dozen previously known ultracool white dwarf candidates. Our optical spectroscopic observations demonstrate that a clean selection of large samples of cool white dwarfs in the SDSS (and the SkyMapper, Pan-STARRS, and the Large Synoptic Survey Telescope datasets) is possible using a reduced proper motion diagram and a tangential velocity cut-off (depending on the proper motion accuracy) of 30 km/s. Our near-infrared observations reveal eight new stars with significant absorption. We use the optical and near-infrared photometry to perform a detailed model atmosphere analysis. More than 80% of the stars in our sample are consistent with either pure hydrogen or pure helium atmospheres. However, the eight stars with significant infrared absorption and the majority of the previously known ultracool white dwarf candidates are best explained with mixed hydrogen and helium atmosphere models. The age distribution of our sample is consistent with a Galactic disk age of 8 Gyr. A few ultracool white dwarfs may be as old as 12-13 Gyr, but our models have problems matching the spectral energy distributions of these objects. There are only two halo white dwarf candidates in our sample. However, trigonometric parallax observations are required for accurate mass and age determinations and to confirm their membership in the halo.

We review at length the longstanding problem in the spectroscopic analysis of cool hydrogen-line (DA) white dwarfs (Teff < 13,000 K) where gravities are significantly higher than those found in hotter DA stars. The first solution that has been proposed for this problem is a mild and systematic helium contamination from convective mixing that would mimic the high gravities. We constrain this scenario by determining the helium abundances in six cool DA white dwarfs using high-resolution spectra from the Keck I 10-m telescope. We obtain no detections, with upper limits as low as He/H = 0.04 in some cases. This allows us to put this scenario to rest for good. We also extend our model grid to lower temperatures using improved Stark profiles with non-ideal gas effects from Tremblay & Bergeron and find that the gravity distribution of cool objects remains suspiciously high. Finally, we find that photometric masses are, on average, in agreement with expected values, and that the high-log g problem is so far unique to the spectroscopic approach.

We present new calculations for the Stark broadening of the hydrogen line profiles in the dense atmospheres of white dwarf stars. Our improved model is based on the unified theory of Stark broadening from Vidal, Cooper & Smith, but it also includes non-ideal gas effects from the Hummer & Mihalas occupation probability formalism directly inside the line profile calculations. This approach improves upon previous calculations that relied on the use of an ad-hoc free parameter to describe the dissolution of the line wing opacity in the presence of high electric microfields in the plasma. We present here the first grid of model spectra for hot Teff >~ 12,000 K DA white dwarfs that has no free parameters. The atmospheric parameters obtained from optical and UV spectroscopic observations using these improved models are shown to differ substantially from those published in previous studies.

We present improved calculations for the Stark broadening of hydrogen lines in dense plasmas typical of white dwarf atmospheres. Our new model is based on the unified theory of Stark broadening from Vidal, Cooper, & Smith. For the first time, we account for the non-ideal effects in a consistent way directly inside the line profile calculations. The Hummer & Mihalas theory is used to describe the non-ideal effects due to perturbations on the absorber from protons and electrons. We use a truncation of the electric microfield distribution in the quasi-static proton broadening to take into account the fact that high electric microfields dissociate the upper state of a transition. This approach represents a significant improvement over previous calculations that relied on the use of an ad hoc parameter to mimic these non-ideal effects. We obtain the first model spectra with line profiles that are consistent with the equation of state. We revisit the properties of DA stars in the range 40,000 K > Teff > 13,000 K by analyzing the optical spectra with our improved models. The updated atmospheric parameters are shown to differ substantially from those published in previous studies, with a mean mass shifted by +0.034 Msun. We also show that these revised atmospheric parameters yield absolute visual magnitudes that remain in excellent agreement with trigonometric parallax measurements.

We determine the ratio of helium- to hydrogen-atmosphere white dwarf stars as a function of effective temperature from a model atmosphere analysis of the infrared photometric data from the Two Micron All Sky Survey combined with available visual magnitudes. Our study surpasses any previous analysis of this kind both in terms of the accuracy of the Teff determinations as well as the size of the sample. We observe that the ratio of helium- to hydrogen-atmosphere white dwarfs increases gradually from a constant value of ~0.25 between Teff = 15,000 K and 10,000 K to a value twice as large in the range 10,000 > Teff > 8000 K, suggesting that convective mixing, which occurs when the bottom of the hydrogen convection zone reaches the underlying convective helium envelope, is responsible for this gradual transition. The comparison of our results with an approximate model used to describe the outcome of this convective mixing process implies hydrogen mass layers in the range log M_H/M_tot = -10 to -8 for about 15% of the DA stars that survived the DA to DB transition near Teff ~ 30,000 K, the remainder having presumably more massive layers above log M_H/M_tot ~ -6.

We review the available near- and mid- infrared photometry data sets for white dwarfs from the Two Micron All-Sky Survey (2MASS) Point Source Catalog and the Spitzer Space Telescope. These data sets have been widely used to search for white dwarfs with an infrared excess as well as to characterize the atmosphere of cool white dwarfs. We evaluate the reliability of the 2MASS photometry by performing a statistical comparison with published JHK CIT magnitudes, and by carrying out a detailed model atmosphere analysis of the available photometry. We then present a critical examination of various results published in the literature including data from the Spitzer Space Telescope.

We review the available near- and mid-infrared photometry for white dwarfs obtained from the Two Micron All-Sky Survey (2MASS) and by the Spitzer Space Telescope. Both data sets have recently been used to seek white dwarfs with infrared excesses due to the presence of unresolved companions or circumstellar disks, and also to derive the atmospheric parameters of cool white dwarfs. We first attempt to evaluate the reliability of the 2MASS photometry by comparing it with an independent set of published JHK CIT magnitudes for 160 cool white dwarf stars, and also by comparing the data with the predictions of detailed model atmosphere calculations. The possibility of using 2MASS to identify unresolved M dwarf companions or circumstellar disks is then discussed. We also revisit the analysis of 46 binary candidates from Wachter et al. using the synthetic flux method and confirm the large near-infrared excesses in most objects. We perform a similar analysis by fitting Spitzer 4.5 and 8 micron photometric observations of white dwarfs with our grid of model atmospheres, and demonstrate the reliability of both the Spitzer data and the theoretical calculations up to 8 micron. Finally, we search for massive disks resulting from the merger of two white dwarfs in a 2MASS sample composed of 57 massive degenerates, and show that massive disks are uncommon in such stars.

We present a quantitative and homogeneous analysis of the broadband multicolor photometric data sets gathered so far on rapidly pulsating hot B subdwarf stars. This concerns seven distinct data sets related to six different stars. Our analysis is carried out within the theoretical framework developed by Randall et al., which includes full nonadiabatic effects. The goal of this analysis is partial mode identification, i.e., the determination of the degree index l of each of the observed pulsation modes. We assume possible values of l from 0 to 5 in our calculations. For each target star, we compute a specific model atmosphere and a specific pulsation model using estimates of the atmospheric parameters coming from time-averaged optical spectroscopy. For every assumed value of l, we use a formal chi-squared approach to model the observed amplitude-wavelength distribution of each mode, and we compute a quality-of-fit Q probability to quantify the derived fit and to discriminate objectively between the various solutions. We find that no completely convincing and unambiguous l identification is possible on the basis of the available data, although partial mode discrimination has been reached for 25 out of the 41 modes studied. A brief statistical study of these results suggests that a majority of the modes must have l values of 0, 1, and 2, but also that modes with l = 4 could very well be present while modes with l = 3 appear to be rarer. This is in line with recent results showing that l = 4 modes in rapidly pulsating B subdwarfs have a higher visibility in the optical domain than modes with l = 3. Although somewhat disappointing in terms of mode discrimination, our results still suggest that the full potential of multicolor photometry for l identification in pulsating subdwarfs is within reach.