Source author record

Aaron Dotter

Aaron Dotter appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

38works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

38 published item(s)

preprint2021arXiv

Predictions for gravity-mode periods and surface abundances in intermediate-mass dwarfs from shear mixing and radiative levitation

The treatment of chemical mixing in the radiative envelopes of intermediate-mass stars has hardly been calibrated so far. Recent asteroseismic studies demonstrated that a constant diffusion coefficient in the radiative envelope is not able to explain the periods of trapped gravity modes in the oscillation spectra of $γ$ Doradus pulsators. We present a new generation of MESA stellar models with two major improvements. First, we present a new implementation for computing radiative accelerations and Rosseland mean opacities that requires significantly less CPU time. Second, the inclusion of shear mixing based on rotation profiles computed with the 2D stellar structure code ESTER is considered. We show predictions for the mode periods of these models covering stellar masses from 1.4 to 3.0${\rm M_\odot}$ across the main sequence (MS), computed for different metallicities. The morphology of the chemical mixing profile resulting from shear mixing in combination with atomic diffusion and radiative levitation does allow for mode trapping, while the diffusion coefficient in the outer envelope is large ($>10^{6}\,{\rm cm^2\,s^{-1}}$). Furthermore, we make predictions for the evolution of surface abundances for which radiative accelerations can be computed. We find that the N/C and C/O abundance ratios correlate with stellar age. We predict that these correlations are observable with precisions $\lesssim 0.1$ dex on these ratios, given that a precise age estimate can be made.

preprint2021arXiv

Probing the progenitors of spinning binary black-hole mergers with long gamma-ray bursts

Long-duration gamma-ray bursts are thought to be associated with the core-collapse of massive, rapidly spinning stars and the formation of black holes. However, efficient angular momentum transport in stellar interiors, currently supported by asteroseismic and gravitational-wave constraints, leads to predominantly slowly-spinning stellar cores. Here, we report on binary stellar evolution and population synthesis calculations, showing that tidal interactions in close binaries not only can explain the observed sub-population of spinning, merging binary black holes but also lead to long gamma-ray bursts at the time of black-hole formation. Given our model calibration against the distribution of isotropic-equivalent energies of luminous long gamma-ray bursts, we find that ~10% of the GWTC-2 reported binary black holes had a luminous long gamma-ray burst associated with their formation, with GW190517 and GW190719 having a probability of ~85% and ~60%, respectively, being among them. Moreover, given an assumption about their average beaming fraction, our model predicts the rate density of long gamma-ray bursts, as a function of redshift, originating from this channel. For a constant beaming fraction $f_\mathrm{B}\sim 0.05$ our model predicts a rate density comparable to the observed one, throughout the redshift range, while, at redshift $z \in [0,2.5]$, a tentative comparison with the metallicity distribution of observed LGRB host galaxies implies that between 20% to 85% of the observed long gamma-ray bursts may originate from progenitors of merging binary black holes. The proposed link between a potentially significant fraction of observed, luminous long gamma-ray bursts and the progenitors of spinning binary black-hole mergers allows us to probe the latter well outside the horizon of current-generation gravitational wave observatories, and out to cosmological distances.

preprint2020arXiv

Asteroseismic modeling of gravity modes in slowly rotating A/F stars with radiative levitation

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.

preprint2020arXiv

MINESweeper: Spectrophotometric Modeling of Stars in the Gaia Era

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.

preprint2019arXiv

Combined Effects of Rotation and Age Spreads on Extended Main Sequence Turn Offs

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.

preprint2016arXiv

GRACES observations of young [alpha/Fe]-rich stars

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.

preprint2016arXiv

MESA Isochrones and Stellar Tracks (MIST) 0: Methods for the construction of stellar isochrones

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.

preprint2016arXiv

MESA Isochrones and Stellar Tracks (MIST). I: Solar-Scaled Models

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

preprint2015arXiv

Capella (alpha Aurigae) revisited: New binary orbit, physical properties, and evolutionary state

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.

preprint2015arXiv

Integrated Light Chemical Tagging Analyses of Seven M31 Outer Halo Globular Clusters from the Pan-Andromeda Archaeological Survey

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.

preprint2015arXiv

Rapidly rotating second-generation progenitors for the blue hook stars of ω Cen

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

preprint2014arXiv

Optimal Integrated Abundances for Chemical Tagging of Extragalactic Globular Clusters

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.

preprint2014arXiv

Stellar Models of Multiple Populations in Globular Clusters. I. The Main Sequence of NGC 6752

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.

preprint2013arXiv

A Dynamical Signature of Multiple Stellar Populations in 47 Tucanae

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.

preprint2013arXiv

A peculiar faint satellite in the remote outer halo of M31

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.

preprint2013arXiv

Comparing the White Dwarf Cooling Sequences in 47 Tuc and NGC 6397

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.

preprint2013arXiv

Modules for Experiments in Stellar Astrophysics (MESA): Giant Planets, Oscillations, Rotation, and Massive Stars

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]

preprint2013arXiv

Spectrum Syntheses of High Resolution Integrated Light Spectra of Galactic Globular Clusters

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.

preprint2013arXiv

The Interior Structure Constants as an Age Diagnostic for Low-Mass, Pre-Main Sequence Detached Eclipsing Binary Stars

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.

preprint2013arXiv

The Second Parameter Problem(s)

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.

preprint2013arXiv

WOCS 40007: A Detached Eclipsing Binary near the Turnoff of the Open Cluster NGC 6819

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.

preprint2012arXiv

Stellar Models with Enhanced Abundances of C, N, O, Ne, Na, Mg, Si, S, Ca, and Ti, in Turn, at Constant Helium and Iron Abundances

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.

preprint2011arXiv

A Deep, Wide-Field, and Panchromatic View of 47 Tuc and the SMC with HST: Observations and Data Analysis Methods

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.

preprint2011arXiv

Accurate Low-Mass Stellar Models of KOI-126

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.

preprint2011arXiv

Detailed Chemical Abundances of Four Stars in the Unusual Globular Cluster, Palomar 1

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 [$α$/Fe] ratios, though in agreement with the Galactic stars within the $1σ$ 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/$α$] 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.

preprint2011arXiv

Determining the Metallicity of Low-Mass Stars and Brown Dwarfs: Tools for Probing Fundamental Stellar Astrophysics, Tracing Chemical Evolution of the Milky Way and Identifying the Hosts of Extrasolar Planets

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.

preprint2011arXiv

Globular clusters in the outer Galactic halo: new HST/ACS imaging of 6 globular clusters and the Galactic globular cluster age-metallicity relation

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.

preprint2011arXiv

PAndAS' progeny: extending the M31 dwarf galaxy cabal

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)

preprint2011arXiv

The ACS Survey of Galactic Globular Clusters XI: The Three-Dimensional Orientation of the Sagittarius Dwarf Spheroidal Galaxy and its Globular Clusters

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.

preprint2010arXiv

Mass transport by buoyant bubbles in galaxy clusters

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

preprint2010arXiv

Modules for Experiments in Stellar Astrophysics (MESA)

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

preprint2010arXiv

The ACS Survey of Galactic Globular Clusters. X. New Determinations of Centers for 65 Clusters

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.

preprint2010arXiv

The photometric properties of a vast stellar substructure in the outskirts of M33

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 $μ_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.

preprint2009arXiv

Comparison of Alpha-Element Enhanced Simple Stellar Population Models with Milky Way Globular Clusters

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.

preprint2009arXiv

Deep 2MASS Photometry of M67 and Calibration of the Main Sequence J-Ks Color Difference as an Age Indicator

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

preprint2009arXiv

Evolutionary implications of the new triple-alpha nuclear reaction rate for low mass stars

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.

preprint2009arXiv

The ACS Survey of Galactic Globular Clusters. IX. Horizontal Branch Morphology and the Second Parameter Phenomenon

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)

preprint2008arXiv

The ACS Survey of Galactic Globular Clusters. VII. Relative Ages

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.