We use recently derived ages for 61 Milky Way (MW) globular clusters (GCs) to show that their age-metallicity relation (AMR) can be divided into two distinct, parallel sequences at [Fe/H] $ga -1.8$. Approximately one-third of the clusters form an off
set sequence that spans the full range in age ($sim 10.5$--13 Gyr), but is more metal rich at a given age by $sim 0.6$ dex in [Fe/H]. All but one of the clusters in the offset sequence show orbital properties that are consistent with membership in the MW disk. They are not simply the most metal-rich GCs, which have long been known to have disk-like kinematics, but they are the most metal-rich clusters at all ages. The slope of the mass-metallicity relation (MMR) for galaxies implies that the offset in metallicity of the two branches of the AMR corresponds to a mass decrement of 2 dex, suggesting host galaxy masses of $M_{*} sim 10^{7-8} msol$ for GCs that belong to the more metal-poor AMR. We suggest that the metal-rich branch of the AMR consists of clusters that formed in-situ in the disk, while the metal-poor GCs were formed in relatively low-mass (dwarf) galaxies and later accreted by the MW. The observed AMR of MW disk stars, and of the LMC, SMC and WLM dwarf galaxies are shown to be consistent with this interpretation, and the relative distribution of implied progenitor masses for the halo GC clusters is in excellent agreement with the MW subhalo mass function predicted by simulations. A notable implication of the bifurcated AMR, is that the identical mean ages and spread in ages, for the metal rich and metal poor GCs are difficult to reconcile with an in-situ formation for the latter population.
Ages have been derived for 55 globular clusters (GCs) from overlays of isochrones onto the turnoff photometry, assuming distances based on fits of zero-age horizontal branch (ZAHB) models to the lower bound of the observed distributions of HB stars.
The error bar arising just from the fitting of ZAHBs and isochrones is ~ +/- 0.25 Gyr, while that associated with distance and chemical abundance uncertainties is ~ +/- 1.5-2 Gyr. Ages vary from mean values of ~12.5 Gyr at [Fe/H] < -1.7 to ~11 Gyr at [Fe/H] > -1.0. At intermediate metallicities, the age-metallicity relation (AMR) appears to be bifurcated: one branch apparently contains clusters with disk-like kinematics, whereas the other branch is populated by clusters with halo-type orbits. There is no apparent dependence of age on Galactocentric distance (R_G) nor is there a clear correlation of HB type with age. Subtle variations in the subgiant branch (SGB) slopes of [Fe/H] < -1.5 GCs are tentatively attributed to helium abundance differences. Curiously, GCs with steep M13-like SGBs tend to be massive systems, located at small R_G, that show the strongest evidence for multiple stellar populations. The others are typically low-mass systems that, at the present time, should not be able to retain the matter lost by mass-losing stars. The apparent separation of the two groups in terms of their present-day gas retention properties is difficult to understand if all GCs were initial ~20 times their current masses. The lowest mass systems may have never been able to retain enough gas to produce a significant population of second-generation stars; in this case, the observed light element abundance variations were presumably present in the gas out of which the observed cluster stars formed.
[Abbreviated] We have investigated the color-magnitude diagram of Omega Centauri and find that the blue main sequence (bMS) can be reproduced only by models that have a of helium abundance in the range Y=0.35-$0.40. To explain the faint subgiant bran
ch of the reddest stars (MS-a/RG-a sequence), isochrones for the observed metallicity ([Fe/H]approx0.7) appear to require both a high age (~13Gyr) and enhanced CNO abundances ([CNO/Fe]approx0.9$). Y~0.35 must also be assumed in order to counteract the effects of high CNO on turnoff colors, and thereby to obtain a good fit to the relatively blue turnoff of this stellar population. This suggest a short chemical evolution period of time (<1Gyr) for Omega Cen. Our intermediate-mass (super-)AGB models are able to reproduce the high helium abundances, along with [N/Fe]~2 and substantial O depletions if uncertainties in the treatment of convection are fully taken into account. These abundance features distinguish the bMS stars from the dominant [Fe/H] $approx1.7$ population. The most massive super-AGB stellar models (M_zams>=6.8M_sun, M_He,core>=1.245M_sun) predict too large N-enhancements, which limits their role in contributing to the extreme populations. We show quantitatively that highly He- and N-enriched AGB ejecta have particularly efficient cooling properties. Based on these results and on the reconstruction of the orbit of Omega Cen with respect to the Milky Way we propose the galactic plane passage gas purging scenario for the chemical evolution of this cluster. Our model addresses the formation and properties of the bMS population (including their central location in the cluster). We follow our model descriptively through four passage events, which could explain not only some key properties of the bMS, but also of the MS-a/RGB-a and the s-enriched stars.
We describe an extensive observational project that has obtained high-quality and homogeneous photometry for a number of different Galactic star clusters (including M 92, M 13, M 3, M 71, and NGC 6791) spanning a wide range in metallicity (-2.3<[Fe/H
]<+0.4), as observed in the ugriz passbands with the MegaCam wide-field imager on the Canada-France-Hawaii Telescope. By employing these purest of stellar populations, fiducial sequences have been defined from color-magnitude diagrams that extend from the tip of the red-giant branch down to approximately 4 magnitudes below the turnoff: these sequences have been accurately calibrated to the standard ugriz system via a set of secondary photometric standards located within these same clusters. Consequently, they can serve as a valuable set of empirical fiducials for the interpretation of stellar populations data in the ugriz system.
Stellar models have been computed for stars having [Fe/H] = 0.0 and -2.0 to determine the effects of using boundary conditions derived from the latest MARCS model atmospheres. The latter were fitted to the interior models at both the photosphere and
at tau = 100, and at least for the 0.8-1.0 solar mass stars considered here, the resultant evolutionary tracks were found to be nearly independent of the chosen fitting point. Particular care was taken to treat the entire star as consistently as possible; i.e., both the interior and atmosphere codes assumed the same abundances and the same treatment of convection. Tracks were also computed using either the classical gray T(tau,T_eff) relation or that derived by Krishna Swamy (1966) to derive the boundary pressure. The latter predict warmer giant branches (by ~150 K) at solar abundances than those based on gray or MARCS atmospheres, which happens to be in good agreement with the inferred temperatures of giants in the open cluster M67 from the latest (V-K)-T_eff relations. Most of the calculations assumed Z=0.0125 (Asplund et al.), though a few models were computed for Z=0.0165 (Grevesse & Sauval) to determine the dependence of the tracks on Z_odot. Grids of scaled solar, differentially corrected (SDC) atmospheres were also computed to try to improve upon theoretical MARCS models. When they were used as boundary conditions, the resultant tracks agreed very well with those based on a standard scaled-solar (e.g., Krishna Swamy) T(tau,T_eff) relation, independently of the assumed metal abundance. Fits of isochrones to the C-M diagram of the [Fe/H] = -2 globular cluster M68 were examined, as was the possibility that the mixing-length parameter varies with stellar parameters.
The mass at which a transition is made between stars that have radiative or convective cores throughout the core H-burning phase is a fairly sensitive function of Z (particularly the CNO abundances). As a consequence, the ~4 Gyr, open cluster M67 pro
vides a constraint on Z_odot (and the solar heavy-element mixture) because (i) high-resolution spectroscopy indicates that this system has virtually the same metal abundances as the Sun, and (ii) its turnoff stars have masses just above the lower limit for sustained core convection on the main sequence. In this study, evolutionary tracks and isochrones using the latest MARCS model atmospheres as boundary conditions have been computed for 0.6-1.4 solar masses on the assumption of a metals mix (implying Z_odot = 0.0125) based on the solar abundances derived by M. Asplund and collaborators using 3-D model atmospheres. These calculations do not predict a turnoff gap where one is observed in M67. No such difficulty is found if the analysis uses isochrones for Z_odot = 0.0165, assuming the Grevesse & Sauval (1998) mix of heavy elements. Our findings, like the inferences from helioseismology, indicate a problem with the Asplund et al. abundances. However, it is possible that low-Z models with diffusive processes taken into account will be less problematic.
