We present our analysis of a large sample (over 150k) of candidate Galactic RR Lyrae (RRL) stars for which we derived high quality photometry at UV, optical and infrared wavelengths, using data from publicly available surveys. For a sub-sample of these stars (~2,400 fundamental mode field RRLs) we have measured their individual metallicity using the Delta S method, resulting in the largest and most homogeneous spectroscopic data set collected for RRLs. We use this sample to study the metallicity distribution in the Galactic Halo, including the long-standing problem of the Oosterhoff dichotomy among Galactic globular clusters. We also analyze the dependence of their pulsation properties, and in particular the shape of their infrared light curves, from their [Fe/H] abundance.

We collected over 6000 high-resolution spectra of four dozen field RR Lyrae (RRL) variables pulsating either in the fundamental (39 RRab) or in the first overtone (9 RRc) mode. We measured radial velocities (RVs) of four strong metallic and four Balmer lines along the entire pulsational cycle and derived RV amplitudes with accuracies better than 1$-$2~\kmsec. The new amplitudes were combined with literature data for 23~RRab and 3~RRc stars (total sample 74 RRLs) which allowed us to investigate the variation of the Bailey diagram (photometric amplitude versus period) when moving from optical to mid-infrared bands and to re-cast the Bailey diagram in terms of RV amplitudes. We found that RV amplitudes for RRab are minimally affected by nonlinear phenomena (shocks) and multi-periodicity (Blazhko effect). The RV slope ($\log P$--A(V$_r$)) when compared with the visual slope ($\log P$--A($V$)) is shallower and the dispersion, at fixed period, decreases by a factor of two. We constructed homogeneous sets of Horizontal Branch evolutionary models and nonlinear, convective pulsation models of RRLs to constrain the impact of evolutionary effects on their pulsation properties. Evolution causes, on the Bailey diagram based on RV amplitudes, a modest variation in pulsation period and a large dispersion in amplitude. The broad dispersion in period of the Bailey diagram is mainly caused by variation in RRL intrinsic parameters (stellar mass, chemical composition). Empirical evidence indicates that RV amplitudes are an optimal diagnostic for tracing the mean effective temperature across the RRab instability strip.

We examine four ancient Large Magellanic Cloud (LMC) globular clusters (GCs) for evidence of multiple stellar populations using the Advanced Camera for Surveys and Wide Field Camera 3 on the Hubble Space Telescope Programme GO-14164. NGC 1466, NGC 1841, and NGC 2257 all show evidence for a redder, secondary population along the main-sequence. Reticulum does not show evidence for the presence of a redder population, but this GC has the least number of stars and Monte Carlo simulations indicate that the sample of main sequence stars is too small to robustly infer whether a redder population exists in this cluster. The second, redder, population of the other three clusters constitutes $\sim30-40\%$ of the total population along the main-sequence. This brings the total number of ancient LMC GCs with known split or broadened main-sequences to five. However, unlike for Hodge 11 and NGC 2210 (see arXiv:1904.01434), none of the clusters show evidence for multiple populations in the horizontal branch. We also do not find evidence of a second population along the Red Giant Branch (RGB).

We collected a large data set of field RR Lyrae stars (RRLs) by using catalogues already available in the literature and Gaia DR2. We estimated the iron abundances for a sub-sample of 2,382 fundamental RRLs ($\Delta$S method: CaIIK, H$\beta$, H$\gamma$ and H$\delta$ lines) for which are publicly available medium-resolution SDSS-SEGUE spectra. We also included similar estimates available in the literature ending up with the largest and most homogeneous spectroscopic data set ever collected for RRLs (2,903). The metallicity scale was validated by using iron abundances based on high resolution spectra for a fundamental field RRL (V~Ind), for which we collected X-shooter spectra covering the entire pulsation cycle. The peak ([Fe/H]=-1.59$\pm$0.01) and the standard deviation ($\sigma$=0.43 dex) of the metallicity distribution agree quite well with similar estimates available in the literature. The current measurements disclose a well defined metal-rich tail approaching Solar iron abundance. The spectroscopic sample plotted in the Bailey diagram (period vs luminosity amplitude) shows a steady variation when moving from the metal-poor ([Fe/H]=-3.0/-2.5) to the metal-rich ([Fe/H]=-0.5/0.0) regime. The smooth transition in the peak of the period distribution as a function of the metallicity strongly indicates that the long-standing problem of the Oosterhoff dichotomy among Galactic globulars is the consequence of the lack of metal-intermediate clusters hosting RRLs. We also found that the luminosity amplitude, in contrast with period, does not show a solid correlation with metallicity. This suggests that period-amplitude-metallicity relations should be cautiously treated.

We present multi-epoch infrared photometry in the $K_s$-band for 74 bright RR Lyrae variable stars tied directly to the 2MASS photometric system. We systematize additional $K$-band photometry from the literature to the 2MASS system and combine it to obtain photometry for 146 RR Lyrae stars on a consistent, modern system. A set of outlier stars in the literature photometry is identified and discussed. Reddening estimates for each star were gathered from the literature and combined to provide an estimate of the interstellar absorption affecting each star, and we find excellent agreement with another source in the literature. We utilize trigonometric parallaxes from the Second Data Release of ESA's Gaia astrometric satellite to determine the absolute magnitude, $M_{K_s}$ for each of these stars, and analyze them using the astrometry based luminosity prescription to obtain a parallax-based calibration of $M_K$(RR). Our period-luminosity-metallicity relationship is $ M_{K_s} = (-2.8\pm 0.2) (\log P +0.27) + (0.12\pm 0.02) ( [\mathrm{Fe/H}] + 1.3) - (0.41\pm 0.03)$ mag. A Gaia global zero-point error of $\pi_{zp} = -0.042\pm 0.013\,$mas is determined for this sample of RR Lyrae stars.

Radial velocities of fifteen double-mode bulge RR Lyrae (RR01) stars are presented, six of which belong to a compact group of RR01 stars in pulsation space, with the ratio of first-overtone period to fundamental-mode period, P_{fo}/P_{f}~0.74, and P_{f}~0.44. It has been suggested that these pulsationally clumped RR01 stars are a relic of a disrupted dwarf galaxy or stellar cluster, as they also appear to be spatially coherent in a vertical strip across the bulge. However, the radial velocities of the stars presented here, along with proper motions from Gaia DR2, show a large range of radial velocities, proper motions and distances for the bulge RR01 stars in the pulsation clump, much larger than the RR01 stars in the Sagittarius dwarf galaxy (Sgr). Therefore, in contrast to the kinematics of the RRL stars belonging to Sgr, and those in and surrounding the bulge globular cluster NGC~6441, there is no obvious kinematic signature within the pulsationally clumped RR01 stars. If the pulsationally clumped RR01 stars belonged to the same system in the past and were accreted, their accretion in the inner Galaxy was not recent, as the kinematic signature of this group has been lost (i.e., these stars are now well-mixed within the inner Galaxy). We show that the apparent spatial coherence reported for these stars could have been caused by small number statistics. The orbits of the RR01 stars in the inner Galaxy suggest they are confined to the innermost ~4~kpc of the Milky Way.

We present a multiple population search in two old Large Magellanic Cloud (LMC) Globular Clusters, Hodge 11 and NGC 2210. This work uses data from the Advanced Camera for Surveys and Wide Field Camera 3 on the Hubble Space Telescope from programme GO-14164 in Cycle 23. Both of these clusters exhibit a broadened main sequence with the second population representing ($20 \pm \! \sim \! 5$)% for NGC 2210 and ($30 \pm \! \sim \! 5$)% for Hodge 11. In both clusters, the smaller population is redder than the primary population, suggesting CNO variations. Hodge 11 also displays a bluer second population in the horizontal branch, which is evidence for helium enhancement. However, even though NGC 2210 shows similarities to Hodge 11 in the main sequence, there does not appear to be a second population on NGC 2210's horizontal branch. This is the first photometric evidence that ancient LMC Globular Clusters exhibit multiple stellar populations.

We provide homogeneous optical (UBVRI) and near-infrared (JHK) time series photometry for 254 cluster (omega Cen, M4) and field RR Lyrae (RRL) variables. We ended up with more than 551,000 measurements. For 94 fundamental (RRab) and 51 first overtones (RRc) we provide a complete optical/NIR characterization (mean magnitudes, luminosity amplitudes, epoch of the anchor point). The NIR light curves of these variables were adopted to provide new and accurate light-curve templates for both RRc (single period bin) and RRab (three period bins) variables. The templates for the J and the H band are newly introduced, together with the use of the pulsation period to discriminate among the different RRab templates. To overcome subtle uncertainties in the fit of secondary features of the light curves we provide two independent sets of analytical functions (Fourier series, Periodic Gaussian functions). The new templates were validated by using 26 omega Cen and Bulge RRLs covering the four period bins. We found that the difference between the measured mean magnitude along the light curve and the mean magnitude estimated by using the template on a single randomly extracted phase point is better than 0.01 mag (sigma=0.04 mag). We also validated the template on variables for which at least three phase points were available, but without information on the phase of the anchor point. The accuracy of the mean magnitudes is ~0.01 mag (sigma=0.04 mag). The new templates were applied to the Large Magellanic Cloud (LMC) globular Reticulum and by using literature data and predicted PLZ relations we found true distance moduli of 18.47+-0.10+-0.03 mag (J) and 18.49+-0.09+-0.05 mag (K). We also used literature optical and mid-infrared data and we found a mean true distance modulus of 18.47+-0.02+-0.06 mag, suggesting that Reticulum is ~1 kpc closer than the LMC.

The main aim of this experiment is to provide a complete census of old (t > 10 Gyr, RR Lyrae, type II Cepheids, red horizontal branch), intermediate age (red clump, Miras) and young (classical Cepheids) stellar tracers across the Galactic Bulge. To fully exploit the unique photometric quality of LSST images, we plan to perform a Shallow minisurvey (ugrizy, -20 < l < 20 deg, -15 < b < 10 deg) and a Deep minisurvey (izy, -20 < l < 20 deg, -3 < b < 3 deg). The former one is aimed at constraining the 3D structure of the galactic Bulge across the four quadrants, and in particular, the transition between inner and outer Bulge. The u,g,r,i,z,y LSST bands provide fundamental diagnostics to constrain the evolutionary properties of low and intermediate-mass stars when moving from a metal-poor to a metal-rich regime. The deep minisurvey is aimed at tracing RR Lyrae, Red Clump stars, Miras and classical Cepheids in highly reddened regions of the Galactic center. These images will allow us to investigate the role that baryonic mass and dark matter played in the early formation and evolution of the MW.

We developed a new approach to provide accurate estimates of metal content, reddening and true distance modulus of RR Lyrae stars (RRLs). The method is based on homogeneous optical ($BVI$) and near-infrared ($JHK$) mean magnitudes and on predicted period--luminosity--metallicity relations ($IJHK$) and absolute mean magnitude--metallicity relations ($BV$). We obtained solutions for three different RRL samples in $\omega$ Cen: first overtone (RRc,~90), fundamental (RRab,~80) and global (RRc+RRab) in which the period of first overtones were fundamentalized. The metallicity distribution shows a well defined peak at [Fe/H]$\sim$--1.98 and a standard deviation of $\sigma$=0.54 dex. The spread is, as expected, metal-poor ([Fe/H]$\le$--2.3) objects. The current metallicity distribution is $\sim$0.3 dex more metal-poor than similar estimates for RRLs available in the literature. The difference vanishes if the true distance modulus we estimated is offset by --0.06/--0.07~mag in true distance modulus. We also found a cluster true distance modulus of $\mu$=13.720$\pm$0.002$\pm$0.030~mag, where the former error is the error on the mean and the latter is the standard deviation. Moreover, we found a cluster reddening of E($B-V$)=0.132$\pm$0.002$\pm$0.028~mag and spatial variations of the order of a few arcmin across the body of the cluster. Both the true distance modulus and the reddening are slightly larger than similar estimates available in the literature, but the difference is within 1$\sigma$. The metallicity dependence of distance diagnostics agree with theory and observations, but firm constraints require accurate and homogeneous spectroscopic measurements.

The bright, nearby binary $\alpha$ Centauri provides an excellent laboratory for testing stellar evolution models, as it is one of the few stellar systems for which we have high-precision classical (mass, radius, luminosity) and asteroseismic ($p$-mode) observations. Stellar models are created and fit to the classical and seismic observations of both stars by allowing for the free variation of convective mixing length parameter $\alpha_{\text{MLT}}$. This system is modeled using five different sets of assumptions about the physics governing the stellar models. There are 31 pairs of tracks (out of ${\sim} 150,000$ generated) which fit the classical, binary, and seismic observational constraints of the system within $3\,\sigma$. Models with each tested choice of input physics are found to be viable, but the optimal mixing lengths for Cen A and Cen B remain the same regardless of the physical prescription. The optimal mixing lengths are $\alpha_{\text{MLT,A}} /\alpha_{\odot}= 0.932$ and $\alpha_{\text{MLT,B}}/\alpha_{\odot} = 1.095$. That Cen A and Cen B require sub- and super-solar mixing lengths, respectively, to fit the observations is a trend consistent with recent findings, such as in Kervella et al. (2017), Joyce and Chaboyer (2018), and Viani et al. (2018). The optimal models find an age for $\alpha$ Centauri of $5.3 \pm 0.3$ Gyr.

Theoretical stellar evolution models are constructed and tailored to the best known, observationally derived characteristics of metal-poor ([Fe/H]$\sim-2.3$) stars representing a range of evolutionary phases: subgiant HD140283, globular cluster M92, and four single, main sequence stars with well-determined parallaxes: HIP46120, HIP54639, HIP106924, and WOLF1137. It is found that the use of a solar-calibrated value of the mixing length parameter $\alpha_{\text{MLT}}$ in models of these objects is ineffective at reproducing their observed properties. Empirically calibrated values of $\alpha_{\text{MLT}}$ are presented for each object, accounting for uncertainties in the input physics employed in the models. It is advocated that the implementation of an adaptive mixing length is necessary in order for stellar evolution models to maintain fidelity in the era of asteroseismic observations.

We analyze Hubble Space Telescope observations of six globular clusters in the Large Magellanic Cloud from program GO-14164 in Cycle 23. These are the deepest available observations of the LMC globular cluster population; their uniformity facilitates a precise comparison with globular clusters in the Milky Way. Measuring the magnitude of the main sequence turnoff point relative to template Galactic globular clusters allows the relative ages of the clusters to be determined with a mean precision of 8.4%, and down to 6% for individual objects. We find that the mean age of our LMC cluster ensemble is identical to the mean age of the oldest metal-poor clusters in the Milky Way halo to 0.2 $\pm$ 0.4 Gyr. This provides the most sensitive test to date of the synchronicity of the earliest epoch of globular cluster formation in two independent galaxies. Horizontal branch magnitudes and subdwarf fitting to the main sequence allow us to determine distance estimates for each cluster, and examine their geometric distribution in the LMC. Using two different methods, we find an average distance to the LMC of 18.52 $\pm$ 0.05.

We obtain high resolution spectra of nine red giant branch stars in NGC 6681 and perform the first detailed abundance analysis of stars in this cluster. We confirm cluster membership for these stars based on consistent radial velocities of 214.5$\pm$3.7 km/s and find a mean [Fe/H] = -1.63$\pm$0.07 dex and [{\alpha}/Fe] = 0.42$\pm$0.11 dex. Additionally, we confirm the existence of a Na-O anti-correlation in NGC 6681 and identify two populations of stars with unique abundance trends. With the use of HST photometry from Sarajedini et al. (2007) and Piotto et al. (2015) we are able to identify these two populations as discrete sequences in the cluster CMD. Although we cannot confirm the nature of the polluter stars responsible for the abundance differences in these populations, these results do help put constraints on possible polluter candidates.

The recent Gaia Data Release 1 of stellar parallaxes provides ample opportunity to find metal-poor main-sequence stars with precise parallaxes. We select 21 such stars with parallax uncertainties better than $\sigma_\pi/\pi\leq0.10$ and accurate abundance determinations suitable for testing metal-poor stellar evolution models and determining the distance to Galactic globular clusters. A Monte Carlo analysis was used, taking into account uncertainties in the model construction parameters, to generate stellar models and isochrones to fit to the calibration stars. The isochrones which fit the calibration stars best were then used to determine the distances and ages of 22 globular clusters with metallicities ranging from -2.4 dex to -0.7 dex. We find distances with an average uncertainty of 0.15 mag and absolute ages ranging from 10.8 - 13.6 Gyr with an average uncertainty of 1.6 Gyr. Using literature proper motion data we calculate orbits for the clusters finding six that reside within the Galactic disk/bulge while the rest are considered halo clusters. We find no strong evidence for a relationship between age and Galactocentric distance, but we do find a decreasing age-[Fe/H] relation.

We have performed a differential, line-by-line, chemical abundance analysis, ultimately relative to the Sun, of nine very metal-poor main sequence halo stars, near [Fe/H]=$-$2 dex. Our abundances range from $-2.66\leq\mathrm{[Fe/H]}\leq-1.40$ dex with conservative uncertainties of 0.07 dex. We find an average [$\alpha$/Fe]$=0.34\pm0.09$ dex, typical of the Milky Way. While our spectroscopic atmosphere parameters provide good agreement with HST parallaxes, there is significant disagreement with temperature and gravity parameters indicated by observed colors and theoretical isochrones. Although a systematic underestimate of the stellar temperature by a few hundred degrees could explain this difference, it is not supported by current effective temperature studies and would create large uncertainties in the abundance determinations. Both 1D and $\langle$3D$\rangle$ hydrodynamical models combined with separate 1D non-LTE effects do not yet account for the atmospheres of real metal-poor MS stars, but a fully 3D non-LTE treatment may be able to explain the ionization imbalance found in this work.

Hubble Space Telescope (HST) fine guidance sensor observations were used to obtain parallaxes of eight metal-poor ([Fe/H] < -1.4) stars. The parallaxes of these stars determined by the revised Hipparcos reduction average 17% accuracy, in contrast to our new HST parallaxes which average 1% accuracy and have errors on the individual parallaxes ranging from 85 to 144 microarcsecond. This parallax data has been combined with HST ACS photometry in the F606W and F814W filters to obtain the absolute magnitudes of the stars with an accuracy of 0.02 to 0.03 magnitudes. Six of these stars are on the main sequence (with -2.7 < [Fe/H] < -1.8), and suitable for testing metal-poor stellar evolution models and determining the distances to metal-poor globular clusters. Using the abundances obtained by O'Malley et al. (2017) we find that standard stellar models using the Vandenberg & Clem (2003) color transformation do a reasonable job of matching five of the main sequence stars, with HD 54639 ([Fe/H] = -2.5) being anomalous in its location in the color-magnitude diagram. Stellar models and isochrones were generated using a Monte Carlo analysis to take into account uncertainties in the models. Isochrones which fit the parallax stars were used to determine the distances and ages of nine globular clusters (with -2.4 <= [Fe/H] <= -1.9$). Averaging together the age of all nine clusters, leads to an absolute age of the oldest, most metal-poor globular clusters of 12.7+/- 1.0 Gyr, where the quoted uncertainty takes into account the known uncertainties in the stellar models and isochrones, along with the uncertainty in the distance and reddening of the clusters.

Synthetic RGBB magnitudes are generated with the most recent theoretical stellar evolution models computed with the Dartmouth Stellar Evolution Program (DSEP) code. They are compared to the observational work of Nataf et al., who present RGBB magnitudes for 72 globular clusters. A DSEP model using a chemical composition with enhanced $\alpha$ capture [$\alpha$/Fe] $ =+0.4$ and an age of 13 Gyr shows agreement with observations over metallicities ranging from [Fe/H] = $0$ to [Fe/H] $\approx-1.5$, with discrepancy emerging at lower metallicities.

We report an age revision for the low-mass detached eclipsing binary CM Draconis and its common proper motion companion, WD 1633+572. An age of 8.5 $\pm$ 3.5 Gyr is found by combining an age estimate for the lifetime of WD 1633+572 and an estimate from galactic space motions. The revised age is greater than a factor of two older than previous estimates. Our results provide consistency between the white dwarf age and the system's galactic kinematics, which reveal the system is a highly probable member of the galactic thick disk. We find the probability that CM Draconis and WD 1633+572 are members of the thick disk is 8500 times greater than the probability that they are members of the thin disk and 170 times greater than the probability they are halo interlopers. If CM Draconis is a member of the thick disk, it is likely enriched in $\alpha$-elements compared to iron by at least 0.2 dex relative to the Sun. This leads to the possibility that previous studies under-estimate the [Fe/H] value, suggesting the system has a near-solar [Fe/H]. Implications for the long-standing discrepancies between the radii of CM Draconis and predictions from stellar evolution theory are discussed. We conclude that CM Draconis is only inflated by about 2% compared to stellar evolution predictions.

We present the current status of an effort to create an updated grid of low-mass stellar evolution mass tracks and isochrones computed using the Dartmouth stellar evolution code. Emphasis is placed on reliably extending the present grid to the pre-main-sequence, where modeling uncertainties have the greatest impact. Revisions to the original code release include: updated surface boundary conditions, the introduction of deuterium burning, and magnetic fields. The mass track grid contains models with a mass above 0.1 $M_{\odot}$ and metallicities in the range of -0.5 to +0.5 dex. Magnetic mass tracks are calculated for surface magnetic field strengths between 0.1 kG and 4.0 kG using two different prescriptions for magneto-convection. Standard and magnetic model isochrones are available for ages older than 1 Myr. Tabulated quantities include the stellar fundamental properties, absolute photometric magnitudes, magnetic field properties, convective turnover times, apsidal motion constants, and lithium abundances. The complete grid of mass tracks and isochrones will be made publicly available.

We present new $BVI$ photometry containing 40,000 stars of the Galactic globular cluster M15 (NGC 7078) covering a $25 \times 25$ arcminute area centered on the cluster with a magnitude range from the tip of the Red Giant Branch to three magnitudes below the main sequence turn-off. Using $\alpha$-enhanced Dartmouth Stellar Evolution Program models, we find an age of $13.0 \ \pm \ 1.0$ Gyr and distance modulus of $(m-M)_V = 15.4 \pm 0.1$ through isochrone fitting. Unlike previous works, we find good agreement between the observed completeness-corrected stellar luminosity function and models.

We examine the hypothesis that magnetic fields are inflating the radii of fully convective main sequence stars in detached eclipsing binaries (DEBs). The magnetic Dartmouth stellar evolution code is used to analyze two systems in particular: Kepler-16 and CM Draconis. Magneto-convection is treated assuming stabilization of convection and also by assuming reductions in convective efficiency due to a turbulent dynamo. We find that magnetic stellar models are unable to reproduce the properties of inflated fully convective main sequence stars, unless strong interior magnetic fields in excess of 10 MG are present. Validation of the magnetic field hypothesis given the current generation of magnetic stellar evolution models therefore depends critically on whether the generation and maintenance of strong interior magnetic fields is physically possible. An examination of this requirement is provided. Additionally, an analysis of previous studies invoking the influence of star spots is presented to assess the suggestion that star spots are inflating stars and biasing light curve analyses toward larger radii. From our analysis, we find that there is not yet sufficient evidence to definitively support the hypothesis that magnetic fields are responsible for the observed inflation among fully convective main sequence stars in DEBs.

In this study we compare the photometric data of 34 Milky Way globular clusters, observed within the ACS Treasury Program (PI: Ata Sarajedini) with the corresponding ground-based data, provided by the Photometric Standard Field Catalogs of Stetson (2000, 2005). We focus on the transformation between the HST/ACS F606W to V-band and F814W to I-band only. The goal is to assess the validity of the filter transformation equations by Sirianni et al.(2005) with respect to their dependence on metallicity, Horizontal Branch morphology, mass and integrated (V-I) colour of the various globular clusters. Such a dependence is expected due to the fact that the transformation equations are based on the observations of only one globular cluster, i.e., NGC 2419. Surprisingly, the correlation between offset and metallicity is found to be weak, with a low level significance. The correlation between offset and Horizontal Branch structure, as well as total cluster mass is still weaker. Based on the available data we do not find the photometric offset to be linked to multiple stellar populations, e.g., as found in NGC 0288, NGC 1851, and NGC 5139. The results of this study show that there are small systematic offsets between the transformed ACS- and observed ground based photometry, and that these are only weakly correlated, if at all, with various cluster parameters and their underlying stellar populations. As a result, investigators wishing to transform globular cluster photometry from the Sirianni et al.(2005) ground-based V, I system onto the Stetson (2000) system simply need to add 0.040 (+/-0.012) to the V-band magnitudes and 0.047 (+/-0.011) to the I-band magnitudes. This in turn means that the transformed ACS (V-I) colours match the ground-based values from Stetson (2000) to within ~0.01 mag.

Magnetic fields are hypothesized to inflate the radii of low-mass stars---defined as less massive than 0.8$M_\odot$---in detached eclipsing binaries (DEBs). We investigate this hypothesis using the recently introduced magnetic Dartmouth stellar evolution code. In particular, we focus on stars thought to have a radiative core and convective outer envelope by studying in detail three individual DEBs: UV Psc, YY Gem, and CU Cnc. The results suggest that the stabilization of thermal convection by a magnetic field is a plausible explanation for the observed model-radius discrepancies. However, surface magnetic field strengths required by the models are significantly stronger than those estimated from the observed coronal X-ray emission. Agreement between model predicted surface magnetic field strengths and those inferred from X-ray observations can be found by assuming that the magnetic field sources its energy from convection. This approach makes the transport of heat by convection less efficient and is akin to reduced convective mixing length methods used in other studies. Predictions for the metallicity and magnetic field strengths of the aforementioned systems are reported. We also develop an expression relating a reduction in the convective mixing length to a magnetic field strength in units of the equipartition value. Our results are compared with those from previous investigations to incorporate magnetic fields to explain the low-mass DEB radius inflation. Finally, we explore how the effects of magnetic fields might affect mass determinations using asteroseismic data and the implication of magnetic fields on exoplanet studies.

We present results obtained using the magnetic Dartmouth stellar evolution code that address the possibility that magnetic fields are inflating low-mass stars in detached eclipsing binaries. While it seems plausible that magnetic fields are inflating stars with radiative cores, the level of inflation observed among fully convective stars appears too large to be explained by magnetic fields. We provide an alternative explanation, stellar metallicity, and propose observations that can help further constrain stellar models.

Magnetic fields have been hypothesized to inflate the radii of low-mass stars---defined as less than 0.8 solar masses---in detached eclipsing binaries (DEBs). We evaluate this hypothesis using the magnetic Dartmouth stellar evolution code. Results suggest that magnetic suppression of thermal convection can inflate low-mass stars that possess a radiative core and convective outer envelope. A scaling relation between X-ray luminosity and surface magnetic flux indicates that model surface magnetic field strength predictions are consistent with observations. This supports the notion that magnetic fields may be inflating these stars. However, magnetic models are unable to reproduce radii of fully convective stars in DEBs. Instead, we propose that model discrepancies below the fully convective boundary are related to metallicity.

The first calibrated broadband UBVI time-series photometry is presented for the RR Lyrae variable stars in NGC6656 (M22), with observations spanning a range of twenty-two years. We have also redetermined the variability types and periods for the RR Lyrae stars identified previously by photographic observations, revising the number of fundamental-mode RR Lyrae variables (RR0) to 10 and the number of first-overtone variables (RR1) to 16. The mean periods of the RR0 and RR1 variables are <P>_{RR0}=0.66 \pm 0.02 d and <P>_{RR1}=0.33 \pm 0.01 d, respectively, supporting an Oosterhoff II classification for the cluster. The number ratio of RR1- to all RR-type variables is N_1/N_{RR}=0.61, also consistent with an Oosterhoff II designation. Both the RR Lyrae stars' minimum light colors and the blue edge of the RR Lyrae instability strip suggest E(B-V)=0.36 \pm 0.02 mag toward M22. Regarding the HB morphology of M22, we find (B-R)/(B+V+R)=+0.97 \pm 0.1 and at least one "gap" located in an unusual part of the blue HB, in the middle of the so-called hot HB stars.

We introduce a new one-dimensional stellar evolution code, based on the existing Dartmouth code, that self-consistently accounts for the presence of a globally pervasive magnetic field. The methods involved in perturbing the equations of stellar structure, the equation of state, and the mixing-length theory of convection are presented and discussed. As a first test of the code's viability, stellar evolution models are computed for the components of a solar-type, detached eclipsing binary (DEB) system, EF Aquarii, shown to exhibit large disagreements with stellar models. The addition of the magnetic perturbation corrects the radius and effective temperature discrepancies observed in EF Aquarii. Furthermore, the required magnetic field strength at the model photosphere is within a factor of two of the magnetic field strengths estimated from the stellar X-ray luminosities measured by ROSAT and those predicted from Ca II K line core emission. These models provide firm evidence that the suppression of thermal convection arising from the presence of a magnetic field is sufficient to significantly alter the structure of solar-type stars, producing noticeably inflated radii and cooler effective temperatures. The inclusion of magnetic effects within a stellar evolution model has a wide range of applications, from DEBs and exoplanet host stars to the donor stars of cataclysmic variables.

We examine the agreement between the observed and theoretical low-mass (< 0.8 solar masses) stellar main sequence mass-radius relationship by comparing detached eclipsing binary (DEB) data with a new, large grid of stellar evolution models. The new grid allows for a realistic variation in the age and metallicity of the DEB population, characteristic of the local galactic neighborhood. Overall, our models do a reasonable job of reproducing the observational data. A large majority of the models match the observed stellar radii to within 4%, with a mean absolute error of 2.3%. These results represent a factor of two improvement compared to previous examinations of the low-mass mass-radius relationship. The improved agreement between models and observations brings the radius deviations within the limits imposed by potential starspot-related uncertainties for 92% of the stars in our DEB sample.

The fraction of binary stars is an important ingredient to interpret globular cluster dynamical evolution and their stellar population. We investigate the properties of main-sequence binaries measured in a uniform photometric sample of 59 Galactic globular clusters that were observed by HST WFC/ACS as a part of the Globular Cluster Treasury project. We measured the fraction of binaries and the distribution of mass-ratio as a function of radial location within the cluster, from the central core to beyond the half-mass radius. We studied the radial distribution of binary stars, and the distribution of stellar mass ratios. We investigated monovariate relations between the fraction of binaries and the main parameters of their host clusters. We found that in nearly all the clusters, the total fraction of binaries is significantly smaller than the fraction of binaries in the field, with a few exceptions only. Binary stars are significantly more centrally concentrated than single MS stars in most of the clusters studied in this paper. The distribution of the mass ratio is generally flat (for mass-ratio parameter q>0.5). We found a significant anti-correlation between the binary fraction in a cluster and its absolute luminosity (mass). Some, less significant correlation with the collisional parameter, the central stellar density, and the central velocity dispersion are present. There is no statistically significant relation between the binary fraction and other cluster parameters. We confirm the correlation between the binary fraction and the fraction of blue stragglers in the cluster.

We present new absolute trigonometric parallaxes and proper motions for seven Pop II variable stars: five RR Lyr variables; RZ Cep, XZ Cyg, SU Dra, RR Lyr, UV Oct; and two type 2 Cepheids; VY Pyx and {\kappa} Pav. We obtained these results with astrometric data from Fine Guidance Sensors, white-light interfer- ometers on Hubble Space Telescope. We find absolute parallaxes in milliseconds of arc: RZ Cep, 2.12\pm0.16 mas; XZ Cyg, 1.67\pm0.17 mas; SU Dra, 1.42\pm0.16 mas; RR Lyr, 3.77\pm0.13 mas; UV Oct, 1.71\pm0.10 mas; VY Pyx, 6.44\pm0.23 mas; and {\kappa} Pav, 5.57\pm0.28 mas; an average {\sigma}{\pi}/{\pi} = 5.4%. With these parallaxes we compute absolute magnitudes in V and K bandpasses corrected for interstellar extinction and Lutz-Kelker-Hanson bias. Using these RRL absolute magnitudes, we then derive zero-points for MV -[Fe/H] and MK -[Fe/H]-Log P relations. The technique of reduced parallaxes corroborates these results. We employ our new results to determine distances and ages of several Galactic globular clusters and the dis- tance of the LMC. The latter is close to that previously derived from Classical Cepheids uncorrected for any metallicity effect, indicating that any such effect is small. We also discuss the somewhat puzzling results obtained for our two type 2 Cepheids.

The recent discovery of an eclipsing hierarchical triple system with two low-mass stars in a close orbit (KOI-126) by Carter et al. (2011) appeared to reinforce the evidence that theoretical stellar evolution models are not able to reproduce the observational mass-radius relation for low-mass stars. We present a set of stellar models for the three stars in the KOI-126 system that show excellent agreement with the observed radii. This agreement appears to be due to the equation of state implemented by our code. A significant dispersion in the observed mass-radius relation for fully convective stars is demonstrated; indicative of the influence of physics currently not incorporated in standard stellar evolution models. We also predict apsidal motion constants for the two M-dwarf companions. These values should be observationally determined to within 1% by the end of the Kepler mission.

We use observations from the ACS study of Galactic globular clusters to investigate the spatial distribution of the inner regions of the disrupting Sagittarius dwarf spheroidal galaxy (Sgr). We combine previously published analyses of four Sgr member clusters located near or in the Sgr core (M54, Arp 2, Terzan 7 and Terzan 8) with a new analysis of diffuse Sgr material identified in the background of five low-latitude Galactic bulge clusters (NGC 6624, 6637, 6652, 6681 and 6809) observed as part of the ACS survey. By comparing the bulge cluster CMDs to our previous analysis of the M54/Sgr core, we estimate distances to these background features. The combined data from four Sgr member clusters and five Sgr background features provides nine independent measures of the Sgr distance and, as a group, provide uniformly measured and calibrated probes of different parts of the inner regions of Sgr spanning twenty degrees over the face of the disrupting dwarf. This allows us, for the first time, to constrain the three dimensional orientation of Sgr's disrupting core and globular cluster system and compare that orientation to the predictions of an N-body model of tidal disruption. The density and distance of Sgr debris is consistent with models that favor a relatively high Sgr core mass and a slightly greater distance (28-30 kpc, with a mean of 29.4 kpc). Our analysis also suggests that M54 is in the foreground of Sgr by ~2 kpc, projected on the center of the Sgr dSph. While this would imply a remarkable alignment of the cluster and the Sgr nucleus along the line of sight, we can not identify any systematic effect in our analysis that would falsely create the measured 2 kpc separation. Finally, we find that the cluster Terzan 7 has the most discrepant distance (25 kpc) among the four Sgr core clusters, which may suggest a different dynamical history than the other Sgr core clusters.

We present a new analysis of the Jupiter+Saturn analog system, OGLE-2006-BLG-109Lb,c, which was the first double planet system discovered with the gravitational microlensing method. This is the only multi-planet system discovered by any method with measured masses for the star and both planets. In addition to the signatures of two planets, this event also exhibits a microlensing parallax signature and finite source effects that provide a direct measure of the masses of the star and planets, and the expected brightness of the host star is confirmed by Keck AO imaging, yielding masses of M_* = 0.51(+0.05-0.04) M_sun, M_b = 231+-19 M_earth, M_c = 86+-7 M_earth. The Saturn-analog planet in this system had a planetary light curve deviation that lasted for 11 days, and as a result, the effects of the orbital motion are visible in the microlensing light curve. We find that four of the six orbital parameters are tightly constrained and that a fifth parameter, the orbital acceleration, is weakly constrained. No orbital information is available for the Jupiter-analog planet, but its presence helps to constrain the orbital motion of the Saturn-analog planet. Assuming co-planar orbits, we find an orbital eccentricity of eccentricity = 0.15 (+0.17-0.10) and an orbital inclination of i = 64 (+4-7) deg. The 95% confidence level lower limit on the inclination of i > 49 deg. implies that this planetary system can be detected and studied via radial velocity measurements using a telescope of >30m aperture.

The open cluster M67 has solar metallicity and an age of about 4Gyr. The turn-off mass is close to the minimum mass for which solar metallicity stars develop a convective core during main sequence evolution as a result of the development of hydrogen burning through the CNO-cycle. The morphology of the color-magnitude-diagram (CMD) of M67 around the turn-off shows a clear hook-like feature, direct sign that stars close to the turn-off have convective cores. VandenBerg et al. investigated the possibility of using the morphology of the M67 turn-off to put constraints on the solar metallicity, particularly CNO elements, for which solar abundances have been revised downwards by more than 30% over the last few years. Here, we extend their work filling in the gaps in their analysis. To this aim, we compute isochrones appropriate for M67 using new (low metallicity) and old (high metallicity) solar abundances and study whether the characteristic turn-off in the CMD of M67 can be reproduced or not. We also study the importance of other constitutive physics on determining the presence of such a hook, particularly element diffusion, overshooting and nuclear reaction rates. We find that using the new solar abundance determinations, with low CNO abundances, makes it more difficult to reproduce the characteristic CMD of M67. This result is in agreement with results by VandenBerg et al. However, changes in the constitutive physics of the models, particularly overshooting, can influence and alter this result to the extent that isochrones constructed with models using low CNO solar abundances can also reproduce the turn-off morphology in M67. We conclude that only if all factors affecting the turn-off morphology are completely under control (and this is not the case), M67 could be used to put constraints on solar abundances.

The horizontal branch (HB) morphology of globular clusters (GCs) is most strongly influenced by metallicity. The second parameter phenomenon acknowledges that metallicity alone is not enough to describe the HB morphology of all GCs. In particular, the outer Galactic halo contains GCs with redder HBs at a given metallicity than are found inside the Solar circle. Thus, at least a second parameter is required to characterize HB morphology. Here we analyze the median color difference between the HB and the red giant branch (RGB), d(V-I), measured from HST ACS photometry of 60 GCs within ~20 kpc of the Galactic Center. Analysis of this homogeneous data set reveals that, after the influence of metallicity has been removed, the correlation between d(V-I) and age is stronger than that of any other parameter considered. Expanding the sample to include HST photometry of the 6 most distant Galactic GCs lends additional support to the correlation between d(V-I) and age. This result is robust with respect to the adopted metallicity scale and the method of age determination, but must bear the caveat that high quality, detailed abundance information is not available for a significant fraction of the sample. When a subset of GCs with similar metallicities and ages are considered, a correlation between d(V-I) and central luminosity density is exposed. With respect to the existence of GCs with anomalously red HBs at a given metallicity, we conclude that age is the second parameter and central density is most likely the third. Important problems related to HB morphology in GCs, notably multi-modal distributions and faint blue tails, remain to be explained. (Abridged)

New R-band observations of 21 local field RR Lyrae variable stars are used to explore the reliability of minimum light (V-R) colors as a tool for measuring interstellar reddening. For each star, R-band intensity mean magnitudes and light amplitudes are presented. Corresponding V-band light curves from the literature are supplemented with the new photometry, and (V-R) colors at minimum light are determined for a subset of these stars as well as for other stars in the literature. Two different definitions of minimum light color are examined, one which uses a Fourier decomposition to the V and R light curves to find (V-R) at minimum V-band light, (V-R)_{min}^F, and the other which uses the average color between the phase interval 0.5-0.8, (V-R)_{min}^{\phi(0.5-0.8)}. From 31 stars with a wide range of metallicities and pulsation periods, the mean dereddened RR Lyrae color at minimum light is (V-R)_{min,0}^F = 0.28 pm 0.02 mag and (V-R)_{min,0}^{\phi(0.5-0.8)} = 0.27 pm 0.02 mag. As was found by Guldenschuh et al. (2005) using (V-I) colors, any dependence of the star's minimum light color on metallicity or pulsation amplitude is too weak to be formally detected. We find that the intrinsic (V-R) of Galactic bulge RR Lyrae stars are similar to those found by their local counterparts and hence that Bulge RR0 Lyrae stars do not have anomalous colors as compared to the local RR Lyrae stars.

We present an analysis of the period-$V$-amplitude plane for RR0 Lyrae stars (fundamental mode pulsators) with "normal" light curves in the bulge using the MACHO bulge fields. Although bulge globular clusters have RR Lyraes that divide into two reasonable distinct groups according to the average period of the RR0 Lyraes (Oosterhoff 1939), there is no evidence of a gap between Oosterhoff I and II stars in the bulge field star sample. The majority of the bulge RR0 Lyrae field star population have a difference in period compared to the Oosterhoff I cluster M3 (Delta log P) that is shifted by about 0.02 days with regard to the Milky Way Oosterhoff I population, and the sample includes stars with Delta log P > 0.06 days, a characteristic hardly seen in Milky Way globular clusters. The metal-rich RR0 Lyrae stars in the Galactic bulge sample have Delta log P values on the other side of the spectrum as those in the metal-rich globular clusters NGC 6388 and NGC 6441. We find that the $V$-amplitude for a given period is a function of Delta log P, and not of metal abundance, similar to the result found by Clement & Shelton (1999) for RR Lyrae stars in Milky Way globular clusters. A comparative study of the bulge field stars with similar metallicities but different Oosterhoff types is carried out. Bulge field RR0 Lyrae variables with Delta log P values similar to Oosterhoff II clusters are about 0.2 mag brighter than RR0 Lyrae variables with Delta log P similar to Oosterhoff I clusters. Reliance upon a single absolute magnitude-[Fe/H] relation may be inappropriate when considering populations with different Delta log P.

We derive the distance to the northern extension of the Sagittarius (Sgr) dwarf spheroidal galaxy from 203 Sgr RR0 Lyrae stars found in the MACHO database. Their distances are determined differentially with respect to 288 Galactic Bulge RR0 Lyrae stars also found in the MACHO data. We find a distance modulus difference of 2.41 mags at $l$ = 5$^{\circ}$ and $b$ = -8$^{\circ}$ and that the extension of the Sgr galaxy towards the galactic plane is inclined toward us. Assuming $\rm R_{GC}$ = 8 kpc, this implies the distance to these stars is $(m-M)_0$ = 16.97 $\pm$ 0.07 mags, which corresponds to D = 24.8 $\pm$ 0.8 kpc. Although this estimate is smaller than previous determinations for this galaxy and agrees with previous suggestions that Sgr's body is truly closer to us, this estimate is larger than studies at comparable galactic latitudes.

Abridged: The golden age of astrophysics is upon us with both grand discoveries (extra-solar planets, dark matter, dark energy) and precision cosmology. Fundamental understanding of the working of stars and galaxies is within reach, thanks to newly available precision measurements. We highlight the importance of distances and model independent distances and masses. Distances are fundamental in astrophysics and their knowledge can change our perception of phenomena dramatically: e.g., in antiquity, the Heliocentric model was rejected because the predicted stellar parallaxes were not observed. Distance measurements are directly related to the history & fate of the uni- verse as they provide 2 of 3 methods available to date the universe. The 1st method is based on the ages of stars, which can be ascertained if their lumi- nosities (distances) are accurately known. The 2nd method relies on cosmolo- gical methods. To 1st order, the age of the universe is the inverse H_0. As stressed by the previous decadal report, "the fundamental goal of ... astrophysics is to understand how the universe ... galaxies [and] stars ... formed, how they evolved, and what their destiny will be." These questions can be answered partly by micro-arcsecond astrometry: 1) Galactic archeology: a detailed reconstruction of the formation history of the Milky Way and other Local Group galaxies, 2) the very oldest stars in the Milky Way and the age of the Universe, and 3) H_0 and concordance cosmology. These goals are achievable by combining muas-arcsecond astrometry from the proposed SIM-Lite mission supplemented with ground-based spectroscopy. The results of our proposed project will force the biggest reassessment of stellar astrophysics in 50 years, which will affect most branches of astrophysics.

The ACS Survey of Galactic Globular Clusters is a Hubble Space Telescope (HST) Treasury program designed to provide a new large, deep and homogeneous photometric database. Based on observations from this program, we have measured precise relative ages for a sample of 64 Galactic globular clusters by comparing the relative position of the clusters' main sequence turn offs, using main-sequence fitting to cross-compare clusters within the sample. This method provides relative ages to a formal precision of 2-7%. We demonstrate that the calculated relative ages are independent of the choice of theoretical model. We find that the Galactic globular cluster sample can be divided into two groups -- a population of old clusters with an age dispersion of ~5% and no age-metallicity relation, and a group of younger clusters with an age-metallicity relation similar to that of the globular clusters associated with the Sagittarius dwarf galaxy. These results are consistent with the Milky Way halo having formed in two phases or processes. The first one would be compatible with a rapid (<0.8 Gyr) assembling process of the halo, in which the clusters in the old group were formed. Also these clusters could have been formed before reionization in dwarf galaxies that would later merge to build the Milky Way halo as predicted by Lambda-CDM cosmology. However, the galactocentric metallicity gradient shown by these clusters seems difficult to reconcile with the latter. As for the younger clusters, it is very tempting to argue that their origin is related to their formation within Milky Way satellite galaxies that were later accreted, but the origin of the age-metallicity relation remains unclear.

The first paper in this series explored the effects of altering the chemical mixture of the stellar population on an element by element basis on stellar evolutionary tracks and isochrones to the end of the red giant branch. This paper extends the discussion by incorporating the fully consistent synthetic stellar spectra with those isochrone models in predicting integrated colors, Lick indices, and synthetic spectra. Older populations display element ratio effects in their spectra at higher amplitude than younger populations. In addition, spectral effects in the photospheres of stars tend to dominate over effects from isochrone temperatures and lifetimes, but, further, the isochrone-based effects that are present tend to fall along the age-metallicity degeneracy vector, while the direct stellar spectral effects usually show considerable orthogonality.

We present metallicities of 2690 RR0 Lyrae stars observed toward the MACHO Survey fields in the Galactic Bulge. These [Fe/H] values are based upon an empirically calibrated relationship that uses the Fourier coefficients of the light curve and are accurate to ~0.2 dex. The majority of the RR0 Lyrae stars in our sample are located in the Galactic Bulge, but 255 RR0 stars are associated with Sagittarius dwarf galaxy. We find that the RR0 Lyrae stars that belong to the Galactic Bulge have average metallicities [Fe/H] = -1.25, with a broad metallicity range from [Fe/H] = -2.26 to -0.15. The RR0 stars from Sagittarius dwarf galaxy have lower average metallicity of [Fe/H] = -1.55 \pm 0.02, with an intrinsic dispersion of 0.25 dex, similar to that in the bulge. A correlation between metallicity and galactocentric distance is found, in a sense that for the metal-poor RR0 Lyrae stars ([Fe/H]< -1.5 dex), the closer a star is to the Galactic center, on average, the more metal rich it is. However, for the metal-rich RR0 Lyrae stars ([Fe/H] > -1.2 dex), this trend is reversed. Using mean magnitudes of MACHO RR Lyrae stars, we searched for the evidence of the Galactic bar, and found marginal evidence of a bar. The absence of a strong bar indicates that the RR Lyrae in the bulge represent a different population than the majority of the bulge stars, which are metal rich and are part of a bar.

The ever-expanding depth and quality of photometric and spectroscopic observations of stellar populations increase the need for theoretical models in regions of age-composition parameter space that are largely unexplored at present. Stellar evolution models that employ the most advanced physics and cover a wide range of compositions are needed to extract the most information from current observations of both resolved and unresolved stellar populations. The Dartmouth Stellar Evolution Database is a collection of stellar evolution tracks and isochrones that spans a range of [Fe/H] from -2.5 to +0.5, [alpha/Fe] from -0.2 to +0.8 (for [Fe/H] <=0) or +0.2 (for [Fe/H] >0), and initial He mass fractions from Y=0.245 to 0.40. Stellar evolution tracks were computed for masses between 0.1 and 4 Msun, allowing isochrones to be generated for ages as young as 250 Myr. For the range in masses where the core He flash occurs, separate He-burning tracks were computed starting from the zero age horizontal branch. The tracks and isochrones have been transformed to the observational plane in a variety of photometric systems including standard UBV(RI)c, Stromgren uvby, SDSS ugriz, 2MASS JHKs, and HST ACS-WFC and WFPC2. The Dartmouth Stellar Evolution Database is accessible through a website at http://stellar.dartmouth.edu/~models/ where all tracks, isochrones, and additional files can be downloaded.

The ACS Survey of Globular Clusters has used HST's Wide-Field Channel to obtain uniform imaging of 65 of the nearest globular clusters to provide an extensive homogeneous dataset for a broad range of scientific investigations. The survey goals required not only a uniform observing strategy, but also a uniform reduction strategy. To this end, we designed a sophisticated software program to process the cluster data in an automated way. The program identifies stars simultaneously in the multiple dithered exposures for each cluster and measures them using the best available PSF models. We describe here in detail the program's rationale, algorithms, and output. The routine was also designed to perform artificial-star tests, and we run a standard set of ~10^5 tests for each cluster in the survey. The catalog described here will be exploited in a number of upcoming papers and will eventually be made available to the public via the world-wide web.

Searches for extrasolar planets have uncovered an astonishing diversity of planetary systems, yet the frequency of solar system analogs remains unknown. The gravitational microlensing planet search method is potentially sensitive to multiple-planet systems containing analogs of all the solar system planets except Mercury. We report the detection of a multiple-planet system with microlensing. We identify two planets with masses of ~0.71 and ~0.27 times the mass of Jupiter and orbital separations of ~2.3 and ~4.6 astronomical units orbiting a primary star of mass ~0.50 solar masses at a distance of ~1.5 kiloparsecs. This system resembles a scaled version of our solar system in that the mass ratio, separation ratio, and equilibrium temperatures of the planets are similar to those of Jupiter and Saturn. These planets could not have been detected with other techniques; their discovery from only six confirmed microlensing planet detections suggests that solar system analogs may be common.

We present mean reddenings toward 3525 RR0 Lyrae stars from the Galactic bulge fields of the MACHO Survey. These reddenings are determined using the color at minimum $V$-band light of the RR0 Lyrae stars themselves and are found to be in general agreement with extinction estimates at the same location obtained from other methods. Using 3256 stars located in the Galactic Bulge, we derive the selective extinction coefficient $R_{V,VR}=A_V/E(V-R) = 4.3 \pm 0.2$. This value is what is expected for a standard extinction law with $R_{V,BV} = 3.1 \pm 0.3$.

Precision astrometry at microarcsecond accuracy has application to a wide range of astrophysical problems. This paper is a study of the science questions that can be addressed using an instrument that delivers parallaxes at about 4 microarcsec on targets as faint as V = 20, differential accuracy of 0.6 microarcsec on bright targets, and with flexible scheduling. The science topics are drawn primarily from the Team Key Projects, selected in 2000, for the Space Interferometry Mission PlanetQuest (SIM PlanetQuest). We use the capabilities of this mission to illustrate the importance of the next level of astrometric precision in modern astrophysics. SIM PlanetQuest is currently in the detailed design phase, having completed all of the enabling technologies needed for the flight instrument in 2005. It will be the first space-based long baseline Michelson interferometer designed for precision astrometry. SIM will contribute strongly to many astronomical fields including stellar and galactic astrophysics, planetary systems around nearby stars, and the study of quasar and AGN nuclei. SIM will search for planets with masses as small as an Earth orbiting in the `habitable zone' around the nearest stars using differential astrometry, and could discover many dozen if Earth-like planets are common. It will be the most capable instrument for detecting planets around young stars, thereby providing insights into how planetary systems are born and how they evolve with time. SIM will observe significant numbers of very high- and low-mass stars, providing stellar masses to 1%, the accuracy needed to challenge physical models. Using precision proper motion measurements, SIM will probe the galactic mass distribution and the formation and evolution of the Galactic halo. (abridged)

Photometry with HST's ACS reveals that the subgiant branch (SGB) of the globular cluster NGC 1851 splits into two well-defined branches. If the split is due only to an age effect, the two SGBs would imply two star formation episodes separated by $\sim$ 1 Gyr. We discuss other anomalies in NGC 1851 which could be interpreted in terms of a double stellar population. Finally, we compare the case of NGC 1851 with the other two globulars known to host multiple stellar populations, and show that all three clusters differ in several important respects.

We present new Hubble Space Telescope photometry of the massive globular cluster M54 (NGC 6715) and the superposed core of the tidally disrupted Sagittarius (Sgr) dSph galaxy as part of the ACS Survey of Galactic Globular Clusters. Our deep (F606W~26.5), high-precision photometry yields an unprecedentedly detailed color-magnitude diagram showing the extended blue horizontal branch and multiple main sequences of the M54+Sgr system. The distance and reddening to M54 are revised usingboth isochrone and main-sequence fitting to (m-M)_0=17.27 and E(B-V)=0.15. Preliminary assessment finds the M54+Sgr field to be dominated by the old metal-poor populations of Sgr and the globular cluster. Multiple turnoffs indicate the presence of at least two intermediate-aged star formation epochs with 4 and 6 Gyr ages and [Fe/H]=-0.4 to -0.6. We also clearly show, for the first time, a prominent, 2.3 Gyr old Sgr population of near-solar abundance. A trace population of even younger (0.1-0.8 Gyr old), more metal-rich ([Fe/H]\sim0.6) stars is also indicated. The Sgr age-metallicity relation is consistent with a closed-box model and multiple (4-5) star formation bursts over the entire life of the satellite, including the time since Sgr began disrupting.

The ACS Survey of Galactic Globular Clusters, an HST Treasury Project, will deliver high quality, homogeneous photometry of 65 globular clusters. This paper introduces a new collection of stellar evolution tracks and isochrones suitable for analyzing the ACS Survey data. Stellar evolution models were computed at [Fe/H]= -2.5, -2.0, -1.5, -1.0, -0.5, and 0; [alpha/Fe]= -0.2, 0, 0.2, 0.4, 0.6, and 0.8; and three initial He abundances for masses from 0.1 to 1.8 Msun and ages from 2 to 15 Gyr. Each isochrone spans a wide range in luminosity from Mv~14 up to the tip of the red giant branch. These are complemented by a set of He-burning tracks that extend from the zero age horizontal branch to the onset of thermal pulsations on the asymptotic giant branch. In addition, a set of computer programs are provided that make it possible to interpolate the isochrones in [Fe/H], generate luminosity functions from the isochrones, and create synthetic horizontal branch models. The tracks and isochrones have been converted to the observational plane with two different color-Teff transformations, one synthetic and one semi-empirical, in ground-based B, V, and I, and F606W and F814W for both ACS-WFC and WFPC2 systems. All models and programs presented in this paper are available from http://stellar.dartmouth.edu/~models/

Integrated light from distant galaxies is often compared to stellar population models via the equivalent widths of spectral features--spectral indices--whose strengths rely on the abundances of one or more elements. Such comparisons hinge not only on the overall metal abundance but also on relative abundances. Studies have examined the influence of individual elements on synthetic spectra but little has been done to address similar issues in the stellar evolution models that underlie most stellar population models. Stellar evolution models will primarily be influenced by changes in opacities. In order to explore this issue in detail, twelve sets of stellar evolution tracks and isochrones have been created at constant heavy element mass fraction Z that self-consistently account for varying heavy element mixtures. These sets include scaled-solar, alpha-enhanced, and individual cases where the elements C, N, O, Ne, Mg, Si, S, Ca, Ti, and Fe have been enhanced above their scaled-solar values. The variations that arise between scaled-solar and the other cases are examined with respect to the H-R diagram and main sequence lifetimes.

The past two Decadal Surveys in Astronomy and Astrophysics recommended the completion of a space-based interferometry mission, known today as SIM PlanetQuest, for its unique ability to detect and characterize nearby rocky planets (Bahcall 1991, McKee & Taylor 2001), as well as contributions to a broad range of problems in astrophysics. Numerous committees of the National Research Council as well as NASA Roadmaps have similarly highlighted SIM as the one technology that offers detection and characterization of rocky planets around nearby stars and which is technically ready. To date, SIM remains the only program with the capability of detecting and confirming rocky planets in the habitable zones of nearby solar-type stars. Moreover, SIM measures masses and three-dimensional orbits of habitable planets around nearby stars (within 25 pc); these are the only stars for which follow-up by other techniques is feasible, such as space-based spectroscopy, ground-based interferometry, and of course TPF.

We present new BVI ground-based photometry and VI space-based photometry for the globular cluster M92 (NGC 6341) and examine luminosity functions in B, V, and I containing over 50,000 stars ranging from the tip of the red giant branch to several magnitudes below the main sequence turn off. Once corrected for completeness, the observed luminosity functions agree very well with theoretical models and do not show stellar excesses in any region of the luminosity function. Using reduced chi squared fitting, the new M92 luminosity function is shown to be an excellent match to the previously published luminosity function for M30. These points combine to establish that the "subgiant excess" found in previously published luminosity functions of Galactic globular clusters are due to deficiencies in the stellar models used at that time. Using up to date stellar models results in good agreement between observations and theory. Several statistical methods are presented to best determine the age of M92. These methods prove to be insensitive to the exact choice of metallicity within the published range. Using [Fe/H]=-2.17 to match recent studies we find an age of 14.2 plus or minus 1.2 Gyr for the cluster.

We are developing a brand new stellar population models with flexible chemistry (isochrones plus stellar colors and spectra) in order to set a new standard of completeness and excellence. Here we present preliminary results to assess the effects of stellar evolution models and stellar model atmosphere to the well-known Lick indices at constant heavy element mass fraction Z that self-consistently account for varying heavy element mixtures. We have enhanced chemical elements one by one. Our ultimate goal is to demonstrate 10% absolute mean ages for a sample of local galaxies derived from an integrated light spectrum.

We present the first results of a large ACS Survey of Galactic globular clusters. This Hubble Space Telescope (HST) Treasury project is designed to obtain photometry with S/N > ~10 for main sequence stars with masses > ~0.2Msun in a sample of globulars using the Advanced Camera for Surveys (ACS) Wide Field Channel. Here we focus on clusters without previous HST imaging data. These include NGC 5466, 6779, 5053, 6144, Palomar 2, E 3, Lynga 7, Palomar 1, and NGC 6366. Our CMDs extend reliably from the horizontal branch to as much as seven magnitudes fainter than the main sequence turnoff and represent the deepest CMDs published to-date for these clusters. Using fiducial sequences for three standard clusters (M92, NGC 6752, and 47 Tuc) with well-known metallicities and distances, we perform main sequence fitting on the target clusters in order to obtain estimates of their distances and reddenings. These comparisons along with fitting the cluster main sequences to theoretical isochrones yield ages for the target clusters. We find that the majority of the clusters have ages that are consistent with the standard clusters at their metallicities. The exceptions are E 3 which appears ~2 Gyr younger than 47 Tuc, and Pal 1, which could be as much as 8 Gyr younger than 47 Tuc.

We present HI observations of 54 galaxies in the Pegasus cluster. The observations include single dish HI measurements, obtained with the Arecibo telescope for all 54 galaxies in the sample, as well as HI images, obtained with the VLA for 10 of these. The Arecibo profiles reveal an overall HI deficiency in the cluster, with ~40% of the galaxies in the core of the cluster showing modest deficiencies of typically a factor of 2-3. The HI morphology of some galaxies shows that the HI disk is smaller than the optical disk and slightly offset from the stars. We find a correlation between HI deficiency and the ratio of the HI disk size to optical disk size. More HI deficient galaxies have relatively smaller HI disks, a configuration that is usually attributed to an interaction between the ISM of the galaxy and the hot ICM. Such a result is surprising since the Pegasus cluster has a low level of X-ray emission, and a low velocity dispersion. The low velocity dispersion, coupled with the lack of a dense hot ICM indicate that ram pressure stripping should not play a significant role in this environment. In addition, two of the galaxies, NGC7604 and NGC7648, are morphologically peculiar. Their peculiarities indicate contradictory scenarios of what is triggering their unusual star formation. Halpha imaging, along with long-slit spectroscopy of NGC7648 reveal morphological features which point to a recent tidal interaction. On the other hand, Halpha imaging,along with VLA HI mapping, of NGC7604 reveal morphological features suggestive of a ram pressure event. Our data indicate that ISM-ICM interactions may play a role in a wider variety of environments than suggested by simple ram pressure arguments.

SDSSJ1049+5103, commonly known as Willman 1, is an extremely low-luminosity Milky Way companion whose properties are intermediate between those of globular clusters and dwarf spheroidals. In this paper, we use deep photometry to show that this object is old and moderately metal-poor, has a distance of 38 +/- 7 kpc, has an M_V of -2.5 mag, and has a half-light radius of 21 +/- 7 pc, consistent with previous estimates. The spatial distribution of Willman 1's main sequence stars shows 1) its total spatial extent exceeds its tidal radius for a range of assumptions about its total mass and its orbit and 2) the presence of prominent multi-directional stellar tails. The tidal interactions causing these tail features may explain the large physical size of Willman 1 relative to low-luminosity globular clusters. It is the most distant Galactic object yet known to display prominent tails, and is the only distant satellite to display multi-directional tails. Although we cannot at present determine the cause of this unusual morphology, preliminary comparisons between the morphology of Willman 1 and published simulations suggest that it may be near the apocenter of its orbit, or that it may have interacted with another halo object. We find a significant difference between the luminosity functions of stars in the center and in the tails of Willman 1, strongly suggesting mass segregation much like that seen in Palomar 5. Although Willman 1 has more pronounced tidal tails than most confirmed Milky Way dwarf galaxies, because of its very low stellar mass we cannot at present rule out the possibility that it has a dark matter halo. (Abridged)

We present new photometry for the Galactic thick disk globular cluster Palomar 11 extending well past the main sequence turn-off in the V and I bands. This photometry shows noticeable red giant and subgiant branches. The difference in magnitude between the red horizontal branch (red clump) and the subgiant branch is used to determine that Palomar 11 has an age of 10.4+/-0.5 Gyr. The red clump is used to derive a distance d_\sun=14.3+/-0.4 kpc, and a mean cluster reddening of E(V-I)=0.40+/-0.03. There is differential reddening across the cluster, of order \delta E(V-I)~0.07. The colour magnitude diagram of Palomar 11 is virtually identically to that of the thick disk globular cluster NGC 5927, implying that these two clusters have a similar age and metallicity. Palomar 11 has a slightly redder red giant branch than 47 Tuc, implying that Palomar 11 is 0.15 dex more metal-rich, or 1 Gyr older than 47 Tuc. Ca II triplet observations (Rutledge et al. 1997) favour the hypothesis that Palomar 11 is the same age as 47 Tuc, but slightly more metal-rich.

Stellar models which incorporate simple diffusion or shear induced mixing are used to describe canonical extra mixing in low mass red giants of low and solar metallicity. These models are able to simultaneously explain the observed Li and CN abundance changes along upper red giant branch (RGB) in field low-metallicity stars and match photometry, rotation and carbon isotopic ratios for stars in the old open cluster M67. The shear mixing model requires that main sequence (MS) progenitors of upper RGB stars possessed rapidly rotating radiative cores and that specific angular momentum was conserved in each of their mass shells during their evolution. We surmise that solar-type stars will not experience canonical extra mixing on the RGB because their more efficient MS spin-down resulted in solid-body rotation, as revealed by helioseismological data for the Sun. Thus, RGB stars in the old, high metallicity cluster NGC 6791 should show no evidence for mixing in their carbon isotopic ratios. We develop the idea that canonical extra mixing in a giant component of a binary system may be switched to its enhanced mode with much faster and somewhat deeper mixing as a result of the giant's tidal spin-up. This scenario can explain photometric and composition peculiarities of RS CVn binaries. The tidally enforced enhanced extra mixing might contribute to the star-to-star abundance variations of O, Na and Al in globular clusters. This idea may be tested with observations of carbon isotopic ratios and CN abundances in RS CVn binaries.

A Monte Carlo simulation exploring uncertainties in standard stellar evolution theory on the red giant branch of metal-poor globular clusters has been conducted. Confidence limits are derived on the absolute V-band magnitude of the bump in the red giant branch luminosity function (M_v,b) and the excess number of stars in thebump, R_b. The analysis takes into account uncertainties in the primordial helium abundance, abundance of alpha-capture elements, radiative and conductive opacities, nuclear reaction rates, neutrino energy losses, the treatments of diffusion and convection, the surface boundary conditions, and color transformations. The uncertainty in theoretical values for the red giant bump magnitude varies with metallicity between +0.13/-0.12 mag at [Fe/H] = -2.4 and +0.23/-0.21 mag at [Fe/H] = -1.0$. The dominant sources of uncertainty are the abundance of the alpha-capture elements, the mixing length, and the low-temperature opacities. The theoretical values of M_v,b are in good agreement with observations. The uncertainty in the theoretical value of R_b is +/-0.01 at all metallicities studied. The dominant sources of uncertainty are the abundance of the alpha-capture elements, the mixing length, and the high-temperature opacities. The median value of R_b varies from 0.44 at [Fe/H] = -2.4$ to 0.50 at [Fe/H] = -1.0. These theoretical values for R_b are in agreement with observations.

Photometric UBVI CCD photometry is presented for NGC 188 and Berkeley 17. Color-magnitude diagrams (CMDs) are constructed and reach well past the main-sequence turn-off for both clusters. Cluster ages are determined by means of isochrone fitting to the cluster CMDs. These fits are constrained to agree with spectroscopic metallicity and reddening estimates. Cluster ages are determined to be 7.0+/-0.5 Gyr for NGC 188, and 10.0+/- 1.0 Gyr for Berkeley 17, where the errors refer to uncertainties in the relative age determinations. These ages are compared to the ages of relatively metal-rich inner halo/thick disk globular clusters and other old open clusters. Berkeley 17 and NGC 6791 are the oldest open clusters with an age of 10 Gyr. They are 2 Gyr younger than the thick disk globular clusters. These results confirm the status of Berkeley 17 as one of the oldest known open cluster in the Milky Way, and its age provides a lower limit to the age of the Galactic disk.

The first CCD UBV(RI) photometric study in the area of the doubtful open cluster NGC 2129 is presented. Photometry of a field offset 15 arcmin northward is also provided, to probe the Galactic disk population toward the cluster. Using star counts, proper motions from the UCAC2 catalog, colour-magnitude and colour-colour diagrams we demonstrate that NGC 2129 is a young open cluster. The cluster radius is 2.5 arcmin, and across this region we find evidence of significant differential reddening, although the reddening law seems to be normal toward its direction. Updated estimates of the cluster fundamental parameters are provided. The mean reddening is found to be E(B-V)=0.80$\pm$0.08 and the distance modulus is $(m-M)_o$= 11.70$\pm0.30$. Hence, NGC 2129 is located at 2.2$\pm$0.2 kpc from the Sun inside the Local spiral arm. The age derived from 37 photometrically selected members is estimated to be approximately 10 million years. These stars are used to provide new estimates of the cluster absolute proper motion components.

An approach is introduced for incorporating the concept of stellar pollution into stellar evolution models. The approach involves enhancing the metal content of the surface layers of stellar models. In addition, the surface layers of stars in the mass range of 0.5-2.0 Solar masses are mixed to an artificial depth motivated by observations of lithium abundance. The behavior of polluted stellar evolution models is explored assuming the pollution occurs after the star has left the fully convective pre main sequence phase. Stellar models polluted with a few Earth masses of iron are significantly hotter than stars of the same mass with an equivalent bulk metallicity. Polluted stellar evolution models can successfully reproduce the metal-rich, parent star tau Bootis and suggest a slightly lower mass than standard evolution models. Finally, the possibility that stars in the Hyades open cluster have accreted an average of 0.5 Earth masses of iron is explored. The results indicate that it is not possible to rule out stellar pollution on this scale from the scatter of Hyades stars on a color-magnitude diagram. The small amount of scatter in the observational data set does rule out pollution on the order of 1.5 Earth masses of iron. Pollution effects at the low level of 0.5 Earth masses of iron do not produce substantial changes in a star's evolution.

The Hyades open cluster presents a unique laboratory for planet formation and stellar pollution studies because all of the stars have essentially the same age and were born from the same cloud of gas. Furthermore, with an age of roughly 650 Myr most of the intermediate and low mass stars are on the main sequence. Given these assumptions, the accretion of metal rich material onto the surface of a star during and shortly after the formation of planetary systems should be evident via the enhanced metallicity of the star. Building on previous work, stellar evolution models which include the effects of stellar pollution are applied to the Hyades. The results of several Monte Carlo simulations, in which the amount of accreted material is drawn at random from a Gaussian distribution with standard deviation equal to half the mean, are presented. An effective temperature-[Fe/H] relation is produced and compared to recent observations. The theoretical predictions presented in this letter will be useful in future searches for evidence of stellar pollution due to planet formation. It is concluded that stellar pollution effects at the mean level of >=2 Earth masses of iron are ruled out by current observational data.

The theoretical uncertainties in the calibration of the relationship between the subgiant mass and age in metal-poor stars are investigated using a Monte Carlo approach. Assuming that the mass and iron abundance of a subgiant star are known exactly, uncertainties in the input physics used to construct stellar evolution models and isochrones lead to a Gaussian 1-sigma uncertainty of +/-2.9% in the derived ages. The theoretical error budget is dominated by the uncertainties in the calculated opacities. Observations of detached double lined eclipsing binary OGLEGC-17 in the globular cluster Omega Cen have found that the primary is on the subgiant branch with a mass of M = 0.809+/-0.012 M_sun and [Fe/H]= -2.29+/-0.15 (Kaluzny et al. 2001). Combining the theoretical uncertainties with the observational errors leads to an age for OGLEGC-17 of 11.10+/-0.67 Gyr. The one-sided, 95% lower limit to the age of OGLEGC-17 is 10.06 Gyr, while the one-sided, 95% upper limit is 12.27 Gyr.

New estimates of globular cluster distances, combined with revised ranges for input parameters in stellar evolution codes and recent estimates of the earliest redshift of cluster formation allow us to derive a new 95% confidence level lower limit on the age of the Universe of 11 Gyr. This is now definitively inconsistent with the expansion age for a flat Universe for the currently allowed range of the Hubble constant unless the cosmic equation of state is dominated by a component that violates the strong energy condition. This solidifies the case for a dark energy-dominated universe, complementing supernova data and direct measurements of the geometry and matter density in the Universe. The best-fit age is consistent with a cosmological constant-dominated (w=pressure/energy density = -1) universe. For the Hubble Key project best fit value of the Hubble Constant our age limits yields the constraints w < -0.4 and Omega_matter < 0.38 at the 68 % confidence level, and w < -0.26 and Omega_matter < 0.58 at the 95 % confidence level.

We present high-precision V, B-V color-magnitude diagrams (CMDs) for the classic second parameter globular clusters M3 and M13 from wide-field deep CCD photometry. The data for the two clusters were obtained during the same photometric nights with the same instrument, allowing us to determine accurate relative ages. Based on a differential comparison of the CMDs using the Delta (B-V) method, an age difference of 1.7 +/- 0.7 Gyr is obtained between these two clusters. We compare this result with our updated horizontal-branch (HB) population models, which confirm that the observed age difference can produce the difference in HB morphology between the clusters. This provides further evidence that age is the dominant second parameter that influences HB morphology.

Stellar evolution models which include the effect of helium and heavy element diffusion have been calculated for initial iron abundances of [Fe/H] = -2.3, -2.1, -1.9 and -1.7. These models were calculated for a large variety of masses and three separate mixing lengths, \alpha = 1.50, 1.75 and 2.00 (with \alpha = 1.75 being the solar calibrated mixing length). The change in the surface iron abundance for stars of different masses was determined for the ages 11, 13 and 15 Gyr. Iron settles out of the surface convection zone on the main sequence; this iron is dredged back up when the convection zone deepens on the giant branch. In all cases, the surface [Fe/H] abundance in the turn-off stars was at least 0.28 dex lower than the surface [Fe/H] abundance in giant-branch stars of the same age. However, recent high dispersion spectra of stars in the globular cluster NGC 6397 found that the turn-off and giant branch stars had identical (within a few percent) iron abundances of [Fe/H] =-2.03 (Gratton et al. 2001). These observations prove that heavy element diffusion must be inhibited in the surface layers of metal-poor stars. When diffusion is inhibited in the outer layers of a stellar model, the predicted temperatures of the models are similar to models evolved without diffusion, while the predicted lifetimes are similar to stars in which diffusion is not inhibited. Isochrones constructed from the models in which diffusion is inhibited fall half-way between isochrones without diffusion, and isochrones with full diffusion. As a result, absolute globular cluster ages which are based upon the absolute magnitude of the turn-off are 4% larger than ages inferred from full diffusion isochrones, and 4% smaller than ages inferred from no diffusion isochrones.

We compare the metallicities of stars with radial velocity planets to the metallicity of a sample of field dwarfs. We confirm recent work indicating that the stars-with-planet sample as a whole is iron rich. However, the lowest mass stars tend to be iron poor, with several having [Fe/H]<-0.2, demonstrating that high metallicity is not required for the formation of short period Jupiter-mass planets. We show that the average [Fe/H] increases with increasing stellar mass (for masses below 1.25 solar masses) in both samples, but that the increase is much more rapid in the stars-with-planet sample. The variation of metallicity with stellar age also differs between the two samples. We examine possible selection effects related to variations in the sensitivity of radial velocity surveys with stellar mass and metallicity, and identify a color cutoff (B-V>0.48) that contributes to but does not explain the mass-metallicity trend in the stars-with-planets sample. We use Monte Carlo models to show that adding an average of 6.5 Earth masses of iron to each star can explain both the mass-metallicity and the age-metallicity relations of the stars-with-planets sample. However, for at least one star, HD 38529, there is good evidence that the bulk metallicity is high. We conclude that the observed metallicities and metallicity trends are the result of the interaction of three effects; accretion of about 6 Earth masses of iron rich material, selection effects, and in some cases, high intrinsic metallicity.

We present our recent revision of model constructions for the horizontal-branch (HB) morphology of globular clusters, which suggests the HB morphology is more sensitive to age compared to our earlier models. We also present our high precision CCD photometry for the classic second parameter pair M3 and M13. The relative age dating based on this photometry indicates that M13 is indeed older than M3 by 1.7 Gyr. This is consistent with the age difference predicted from our new models, which provides a further support that the HB morphology is a reliable age indicator in most population II stellar systems.

We study spectroscopically determined iron abundances of 642 solar-type stars to search for the signature of accreted iron-rich material. We find that the metallicity [Fe/H] of a subset of 466 main sequence stars, when plotted as a function of stellar mass, mimics the pattern seen in lithium abundances in open clusters. Using Monte Carlo models we find that, on average, these stars have accreted about 0.4 Earth masses of iron while on the main sequence. A much smaller sample of 19 stars in the Hertzsprung gap, which are slightly evolved and whose convection zones are significantly more massive, have lower average [Fe/H], and their metallicity shows no clear variation with stellar mass. These findings suggest that terrestrial-type material is common around solar type stars.

We present long-slit spectroscopy, B and R bandpass imaging, and 21 cm observations of a sample of early-type galaxies in nearby clusters which are known to be either in a star-forming phase or to have had star formation which recently terminated. From the long-slit spectra, obtained with the Blanco 4-m telescope, we find that emission lines in the star-forming cluster galaxies are significantly more centrally concentrated than in a sample of field galaxies. The broadband imaging reveals that two currently star-forming early-type galaxies in the Pegasus I cluster have blue nuclei, again indicating that recent star formation has been concentrated. In contrast, the two galaxies for which star formation has already ended show no central color gradient. The Pegasus I galaxy with the most evident signs of ongoing star formation (NGC7648), exhibits signatures of a tidal encounter. Neutral hydrogen observations of that galaxy with the Arecibo radiotelescope reveal the presence of ~4 x 10^8 solar masses of HI. Arecibo observations of other current or recent star-forming early-type galaxies in Pegasus I indicate smaller amounts of gas in one of them, and only upper limits in others.

HST (V,I) photometry has been obtained for the inner halo globular cluster NGC 6652. The photometry reaches approximately 4 mag below the turn-off and includes a well populated horizontal branch. This cluster is located close to the Galactic center at a galactocentric distance of approximately 2.0 kpc with a reddening of E(V-I) = 0.15 +/- 0.02 and has a metallicity of [Fe/H] approximately -0.85. Based upon Delta(V) between the point on the sub-giant branch which is 0.05 mag redder than the turn-off and the horizontal branch, NGC 6652 is 11.7 +/- 1.6 Gyr old. Using this same Delta(V), precise differential ages for 47 Tuc (a thick disk globular), M107 and NGC 1851 (both halo clusters) were obtained. NGC 6652 appears to be the same age as 47 Tuc and NGC 1851 (within +/- 1.2 Gyr), while there is a slight suggestion that M107 is older than NGC 6652 by 2.3 +/- 1.5 Gyr. As this is a less than 2-sigma result, this issue needs to be investigated further before a definitive statement regarding the relative age of M107 and NGC 6652 may be made.

We determine the relative ages of the open cluster NGC 188 and selected Hipparcos field stars by isochrone fitting, and compare them to the age of the thick disk globular cluster 47 Tuc. The best fit age for NGC 188 was determined to be $6.5 \pm 1.0$ Gyr. The solar metallicity Hipparcos field stars yielded a slightly older thin disk age, $7.5 \pm 0.7$ Gyr. Two slightly metal-poor ($\feh = -0.22$) field stars whose kinematic and orbital parameters indicate that they are members of the thin disk were found to have an age of $9.7\pm 0.6$ Gyr. The age for 47 Tuc was determined to be $12.5 \pm 1.5$ Gyr. All errors are internal errors due to the uncertainty in the values of metallicity and reddening. Thus, the oldest stars dated in the thin disk are found to be $2.8\pm 1.6$ Gyr younger than 47 Tuc. Furthermore, as discussed by \citet{Chb99} 47 Tuc has a similar age to three globular clusters located in the inner part of the Galactic halo, implying that star formation in the thin disk started within $2.8\pm 1.6$ Gyr of star formation in the halo.

A grid of stellar evolution models for Procyon A has been calculated. These models include the best physics available to us (including the latest opacities and equation of state) and are based on the revised astrometric mass of Girard et al (1996). Models were calculated with helium diffusion and with the combined effects of helium and heavy element diffusion. Oscillation frequencies for l=0,1,2 and 3 p-modes and the characteristic period spacing for the g-modes were calculated for these models. We find that g-modes are sensitive to model parameters which effect the structure of the core, such as convective core overshoot, the heavy element abundance and the evolutionary state (main sequence or shell hydrogen burning) of Procyon A. The p-modes are relatively insensitive to the details of the physics used to model Procyon A, and only depend on the evolutionary state of Procyon A. Hence, observations of p-mode frequencies on Procyon A will serve as a robust test of stellar evolution models.

An extensive grid of metal-rich isochrones utilizing the latest available input physics has been calculated for comparison with the old, metal-rich open cluster NGC 6791. The isochrones have been simultaneously fit to BV and VI color magnitude diagrams, with the same composition, reddening and distance modulus required for both colors. Our best fitting isochrone assumes [Fe/H] = +0.4, scaled solar abundance ratios, and dY/dZ = 2 (Y = 0.31), yielding an excellent fit to the data at all points along the major sequences. The resulting age is 8 Gyr, with E(B-V) = 0.10 and (m-M)_v = 13.42. The derived cluster parameters are fairly robust to variations in the isochrone [Fe/H] and helium abundances. All of the acceptable fits indicate that 0.07 < E(B-V) < 0.14$, 13.29 < (m-M)_v < 13.46, and that NGC 6791 has an age of 8.0+/- 0.5 Gyr. The fits also suggest that dY/dZ lies between 1 and 3. A metallicity as low as solar is clearly ruled out, as is dY/dZ = 0. Comparison with previous isochrone studies indicates that the derived reddening is primarily due to our use of the most recent color transformations, whereas the age depends upon both the colors and the input physics. Our isochrones provide an excellent fit to the Hyades zero-age main sequence as determined by Hipparcos, providing evidence that our derived reddening and distance modulus are reliable.

The present status of the distance scale to Galactic globular clusters is reviewed. Six distance determination techniques which are deemed to be most reliable are discussed in depth. These different techniques are used to calibrate the absolute magnitude of the RR Lyrae stars. The various calibrations fall into three groups. Main sequence fitting using Hipparcos parallaxes, theoretical HB models and the RR Lyrae in the LMC all favor a bright calibration, implying a `long' globular cluster distance scale. White dwarf fitting and the astrometric distances yield a somewhat fainter RR Lyrae calibration, while the statistical parallax solution yields faint RR Lyrae stars implying a `short' distance scale to globular clusters. Various secondary distance indicators discussed all favor the long distance scale. The `long' and `short' distance scales differ by (0.31+/-0.16) mag. Averaging together all of the different distance determinations yields Mv(RR) = (0.23+/-0.04)([Fe/H] + 1.6) + (0.56+/-0.12) mag.

A minimum age of the universe can be estimated directly by determining the age of the oldest objects in the our Galaxy. These objects are the metal-poor stars in the halo of the Milky Way. Recent work on nucleochronology finds that the oldest stars are 15.2+/-3.7 Gyr old. White dwarf cooling curves have found a minimum age for the oldest stars of 8 Gyr. Currently, the best estimate for the age of the oldest stars is based upon the absolute magnitude of the main sequence turn-off in globular clusters. The oldest globular clusters are 11.5+/-1.3 Gyr old, implying a minimum age of the universe of t_universe > 9.5 Gyr (95% confidence level).

Lithium is an excellent tracer of mixing in stars as it is destroyed (by nuclear reactions) at a temperature around $\sim 2.5\times 10^6$ K. The lithium destruction zone is typically located in the radiative region of a star. If the radiative regions are stable, the observed surface value of lithium should remain constant with time. However, comparison of the meteoritic and photospheric Li abundances in the Sun indicate that the surface abundance of Li in the Sun has been depleted by more than two orders of magnitude. This is not predicted by solar models and is a long standing problem. Observations of Li in open clusters indicate that Li depletion is occurring on the main sequence. Furthermore, there is now compelling observational evidence that a spread of lithium abundances is present in nearly identical stars. This suggests that some transport process is occurring in stellar radiative regions. Helioseismic inversions support this conclusion, for they suggest that standard solar models need to be modified below the base of the convection zone. There are a number of possible theoretical explanations for this transport process. The relation between Li abundances, rotation rates and the presence of a tidally locked companion along with the observed internal rotation in the Sun indicate that the mixing is most likely induced by rotation. The current status of non-standard (particularly rotational) stellar models which attempt to account for the lithium observations are reviewed.

We present a review of the present state of knowledge regarding the relative ages of Galactic globular clusters. First, we discuss the relevant galaxy formation models and describe the detailed predictions they make with respect to the formation timescale and chemical evolution of the globular clusters. Next, the techniques used to estimate globular cluster ages are described and evaluated with particular emphasis on the advantages and disadvantages of each method. With these techniques as a foundation, we present arguments in favor of the following assertions: 1) The age of a globular cluster is the likeliest candidate to be the global second parameter, which along with metal abundance, controls the morphology of the horizonal branch. 2) A total age range of as much as $\sim$5 Gyr exists among the bulk of the Galactic globulars. 3) There is a significant relation between age and metallicity among the Galactic globular clusters if the slope of the \mvrr-\feh relation is less than $\sim$0.23. These conclusions along with other supporting evidence favor a formation scenario in which the inner regions of the Galactic halo collapsed in a monotonic fashion over a short time period much less than 1 Gyr. In contrast, the outer regions of the halo fragmented and collapsed in a chaotic manner over several Gyrs.

A grid of stellar evolution models for Procyon A has been calculated. These models include the best physics available to us (including the latest opacities and equation of state) and are based on the revised astrometric mass of Girard et al. (1996). The long standing discrepancy between the evolutionary mass and the astrometric mass is now resolved, a result of the newly determined astrometric mass. Models were calculated with helium diffusion and with the combined effects of helium and heavy element diffusion. Oscillation frequencies for $\ell=0,1,2$ and 3 $p$-modes (and $g$-modes) were calculated for these models. The predicted $p$-mode frequencies are relatively unaffected by heavy element diffusion and convective core overshoot. The inclusion of a modest stellar wind which effectively suppresses the helium diffusion in the surface layers has a modest effect on the $p$-mode frequencies. The evolutionary state (main sequence or shell hydrogen burning) of Procyon A has the largest effect on the predicted $p$-mode frequencies. The $g$-modes show a greater sensitivity to the various model parameters.

We review five independent techniques which are used to set the distance scale to globular clusters, including subdwarf main sequence fitting utilizing the recent Hipparcos parallax catalogue. These data together all indicate that globular clusters are farther away than previously believed, implying a reduction in age estimates. This new distance scale estimate is combined with a detailed numerical Monte Carlo study designed to assess the uncertainty associated with the theoretical age-turnoff luminosity relationship in order to estimate both the absolute age and uncertainty in age of the oldest globular clusters. Our best estimate for the mean age of the oldest globular clusters is now $11.5\pm 1.3 $Gyr, with a one-sided, 95% confidence level lower limit of 9.5 Gyr. This represents a systematic shift of over 2 $\sigma$ compared to our earlier estimate, due completely to the new distance scale---which we emphasize is not just due to the Hipparcos data. This now provides a lower limit on the age of the universe which is consistent with either an open universe, or a flat, matter dominated universe (the latter requiring $H_0 \le 67 \kmsmpc$). Our new study also explicitly quantifies how remaining uncertainties in the distance scale and stellar evolution models translate into uncertainties in the derived globular cluster ages. Simple formulae are provided which can be used to update our age estimate as improved determinations for various quantities become available.

Globular clusters are the oldest objects in the Galaxy whose age may be accurately determined. As such globular cluster ages provide the best estimate for the age of the universe. The age of a globular cluster is determined by a comparison between theoretical stellar evolution models and observational data. Current uncertainties in the stellar models and age dating process are discussed in detail. The best estimate for the absolute age of the globular clusters is $14.6\pm 1.7\,$Gyr. The one-sided, 95\% confidence limit on the lower age of the universe is $12.2\,$Gyr.

Globular cluster age estimates based on the absolute magnitude of the main sequence turn-off (\mvto) are generally considered to be the most reliable from a theoretical viewpoint. However, the difficulty in determining \mvto\ in observed colour-magnitude diagrams leads to a large error in the derived age. In this paper, we advocate the use of the absolute magnitude of the point which is brighter than the turn-off and 0.05 mag redder (\mvb) as a precision age indicator. It is easy to measure this point on observed colour-magnitude diagrams, leading to small observational error bars. Furthermore, an extensive Monte Carlo calculation indicates that the theoretical uncertainty in \mvb\ is similar to \mvto. As a result, ages derived using \mvb\ are at least a factor of 2 more precise than those derived using \mvto. This technique is applied to the globular cluster M68 and an age of $12.8\pm 0.3\,$Gyr is derived (assuming $\mvrr = 0.20\,\feh + 0.98$), indicating that M68 is a `young' globular cluster. A homogeneous set of globular cluster age estimates with this precision would provide unprecedented insight into the formation of the Galactic halo.

Using the redshift evolution of the neutral hydrogen density, as inferred from observations of damped Ly$\alpha$ clouds, we calculate the evolution of star formation rates and elemental abundances in the universe. For most observables our calculations are in rough agreement with previous results based on the instantaneous re-cycling approximation (IRA). However, for the key metallicity tracer Zn, we find a better match to the observed abundance at high redshift than that given by the constant-yield IRA model. We investigate whether the redshift evolution of deuterium, depressions in the diffuse extragalactic gamma-ray background, and measurement of the MeV neutrino background may help determine if observational bias due to dust obscuration is important. We also indicate how the importance of dust on the calculations can be significantly reduced if correlations of the HI column density with metallicity are present. The possibilities for measuring $q_o$ with observations of elemental abundances in damped Ly$\alpha$ systems are discussed.

We report the results of a detailed numerical study designed to estimate both the absolute age and the uncertainty in age (with confidence limits) of the oldest globular clusters. Such an estimate is essential if a comparison with the Hubble age of the universe is to be made to determine the consistency, or lack thereof, of various cosmological models. Utilizing estimates of the uncertainty range (and distribution) in the input parameters of stellar evolution codes we produced 1000 Monte Carlo realizations of stellar isochrones, with which we could fit the ages of the 18 oldest globular clusters. Incorporating the observational uncertainties in the measured color-magnitude diagrams for these systems and the predicted isochrones, we derived a probability distribution for the mean age of these systems. The one-sided 95$\%$ C.L.~ lower bound for this distribution occurs at an age of 12.07 Gyr. This puts interesting constraints on cosmology which we discuss. Further details, including a description of the distributions, covariance matrices, dependence upon individual input parameters, etc.\ will appear in a future article.

Main sequence turnoff magnitudes from the recent set of Yale isochrones (Chaboyer \ea 1995) have been combined with a variety of relations for the absolute magnitude of RR Lyr stars (\mvrr) to calibrate age as a function of the difference in magnitude between the main sequence turn-off and the horizontal branch (\dv). A best estimate for the calibration of \mvrr is derived from a survey of the current literature: $\mvrr = 0.20 \feh + 0.98$. This estimate, together with other calibrations (with slopes ranging from 0.15 to 0.30) has been used to derive \dv ages for 43 Galactic globular clusters. Independent of the choice of \mvrr, there is no strong evidence for an age-Galacto\-centric distance relationship among the 43 globular clusters. However, an age-metallicity relation exists, with the metal-poor clusters being the oldest. A study of the age distribution reveals that an age range of 5 Gyr exists among the bulk of the globular clusters. In addition, about 10\% of the sample are substantially younger, and including them in the analysis increases the age range to 9 Gyr. Once again, these statements are independent of the \mvrr ~relation. Evidence for age being the second parameter governing horizontal branch morphology is found by comparing the average \dv ~age of the second parameter clusters to the normal clusters. The second parameter clusters are found to be on average 2 - 3 Gyr younger than the other clusters, which is consistent with age being the second parameter. These results suggest that globular clusters were formed over an extended period of time, with progressively more metal-rich globular clusters ($\feh \ga -1.7$) being formed at later times.

In this review, I will discuss a few problems which point to the need for improved stellar evolution models of halo stars. Current stellar evolution models do not match the observed $^7$Li abundance patterns, suggesting that the input physics and/or the assumptions used in constructing the models are in need of revision. It appears that all halo stars have suffered some $^7$Li depletion, implying that the primordial $^7$Li abundance is higher than that presently observed in hot halo stars. Observations of abundances of various elements in globular cluster giant branch stars have suggested for some time now that some form of deep mixing, which is not present in theoretical models, occurs in halo stars. The driving mechanism for this mixing, and its incorporation into stellar models remain one of the key problems in stellar modeling. Current theoretical isochrones are able to provide a good match to observed colour-magnitude diagrams. However, there is some evidence that the theoretical luminosity functions are in disagreement with observations. This is an area which requires further study, as it suggests that the relative main sequence/giant branch lifetimes predicted by the models are incorrect. A discussion of some of the uncertainties involved in determining the ages of globular clusters is presented. The absolute ages of globular clusters provide a lower bound to the age of the universe, and so are of great interest to cosmologists. Unfortunately, present uncertainties in stellar models lead to a rather large range in the inferred ages of globular clusters of 11 -- 18 Gyr.

The Yale stellar evolution code has been modified to use the OPAL equation of state tables (Rogers 1994). Stellar models and isochrones were constructed for low metallicity systems ($-2.8 \le [Fe/H] \le -0.6$). Above $M\sim 0.7\,\msun$, the isochrones are very similar to those which are constructed using an equation of state which includes the analytical Debye-Huckel correction at high temperatures. The absolute magnitude of the main sequence turn-off (\mvto) with the OPAL or Debye-Huckel isochrones is about 0.06 magnitudes fainter, at a given age, than \mvto derived from isochrones which do not include the Debye-Huckel correction. As a consequence, globular clusters ages derived using \mvto are reduced by 6 -- 7\% as compared to the ages determined from the standard isochrones. Below $M\sim 0.7\,\msun$, the OPAL isochrones are systematically hotter (by approximately 0.04 in B-V) at a given magnitude as compared to the standard, or Debye-Huckel isochrones. However, the lower mass models fall out of the OPAL table range, and this could be the cause of the differences in the location of the lower main-sequences.

We have studied the importance of the combined effects of rotation, diffusion, and convective overshoot on the p-mode oscillation spectrum and the neutrino flux of the standard solar model. To isolate the various physical affects included in the new rotation plus diffusion models we also constructed solar models to test the significance of diffusion and of overshoot by themselves. In previous studies, models that include helium diffusion during solar evolution were found to improve the predicted p-mode frequencies for some modes and worsen the agreement for others (Guenther \ea 1993). Here we verify this result for both the Bahcall and Loeb (1990) formulation of diffusion and the Proffitt and Michaud (1991) formulation of diffusion. We find that the effects of rotation on the Sun's structure in the outer layers perturbs the $p$-mode frequencies only slightly when compared to the more substantial effects due to diffusion. In the thin overshoot layer (taken here to be $0.1\, H_p$), we have compared the effect of overmixing in a radiative layer versus convective (adiabatic) penetration. Neither radiative overmixing nor adiabatic penetration has any significant effect on the $p$-modes, probably in part because the overshoot layer is constrained to be thin. The predicted neutrino flux in our diffusion plus rotation model is 7.12 SNU for Cl detectors, 127 SNU for Ga detectors and $5.00\times 10^6\,{\rm erg \, cm^{-2} }$ for the $^8$B neutrinos; this is approximately half-way between the standard solar model without diffusion, and the standard solar model with diffusion alone.

We present UBR and H$\alpha$ imaging of NGC 1620, a highly inclined spiral galaxy that contains a large scale, arc-like feature of radius 3 kpc in its outer disk at a distance of $\sim$ 11 kpc from the center. What is unusual about this arc-like feature is its stellar nature and the presence of a luminous star cluster at its center. The arc is fragmented into HII region complexes and OB star clusters and shows two kinks in optical continuum light. It spans an angle of 220$^{\circ}$ on our U image and a full, though fragmented, circle on an unsharp masked R image. It is centered on a young star cluster that is the most luminous clump in blue optical continuum light besides the nucleus of the galaxy. This central star cluster has UBR colors and a surface brightness similar to those of other HII regions, but is a relatively weak H$\alpha$ emitter. It consists of at least three unresolved condensations in optical continuum light. Its location at the center of the arc and its prominence within the galaxy suggests that it has been the site of several generations of supernova explosions that swept up the surrounding gas into a supershell. When it attained a radius of $0.5-1$ kpc, this shell became gravitationally unstable and formed the stars which now delineate the arc. The constraints imposed by the survival of the expanding arc against random stellar motions and the age of the stars in the arc yield a required energy input by a minimum of 400 and a maximum of 6500 supernovae. In this scenario the asymmetry in surface brightness of the arc reflects the radial gradient of the gas density in the disk of NGC 1620, while the kinks reflect inhomogeneities in the original gas distribution with respect to the central star cluster. The supernova superbubble formed at least $5 \times 10^7$ yr ago so that, unless

The main sequence turnoff luminosity is the best stellar `clock' which can be used to determine the absolute ages of globular clusters. This is due to the fact that it is generally assumed that the luminosity and lifetimes of main sequence globular cluster stars are independent of the properties of stellar convection and atmospheres, two areas of stellar evolution which are poorly understood. Several possible sources of error in this stellar clock are discussed, and isochrones are constructed using a variety of different physical assumptions. The mean age of the oldest globular clusters are determined from these isochrones and it is found that the uncertainties in the input physics can lead to changes in the derived age of $\pm 15\%$. Surprisingly the largest source of error is the mixing length theory of convection. It is well known that uncertainties in the distance scale and chemical composition of globular cluster stars lead to changes of order $\sim 22\%$ in the determination of absolute ages. Combining the various sources of error, the absolute age of the oldest globular clusters are found to lie in the range 11 --- 21 Gyr. This is meant to be a total theoretical range. For the standard inflationary model ($\Omega = 1, \Lambda = 0$), a minimum age of the universe of 11 Gyr requires $H_o \la 60~\kmsmpc$.

A new grid of theoretical isochrones based on the Yale stellar evolution code using the OPAL and Kurucz opacities has been constructed. The grid of isochrones spans a wide range of metallicities, helium abundances and masses. The construction of the isochrones is described and the isochrones are compared to galactic globular cluster observations. A solar calibrated mixing length ($\alpha = 1.7$) yields a good fit to globular cluster colour-magnitude diagrams. Ages for 40 globular clusters are determined using the \dv method and the formation of the halo is discussed.

Stellar models with masses ranging from 0.5 to $1.3~M_\odot$ were constructed in order to compare to young cluster observations of Li and of rotation velocities. The amount of Li depletion in cool stars is sensitive to the amount of overshoot at the base of the surface convection zone, and the exact metallicity of the models. Even when this uncertainty is taken into account, the Li observations are a severe constraint for the models and rule out standard models, and pure diffusion models. Stellar models which include diffusion and rotational mixing in the radiative regions of stars are able to simultaneously match the Li abundances observed in the Pleiades, UMaG, Hyades, Praesepe, NGC 752 and M67. They also match the observed rotation periods in the Hyades. However, these models are unable to simultaneously explain the presence of the rapidly rotating late G and K stars in the Pleiades and the absence of rapidly rotating late F and early G stars.

The Yale stellar evolution code has been modified to include the combined effects of diffusion and rotational mixing on $^1$H, $^4$He and the trace elements $^3$He, $^6$Li, $^7$Li, and $^9$Be. The interaction between rotational mixing and diffusion is studied by calculating a number of calibrated solar models. The rotational mixing inhibits the diffusion in the outer parts of the models, leading to a decrease in the envelope diffusion by 25 -- 50\%. Conversely, diffusion leads to gradients in mean molecular weight which can inhibit the rotational mixing. The degree to which gradients in mean molecular weight inhibits the rotational mixing is somewhat uncertain. A comparison to the observed solar oblateness suggests that gradients in the mean molecular weight play a smaller role in inhibiting the rotational mixing previously believed. This is reinforced by the fact that the model with the standard value for the inhibiting effect of mean molecular weight on the rotational mixing depletes no Li on the main sequence. This is in clear contrast to the observations. A reduction in the inhibiting effect of mean molecular weight gradients by a factor of ten leads to noticeable main sequence Li depletion.

Light element ($^6$Li, $^7$Li and $^9$Be) depletion isochrones for halo stars have been calculated with standard stellar evolution models. These models include the latest available opacities and are computed through the sub-giant branch. If $^6$Li is not produced in appreciable amounts by stellar flares, then the detection of $^6$Li in HD 84937 by Smith, Lambert \& Nissen (1993) is compatible with standard stellar evolution and standard big bang nucleosynthesis only if HD 84937 is a sub-giant. The present parallax is inconsistent with HD 84937 being a sub-giant star at the $2.5\, \sigma$ level. The most metal poor star with a measured $^9$Be abundance is HD 140283, which is a relatively cool sub-giant. Standard stellar evolution predict that $^9$Be will have been depleted in this star by $\sim 0.3$ dex (for ${\rm T_{eff}} = 5640$ K). Revising the abundance upward changes the oxygen-beryllium relation, suggesting incompatible with standard comic ray production models, and hence, standard big bang nucleosynthesis. However, an increase in the derived temperature of HD 140283 to 5740 K would result in no depletion of $^9$Be and agreement with standard big bang nucleosynthesis.

A large number of stellar evolution models with [Fe/H] = -2.3 and -3.3 have been calculated in order to determine the primordial $^7$Li abundance and to test current stellar evolution models by a comparison to the extensive database of Li abundances in extremely metal poor halo stars observed by Thorburn (1994). Standard models do a good job of fitting the observed Li abundances in stars hotter than 5600 K. They predict a primordial $^7$Li abundance of Log N(Li) = 2.24\pm 0.03$. Models which include microscopic diffusion predict a downward curvature in the $^7$Li destruction isochrones at hot temperatures which is not present in the observations. Thus, the observations clearly rule out models which include uninhibited microscopic diffusion of $^7$Li from the surface of the star. The [Fe/H] = -2.3 stellar models which include both diffusion and rotational mixing provide an excellent match to the observations. Both the plateau stars and the heavily depleted cool stars are well fit by these models. The rotational mixing leads to considerable $^7$Li depletion in these models and the primordial $^7$Li abundance inferred from these models is Log N(Li) = $3.08\pm 0.1$.