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$.

Aims: Hydrogen deficient white dwarfs are characterized by very dense, fluid-like atmospheres of complex physics and chemistry that are still poorly understood. The incomplete description of these atmospheres by the models results in serious problems with the description of spectra of these stars and subsequent difficulties in derivation of their surface parameters. Here, we address the problem of infrared (IR) opacities in the atmospheres of cool white dwarfs by direct $ab$ $initio$ simulations of IR absorption of dense helium. Methods: We applied state-of-the-art density functional theory-based quantum molecular dynamics simulations to obtain the time evolution of the induced dipole moment. The IR absorption coefficients were obtained by the Fourier transform of the dipole moment time autocorrelation function. Results: We found that a dipole moment is induced due to three- and more-body simultaneous collisions between helium atoms in highly compressed helium. This results in a significant IR absorption that is directly proportional to $\rm \rho_{\rm He}^3$, where $\rho_{\rm He}$ is the density of helium. To our knowledge, this absorption mechanism has never been measured or computed before and is therefore not accounted for in the current atmosphere models. It should dominate the other collisionally induced absorptions (CIA), arising from $\rm H-He$ and $\rm H_2-He$ pair collisions, and therefore shape the IR spectra of helium-dominated and pure helium atmosphere cool white dwarfs for $\rm He/H>10^4$. Conclusions: Our work shows that there exists an unaccounted IR absorption mechanism arising from the multi-collisions between He atoms in the helium-rich atmospheres of cool white dwarfs, including pure helium atmospheres. This absorption may be responsible for a yet unexplained frequency dependence of near- and mid- IR spectra of helium-rich stars.

The atmospheres of very cool, hydrogen-rich white dwarfs (Teff <6000 K) are challenging to models because of the increased complexity of the equation of state, chemical equilibrium, and opacity sources in a low-temperature, weakly ionized dense gas. In particular, many models that assume relatively simple models for the broadening of atomic levels and mostly ideal gas physics overestimate the flux in the blue part of their spectra. A solution to this problem that has met with some success is that additional opacity at short wavelengths comes for the extreme broadening of the Lyman alpha line of atomic H by collisions primarily with H2. For the purpose of validating this model more rigorously, we acquired Hubble Space Telescope STIS spectra of 8 very cool white dwarfs (5 DA and 3 DC stars). Combined with their known parallaxes, BVRIJHK and Spitzer IRAC photometry, we analyze their entire spectral energy distribution (from 0.24 to 9.3 micron) with a large grid of model atmospheres and synthetic spectra. We find that the red wing of the Lyman alpha line reproduces the rapidly decreasing near-UV flux of these very cool stars very well. We determine better constrained values of Teff and gravity as well as upper limits to the helium abundance in their atmospheres.

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.

Cool white dwarfs with Teff < 6000 K are the remnants of the oldest stars that existed in our Galaxy. Their atmospheres, when properly characterized, can provide valuable information on white dwarf evolution and ultimately star formation through the history of the Milky Way. Understanding the atmospheres of these stars requires joined observational effort and reliable atmosphere modeling. We discuss and analyze recent observations of the near-ultraviolet (UV) and near-infrared (IR) spectrum of several cool white dwarfs including DQ/DQp stars showing carbon in their spectra. We present fits to the entire spectral energy distribution (SED) of selected cool stars, showing that the current pure-hydrogen atmosphere models are quite reliable, especially in the near-UV spectral region. Recently, we also performed an analysis of the coolest known DQ/DQp stars investigating further the origin of the C2 Swan bands-like spectral features that characterize the DQp stars. We show that the carbon abundances derived for DQp stars fit the trend of carbon abundance with Teff seen in normal cool DQ stars. This further supports the recent conclusion of Kowalski A&A (2010) that DQp stars are DQ stars with pressure distorted Swan bands. However, we encounter some difficulties in reproducing the IR part of the SED of stars having a mixed He/H atmosphere. This indicates limitations in current models of the opacity in dense He/H fluids.

Since there are several ways planets can survive the giant phase of the host star, we examine the habitability and detection of planets orbiting white dwarfs. As a white dwarf cools from 6000 K to 4000 K, a planet orbiting at 0.01 AU would remain in the Continuous Habitable Zone (CHZ) for ~8 Gyr. We show that photosynthetic processes can be sustained on such planets. The DNA-weighted UV radiation dose for an Earth-like planet in the CHZ is less than the maxima encountered on Earth, hence non-magnetic white dwarfs are compatible with the persistence of complex life. Polarisation due to a terrestrial planet in the CHZ of a cool white dwarf is 10^2 (10^4) times larger than it would be in the habitable zone of a typical M-dwarf (Sun-like star). Polarimetry is thus a viable way to detect close-in rocky planets around white dwarfs. Multi-band polarimetry would also allow reveal the presence of a planet atmosphere, providing a first characterisation. Planets in the CHZ of a 0.6 M_sun white dwarf will be distorted by Roche geometry, and a Kepler-11d analogue would overfill its Roche lobe. With current facilities a Super-Earth-sized atmosphereless planet is detectable with polarimetry around the brightest known cool white dwarf. Planned future facilities render smaller planets detectable, in particular by increasing the instrumental sensitivity in the blue.

We present a detailed model atmosphere analysis of two of the oldest stars known in the solar neighborhood, the high proper motion white dwarfs SDSS J110217.48+411315.4 (hereafter J1102) and WD 0346+246 (hereafter WD0346). We present trigonometric parallax observations of J1102, which places it at a distance of only 33.7 +- 2.0 pc. Based on the state of the art model atmospheres, optical, near-, mid-infrared photometry, and distances, we constrain the temperatures, atmospheric compositions, masses, and ages for both stars. J1102 is an 11 Gyr old (white dwarf plus main-sequence age), 0.62 Msol white dwarf with a pure H atmosphere and Teff = 3830 K. WD0346 is an 11.5 Gyr old, 0.77 Msol white dwarf with a mixed H/He atmosphere and Teff = 3650 K. Both stars display halo kinematics and their ages agree remarkably well with the ages of the nearest globular clusters, M4 and NGC 6397. J1102 and WD0346 are the closest examples of the oldest halo stars that we know of.

We present extensive spectral and photometric observations of the recycled pulsar/white-dwarf binary containing PSR J0437-4715, which we analyzed together with archival X-ray and gamma-ray data, to obtain the complete mid-infrared to gamma-ray spectrum. We first fit each part of the spectrum separately, and then the whole multi-wavelength spectrum. We find that the optical-infrared part of the spectrum is well fit by a cool white dwarf atmosphere model with pure hydrogen composition. The model atmosphere (Teff = 3950pm150K, log g=6.98pm0.15, R_WD=(1.9pm0.2)e9 cm) fits our spectral data remarkably well for the known mass and distance (M=0.25pm0.02Msun, d=156.3pm1.3pc), yielding the white dwarf age (tau=6.0pm0.5Gyr). In the UV, we find a spectral shape consistent with thermal emission from the bulk of the neutron star surface, with surface temperature between 1.25e5 and 3.5e5K. The temperature of the thermal spectrum suggests that some heating mechanism operates throughout the life of the neutron star. The temperature distribution on the neutron star surface is non-uniform. In the X-rays, we confirm the presence of a high-energy tail which is consistent with a continuation of the cut-off power-law component (Gamma=1.56pm0.01, Ecut=1.1pm0.2GeV) that is seen in gamma-rays and perhaps even extends to the near-UV.

The mass-dependent equilibrium stable isotope fractionation between different materials is an important geochemical process. Here we present an efficient method to compute the isotope fractionation between complex minerals and fluids at high pressure, P, and temperature, T, representative for the Earth's crust and mantle. The method is tested by computation of the equilibrium fractionation of lithium isotopes between aqueous fluids and various Li bearing minerals such as staurolite, spodumene and mica. We are able to correctly predict the direction of the isotope fractionation as observed in the experiments. On the quantitative level the computed fractionation factors agree within 1.0 permil with the experimental values indicating predictive power of ab initio methods. We show that with ab initio methods we are able to investigate the underlying mechanisms driving the equilibrium isotope fractionation process, such as coordination of the fractionating elements, their bond strengths to the neighboring atoms, compression of fluids and thermal expansion of solids. This gives valuable insight into the processes governing the isotope fractionation mechanisms on the atomic scale. The method is applicable to any state and does not require different treatment of crystals and fluids.

Reliable modeling of the atmospheres of cool white dwarfs is crucial for understanding the atmospheric evolution of these stars and for accurate white dwarfs cosmochronology. Over the last decade {\it ab initio} modeling entered many research fields and has been successful in predicting properties of various materials under extreme conditions. In many cases the investigated physical regimes are difficult or even impossible to access by experimental methods, and first principles quantum mechanical calculations are the only tools available for investigation. Using modern methods of computational chemistry and physics we investigate the atmospheres of helium-rich, old white dwarfs. Such atmospheres reach extreme, fluid like densities (up to grams per cm$^3$) and represent an excellent laboratory for high temperature and pressure physics and chemistry. We show our results for the stability and opacity of $\rm H^-$ and $\rm C_2$ in dense helium and the implications of our work for understanding cool white dwarfs.

Over the last decade {\it ab initio} modeling of material properties has become widespread in diverse fields of research. It has proved to be a powerful tool for predicting various properties of matter under extreme conditions. We apply modern computational chemistry and materials science methods, including density functional theory (DFT), to solve lingering problems in the modeling of the dense atmospheres of cool white dwarfs ($T_{\rm eff}\rm <7000 \, K$). Our work on the revision and improvements of the absorption mechanisms in the hydrogen and helium dominated atmospheres resulted in a new set of atmosphere models. By inclusion of the Ly-$\rm \alpha$ red wing opacity we successfully fitted the entire spectral energy distributions of known cool DA stars. In the subsequent work we fitted the majority of the coolest stars with hydrogen-rich models. This finding challenges our understanding of the spectral evolution of cool white dwarfs. We discuss a few examples, including the cool companion to the pulsar PSR J0437-4715. The two problems important for the understanding of cool white dwarfs are the behavior of negative hydrogen ion and molecular carbon in a fluid-like, helium dominated medium. Using {\it ab initio} methods we investigate the stability and opacity of these two species in dense helium. Our investigation of $\rm C_2$ indicates that the absorption features observed in the ``peculiar'' DQp white dwarfs resemble the absorption of perturbed $\rm C_2$ in dense helium.

The DQ white dwarfs are stars whose atmosphere is enriched with carbon, which for cool stars ($T_{\rm eff}<8000\rm \, K$) is indicated by the Swan bands of $\rm C_2$ in the optical part of their spectra. With decreasing effective temperature these molecular bands undergo a significant blueshift ($\sim 100-300 \AA$). The origin of this phenomenon has been disputed over the last two decades and has remained unknown. We attempt to address this problem by investigating the impact of dense helium on the spectroscopic properties of molecular carbon under the physical conditions encountered inside helium-rich, fluid-like atmospheres of cool DQ white dwarfs. We found that the electronic transition energy $T_e$ increases monotonically with the helium density ($\Delta T_{\rm e}\rm\, (eV)\sim1.6 \, \it \rho \rm \, (g/cm^3)$). This causes the Swan absorption to occur at shorter wavelengths compared with unperturbed $\rm C_2$. On the other hand the pressure-induced increase in the vibrational frequency is insufficient to account for the observed Swan bands shifts. This is consistent with the observations and indicates that the observed Swan-like molecular bands are most likely the pressure-shifted bands of $\rm C_2$.

We present a new set of cooling models and isochrones for both H- and He-atmosphere white dwarfs, incorporating accurate boundary conditions from detailed model atmosphere calculations, and carbon-oxygen chemical abundance profiles based on updated stellar evolution calculations from the BaSTI stellar evolution archive - a theoretical data center for the Virtual Observatory. We discuss and quantify the uncertainties in the cooling times predicted by the models, arising from the treatment of mixing during the central H- and He-burning phases, number of thermal pulses experienced by the progenitors, progenitor metallicity and the $^{12}C(\alpha,\gamma)^{16}O$ reaction rate. The largest sources of uncertainty turn out to be related to the treatment of convection during the last stages of the progenitor central He-burning phase, and the $^{12}C(\alpha,\gamma)^{16}O$ reaction rate. We compare our new models to previous calculations performed with the same stellar evolution code, and discuss their application to the estimate of the age of the solar neighborhood, and the interpretation of the observed number ratios between H- and He-atmosphere white dwarfs. The new white dwarf sequences and an extensive set of white dwarf isochrones that cover a large range of ages and progenitor metallicities are made publicly available at the official BaSTI website.

We report the discovery of three nearby old halo white dwarf candidates in the Sloan Digital Sky Survey (SDSS), including two stars in a common proper motion binary system. These candidates are selected from our 2800 square degree proper motion survey on the Bok and U.S. Naval Observatory Flagstaff Station 1.3m telescopes, and they display proper motions of 0.4-0.5 arcsec/yr. Follow-up MMT spectroscopy and near-infrared photometry demonstrate that all three objects are hydrogen-dominated atmosphere white dwarfs with Teff = 3700 - 4100 K. For average mass white dwarfs, these temperature estimates correspond to cooling ages of 9-10 Gyr, distances of 70-80 pc, and tangential velocities of 140-200 km/s. Based on the UVW space velocities, we conclude that they most likely belong to the halo. Furthermore, the combined main-sequence and white dwarf cooling ages are 10-11 Gyr. Along with SDSS J1102+4113, they are the oldest field white dwarfs currently known. These three stars represent only a small fraction of the halo white dwarf candidates in our proper motion survey, and they demonstrate that deep imaging surveys like the Pan-STARRS and Large Synoptic Survey Telescope should find many old thick disk and halo white dwarfs that can be used to constrain the age of the Galactic thick disk and halo.

We present Spitzer 5-15 micron spectroscopy of one cool white dwarf and 3.6-8 micron photometry of 51 cool white dwarfs with T_eff < 6000 K. The majority of our targets have accurate BVRIJHK photometry and trigonometric parallax measurements available, which enables us to perform a detailed model atmosphere analysis using their optical, near- and mid-infrared photometry with state- of-the-art model atmospheres. We demonstrate that the optical and infrared spectral energy distributions of cool white dwarfs are well reproduced by our grid of models. Our best fit models are consistent with the observations within 5% in all filters except the IRAC 8 micron band, which has the lowest signal- to-noise ratio photometry. Excluding the ultracool white dwarfs, none of the stars in our sample show significant mid-infrared flux deficits or excesses. The non-detection of mid-infrared excess flux around our 2-9 Gyr old targets constrain the fraction of cool white dwarfs with warm debris disks to 0.8% (+1.5% -0.8%).

We report the discovery of a nearby, old, halo white dwarf candidate from the Sloan Digital Sky Survey. SDSS J110217.48+411315.4 has a proper motion of 1.75 arcsec/year and redder optical colors than all other known featureless (type DC) white dwarfs. We present SDSS imaging and spectroscopy of this object, along with near-infrared photometry obtained at the United Kingdom Infra-Red Telescope. Fitting its photometry with up-to-date model atmospheres, we find that its overall spectral energy distribution is fit reasonably well with a pure hydrogen composition and T_eff~3800 K (assuming log g=8). That temperature and gravity would place this white dwarf at 35 pc from the Sun with a tangential velocity of 290 km/s and space velocities consistent with halo membership; furthermore, its combined main sequence and white dwarf cooling age would be ~11 Gyr. However, if this object is a massive white dwarf, it could be a younger object with a thick disk origin. Whatever its origin, the optical colors of this object are redder than predicted by any current pure hydrogen, pure helium or mixed hydrogen-helium atmospheric model, indicating that there remain problems in our understanding of the complicated physics of the dense atmospheres of cool white dwarfs.

We present the first mid-infrared spectra of two cool white dwarfs obtained with the Spitzer Space Telescope. We also present 3.5-8 micron photometry for 19 cool white dwarfs with 5000K < Teff < 9000K. We perform a detailed model atmosphere analysis of these white dwarfs by fitting their UBVRIJHK and Spitzer photometry with state-of-the-art model atmospheres, and demonstrate that the optical and infrared spectral energy distributions of cool white dwarfs are well reproduced by our grid of models. Our mid-IR photometry and 7.5-14.5 micron spectrum of WD0018-267 are consistent with a Teff = 5720K, pure hydrogen white dwarf model atmosphere. On the other hand, LHS 1126 remains peculiar with significant mid-IR flux deficits in all IRAC bands and a featureless spectrum in the 5.2-7.5 micron range. Even though this deficit is attributed to collision induced absorption (CIA) due to molecular hydrogen, the shape of the deficit cannot be explained with current CIA opacity calculations. The infrared portion of the LHS 1126 spectral energy distribution is best-fit with a power law index of -1.99; identical to a Rayleigh-Jeans spectrum. This argues that the deficit may be due to an unrecognized grey-like opacity source in the infrared.

We discuss the importance of pure hydrogen white dwarf atmosphere models with Ly-$\rm \alpha$ far red wing opacity in the analysis of the white dwarf cooling sequence of the globular cluster NGC 6397. Our recently improved atmosphere models account for the previously missing opacity from the Ly-$\rm \alpha$ hydrogen line broadened by collisions of the absorbing hydrogen atoms with molecular and atomic hydrogen. These models are the first that well reproduce the UV colors and spectral energy distributions of cool white dwarfs with $T_{\rm eff}<6000 \rm K$ observed in the Galactic Disk. Fitting the observed $F814W$ magnitude and $F606W-F814W$ color we obtained a value for the true distance modulus, $\mu=12.00 \pm 0.02$, that is in agreement with recent analyses. We show that the stars at the end of the cooling sequence appear to be $\rm \sim 160 K$ cooler when models that account for Ly-$\rm \alpha$ opacity are used. This indicates that the age of NGC 6397 derived from the white dwarf cooling sequence using atmosphere models that do not include the correct Ly-$\alpha$ opacity is underestimated by $\sim 0.5$ Gyr. Our analysis shows that it is essential to use white dwarf atmosphere models with Ly-$\rm \alpha$ opacity for precise dating of old stellar populations from white dwarf cooling sequences.

We investigate the importance of the far red wing of the Lyman $ \alpha$ line of hydrogen in the atmospheres of cool white dwarfs of pure hydrogen composition. We find that this absorption process dominates all important sources of opacity in the blue part of the optical spectrum of these stars. Our successful fits to the spectra of cool DA/DC white dwarfs indicate that the far red wing of the $\rm Ly \alpha$ line is the source of opacity that had been missing in the models. The observed sequence of cool white dwarfs in color-color diagrams is very well reproduced by our new pure hydrogen atmosphere models, suggesting that the atmospheric composition of the coolest DC white dwarfs must be revisited.

We investigate the importance of the pseudo-continuum bound-free opacity from hydrogen atoms in the atmospheres of cool white dwarfs. This source of absorption, when calculated by the occupation probability formalism applied in the modeling of white dwarf atmospheres with $T_{\rm eff}\rm <17000 K$, dominates all other sources of opacity at optical wavelengths. This is unrealistic and not observed. On the other hand, a significant flux suppression in the blue part of the spectra of cool white dwarfs has been reported, and mainly interpreted as a result of the pseudo-continuum absorption from atomic hydrogen. We investigate this problem by proposing a new, more realistic approach to calculating this source of opacity. We show that this absorption is orders of magnitude smaller than that predicted by current methods. Therefore, we rule out the pseudo-continuum opacity as a source of the flux deficiency observed in the spectra of cool white dwarfs.

We investigate the dissociation equilibrium of $\rm H_2$ in very cool, helium-rich white dwarf atmospheres. We present the solution of the non-ideal chemical equilibrium for the dissociation of molecular hydrogen in a medium of dense helium. We find that at the photosphere of cool white dwarfs of $T_{\rm eff}\rm=4000 K$, the non-ideality results in an increase of the mole fraction of molecular hydrogen by up to a factor of $\sim 10$, compared to the equilibrium value for the ideal gas. This increases the $\rm H_{2}-He$ CIA opacity by an order of magnitude and will affect the determination of the abundance of hydrogen in very cool, helium-rich white dwarfs.

The atmospheres of cool, helium-rich white dwarfs constitute an exotic and poorly explored physical regime of stellar atmospheres. Under physical conditions where the temperature varies from $\rm 1000K$ to $\rm 10000K$, the density can reach values as large as $\rm 2 g/cm^{3}$, and the pressure is as high as $\rm 1 Mbar$, the atmosphere is no longer an ideal gas and must be treated as a dense fluid. Helium atoms become strongly correlated and refraction effects are present. Opacity sources such as $\rm He^{-}$ free-free absorption must be calculated with a formalism that has never been applied to astrophysical opacities. These effects have been ignored in previous models of cool white dwarf atmospheres.

We consider the problem of radiative transfer in stellar atmospheres where the index of refraction departs from unity and is a function of density and temperature. We present modified Feautrier and Lambda-iteration methods to solve the equation of radiative transfer with refraction in a plane parallel atmosphere. These methods are general and can be used in any problem with 1-D geometry where the index of refraction is a monotonically varying function of vertical optical depth. We present an application to very cool white dwarf atmospheres where the index of refraction departs significantly from unity. We investigate how ray curvature and total internal reflection affect the limb darkening and the pressure-temperature structure of the atmosphere. Refraction results in a much weakened limb darkening effect. We find that through the constraint of radiative equilibrium, total internal reflection warms the white dwarf atmosphere near the surface ($\tau \wig< 1$). This effect may have a significant impact on studies of very cool white dwarf stars.