The Kilodegree Extremely Little Telescope (KELT) project has been conducting a photometric survey for transiting planets orbiting bright stars for over ten years. The KELT images have a pixel scale of ~23"/pixel---very similar to that of NASA's Transiting Exoplanet Survey Satellite (TESS)---as well as a large point spread function, and the KELT reduction pipeline uses a weighted photometric aperture with radius 3'. At this angular scale, multiple stars are typically blended in the photometric apertures. In order to identify false positives and confirm transiting exoplanets, we have assembled a follow-up network (KELT-FUN) to conduct imaging with higher spatial resolution, cadence, and photometric precision than the KELT telescopes, as well as spectroscopic observations of the candidate host stars. The KELT-FUN team has followed-up over 1,600 planet candidates since 2011, resulting in more than 20 planet discoveries. Excluding ~450 false alarms of non-astrophysical origin (i.e., instrumental noise or systematics), we present an all-sky catalog of the 1,128 bright stars (6<V<10) that show transit-like features in the KELT light curves, but which were subsequently determined to be astrophysical false positives (FPs) after photometric and/or spectroscopic follow-up observations. The KELT-FUN team continues to pursue KELT and other planet candidates and will eventually follow up certain classes of TESS candidates. The KELT FP catalog will help minimize the duplication of follow-up observations by current and future transit surveys such as TESS.

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.

Recent observations and a photometric search for variable stars in the Ursa Minor dwarf spheroidal galaxy (UMi dSph) are presented. Our observations were taken at Apache Point Observatory in 2014 and 2016 using the 0.5m ARCSAT telescope and the West Mountain Observatory 0.9m telescope of Brigham Young University in 2016. Previously known RR Lyrae stars in our field of view of the UMi dSph are identified, and we also catalog new variable star candidates. Tentative classifications are given for some of the new variable stars. We have conducted period searches with the data collected with the WMO telescope. Our ultimate goal is to create an updated catalog of variable stars in the UMi dSph and to compare the RR Lyrae stellar characteristics to other RR Lyrae stars found in the Local Group dSph galaxies.

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.

We have obtained 53 images with the $g$ filter and 19 images with the $i$ filter, each with 600-second exposures of the super metal rich open cluster NGC 6253 with the Gemini-South telescope to create deep images of the cluster to observe the cluster white dwarfs for the first time. We will analyze the white dwarf luminosity function to measure the cluster's white dwarf age, search for any anomalous features (as has been seen in the similarly metal rich cluster NGC 6791), and constrain the initial-final mass relation at high metallicities. We present an update on these observations and our program to study the formation of white dwarfs in super high metallicity environments.

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.

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.

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 present HST/ACS observations of RR Lyrae variable stars in six ultra deep fields of the Andromeda galaxy (M31), including parts of the halo, disk, and giant stellar stream. Past work on the RR Lyrae stars in M31 has focused on various aspects of the stellar populations that make up the galaxy's halo, including their distances and metallicities. This study builds upon this previous work by increasing the spatial coverage (something that has been lacking in previous studies) and by searching for these variable stars in constituents of the galaxy not yet explored. Besides the 55 RR Lyrae stars we found in our initial field located 11kpc from the galactic nucleus, we find additional RR Lyrae stars in four of the remaining five ultra deep fields as follows: 21 in the disk, 24 in the giant stellar stream, 3 in the halo field 21kpc from the galactic nucleus, and 5 in one of the halo fields at 35kpc. No RR Lyrae were found in the second halo field at 35kpc. The RR Lyrae populations of these fields appear to mostly be of Oosterhoff I type, although the 11kpc field appears to be intermediate or mixed. We will discuss the properties of these stars including period and reddening distributions. We calculate metallicities and distances for the stars in each of these fields using different methods and compare the results, to an extent that has not yet been done. We compare these methods not just on RR Lyrae stars in our M31 fields, but also on a data set of Milky Way field RR Lyrae stars.

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 new UV-to-IR stellar photometry of four low-extinction windows in the Galactic bulge, obtained with the Wide Field Camera 3 on the Hubble Space Telescope (HST). Using our five bandpasses, we have defined reddening-free photometric indices sensitive to stellar effective temperature and metallicity. We find that the bulge populations resemble those formed via classical dissipative collapse: each field is dominated by an old (~10 Gyr) population exhibiting a wide metallicity range (-1.5 < [Fe/H] < 0.5). We detect a metallicity gradient in the bulge population, with the fraction of stars at super-solar metallicities dropping from 41% to 35% over distances from the Galactic center ranging from 0.3 to 1.2 kpc. One field includes candidate exoplanet hosts discovered in the SWEEPS HST transit survey. Our measurements for 11 of these hosts demonstrate that exoplanets in the distinct bulge environment are preferentially found around high-metallicity stars, as in the solar neighborhood, supporting the view that planets form more readily in metal-rich environments.

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

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.

We report the discovery of several optical burst-like events from the low-mass X-ray binary MS1603.6+2600 (UW CrB). The events last for a few tens of seconds, exhibit a very fast rise and slow decay, and involve optical brightening of a factor of 2-3. The flares appear distinct from the lower level flickering and instead strongly resemble reprocessed type-I X-ray bursts as seen in a number of other neutron star low-mass X-ray binaries. In conjunction with the previously reported candidate X-ray burst, these confirm that the compact object in UW CrB is a neutron star. We examine the optical burst brightness and recurrence times and discuss how the nature of the system can be constrained. We conclude that the source is most likely an accretion disk corona source at an intermediate distance, rather than a nearby quiescent system or very distant dipper.

During the course of an ongoing CCD monitoring program to investigate low-level light variations in subdwarf B (sdB) stars, we have serendipitously discovered a new class of low amplitude, multimode sdB pulsators with periods of the order of an hour. These periods are more than a factor of ten longer than those of previously known multimode sdB pulsators (EC 14026 stars), implying that they are due to gravity modes rather than pressure modes. The longer period pulsators are found only among cooler sdB stars, where they are surprisingly common. The iron opacity instability that drives the short period EC 14026 stars is effective only in hot sdB's, leaving the driving mechanism for the deeper gravity modes in cool sdB's currently unknown. We present the first observational results for our newly identified sdB variables, and discuss possible implications.