PUBLISHED IN APJ, 797, 35 (2014) Preprint typeset using LATEX style emulateapj v. 08/22/09 EXTENDED MAIN SEQUENCE TURNOFFS IN INTERMEDIATE-AGE STAR CLUSTERS: A CORRELATION BETWEEN TURNOFF WIDTH AND EARLY ESCAPE VELOCITY1 PAUL GOUDFROOIJ2,LÉO GIRARDI3,VERA KOZHURINA-PLATAIS2, JASON S. KALIRAI2,4,IMANTS PLATAIS5,THOMAS H. PUZIA6, MATTEO CORRENTI2,ALESSANDRO BRESSAN7,RUPALI CHANDAR8,LEANDRO KERBER9,PAOLA MARIGO10, AND STEFANO RUBELE3 2 Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA; goudfroo, verap, jkalirai, [email protected] 3 Osservatorio Astronomico di Padova – INAF, Vicolo dell’Osservatorio 5, I-35122 Padova, Italy; [email protected] 4 Center for Astrophysical Sciences, Johns Hopkins University, 3400 N. Charles Street, Baltimore, MD 21218, USA 5 Department of Physics and Astronomy, Johns Hopkins University, 3400 N. Charles Street, Baltimore, MD 21218, USA; [email protected] 6 Institute of Astrophysics, Pontificia Universidad Católica de Chile, Avenida Vicuña Mackenna 4860, Macul, Santiago, Chile; [email protected] 7 SISSA, via Bonomea 365, I-34136 Trieste, Italy; [email protected] 8 Department of Physics and Astronomy, The University of Toledo, 2801 West Bancroft Street, Toledo, OH 43606, USA; [email protected] 9 Universidade Estadual de Santa Cruz, Rodovia Ilhéus-Itabuna, km 16, 45662-000 Ilhéus, Bahia, Brazil; [email protected] 10 Dipartimento di Fisica e Astronomia Galileo Galilei, Università di Padova, Vicolo dell’Osservatorio 3, I-35122 Padova, Italy; [email protected] Published in ApJ, 797, 35 (2014) ABSTRACT We present color-magnitude diagram analysis of deep Hubble Space Telescope imaging of a mass-limited sample of 18 intermediate-age (1–2 Gyr old) star clusters in the Magellanic Clouds, including 8 clusters for which new data was obtained. We find that all star clusters in our sample feature extended main sequence turnoff (eMSTO) regions that are wider than can be accounted for by a simple stellar population (including unresolved binary stars). FWHM widths of the MSTOs indicate age spreads of 200–550 Myr. We evaluate dynamical evolution of clusters with and without initial mass segregation. Our main results are: (1) the fraction of red clump (RC) stars in secondary RCs in eMSTO clusters scales with the fraction of MSTO stars having pseudo-ages . 1.35 Gyr; (2) the width of the pseudo-age distributions of eMSTO clusters is correlated with their central escape velocity vesc, both currently and at an age of 10 Myr. We find that these two results are unlikely to be reproduced by the effects of interactive binary stars or a range of stellar rotation velocities. We therefore argue that the eMSTO phenomenon is mainly caused by extended star formation within the clusters; −1 (3) we find that vesc 15 kms out to ages of at least 100 Myr for all clusters featuring eMSTOs, while −1 ≥ vesc 12 kms at all ages for two lower-mass clusters in the same age range that do not show eMSTOs. We argue≤ that eMSTOs only occur for clusters whose early escape velocities are higher than the wind velocities of stars that provide material from which second-generationstars can form. The threshold of 12–15 kms−1 is consistent with wind velocities of intermediate-mass AGB stars and massive binary stars in the literature. Subject headings: globular clusters: general — Magellanic Clouds 1. INTRODUCTION ing along the RGB. Instead, their origin must be primordial, For almost a century and counting, the study of globular being imprinted on the stars during their formation process. clusters (GCs) has contributed enormously to our understand- The chemical processes involved in causing the light-element ing of stellar evolution. Until recently, this was especially true abundance variations have largely been identified as proton capture reactions at T & 2 107 K, such as the CNO and because they were thought to be simple objects consisting of × thousands to millions of coeval stars with the same chemical NeNa cycles. Currently, the leading candidates for “polluter” composition. However, this notion has had to face serious sources are stars in which such reactions occur readily and challenges over the last dozen years. It is now commonly which feature slow stellar winds so that their ejecta are rela- arXiv:1410.3840v3 [astro-ph.SR] 24 Nov 2014 recognized that GCs typically∼ harbor multiple stellar popula- tively easy to retain within the potential well of massive clus- ters: (i) intermediate-mass AGB stars (4 . /M⊙ . 8, here- tions featuring several unexpected characteristics (for recent M reviews, see Renzini 2008; Gratton et al. 2012). after IM-AGB; e.g., D’Antona & Ventura 2007 and references Recent spectroscopic surveys established that light ele- therein), (ii) rapidly rotating massive stars (often referred to as ments like C, N, O, Na, Mg, and Al show large star-to-star “FRMS”; Decressin et al. 2007) and (iii) massive binary stars abundance variations (often dubbed “Na-O anticorrelations”) (de Mink et al. 2009). within virtually all Galactic GCs studied to date in sufficient In the two currently favored formation scenarios, the abun- detail (Carretta et al. 2010, and references therein). These dance variations are due to stars having either formed from abundance variations have been found among both red giant or polluted by gas that is a mixture of pristine material branch (RGB) stars and main sequence (MS) stars in sev- and material shed by such “polluters”. In the “in situ star eral GCs (Gratton et al. 2004). This clarified that the vari- formation” scenario (see, e.g., D’Ercole et al. 2008, 2010; ations cannot be due to internal mixing within stars evolv- Conroy & Spergel 2011), the abundance variations are due to a second generation of stars that formed out of gas clouds 1 Based on observations with the NASA/ESA Hubble Space Telescope, that were polluted by winds of first-generation stars to vary- obtained at the Space Telescope Science Institute, which is operated by the ing extents, during a period spanning up to a few hundreds of Association of Universities for Research in Astronomy, Inc., under NASA Myr, depending on the nature of the polluters. In the alter- contract NAS5-26555 native “early disc accretion” scenario (Bastian et al. 2013a), 2 Goudfrooij et al. the polluted gas is instead accreted by low-mass pre-main- further study potential effects of dynamical properties on the sequence stars during the first 20 Myr after the formation morphology of MSTOs, we present MSTO properties of 20 of the star cluster. Note that in≈ the latter scenario, the chemi- intermediate-age star clusters in the Magellanic Clouds in this cal enrichment that causes the abundance variations currently paper. This includes 8 such clusters for which new imaging seen among RGB and MS stars in ancient GCs would only data was obtained using the Wide Field Camera #3 (WFC3) have occurred by FRMS and massive binary stars, given the on board the Hubble Space Telescope (HST). time scales involved. This paper is set up as follows. Sect. 2 describes the data An unfortunate issue in distinguishing between these two used in this paper and the star cluster sample. Isochronefitting distinct scenarios for the formation of GCs is the ancient age is described in Sect. 3. Dynamical properties and dynamical of Galactic GCs ( 12–13 Gyr), which preventsa direct mea- evolution of the star clusters in our sample are discussed in surementof the shorttime∼ scales (and hence the types of stars) Sect. 4. Sect. 5 describes pseudo-age distributions of the sam- involved in the chemical enrichment of the “polluted” stars. ple star clusters as derived from the MSTO morphology and In the context of the nature of Na-O anticorrelations in presents a correlation between the MSTO widths and the es- Galactic GCs, the recent discoveryof extended main sequence cape velocities of the clusters at early times. Sect. 6 discusses turn-offs (hereafter eMSTOs) in intermediate-age (1–2 Gyr our findings in the context of predictions of currently popular old) star clusters in the Magellanic Clouds (Mackey et al. scenarios on the nature of eMSTOs, and Sect. 7 summarizes 2008a; Glatt et al. 2008; Milone et al. 2009; Goudfrooij et al. our conclusions. 2009) has generated much interest in the literature, espe- cially since many investigations concluded that the simplest 2. CLUSTER SAMPLE AND NEW DATA viable interpretation of the eMSTOs is the presence of mul- The selection procedure for our “full” target cluster sam- tiple stellar populations spanning an age interval of several ple is based on integrated-light photometry in the literature: 108 yr within these clusters (see also Rubele et al. 2010, Clusters are selected to have a “S parameter” (Girardi et al. 2011; Goudfrooijet al. 2011a,b; Conroy & Spergel 2011; 1995; Pessev et al. 2008) in the range 35–40 along with an Keller et al. 2011; Mackey et al. 2013). However, the eMSTO unreddened integrated-light V-band magnitude 12.5. These phenomenon has been interpreted in two other main ways: criteria translate to cluster ages between roughly≤ 1.0 and 2.0 4 spreads in rotation velocity among turnoff stars (hereafter the Gyr and masses & 3 10 M⊙. The global properties of the “stellar rotation” scenario: Bastian & de Mink 2009; Li et al. star clusters in our “full”× sample are listed in Table 1. At 2012; Yang et al. 2013, butsee Girardi et al. 2011), anda pho- the onset of this study, data of adequate quality was already tometric feature of interacting binaries within a simple stel- available in the HST archive for several star clusters in this lar population (the “interacting binaries” scenario: Yang et al. sample. This includes star clusters in the HST programs GO- 2011; Li et al.