We present a detailed analysis of the white dwarf luminosity functions derived from the local 40 pc sample and the deep proper motion catalog of Munn et al (2014, 2017). Many of the previous studies ignored the contribution of thick disk white dwarfs to the Galactic disk luminosity function, which results in an erronous age measurement. We demonstrate that the ratio of thick/thin disk white dwarfs is roughly 20\% in the local sample. Simultaneously fitting for both disk components, we derive ages of 6.8-7.0 Gyr for the thin disk and 8.7 $\pm$ 0.1 Gyr for the thick disk from the local 40 pc sample. Similarly, we derive ages of 7.4-8.2 Gyr for the thin disk and 9.5-9.9 Gyr for the thick disk from the deep proper motion catalog, which shows no evidence of a deviation from a constant star formation rate in the past 2.5 Gyr. We constrain the time difference between the onset of star formation in the thin disk and the thick disk to be $1.6^{+0.3}_{-0.4}$ Gyr. The faint end of the luminosity function for the halo white dwarfs is less constrained, resulting in an age estimate of $12.5^{+1.4}_{-3.4}$ Gyr for the Galactic inner halo. This is the first time ages for all three major components of the Galaxy are obtained from a sample of field white dwarfs that is large enough to contain significant numbers of disk and halo objects. The resultant ages agree reasonably well with the age estimates for the oldest open and globular clusters.

A catalog of 8472 white dwarf (WD) candidates is presented, selected using reduced proper motions from the deep proper motion catalog of Munn et al. 2014. Candidates are selected in the magnitude range 16 < r < 21.5 over 980 square degrees, and 16 < r < 21.3 over an additional 1276 square degrees, within the Sloan Digital Sky Survey (SDSS) imaging footprint. Distances, bolometric luminosities, and atmospheric compositions are derived by fitting SDSS ugriz photometry to pure hydrogen and helium model atmospheres (assuming surface gravities log g = 8). The disk white dwarf luminosity function (WDLF) is constructed using a sample of 2839 stars with 5.5 < M_bol < 17, with statistically significant numbers of stars cooler than the turnover in the luminosity function. The WDLF for the halo is also constructed, using a sample of 135 halo WDs with 5 < M_bol < 16. We find space densities of disk and halo WDs in the solar neighborhood of 5.5 +- 0.1 x 10^-3 pc^-3 and 3.5 +- 0.7 x 10^-5 pc^-3, respectively. We resolve the bump in the disk WDLF due to the onset of fully convective envelopes in WDs, and see indications of it in the halo WDLF as well.

BASE-9 is a Bayesian software suite that recovers star cluster and stellar parameters from photometry. BASE-9 is useful for analyzing single-age, single-metallicity star clusters, binaries, or single stars, and for simulating such systems. BASE-9 uses Markov chain Monte Carlo and brute-force numerical integration techniques to estimate the posterior probability distributions for the age, metallicity, helium abundance, distance modulus, and line-of-sight absorption for a cluster, and the mass, binary mass ratio, and cluster membership probability for every stellar object. BASE-9 is provided as open source code on a version-controlled web server. The executables are also available as Amazon Elastic Compute Cloud images. This manual provides potential users with an overview of BASE-9, including instructions for installation and use.

We present deep photometric observations of the open cluster NGC 2477 using HST/WFPC2. By identifying seven cluster white dwarf candidates, we present an analysis of the white dwarf age of this cluster, using both the traditional method of fitting isochrones to the white dwarf cooling sequence, and by employing a new Bayesian statistical technique that has been developed by our group. This new method performs an objective, simultaneous model fit of the cluster and stellar parameters (namely age, metallicity, distance, reddening, as well as individual stellar masses, mass ratios, and cluster membership) to the photometry. Based on this analysis, we measure a white dwarf age of 1.035 +/- 0.054 +/- 0.087 Gyr (uncertainties represent the goodness of model fits and discrepancy among models, respectively), in good agreement with the cluster's main sequence turnoff age. This work is part of our ongoing work to calibrate main sequence turnoff and white dwarf ages using open clusters, and to improve the precision of cluster ages to the ~5% level.

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 infrared observations in search of a planet around the white dwarf, GD66. Time-series photometry of GD66 shows a variation in the arrival time of stellar pulsations consistent with the presence of a planet with mass > 2.4Mj. Any such planet is too close to the star to be resolved, but the planet's light can be directly detected as an excess flux at 4.5um. We observed GD66 with the two shorter wavelength channels of IRAC on Spitzer but did not find strong evidence of a companion, placing an upper limit of 5--7Mj on the mass of the companion, assuming an age of 1.2--1.7Gyr.

We present the results of a search for pulsations in six of the recently discovered carbon-atmosphere white dwarf ("hot DQ") stars. Based on our theoretical calculations, the star SDSS J142625.71+575218.3 is the only object expected to pulsate. We observe this star to be variable, with significant power at 417.7 s and 208.8 s (first harmonic), making it a strong candidate as the first member of a new class of pulsating white dwarf stars, the DQVs. Its folded pulse shape, however, is quite different from that of other white dwarf variables, and shows similarities with that of the cataclysmic variable AM CVn, raising the possibility that this star may be a carbon-transferring analog of AM CVn stars. In either case, these observations represent the discovery of a new and exciting class of object.

We present limits on planetary companions to pulsating white dwarf stars. A subset of these stars exhibit extreme stability in the period and phase of some of their pulsation modes; a planet can be detected around such a star by searching for periodic variations in the arrival time of these pulsations. We present limits on companions greater than a few Jupiter masses around a sample of 15 white dwarf stars as part of an on-going survey. One star shows a variation in arrival time consistent with a 2 M_J planet in a 4.5 year orbit. We discuss other possible explanations for the observed signal and conclude that a planet is the most plausible explanation based on the data available.

We present the first phase in our ongoing work to use Sloan Digital Sky Survey (SDSS) data to create separate white dwarf (WD) luminosity functions for two or more different mass ranges. In this paper, we determine the completeness of the SDSS spectroscopic white dwarf sample by comparing a proper-motion selected sample of WDs from SDSS imaging data with a large catalog of spectroscopically determined WDs. We derive a selection probability as a function of a single color (g-i) and apparent magnitude (g) that covers the range -1.0 < g-i < 0.2 and 15 < g < 19.5. We address the observed upturn in log g for white dwarfs with Teff <~ 12,000K and offer arguments that the problem is limited to the line profiles and is not present in the continuum. We offer an empirical method of removing the upturn, recovering a reasonable mass function for white dwarfs with Teff < 12,000K. Finally, we present a white dwarf luminosity function with nearly an order of magnitude (3,358) more spectroscopically confirmed white dwarfs than any previous work.

Currently there are two main techniques for independently determining the ages of stellar populations: main sequence evolution theory (via cluster isochrones) and white dwarf cooling theory. Open clusters provide the ideal environment for the calibration of these two clocks. Because current techniques to derive cluster ages from white dwarfs are observationally challenging, we discuss the feasibility of determining white dwarf ages from the brighter white dwarfs alone. This would eliminate the requirement of observing the coolest (i.e., faintest) white dwarfs. We discuss our method for testing this new idea, as well as the required photometric precision and prior constraints on metallicity, distance, and reddening. We employ a new Bayesian statistical technique to obtain and interpret results.

We demonstrate a new Bayesian technique to invert color-magnitude diagrams of main sequence and white dwarf stars to reveal the underlying cluster properties of age, distance, metallicity, and line-of-sight absorption, as well as individual stellar masses. The advantages our technique has over traditional analyses of color-magnitude diagrams are objectivity, precision, and explicit dependence on prior knowledge of cluster parameters. Within the confines of a given set of often-used models of stellar evolution, the initial-final mass relation, and white dwarf cooling, and assuming photometric errors that one could reasonably achieve with the Hubble Space Telescope, our technique yields exceptional precision for even modest numbers of cluster stars. For clusters with 50 to 400 members and one to a few dozen white dwarfs, we find typical internal errors of sigma([Fe/H]) <= 0.03 dex, sigma(m-M) <= 0.02 mag, and sigma(Av) <= 0.01 mag. We derive cluster white dwarf ages with internal errors of typically only 10% for clusters with only three white dwarfs and almost always <= 5% with ten white dwarfs. These exceptional precisions will allow us to test white dwarf cooling models and standard stellar evolution models through observations of white dwarfs in open and globular clusters.