We present MINESweeper, a tool to measure stellar parameters by jointly fitting observed spectra and broadband photometry to model isochrones and spectral libraries. This approach enables the measurement of spectrophotometric distances, in addition to stellar parameters such as Teff, log(g), [Fe/H], [a/Fe], and radial velocity. MINESweeper employs a Bayesian framework and can easily incorporate a variety of priors, including Gaia parallaxes. Mock data are fit in order to demonstrate how the precision of derived parameters depends on evolutionary phase and SNR. We then fit a selection of data in order to validate the model outputs. Fits to a variety of benchmark stars including Procyon, Arcturus, and the Sun result in derived stellar parameters that are in good agreement with the literature. We then fit combined spectra and photometry of stars in the open and globular clusters M92, M13, M3, M107, M71, and M67. Derived distances, [Fe/H], [a/Fe], and log(g)-Teff, relations are in overall good agreement with literature values, although there are trends between metallicity and log(g), within clusters that point to systematic uncertainties at the ~0.1 dex level. Finally, we fit a large sample of stars from the H3 Spectroscopic Survey in which high quality Gaia parallaxes are also available. These stars are fit without the Gaia parallaxes so that the geometric parallaxes can serve as an independent test of the spectrophotometric distances. Comparison between the two reveals good agreement within their formal uncertainties after accounting for the Gaia zero point uncertainties.

It has been known for several decades that transport of chemical elements is induced by the process of microscopic atomic diffusion. Yet, the effect of atomic diffusion, including radiative levitation, has hardly been studied in the context of gravity mode pulsations of core-hydrogen burning stars. In this paper, we study the difference in the properties of such modes for models with and without atomic diffusion. We perform asteroseismic modeling of two slowly rotating A- and F-type pulsators, KIC11145123 ($f_{\rm rot} \approx0.010~{\rm d}^{-1}$) and KIC9751996 ($f_{\rm rot} \approx0.0696~{\rm d}^{-1}$), respectively, based on the periods of individual gravity modes. For both stars, we find models whose g-mode periods are in very good agreement with the {\it Kepler\/} asteroseismic data, keeping in mind that the theoretical/numerical precision of present-day stellar evolution models is typically about two orders of magnitude lower than the measurement errors. Using the Akaike Information Criterion (AIC) we have made a comparison between our best models with and without diffusion, and found very strong evidence for signatures of atomic diffusion in the pulsations of KIC11145123. In the case of KIC9751996 the models with atomic diffusion are not able to explain the data as well as the models without it. Furthermore, we compare the observed surface abundances with those predicted by the best fitting models. The observed abundances are inconclusive for KIC9751996, while those of KIC11145123 from the literature can better be explained by a model with atomic diffusion.

We use photometric and spectroscopic observations of the eclipsing binary E32 in the globular cluster 47 Tuc to derive the masses, radii, and luminosities of the component stars. The system has an orbital period of 40.9 d, a markedly eccentric orbit with e = 0.24, and is shown to be a member of or a recent escaper from the cluster. We obtain Mp = 0.862(5) Msun , Rp = 1.183(3) Rsun , Lp = 1.65(5) Lsun for the primary and Ms = 0.827(5) Msun , Rs = 1.004(4) Rsun , Ls = 1.14(4) Lsun for the secondary. Based on these data and on an earlier analysis of the binary V69 in 47 Tuc we measure the distance to the cluster from the distance moduli of the component stars, and, independently, from a color - surface brightness calibration. We obtain 4.55(3) and 4.50(7) kpc, respectively - values compatible within 1 sigma with recent estimates based on Gaia DR2 parallaxes. By comparing the M - R diagram of the two binaries and the color-magnitude diagram of 47 Tuc to Dartmouth model isochrones we estimate the age of the cluster to be 12.0 pm 0.5 Gyr, and the helium abundance of the cluster to be Y approx 0.25.

The extended main sequence turn offs (eMSTOs) of several young to intermediate age clusters are examined in the Magellanic Clouds and the Milky Way. We explore the effects of extended star formation (eSF) and a range of stellar rotation rates on the behavior of the color-magnitude diagram (CMD), paying particular attention to the MSTO. We create synthetic stellar populations based on MESA stellar models to simulate observed Hubble Space Telescope and Gaia star cluster data. We model the effect of rotation as a non-parametric distribution, allowing for maximum flexibility. In our models the slow rotators comprise the blueward, and fast rotators the redward portion of the eMSTO. We simulate data under three scenarios: non-rotating eSF, a range of rotation rates with a single age, and a combination of age and rotation effects. We find that two of the five clusters (the youngest and oldest) favor an age spread, but these also achieve the overall worst fits of all clusters. The other three clusters show comparable statistical evidence between rotation and an age spread. In all five cases, a rotation rate distribution alone is capable of qualitatively matching the observed eMSTO structure. In future work, we aim to compare our predicted Vsin(i) with observations in order to better constrain the physics related to stellar rotation.

The Asymptotic Giant Branch (AGB) scenario ascribes the multiple populations in old Galactic Globular Clusters (GGC) to episodes of star formation in the gas contaminated by the ejecta of massive AGBs and super-AGBs of a first stellar population. The mass of these AGBs (4-8Msun) today populate the Young Massive Clusters (YMC) of the Magellanic Clouds, where rapid rotation and its slowing down play an important role in shaping the color magnitude diagram features. Consequently, we must reconsider whether the rotational evolution of these masses affects the yields, and whether the resulting abundances are compatible with the chemical patterns observed in GGC. We show the first results of a differential analysis, by computing the hot bottom burning evolution of non rotating models with increased CNO-Na abundances at the 2DU, following the results of MESA rotational models.

The recent measurements of internal variations of helium in Galactic and extragalactic Globular Clusters (GCs) set binding constraints to the models of formation of Multiple Populations (MPs) in GCs, and gave rise, at the same time, to crucial questions related with the influence of the environment on MP formation as well as with the role played by GCs in the early galactic formation. We present the most recent estimates of helium enrichment in the main populations of a large sample of Galactic and extragalactic GCs.

Context. The benchmark open cluster M67 is known to have solar metallicity and similar age as the Sun. It thus provides us a great opportunity to study the properties of solar twins, as well as the evolution of Sun-like stars. Aims. Previous spectroscopic studies reported to detect possible subtle changes in stellar surface abundances throughout the stellar evolutionary phase, namely the effect of atomic diffusion, in M67. In this study we attempt to confirm and quantify more precisely the effect of atomic diffusion, as well as to explore the level of chemical (in)homogeneity in M67. Methods. We presented a strictly line-by-line differential chemical abundance analysis of two groups of stars in M67: three turn-off stars and three sub-giants. Stellar atmospheric parameters and elemental abundances were obtained with very high precision using the Keck/HIRES spectra. Results. The sub-giants in our sample show negligible abundance variations ($\le$ 0.02 dex), which implies that M67 was born chemically homogeneous. We note there is a significant abundance difference ($\sim$ 0.1 - 0.2 dex) between sub-giants and turn-off stars, which can be interpreted as the signature of atomic diffusion. Qualitatively stellar models with diffusion agree with the observed abundance results. Some turn-off stars do not follow the general pattern, which suggests that in some cases diffusion can be inhibited, or they might suffered some sort of mixing event related to planets. Conclusions. Our results pose additional challenges for chemical tagging when using turn-off stars. In particular, the effects of atomic diffusion, which could be as large as 0.1 - 0.2 dex, must be taken into account in order for chemical tagging to be successfully applied.

We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (MESA). RSP is a new functionality in MESAstar that models the non-linear radial stellar pulsations that characterize RR Lyrae, Cepheids, and other classes of variable stars. We significantly enhance numerical energy conservation capabilities, including during mass changes. For example, this enables calculations through the He flash that conserve energy to better than 0.001 %. To improve the modeling of rotating stars in MESA, we introduce a new approach to modifying the pressure and temperature equations of stellar structure, and a formulation of the projection effects of gravity darkening. A new scheme for tracking convective boundaries yields reliable values of the convective-core mass, and allows the natural emergence of adiabatic semiconvection regions during both core hydrogen- and helium-burning phases. We quantify the parallel performance of MESA on current generation multicore architectures and demonstrate improvements in the computational efficiency of radiative levitation. We report updates to the equation of state and nuclear reaction physics modules. We briefly discuss the current treatment of fallback in core-collapse supernova models and the thermodynamic evolution of supernova explosions. We close by discussing the new MESA Testhub software infrastructure to enhance source-code development.

Early type galaxies (ETGs) frequently show emission from warm ionized gas. These Low Ionization Emission Regions (LIERs) were originally attributed to a central, low-luminosity active galactic nuclei. However, the recent discovery of spatially-extended LIER emission suggests ionization by both a central source and an extended component that follows a stellar-like radial distribution. For passively-evolving galaxies with old stellar populations, hot post-Asymptotic Giant Branch (AGB) stars are the only viable extended source of ionizing photons. In this work, we present the first prediction of LIER-like emission from post-AGB stars that is based on fully self-consistent stellar evolution and photoionization models. We show that models where post-AGB stars are the dominant source of ionizing photons reproduce the nebular emission signatures observed in ETGs, including LIER-like emission line ratios in standard optical diagnostic diagrams and H$\alpha$ equivalent widths of order 0.1-3 angstroms. We test the sensitivity of LIER-like emission to the details of post-AGB models, including the mass loss efficiency and convective mixing efficiency, and show that line strengths are relatively insensitive to post-AGB timescale variations. Finally, we examine the UV-optical colors of the models and the stellar populations responsible for the UV-excess observed in some ETGs. We find that allowing as little as 3% of the HB population to be uniformly distributed to very hot temperatures (30,000 K) produces realistic UV colors for old, quiescent ETGs.

We present an analysis of the second epoch HST WFC3 F110W near-Infrared (NIR) imaging data of the globular cluster M4. The new dataset suggests that one of the previously suggested four brown dwarf candidates in this cluster is indeed a high-probability cluster member. The position of this object in the NIR colour magnitude diagrams (CMDs) is in the white dwarf/brown dwarf area. The source is too faint to be a low-mass main sequence star, but, according to theoretical considerations, also most likely somewhat too bright to be a bona-fide brown dwarf. Since we know that the source is a cluster member, we determined a new optical magnitude estimate at the position the source should have in the optical image. This new estimate places the source closer to the white dwarf sequence in the optical-NIR CMD and suggests that it might be a very cool (T_eff < 4500 K) white dwarf at the bottom of the white dwarf cooling sequence in M4, or a white dwarf/brown dwarf binary. We cannot entirely exclude the possibility that the source is a very massive, bright brown dwarf, or a very low-mass main sequence star, however, we conclude that we still have not convincingly detected a brown dwarf in a globular cluster, but we expect to be very close to the start of the brown dwarf cooling sequence in this cluster. We also note that the main sequence ends at F110W approx.22.5 mag in the proper-motion cleaned CMDs, where completeness is still high.

The multiple stellar populations (MPs) of $\sim$11-13 Gyr-old Globular Clusters (GCs) in our Galaxy are characterized by different content of several light elements. These elements describe well-defined patterns like the C-N and the Na-O anticorrelations and the He-N and Na-N correlations. The discovery of the MPs in Magellanic Cloud GCs opened up new paths for the investigation of chemical anomalies in clusters with different age and physical properties. In this context, we used Hubble Space Telescope photometry to investigate the MPs and constrain their chemical composition of four $\sim$6-11 Gyr extragalactic GCs, namely NGC 121,NGC 339, NGC 416 and Lindsay 1 in the Small Magellanic Cloud. The comparison of the stellar colors with synthetic spectra suggests that second-population stars of NGC 121, NGC 339, NGC 416 are slightly enhanced in helium by $\delta$Y = 0.009$\pm$0.006, 0.007$\pm$0.004 and 0.010$\pm$0.003, respectively, with respect to the first population, while we find no significant helium variation in Lindsay 1 ($\delta$Y = 0.000$\pm$0.004). Moreover, second-population stars of all the clusters are, on average, enhanced in nitrogen and depleted in carbon and oxygen, in close analogy with what we observe in Galactic GCs.

The overlap between the spectroscopic Galactic Archaeology with HERMES (GALAH) survey & $Gaia$ provides a high-dimensional chemodynamical space of unprecedented size. We present a first analysis of a subset of this overlap, of 7066 dwarf, turn-off, & sub-giant stars. [...] We investigate correlations between chemical compositions, ages, & kinematics for this sample. Stellar parameters & elemental abundances are derived from the GALAH spectra with the spectral synthesis code SME. [...] We report Li, C, O, Na, Mg, Al, Si, K, Ca, Sc, Ti, V, Cr, Mn, Co, Ni, Cu, Zn, Y, as well as Ba & we note that we employ non-LTE calculations for Li, O, Al, & Fe. We show that the use of astrometric & photometric data improves the accuracy of the derived spectroscopic parameters, especially $\log g$. [...] we recover the result that stars of the high-$\alpha$ sequence are typically older than stars in the low-$\alpha$ sequence, the latter spanning $-0.7<$[Fe/H]$<+0.5$. While these two sequences become indistinguishable in [$\alpha$/Fe] vs. [Fe/H] at the metal-rich regime, we find that age can be used to separate stars from the extended high-$\alpha$ & the low-$\alpha$ sequence even in this regime. [...] we find that the old stars ($>8$ Gyr have lower angular momenta $L_z$ than the Sun, which implies that they are on eccentric orbits & originate from the inner disk. Contrary to some previous smaller scale studies we find a continuous evolution in the high-$\alpha$-sequence up to super-solar [Fe/H] rather than a gap, which has been interpreted as a separate "high-$\alpha$ metal-rich" population. Stars in our sample that are younger than 10 Gyr, are mainly found on the low $\alpha$-sequence & show a gradient in $L_z$ from low [Fe/H] ($L_z>L_{z,\odot}$) towards higher [Fe/H] ($L_z<L_{z,\odot}$), which implies that the stars at the ends of this sequence are likely not originating from the close solar vicinity.

We present photometric analysis of twelve Galactic open clusters and show that the same multiple-population phenomenon observed in Magellanic Clouds (MCs) is present in nearby open clusters. Nearly all the clusters younger than $\sim$2.5 Gyr of both MCs exhibit extended main-sequence turnoffs (eMSTOs) and all the cluster younger than $\sim$700 Myr show broadened/split main sequences (MSs). High-resolution spectroscopy has revealed that these clusters host stars with a large spread in the observed projected rotations. In addition to rotation, internal age variation is indicated as a possible responsible for the eMSTOs, making these systems the possible young counterparts of globular clusters with multiple populations. Recent work has shown that the eMSTO+broadened MSs are not a peculiarity of MCs clusters. Similar photometric features have been discovered in a few Galactic open clusters, challenging the idea that the color-magnitude diagrams (CMDs) of these systems are similar to single isochrones and opening new windows to explore the eMSTO phenomenon. We exploit photometry+proper motions from Gaia DR2 to investigate the CMDs of open clusters younger than $\sim$1.5 Gyr. Our analysis suggests that: (i) twelve open clusters show eMSTOs and/or broadened MSs, that cannot be due neither to field contamination, nor binaries; (ii) split/broadened MSs are observed in clusters younger than $\sim$700 Myr, while older objects display only an eMSTO, similarly to MCs clusters; (iii) the eMSTO, if interpreted as a pure age spread, increases with age, following the relation observed in MCs clusters and demonstrating that rotation is the responsible for this phenomenon.

The Pan-Andromeda Archaeological Survey is a survey of $>400$ square degrees centered on the Andromeda (M31) and Triangulum (M33) galaxies that has provided the most extensive panorama of a $L_\star$ galaxy group to large projected galactocentric radii. Here, we collate and summarise the current status of our knowledge of the substructures in the stellar halo of M31, and discuss connections between these features. We estimate that the 13 most distinctive substructures were produced by at least 5 different accretion events, all in the last 3 or 4 Gyrs. We suggest that a few of the substructures furthest from M31 may be shells from a single accretion event. We calculate the luminosities of some prominent substructures for which previous estimates were not available, and we estimate the stellar mass budget of the outer halo of M31. We revisit the problem of quantifying the properties of a highly structured dataset; specifically, we use the OPTICS clustering algorithm to quantify the hierarchical structure of M31's stellar halo, and identify three new faint structures. M31's halo, in projection, appears to be dominated by two `mega-structures', that can be considered as the two most significant branches of a merger tree produced by breaking M31's stellar halo into smaller and smaller structures based on the stellar spatial clustering. We conclude that OPTICS is a powerful algorithm that could be used in any astronomical application involving the hierarchical clustering of points. The publication of this article coincides with the public release of all PAndAS data products.

Galactic globular clusters (GCs) show overwhelming photometric and spectroscopic evidence for the existence of multiple stellar populations. The question of whether or not there exists a GC that represents a true 'simple stellar population' remains open. Here we focus on Ruprecht 106 (R106), a halo GC with [Fe/H]=-1.5 and [alpha/Fe]~0. A previous spectroscopic study found no sign of the Na-O anticorrelation among 9 of its brightest red giants, which led to the conclusion that R106 is a true simple stellar population GC. Here we present new Hubble Space Telescope (HST) Wide Field Camera 3 photometry of R106 that, when combined with archival HST images spanning a 6-year baseline, allows us to create proper motion cleaned color-magnitude diagrams spanning the ultraviolet (F336W) to the near-infrared (F814W). These data allow us to construct the pseudo-color C_{U,B,I} that is sensitive to the presence of light-element abundance spreads. We find no evidence of a split along the red giant branch (RGB) in the C_{U,B,I} diagram but the width of the RGB (sigma_CUBI = 0.015) is marginally broader than expected from artificial star tests (sigma_CUBI = 0.009). The observed spread in C_{U,B,I} is smaller than any other Galactic GC studied to date. Our results raise important questions about the role of formation environment and primordial chemical composition in the formation of multiple stellar populations in GCs.

We use Hubble Space Telescope multicolour photometry of the globular cluster 47 Tucanae to uncover a population of 24 objects with no previous classification that are outliers from the single-star model tracks in the colour-magnitude diagram and yet are likely cluster members. By comparing those sources with evolutionary models and X-ray source catalogues, we were able to show that the majority of those sources are likely binary systems that do not have any X-ray source detected nearby, most possibly formed by a white dwarf and a main-sequence star and a small number of possible double-degenerate systems.

Open cluster members are coeval and share the same initial bulk chemical compositions. Consequently, differences in surface abundances between members of a cluster that are at different evolutionary stages can be used to study the effects of mixing and internal chemical processing. We carry out an abundance analysis of seven elements (Li, O, Na, Mg, Al, Si, Fe) in 66 stars belonging to the open cluster M67, based on high resolution GALAH spectra, 1D MARCS model atmospheres, and, for the first time for a large sample of M67 stars, non-local thermodynamic equilibrium (non-LTE) radiative transfer. From the non-LTE analysis, we find a typical star-to-star scatter in the abundance ratios of around 0.05 dex; this scatter is slightly but systematically larger when LTE is assumed instead. We find trends in the abundance ratios with effective temperature, indicating systematic differences in the surface abundances between turn-off and giant stars; these trends are more pronounced when LTE is assumed. However, in the non-LTE analysis, most of the element trends have been flattened. Two species are exceptions to this behaviour, namely Al and Si, which both clearly display remaining trends in the non-LTE analysis. We comment on the possible origin of these trends, by comparing them with recent stellar models that include atomic diffusion.

We present the first results from the GeMS/GSAOI Galactic Globular Cluster Survey (G4CS) of the Milky-Way globular clusters (GCs) NGC 3201 and NGC 2298. Using the Gemini South Adaptive Optics Imager (GSAOI), in tandem with the Gemini Multi-conjugate adaptive optics System (GeMS) on the 8.1-meter Gemini-South telescope, we collected deep near-IR observations of both clusters, resolving their constituent stellar populations down to $K_s\simeq21$ Vega mag. Point spread function (PSF) photometry was performed on the data using spatially-variable PSFs to generate $JHK_{s}$ photometric catalogues for both clusters. These catalogues were combined with Hubble Space Telescope (HST) data to augment the photometric wavelength coverage, yielding catalogues that span the near-ultraviolet (UV) to near-infrared (near-IR). We then applied 0.14 mas/year accurate proper-motion cleaning, differential-reddening corrections and chose to anchor our isochrones using the lower main-sequence knee (MSK) and the main-sequence turn-off (MSTO) prior to age determination. As a result of the data quality, we found that the $K_{s}$ vs. F606W$-K_{s}$ and F336W vs. F336W$-K_{s}$ color-magnitude diagrams (CMDs) were the most diagnostically powerful. We used these two color combinations to derive the stellar-population ages, distances and reddening values for both clusters. Following isochrone-fitting using three different isochrone sets, we derived best-fit absolute ages of $12.2\pm0.5$ Gyr and $13.2\pm0.4$ Gyr for NGC 3201 and NGC 2298, respectively. This was done using a weighted average over the two aforementioned color combinations, following a pseudo-$\chi^2$ determination of the best-fit isochrone set. Our derived parameters are in good agreement with recent age determinations of the two clusters, with our constraints on the ages being or ranking among the most statistically robust.

Stellar photometric variability offers a novel probe of the interior structure and evolutionary state of stars. Here we present a census of stellar variability on day to decade timescales across the color-magnitude diagram for 73,000 stars brighter than $M_I$=-5 in the Whirlpool Galaxy (M51). Our Cycle 24 HST program acquired V and I-band images over 34 epochs spanning one year with pseudo-random cadences enabling sensitivity to periods from days to months. We supplement these data with archival V and I-band HST data obtained in 1995 and 2005 providing sensitivity to variability on decade timescales. At least 50% of stars brighter than $M_I$=-7 show strong evidence for variability within our Cycle 24 data; among stars with V-I>2 the variability fraction rises to ~100%. Large amplitude variability (>0.3 mag) on decade timescales is restricted to red supergiants and very luminous blue stars. Both populations display fairly smooth variability on month-year timescales. The Cepheid instability strip is clearly visible in our data, although the variability fraction within this region never exceeds ~10%. The location of variable stars across the color magnitude diagram broadly agrees with theoretical sources of variability, including the instability strip, red supergiant pulsational instabilities, long-period fundamental mode pulsations, and radiation-dominated envelopes in massive stars. Our data can be used to place stringent constraints on the precise onset of these various instabilities and their lifetimes and growth rates.

The Hyades, Praesepe, and Pleiades are well studied stellar clusters that anchor important secondary stellar age indicators. Recent studies have shown that main sequence turn off-based ages for these clusters may depend on the degree of rotation in the underlying stellar models. Rotation induces structural instabilities that can enhance the chemical mixing of a star, extending its fuel supply. In addition, rotation introduces a modulation of the star's observed magnitude and color due to the effects of gravity darkening. We aim to investigate the extent to which stellar rotation affects the age determination of star clusters. We utilize the MESA stellar evolution code to create models that cover a range of rotation rates corresponding to $\Omega/\Omega_c=0.0$ to $0.6$ in $0.1$ dex steps, allowing the assessment of variations in this dimension. The statistical analysis package, MATCH, is employed to derive ages and metallicities by fitting our MESA models to Tycho $B_T$, $V_T$ and 2MASS $J$, $K_s$ color-magnitude diagrams. We find that the derived ages are relatively insensitive to the effects of rotation. For the Hyades, Praesepe, and Pleiades, we derive ages based on synthetic populations that model a distribution of rotation rates or a fixed rate. Across each case, derived ages tend to agree roughly within errors, near $680$, $590$, and $110-160$ Myr for the Hyades, Praesepe, and Pleiades, respectively. These ages are in agreement with Li depletion boundary-based ages and previous analyses that used non-rotating isochrones. Our methods do not provide a strong constraint on the metallicities of these clusters.

We present a spectroscopic analysis of main sequence (MS) stars in the young globular cluster NGC1818 (age~40 Myrs) in the Large Magellanic Cloud. Our photometric survey on Magellanic Clouds clusters has revealed that NGC1818, similarly to the other young objects with age 600 Myrs, displays not only an extended MS Turn-Off (eMSTO), as observed in intermediate-age clusters (age~1-2 Gyrs), but also a split MS. The most straightforward interpretation of the double MS is the presence of two stellar populations: a sequence of slowly-rotating stars lying on the blue-MS and a sequence of fast rotators, with rotation close to the breaking speed, defining a red-MS. We report the first direct spectroscopic measurements of projected rotational velocities vsini for the double MS, eMSTO and Be stars of a young cluster. The analysis of line profiles includes non-LTE effects, required for correctly deriving v sini values. Our results suggest that: (i) the mean rotation for blue- and red-MS stars is vsini=71\pm10 km/s (sigma=37 km/s) and vsini=202\pm23 km/s (sigma=91 km/s), respectively; (ii) eMSTO stars have different vsini, which are generally lower than those inferred for red-MS stars, and (iii) as expected, Be stars display the highest vsini values. This analyis supports the idea that distinct rotational velocities play an important role in the appearence of multiple stellar populations in the color-magnitude diagrams of young clusters, and poses new constraints to the current scenarios.

We use the framework developed as part of the MESA Isochrones and Stellar Tracks (MIST) project to assess the utility of several types of observables in jointly measuring the age and 1D stellar model parameters in star clusters. We begin with a pedagogical overview summarizing the effects of stellar model parameters, such as the helium abundance, mass-loss efficiency, and the mixing length parameter, on observational diagnostics including the color-magnitude diagram, mass-radius relation, and surface abundances, amongst others. We find that these parameters and the stellar age influence observables in qualitatively distinctive, degeneracy-breaking ways. To assess the current state of affairs, we use the recent Gaia Data Release 2 (DR2) along with data from the literature to investigate three well-studied old open clusters---NGC6819, M67, NGC6791---as case studies. Although there is no obvious tension between the existing observations and the MIST models for NGC6819, there are interesting discrepancies in the cases of M67 and NGC6791. At this time, parallax zero point uncertainties in Gaia DR2 remain one of the limiting factors in the analysis of these clusters. With a combination of exquisite photometry, parallax distances, and cluster memberships from Gaia at the end of its mission, we anticipate precise and accurate ages for these and other star clusters in the Galaxy.

We present a color-based method for identifying red giants based on Pan-STARRS grz and WISE W1 and W2 photometry. We utilize a subsample of bright stars with precise parallaxes from Gaia DR2 to verify that the color-based selection reliably separates dwarfs from giants. The selection is conservative in the sense that contamination is small (~30%) but not all giants are included (the selection primarily identifies K giants). The color-based selection can be applied to stars brighter than $W1\approx16$, more than two magnitudes fainter than techniques relying on shallower 2MASS photometry. Many streams and clouds are visible in the resulting sky maps, especially when binned by Gaia DR2 proper motions, including the Sagittarius stream, the Hercules-Aquila Cloud, the Eastern Banded Structure, Monoceros, and the Virgo Overdensity. In addition to the characterization of new and known stellar streams, we expect that this method for selecting red giants will enable detailed analysis of the diffuse stellar halo to distances exceeding 100 kpc.

Observations of young open clusters show a bimodal distribution of rotation periods that has been difficult to explain with existing stellar spin-down models. Detailed MHD stellar wind simulations have demonstrated that surface magnetic field morphology has a strong influence on wind-driven angular momentum loss. Observations suggest that faster rotating stars store a larger fraction of their magnetic flux in higher-order multipolar components of the magnetic field. In this work, we present a new model for stellar spin-down that, for the first time, accounts for the stellar surface magnetic field configuration. We show how a magnetic complexity that evolves from complex toward simple configurations as a star spins down can explain the salient features of stellar rotation evolution, including the bimodal distribution of both slow and fast rotators seen in young open clusters.

The effective temperature (Teff) distribution of stellar evolution models along the red giant branch (RGB) is sensitive to a number of parameters including the overall metallicity, elemental abundance patterns, the efficiency of convection, and the treatment of the surface boundary condition. Recently there has been interest in using observational estimates of the RGB Teff to place constraints on the mixing length parameter, a_MLT, and possible variation with metallicity. Here we use 1D MESA stellar evolution models to explore the sensitivity of the RGB Teff to the treatment of the surface boundary condition. We find that different surface boundary conditions can lead to +/- 100 K metallicity-dependent offsets on the RGB relative to one another in spite of the fact that all models can reproduce the properties of the Sun. Moreover, for a given atmosphere T-tau relation, we find that the RGB Teff is also sensitive to the optical depth at which the surface boundary condition is applied in the stellar model. Nearly all models adopt the photosphere as the location of the surface boundary condition but this choice is somewhat arbitrary. We compare our models to stellar parameters derived from the APOGEE-Kepler sample of first ascent red giants and find that systematic uncertainties in the models due to treatment of the surface boundary condition place a limit of ~100 K below which it is not possible to make firm conclusions regarding the fidelity of the current generation of stellar models.

We present a Bayesian implementation of isochrone fitting in deriving stellar ages and masses, incorporating absolute K magnitude ($\rm M_K$) derived from 2MASS photometry and Gaia DR1 parallax and differentiation between initial bulk metallicity and present day surface metallicity, with allowance for incorporating further constraints (e.g., asteroseismology) when available. As a test, we re-computed stellar ages and masses of $\sim 4000$ stars in the solar neighbourhood from six well-studied literature samples using both Hipparcos and TGAS parallaxes. Our ages are found to be compatible with literature values but with reduced uncertainties in general. The inclusion of parallax-based $\rm M_K$ serves as an additional constraint on the derived quantities, especially when systematic errors in stellar parameters are underestimated. We reconstructed the age-metallicity relationship in the solar neighbourhood by re-analysing the Geneva-Copenhagen Survey with the inclusion of TGAS-parallaxes and initial bulk metallicity sampling. We found a flat trend for disk stars with ages $<$11\,Gyr but with smaller scatter at all ages compared to literature.

The split main sequences (MSs) and extended MS turnoffs (eMSTOs) detected in a few young clusters have demonstrated that these stellar systems host multiple populations differing in a number of properties such as rotation and, possibly, age.We analyze Hubble Space Telescope photometry for thirteen clusters with ages between ~40 and ~1000 Myrs and of different masses. Our goal is to investigate for the first time the occurrence of multiple populations in a large sample of young clusters. We find that all the clusters exhibit the eMSTO phenomenon and that MS stars more massive than ~1.6 solar masses define a blue and red MS, with the latter hosting the majority of MS stars. The comparison between the observations and isochrones suggests that the blue MSs are made of slow-rotating stars, while the red MSs host stars with rotational velocities close to the breakup value. About half of the bright MS stars in the youngest clusters are H-alpha emitters. These Be stars populate the red MS and the reddest part of the eMSTO thus supporting the idea that the red MS is made of fast rotators. We conclude that the split MS and the eMSTO are a common feature of young clusters in both Magellanic Clouds. The phenomena of a split MS and an eMSTO occur for stars that are more massive than a specific threshold which is independent of the host-cluster mass. As a by-product, we report the serendipitous discovery of a young SMC cluster, GALFOR1.

High-resolution spectroscopic observations were taken of 29 extended main sequence turn-off (eMSTO) stars in the young ($\sim$200 Myr) LMC cluster, NGC 1866 using the Michigan/Magellan Fiber System and MSpec spectrograph on the Magellan-Clay 6.5-m telescope. These spectra reveal the first direct detection of rapidly rotating stars whose presence has only been inferred from photometric studies. The eMSTO stars exhibit H-alpha emission (indicative of Be-star decretion disks), others have shallow broad H-alpha absorption (consistent with rotation $\gtrsim $150 km s$^{-1}$), or deep H-alpha core absorption signaling lower rotation velocities ($ \lesssim $150 km s$^{-1}$ ). The spectra appear consistent with two populations of stars - one rapidly rotating, and the other, younger and slowly rotating.

Chemical abundances are presented for 19 elements in a sample of 63 red giants in the Carina dwarf spheroidal galaxy (dSph), based on homogeneous 1D/LTE model atmosphere analyses of our own observations (32 stars) and data available in the literature (a further 31 independent stars). The (Fe) metallicity and [$\alpha$/Fe] distribution functions have mean values and dispersions of -1.59 and 0.33 dex ([Fe/H] range: -2.68 to -0.64), and 0.07 and 0.13 dex ([$\alpha$/Fe] range: -0.27 to 0.25), respectively. We confirm the finding of Venn et al. (2012) that a small percentage (some 10% in the present investigation) of the sample show clear evidence for significant enrichment by Type Ia supernovae ejecta. Calcium, with the most accurately determined abundance of the alpha-elements, shows an asymmetric distribution towards smaller values of [Ca/Fe] at all [Fe/H], most significantly over -2.0 < [Fe/H] < -1.0, suggestive of incomplete mixing of the ejecta of Type Ia SNe with the ambient medium of each of Carina's generations. Approximate color-magnitude-diagram age estimates are presented for the sample and, together with our chemical abundances, compared with the results of our previous synthetic CMD analysis, which reported the details of Carina's four well-defined populations. We searched for the Na-O anti-correlation universally reported in the Galaxy's globular clusters, and confirm that this phenomenon does not exist in Carina. We also found that one of the 32 stars in our sample has an extremely enhanced lithium abundance -- A(Li)$_{\text{NLTE}}$ = +3.36, consistent with membership of the ~1% group of Li-rich stars in dSph described by Kirby et al.

We investigate the color-magnitude diagram (CMD) of the Carina dwarf spheroidal galaxy using data of Stetson et al. (2011) and synthetic CMDs based on isochrones of Dotter et al. (2008), in terms of the parameters [Fe/H], age, and [alpha/Fe], for the cases when (i) [alpha/Fe] is held constant and (ii) [alpha/Fe] is varied. The data are well described by four basic epochs of star formation, having [Fe/H] = -1.85, -1.5, -1.2, and ~-1.15 and ages ~13, 7, ~3.5, and ~1.5 Gyr, respectively (for [alpha/Fe] = 0.1 (constant [alpha/Fe]) and [alpha/Fe] = 0.2, 0.1, -0.2, -0.2 (variable [alpha/Fe])), with small spreads in [Fe/H] and age of order 0.1 dex and 1 - 3 Gyr. Within an elliptical radius 13.1 arcmin, the mass fractions of the populations, at their times of formation, were (in decreasing age order) 0.34, 0.39, 0.23, and 0.04. This formalism reproduces five observed CMD features (two distinct subgiant branches of old and intermediate-age populations, two younger, main-sequence components, and the small color dispersion on the red giant branch (RGB)). The parameters of the youngest population are less certain than those of the others, and given it is less centrally concentrated it may not be directly related to them. High-resolution spectroscopically analyzed RGB samples appear statistically incomplete compared with those selected using radial velocity, which contain bluer stars comprising ~5 - 10% of the samples. We conjecture these objects may, at least in part, be members of the youngest population. We use the CMD simulations to obtain insight into the population structure of Carina's upper RGB.

In the era of large stellar spectroscopic surveys, there is emphasis on deriving not only stellar abundances but also ages for millions of stars. In the context of Galactic archeology, stellar ages provide a direct probe of the formation history of the Galaxy. We use the stellar evolution code MESA to compute models with atomic diffusion--with and without radiative acceleration--and extra mixing in the surface layers. The extra mixing consists of both density-dependent turbulent mixing and envelope overshoot mixing. Based on these models we argue that it is important to distinguish between initial, bulk abundances (parameters) and current, surface abundances (variables) in the analysis of individual stellar ages. In stars that maintain radiative regions on evolutionary timescales, atomic diffusion modifies the surface abundances. We show that when initial, bulk metallicity is equated with current, surface metallicity in isochrone age analysis the resulting stellar ages can be systematically over-estimated by up to 20%. The change of surface abundances with evolutionary phase also complicates chemical tagging, the concept that dispersed star clusters can be identified through unique, high-dimensional chemical signatures. Stars from the same cluster, but in different evolutionary phases, will show different surface abundances. We speculate that calibration of stellar models may allow us to estimate not only stellar ages but also initial abundances for individual stars. In the meantime, analyzing the chemical properties of stars in similar evolutionary phases is essential to minimize the effects of atomic diffusion in the context of chemical tagging.

As part of a large investigation with Hubble Space Telescope to study the faintest stars within the globular cluster Omega Centauri, in this work we present early results on the multiplicity of its main sequence (MS) stars, based on deep optical and near-infrared observations. By using appropriate color-magnitude diagrams we have identified, for the first time, the two main stellar populations I, and II along the entire MS, from the turn-off towards the hydrogen-burning limit. We have compared the observations with suitable synthetic spectra of MS stars and conclude that the two MSs are consistent with stellar populations with different metallicity, helium, and light-element abundance. Specifically, MS-I corresponds to a metal-poor stellar population ([Fe/H]~-1.7) with Y~ 0.25 and [O/Fe]~0.30. The MS-II hosts helium-rich (Y~0.37-0.40) stars with metallicity ranging from [Fe/H]~-1.7 to -1.4. Below the MS knee (mF160W~19.5, our photometry reveals that each of the two main MSs hosts stellar subpopulations with different oxygen abundances, with very O-poor stars ([O/Fe]~-0.5) populating the MS-II. Such a complexity has never been observed in previous studies of M-dwarfs in globular clusters. A few months before the lunch of the James Webb Space Telescope, these results demonstrate the power of optical and near-infrared photometry in the study of multiple stellar populations in globular clusters.

One of the most unexpected results in the field of stellar populations of the last few years, is the discovery that some Magellanic-Cloud globular clusters younger than ~400 Myr, exhibit bimodal main sequences (MSs) in their color-magnitude diagrams (CMDs). Moreover, these young clusters host an extended main sequence turn off (eMSTO) in close analogy with what is observed in most ~1-2 Gyr old clusters of both Magellanic Clouds. We use high-precision Hubble-Space-Telescope photometry to study the young star cluster NGC1866 in the Large Magellanic Cloud. We discover an eMSTO and a split MS. The analysis of the CMD reveals that (i) the blue MS is the less populous one, hosting about one-third of the total number of MS stars; (ii) red-MS stars are more centrally concentrated than blue-MS stars; (iii) the fraction of blue-MS stars with respect to the total number of MS stars drops by a factor of ~2 in the upper MS with F814W <~19.7. The comparison between the observed CMDs and stellar models reveals that the observations are consistent with ~200 Myr old highly-rotating stars on the red-MS, with rotation close to critical value, plus a non-rotating stellar population spanning an age interval between ~140 and 220 Myr, on the blue-MS. Noticeable, neither stellar populations with different ages only, nor coeval stellar models with different rotation rates, properly reproduce the observed split MS and eMSTO. We discuss these results in the context of the eMSTO and multiple MS phenomenon.

We provide a set of stellar evolution and nucleosynthesis calculations that applies established physics assumptions simultaneously to low- and intermediate-mass and massive star models. Our goal is to provide an internally consistent and comprehensive nuclear production and yield data base for applications in areas such as pre-solar grain studies. Our non-rotating models assume convective boundary mixing where it has been adopted before. We include 8 (12) initial masses for $Z = 0.01$ ($0.02$). Models are followed either until the end of the asymptotic giant branch phase or the end of Si burning, complemented by a simple analytic core-collapse supernova models with two options for fallback and shock velocities. The explosions show which pre-supernova yields will most strongly be effected by the explosive nucleosynthesis. We discuss how these two explosion parameters impacts the light elements and the $s$ and $p$ process. For low- and intermediate-mass models our stellar yields from H to Bi include the effect of convective boundary mixing at the He-intershell boundaries and the stellar evolution feedback of the mixing process that produces the $^{13}$C pocket. All post-processing nucleosynthesis calculations use the same nuclear reaction rate network and nuclear physics input. We provide a discussion of the nuclear production across the entire mass range organized by element group. All our stellar nucleosynthesis profile and time evolution output is available electronically, and tools to explore the data on the NuGrid VOspace hosted by the Canadian Astronomical Data Centre are introduced.

This is the first of a series of papers presenting the Modules for Experiments in Stellar Astrophysics (MESA) Isochrones and Stellar Tracks (MIST) project, a new comprehensive set of stellar evolutionary tracks and isochrones computed using MESA, a state-of-the-art open-source 1D stellar evolution package. In this work, we present models with solar-scaled abundance ratios covering a wide range of ages ($5 \leq \rm \log(Age)\;[yr] \leq 10.3$), masses ($0.1 \leq M/M_{\odot} \leq 300$), and metallicities ($-2.0 \leq \rm [Z/H] \leq 0.5$). The models are self-consistently and continuously evolved from the pre-main sequence to the end of hydrogen burning, the white dwarf cooling sequence, or the end of carbon burning, depending on the initial mass. We also provide a grid of models evolved from the pre-main sequence to the end of core helium burning for $-4.0 \leq \rm [Z/H] < -2.0$. We showcase extensive comparisons with observational constraints as well as with some of the most widely used existing models in the literature. The evolutionary tracks and isochrones can be downloaded from the project website at http://waps.cfa.harvard.edu/MIST/.

We measure chemical abundance ratios and radial velocities in four massive (i.e., young) [alpha/Fe]-rich red giant stars using high-resolution high-S/N spectra from ESPaDOnS fed by Gemini-GRACES. Our differential analysis ensures that our chemical abundances are on the same scale as the Alves-Brito et al. (2010) study of bulge, thin and thick disk red giants. We confirm that the program stars have enhanced [alpha/Fe] ratios and are slightly metal poor. Aside from lithium enrichment in one object, the program stars exhibit no chemical abundance anomalies when compared to giant stars of similar metallicity throughout the Galaxy. This includes the elements Li, O, Si, Ca, Ti, Cr, Ni, Cu, Ba, La, and Eu. Therefore, there are no obvious chemical signatures that can help to reveal the origin of these unusual stars. While our new observations show that only one star (not the Li-rich object) exhibits a radial velocity variation, simulations indicate that we cannot exclude the possibility that all four could be binaries. In addition, we find that two (possibly three) stars show evidence for an infrared excess, indicative of a debris disk. This is consistent with these young [alpha/Fe]-rich stars being evolved blue stragglers, suggesting their apparent young age is a consequence of a merger or mass transfer. We would expect a binary fraction of ~50% or greater for the entire sample of these stars, but the signs of the circumbinary disk may have been lost since these features can have short timescales. Radial velocity monitoring is needed to confirm the blue straggler origin.

The globular cluster 47 Tuc exhibits a complex sub-giant branch (SGB) with a faint-SGB comprising only about the 10% of the cluster mass and a bright-SGB hosting at least two distinct populations.We present a spectroscopic analysis of 62 SGB stars including 21 faint-SGB stars. We thus provide the first chemical analysis of the intriguing faint-SGB population and compare its abundances with those of the dominant populations. We have inferred abundances of Fe, representative light elements C, N, Na, and Al, {\alpha} elements Mg and Si for individual stars. Oxygen has been obtained by co-adding spectra of stars on different sequences. In addition, we have analysed 12 stars along the two main RGBs of 47 Tuc. Our principal results are: (i) star-to-star variations in C/N/Na among RGB and bright-SGB stars; (ii) substantial N and Na enhancements for the minor population corresponding to the faint-SGB; (iii) no high enrichment in C+N+O for faint-SGB stars. Specifically, the C+N+O of the faint-SGB is a factor of 1.1 higher than the bright-SGB, which, considering random (+-1.3) plus systematic errors (+-0.3), means that their C+N+O is consistent within observational uncertainties. However, a small C+N+O enrichment for the faint-SGB, similar to what predicted on theoretical ground, cannot be excluded. The N and Na enrichment of the faint-SGB qualitatively agrees with this population possibly being He-enhanced, as suggested by theory. The iron abundance of the bright and faint-SGB is the same to a level of ~0.10 dex, and no other significant difference for the analysed elements has been detected.

I describe a method to transform a set of stellar evolution tracks onto a uniform basis and then interpolate within that basis to construct stellar isochrones. The method accommodates a broad range of stellar types, from substellar objects to high-mass stars, and phases of evolution, from the pre-main sequence to the white dwarf cooling sequence. I discuss situations in which stellar physics leads to departures from the otherwise monotonic relation between initial stellar mass and lifetime and how these may be dealt with in isochrone construction. I close with convergence tests and recommendations for the number of points in the uniform basis and the mass between tracks in the original grid required in order to achieve a certain level of accuracy in the resulting isochrones. The programs that implement these methods are free and open-source; they may be obtained from the project webpage.

We present an analysis of deep HST/WFC3 near-IR (NIR) imaging data of the globular cluster M4. The best-photometry NIR colour-magnitude diagram (CMD) clearly shows the main sequence extending towards the expected end of the Hydrogen-burning limit and going beyond this point towards fainter sources. The white dwarf sequence can be identified. As such, this is the deepest NIR CMD of a globular cluster to date. Archival HST optical data were used for proper-motion cleaning of the CMD and for distinguishing the white dwarfs (WDs) from brown dwarf (BD) candidates. Detection limits in the NIR are around F110W approx 26.5 mag and F160W approx27 mag, and in the optical around F775W approx 28 mag. Comparing our observed CMDs with theoretical models, we conclude that we have reached beyond the H-burning limit in our NIR CMD and are probably just above or around this limit in our optical-NIR CMDs. Thus, any faint NIR sources that have no optical counterpart are potential BD candidates, since the optical data are not deep enough to detect them. We visually inspected the positions of NIR sources which are fainter than the H-burning limit in F110W and for which the optical photometry did not return a counterpart. We found in total five sources for which we did not get an optical measurement. For four of these five sources, a faint optical counterpart could be visually identified, and an upper optical magnitude was estimated. Based on these upper optical magnitude limits, we conclude that one source is likely a WD, one source could either be a WD or BD candidate, and the remaining two sources agree with being BD candidates. For only one source no optical counterpart could be detected, which makes this source a good BD candidate. We conclude that we found in total four good BD candidates.

We analyse color-magnitude diagrams of eight Globular Clusters (GCs) in the outer Galactic Halo. Images were taken with the Wide Field Channel of the Advanced Camera for Survey and the Ultraviolet and Visual Channel of the Wide Field Camera 3 on board of the Hubble Space Telescope. We have determined the fraction of binary stars along the main sequence and combined results with those of a recent paper where some of us have performed a similar analysis on 59 Galactic GCs. In total, binaries have been now studied homogeneously in 67 GCs. We studied the radial and luminosity distributions of the binary systems, the distribution of their mass-ratios and investigated univariate relations with several parameters of the host GCs. We confirm the anti-correlation between the binary fraction and the luminosity of the host cluster, and find that low-luminosity clusters can host a large population in excess of ~40% in the cluster core. However, our results do not support a significant correlation with the cluster age as suggested in the literature. In most GCs, binaries are more centrally concentrated than single stars. If the fraction of binaries is normalised to the core binary fraction the radial density profiles follow a common trend. It has a maximum in the center and declines by a factor of two at a distance of about two core radii from the cluster center. After dropping to its minimum at a radial distance of $\sim$5 core radii it stays approximately constant at larger radii. We also find that the mass-ratio and the distribution of binaries as a function of the mass of the primary star is almost flat.

We use photometric and spectroscopic observations of the detached eclipsing binaries V40 and V41 in the globular cluster NGC 6362 to derive masses, radii, and luminosities of the component stars. The orbital periods of these systems are 5.30 and 17.89 d, respectively. The measured masses of the primary and secondary components ($M_p$, $M_s$) are (0.8337$\pm$0.0063, 0.7947$\pm$0.0048) M$_\odot$ for V40 and (0.8215$\pm$0.0058, 0.7280$\pm$0.0047) M$_\odot$ for V41. The measured radii ($R_p$, $R_s$) are (1.3253$\pm$0.0075, 0.997$\pm$0.013) R$_\odot$ for V40 and (1.0739$\pm$0.0048, 0.7307$\pm$0.0046) R$_\odot$ for V41. Based on the derived luminosities, we find that the distance modulus of the cluster is 14.74$\pm$0.04 mag -- in good agreement with 14.72 mag obtained from CMD fitting. We compare the absolute parameters of component stars with theoretical isochrones in mass-radius and mass-luminosity diagrams. For assumed abundances [Fe/H] = -1.07, [$\alpha$/Fe] = 0.4, and Y = 0.25 we find the most probable age of V40 to be 11.7$\pm$0.2 Gyr, compatible with the age of the cluster derived from CMD fitting (12.5$\pm$0.5 Gyr). V41 seems to be markedly younger than V40. If independently confirmed, this result will suggest that V41 belongs to the younger of the two stellar populations recently discovered in NGC 6362. The orbits of both systems are eccentric. Given the orbital period and age of V40, its orbit should have been tidally circularized some $\sim$7 Gyr ago. The observed eccentricity is most likely the result of a relatively recent close stellar encounter.

We use images acquired with the Hubble Space Telescope Wide Field Camera 3 and new models to probe the Horizontal Branch (HB) population of the We use images acquired with the Hubble Space Telescope Wide Field Camera 3 and new models to probe the horizontal branch (HB) population of the Galactic globular cluster (GC) NGC 2419. A detailed analysis of the composite HB highlights three populations:(1) the blue luminous HB, hosting standard helium stars (Y=0.25) with a very small spread of mass, (2) a small population of stars with intermediate helium content (0.26<Y<=0.29), and (3) the well-populated extreme HB. We can fit the last group with models having high helium abundance (Y \sim 0.36), half of which (the hottest part, 'blue hook' stars) are identified as possible 'late flash mixed stars'. The initial helium abundance of this extreme population is in nice agreement with the predicted helium abundance in the ejecta of massive asymptotic giant branch (AGB) stars of the same metallicity as NGC 2419. This result further supports the hypothesis that second-generation stars in GCs formed from the ashes of intermediate-mass AGB stars. We find that the distribution in magnitude of the blue hook stars is larger than that predicted by theoretical models. We discuss the possible uncertainties in the magnitude scales and different attempts to model this group of stars. Finally, we suggest that consistency can be better achieved if we assume core masses larger than predicted by our models. This may be possible if the progenitors were fast rotators on the main sequence. If further study confirms this interpretation, a fast initial rotation would be a strong signature of the peculiarity of extreme second-generation stars in GCs.

Horizontal Branch stars belong to an advanced stage in the evolution of the oldest stellar galactic population, occurring either as field halo stars or grouped in globular clusters. The discovery of multiple populations in these clusters, that were previously believed to have single populations gave rise to the currently accepted theory that the hottest horizontal branch members (the blue hook stars, which had late helium-core flash ignition, followed by deep mixing) are the progeny of a helium-rich "second generation" of stars. It is not known why such a supposedly rare event (a late flash followed by mixing) is so common that the blue hook of {\omega} Cen contains \sim 30% of horizontal branch stars 10 , or why the blue hook luminosity range in this massive cluster cannot be reproduced by models. Here we report that the presence of helium core masses up to \sim 0.04 solar masses larger than the core mass resulting from evolution is required to solve the luminosity range problem. We model this by taking into account the dispersion in rotation rates achieved by the progenitors, whose premain sequence accretion disc suffered an early disruption in the dense environment of the cluster's central regions where second-generation stars form. Rotation may also account for frequent late-flash-mixing events in massive globular clusters.

Knowledge of the chemical composition and absolute masses of Capella are key to understanding the evolutionary state of this benchmark binary system comprising two giant stars. Previous efforts, including our own 2009 study, have largely failed to reach an acceptable agreement between the observations and current stellar evolution models, preventing us from assessing the status of the primary. Here we report a revision of the physical properties of the components incorporating recently published high-precision radial velocity measurements, and a new detailed chemical analysis providing abundances for more than 20 elements in both stars. We obtain highly precise (to about 0.3%) masses of 2.5687 +/- 0.0074 and 2.4828 +/- 0.0067 solar masses, radii of 11.98 +/- 0.57 and 8.83 +/- 0.33 solar radii, effective temperatures of 4970 +/- 50 K and 5730 +/- 60 K, and independently measured luminosities based on the orbital parallax (78.7 +/- 4.2 and 72.7 +/- 3.6 solar luminosities). We find an excellent match to stellar evolution models at the measured composition of [Fe/H] = -0.04 +/- 0.06. Three different sets of models place the primary star firmly at the end of the core helium-burning phase (clump), while the secondary is known to be evolving rapidly across the Hertzprung gap. The measured lithium abundance, the C/N ratio, and the 12C/13C isotopic carbon abundance ratio, which change rapidly in the giant phase, are broadly in agreement with expectations from models. Predictions from tidal theory for the spin rates, spin-orbit alignment, and other properties do not fare as well, requiring a 40-fold increase in the efficiency of the dissipation mechanisms in order to match the observations.

The GALAH survey is a large high-resolution spectroscopic survey using the newly commissioned HERMES spectrograph on the Anglo-Australian Telescope. The HERMES spectrograph provides high-resolution (R ~28,000) spectra in four passbands for 392 stars simultaneously over a 2 degree field of view. The goal of the survey is to unravel the formation and evolutionary history of the Milky Way, using fossil remnants of ancient star formation events which have been disrupted and are now dispersed throughout the Galaxy. Chemical tagging seeks to identify such dispersed remnants solely from their common and unique chemical signatures; these groups are unidentifiable from their spatial, photometric or kinematic properties. To carry out chemical tagging, the GALAH survey will acquire spectra for a million stars down to V~14. The HERMES spectra of FGK stars contain absorption lines from 29 elements including light proton-capture elements, alpha-elements, odd-Z elements, iron-peak elements and n-capture elements from the light and heavy s-process and the r-process. This paper describes the motivation and planned execution of the GALAH survey, and presents some results on the first-light performance of HERMES.

Detailed chemical abundances are presented for seven M31 outer halo globular clusters (with projected distances from M31 greater than 30 kpc), as derived from high resolution integrated light spectra taken with the Hobby Eberly Telescope. Five of these clusters were recently discovered in the Pan-Andromeda Archaeological Survey (PAndAS)---this paper presents the first determinations of integrated Fe, Na, Mg, Ca, Ti, Ni, Ba, and Eu abundances for these clusters. Four of the target clusters (PA06, PA53, PA54, and PA56) are metal-poor ([Fe/H] < -1.5), alpha-enhanced (though they are possibly less alpha-enhanced than Milky Way stars at the 1 sigma level), and show signs of star-to-star Na and Mg variations. The other three globular clusters (H10, H23, and PA17) are more metal rich, with metallicities ranging from [Fe/H] = -1.4 to -0.9. While H23 is chemically similar to Milky Way field stars, Milky Way globular clusters, and other M31 clusters, H10 and PA17 have moderately low [Ca/Fe], compared to Milky Way field stars and clusters. Additionally, PA17's high [Mg/Ca] and [Ba/Eu] ratios are distinct from Milky Way stars, and are in better agreement with the stars and clusters in the Large Magellanic Cloud (LMC). None of the clusters studied here can be conclusively linked to any of the identified streams from PAndAS; however, based on their locations, kinematics, metallicities, and detailed abundances, the most metal-rich PAndAS clusters H23 and PA17 may be associated with the progenitor of the Giant Stellar Stream, H10 may be associated with the SW Cloud, and PA53 and PA56 may be associated with the Eastern Cloud.

We present stellar atmosphere and evolution models of main sequence stars in two stellar populations of the Galactic globular cluster NGC 6752. These populations represent the two extremes of light-element abundance variations in the cluster. NGC 6752 is a benchmark cluster in the study of multiple stellar populations because of the rich array of spectroscopic abundances and panchromatic Hubble Space Telescope photometry. The spectroscopic abundances are used to compute stellar atmosphere and evolution models. The synthetic spectra for the two populations show significant differences in the ultraviolet and, for the coolest temperatures, in the near-infrared. The stellar evolution models exhibit insignificant differences in the H-R diagram except on the lower main sequence. The appearance of multiple sequences in the colour-magnitude diagrams (CMDs) of NGC 6752 is almost exclusively due to spectral effects caused by the abundance variations. The models reproduce the observed splitting and/or broadening of sequences in a range of CMDs. The ultraviolet CMDs are sensitive to variations in carbon, nitrogen, and oxygen but the models are not reliable enough to directly estimate abundance variations from photometry. On the other hand, the widening of the lower main sequence in the near-infrared CMD, driven by oxygen-variation via the water molecule, is well-described by the models and can be used to estimate the range of oxygen present in a cluster from photometry. We confirm that it is possible to use multiband photometry to estimate helium variations among the different populations, with the caveat that the estimated amount of helium-enhancement is model-dependent.

We present the first detailed analysis of the detached eclipsing binary V15 in the super-metal rich open cluster NGC 6253. We obtain the following absolute parameters: M_p=1.303+-0.006 Msun, R_p=1.71+-0.03 Rsun, L_p=2.98+-0.10 Lsun for the primary, and M_s=1.225+-0.006 Msun, R_s=1.44+-0.02 Rsun, L_s=2.13+-0.06 Lsun for the secondary. Based on Dartmouth isochrones, the age of NGC 6253 is estimated to be 3.80 - 4.25 Gyr from the mass-radius diagram and 3.9 - 4.6 Gyr from color-magnitude diagram (CMD) fitting. Both of these estimates are significantly higher than those reported so far. The derived apparent distance modulus of 11.65 mag agrees well with the range of 10.9 - 12.2 mag derived by other authors; however our estimated reddening (0.113 mag) is lower than the lowest published value (0.15 mag). We confirm earlier observations that model atmospheres are not accurate enough to account for the whole CMD of the cluster, with the largest discrepancies appearing on the subgiant and giant branches. Although age estimation from the mass-radius diagram is a relatively safe, distance- and reddening-independent procedure, our results should be verified by photometric and spectroscopic observations of additional detached eclipsing binaries which we have discovered, at least two of which are proper-motion members of NGC 6253.

High resolution integrated light (IL) spectroscopy provides detailed abundances of distant globular clusters whose stars cannot be resolved. Abundance comparisons with other systems (e.g. for chemical tagging) require understanding the systematic offsets that can occur between clusters, such as those due to uncertainties in the underlying stellar population. This paper analyses high resolution IL spectra of the Galactic globular clusters 47 Tuc, M3, M13, NGC 7006, and M15 to (1) quantify potential systematic uncertainties in Fe, Ca, Ti, Ni, Ba, and Eu and (2) identify the most stable abundance ratios that will be useful in future analyses of unresolved targets. When stellar populations are well-modelled, uncertainties are ~0.1-0.2 dex based on sensitivities to the atmospheric parameters alone; in the worst case scenarios, uncertainties can rise to 0.2-0.4 dex. The [Ca I/Fe I] ratio is identified as the optimal integrated [alpha/Fe] indicator (with offsets <0.1 dex), while [Ni I/Fe I] is also extremely stable to within <0.1 dex. The [Ba II/Eu II] ratios are also stable when the underlying populations are well modelled and may also be useful for chemical tagging.

NGC1851 is surrounded by a stellar component that extends more than ten times beyond the tidal radius. Although the nature of this stellar structure is not known, it has been suggested to be a sparse halo of stars or associated with a stellar stream. We analyse the nature of this intriguing stellar component surrounding NGC1851 by investigating its radial velocities and chemical composition, in particular in comparison with those of the central cluster analysed in a homogeneous manner. In total we observed 23 stars in the halo with radial velocities consistent with NGC1851, and for 15 of them we infer [Fe/H] abundances. Our results show that: (i) stars dynamically linked to NGC1851 are present at least up to ~2.5 tidal radii, supporting the presence of a halo of stars surrounding the cluster; (ii) apart from the NGC1851 radial velocity-like stars, our observed velocity distribution agrees with that expected from Galactic models, suggesting that no other sub-structure (such as a stream) at different radial velocities is present in our field; (iii) the chemical abundances for the s-process elements Sr and Ba are consistent with the s-normal stars observed in NGC1851; (iv) all halo stars have metallicities, and abundances for the other studied elements Ca, Mg and Cr, consistent with those exhibited by the cluster. The complexity of the whole NGC1851 cluster+halo system may agree with the scenario of a tidally-disrupted dwarf galaxy in which NGC1851 was originally embedded.

We present an analysis of the detached eclipsing binaries V44 and V54 belonging to the globular cluster M55. For V54 we obtain the following absolute parameters: $M_p=0.726\pm 0.015\,M_\odot$, $R_p=1.006\pm 0.009\,R_\odot$, $L_p=1.38\pm 0.07\,L_\odot$ for the primary, and $M_s=0.555\pm 0.008\,M_\odot$, $R_s=0.528\pm 0.005\,R_\odot$, $L_s=0.16\pm0.01\,L_\odot$ for the secondary. The age and apparent distance modulus of V54 are estimated at 13.3 - 14.7 Gyr and $13.94\pm0.05$ mag, respectively. This age is substantially larger than ages we have derived from the analysis of binary systems in 47 Tuc and M4. The secondary of V44 is so weak in the optical domain that only mass function and relative parameters are obtained for the components of this system. However, there is a good chance that the velocity curve of the secondary could be derived from near-IR spectra. As the primary of V44 is more evolved than that of V54, such data would impose much tighter limits on the age and distance of M55.

Asteroseismology has the capability of precisely determining stellar properties which would otherwise be inaccessible, such as radii, masses and thus ages of stars. When coupling this information with classical determinations of stellar parameters, such as metallicities, effective temperatures and angular diameters, powerful new diagnostics for Galactic studies can be obtained. The ongoing Stroemgren survey for Asteroseismology and Galactic Archaeology (SAGA) has the goal of transforming the Kepler field into a new benchmark for Galactic studies, similarly to the solar neighborhood. Here we present first results from a stripe centred at Galactic longitude 74deg and covering latitude from about 8 to 20deg, which includes almost 1000 K-giants with seismic information and the benchmark open cluster NGC 6819. We describe the coupling of classical and seismic parameters, the accuracy as well as the caveats of the derived effective temperatures, metallicities, distances, surface gravities, masses, and radii. Confidence in the achieved precision is corroborated by the detection of the first and secondary clump in a population of field stars with a ratio of 2 to 1, and by the negligible scatter in the seismic distances among NGC 6819 member stars. An assessment of the reliability of stellar parameters in the Kepler Input Catalogue is also performed, and the impact of our results for population studies in the Milky Way is discussed, along with the importance of an all-sky Stroemgren survey.

We present new Stromgren and Washington data sets for the Bootes I dwarf galaxy, and combine them with the available SDSS photometry. The goal of this project is to refine a ground-based, practical, accurate method to determine age and metallicity for individual stars in Bootes I that can be selected in an unbiased imaging survey, without having to take spectra. We produce photometric metallicities from Stromgren and Washington photometry, for stellar systems with a range of $-1.0>[Fe/H]>-3.5$. To avoid the decrease in sensitivity of the Stromgren metallicity index on the lower red-giant branch, we replace the Stromgren v-filter with the broader Washington C-filter; we find that $CT_1by$ is the most successful filter combination, for individual stars with $[Fe/H]<-2.0$, to maintain ~0.2 dex $[Fe/H]$-resolution over the whole red-giant branch. We demonstrate that we can break the isochrones' age-metallicity degeneracy with these filters, using stars with log g=2.5-3.0, which have less than a 2% change in their $(C-T_1)$-colour due to age, over a range of 11-14 Gyr.

The M4 Core Project with HST is designed to exploit the Hubble Space Telescope to investigate the central regions of M4, the Globular Cluster closest to the Sun. In this paper we combine optical and near-infrared photometry to study multiple stellar populations in M4. We detected two sequences of M-dwarfs containing ~38% (MS_I) and ~62% (MS_II) of MS stars below the main-sequence (MS) knee. We compare our observations with those of NGC2808, which is the only other GCs where multiple MSs of very low-mass stars have been studied to date. We calculate synthetic spectra for M-dwarfs, assuming the chemical composition mixture inferred from spectroscopic studies of stellar populations along the red giant branch, and different Helium abundances, and we compare predicted and observed colors. Observations are consistent with two populations, one with primordial abundance and another with enhanced nitrogen and depleted oxygen.

The horizontal branch (HB) morphology of globular clusters (GCs) is mainly determined by metallicity. However, the fact that GCs with almost the same metallicity exhibit different HB morphologies demonstrates that at least one more parameter is needed to explain the HB morphology. It has been suggested that one of these should be a global parameter that varies from GC to GC, and the other a non-global parameter that varies within the GC. In this study we provide empirical evidence corroborating this idea. We used the photometric catalogs obtained with the Advanced Camera for Surveys (ACS) of the Hubble Space Telescope (HST) and analyse the CMDs of 74 GCs. The HB morphology of our sample of GCs has been investigated on the basis of the two new parameters L1 and L2 that measure the distance between the RGB and the coolest part of the HB, and the color extension of the HB, respectively. We find that L1 correlates with both metallicity and age, whereas L2 most strongly correlates with the mass of the hosting GC. The range of helium abundance among the stars in a GC, characterised by Delta Y and associated with the presence of multiple stellar populations, has been estimated in a few GCs to date. In these GCs we find a close relationship among Delta Y, GC mass, and L2. We conclude that age and metallicity are the main global parameters while the range of helium abundance within a GC is the main non-global parameter defining the HB morphology of Galactic GCs.

Using deep Hubble Space Telescope imaging, color-magnitude diagrams are constructed for the globular clusters 47 Tuc and NGC 6397. As expected, because of its lower metal abundance, the main sequence of NGC 6397 lies well to the blue of that of 47 Tuc. A comparison of the white dwarf cooling sequences of the two clusters, however, demonstrates that these sequences are indistinguishable over most of their loci - a consequence of the settling out of heavy elements in the dense white dwarf atmosphere and the near equality of their masses. Lower quality data on M4 continues this trend to a third cluster whose metallicity is intermediate between these two. While the path of the white dwarfs in the color-magnitude diagram is nearly identical in 47 Tuc and NGC 6397, the numbers of white dwarfs along the path are not. This results from the relatively rapid relaxation in NGC 6397 compared to 47 Tuc and provides a cautionary note that simply counting objects in star clusters in random locations as a method of testing stellar evolutionary theory is likely dangerous unless dynamical considerations are included.

Globular clusters trace the formation history of the spheroidal components of both our Galaxy and others, which represent the bulk of star formation over the history of the universe. They also exhibit a range of metallicities, with metal-poor clusters dominating the stellar halo of the Galaxy, and higher metallicity clusters found within the inner Galaxy, associated with the stellar bulge, or the thick disk. Age differences between these clusters can indicate the sequence in which the components of the Galaxy formed, and in particular which clusters were formed outside the Galaxy and later swallowed along with their original host galaxies, and which were formed in situ. Here we present an age determination of the metal-rich globular cluster 47 Tucanae by fitting the properties of the cluster white dwarf population, which implies an absolute age of 9.9 (0.7) Gyr at 95% confidence. This is about 2.0 Gyr younger than inferred for the metal-poor cluster NGC 6397 from the same models, and provides quantitative evidence that metal-rich clusters like 47 Tucanae formed later than the metal-poor halo clusters like NGC 6397.

The Second Parameter (2ndP) Problem recognizes the remarkable role played by horizontal branch (HB) morphology in the development of our understanding of globular clusters, and the Galaxy, over the last 50 years. I will describe the historical development of the 2ndP and discuss recent advances that are finally providing some answers. I will discuss how the controversies surrounding the nature of the 2ndP can be reconciled if we acknowledge that there are actually two distinct problems with entirely different solutions.

We analyze extensive BVR_cI_c time-series photometry and radial-velocity measurements for WOCS 40007 (Auner 259; KIC 5113053), a double-lined detached eclipsing binary and a member of the open cluster NGC 6819. Utilizing photometric observations from the 1-meter telescope at Mount Laguna Observatory and spectra from the WIYN 3.5-meter telescope, we measure precise and accurate masses (~1.6% uncertainty) and radii (~0.5%) for the binary components. In addition, we discover a third star orbiting the binary with a period greater than 3000 days using radial velocities and Kepler eclipse timings. Because the stars in the eclipsing binary are near the cluster turnoff, they are evolving rapidly in size and are sensitive to age. With a metallicity of [Fe/H]=+0.09+/-0.03, we find the age of NGC 6819 to be about 2.4 Gyr from CMD isochrone fitting and 3.1+/-0.4 Gyr by analyzing the mass-radius (M-R) data for this binary. The M-R age is above previous determinations for this cluster, but consistent within 1 sigma uncertainties. When the M-R data for the primary star of the additional cluster binary WOCS 23009 is included, the weighted age estimate drops to 2.5+/-0.2 Gyr, with a systematic uncertainty of at least 0.2 Gyr. The age difference between our CMD and M-R findings may be the result of systematic error in the metallicity or helium abundance used in models, or due to slight radius inflation of one or both stars in the WOCS 40007 binary.

We substantially update the capabilities of the open source software package Modules for Experiments in Stellar Astrophysics (MESA), and its one-dimensional stellar evolution module, MESA Star. Improvements in MESA Star's ability to model the evolution of giant planets now extends its applicability down to masses as low as one-tenth that of Jupiter. The dramatic improvement in asteroseismology enabled by the space-based Kepler and CoRoT missions motivates our full coupling of the ADIPLS adiabatic pulsation code with MESA Star. This also motivates a numerical recasting of the Ledoux criterion that is more easily implemented when many nuclei are present at non-negligible abundances. This impacts the way in which MESA Star calculates semi-convective and thermohaline mixing. We exhibit the evolution of 3-8 Msun stars through the end of core He burning, the onset of He thermal pulses, and arrival on the white dwarf cooling sequence. We implement diffusion of angular momentum and chemical abundances that enable calculations of rotating-star models, which we compare thoroughly with earlier work. We introduce a new treatment of radiation-dominated envelopes that allows the uninterrupted evolution of massive stars to core collapse. This enables the generation of new sets of supernovae, long gamma-ray burst, and pair-instability progenitor models. We substantially modify the way in which MESA Star solves the fully coupled stellar structure and composition equations, and we show how this has improved MESA's performance scaling on multi-core processors. Updates to the modules for equation of state, opacity, nuclear reaction rates, and atmospheric boundary conditions are also provided. We describe the MESA Software Development Kit (SDK) that packages all the required components needed to form a unified and maintained build environment for MESA. [Abridged]

Spectrum syntheses for three elements (Mg, Na, and Eu) in high-resolution integrated light spectra of the Galactic globular clusters 47 Tuc, M3, M13, NGC 7006, and M15 are presented, along with calibration syntheses of the Solar and Arcturus spectra. Iron abundances in the target clusters are also derived from integrated light equivalent width analyses. Line profiles in the spectra of these five globular clusters are well fit after careful consideration of the atomic and molecular spectral features, providing levels of precision that are better than equivalent width analyses of the same integrated light spectra, and that are comparable to the precision in individual stellar analyses. The integrated light abundances from the 5528 and 5711 A Mg I lines, the 6154 and 6160 A Na I lines, and the 6645 A Eu II line fall within the observed ranges from individual stars; however, these integrated light abundances do not always agree with the average literature abundances. Tests with the second parameter clusters M3, M13, and NGC 7006 show that assuming an incorrect horizontal branch morphology is likely to have only a small (< 0.06 dex) effect on these Mg, Na, and Eu abundances. These tests therefore show that integrated light spectrum syntheses can be applied to unresolved globular clusters over a wide range of metallicities and horizontal branch morphologies. Such high precision in integrated light spectrum syntheses is valuable for interpreting the chemical abundances of globular cluster systems around other galaxies.

Based on the width of its main sequence, and an actual observed split when viewed through particular filters, it is widely accepted that 47 Tucanae contains multiple stellar populations. In this contribution, we divide the main-sequence of 47 Tuc into four color groups, which presumably represent stars of various chemical compositions. The kinematic properties of each of these groups is explored via proper-motions, and a strong signal emerges of differing proper-motion anisotropies with differing main-sequence color; the bluest main-sequence stars exhibit the largest proper-motion anisotropy which becomes undetectable for the reddest stars. In addition, the bluest stars are also the most centrally concentrated. A simiilar analysis for SMC stars, which are located in the background of 47 Tuc on our frames, yields none of the anisotropy exhibited by the 47 Tuc stars. We discuss implications of these results for possible formation scenarios of the various populations.

We present Hubble Space Telescope imaging of a newly-discovered faint stellar system, PAndAS-48, in the outskirts of the M31 halo. Our photometry reveals this object to be comprised of an ancient and very metal-poor stellar population with age > 10 Gyr and [Fe/H] < -2.3. Our inferred distance modulus of 24.57 +/- 0.11 confirms that PAndAS-48 is most likely a remote M31 satellite with a 3D galactocentric radius of 149 (+19 -8) kpc. We observe an apparent spread in color on the upper red giant branch that is larger than the photometric uncertainties should allow, and briefly explore the implications of this. Structurally, PAndAS-48 is diffuse, faint, and moderately flattened, with a half-light radius rh = 26 (+4 -3) pc, integrated luminosity Mv = -4.8 +/- 0.5, and ellipticity = 0.30 (+0.08 -0.15). On the size-luminosity plane it falls between the extended globular clusters seen in several nearby galaxies, and the recently-discovered faint dwarf satellites of the Milky Way; however, its characteristics do not allow us to unambiguously class it as either type of system. If PAndAS-48 is a globular cluster then it is the among the most elliptical, isolated, and metal-poor of any seen in the Local Group, extended or otherwise. Conversely, while its properties are generally consistent with those observed for the faint Milky Way dwarfs, it would be a factor ~2-3 smaller in spatial extent than any known counterpart of comparable luminosity.

We propose a novel method for determining the ages of low-mass, pre-main sequence stellar systems using the apsidal motion of low-mass detached eclipsing binaries. The apsidal motion of a binary system with an eccentric orbit provides information regarding the interior structure constants of the individual stars. These constants are related to the normalized stellar interior density distribution and can be extracted from the predictions of stellar evolution models. We demonstrate that low-mass, pre-main sequence stars undergoing radiative core contraction display rapidly changing interior structure constants (greater than 5% per 10 Myr) that, when combined with observational determinations of the interior structure constants (with 5 -- 10% precision), allow for a robust age estimate. This age estimate, unlike those based on surface quantities, is largely insensitive to the surface layer where effects of magnetic activity are likely to be most pronounced. On the main sequence, where age sensitivity is minimal, the interior structure constants provide a valuable test of the physics used in stellar structure models of low-mass stars. There are currently no known systems where this technique is applicable. Nevertheless, the emphasis on time domain astronomy with current missions, such as Kepler, and future missions, such as LSST, has the potential to discover systems where the proposed method will be observationally feasible.

We use photometric and spectroscopic observations of the eclipsing binaries V65, V66 and V69 in the field of the globular cluster M4 to derive masses, radii, and luminosities of their components. The orbital periods of these systems are 2.29, 8.11 and 48.19 d, respectively. The measured masses of the primary and secondary components (Mp and Ms) are 0.8035+-0.0086 and 0.6050+-0.0044 Msun for V65, 0.7842+-0.0045 and 0.7443+-0.0042 Msun for V66, and 0.7665+-0.0053 and 0.7278+-0.0048 Msun for V69. The measured radii (Rp and Rs) are 1.147+_0.010 and 0.6110+-0.0092 Rsun for V66, 0.9347+_0.0048 and 0.8298+-0.0053 Rsun for V66, and 0.8655+-0.0097 and 0.8074+-0.0080 Rsun for V69. The orbits of V65 and V66 are circular, whereas that of V69 has an eccentricity of 0.38. Based on systemic velocities and relative proper motions, we show that all the three systems are members of the cluster. We find that the distance to M4 is 1.82+-0.04 kpc - in good agreement with recent estimates based on entirely different methods. We compare the absolute parameters of V66 and V69 with two sets of theoretical isochrones in mass-radius and mass-luminosity diagrams, and for an assumed [Fe/H] = -1.20, [alpha/Fe] = 0.4, and Y = 0.25 we find the most probable age of M4 to be between 11.2 and 11.3 Gyr. CMD-fitting with the same parameters yields an age close to, or slightly in excess of, 12 Gyr. However, considering the sources of uncertainty involved in CMD fitting, these two methods of age determination are not discrepant. Age and distance determinations can be further improved when infrared eclipse photometry is obtained.

We present a new measurement of the stellar initial mass function (IMF) based on ultra-deep, high-resolution photometry of >5,000 stars in the outskirts of the Small Magellanic Cloud (SMC) galaxy. The Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS) observations reveal this rich, co-spatial population behind the foreground globular cluster 47 Tuc, which we targeted for 121 HST orbits. The stellar main sequence of the SMC is measured in the F606W, F814W color-magnitude diagram (CMD) down to ~30th magnitude, and is cleanly separated from the foreground star cluster population using proper motions. We simulate the SMC population by extracting stellar masses (single and unresolved binaries) from specific IMFs, and converting those masses to luminosities in our bandpasses. The corresponding photometry for these simulated stars is drawn directly from a rich cloud of 4 million artificial stars, thereby accounting for the real photometric scatter and completeness of the data. Over a continuous and well populated mass range of M = 0.37 - 0.93 Msun (i.e., down to a ~75% completeness limit at F606W = 28.7), we demonstrate that the IMF is well represented by a single power-law form with slope \alpha = -1.90 (^{+0.15}_{-0.10}) (3 sigma error) (i.e., dN/dM \propto M^{\alpha}). This is shallower than the Salpeter slope of \alpha = -2.35, which agrees with the observed stellar luminosity function at higher masses. Our results indicate that the IMF does {\it not} turn over to a more shallow power-law form within this mass range. We discuss implications of this result for the theory of star formation, the inferred masses of galaxies, and the (lack of a) variation of the IMF with metallicity.

Multi-epoch observations with ACS on HST provide a unique and comprehensive probe of stellar dynamics within NGC 6397. We are able to confront analytic models of the globular cluster with the observed stellar proper motions. The measured proper motions probe well along the main sequence from 0.8 to below 0.1 M$_\odot$ as well as white dwarfs younger than one gigayear. The observed field lies just beyond the half-light radius where standard models of globular cluster dynamics (e.g. based on a lowered Maxwellian phase-space distribution) make very robust predictions for the stellar proper motions as a function of mass. The observed proper motions show no evidence for anisotropy in the velocity distribution; furthermore, the observations agree in detail with a straightforward model of the stellar distribution function. We do not find any evidence that the young white dwarfs have received a natal kick in contradiction with earlier results. Using the observed proper motions of the main-sequence stars, we obtain a kinematic estimate of the distance to NGC 6397 of $2.2^{+0.5}_{-0.7}$ kpc and a mass of the cluster of $1.1 \pm 0.1 \times 10^5 \mathrm{M}_\odot$ at the photometric distance of 2.53 kpc. One of the main-sequence stars appears to travel on a trajectory that will escape the cluster, yielding an estimate of the evaporation timescale, over which the number of stars in the cluster decreases by a factor of e, of about 3 Gyr. The proper motions of the youngest white dwarfs appear to resemble those of the most massive main-sequence stars, providing the first direct constraint on the relaxation time of the stars in a globular cluster of greater than or about 0.7 Gyr.

We present an empirical determination of the white dwarf cooling sequence in the globular cluster 47 Tucanae. Using spectral models, we determine temperatures for 887 objects from Wide Field Camera 3 data, as well as 292 objects from data taken with the Advanced Camera for Surveys. We make the assumption that the rate of white dwarf formation in the cluster is constant. Stellar evolution models are then used to determine the rate at which objects are leaving the main sequence, which must be the same as the rate at which objects are arriving on the white dwarf sequence in our field. The result is an empirically derived relation between temperature ($T_{eff}$) and time ($t$) on the white dwarf cooling sequence. Comparing this result to theoretical cooling models, we find general agreement with the expected slopes between 20,000K and 30,000K and between 6,000K and 20,000K, but the transition to the Mestel cooling rate of $T_{eff} \propto t^{-0.4}$ is found to occur at hotter temperatures, and more abruptly than is predicted by any of these models.

By means of high-resolution cosmological hydrodynamical simulations of Milky Way-like disc galaxies, we conduct an analysis of the associated stellar metallicity distribution functions (MDFs). After undertaking a kinematic decomposition of each simulation into spheroid and disc sub-components, we compare the predicted MDFs to those observed in the solar neighbourhood and the Galactic bulge. The effects of the star formation density threshold are visible in the star formation histories, which show a modulation in their behaviour driven by the threshold. The derived MDFs show median metallicities lower by 0.2-0.3 dex than the MDF observed locally in the disc and in the Galactic bulge. Possible reasons for this apparent discrepancy include the use of low stellar yields and/or centrally-concentrated star formation. The dispersions are larger than the one of the observed MDF; this could be due to simulated discs being kinematically hotter relative to the Milky Way. The fraction of low metallicity stars is largely overestimated, visible from the more negatively skewed MDF with respect to the observational sample. For our fiducial Milky Way analog, we study the metallicity distribution of the stars born "in situ" relative to those formed via accretion (from disrupted satellites), and demonstrate that this low-metallicity tail to the MDF is populated primarily by accreted stars. Enhanced supernova and stellar radiation energy feedback to the surrounding interstellar media of these pre-disrupted satellites is suggested as an important regulator of the MDF skewness.

Recent work has shown that most globular clusters have at least two chemically distinct components, as well as cluster-to-cluster differences in the mean [O/Fe], [Mg/Fe], and [Si/Fe] ratios at similar [Fe/H] values. In order to investigate the implications of variations in the abundances of these and other metals for H-R diagrams and predicted ages, grids of evolutionary sequences have been computed for scaled solar and enhanced alpha-element mixtures, and for mixtures in which the assumed [m/Fe] value for each of the metals C, N, O, Ne, Na, Mg, Si, S, Ca, and Ti has been increased, in turn, by 0.4 dex at constant [Fe/H]. These tracks, together with isochrones for ages from 6 to 14 Gyr, have been computed for -3.0 < [Fe/H] < -0.6, with helium abundances Y = 0.25, 0.29, and 0.33 at each [Fe/H] value, using upgraded versions of the Victoria stellar structure program and the Regina interpolation code, respectively. Turnoff luminosity versus age relations from isochrones are found to depend almost entirely on the importance of the CNO-cycle, and thereby mainly on the abundance of oxygen. Since C, N, and O, as well as Ne and S, do not contribute significantly to the opacities at low temperatures and densities, variations in their abundances do not impact the Teff scale of red giants. The latter is a strong function of the abundances of only Mg and Si (and Fe, possibly to a lesser extent), because they are so abundant and because they are strong sources of opacity at low temperatures. For these reasons, Mg and Si also have important effects on the temperatures of main-sequence stars. Due to their low abundances, Na, Ca, and Ti are of little consequence for stellar models. The effects of varying the adopted solar metals mix and the helium abundance at a fixed [Fe/H] are also briefly discussed.

Models of stellar structure and evolution can be constrained by measuring accurate parameters of detached eclipsing binaries in open clusters. Multiple binary stars provide the means to determine helium abundances in these old stellar systems, and in turn, to improve estimates of their age. In the first paper of this series, we demonstrated how measurements of multiple eclipsing binaries in the old open cluster NGC6791 sets tighter constraints on the properties of stellar models than has previously been possible, thereby potentially improving both the accuracy and precision of the cluster age. Here we add additional constraints and perform an extensive model comparison to determine the best estimates of the cluster age and helium content, employing as many observational constraints as possible. We improve our photometry and correct empirically for differential reddening effects. We then perform an extensive comparison of the CMDs and eclipsing binary measurements to Victoria and DSEP isochrones to estimate cluster parameters. We also reanalyse a spectrum of the star 2-17 to improve [Fe/H] constraints. We find a best estimate of the age of ~8.3 Gyr while demonstrating that remaining age uncertainty is dominated by uncertainties in the CNO abundances. The helium mass fraction is well constrained at Y = 0.30 \pm 0.01 resulting in dY/dZ ~ 1.4 assuming that such a relation exists. During the analysis we firmly identify blue straggler stars, including the star 2-17, and find indications for the presence of their evolved counterparts. Our analysis supports the RGB mass-loss found from asteroseismology and we determine precisely the absolute mass of stars on the lower RGB, 1.15\pm0.02Msun. This will be an important consistency check for the detailed asteroseismology of cluster stars.

We present a detailed analysis of the upper main sequence of the 1.3 Gyr old open cluster Trumpler 20. High accuracy BV photometry combined with the Very Large Telescope/FLAMES medium-resolution spectroscopy of 954 stars is essential to understanding the unusual appearance of the color-magnitude diagram (CMD), initially suggesting multiple populations in Trumpler 20. We show that differential reddening is a dominant contributor to the apparent splitting/widening of the main-sequence turnoff region. At its extreme, the excess differential reddening reaches Delta(B-V)=0.1 while the adopted minimum reddening for the cluster is E(B-V)=0.36. A unique sample of measured projected rotational velocities indicates that stellar rotation is high near the main-sequence turnoff, reaching vsin i=180 km/s. By dividing the upper main-sequence stars into equal groups of slow and fast rotators, we find that fast rotators have a marginal blueshift of delta(V-I)=-0.01, corresponding to a difference in the median vsin i of 60 km/s between these subsamples. We conclude that stellar rotation has an insignificant effect on the morphology of the upper main sequence of this intermediate-age open cluster. Trumpler 20 appears to contain a single coeval population of stars but there is evidence that the red clump is extended.

We present a new distance determination to the Galactic globular cluster 47 Tucanae by fitting the spectral energy distributions of its white dwarfs to pure hydrogen atmosphere white dwarf models. Our photometric dataset is obtained from a 121 orbit Hubble Space Telescope program using the Wide Field Camera 3 UVIS/IR channels, capturing F390W, F606W, F110W, and F160W images. These images cover more than 60 square arcmins and extend over a radial range of 5-13.7 arcmin (6.5-17.9 pc) within the globular cluster. Using a likelihood analysis, we obtain a best fitting unreddened distance modulus of (m - M)o=13.36+/-0.02+/-0.06 corresponding to a distance of 4.70+/-0.04+/-0.13 kpc, where the first error is random and the second is systematic. We also search the white dwarf photometry for infrared excess in the F160W filter, indicative of debris disks or low mass companions, and find no convincing cases within our sample.

In HST Cycle 17, we imaged the well known globular star cluster 47 Tucanae for 121 orbits using ACS and both the UVIS and IR channels of WFC3 (GO-11677, PI - H. Richer). This unique data set was obtained to address many scientific questions that demand a very deep, panchromatic, and panoramic view of the cluster's stellar populations. In total, the program obtained over 0.75 Ms of imaging exposure time with the three HST cameras, over a time span of 9 months in 2010. The primary ACS field was imaged in the two broadband filters F606W and F814W filters, at 13 orientations, for all 121 orbits. The parallel WFC3 imaging provides a panchromatic (0.4 - 1.7 micron) and contiguous imaging swath over a 250 degree azimuthal range at impact radii of 6.5 -- 17.9 pc in 47 Tuc. This imaging totals over 60 arcmin^2 in area and utilizes the F390W and F606W broadband filters on WFC3/UVIS and the F110W and F160W broadband filters on WFC3/IR. In this paper, we describe the observational design of the new survey and one of the methods used to analyze all of the imaging data. This analysis combines over 700 full-frame images taken with the three HST cameras into a handful of ultra-deep, well-sampled combined images in each of the six filters. The results reveal unprecedented color-magnitude diagrams (CMDs) of the cluster extending to >30th magnitude in the optical, 29th magnitude in the UV, and 27th magnitude in the IR. The data set provides a characterization of the complete stellar populations of 47 Tuc, extending from the faintest hydrogen burning dwarfs through the main-sequence and giant branches, down to very cool white dwarf remnants in the cluster. The imaging also provides the deepest probe of the stellar populations of the background Small Magellanic Cloud (SMC) galaxy, resolving low mass main-sequence dwarfs with M < 0.2 Msun.

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.

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.

Detailed chemical abundances for twenty one elements are presented for four red giants in the anomalous outer halo globular cluster Palomar 1 ($R_{\rm{GC}} = 17.2$ kpc, $Z=3.6$ kpc) using high-resolution (R=36000) spectra from the High Dispersion Spectrograph (HDS) on the Subaru Telescope. Pal 1 has long been considered unusual because of its low surface brightness, sparse red giant branch, young age, and its possible association with two extragalactic streams of stars---this paper shows that its chemistry further confirms its unusual nature. The mean metallicity of the four stars, $[\rm{Fe/H}] = -0.60 \pm 0.01$, is high for a globular cluster so far from the Galactic center, but is low for a typical open cluster. The [$\alpha$/Fe] ratios, though in agreement with the Galactic stars within the $1\sigma$ errors, agree best with the lower values in dwarf galaxies. No signs of the Na/O anticorrelation are detected in Pal 1, though Na appears to be marginally high in all four stars. Pal 1's neutron capture elements are also unusual: its high [Ba/Y] ratio agrees best with dwarf galaxies, implying an excess of second-peak over first-peak s-process elements, while its [Eu/$\alpha$] and [Ba/Eu] ratios show that Pal 1's contributions from the r-process must have differed in some way from normal Galactic stars. Therefore, Pal 1 is chemically unusual, as well in its other properties. Pal 1 shares some of its unusual abundance characteristics with the young clusters associated with the Sagittarius dwarf galaxy remnant and the intermediate-age LMC clusters, and could be chemically associated with the Canis Majoris overdensity; however it does not seem to be similar to the Monoceros/Galactic Anticenter Stellar Stream.

Color-magnitude diagrams (CMDs) derived from Hubble Space Telescope (HST) Advanced Camera for Surveys F606W,F814W photometry of 6 globular clusters (GCs) are presented. The six GCs form two loose groupings in Galactocentric distance (Rgc): IC 4499, NGC 6426, and Ruprecht 106 at ~15-20 kpc and NGC 7006, Palomar 15, and Pyxis at ~40 kpc. The CMDs allow the ages to be estimated from the main sequence turnoff in every case. In addition, the age of Palomar 5 (Rgc ~ 18 kpc) is estimated using archival HST Wide Field Planetary Camera 2 V,I photometry. The age analysis reveals the following: IC 4499, Ruprecht 106, and Pyxis are 1-2 Gyr younger than inner halo GCs with similar metallicities; NGC 7006 and Palomar 5 are marginally younger than their inner halo counterparts; NGC 6426 and Palomar 15, the two most metal-poor GCs in the sample, are coeval with all the other metal-poor GCs within the uncertainties. Combined with our previous efforts, the current sample provides strong evidence that the Galactic GC age-metallicity relation consists of two distinct branches. One suggests a rapid chemical enrichment in the inner Galaxy while the other suggests prolonged GC formation in the outer halo. The latter is consistent with the outer halo GCs forming in dwarf galaxies and later being accreted by the Milky Way.

The Pan Andromeda Archeological Survey (PAndAS) CFHT Megaprime survey of the M31-M33 system has found a star stream which extends about 120 kpc NW from the center of M31. The great length of the stream, and the likelihood that it does not significantly intersect the disk of M31, means that it is unusually well suited for a measurement of stream gaps and clumps along its length as a test for the predicted thousands of dark matter sub-halos. The main result of this paper is that the density of the stream varies between zero and about three times the mean along its length on scales of 2 to 20 kpc. The probability that the variations are random fluctuations in the star density is less than 10^-5. As a control sample we search for density variations at precisely the same location in stars with metallicity higher than the stream, [Fe/H]=[0, -0.5] and find no variations above the expected shot noise. The lumpiness of the stream is not compatible with a low mass star stream in a smooth galactic potential, nor is it readily compatible with the disturbance caused by the visible M31 satellite galaxies. The stream's density variations appear to be consistent with the effects of a large population of steep mass function dark matter sub-halos, such as found in LCDM simulations, acting on an approximately 10Gyr old star stream. The effects of a single set of halo substructure realizations are shown for illustration, reserving a statistical comparison for another study.

We present the discovery of five new dwarf galaxies, Andromeda XXIII-XXVII, located in the outer halo of M31. These galaxies were found in the second year of data from the Pan-Andromeda Archaeological Survey (PAndAS) of the M31/M33 subgroup. This survey now provides an almost complete panoramic view of the M31 halo out to an average projected radius of ~150 kpc. The metal-poor stellar density map for this whole region serves, not only as an illustration of the discovery space for satellite galaxies, but also gives a birds-eye view of the ongoing assembly process of an L* disk galaxy. Four of the new dwarfs appear as well-defined spatial over-densities of stars lying on the expected locus of metal-poor red giant branch stars at the distance of M31. The fifth over-density, And XXVII, is embedded in an extensive stream of such stars and is possibly the remnant of a strong tidal disruption event. All five satellites have metallicities and luminosities typical of dwarf spheroidal galaxies and continue the trend whereby the brighter dwarf spheroidal satellites of M31 generally have much larger half-light radii than their Milky Way counterparts. With an extended sample of M31 satellite galaxies we also revisit the spatial distribution of this population and find that, within the current limits of the PAndAS survey, the surface density of satellites is essentially constant out to 150 kpc. This corresponds to a radial density distribution of satellites varying as 1/r, a result seemingly in conflict with the predictions of cosmological simulations. (Abridged)

We present a brief overview of a splinter session on determining the metallicity of low-mass dwarfs that was organized as part of the Cool Stars 16 conference. We review contemporary spectroscopic and photometric techniques for estimating metallicity in low-mass dwarfs and discuss the importance of measuring accurate metallicities for studies of Galactic and chemical evolution using subdwarfs, creating metallicity benchmarks for brown dwarfs, and searching for extrasolar planets that are orbiting around low--mass dwarfs. In addition, we present the current understanding of the effects of metallicity on stellar evolution and atmosphere models and discuss some of the limitations that are important to consider when comparing theoretical models to data.

We have surveyed $\sim40$sq.degrees surrounding M33 with CFHT MegaCam in the g and i filters, as part of the Pan-Andromeda Archaeological Survey. Our observations are deep enough to resolve the top 4mags of the red giant branch population in this galaxy. We have previously shown that the disk of M33 is surrounded by a large, irregular, low-surface brightness substructure. Here, we quantify the stellar populations and structure of this feature using the PAndAS data. We show that the stellar populations of this feature are consistent with an old population with $<[Fe/H]>\sim-1.6$dex and an interquartile range in metallicity of $\sim0.5$dex. We construct a surface brightness map of M33 that traces this feature to $\mu_V\simeq33$mags\,arcsec$^{-2}$. At these low surface brightness levels, the structure extends to projected radii of $\sim40$kpc from the center of M33 in both the north-west and south-east quadrants of the galaxy. Overall, the structure has an "S-shaped" appearance that broadly aligns with the orientation of the HI disk warp. We calculate a lower limit to the integrated luminosity of the structure of $-12.7\pm0.5$mags, comparable to a bright dwarf galaxy such as Fornax or AndII and slightly less than $1\$ of the total luminosity of M33. Further, we show that there is tentative evidence for a distortion in the distribution of young stars near the edge of the HI disk that occurs at similar azimuth to the warp in HI. The data also hint at a low-level, extended stellar component at larger radius that may be a M33 halo component. We revisit studies of M33 and its stellar populations in light of these new results, and we discuss possible formation scenarios for the vast stellar structure. Our favored model is that of the tidal disruption of M33 in its orbit around M31.

Stellar physics and evolution calculations enable a broad range of research in astrophysics. Modules for Experiments in Stellar Astrophysics (MESA) is a suite of open source libraries for a wide range of applications in computational stellar astrophysics. A newly designed 1-D stellar evolution module, MESA star, combines many of the numerical and physics modules for simulations of a wide range of stellar evolution scenarios ranging from very-low mass to massive stars, including advanced evolutionary phases. MESA star solves the fully coupled structure and composition equations simultaneously. It uses adaptive mesh refinement and sophisticated timestep controls, and supports shared memory parallelism based on OpenMP. Independently usable modules provide equation of state, opacity, nuclear reaction rates, and atmosphere boundary conditions. Each module is constructed as a separate Fortran 95 library with its own public interface. Examples include comparisons to other codes and show evolutionary tracks of very low mass stars, brown dwarfs, and gas giant planets; the complete evolution of a 1 Msun star from the pre-main sequence to a cooling white dwarf; the Solar sound speed profile; the evolution of intermediate mass stars through the thermal pulses on the He-shell burning AGB phase; the interior structure of slowly pulsating B Stars and Beta Cepheids; evolutionary tracks of massive stars from the pre-main sequence to the onset of core collapse; stars undergoing Roche lobe overflow; and accretion onto a neutron star. Instructions for downloading and installing MESA can be found on the project web site (http://mesa.sourceforge.net/).

We present new measurements of the centers for 65 Milky Way globular clusters. Centers were determined by fitting ellipses to the density distribution within the inner $2\arcmin$ of the cluster center, and averaging the centers of these ellipses. The symmetry of clusters was also analyzed by comparing cumulative radial distributions on opposite sides of the cluster across a grid of trial centers. All of the determinations were done with stellar positions derived from a combination of two single-orbit ACS images of the core of the cluster in $F606W$ and $F814W$. We find that the ellipse-fitting method provides remarkable accuracy over a wide range of core sizes and density distributions, while the symmetry method is difficult to use on clusters with very large cores, or low density. The symmetry method requires a larger field, or a very sharply peaked density distribution.

We investigate the effect of three important processes by which AGN-blown bubbles transport material: drift, wake transport and entrainment. The first of these, drift, occurs because a buoyant bubble pushes aside the adjacent material, giving rise to a net upward displacement of the fluid behind the bubble. For a spherical bubble, the mass of upwardly displaced material is roughly equal to half the mass displaced by the bubble, and should be ~ 10^{7-9} solar masses depending on the local ICM and bubble parameters. We show that in classical cool core clusters, the upward displacement by drift may be a key process in explaining the presence of filaments behind bubbles. A bubble also carries a parcel of material in a region at its rear, known as the wake. The mass of the wake is comparable to the drift mass and increases the average density of the bubble, trapping it closer to the cluster centre and reducing the amount of heating it can do during its ascent. Moreover, material dropping out of the wake will also contribute to the trailing filaments. Mass transport by the bubble wake can effectively prevent the build-up of cool material in the central galaxy, even if AGN heating does not balance ICM cooling. Finally, we consider entrainment, the process by which ambient material is incorporated into the bubble. Abridged

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)

We use photometric and spectroscopic observations of the eclipsing binary V69-47 Tuc to derive the masses, radii, and luminosities of the component stars. Based on measured systemic velocity, distance, and proper motion, the system is a member of the globular cluster 47 Tuc. The system has an orbital period of 29.5 d and the orbit is slightly eccentric with e=0.056. We obtain Mp=0.8762 +- 0.0048 M(Sun), Rp=1.3148 +-0.0051 R(Sun), Lp=1.94 +- 0.21 L(Sun) for the primary and Ms=0.8588 +- 0.0060 M(Sun), Rs=1.1616 +- 0.0062 R(Sun), Ls=1.53 +- 0.17 L(Sun) for the secondary. These components of V69 are the first Population II stars with masses and radii derived directly and with an accuracy of better than 1%. We measure an apparent distance modulus of (m-M)v=13.35 +- 0.08 to V69. We compare the absolute parameters of V69 with five sets of stellar evolution models and estimate the age of V69 using mass-luminosity-age, mass-radius-age, and turnoff mass - age relations. The masses, radii, and luminosities of the component stars are determined well enough that the measurement of ages is dominated by systematic differences between the evolutionary models, in particular, the adopted helium abundance. By comparing the observations to Dartmouth model isochrones we estimate the age of V69 to be 11.25 +- 0.21(random) +- 0.85(systematic) Gyr assuming [Fe/H]=-0.70, [alpha/Fe]=0.4, and Y=0.255. The determination of the distance to V69, and hence to 47Tuc, can be further improved when infrared eclipse photometry is obtained for the variable.

We present simple stellar population (SSP) models with scaled-solar and alpha-element enhanced abundances. The SSP models are based on the Dartmouth Stellar Evolution Database, our library of synthetic stellar spectra, and a detailed systematic variation of horizontal-branch (HB) morphology with age and metallicity. In order to test the relative importance of a variety of SSP model ingredients, we compare our SSP models with integrated spectra of 41 Milky Way Globular Clusters (MWGCs) from Schiavon et al. (2005). Using the Mg b and Ca4227 indices, we confirm that Mg and Ca are enhanced by about +0.4 and +0.2 dex, respectively, in agreement with results from high resolution spectra of individual stars in MWGCs. Balmer lines, particularly Hgamma and Hdelta, of MWGCs are reproduced by our alpha-enhanced SSP models not only because of the combination of isochrone and spectral effects but also because of our reasonable HB treatment. Moreover, it is shown that the Mg abundance significantly influences Balmer and iron line indices. Finally, the investigation of power-law initial mass function (IMF) variations suggests that an IMF much shallower than Salpeter is unrealistic because the Balmer lines are too strong on the metal-poor side to be compatible with observations.

In hierarchical cosmological models, galaxies grow in mass through the continual accretion of smaller ones. The tidal disruption of these systems is expected to result in loosely bound stars surrounding the galaxy, at distances that reach $10 - 100$ times the radius of the central disk. The number, luminosity and morphology of the relics of this process provide significant clues to galaxy formation history, but obtaining a comprehensive survey of these components is difficult because of their intrinsic faintness and vast extent. Here we report a panoramic survey of the Andromeda galaxy (M31). We detect stars and coherent structures that are almost certainly remnants of dwarf galaxies destroyed by the tidal field of M31. An improved census of their surviving counterparts implies that three-quarters of M31's satellites brighter than $M_V < -6$ await discovery. The brightest companion, Triangulum (M33), is surrounded by a stellar structure that provides persuasive evidence for a recent encounter with M31. This panorama of galaxy structure directly confirms the basic tenets of the hierarchical galaxy formation model and reveals the shared history of M31 and M33 in the unceasing build-up of galaxies.

Age constraints are most often placed on globular clusters by comparing their CMDs with theoretical isochrones. The recent discoveries of detached, eclipsing binaries in such systems by the Cluster AgeS Experiment (CASE) provide new insights into their ages and, at the same time, much needed tests of stellar evolution models. We describe efforts to model the properties of the detached, eclipsing binary V69 in 47 Tuc and compare age constraints derived from stellar evolution models of V69A and B with ages obtained from fitting isochrones to the cluster CMD. We determine, under reasonable assumptions of distance, reddening, and metallicity, the extent to which it is possible to simultaneously constrain the age and He content of 47 Tuc.

Context: Ogata et al. (2009; hereafter OKK) presented a theoretical determination of the triple-alpha nuclear reaction rate. Their rate differs from the NACRE rate by many orders of magnitude at temperatures relevant for low mass stars. Aims: We explore the evolutionary implications of adopting the OKK triple-alpha reaction rate in low mass stars and compare the results with those obtained using the NACRE rate. Methods: The triple-alpha reaction rates are compared by following the evolution of stellar models at 1 and 1.5 Msol with Z=0.0002 and Z=0.02. Results: Results show that the OKK rate has severe consequences for the late stages of stellar evolution in low mass stars. Most notable is the shortening--or disappearance--of the red giant phase. Conclusions: The OKK triple-alpha reaction rate is incompatible with observations of extended red giant branches and He burning stars in old stellar systems.

We present an analysis of Two Micron All Sky Survey (2MASS) calibration photometry of the old open cluster M67 (NGC 2682). The proper motion-cleaned color-magnitude diagram (CMD) resulting from these data extends ~3 magnitudes deeper than one based on data from the point source catalog. The CMD extends from above the helium-burning red clump to a faint limit that is more than 7 magnitudes below the main sequence turnoff in the Ks band. After adopting a reddening of E(B-V) = 0.041 +/- 0.004 and a metal abundance of [Fe/H] = -0.009 +/- 0.009 based on a survey of published values, we fit the unevolved main sequence of M67 to field main sequence stars with 2MASS photometry and Hipparcos parallaxes. This analysis yields distance moduli of (m-M)Ks = 9.72 +/- 0.05 and (m-M)o = 9.70 +/- 0.05, which are consistent with published values. We compare the theoretical isochrones of Girardi et al. and Dotter et al. to the CMD of M67 and comment on the relative merits of each set of models. These comparisons suggest an age between 3.5 and 4.0 Gyr for M67. The depth of the M67 data make them ideal for the calibration of a new age indicator that has recently been devised by Calamida et al.- the difference in (J-Ks) color between the main sequence turnoff (TO) and the point on the lower main sequence where it turns down (TD) and becomes nearly vertical [D(J-Ks)]. Coupled with deep 2MASS photometry for three other open clusters, NGC 2516, M44, and NGC 6791, we calibrate D(J-Ks) in terms of age and find D(J-Ks) = (3.017 +/- 0.347) - (0.259 +/- 0.037)*Log Age (yrs).

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.

The connection between mass loss on the red giant branch (RGB) and horizontal branch (HB) morphology in globular clusters (GCs) has long been acknowledged but the mechanisms governing mass loss remains poorly understood from a theoretical perspective. The present study uses synthetic HB models to demonstrate for the first time that alpha-enhancement and a simple relation between mass loss and metallicity can explain the entire range of HB morphology (characterized by the HB type index) observed in old, coeval GCs. The mass loss-metallicity relation accounts naturally for the fact that the most metal poor GCs ([Fe/H] < -2) have redder HBs than is typical of GCs with -2 < [Fe/H] < -1.5 without invoking younger ages. These results may prove useful in studying the contribution of HB stars to integrated light via stellar population synthesis.

AM-1, at ~120 kpc, and Pal 14, at ~70 kpc, are two of the most distant Galactic globular clusters known. We present Hubble Space Telescope WFPC2 photometry of AM-1 and Pal 14 that reveals unprecedented depth and detail in the color-magnitude diagrams of these two clusters. Absolute and relative age measurements confirm that both are younger than the inner halo globular cluster M 3 by 1.5-2 Gyr, assuming all three clusters have similar compositions. Thus AM-1 and Pal 14 join Pal 3, Pal 4, and Eridanus (studied by Stetson et al.) as distant Galactic globular clusters with red horizontal branches and young ages relative to the inner halo. Within the context of the entire body of research on the ages of second parameter globular clusters, the observed correlation between age and horizontal branch morphology suggests that age is the best candidate for the second parameter. However, this conclusion is tempered by the lack of precise chemical abundance determinations for a significant fraction of second parameter globular clusters.

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.

We present N-body models to complement deep imaging of the metal-poor core-collapsed cluster NGC6397 obtained with the Hubble Space Telescope. All simulations include stellar and binary evolution in-step with the stellar dynamics and account for the tidal field of the Galaxy. We focus on the results of a simulation that began with 100000 objects (stars and binaries), 5% primordial binaries and Population II metallicity. After 16 Gyr of evolution the model cluster has about 20% of the stars remaining and has reached core-collapse. We compare the color-magnitude diagrams of the model at this age for the central region and an outer region corresponding to the observed field of NGC6397 (about 2-3 half-light radii from the cluster centre). This demonstrates that the white dwarf population in the outer region has suffered little modification from dynamical processes - contamination of the luminosity function by binaries and white dwarfs with non-standard evolution histories is minimal and should not significantly affect measurement of the cluster age. We also show that the binary fraction of main-sequence stars observed in the NGC6397 field can be taken as representative of the primordial binary fraction of the cluster. For the mass function of the main-sequence stars we find that although this has been altered significantly by dynamics over the cluster lifetime, especially in the central and outer regions, that the position of the observed field is close to optimal for recovering the initial mass function of the cluster stars (below the current turn-off mass). More generally we look at how the mass function changes with radius in a dynamically evolved stellar cluster and suggest where the best radial position to observe the initial mass function is for clusters of any age.

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 the CMD from deep HST imaging in the globular cluster NGC 6397. The ACS was used for 126 orbits to image a single field in two colors (F814W, F606W) 5 arcmin SE of the cluster center. The field observed overlaps that of archival WFPC2 data from 1994 and 1997 which were used to proper motion (PM) clean the data. Applying the PM corrections produces a remarkably clean CMD which reveals a number of features never seen before in a globular cluster CMD. In our field, the main sequence stars appeared to terminate close to the location in the CMD of the hydrogen-burning limit predicted by two independent sets of stellar evolution models. The faintest observed main sequence stars are about a magnitude fainter than the least luminous metal-poor field halo stars known, suggesting that the lowest luminosity halo stars still await discovery. At the bright end the data extend beyond the main sequence turnoff to well up the giant branch. A populous white dwarf cooling sequence is also seen in the cluster CMD. The most dramatic features of the cooling sequence are its turn to the blue at faint magnitudes as well as an apparent truncation near F814W = 28. The cluster luminosity and mass functions were derived, stretching from the turn off down to the hydrogen-burning limit. It was well modeled with either a very flat power-law or a lognormal function. In order to interpret these fits more fully we compared them with similar functions in the cluster core and with a full N-body model of NGC 6397 finding satisfactory agreement between the model predictions and the data. This exercise demonstrates the important role and the effect that dynamics has played in altering the cluster IMF.

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.

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

We present the results of a deep Hubble Space Telescope (HST) exposure of the nearby globular cluster NGC6397, focussing attention on the cluster's white dwarf cooling sequence. This sequence is shown to extend over 5 magnitudes in depth, with an apparent cutoff at magnitude F814W=27.6. We demonstrate, using both artificial star tests and the detectability of background galaxies at fainter magnitudes, that the cutoff is real and represents the truncation of the white dwarf luminosity function in this cluster. We perform a detailed comparison between cooling models and the observed distribution of white dwarfs in colour and magnitude, taking into account uncertainties in distance, extinction, white dwarf mass, progenitor lifetimes, binarity and cooling model uncertainties. After marginalising over these variables, we obtain values for the cluster distance modulus and age of \mu_0 = 12.02 \pm 0.06 and T_c = 11.47 \pm 0.47Gyr (95% confidence limits). Our inferred distance and white dwarf initial-final mass relations are in good agreement with other independent determinations, and the cluster age is consistent with, but more precise than, prior determinations made using the main sequence turnoff method. In particular, within the context of the currently accepted \Lambda CDM cosmological model, this age places the formation of NGC6397 at a redshift z=3, at a time when the cosmological star formation rate was approaching its peak.

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.

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)

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.

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.