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.

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.

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

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

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

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.

Bayesian statistical calculations and linear-bisector calculations for obtaining Cepheid distances and radii by the infrared surface brightness method have been compared for a set of 38 Cepheids. The distances obtained by the two techniques agree to 1.5%+/-0.6% and the radii agree to 1.1%+/-0.7%. Thus the two methods yield the same distances and radii at the 2 sigma level. This implies that the short distance to the LMC found in recent linearbisector studies of Cepheids is not a result of simplifications in the mathematical approach. The computed uncertainties in distance and radius are larger in the Bayesian calculation typically by a factor of three.

We have compared the results of Bayesian statistical calculations and linear-bisector calculations for obtaining Cepheid distances and radii by the infrared surface brightness method. We analyzed a set of 38 Cepheids using a Bayesian Markov Chain Monte Carlo method that had been recently studied with a linear-bisector method. The distances obtained by the two techniques agree to 1.5 \pm 0.6% with the Bayesian distances being larger. The radii agree to 1.1% \pm 0.7% with the Bayesian determinations again being larger. We interpret this result as demonstrating that the two methods yield the same distances and radii. This implies that the short distance to the LMC found in recent linear-bisector studies of Cepheids is not caused by deficiencies in the mathematical treatment. However, the computed uncertainties in distance and radius for our dataset are larger in the Bayesian calculation by factors of 1.4-6.7. We give reasons to favor the Bayesian computations of the uncertainties. The larger uncertainties can have a significant impact upon interpretation of Cepheid distances and radii obtained from the infrared surface brightness method.

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 present an absolute parallax and relative proper motion for the central star of the planetary nebula NGC 6853 (The Dumbell). We obtain these with astrometric data from FGS 3, a white-light interferometer on {\it HST}. Spectral classifications and VRIJHKT$_2$M and DDO51 photometry of the stars comprising the astrometric reference frame provide spectrophotometric estimates of their absolute parallaxes. Introducing these into our model as observations with error, we find $\pi_{abs} = 2.10 \pm 0.48$ mas for the DAO central star of NGC 6853. A weighted average with a previous ground-based USNO determination yields $\pi_{abs} = 2.40 \pm 0.32$. We assume that the extinction suffered by the reference stars nearest (in angular separation and distance) to the central star is the same as for the central star. Correcting for color differences, we find $<$A$_V>$ = 0.30 $ \pm $ 0.06 for the central star, hence, an absolute magnitude M$_V = 5.48^{-0.16}_{+0.15}$. A recent determination of the central star effective temperature aided in estimating the central star radius, $R_{*}=0.055 \pm 0.02R_{\sun}$, a star that may be descending to the white dwarf cooling track.

We develop and describe a Bayesian statistical analysis to solve the surface brightness equations for Cepheid distances and stellar properties. Our analysis provides a mathematically rigorous and objective solution to the problem, including immunity from Lutz-Kelker bias. We discuss the choice of priors, show the construction of the likelihood distribution, and give sampling algorithms in a Markov Chain Monte Carlo approach for efficiently and completely sampling the posterior probability distribution. Our analysis averages over the probabilities associated with several models rather than attempting to pick the `best model' from several possible models. Using a sample of thirteen Cepheids we demonstrate the method. We discuss diagnostics of the analysis and the effects of the astrophysical choices going into the model. We show that we can objectively model the order of Fourier polynomial fits to the light and velocity data. By comparison with theoretical models of Bono et al. (2001) we find that EU Tau and SZ Tau are overtone pulsators, most likely without convective overshoot. The period-radius and period-luminosity relations we obtain are shown to be compatible with those in the recent literature.

We report the first astrometrically determined mass of an extrasolar planet, a companion previously detected by Doppler spectroscopy. Radial velocities first provided an ephemeris with which to schedule a significant fraction of the {\it HST} observations near companion peri- and apastron. The astrometry residuals at these orbital phases exhibit a systematic deviation consistent with a perturbation due to a planetary mass companion. Combining {\it HST} astrometry with radial velocities, we solve for the proper motion, parallax, perturbation size, inclination, and position angle of the line of nodes, while constraining period, velocity amplitude, longitude of periastron, and eccentricity to values determined from radial velocities. We find a perturbation semi-major axis and inclination, $\alpha$ = 0.25 $\pm$ 0.06 mas, $i$ = 84\arcdeg $\pm$6\arcdeg, and Gl 876 absolute parallax, $\pi_{abs}= 214.6 \pm$ 0.2 mas. Assuming that the mass of the primary star is $M_* = 0.32M_{\sun}$, we find the mass of the planet, Gl 876b, $M_b = 1.89\pm0.34M_{Jup}$.

We use data from the BIMA Survey of Nearby Galaxies (SONG) to investigate the relationship between ellipticity and central mass concentration in barred spirals. Existing simulations predict that bar ellipticity decreases as inflowing mass driven by the bar accumulates in the central regions, ultimately destroying the bar. Using the ratio of the bulge mass to the mass within the bar radius as an estimate of the central mass concentration, we obtain dynamical mass estimates from SONG CO 1-0 rotation curve data. We find an inverse correlation between bar ellipticity and central mass concentration, consistent with simulations of bar dissolution.

We present an absolute parallax and relative proper motion for the fundamental distance scale calibrator, delta Cep. We obtain these with astrometric data from FGS 3, a white-light interferometer on HST. Utilizing spectrophotometric estimates of the absolute parallaxes of our astrometric reference stars and constraining delta Cep and reference star HD 213307 to belong to the same association (Cep OB6, de Zeeuw et al. 1999), we find pi_{abs} = 3.66 +/- 0.15 mas. The larger than typical astrometric residuals for the nearby astrometric reference star HD 213307 are found to satisfy Keplerian motion with P = 1.07 +/- 0.02 years, a perturbation and period that could be due to a F0V companion ~7 mas distant from and ~4 magnitudes fainter than the primary. Spectral classifications and VRIJHKT$_2$M and DDO51 photometry of the astrometric reference frame surrounding delta Cep indicate that field extinction is high and variable along this line of sight. However the extinction suffered by the reference star nearest (in angular separation and distance) to delta Cep, HD 213307, is lower and nearly the same as for delta Cep. Correcting for color differences, we find <A_V> = 0.23 +/- 0.03 for delta Cep, hence, an absolute magnitude M_V = -3.47 +/- 0.10. Adopting an average V magnitude, <V> = 15.03 +/- 0.03, for Cepheids with log P = 0.73 in the LMC from Udalski et al. (1999), we find a V-band distance modulus for the LMC, m-M = 18.50 +/- 0.13 or, 18.58 +/- 0.15, where the latter value results from a highly uncertain metallicity correction (Freedman et al. 2001). These agree with our previous RR Lyr HST parallax-based determination of the distance modulus of the LMC.

We present an absolute parallax and relative proper motion for the fundamental distance scale calibrator, RR Lyr. We obtain these with astrometric data from FGS 3, a white-light interferometer on HST. We find $\pi_{abs} = 3.82 \pm 0.2$ mas. Spectral classifications and VRIJHKT$_2$M and DDO51 photometry of the astrometric reference frame surrounding RR Lyr indicate that field extinction is low along this line of sight. We estimate <A_V>=0.07\pm0.03 for these reference stars. The extinction suffered by RR Lyr becomes one of the dominant contributors to the uncertainty in its absolute magnitude. Adopting the average field absorption, <A_V>=0.07 \pm 0.03, we obtain M_V^{RR} = 0.61 ^{-0.11}_{+0.10}. This provides a distance modulus for the LMC, m-M = 18.38 - 18.53^{-0.11}_{+0.10} with the average extinction-corrected magnitude of RR Lyr variables in the LMC, <V(RR)>, remaining a significant uncertainty. We compare this result to more than 80 other determinations of the distance modulus of the LMC.

TV Columbae (TV Col) is a 13th magnitude Intermediate Polar (IP) Cataclysmic Variable (CV), with multiple periods found in the light curves. Past estimates predicted a distance of 400 parsec to greater than 500 parsec. Recently completed Hubble Space Telescope (HST) Fine Guidance Sensor (FGS) interferometric observations allow us to determine the first trigonometric parallax to TV Col. This determination puts the distance of TV Col at 368 -15+17 parsecs. CD-32 2376, a 10th magnitude Tycho Catalog star, is a reference star in the TV Col frame. We find a distance of $127.7 -1+1 parsecs.

We present an analysis of astrometric data from FGS 3, a white-light interferometer on {\it HST}, and of radial velocity data from two ground-based campaigns. We model the astrometric and radial velocity measurements simultaneously to obtain parallax, proper motion and component masses for Wolf 1062 = Gl 748 AB (M3.5V). To derive the mass fraction, we relate FGS 3 fringe scanning observations of the science target to a reference frame provided by fringe tracking observations of a surrounding star field. We obtain an absolute parallax $\pi_{abs} = 98.0 \pm 0.4$ milliseconds of arc, yielding ${\cal M}_A = 0.379 \pm 0.005{\cal M}_{\sun}$ and ${\cal M}_B= 0.192 \pm 0.003 {\cal M}_{\sun}$, high quality component masses with errors of only 1.5%.

With HST FGS 3 we have determined a parallax for the white dwarf - M dwarf interacting binary, Feige 24. The white dwarf (DA) component has an effective temperature, T_eff of approximately 56,000 K degrees. A weighted average with past parallax determinations (pi_abs = 14.6 +- 0.4 milliseconds of arc) narrows the range of possible radius values, compared to past estimates. We obtain R_DA = 0.0185 +- 0.0008 R_solar with uncertainty in the temperature and bolometric correction the dominant contributors to the error. FGS photometry provides a light curve entirely consistent with reflection effects. A recently refined model Mass-Luminosity Relation (Baraffe et al. 1998) for low mass stars provides a mass estimate for the M dwarf companion, M_dM = 0.37 +- 0.20 M_solar, where the mass range is due to metallicity and age uncertainties. Radial velocities from Vennes and Thorstensen (1994) provide a mass ratio from which we obtain M_DA = 0.49+0.19-0.05 ~ M_solar. Independently, our radius and recent log_g determinations yield 0.44 < M_DA < 0.47 M_solar. In each case the minimum DA mass is that derived by Vennes & Thorstensen from their radial velocities and Keplerian circular orbits with inclination <= 90 degrees. Locating Feige 24 on an M - R plane suggests a carbon core. Our radius and these mass estimates yield a gamma_grav inconsistent with that derived by Vennes & Thorstensen. We speculate on the nature of a third component whose existence would resolve the discrepancy.

We report on a sub-stellar companion search utilizing interferometric fringe-tracking astrometry acquired with Fine Guidance Sensor 3 (FGS 3) on the Hubble Space Telescope. Our targets were Proxima Centauri and Barnard's Star. We obtain absolute parallax values for Proxima Cen pi_{abs} = 0.7687 arcsecond and for Barnard's Star pi_{abs} = 0.5454 arcsecond. Once low-amplitude instrumental systematic errors are identified and removed, our companion detection sensitivity is less than or equal to one Jupiter mass for periods longer than 60 days for Proxima Cen. Between the astrometry and the radial velocity results we exclude all companions with M > 0.8M_{Jup} for the range of periods 1 < P < 1000 days. For Barnard's Star our companion detection sensitivity is less than or equal to one Jupiter mass for periods long er than 150 days. Our null results for Barnard's Star are consistent with those of Gatewood (1995).

RW Triangulum (RW Tri) is a 13th magnitude Nova-like Cataclysmic Variable star with an orbital period of 0.2319 days (5.56 hours). Infrared observations of RW Tri indicate that its secondary is most likely a late K-dwarf. Past analyses predicted a distance of 270 parsec, derived from a black-body fit to spectrum of the central part of the disk. Recently completed Hubble Space Telescope Fine Guidance Sensor interferometric observations allow us to determine the first trigonometric parallax to RW Tri. This determination puts the distance of RW Tri at 341, one of the most distant objects with a direct parallax measurement.

We have observed Proxima Centauri and Barnard's Star with Hubble Space Telescope Fine Guidance Sensor 3. Proxima Centauri exhibits small-amplitude, periodic photometric variations. Once several sources of systematic photometric error are corrected, we obtain 2 milli-magnitude internal photometric precision. We identify two distinct behavior modes over the past four years: higher amplitude, longer period; smaller amplitude, shorter period. Within the errors one period (P ~ 83d) is twice the other. Barnard's Star shows very weak evidence for periodicity on a timescale of approximately 130 days. If we interpret these periodic phenomena as rotational modulation of star spots, we identify three discrete spots on Proxima Cen and possibly one spot on Barnard's Star. We find that the disturbances change significantly on time scales as short as one rotation period.