We have made high precision polarimetric observations of the polluted white dwarf G29-38 with the HIgh Precision Polarimetric Instrument 2. The observations were made at two different observatories -- using the 8.1-m Gemini North Telescope and the 3.9-m Anglo AustralianTelescope -- and are consistent with each other. After allowing for a small amount of interstellar polarization, the intrinsic linear polarization of the system is found to be 275.3 +/- 31.9 parts-per-million at a position angle of 90.8 +/- 3.8 degrees in the SDSS g' band. We compare the observed polarization with the predictions of circumstellar disc models. The measured polarization is small in the context of the models we develop which only allows us to place limits on disc inclination and Bond albedo for optically thin disc geometries. In this case either the inclination is near face-on or the albedo is small -- likely in the range 0.05 to 0.15 -- which is in line with other debris disc measurements. A preliminary search for the effects of G29-38's pulsations in the polarization signal produced inconsistent results. This may be caused by beating effects, indicate a clumpy dust distribution, or be a consequence of measurement systematics.

The inwards scattering of planetesimals towards white dwarfs is expected to be a stochastic process with variability on human time-scales. The planetesimals tidally disrupt at the Roche radius, producing dusty debris detectable as excess infrared emission. When sufficiently close to the white dwarf, this debris sublimates and accretes on to the white dwarf and pollutes its atmosphere. Studying this infrared emission around polluted white dwarfs can reveal how this planetary material arrives in their atmospheres. We report a near-infrared monitoring campaign of 34 white dwarfs with infrared excesses with the aim to search for variability in the dust emission. Time series photometry of these white dwarfs from the United Kingdom Infrared Telescope (Wide Field Camera) in the J, H and K bands were obtained over baselines of up to three years. We find no statistically significant variation in the dust emission in all three near-infrared bands. Specifically, we can rule out variability at ~1.3% for the 13 white dwarfs brighter than 16th mag in K band, and at ~10% for the 32 white dwarfs brighter than 18th mag over time-scales of three years. Although to date two white dwarfs, SDSS J095904.69-020047.6 and WD 1226+110, have shown K band variability, in our sample we see no evidence of new K band variability at these levels. One interpretation is that the tidal disruption events which lead to large variabilities are rare, occur on short time-scales, and after a few years the white dwarfs return to being stable in the near-infrared.

The binary fractions of open and globular clusters yield powerful constraints on their dynamical state and evolutionary history. We apply publicly available Bayesian analysis code to a $UBVRIJHK_{S}$ photometric catalog of the old open cluster NGC 188 to detect and characterize photometric binaries along the cluster main sequence. This technique has the advantage that it self-consistently handles photometric errors, missing data in various bandpasses, and star-by-star prior constraints on cluster membership. Simulations are used to verify uncertainties and quantify selection biases in our analysis, illustrating that among binaries with mass ratios >0.5, we recover the binary fraction to better than 7% in the mean, with no significant dependence on binary fraction and a mild dependence on assumed mass ratio distribution. Using our photometric catalog, we recover the majority (65%$\pm$11%) of spectroscopically identified main sequence binaries, including 8 of the 9 with spectroscopically measured mass ratios. Accounting for incompleteness and systematics, we derive a mass ratio distribution that rises toward lower mass ratios (within our $q >$0.5 analysis domain). We observe a raw binary fraction for solar-type main sequence stars with mass ratios $q >$0.5 of 42%$\pm$4%, independent of the assumed mass ratio distribution to within its uncertainties, consistent with literature values for old open clusters but significantly higher than the field solar-type binary fraction. We confirm that the binaries identified by our method are more concentrated than single stars, in agreement with previous studies, and we demonstrate that the binary nature of those candidates which remain unidentified spectroscopically is strongly supported by photometry from Gaia DR2.

The infrared dust emission from the white dwarf GD 56 is found to rise and fall by 20% peak-to-peak over 11.2 yr, and is consistent with ongoing dust production and depletion. It is hypothesized that the dust is produced via collisions associated with an evolving dust disk, temporarily increasing the emitting surface of warm debris, and is subsequently destroyed or assimilated within a few years. The variations are consistent with debris that does not change temperature, indicating that dust is produced and depleted within a fixed range of orbital radii. Gas produced in collisions may rapidly re-condense onto grains, or may accrete onto the white dwarf surface on viscous timescales that are considerably longer than Poynting-Robertson drag for micron-sized dust. This potential delay in mass accretion rate change is consistent with multi-epoch spectra of the unchanging Ca II and Mg II absorption features in GD 56 over 15 yr, although the sampling is sparse. Overall these results indicate that collisions are likely to be the source of dust and gas, either inferred or observed, orbiting most or all polluted white dwarfs.

The initial-final mass relation (IFMR) of white dwarfs (WDs) plays an important role in stellar evolution. To derive precise estimates of IFMRs and explore how they may vary among star clusters, we propose a Bayesian hierarchical model that pools photo- metric data from multiple star clusters. After performing a simulation study to show the benefits of the Bayesian hierarchical model, we apply this model to five star clus- ters: the Hyades, M67, NGC 188, NGC 2168, and NGC 2477, leading to reasonable and consistent estimates of IFMRs for these clusters. We illustrate how a cluster-specific analysis of NGC 188 using its own photometric data can produce an unreasonable IFMR since its WDs have a narrow range of zero-age main sequence (ZAMS) masses. However, the Bayesian hierarchical model corrects the cluster-specific analysis by bor- rowing strength from other clusters, thus generating more reliable estimates of IFMR parameters. The data analysis presents the benefits of Bayesian hierarchical modelling over conventional cluster-specific methods, which motivates us to elaborate the pow- erful statistical techniques in this article.

In astrophysics, we often aim to estimate one or more parameters for each member object in a population and study the distribution of the fitted parameters across the population. In this paper, we develop novel methods that allow us to take advantage of existing software designed for such case-by-case analyses to simultaneously fit parameters of both the individual objects and the parameters that quantify their distribution across the population. Our methods are based on Bayesian hierarchical modelling which is known to produce parameter estimators for the individual objects that are on average closer to their true values than estimators based on case-by-case analyses. We verify this in the context of estimating ages of Galactic halo white dwarfs (WDs) via a series of simulation studies. Finally, we deploy our new techniques on optical and near-infrared photometry of ten candidate halo WDs to obtain estimates of their ages along with an estimate of the mean age of Galactic halo WDs of [11.25, 12.96] Gyr. Although this sample is small, our technique lays the ground work for large-scale studies using data from the Gaia mission.

Observational tests of stellar and Galactic chemical evolution call for the joint knowledge of a star's physical parameters, detailed element abundances, and precise age. For cool main-sequence (MS) stars the abundances of many elements can be measured from spectroscopy, but ages are very hard to determine. The situation is different if the MS star has a white dwarf (WD) companion and a known distance, as the age of such a binary system can then be determined precisely from the photometric properties of the cooling WD. As a pilot study for obtaining precise age determinations of field MS stars, we identify nearly one hundred candidates for such wide binary systems: a faint WD whose GPS1 proper motion matches that of a brighter MS star in Gaia/TGAS with a good parallax ($\sigma_\varpi/\varpi\le 0.05$). We model the WD's multi-band photometry with the BASE-9 code using this precise distance (assumed to be common for the pair) and infer ages for each binary system. The resulting age estimates are precise to $\le 10\%$ ($\le 20\%$) for $42$ ($67$) MS-WD systems. Our analysis more than doubles the number of MS-WD systems with precise distances known to date, and it boosts the number of such systems with precise age determination by an order of magnitude. With the advent of the Gaia DR2 data, this approach will be applicable to a far larger sample, providing ages for many MS stars (that can yield detailed abundances for over 20 elements), especially in the age range 2 to 8\,\Gyr, where there are only few known star clusters.

The rapidly evolving dust and gas extinction observed towards WD 1145+017 has opened a real-time window onto the mechanisms for destruction-accretion of planetary bodies onto white dwarf stars, and has served to underline the importance of considering the dynamics of dust particles around such objects. Here it is argued that the interaction between (charged) dust grains and the stellar magnetic field is an important ingredient in understanding the physical distribution of infrared emitting particles in the vicinity of such white dwarfs. These ideas are used to suggest a possible model for WD 1145+017 in which the unusual transit shapes are caused by opaque clouds of dust trapped in the stellar magnetosphere. The model can account for the observed transit periodicities if the stellar rotation is near 4.5 h, as the clouds of trapped dust are then located near or within the co-rotation radius. The model requires the surface magnetic field to be at least around some tens of kG. In contrast to the eccentric orbits expected for large planetesimals undergoing tidal disintegration, the orbits of magnetospherically-trapped dust clouds are essentially circular, consistent with the observations.

We use Hubble Space Telescope (HST) imaging from the ACS Treasury Survey to determine fits for single population isochrones of 69 Galactic globular clusters. Using robust Bayesian analysis techniques, we simultaneously determine ages, distances, absorptions, and helium values for each cluster under the scenario of a "single" stellar population on model grids with solar ratio heavy element abundances. The set of cluster parameters is determined in a consistent and reproducible manner for all clusters using the Bayesian analysis suite BASE-9. Our results are used to re-visit the age-metallicity relation. We find correlations with helium and several other parameters such as metallicity, binary fraction, and proxies for cluster mass. The helium abundances of the clusters are also considered in the context of CNO abundances and the multiple population scenario.

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.

In this paper we apply a Bayesian technique to determine the best fit of stellar evolution models to find the main sequence turn off age and other cluster parameters of four intermediate-age open clusters: NGC 2360, NGC 2477, NGC 2660, and NGC 3960. Our algorithm utilizes a Markov chain Monte Carlo technique to fit these various parameters, objectively finding the best-fit isochrone for each cluster. The result is a high-precision isochrone fit. We compare these results with the those of traditional "by-eye" isochrone fitting methods. By applying this Bayesian technique to NGC 2360, NGC 2477, NGC 2660, and NGC 3960, we determine the ages of these clusters to be 1.35 +/- 0.05, 1.02 +/- 0.02, 1.64 +/- 0.04, and 0.860 +/- 0.04 Gyr, respectively. The results of this paper continue our effort to determine cluster ages to higher precision than that offered by these traditional methods of isochrone fitting.

We use Cycle 21 Hubble Space Telescope (HST) observations and HST archival ACS Treasury observations of 30 Galactic Globular Clusters to characterize two distinct stellar populations. A sophisticated Bayesian technique is employed to simultaneously sample the joint posterior distribution of age, distance, and extinction for each cluster, as well as unique helium values for two populations within each cluster and the relative proportion of those populations. We find the helium differences among the two populations in the clusters fall in the range of ~0.04 to 0.11. Because adequate models varying in CNO are not presently available, we view these spreads as upper limits and present them with statistical rather than observational uncertainties. Evidence supports previous studies suggesting an increase in helium content concurrent with increasing mass of the cluster and also find that the proportion of the first population of stars increases with mass as well. Our results are examined in the context of proposed globular cluster formation scenarios. Additionally, we leverage our Bayesian technique to shed light on inconsistencies between the theoretical models and the observed data.

We present optical spectroscopy and near-infrared photometry of 57 faint ($g= 19-22$) high proper motion white dwarfs identified through repeat imaging of $\approx3100$ square degrees of the Sloan Digital Sky Survey footprint by Munn et al. (2014). We use $ugriz$ and $JH$ photometry to perform a model atmosphere analysis, and identify ten ultracool white dwarfs with $T_{\rm eff}<4000$ K, including the coolest pure H atmosphere white dwarf currently known, J1657+2638, with $T_{\rm eff}= 3550 \pm 100$K. The majority of the objects with cooling ages larger than 9 Gyr display thick disc kinematics and constrain the age of the thick disc to $\geq11$ Gyr. There are four white dwarfs in our sample with large tangential velocities ($v_{tan} > 120$ km $s^{-1}$) and UVW velocities that are more consistent with the halo than the Galactic disc. For typical $0.6 M_{\odot}$ white dwarfs, the cooling ages for these halo candidates range from 2.3 to 8.5 Gyr. However, the total main-sequence + white dwarf cooling ages of these stars would be consistent with the Galactic halo if they are slightly undermassive. Given the magnitude limits of the current large scale surveys, many of the coolest and oldest white dwarfs remain undiscovered in the solar neighborhood, but upcoming surveys such as GAIA and the Large Synoptic Survey Telescope (LSST) should find many of these elusive thick disc and halo white dwarfs.

It takes years of effort employing the best telescopes and instruments to obtain high-quality stellar photometry, astrometry, and spectroscopy. Stellar evolution models contain the experience of lifetimes of theoretical calculations and testing. Yet most astronomers fit these valuable models to these precious datasets by eye. We show that a principled Bayesian approach to fitting models to stellar data yields substantially more information over a range of stellar astrophysics. We highlight advances in determining the ages of star clusters, mass ratios of binary stars, limitations in the accuracy of stellar models, post-main-sequence mass loss, and the ages of individual white dwarfs. We also outline a number of unsolved problems that would benefit from principled Bayesian analyses.

We develop a Bayesian model for globular clusters composed of multiple stellar populations, extending earlier statistical models for open clusters composed of simple (single) stellar populations (vanDyk et al. 2009, Stein et al. 2013). Specifically, we model globular clusters with two populations that differ in helium abundance. Our model assumes a hierarchical structuring of the parameters in which physical properties---age, metallicity, helium abundance, distance, absorption, and initial mass---are common to (i) the cluster as a whole or to (ii) individual populations within a cluster, or are unique to (iii) individual stars. An adaptive Markov chain Monte Carlo (MCMC) algorithm is devised for model fitting that greatly improves convergence relative to its precursor non-adaptive MCMC algorithm. Our model and computational tools are incorporated into an open-source software suite known as BASE-9. We use numerical studies to demonstrate that our method can recover parameters of two-population clusters, and also show model misspecification can potentially be identified. As a proof of concept, we analyze the two stellar populations of globular cluster NGC 5272 using our model and methods. (BASE-9 is available from GitHub: https://github.com/argiopetech/base/releases).

We use Cycle 21 Hubble Space Telescope (HST) observations and HST archival ACS Treasury observations of Galactic Globular Clusters to find and characterize two stellar populations in NGC 5024 (M53), NGC 5272 (M3), and NGC 6352. For these three clusters, both single and double-population analyses are used to determine a best fit isochrone(s). We employ a sophisticated Bayesian analysis technique to simultaneously fit the cluster parameters (age, distance, absorption, and metallicity) that characterize each cluster. For the two-population analysis, unique population level helium values are also fit to each distinct population of the cluster and the relative proportions of the populations are determined. We find differences in helium ranging from $\sim$0.05 to 0.11 for these three clusters. Model grids with solar $\alpha$-element abundances ([$\alpha$/Fe] =0.0) and enhanced $\alpha$-elements ([$\alpha$/Fe]=0.4) are adopted.

Over the last decade proposal success rates in the fundamental sciences have dropped significantly. Astronomy and related fields funded by NASA and NSF are no exception. Data across agencies show that this is not principally the result of a decline in proposal merit (the proportion of proposals receiving high rankings is largely unchanged), nor of a shift in proposer demographics (seniority, gender, and institutional affiliation have all remained unchanged), nor of an increase (beyond inflation) in the average requested funding per proposal, nor of an increase in the number of proposals per investigator in any one year. Rather, the statistics are consistent with a scenario in which agency budgets for competed research are flat or decreasing in inflation-adjusted dollars, the overall population of investigators has grown, and a larger proportion of these investigators are resubmitting meritorious but unfunded proposals. This White Paper presents statistics which support this conclusion, as well as recent research on the time cost of proposal writing versus that of producing publishable results. We conclude that an aspirational proposal success rate of 30-35% would still provide a healthily competitive environment for researchers, would more fully utilize the scientific capacity of the community's facilities and missions, and provide relief to the funding agencies who face the logistics of ever-increasing volumes of proposals.

We surveyed 113 astronomers and 82 psychologists active in applying for federally funded research on their grant-writing history between January, 2009 and November, 2012. We collected demographic data, effort levels, success rates, and perceived non-financial benefits from writing grant proposals. We find that the average proposal takes 116 PI hours and 55 CI hours to write; although time spent writing was not related to whether the grant was funded. Effort did translate into success, however, as academics who wrote more grants received more funding. Participants indicated modest non-monetary benefits from grant writing, with psychologists reporting a somewhat greater benefit overall than astronomers. These perceptions of non-financial benefits were unrelated to how many grants investigators applied for, the number of grants they received, or the amount of time they devoted to writing their proposals. We also explored the number of years an investigator can afford to apply unsuccessfully for research grants and our analyses suggest that funding rates below approximately 20%, commensurate with current NIH and NSF funding, are likely to drive at least half of the active researchers away from federally funded research. We conclude with recommendations and suggestions for individual investigators and for department heads.

This paper provides a detailed comparison of the differences in parameters derived for a star cluster from its color-magnitude diagrams depending on the filters and models used. We examine the consistency and reliability of fitting three widely-used stellar evolution models to fifteen combinations of optical and near-IR photometry for the old open cluster NGC 188. The optical filter response curves match those of the theoretical systems and are thus not the source of fit inconsistencies. NGC 188 is ideally suited to the present study thanks to a wide variety of high-quality photometry and available proper motions and radial velocities which enable us to remove non-cluster members and many binaries. Our Bayesian fitting technique yields inferred values of age, metallicity, distance modulus, and absorption as a function of the photometric band combinations and stellar models. We show that the historically-favored three band combinations of UBV and VRI can be meaningfully inconsistent with each other and with longer baseline datasets such as UBVRIJHKs. Differences among model sets can also be substantial. For instance, fitting Yi et al. (2001) and Dotter et al. (2008) models to UBVRIJHKs photometry for NGC 188 yields the following cluster parameters: age={5.78+ 0.03, 6.45+-0.04} Gyr, [Fe/H]={+0.125+-0.003, -0.077+-0.003} dex, m-M={11.441+-0.007, 11.525+-0.005} mag, and Av={0.162+-0.003, 0.236+-0.003} mag, respectively. Within the formal fitting errors, these two fits are substantially and statistically different. Such differences amongst fits using different filters and models are a cautionary tale regarding our current ability to fit star cluster color-magnitude diagrams. Additional modeling of this kind, with more models and star clusters, and future GAIA parallaxes are critical for isolating and quantifying the most relevant uncertainties in stellar evolutionary models.

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.

A new proper motion catalog is presented, combining the Sloan Digital Sky Survey (SDSS) with second epoch observations in the r band within a portion of the SDSS imaging footprint. The new observations were obtained with the 90prime camera on the Steward Observatory Bok 90 inch telescope, and the Array Camera on the U.S. Naval Observatory, Flagstaff Station, 1.3 meter telescope. The catalog covers 1098 square degrees to r = 22.0, an additional 1521 square degrees to r = 20.9, plus a further 488 square degrees of lesser quality data. Statistical errors in the proper motions range from 5 mas/year at the bright end to 15 mas/year at the faint end, for a typical epoch difference of 6 years. Systematic errors are estimated to be roughly 1 mas/year for the Array Camera data, and as much as 2 - 4 mas/year for the 90prime data (though typically less). The catalog also includes a second epoch of r band photometry.

We apply a self-consistent and robust Bayesian statistical approach to determining the ages, distances, and ZAMS masses of 28 field DA white dwarfs with ages of approximately 4 to 8 Gyrs. Our technique requires only quality optical and near-IR photometry to derive ages with < 15% uncertainties, generally with little sensitivity to our choice of modern initial-final mass relation. We find that age, distance, and ZAMS mass are correlated in a manner that is too complex to be captured by traditional error propagation techniques. We further find that the posterior distributions of age are often asymmetric, indicating that the standard approach to deriving WD ages can yield misleading results.

We identify proper motion objects in the Hubble Ultra Deep Field (UDF) using the optical data from the original UDF program in 2004 and the near-infrared data from the 128-orbit UDF 2012 campaign. There are 12 sources brighter than I=27 mag that display >3sigma significant proper motions. We do not find any proper motion objects fainter than this magnitude limit. Combining optical and near-infrared photometry, we model the spectral energy distribution of each point-source using stellar templates and state-of-the-art white dwarf models. For I<27 mag, we identify 23 stars with K0-M6 spectral types and two faint blue objects that are clearly old, thick disk white dwarfs. We measure a thick disk white dwarf space density of 0.1-1.7 E-3 per cubic parsec from these two objects. There are no halo white dwarfs in the UDF down to I=27 mag. Combining the Hubble Deep Field North, South, and the UDF data, we do not see any evidence for dark matter in the form of faint halo white dwarfs, and the observed population of white dwarfs can be explained with the standard Galactic models.

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

We 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 constrain the distribution of calcium across the surface of the white dwarf star G29-38 by combining time series spectroscopy from Gemini-North with global time series photometry from the Whole Earth Telescope. G29-38 is actively accreting metals from a known debris disk. Since the metals sink significantly faster than they mix across the surface, any inhomogeneity in the accretion process will appear as an inhomogeneity of the metals on the surface of the star. We measure the flux amplitudes and the calcium equivalent width amplitudes for two large pulsations excited on G29-38 in 2008. The ratio of these amplitudes best fits a model for polar accretion of calcium and rules out equatorial accretion.

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 model the mineralogy and distribution of dust around the white dwarf G29-39 using the infrared spectrum from 1-35 microns and combining a wide range of materials based on spectral studies of comets and debris disks. In order of their contribution to the mid-infrared emission, the most abundant minerals for G29-38 are amorphous carbon, amorphous and crystalline silicates, water ice, and metal sulfides. The amorphous C can be equivalently replaced by other materials (like metallic Fe) with featureless infrared spectra. The best-fitting crystalline silicate is Fe-rich pyroxene. In order to absorb enough starlight to power the observed emission, the disk must either be much thinner than the stellar radius (so that it can be heated from above and below) or it must have an opening angle wider than 2 degrees. A `moderately optically thick' torus model with mass 2x10^19 g fits the spectrum well. A physically thin (less than the white dwarf radius) and optically thick disk can contribute to the near-infrared continuum only; such a disk cannot explain the longer-wavelength continuum or strong silicate features. The silicate composition contains minerals found from cometary spectra and meteorites, but Fe-rich pyroxene is more abundant than enstatite (Mg-rich pyroxene) or forsterite (Mg-rich olivine) in G29-38, in contrast to what is found in most comet or meteorite mineralogies. Enstatite meteorites may be the most similar solar system materials to the G29-38 dust. Finally, we suggest the surviving core of a `hot jupiter' as an alternative (neither cometary nor asteroidal) origin for the debris, though further theoretical work is needed to determine if this hypothesis is viable.

Due to the short settling times of metals in DA white dwarf atmospheres, any white dwarfs with photospheric metals must be actively accreting. It is therefore natural to expect that the metals may not be deposited uniformly on the surface of the star. We present calculations showing how the temperature variations associated with white dwarf pulsations lead to an observable diagnostic of the surface metal distribution, and we show what constraints current data sets are able to provide. We also investigate the effect that time-variable accretion has on the metal abundances of different species, and we show how this can lead to constraints on the gravitational settling times.

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

We present the result of a wide-field survey of globular clusters (GCs) in M31 covering a 3deg x 3deg field c. We have searched for GCs on CCD images taken with Washington CMT1 filters at the KPNO 0.9 m telescope using steps: (1) inspection of morphological parameters given by the SExtractor package such as stellarity, full maximum, and ellipticity; (2) consulting the spectral types and radial velocities obtained from spectra takena spectrograph at the WIYN 3.5 m telescope; and (3) visual inspection of the images of each object. We have and GC candidates, of which 605 are newly found GCs and GC candidates and 559 are previously known GCs. Amoects there are 113 genuine GCs, 258 probable GCs, and 234 possible GCs, according to our classification critee known objects there are 383 genuine GCs, 109 probable GCs, and 67 possible GCs. In total there are 496 genprobable GCs and 301 possible GCs. Most of these newly found GCs have T1 magnitudes of 17.5 - 19.5 mag, [17.9 < V < 19.9 mag assuming (C-T1) ~ 1.5], and (C-T1) colors in the range 1 - 2.

Our mid-infrared survey of 124 white dwarfs with the Spitzer Space Telescope and the IRAC imager has revealed an infrared excess associated with the white dwarf WD 2115-560 naturally explained by circumstellar dust. This object is the fourth white dwarf observed to have circumstellar dust. All four are DAZ white dwarfs, i.e. they have both photospheric Balmer lines and photospheric metal lines. We discuss these four objects as a class, which we abbreviate "DAZd", where the "d" stands for "dust". Using an optically-thick, geometrically-thin disk model analogous to Saturn's rings, we find that the inner disk edges are at >~0.1 to 0.2 Ro and that the outer disk edges are ~0.3 to 0.6 Ro. This model naturally explains the accretion rates and lifetimes of the detected WD disks and the accretion rates inferred from photospheric metal abundances.

Metals in the photospheres of white dwarfs with Teff between 12,000 and 25,000 K should gravitationally settle out of these atmospheres in 1-2 weeks. Temporal variations in the line strengths of these metals could provide a direct measurement of episodic metal accretion. Using archival VLT and Keck spectroscopy, we find evidence that the DAZd white dwarf G29-38 shows significant changes in its Ca II K line strength. At the two best-observed epochs, we find that the Ca line equivalent width (EW) = 165 +- 4 mA (in 1996.885) and 280 +- 8 mA (in 1999.653), which is an increase of 70%. We consider the effect that pulsation has on the Ca EWs for this known variable star, and find that it adds an error of < 1 mA to these measurements. Calcium line strengths at other observational epochs support variations with timescales as short as two weeks. These Ca EW variations indicate that the metal accretion process in G29-38, presumably from its debris disk, is episodic on timescales of a few weeks or less, and thus the accretion is not dominated by Poynting-Robertson drag from an optically thin, continuous disk, which has a timescale of ~1 year.

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 present mid-infrared photometry of 124 white dwarf stars with Spitzer Space Telescope. Objects were observed simultaneously at 4.5 and 8.0um with sensitivities better than 1 mJy. This data can be used to test models of white dwarf atmospheres in a new wavelength regime, as well as to search for planetary companions and debris disks.

We report on the discovery of a faint (M_V ~ -10.6 +/- 0.2) dwarf spheroidal galaxy on deep F606W and F814W Hubble Space Telescope images of a Virgo intracluster field. The galaxy is easily resolved in our images, as our color magnitude diagram (CMD) extends > 1 magnitude beyond the tip of the red giant branch (RGB). Thus, it is the deepest CMD for a small dwarf galaxy inside a cluster environment. Using the colors of the RGB stars, we derive a metal abundance for the dwarf of [M/H]= -2.3 +/- 0.3, and show that the metallicity dispersion is less than 0.6 dex at 95% confidence. We also use the galaxy's lack of AGB stars and the absence of objects brighter than M_bol ~ -4.1 +/- 0.2 to show that the system is old (t >~10 Gyr). Finally, we derive the object's structural parameters, and show that the galaxy displays no obvious evidence of tidal threshing. Since the tip of the red giant branch distance ((m-M)_0 = 31.23 +/- 0.17 or D = 17.6 +/- 1.4 Mpc) puts the galaxy near the core of the Virgo cluster, one might expect the object to have undergone some tidal processing. Yet the chemical and morphological similarity between the dwarf and the dSph galaxies of the Local and M81 Group demonstrates that the object is indeed pristine, and not the shredded remains of a much larger galaxy. We discuss the possible origins of this galaxy, and suggest that it is just now falling into Virgo for the first time.

We present preliminary limits on the presence of planets around white dwarf stars using the IRAC photometer on the Spitzer space telescope. Planets emit strongly in the mid-infrared which allows their presence to be detected as an excess at these wavelengths. We place limits of $5 M_J$ for 8 stars assuming ages of $1 Gyr$, and $10 M_J$ for 23 stars.We describe our survey, present our results and comment on approaches to improve our methodology.

We have used the Hubble Space Telescope's Advanced Camera for Surveys (ACS) to detect and measure ~5300 stars in a single intracluster field in the Virgo Cluster. By performing F606W and F814W photometry on these stars, we have determined their metallicity distribution function, and constrained the types of stars present in this portion of Virgo's intracluster space. Based on the small number of stars detected brighter than the red giant branch (RGB) tip, we suggest that in this region, Virgo's intracluster stars are mostly old (>~10 Gyr). Through analysis of the RGB stars themselves, we determine that the population contains the full range of metallicities probed (-2.3<[M/H]<0.0). We also present evidence that the younger (<10 Gyr) component of the population is more metal-rich, with [M/H]>-0.5. The spatial distribution of the most metal-poor stars in the field shows significantly more structure than that of the metal-rich stars, indicating that the intracluster population is not well-mixed. We discuss the implications these observations have for the production of intracluster stars and the dynamical evolution of the Virgo Cluster.

We report the discovery of 4 candidate intracluster globular clusters (IGCs) in a single deep HST ACS field of the Virgo Cluster. We show that each cluster is roughly spherical, has a magnitude near the peak of the Virgo globular cluster luminosity function, has a radial profile that is best-fit by a King model, and is surrounded by an excess of point sources which have the colors and magnitudes of cluster red giant stars. Despite the fact that two of our IGC candidates have integrated colors redder than the mean of the M87 globular cluster system, we propose that all of the objects are metal-poor with [M/H] < -1. We show that the tidal radii of our intracluster globulars are all larger than the mean for Milky Way clusters, and suggest that the clusters have undergone less tidal stress than their Galactic counterparts. Finally, we normalize our globular cluster observations to the luminosity of intracluster stars, and derive a value of S_N ~ 6 for the specific frequency of Virgo intracluster globular clusters. We use these data to constrain the origins of Virgo's intracluster population, and suggest that globular clusters in our intracluster field have a different origin than globular clusters in the vicinity of M87. In particular, we argue that dwarf elliptical galaxies may be an important source of intracluster stars.

We present near-infrared spectroscopic observations of 20 previously known DAZ white dwarfs obtained at the NASA Infrared Telescope Facility. Two of these white dwarfs (G29-38 and GD362) are known to display significant K-band excesses due to circumstellar debris disks. Here we report the discovery of excess K-band radiation from another DAZ white dwarf, WD0408-041 (GD56). Using spectroscopic observations, we show that the excess radiation cannot be explained by a stellar or substellar companion, and is likely to be caused by a warm debris disk. Our observations strengthen the connection between the debris disk phenomena and the observed metal abundances in cool DAZ white dwarfs. However, we do not find any excess infrared emission from the most metal rich DAZs with $T_{\rm eff}=$ 16000 -- 20000 K. This suggests that the metal abundances in warmer DAZ white dwarfs may require another explanation.

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.

We present 4.5$\mu$m and 8$\mu$m photometric observations of 18 cool white dwarfs obtained with the Spitzer Space Telescope. Our observations demonstrate that four white dwarfs with T_eff< 6000 K show slightly depressed mid-infrared fluxes relative to white dwarf models. In addition, another white dwarf with a peculiar optical and near-infrared spectral energy distribution (LHS 1126) is found to display significant flux deficits in Spitzer observations. These mid-infrared flux deficits are not predicted by the current white dwarf models including collision induced absorption due to molecular hydrogen. We postulate that either the collision induced absorption calculations are incomplete or there are other unrecognized physical processes occuring in cool white dwarf atmospheres. The spectral energy distribution of LHS 1126 surprisingly fits a Rayleigh-Jeans spectrum in the infrared, mimicking a hot white dwarf with effective temperature well in excess of 10$^5$ K. This implies that the source of this flux deficit is probably not molecular absorption but some other process.

We present new observations of the white dwarf G 29-38 with the camera (4.5 and 8 microns), photometer (24 microns), and spectrograph (5.5-14 microns) of the Spitzer Space Telescope. This star has an exceptionally large infrared excess amounting to 3% of the bolometric luminosity. The spectral energy distribution has a continuum peak around 4.5 micros and a 9-11 micron emission feature 1.25 times brighter than the continuum. A mixture of amorphous olivine and a small amount of forsterite in an emitting region 1-5 Rsun from the star can reproduce the shape of the 9-11 micron feature. The spectral energy distribution also appears to require amorphous carbon to explain the hot continuum. Our new measurements support the idea that a relatively recent disruption of a comet or asteroid created the cloud.

A sample of white dwarfs is selected from SDSS DR3 imaging data using their reduced proper motions, based on improved proper motions from SDSS plus USNO-B combined data. Numerous SDSS and followup spectra (Kilic et al. 2005) are used to quantify completeness and contamination of the sample; kinematic models are used to understand and correct for velocity-dependent selection biases. A luminosity function is constructed covering the range 7 < M_bol < 16, and its sensitivity to various assumptions and selection limits is discussed. The white dwarf luminosity function based on 6000 stars is remarkably smooth, and rises nearly monotonically to M_bol = 15.3. It then drops abruptly, although the small number of low-luminosity stars in the sample and their unknown atmospheric composition prevent quantitative conclusions about this decline. Stars are identified that may have high tangential velocities, and a preliminary luminosity function is constructed for them.

We report the discovery of excess K-band radiation from a massive DAZ white dwarf star, GD362. Combining infrared photometric and spectroscopic observations, we show that the excess radiation cannot be explained by a stellar or substellar companion, and is likely to be caused by a debris disk. This would be only the second such system known, discovered 18 years after G29-38, the only single white dwarf currently known to be orbited by circumstellar dust. Both of these systems favor a model with accretion from a surrounding debris disk to explain the metal abundances observed in DAZ white dwarfs. Nevertheless, observations of more DAZs in the mid-infrared are required to test if this model can explain all DAZs.

The study of open clusters has a classic feel to it since the subject predates anyone alive today. Despite the age of this topic, I show via an ADS search that its relevance and importance in astronomy has grown faster in the last few decades than astronomy in general. This is surely due to both technical reasons and the interconnection of the field of stellar evolution to many branches of astronomy. In this review, I outline what we know today about open clusters and what they have taught us about a range of topics from stellar evolution to Galactic structure to stellar disk dissipation timescales. I argue that the most important astrophysics we have learned from open clusters is stellar evolution and that its most important product has been reasonably precise stellar ages. I discuss where open cluster research is likely to go in the next few years, as well as in the era of 20m telescopes, SIM, and GAIA. Age will continue to be of wide relevance in astronomy, from cosmology to planet formation timescales, and with distance errors soon no longer a problem, improved ages will be critically important to many of the most fascinating astrophysical questions.

We explore the nature of the faint blue objects in the Hubble Deep Field South. We have derived proper motions for the point sources in the Hubble Deep Field South using a 3 year baseline. Combining our proper motion measurements with spectral energy distribution fitting enabled us to identify 4 quasars and 42 stars, including 3 white dwarf candidates. Two of these white dwarf candidates, HDFS 1444 and 895, are found to display significant proper motion, 21.1 $\pm$ 7.9 mas/yr and 34.9 $\pm$ 8.0 mas/yr, and are consistent with being thick disk or halo white dwarfs located at ~2 kpc. The other faint blue objects analyzed by Mendez & Minniti do not show any significant proper motion and are inconsistent with being halo white dwarfs; they do not contribute to the Galactic dark matter. The observed population of stars and white dwarfs is consistent with standard Galactic models.

We report our measurement of the rate of change of period with time dP/dt for the 215 s periodicity in the pulsating white dwarf G 117-B15A, the most stable optical clock known. After 31 years of observations, we have finally obtained a 4 sigma measurement dP/dt_observed = (4.27 +/- 0.80) x 10^{-15} s/s. Taking into account the proper-motion effect of dP/dt_pm = (7.0 +/- 2.0) x 10^{-16} s/s, we obtain a rate of change of period with time of dP/dt = (3.57 +/- 0.82) x 10^{-15} s/s. This value is consistent with the cooling rate in our white dwarf models only for cores of C or C/O. With the refinement of the models, the observed rate of period change can be used to accurately measure the ratio of C/O in the core of the white dwarf.

We present deep VI photometry of stars in the globular cluster M5 (NGC 5904) based on images taken with the Hubble Space Telescope. The resulting color-magnitude diagram reaches below V ~ 27 mag, revealing the upper 2-3 magnitudes of the white dwarf cooling sequence, and main sequence stars eight magnitudes and more below the turn-off. We fit the main sequence to subdwarfs of known parallax to obtain a true distance modulus of (m-M)_0 = 14.45 +/- 0.11 mag. A second distance estimate based on fitting the cluster white dwarf sequence to field white dwarfs with known parallax yielded (m-M)_0 = 14.67 +/- 0.18 mag. We couple our distance estimates with extensive photometry of the cluster's RR Lyrae variables to provide a calibration of the RR Lyrae absolute magnitude yielding M_V(RR) = 0.42 +/- 0.10 mag at [Fe/H] = -1.11 dex. We provide another luminosity calibration in the form of reddening-free Wasenheit functions. Comparison of our calibrations with predictions based on recent models combining stellar evolution and pulsation theories shows encouraging agreement. (Abridged)

A reduced proper motion diagram utilizing Sloan Digital Sky Survey (SDSS) photometry and astrometry and USNO-B plate astrometry is used to separate cool white dwarf candidates from metal-weak, high-velocity main sequence Population II stars (subdwarfs) in the SDSS Data Release 2 imaging area. Follow-up spectroscopy using the Hobby-Eberly Telescope, the MMT, and the McDonald 2.7m Telescope is used to demonstrate that the white dwarf and subdwarf loci separate cleanly in the reduced proper motion diagram, and that the contamination by subdwarfs is small near the cool white dwarf locus. This enables large statistically complete samples of white dwarfs, particularly the poorly understood cool white dwarfs, to be created from the SDSS imaging survey, with important implications for white dwarf luminosity function studies. SDSS photometry for our sample of cool white dwarfs is compared to current white dwarf models.

I explore the current ability of both white dwarf cooling theory and main sequence stellar evolution theory to accurately determine stellar population ages by comparing ages derived using both techniques for open clusters ranging from 0.1 to 4 Gyr. I find good agreement between white dwarf and main sequence evolutionary ages over the entire age range currently available for study. I also find that directly comparing main sequence turn-off ages to white dwarf ages is only weakly sensitive to realistic levels of errors in cluster distance, metallicity, and reddening. Additional detailed comparisons between white dwarf and main sequence ages have tremendous potential to refine and calibrate both of these important clocks, and I present new simulations of promising open cluster targets. The most demanding requirement for these white dwarf studies are very deep (V > 25-28) cluster observations made necessary by the faintness of the oldest white dwarfs.

We present new observations of very faint white dwarfs (WDs) in the rich open star cluster NGC 2099 (M37). Following deep, wide field imaging of the cluster using CFHT, we have now obtained spectroscopic observations of candidate WDs using both GMOS on Gemini and LRIS on Keck. Of our 24 WD candidates (all fainter than V = 22.4), 21 are spectroscopically confirmed to be bona fide WDs, 4-5 of which are most likely field objects. Fitting 18 of the 21 WD spectra with model atmospheres, we find that most WDs in this cluster are quite massive (0.7-0.9 Msun), as expected given the cluster's young age (650 Myr) and, hence, high turnoff mass (~2.4 Msun). We determine a new initial-final mass relationship and almost double the number of existing data points from previous studies. The results indicate that stars with initial masses between 2.8 and 3.4 Msun lose 70-75% of their mass through stellar evolution. For the first time, we find some evidence of a metallicity dependence on the initial-final mass relationship.

We present an update of our on-going effort to improve the precision of white dwarf cosmochronology via careful analyses of white dwarf photometry in open clusters. To improve the precision of white dwarf and main sequence age analysis, we are developing a new interpretative scheme using a Bayesian statistical approach that matches observations to simulated stellar clusters. Here we present our first tests of the Bayesian approach with simulated stellar clusters with ages of 1, 2, and 4 billion years.

We describe a new survey for cool white dwarfs that supplements Sloan Digital Sky Survey photometry with USNO proper motions and follow-up spectroscopy. To date we have discovered and spectroscopically confirmed 80 new moderate temperature and cool white dwarfs. We have also found a handful of high-velocity white dwarfs and we expect a sizable fraction of these to be thick disk or possibly halo objects. Our survey is designed to find ~10,000 new white dwarfs, although only ~60 will be among the faintest white dwarfs (M(V) >= 16), where most of the age-sensitivity resides. We discuss an extension of our survey to V ~ 22.

We investigate the use of a narrow-band DDO51 filter for photometric identification of cool white dwarfs. We report photometric observations of 30 known cool white dwarfs with temperatures ranging from 10,000 K down to very cool temperatures (<3500 K). Follow-up spectroscopic observations of a sample of objects selected using this filter and our photometric observations show that DDO51 filter photometry can help select cool white dwarf candidates for follow-up multi--object spectroscopy by rejecting 65% of main sequence stars with the same broad--band colors as the cool white dwarfs. This technique is not selective enough to efficiently feed single--object spectrographs. We present the white dwarf cooling sequence using this filter. Our observations show that very cool white dwarfs form a sequence in the r-DDO vs. r-z color--color diagram and demonstrate that significant improvements are needed in white dwarf model atmospheres.

We use the new ZZ Ceti stars (hydrogen atmosphere white dwarf variables; DAVs) discovered within the Sloan Digital Sky Survey (Mukadam et al. 2004) to re-define the empirical ZZ Ceti instability strip. This is the first time since the discovery of white dwarf variables in 1968 that we have a homogeneous set of spectra acquired using the same instrument on the same telescope, and with consistent data reductions, for a statistically significant sample of ZZ Ceti stars. The homogeneity of the spectra reduces the scatter in the spectroscopic temperatures and we find a narrow instability strip of width ~950K, from 10850--11800K. We question the purity of the DAV instability strip as we find several non-variables within. We present our best fit for the red edge and our constraint for the blue edge of the instability strip, determined using a statistical approach.

Using the deepest and finest resolution images of the Universe acquired with the Hubble Space Telescope and a similar image taken 7 years later for the Great Observatories Origins Deep Survey, we have derived proper motions for the point sources in the Hubble Deep Field--North. Two faint blue objects,HDF2234 and HDF3072, are found to display significant proper motion, 10.0 $\pm$ 2.5 and 15.5 $\pm$ 3.8 mas yr$^{-1}$. Photometric distances and tangential velocities for these stars are consistent with disk white dwarfs located at $\sim$ 500 pc. The faint blue objects analyzed by Ibata et al. (1999) and Mendez & Minniti (2000) do not show any significant proper motion; they are not halo white dwarfs and they do not contribute to the Galactic dark matter. These objects are likely to be distant AGN.

We confirm and extend our previous detection of a population of intergalactic globular clusters in Abell 1185, and report the first discovery of an intergalactic globular cluster in the nearby Virgo cluster of galaxies. The numbers, colors and luminosities of these objects can place constraints on their origin, which in turn may yield new insights to the evolution of galaxies in dense environments.

We obtained new VLT/ISAAC H-band observations for five MACHO LMC source stars and adjacent LMC field regions. After combining our near-IR photometry with HST/PC BVRI optical photometry, we compared the MACHO objects to the adjacent field stars in a variety of color-magnitude and color-color diagrams. These diagnostic diagrams were chosen to be sensitive to our hypothesis that at least some of the MACHO lenses were foreground Galactic disk or thick disk M dwarfs. For the five lensed objects we studied, our hypothesis could be ruled out for main sequence lens masses >= 0.1 Mo for distances out to 4 kpc. On the other hand, the fact that LMC-MACHO-5, an object not in our study, has been recently found to have just such a foreground lens, highlights that the remainder of the LMC MACHO objects should be searched for the signature of their lenses using our photometric technique, or via near-IR spectroscopy. We also constructed diagnostic color-color diagrams sensitive to determining reddening for the individual MACHO source stars and found that these five objects did not show evidence for significant additional reddening. At least these five MACHO objects are thus also inconsistent with the LMC self-lensing hypothesis.

We report our search for the cool white dwarfs belonging to the Galactic disk by extending the NOAO Deep Wide-Field Survey. Narrow-band DDO51 photometry of the Deep Wide-Field Survey's northern field was obtained using the 4m-Mayall Telescope and the MOSAIC imager to separate cool white dwarfs from other stellar types of similar T_eff. Follow-up spectroscopy of four white dwarf candidates from our photometric search resulted in the discovery of two new cool white dwarfs as companions to M dwarfs.

We report new infrared spectroscopic observations of cool DQ white dwarfs by using Coolspec on the 2.7m Harlan-Smith Telescope. DQs have helium-rich atmospheres with traces of molecular carbon thought to be the result of convective dredge-up from their C/O interiors. Recent model calculations predict that oxygen should also be present in DQ atmospheres in detectable amounts. Our synthetic spectra calculations for He-rich white dwarfs with traces of C and O indicate that CO should be easily detected in the cool DQ atmospheres if present in the expected amounts. Determination of the oxygen abundance in the atmosphere will reveal the C/O ratio at the core/envelope boundary, constraining the important and uncertain ^{12}C(alpha,gamma)^{16}O reaction rate.

Stellar spectroscopic classification has been successfully automated by a number of groups. Automated classification and parameterization work best when applied to a homogeneous data set, and thus these techniques primarily have been developed for and applied to large surveys. While most ongoing large spectroscopic surveys target extragalactic objects, many stellar spectra have been and will be obtained. We briefly summarize past work on automated classification and parameterization, with emphasis on the work done in our group. Accurate automated classification in the spectral type domain and parameterization in the temperature domain have been relatively easy. Automated parameterization in the metallicity domain, formally outside the MK system, has also been effective. Due to the subtle effects on the spectrum, automated classification in the luminosity domain has been somewhat more difficult, but still successful. In order to extend the use of automated techniques beyond a few surveys, we present our current efforts at building a web-based automated stellar spectroscopic classification and parameterization machine. Our proposed machinery would provide users with MK classifications as well as the astrophysical parameters of effective temperature, surface gravity, mean abundance, abundance anomalies, and microturbulence.

We present deep BVI observations of the core of M35 and a nearby comparison field obtained at the WIYN 3.5m telescope under excellent seeing. These observations display the lower main sequence in BV and VI CMDs down to V = 23.3 and 24.6, respectively. At these faint magnitudes background field stars are far more numerous than the cluster stars, yet by using a smoothing technique and CMD density distribution subtraction we recover the cluster fiducial main sequence and luminosity function to V = 24.6. We find the location of the main sequence in these CMDs to be consistent with earlier work on other open clusters, specifically NGC 188, NGC 2420, and NGC 2477. We compare these open cluster fiducial sequences to stellar models by Baraffe et al. (1998), Siess et al. (2000), Girardi et al. (2000), and Yi et al. (2001) and find that the models are too blue in both B-V and V-I for stars below ~0.4 Mo. M35 contains stars to the limit of the extracted main sequence, at M ~ 0.10-0.15 Mo, suggesting that M35 may harbor a large number of brown dwarfs, which should be easy targets for near-IR instrumentation on 8-10m telescopes. We also identify a new candidate white dwarf in M35 at V = 21.36 +- 0.01. Depending on which WD models are used to interpret this cluster candidate, it is either a very high mass WD (1.05 +- 0.05 Mo) somewhat older (0.19-0.26 Gyr, 3-4 sigma) than our best isochrone age (150 Myr), or it is a modestly massive WD (0.67-0.78 Mo) much too old (0.42-0.83 Gyr) to belong to the cluster.

We report medium-resolution (1-2 A) spectroscopy and broadband (UBV) photometry for a sample of 39 bright stars (the majority of which are likely to be giants) selected as metal-deficient candidates from an objective-prism survey concentrating on Galactic latitudes below |b| = 30 deg, the LSE survey of Drilling & Bergeron. Although the primary purpose of the LSE survey was to select OB stars (hence the concentration on low latitudes), the small number of bright metal-deficient giant candidates noted during this survey provide interesting information on the metal-weak thick disk (MWTD) population. The kinematics of the LSE giants indicate the presence of a rapidly rotating population, even at quite low metallicity. We consider the distribution of orbital eccentricity of the LSE giants as a function of [Fe/H], and conclude that the local fraction (i.e., within 1 kpc from the Sun) of metal-poor stars that might be associated with the MWTD is on the order of 30%-40% at abundances below [Fe/H] = -1.0. Contrary to recent analyses of previous (much larger) samples of non-kinematically selected metal-poor stars, we find that this relatively high fraction of local metal-poor stars associated with the MWTD may extend to metallicities below [Fe/H] = -1.6, much lower than had been considered before. We identify a subsample of 11 LSE stars that are very likely to be members of the MWTD, based on their derived kinematics; the lowest metallicity among these stars is [Fe/H] = -2.35. Implications of these results for the origin of the MWTD and for the formation of the Galaxy are considered. (abridged)

We explore the application of artificial neural networks (ANNs) for the estimation of atmospheric parameters (Teff, logg, and [Fe/H]) for Galactic F- and G-type stars. The ANNs are fed with medium-resolution (~ 1-2 A) non flux-calibrated spectroscopic observations. From a sample of 279 stars with previous high-resolution determinations of metallicity, and a set of (external) estimates of temperature and surface gravity, our ANNs are able to predict Teff with an accuracy of ~ 135-150 K over the range 4250 <= Teff <= 6500 K, logg with an accuracy of ~ 0.25-0.30 dex over the range 1.0 <= logg <= 5.0 dex, and [Fe/H] with an accuracy ~ 0.15-0.20 dex over the range -4.0 <= [Fe/H] <= +0.3. Such accuracies are competitive with the results obtained by fine analysis of high-resolution spectra. It is noteworthy that the ANNs are able to obtain these results without consideration of photometric information for these stars. We have also explored the impact of the signal-to-noise ratio (S/N) on the behavior of ANNs, and conclude that, when analyzed with ANNs trained on spectra of commensurate S/N, it is possible to extract physical parameter estimates of similar accuracy with stellar spectra having S/N as low as 13. Taken together, these results indicate that the ANN approach should be of primary importance for use in present and future large-scale spectroscopic surveys.

We have used deep HST WFPC2 observations of two fields in NGC 2420 to produce a cluster CMD down to V = 27. After imposing morphological selection criteria we find eight candidate white dwarfs in NGC 2420. Our completeness estimates indicate that we have found the terminus of the WD cooling sequence. We argue that the cluster distance modulus is likely to be close to 12.10 with E(B-V) = 0.04. With these parameters we find a white dwarf cooling age for NGC 2420 of 2.0 +- 0.20 (1 sigma) Gyrs. The 0.20 Gyr uncertainty includes errors in the photometry, sequence fitting, precursor time scales, and theoretical white dwarf cooling time scales. Comparing the cluster white dwarf cooling age to ages derived from stellar isochrone fitting we find a preference for ages derived from models incorporating convective overshoot.

New generation large-aperture telescopes, multi-object spectrographs, and large format detectors are making it possible to acquire very large samples of stellar spectra rapidly. In this context, traditional star-by-star spectroscopic analysis are no longer practical. New tools are required that are capable of extracting quickly and with reasonable accuracy important basic stellar parameters coded in the spectra. Recent analyses of Artificial Neural Networks (ANNs) applied to the classification of astronomical spectra have demonstrated the ability of this concept to derive estimates of temperature and luminosity. We have adapted the back-propagation ANN technique developed by von Hippel et al. (1994) to predict effective temperatures, gravities and overall metallicities from spectra with resolving power ~ 2000 and low signal-to-noise ratio. We show that ANN techniques are very effective in executing a three-parameter (Teff,log g,[Fe/H]) stellar classification. The preliminary results show that the technique is even capable of identifying outliers from the training sample.

We have employed precise V and I photometry of NGC 188 at WIYN to explore the cluster luminosity function (LF) and study the cluster white dwarfs (WDs). Our photometry is offset by V = 0.052 (fainter) from Sandage (1962) and Eggen & Sandage (1969). All published photometry for the past three decades have been based on these two calibrations, which are in error by 0.05 +- 0.01. We employ the Pinsonneault etal (1998) fiducial main sequence to derive a cluster distance modulus of 11.43 +- 0.08. We report observations that are >= 50% complete to V = 24.6 and find that the cluster central-field LF peaks at M_I ~ 3 to 4. This is unlike the solar neighborhood LF and unlike the LFs of dynamically unevolved portions of open and globular clusters, which rise continuously until M_I ~ 9.5. Although we find that >= 50% of the unresolved cluster objects are multiple systems, their presence cannot account for the shape of the NGC 188 LF. For theoretical reasons (Terlevich 1987; Vesperini & Heggie 1997) having to do with the survivability of NGC 188 we believe the cluster is highly dynamically evolved and that the missing low luminosity stars are either in the cluster outskirts or have left the cluster altogether. We identify nine candidate WDs, of which we expect three to six are bona fide cluster WDs. The luminosities of the faintest likely WD indicates an age (Bergeron, Wesemael, & Beauchamp 1995) of 1.14 +- 0.09 Gyrs. This is a lower limit to the cluster age and observations probing to V = 27 or 28 will be necessary to find the faintest cluster WDs and independently determine the cluster age. While our age limit is not surprising for this ~6 Gyr old cluster, our result demonstrates the value of the WD age technique with its very low internal errors. (abridged)

We describe a method for the extraction of spectra from high dispersion objective prism plates. Our method is a catalogue driven plate solution approach, making use of the Right Ascension and Declination coordinates for the target objects. In contrast to existing methods of photographic plate reduction, we digitize the entire plate and extract spectra off-line. This approach has the advantages that it can be applied to CCD objective prism images, and spectra can be re-extracted (or additional spectra extracted) without having to re-scan the plate. After a brief initial interactive period, the subsequent reduction procedure is completely automatic, resulting in fully-reduced, wavelength justified spectra. We also discuss a method of removing stellar continua using a combination of non-linear filtering algorithms. The method described is used to extract over 12,000 spectra from a set of 92 objective prism plates. These spectra are used in an associated project to develop automated spectral classifiers based on neural networks.

We investigate the application of neural networks to the automation of MK spectral classification. The data set for this project consists of a set of over 5000 optical (3800-5200 AA) spectra obtained from objective prism plates from the Michigan Spectral Survey. These spectra, along with their two-dimensional MK classifications listed in the Michigan Henry Draper Catalogue, were used to develop supervised neural network classifiers. We show that neural networks can give accurate spectral type classifications (sig_68 = 0.82 subtypes, sig_rms = 1.09 subtypes) across the full range of spectral types present in the data set (B2-M7). We show also that the networks yield correct luminosity classes for over 95% of both dwarfs and giants with a high degree of confidence. Stellar spectra generally contain a large amount of redundant information. We investigate the application of Principal Components Analysis (PCA) to the optimal compression of spectra. We show that PCA can compress the spectra by a factor of over 30 while retaining essentially all of the useful information in the data set. Furthermore, it is shown that this compression optimally removes noise and can be used to identify unusual spectra.

I present a literature search through 31 July 1997 of white dwarfs (WDs) in open and globular clusters. There are 36 single WDs and 5 WDs in binaries known among 13 open clusters, and 340 single WDs and 11 WDs in binaries known among 11 globular clusters. From these data I have calculated WD mass fractions for four open clusters (the Pleiades, NGC 2168, NGC 3532, and the Hyades) and one globular cluster (NGC 6121). I develop a simple model of cluster evolution that incorporates stellar evolution but not dynamical evolution to interpret the WD mass fractions. I augment the results of my simple model with N-body simulations incorporating stellar evolution (Terlevich 1987; de la Feunte Marcos 1996; Vesperini & Heggie 1997). I find that even though these clusters undergo moderate to strong kinematical evolution the WD mass fraction is relatively insensitive to kinematical evolution. By comparing the cluster mass functions to that of the Galactic disk, and incorporating plausibility arguments for the mass function of the Galactic halo, I estimate the WD mass fraction in these two populations. I assume the Galactic disk is ~10 Gyrs old (Winget et al. 1987; Liebert, Dahn, & Monet 1988; Oswalt et al. 1996) and that the Galactic halo is ~12 Gyrs old (Reid 1997b; Gratton et al. 1997; Chaboyer et al. 1998), although the WD mass fraction is insensitive to age in this range. I find that the Galactic halo should contain 8 to 9% (alpha = -2.35) or perhaps as much as 15 to 17% (alpha = -2.0) of its stellar mass in the form of WDs. The Galactic disk WD mass fraction should be 6 to 7% (alpha = -2.35), consistent with the empirical estimates of 3 to 7% (Liebert, Dahn, & Monet 1988; Oswalt et al. 1996). (abridged)

It has been suspected for nearly 50 years that clusters of galaxies contain a population of intergalactic stars, ripped from galaxies during cluster formation or when the galaxies' orbits take them through the cluster center. Support for the existence of such a population of free-floating stars comes from measurements of the diffuse light in clusters, and from recent detections of planetary nebulae with positions and/or velocities far removed from any observed cluster galaxy. But estimates for the mass of the diffuse population and its distribution relative to the galaxies are still highly uncertain. Here we report the direct detection of intergalactic stars in deep images of a blank field in the Virgo Cluster. The data suggest that approximately 10% of the stellar mass of the cluster is in intergalactic stars. We observe a relatively homogeneous distribution of stars, with evidence of a slight gradient toward M87.

We present a technique which employs artificial neural networks to produce physical parameters for stellar spectra. A neural network is trained on a set of synthetic optical stellar spectra to give physical parameters (e.g. T_eff, log g, [M/H]). The network is then used to produce physical parameters for real, observed spectra. Our neural networks are trained on a set of 155 synthetic spectra, generated using the SPECTRUM program written by Gray (Gray & Corbally 1994, Gray & Arlt 1996). Once trained, the neural network is used to yield T_eff for over 5000 B-K spectra extracted from a set of photographic objective prism plates (Bailer-Jones, Irwin & von Hippel 1997a). Using the MK classifications for these spectra assigned by Houk (1975, 1978, 1982, 1988) we have produced a temperature calibration of the MK system based on this set of 5000 spectra. It is demonstrated through the metallicity dependence of the derived temperature calibration that the neural networks are sensitive to the metallicity signature in the real spectra. With further work it is likely that neural networks will be able to yield reliable metallicity measurements for stellar spectra.

We discuss the luminosity function of SNe Ia under the assumption that recent evidence for dispersion in this standard candle is related to variations in the white dwarf mass function (WDMF) in the host galaxies. We develop a simple parameterization of the WDMF as a function of age of a stellar population and apply this to galaxies of different morphological types. We show that this simplified model is consistent with the observed WDMF of Bergeron et al. (1992) for the solar neighborhood. Our simple models predict that WDMF variations can produce a range of more than 1.8 mag in M$_B$(SN Ia), which is comparable to the observed value using the data of Phillips (1993) and van den Bergh (1996). We also predict a galaxy type dependence of M$_B$(SN Ia) under standard assumptions of the star formation history in these galaxies and show that M$_B$(SN Ia) can evolve with redshift. In principle both evolutionary and galaxy type corrections should be applied to recover the intrinsic range of M$_B$(SN Ia) from the observed values. Our current inadequate knowledge of the star formation history of galaxies coupled with poor physical understanding of the SN Ia mechanism makes the reliable estimation of these corrections both difficult and controversial. The predictions of our models combined with the observed galaxy and redshift correlations may have the power to discriminate between the Chandrasekhar and the sub-Chandrasekhar progenitor scenarios for SNe Ia.

Using HST and the WFPC2 we have acquired very deep V- and I-band photometry of stars in NGC 2420 and NGC 2477 to study cluster luminosity functions at approximately solar metallicity. We have determined these cluster luminosity functions down to $M_I$ = 10.5 (0.2 M$_{\odot}$) and find that the luminosity function of NGC 2420 turns over at $M_I$ $\approx$ 9.0, and possibly stops altogether by $M_I$ $\approx$ 9.5. The luminosity function of NGC 2477 may flatten at $M_I$ $\geq$ 9.5. We compare our open cluster luminosity functions to the solar neighborhood field star luminosity function of Kroupa, Tout \& Gilmore (1993) and the four published HST globular cluster luminosity functions: $\omega$ Cen (Elson {\it et al.}\ 1995), 47 Tuc (De Marchi \& Paresce 1995b), M 15 (De Marchi \& Paresce 1995a), and NGC 6397 (Paresce, De Marchi \& Romaniello 1995). We find a smooth relation between the location of the luminosity function turn-over and the metallicity for all these low mass star samples which matches the expected $M_I$ versus [Fe/H] trend for a model star of $\approx$ 0.27 M$_{\odot}$ (Saumon 1995; Alexander {\it et al.}\ 1996). We interpret this smooth and systematic behavior in the cluster luminosity functions as strong evidence in favor of an invariant initial mass function and a metallicity-dependent mass-luminosity relation.

Modern instrumentation makes it possible to measure the mass to radius ratio for main sequence stars in open clusters from gravitational redshifts. For stars where independent information is available for either the mass or the radius, this application of general relativity directly determines the other quantity. Applicable examples are: 1) measuring the radii of solar metallicity main sequence stars for which the mass - luminosity relation is well known, 2) measuring the radii for stars where model atmospheres can be used to determine the surface gravity (the mass to radius squared ratio), 3) refining the mass - radius relation for main sequence stars, and 4) measuring the change in radius as stars evolve off the main sequence and up the giant branch.

Since the WFPC-2 undersamples the PSF, aperture photometry can produce results which are competetive with profile fitting in many situations. This article reports and investigation of aperture corrections using both real data and PSF models.

The white dwarf luminosity function of a stellar cluster will have a sharp truncation at a luminosity which is determined by the time since formation of the first white dwarfs in that cluster. Calculation of the dependence of this limiting luminosity on age requires relatively well-understood physics and is independent of stellar evolutionary models. Thus, measurement of the termination of the white dwarf luminosity function provides an independent method to determine the age of a cluster, and thereby to calibrate stellar evolutionary ages. We have obtained HST WFPC2 data in two open clusters, identified the white dwarf sequence, and proved the feasibility of this approach, by detecting white dwarfs to V=25. Much deeper data are feasible. From our present limited data, we show that degenerate cooling ages are not consistent with some published isochrone ages for clusters with ages of order 1Gyr.