We present a new set of solar metallicity atmosphere and evolutionary models for very cool brown dwarfs and self-luminous giant exoplanets, which we term ATMO 2020. Atmosphere models are generated with our state-of-the-art 1D radiative-convective equilibrium code ATMO, and are used as surface boundary conditions to calculate the interior structure and evolution of $0.001-0.075\,\mathrm{M_{\odot}}$ objects. Our models include several key improvements to the input physics used in previous models available in the literature. Most notably, the use of a new H-He equation of state including ab initio quantum molecular dynamics calculations has raised the mass by $\sim1-2\%$ at the stellar-substellar boundary and has altered the cooling tracks around the hydrogen and deuterium burning minimum masses. A second key improvement concerns updated molecular opacities in our atmosphere model ATMO, which now contains significantly more line transitions required to accurately capture the opacity in these hot atmospheres. This leads to warmer atmospheric temperature structures, further changing the cooling curves and predicted emission spectra of substellar objects. We present significant improvement for the treatment of the collisionally broadened potassium resonance doublet, and highlight the importance of these lines in shaping the red-optical and near-infrared spectrum of brown dwarfs. We generate three different grids of model simulations, one using equilibrium chemistry and two using non-equilibrium chemistry due to vertical mixing, all three computed self-consistently with the pressure-temperature structure of the atmosphere. We show the impact of vertical mixing on emission spectra and in colour-magnitude diagrams, highlighting how the $3.5-5.5\,\mathrm{\mu m}$ flux window can be used to calibrate vertical mixing in cool T-Y spectral type objects.

The spectra of helium-dominated white dwarf stars with hydrogen in their atmosphere present a distinctive broad feature centered around 1160~\AA\/ in the blue wing of the Lyman-$\alpha$ line. It is extremely apparent in WD 1425+540 recently observed with HST COS. With new theoretical line profiles based on ab initio atomic interaction potentials we show that this feature is a signature of a collision-induced satellite due to an asymptotically forbidden transition. This quasi-molecular spectral satellite is crucial to understanding the asymmetrical shape of Lyman-$\alpha$ seen in this and other white dwarf spectra. Our previous work predicting this absorption feature was limited by molecular potentials that were not adequate to follow the atomic interactions with spectroscopic precision to the asymptotic limit of large separation. A new set of potential energy curves and electronic dipole transition moments for the lowest electronic states of the H--He system were developed to account accurately for the behavior of the atomic interactions at all distances, from the chemical regime within 1~\AA\/ out to where the radiating H atoms are not significantly perturbed by their neighbors. We use a general unified theory of collision-broadened atomic spectral lines to describe a rigorous treatment of hydrogen Lyman-$\alpha$ with these potentials and present a new study of its broadening by radiative collisions of hydrogen and neutral helium. These results enable ab initio modeling of radiative transport in DBA white dwarf atmospheres.

Atmospheric compositions can provide powerful diagnostics of formation and migration histories of planetary systems. We investigate constraints on atmospheric abundances of H$_2$O, Na, and K, in a sample of transiting exoplanets using latest transmission spectra and new H$_2$ broadened opacities of Na and K. Our sample of 19 exoplanets spans from cool mini-Neptunes to hot Jupiters, with equilibrium temperatures between $\sim$300 and 2700 K. Using homogeneous Bayesian retrievals we report atmospheric abundances of Na, K, and H$_2$O, and their detection significances, confirming 6 planets with strong Na detections, 6 with K, and 14 with H$_2$O. We find a mass-metallicity trend of increasing H$_2$O abundances with decreasing mass, spanning generally substellar values for gas giants and stellar/superstellar for Neptunes and mini-Neptunes. However, the overall trend in H$_2$O abundances, from mini-Neptunes to hot Jupiters, is significantly lower than the mass-metallicity relation for carbon in the solar system giant planets and similar predictions for exoplanets. On the other hand, the Na and K abundances for the gas giants are stellar or superstellar, consistent with each other, and generally consistent with the solar system metallicity trend. The H$_2$O abundances in hot gas giants are likely due to low oxygen abundances relative to other elements rather than low overall metallicities, and provide new constraints on their formation mechanisms. The differing trends in the abundances of species argue against the use of chemical equilibrium models with metallicity as one free parameter in atmospheric retrievals, as different elements can be differently enhanced.

The opacity of alkali atoms, most importantly of Na and K, plays a crucial role in the atmospheres of brown dwarfs and exoplanets. We present a comprehensive study of NaH2 collisional profiles at temperatures from 500 to 3000 K, the temperatures prevailing in the atmosphere of brown dwarfs and Jupiter-mass planets.The relevant H2 perturber densities reach several 10^19 cm^-3 in hot Teff > 1500 K Jupiter-mass planets and can exceed 10^20 cm^-3 for more massive or cooler objects. Accurate pressure-broadened profiles that are valid at high densities of H2 should be incorporated into spectral models. Unified profiles of sodium perturbed by molecular hydrogen were calculated in the semi-classical approach using up-to-date molecular data.New NaH2 collisional profiles and their effects on the synthetic spectra of brown dwarfs and hot Jupiters computed with petitCODE are presented.

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.

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.

The hot Jupiter HD 209458b was observed during primary transit at 3.6, 4.5, 5.8 and 8.0 microns using the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. We detail here the procedures we adopted to correct for the systematic trends present in the IRAC data. The light curves were fitted including limb darkening effects and fitted using Markov Chain Monte Carlo and prayer-bead Monte Carlo techniques, finding almost identical results. The final depth measurements obtained by a combined Markov Chain Monte Carlo fit are at 3.6 microns, 1.469 +- 0.013 % and 1.448 +- 0.013 %; at 4.5 microns, 1.478 +- 0.017 % ; at 5.8 microns, 1.549 +- 0.015 % and at 8.0 microns 1.535 +- 0.011 %. Our results clearly indicate the presence of water in the planetary atmosphere. Our broad band photometric measurements with IRAC prevent us from determining the additional presence of other other molecules such as CO, CO2 and methane for which spectroscopy is needed. While water vapour with a mixing ratio of 10^-4-10^-3 combined with thermal profiles retrieved from the day-side may provide a very good fit to our observations, this data set alone is unable to resolve completely the degeneracy between water abundance and atmospheric thermal profile.

T-type dwarfs present a broad and shallow absorption feature centred around 6950 A in the blue wing of the K doublet at 0.77 micron which resembles in depth and shape the satellite absorption predicted by detailed collisional broadening profiles. In our previous work, the predicted line satellite position was however somewhat too blue compared to the observed feature. In this paper we investigate whether new calculations of the energy surfaces of the potentials in the K-H_2 system, including spin-orbit coupling, result in a closer coincidence of the satellite with the observed position. We also investigate the extent to which CaH absorption bands contribute to the feature. We present model atmospheres and synthetic spectra, including gravitational settling for an improved description of depth-dependent abundances of refractory elements, and based on new K-H_2 line profiles using improved interaction potentials. By comparison with a high signal-to-noise optical spectrum of the T1 dwarf epsilon Indi Ba, we find that these new models do reproduce the observed feature, while CaH does not contribute for the atmospheric parameters considered. We also find that CaH is settled out so deep into the atmosphere that even turbulent vertical mixing would appear insufficient to bring significant amounts of CaH to the photosphere in dwarfs later than ~L5. We conclude that previous identification of the feature at this location in T and late L dwarf spectra with CaH was erroneous, as expected on physical grounds: calcium condenses onto grains in early L dwarfs and thus should have settled out of the photosphere in cooler brown dwarfs. This finding revokes one observational verification for the cloud-clearing theory: a gradual clearing of the cloud cover in early T dwarfs.

We present PHOENIX atmosphere models for metal-rich cool white dwarfs using improved line shapes for the Na I and K I resonance doublets. Profiles for collisional broadening due to H2 and He based on the adiabatic representation show strong deviations from Van der Waals interaction at short distances. Comparison with observed spectra that show extremely broadened Na I lines indicates that a He-rich atmospheric composition is required to explain the line strengths and spectral energy distributions. Our current synthetic spectra, using an expansion in powers of density to the third order optimised for brown dwarf atmosphere conditions, significantly underestimate the observed absorption in the far wings, even predicting smaller total line strength than a Lorentzian profile. This is due to the handling of multiple perturber interactions becoming inadequate for the extreme densities of the coolest white dwarfs. The density expansion would have to be extended at least to the 7th order for an accurate treatment of such conditions and might break down altogether in the densest objects. The results of a direct calculation of the unified profile should therefore be used for model atmospheres of cool metal-rich white dwarfs. Qualitative comparison of the full adiabatic profile to the spectrum of WD2356-209 indicates good agreement with the observed line shape. Observations of the coolest white dwarfs may therefore serve as a laboratory for testing the physics of the deeper atmospheres and interiors of brown dwarfs and giant planets.

Ultracool stellar atmospheres show absorption by alkali resonance lines severely broadened by collisions with neutral perturbers. In the coolest and densest atmospheres, such as those of T dwarfs, Na I and K I broadened by molecular hydrogen and helium can come to dominate the entire optical spectrum. Their profiles have been successfully modelled with accurate interaction potentials in the adiabatic theory, computing line profiles from the first few orders of a density expansion of the autocorrelation function. The line shapes in the emergent spectrum also depend on the distribution of absorbers as a function of depth, which can be modelled with improved accuracy by new models of dust condensation and settling. The far red K I wings of the latest T dwarfs still show missing opacity in these models, a phenomenon similar to what has been found for the Na I line profiles observed in extremely cool, metal-rich white dwarfs. We show that the line profile in both cases is strongly determined by multiple-perturber interactions at short distances and can no longer be reproduced by a density expansion, but requires calculation of the full profile in a unified theory. Including such line profiles in stellar atmosphere codes will further improve models for the coolest and densest dwarfs as well as for the deeper atmosphere layers of substellar objects in general.

Water is predicted to be among, if not the most abundant molecular species after hydrogen in the atmospheres of close-in extrasolar giant planets (hot-Jupiters) Several attempts have been made to detect water on an exoplanet, but have failed to find compelling evidence for it or led to claims that should be taken with caution. Here we report an analysis of recent observations of the hot-Jupiter HD189733b taken during the transit, where the planet passed in front of its parent star. We find that absorption by water vapour is the most likely cause of the wavelength-dependent variations in the effective radius of the planet at the infrared wavelengths 3.6, 5.8 and 8 microns. The larger effective radius observed at visible wavelengths may be due to either star variability or the presence of clouds/hazes. We explain the most recent thermal infrared observations of the planet during secondary transit behind the star, reporting a non-detection of water on HD189733b, as being a consequence of the nearly isothermal vertical profile of the planet.s atmosphere. Our results show that water is detectable on extrasolar planets using the primary transit technique and that the infrared should be a better wavelength region than the visible, for such searches.

Two peculiar stars showing an apparent extremely broadened and strong NaI D absorption have been discovered in surveys for cool white dwarfs by Oppenheimer et al. (2001) and Harris et al. (SDSS, 2003). We discuss the nature of these objects using PHOENIX atmosphere models for metal-poor brown dwarfs/very low mass stars, and new white dwarf LTE and NLTE models for hydrogen- and helium-dominated atmospheres with metals. These include complete molecular formation in chemical equilibrium and a model for the alkali resonance line broadening based on the damping profiles of Allard et al. (2003), as well as new molecular line opacities for metal hydrides. First results of our calculations indicate good agreement with a hydrogen-dominated WD atmosphere with a Na abundance roughly consistent with a state of high accretion. We analyse deviations of the abundances of Na, K, Mg and Ca from the cosmic pattern and comment on implications of these results for standard accretion scenarios.

We present new FUV observations of the pulsating DA white dwarf G226-29 obtained with the Far Ultraviolet Spectroscopic Explorer (FUSE). This ZZ Ceti star is the brightest one of its class and the coolest white dwarf observed by FUSE. We report the first detection of the broad quasi-molecular collision-induced satellite of Ly-beta at 1150 A, an absorption feature that is due to transitions which take place during close collisions of hydrogen atoms. The physical interpretation of this feature is based on recent progress of the line broadening theory of the far wing of Ly-beta. This predicted feature had never been observed before, even in laboratory spectra.

We present an analysis of the Far Ultraviolet Spectroscopic Explorer (FUSE) spectrum of the ultramassive (M = 1.31 solar mass), magnetic (B_s = 2.3 MG) white dwarf PG 1658+441. The far ultraviolet (FUV) spectrum exhibits very broad Lyman lines and quasi-molecular Lyman beta satellites, but weak Lyman gamma satellites may also be present. PG 1658+441 is the hottest white dwarf known to show these satellite features. We fit the Lyman lines with stellar models and obtain atmospheric parameters consistent with a published analysis of the Balmer lines. By averaging results obtained for the different FUSE segments, we determine Teff = 29,620 +/- 500K and log g = 9.31 +/- 0.07. The models match the data over large portions of the spectrum but discrepancies remain near the satellite features. Finally, no trace elements have been identified in the FUV spectrum, and we provide abundance upper limits for C, N, Si, and P.

We present new theoretical calculations of the line profile of Lyman gamma that include transitions in which a photon is absorbed by a neutral hydrogen atom while it interacts with a proton. Models show that two absorption features located near 992 A and 996 A are due to H-H^+ collisions. These quasi-molecular satellites are similar to those that were identified in the wings of Lyman alpha and Lyman beta lines of hydrogen-rich white dwarfs. We compute synthetic spectra that take account of these new theoretical profiles and compare them to the spectra of four DA white dwarfs that were observed with FUSE. The models predict the absorption features that are observed in the wing of Lyman gamma near 995 A, and confirm that these features are quasi-molecular satellites.

We present new theoretical calculations of the total line profiles of Lyman alpha and Lyman beta which include perturbations by both neutral hydrogen AND protons and all possible quasi-molecular states of H_2 and H_2^+. They are used to improve theoretical modeling of synthetic spectra for cool DA white dwarfs. We compare them with FUSE observation of G231-40. The appearance of the line wings between Lyman alpha and Lyman beta is shown to be sensitive to the relative abundance of hydrogen ions and neutral atoms, and thereby to provide a temperature diagnostic for stellar atmospheres and laboratory plasmas.

We have analyzed FUSE observations of six hot white dwarf stars: four DA white dwarfs with T_eff >= 45000K, the DAO Feige55 (T_eff = 55000K), and the DA CD -38 10980 (T_eff = 24000K). Photospheric lines from SiIV, PV, and SVI can be observed in the majority of the five hotter objects. Feige55 shows also several other heavier elements. The measured abundances agree only partly with the predictions of the radiative levitation theory. We attribute this to current limitations of the models and the probable presence of mass loss. In the spectrum of CD -38 10980, we have observed the quasi-molecular satellites of Lbeta. This confirms theoretical predictions about the visibility range for these features.

We present new FUV/UV observations of the DA white dwarf Wolf 1346 obtained with the Hopkins Ultraviolet Telescope. The atmospheric parameters of this object are estimated from a fit of model atmospheres to several optical spectra to be Teff = 20000 K, log g = 7.90. From the optical spectrum this star is a normal DA without any indications for chemical elements other than hydrogen. The hydrogen line L beta, however, shows a very unusual shape, with a steep red wing and two absorption features on this wing. The shape is reminiscent of the effects of quasi-molecular line broadening, as observed in L alpha in cooler DA white dwarfs. We show that this is indeed the correct explanation, by identifying 4 quasi-molecular satellites caused through perturbations by the H+ ion (H2+ quasi-molecule). The steep red wing is caused by the exponential decline of the line profile beyond the satellite most distant from the line center at 1078 A.