Thanks to their continuous cooling and relative simplicity, white dwarf stars are routinely used to measure the ages of stellar populations. The usefulness of white dwarfs as cosmochronometers depends on the availability of accurate cooling models. A key ingredient of those models are the conductive opacities, which largely govern the cooling rate. In this work, we present improved conductive opacities for the regime of moderate coupling and moderate degeneracy that characterizes an important portion of the envelopes of DA and DB white dwarfs. We find differences of up to a factor 3 between our calculations and the commonly used opacities of Cassisi et al. (2007), which we attribute to an improved account of electron-electron scattering. The cooling models are strongly affected by those changes in the conductive opacities: the age of a 4000 K white dwarf can be reduced by as much as 2 Gyr. We provide analytical fits to our new opacities to facilitate the implementation of this important effect in white dwarf evolution codes.

We report distinctive core profiles in the strongest optical helium line, He I 5876, from high-resolution high-sensitivity observations of spectral type DB white dwarfs. By analyzing a sample of 40 stars from Keck/HIRES and VLT/UVES, we find the core appearance to be related to the degree of hydrogen and heavy element content in the atmosphere. New Ca K-line measurements or upper limits are reported for about half the sample stars. He I 5876 emission cores with a self-reversed central component are present for those stars with relatively low hydrogen abundance, as well as relatively low atmospheric heavy element pollution. This self-reversed structure disappears for stars with higher degrees of pollution and/or hydrogen abundance, giving way to a single absorption core. From our model atmospheres, we show that the self-reversed emission cores can be explained by temperature inversions in the upper atmosphere. We propose that the transition to a single absorption core is due to the additional opacity from hydrogen and heavy elements that inhibits the temperature inversions. Our current models do not exactly match the effective temperature range of the phenomenon nor the amplitude of the self-reversed structure, which is possibly a result of missing physics such as 3D treatment, convective overshoot, and/or non-LTE effects. The He I 5876 line structure may prove to be a useful new diagnostic for calibrating temperature profiles in DB atmosphere models.

The accretion of rocky material is responsible for the presence of heavy elements in the atmospheres of a large fraction of white dwarf stars. Those objects represent a unique opportunity to infer the bulk composition of exoplanetesimals. This chemical characterization requires the use of detailed atmosphere models to determine the elemental abundances at the photospheres of white dwarfs. In this work, we use a state-of-the-art model atmosphere code to reanalyse the first large survey of metal-polluted white dwarfs for which abundances are found for multiple elements. We show that the improved constitutive physics of our models lead to systematically higher Mg abundances than previous analyses. We find an average log Mg/Ca number abundance ratio of 1.5. This value is significantly above the reference abundance for chondrites, which is expected as current diffusion models predict that for the cool helium-atmosphere white dwarfs of our sample Mg should remain in the atmosphere longer than Ca. This helps resolve a recently identified Mg depletion problem, where the planetesimals accreted by white dwarfs were reported to be Mg-deficient compared to the expected composition of their planetary systems.

We present follow-up spectroscopy of 711 white dwarfs within 100 pc, and present a detailed model atmosphere analysis of the 100 pc white dwarf sample in the SDSS footprint. Our spectroscopic follow-up is complete for 83% of the white dwarfs hotter than 6000 K, where the atmospheric composition can be constrained reliably. We identify 1508 DA white dwarfs with pure hydrogen atmospheres. The DA mass distribution has an extremely narrow peak at $0.59~M_{\odot}$, and reveals a shoulder from relatively massive white dwarfs with $M=0.7$-$0.9~M_{\odot}$. Comparing this distribution with binary population synthesis models, we find that the contribution from single stars that form through mergers cannot explain the over-abundance of massive white dwarfs. In addition, the mass distribution of cool DAs shows a near absence of $M>1~M_{\odot}$ white dwarfs. The pile-up of 0.7-$0.9~M_{\odot}$ and the disappearance of $M>1~M_{\odot}$ white dwarfs is consistent with the effects of core crystallization. Even though the evolutionary models predict the location of the pile-up correctly, the delay from the latent heat of crystallization by itself is insufficient to create a significant pile-up, and additional cooling delays from related effects like phase separation are necessary. We also discuss the population of infrared-faint (ultracool) white dwarfs, and demonstrate for the first time the existence of a well defined sequence in color and magnitude. Curiously, this sequence is connected to a region in the color-magnitude diagrams where the number of helium-dominated atmosphere white dwarfs is low. This suggests that the infrared-faint white dwarfs likely have mixed H/He atmospheres.

Ultra-cool white dwarfs are among the oldest stellar remnants in the Universe. Their efficient gravitational settling and low effective temperatures are responsible for the smooth spectra they exhibit. For that reason, it is not possible to derive their radial velocities or to find the chemistry of the progenitors. The best that can be done is to infer such properties from associated sources, which are coeval. The simplest form of such a system is a common proper motion pair where one star is an evolved degenerate and the other a main-sequence star. In this work, we present the discovery of the first of such a system, the M dwarf LHS 6328 and the ultra-cool white dwarf PSO J1801+625, from the Pan-STARRS 1 3$\pi$ survey and the $Gaia$ Data Release 2. Follow-up spectra were collected covering a usable wavelength range of 3750-24500 {\AA}. Their spectra show that the white dwarf has an effective temperature of 3550K and surface gravity of $\log g = 7.45 \pm 0.13$ or $\log g = 7.49 \pm 0.13$ for a CO or He core, respectively, when compared against synthetic spectra of ultra-cool white dwarf atmosphere models. The system has slightly subsolar metallicity with -0.25 < [Fe/H] < 0.0, and a spatial velocity of (U, V, W) = ($-114.26 \pm 0.24$, $222.94 \pm 0.60$, $10.25 \pm 0.34$) km s $^{-1}$ the first radial velocity and metallicity measurements of an ultra-cool white dwarf. This makes it the first and only benchmark of its kind to date.

Taking advantage of the Gaia Data Release 2, recent studies have revisited the evolution of carbon-polluted white dwarfs (DQs) across a large range of effective temperatures. These analyses have clearly confirmed the existence of two distinct DQ evolutionary sequences: one with normal-mass white dwarfs and one with heavily polluted and generally more massive objects. The first sequence is thought to result from the dredge-up of carbon from the core, while the second could at least partially be made of descendants of Hot DQs. However, the evolution of carbon-polluted white dwarfs below 6500 K remains unexplored, mainly due to the theoretical difficulties associated with modelling their dense atmospheres. In this work, we present a detailed star-by-star analysis of cool carbon-polluted white dwarfs. Our recently improved atmosphere models allow us to obtain good fits to most objects, including very cool DQpec white dwarfs with strongly shifted C$_2$ molecular bands. We show that cool carbon-polluted white dwarfs keep following the two distinct evolutionary tracks previously identified at higher temperatures. We also find that most DQ white dwarfs transform into DQpec when their photospheric densities exceed $\approx$ 0.15 g/cm$^3$. However, we identify stars for which the DQ$\rightarrow$DQpec transition occurs at lower photospheric densities, possibly due to the presence of a strong magnetic field.

We present a homogeneous analysis of 1023 DBZ/DZ(A) and 319 DQ white dwarf stars taken from the Montreal White Dwarf Database. This represents a significant increase over the previous comprehensive studies on these types of objects. We use new trigonometric parallax measurements from the Gaia second data release, together with photometry from the Sloan Digital Sky Survey, Pan-STARRS, Gaia, or BVRI from the literature, which allow the determination of the mass for the majority of the objects in our sample. We use the photometric and spectroscopic techniques with the most recent model atmospheres available, which include high-density effects, to accurately determine the effective temperature, surface gravity, and heavy element abundances for each object. We study the abundance of hydrogen in DBZ/DZ white dwarfs and the properties of the accreted planetesimals. We explore the nature of the second sequence of DQ stars using proper motions from Gaia, and highlight evidence of crystallization in massive DQ stars. We also present mass distributions for both spectral types. Finally, we discuss the implications of our findings in the context of the spectral evolution of white dwarfs, and provide the atmospheric parameters for each star.

As a result of competing physical mechanisms, the atmospheric composition of white dwarfs changes throughout their evolution, a process known as spectral evolution. Because of the ambiguity of their atmospheric compositions and the difficulties inherent to the modeling of their dense atmospheres, no consensus exists regarding the spectral evolution of cool white dwarfs ($T_{\rm eff}<6000$ K). In the previous papers of this series, we presented and observationally validated a new generation of cool white dwarf atmosphere models that include all the necessary constitutive physics to accurately model those objects. Using these new models and a homogeneous sample of 501 cool white dwarfs, we revisit the spectral evolution of cool white dwarfs. Our sample includes all spectroscopically identified white dwarfs cooler than 8300 K for which a parallax is available in Gaia DR2 and photometric observations are available in Pan-STARRS1 and 2MASS. Except for a few cool carbon-polluted objects, our models allow an excellent fit to the spectroscopic and photometric observations of all objects included in our sample. We identify a decrease of the ratio of hydrogen to helium-rich objects between 7500 K and 6250 K, which we interpret as the signature of convective mixing. After this decrease, hydrogen-rich objects become more abundant up to 5000 K. This puzzling increase, reminiscent of the non-DA gap, has yet to be explained. At lower temperatures, below 5000 K, hydrogen-rich white dwarfs become rarer, which rules out the scenario according to which accretion of hydrogen from the interstellar medium dominates the spectral evolution of cool white dwarfs.

We present a critical review of the determination of fundamental parameters of white dwarfs discovered by the Gaia mission. We first reinterpret color-magnitude and color-color diagrams using photometric and spectroscopic information contained in the Montreal White Dwarf Database (MWDD), combined with synthetic magnitudes calculated from a self-consistent set of model atmospheres with various atmospheric compositions. The same models are then applied to measure the fundamental parameters of white dwarfs using the so-called photometric technique, which relies on the exquisite Gaia trigonometric parallax measurements, and photometric data from Pan-STARRS, SDSS, and Gaia. In particular, we discuss at length the systematic effects induced by these various photometric systems. We then study in great detail the mass distribution as a function of effective temperature for the white dwarfs spectroscopically identified in the MWDD, as well as for the white dwarf candidates discovered by Gaia. We pay particular attention to the assumed atmospheric chemical composition of cool, non-DA stars. We also briefly revisit the validity of the mass-radius relation for white dwarfs, and the recent discovery of the signature of crystallization in the Gaia color-magnitude diagram for DA white dwarfs. We finally present evidence that the core composition of most of these white dwarfs is, in bulk, a mixture of carbon and oxygen, an expected result from stellar evolution theory, but never empirically well established before.

Metal-polluted white dwarfs (DZ stars) are characterized by a helium-rich atmosphere contaminated by heavy elements traces originating from accreted rocky planetesimals. As a detailed spectroscopic analysis of those objects can reveal the composition of the accreted debris, there is a great interest in developing accurate DZ atmosphere models. However, the coolest DZ white dwarfs are challenging to model due to the fluidlike density of their atmospheres. Under such extreme conditions, spectral absorption lines are heavily broadened by interactions with neutral helium and it is no longer justified to use the conventional Lorentzian profiles. In this work, we determine the theoretical profiles of the Ca I resonance line (the most prominent spectral line for the coolest DZ white dwarfs) in the dense atmospheres of cool DZ white dwarfs. To do so, we use a unified theory of collisional line profiles and accurate ab initio potential energies and transition dipole moments for the CaHe molecule. We present the resulting profiles for the full range of temperatures and helium densities relevant for the modeling of cool, metal-polluted white dwarfs (from 3000 to 6000 K and from $10^{21}$ to $10^{23}\,{\rm cm}^{-3}$). We also implement these new profiles in our atmosphere models and show that they lead to improved fits to the Ca I resonance line of the coolest DZ white dwarfs.

WD J2356$-$209 is a cool metal-polluted white dwarf whose visible spectrum is dominated by a strong and broad sodium feature. Although discovered nearly two decades ago, no detailed and realistic analysis of this star had yet been realized. In the absence of atmosphere models taking into account the nonideal high-density effects arising at the photosphere of WD J2356$-$209, the origin of its unique spectrum had remained nebulous. We use the cool white dwarf atmosphere code presented in the first paper of this series to finally reveal the secrets of this peculiar object and details about the planetesimal that polluted its atmosphere. Thanks to the improved input physics of our models, we find a solution that is in excellent agreement with the photometric observations and the visible spectrum. Our solution reveals that the photosphere of WD J2356$-$209 has a number density ratio of $\log\,{\rm Na/Ca}= 1.0 \pm 0.2$, which is the highest ever found in a white dwarf. Since we do not know how long ago the accretion episode stopped (if it has), we cannot precisely determine the composition nor the mass of the accreted planetesimal. Nevertheless, all scenarios considered indicate that its composition is incompatible with that of chondrite-like material and that its mass was at least $10^{21}\,{\rm g}$.

In the first paper of this series (Blouin et al. 2018), we presented our upgraded cool white dwarf atmosphere code. In this second paper, we use our new models to analyze SDSS J080440.63+223948.6 (J0804+2239), the first DZ star to show collision-induced absorption (CIA). This object provides a crucial test for our models, since previous versions of our code were unable to simultaneously fit the metal absorption lines and the CIA. We find an excellent fit to both the spectroscopic and photometric data, which further validates the improved constitutive physics of our models. We also show that the presence of metal lines allows to lift the degeneracy between high and low hydrogen abundances that usually affects the fits of white dwarfs with CIA. Finally, we investigate the potential impact of spectroscopically undetected metals on the photometric solutions of DC stars.

Line shapes of the magnesium resonance lines in white dwarf spectra are determined by the properties of magnesium atoms and the structure of the white dwarf atmosphere. Through their blanketing effect, these lines have a dominant influence on the model structure and thus on the determination from the spectra of other physical parameters that describe the stellar atmosphere and elemental abundances. In continuation of previous work on Mg$^+$He lines in the UV, we present theoretical profiles of the resonance line of neutral Mg perturbed by He at the extreme density conditions found in the cool largely transparent atmosphere of DZ white dwarfs. We accurately determined the broadening of Mg by He in a unified theory of collisional line profiles using ab initio calculations of MgHe potential energies and transition matrix elements among the singlet electronic states that are involved for the observable spectral lines. We computed the shapes and line parameters of the Mg lines and studied their dependence on helium densities and temperatures. We present results over the full range of temperatures from 4000 to 12000 K needed for input to stellar spectra models. Atmosphere models were constructed for a range of effective temperatures and surface gravities typical for cool DZ white dwarfs. We present synthetic spectra tracing the behavior of the Mg resonance line profiles under the low temperatures and high gas pressures prevalent in these atmospheres. The determination of accurate opacity data of magnesium resonance lines together with an improved atmosphere model code lead to a good fit of cool DZ white dwarf stars. The broadening of spectral lines by helium needs to be understood to accurately determine the H/He and Mg/He abundance ratio in DZ white dwarf atmospheres. We emphasize that no free potential parameters or ad hoc adjustments were used to calculate the line profiles.

The photospheres of the coolest helium-atmosphere white dwarfs are characterized by fluid-like densities. Under those conditions, standard approximations used in model atmosphere codes are no longer appropriate. Unfortunately, the majority of cool He-rich white dwarfs show no spectral features, giving us no opportunities to put more elaborate models to the test. In the few cases where spectral features are observed (such as in cool DQ or DZ stars), current models completely fail to reproduce the spectroscopic data, signaling shortcomings in our theoretical framework. In order to fully trust parameters derived solely from the energy distribution, it is thus important to at least succeed in reproducing the spectra of the few coolest stars exhibiting spectral features, especially since such stars possess even less extreme physical conditions due to the presence of heavy elements. In this paper, we revise every building block of our model atmosphere code in order to eliminate low-density approximations. Our updated white dwarf atmosphere code incorporates state-of-the-art constitutive physics suitable for the conditions found in cool helium-rich stars (DC and DZ white dwarfs). This includes new high-density metal line profiles, nonideal continuum opacities, an accurate equation of state and a detailed description of the ionization equilibrium. In particular, we present new ab initio calculations to assess the ionization equilibrium of heavy elements (C, Ca, Fe, Mg and Na) in a dense helium medium and show how our improved models allow us to achieve better spectral fits for two cool DZ stars, Ross 640 and LP 658-2.

Collision-induced absorption (CIA) from molecular hydrogen is a dominant opacity source in the atmosphere of cool white dwarfs. It results in a significant flux depletion in the near-IR and IR parts of their spectra. Because of the extreme conditions of helium-rich atmospheres (where the density can be as high as a few g/cm$^3$), this opacity source is expected to undergo strong pressure distortion and the currently used opacities have not been validated at such extreme conditions. To check the distortion of the CIA opacity we applied state-of-the-art ab initio methods of computational quantum chemistry to simulate the CIA opacity at high densities. The results show that the CIA profiles are significantly distorted above densities of $0.1\,{\rm g/cm}^3$ in a way that is not captured by the existing models. The roto-translational band is enhanced and shifted to higher frequencies as an effect of the decrease of the interatomic separation of the H$_2$ molecule. The vibrational band is blueward shifted and split into $Q_R$ and $Q_P$ branches, separated by a pronounced interference dip. Its intensity is also substantially reduced. The distortions result in a shift of the maximum of the absorption from $2.3\,\mu{\rm m}$ to $3-7 \mu{\rm m}$, which could potentially explain the spectra of some very cool, helium-rich white dwarfs.

Dense, He-rich atmospheres of cool white dwarfs represent a challenge to the modeling. This is because these atmospheres are constituted of a dense fluid in which strong multi-atomic interactions determine their physics and chemistry. Therefore, the ideal-gas-based description of absorption is no longer adequate, which makes the opacities of these atmospheres difficult to model. This is illustrated with severe problems in fitting the spectra of cool, He-rich stars. Good description of the infrared (IR) opacity is essential for proper assignment of the atmospheric parameters of these stars. Using methods of computational quantum chemistry we simulate the IR absorption of dense He/H media. We found a significant IR absorption from He atoms (He-He-He CIA opacity) and a strong pressure distortion of the H$_2$-He collision-induced absorption (CIA). We discuss the implication of these results for interpretation of the spectra of cool stars.

Due to their high photospheric density, cool helium-rich white dwarfs (particularly DZ, DQpec and ultracool) are often poorly described by current atmosphere models. As part of our ongoing efforts to design atmosphere models suitable for all cool white dwarfs, we investigate how the ionization ratio of heavy elements and the H$_2$-He collision-induced absorption (CIA) spectrum are altered under fluid-like densities. For the conditions encountered at the photosphere of cool helium-rich white dwarfs, our ab initio calculations show that the ionization of most metals is inhibited and that the H$_2$-He CIA spectrum is significantly distorted for densities higher than 0.1 g/cm$^3$.

We present the "Montreal White Dwarf Database" (MWDD), an accessible database with sortable/filterable table and interactive plots that will, when fully completed, allow the community to explore the physical properties of all white dwarfs ever analyzed by the Montreal group, as well as display data and analyses from the literature. We present its current capability and show how it will continuously be updated to instantly reflect improvements made on both the theoretical and observational fronts.