arXiv:1206.0733v1 [astro-ph.GA] 4 Jun 2012 trfrainhsoyi w ed fteSalMglai C Magellanic Small the of fields two in history formation neIsuue hc soeae yAR n. under Inc., GO-10396. asso- Sci- program AURA are with observations Telescope by ciated These operated Space 5-26555. is NAS the NASAcontract which at Instutute, obtained ence Telescope, Space 02,Italy) 40127, ai,PiaeBg3,Hbr,Tsai 01 Australia 7001, Tasmania Hobart, 37, Bag Private mania, 92 edleg Germany Heidelberg, 69120 M¨onchhofstr. 12–1 Universit¨at Heidelberg, der tronomie Irsadams bten1ad5 and 1 of (between majority mass the a of and that dIrr’s to similar current ) young a and has clos- systems, of the (Z class metallicity is this SMC of The member est dwarf (dIrr’s). of evolution the irregular study to laboratory mental INTRODUCTION 1. 1 2 4 3 7 6 5 h ml aelncCod(M)i funda- a is (SMC) Cloud Magellanic Small The .Cignoni M. ae nosrain ihteNS/S Hubble NASA/ESA the with observations on Based srnm et nvriyo oon,Bologna, Bologna, of University Dept. Astronomy colo ahmtc hsc,Uiest fTas- of University Physics, & Mathematics of School Italy 40127, Observatory, Bologna INAF srnmshsRce-nttt etu ¨rAs- f¨ur Zentrum Rechen-Institut, USA Astronomisches WI, Madison, Wisconsin, of University USA 21218, MD, Baltimore, STScI, h iigo h ao psd nteLreMglai lu ssuggest is headings: Cloud Subject Magellanic o Large timing the in the in of formation. between episode level star similarity major low triggering striking the s very The of of the timing rate ago. share the the Gyr but in 5 activity, both around past differ increase Gy over 4 fields recent to examined of 6 two from ratio The enhancement dramatic then. negligib a a since by experienced activity followed indepe Bar Gyr, two the few star using that first find Bar the the We the investigate photometry. stellar of HST their we portion of deep NE diagrams paper of color-magnitude the this the in to In and applied SW procedures the one. in studied located least fields the also but content h a stems rdciergo fteSalMglai lu nte in Cloud Magellanic Small the of region productive most the is Bar The ≃ 0 2 . , 3 0 sdrvdfo I regions HII from derived as 004 .A Cole A. A. , aais vlto aais niiul ml Magella Small individual: galaxies: - evolution galaxies: 4 .Tosi M. , × 10 [email protected] 9 M Grebel ⊙ e.g. , 3 ABSTRACT .S Gallagher S. J. , 4, Bar 7 .Nota A. , 1 1 aeafrSres(C)o ee intermediate- Advanced seven of HST’s (ACS) Surveys with for Camera photometry deep sented gradients. the age as well and as chemical SFH of global existence the will constrain data to These us SMC. allow the field across of clusters sample and large stars a observe to Survey VLT (VST) Telescope the the We and (HST), (VLT), Telescope Telescope Large space. Space Very Hubble and the com- ground using by the to are view from of is datasets field strategy large bining a Our over power spectroscopic and resolving and 2008). space depth in photometric al. high SMC et achieve Tosi the (see of evolution time the study project to long-term international SMC an the started of have we (SFH) history formation star resolved star the circumstance interactions. of a by driven modulation system, activities the triple formation studying a favors of SMC that the member Moreover a range. their is at of therein) limit references upper and the 2006, al. et Kallivayalil o h lse nlss ehv led pre- already have we analysis, cluster the For spatially- detailed, the derive to order In 5 6 .Sabbi E. , ouain eovdb means by resolved populations ta ciiyadissudden its and activity itial esa omto ciiyin activity formation star le g n eryconstant nearly a and ago r i lu,glxe:stellar galaxies: Cloud, nic omto itr ftwo of history formation a omto n nthe in and formation tar 5 v facoeencounter close a of ive h nacmn and enhancement the f .Anderson J. , dn n eltested well and ndent m fsa formation star of rms 5 .K. E. , loud age and old populous clusters (Glatt et al. 2008a,b, Chiosi & Vallenari (2007) derived the SFH in the 2009, 2011). In combination with our VLT data vicinity of a few SMC clusters. Finally, No¨el et al. we find a complex age-metallicity relation for these (2007) and No¨el et al. (2009) presented a deep clusters with a considerable spread in metallicity ground-based study of 12 fields of the SMC, avoid- at any given age (see e.g. Glatt et al. 2008b). ing the densest regions, because of their high Concerning the SMC field analysis, our plan crowding conditions. is to have Color-Magnitude diagrams (CMDs) In this paper we present the SFH of SFH1 and several magnitudes fainter than the oldest main- SFH4, the two most central fields of the six SMC sequence (MS) turn-off (TO) for the entire . regions we observed with HST/ACS (Sabbi et al. To this purpose we have observed six fields with 2009). The apparent distances from the SMC op- the ACS (Sabbi et al. 2009), sampling regions tical center are about 24′ and 1◦ 52′ respectively. characterized by different stellar and gas densi- SFH1 is located in the SW portion of the SMC ties in the SMC Bar, in the Wing in the direction Bar, where the stellar density, gas and dust con- of the LMC, and in the outskirts (see Fig. 1). A tents are highest, while SFH4 is located to the NE preliminary SFH analysis of such fields has been of the SMC center, at 24′ south of NGC 346, the presented in Sabbi et al. (2009). We also have most active star-forming region in the SMC. an ongoing Guaranteed Time Observation pro- For a better assessment of the intrinsic theo- gram at the VST (Ripepi et al. 2006) designed retical uncertainties, the SFH is derived using two to cover with deep photometry the whole SMC completely independent procedures for the appli- and the Bridge connecting it to the LMC. These cation of the synthetic CMD method. We com- CMDs will allow us to infer for the first time the pare the two methods here and discuss the corre- SFH of the whole SMC over the entire Hubble sponding results. The other ACS fields observed time, covering a much larger area with consid- by us and described by Sabbi et al. (2009) are be- erably better image quality than Zaritsky et al.’s ing treated in the same way, and their SFH will (2002) data. We will derive the SFHs from the be presented in a forthcoming paper. CMDs using the synthetic CMD technique (see We briefly describe our data in Section 2. The e.g. Tolstoy et al. 2009; Cignoni & Tosi 2010, and two procedures for the SFH derivation are sum- references therein). SFHs of some SMC fields marized in Section 3, together with the results of have been derived and presented by other au- their application to SFH1 and SFH4. Similarities thors, based on ground-based observations or HST and differences between the resulting SFHs are dis- studies of small individual regions (see Fig. 1). cussed in Section 4, while in Section 5 we compare Harris & Zaritsky (2004) derived the SFH of the ◦ ◦ our findings with previous literature. Concluding SMC over 4 ×4.5 to a depth of V . 21 using the remarks follow in Section 6. Magellanic Cloud Photometric Survey (MCPS) UBVI catalog Zaritsky et al. (2002). This is cur- 2. Photometry and CMDs of SFH1 and rently the most spatially extended study of the SFH4 galaxy, but does not reach the oldest MSTO. The most comprehensive study of the old population of The data for our six SMC fields were ac- the SMC to date was carried out by Haschke et al. quired with the ACS Wide Field Channel between (2012) using RR Lyrae stars from the Optical November 2005 and January 2006 (GO-10396; P.I. Gravitational Lensing Experiment (OGLE-III; Gallagher) with the F555W and F814W filters. Udalski et al. 2008). They find a uniform metal- The data reduction was performed with the pro- licity distribution across the SMC with a spread of gram img2xym WFC.09x10 (Anderson & King more than 1 dex in [Fe/H]. Dolphin et al. (2001) 2006), and the resulting magnitudes were cali- analyzed the stellar content of the SMC halo, in brated in the Vegamag photometric system using a region close to the globular cluster NGC 121, Sirianni et al. (2005) recipes. For sake of simplic- using both HST Wide Field Planetary Camera 2 ity, from now on we will refer to the mF 555W and (WFPC2) and ground based data. Again with mF 814W magnitudes calibrated in the Vegamag WFPC2 McCumber et al. (2005) studied the stel- system as V and I, respectively. lar content of a small portion of the SMC Wing. Extensive artificial star experiments were per-

2 Fig. 1.— Spatial distribution of the six observed fields (red symbols) together with the observations from Dolphin et al. (2001) (blue symbols), McCumber et al. (2005) (magenta symbol), No¨el et al. (2007) (green symbols), Chiosi & Vallenari (2007) (yellow symbols), superimposed on the DSS image of the SMC.

3 formed to test the level of completeness and the 3. SFH derivation: the two synthetic photometric errors of the data. They followed the CMD approaches approach described in Anderson et al. (2008), and the artificial stars were searched for with exactly The SFH of SFH1 and SFH4 has been de- the same procedure adopted for the real stars. We rived with the synthetic CMD method. As considered an artificial star as recovered if its in- recently reviewed by Tolstoy et al. (2009) and put and output positions agree to within 0.5 pixels Cignoni & Tosi (2010), this approach was first and the fluxes agree to within 0.75 mag. As done applied to nearby resolved galaxies twenty years in the photometric analysis we also required that ago (Tosi et al. 1991) and is now recognized as each star was found in at least three exposures a powerful tool to disentangle the complex his- with a positional error < 0.1 pixels per filter. tory of resolved galaxies. To both fields we have applied it following two independent procedures: Details on both the data reduction and the arti- Cole’s (e.g. Cole et al. 2007) and a combination ficial star tests can be found in Sabbi et al. (2009). of Cignoni’s (e.g. Cignoni et al. 2006) and the The final SFH1 and SFH4 catalogs contain Bologna (e.g. Tosi et al. 1991; Greggio et al. 1998; about 29200 and 17300 stars respectively, and the Angeretti et al. 2005) codes. These tecniques su- corresponding CMDs are shown in Fig.2, where 1 persede the classical isochrone fitting approach, the photometric errors are also plotted . These allowing us to explore a much wider parameter CMDs reach almost four magnitudes fainter than space and to incorporate statistical and obser- the oldest MSTO, thus allowing to study the evo- vational uncertainties. Nonetheless their use is lution of the regions over the whole Hubble time. limited by the reliability of the adopted stellar The sub-giant branch (SGB) is well populated and tracks and color transformations, as well as the its brightness extension, much larger than the pho- nature of the IMF (e.g., §4.1) and quality of the tometric error at that magnitude level, suggests a data. The recovered history is the best of all pos- prolonged star formation between 3 and 12 Gyr sibilities, but this does not necessarily imply that ago. The red giant branch (RGB) and the red it is the actual solution nor that it is unique. clump (RC) also are well populated. In addition to The procedures used here have already been these intermediate-age and old components, both tested in comparisons with each other and similar CMDs show a MS blue plume and a blue loop (BL) methods. Both groups participated in the Coim- sequence typical of late-type dwarf galaxies, cor- bra experiment in 2001 (see Skillman & Gallart responding to young high- and intermediate-mass 2002, and references therein), Cignoni’s method stars. has been applied to the derivation of the Solar Finally, the CMDs do not show any signifi- Neighborhood SFH (Cignoni et al. 2006; Cignoni cant population of stars on the right of the lower 2006) and in combination with the Bologna MS, the CMD locus occupied by pre-main se- code to the analysis of the SMC young clusters quence (PMS) stars. This, in turn, suggests a NGC 346 (Cignoni et al. 2010, 2011) and NGC 602 negligible activity in the last 50 Myr (the aver- (Cignoni et al. 2009), Cole’s method has been age time spent by a solar-like star in PMS; see, compared with two different codes in the anal- e.g., Cignoni et al. 2010) compared to more active ysis of the dwarf irregular IC 1613 (Skillman et al. regions of the SMC like NGC 346 or NGC 602 2003) and the Cetus dwarf spheroidal (Monelli et al. (see e.g. Nota et al. 2006; Cignoni et al. 2009; 2010). Since the procedures are in continuous evo- Carlson et al. 2011; Cignoni et al. 2011). lution, to improve both the reliability of the results and the efficiency of the computations, it is useful 1 To be conservative, the plotted error bars correspond at to perform further comparisons on deep and tight each magnitude level to the larger value of the error result- ing from the photometric package and from the artificial CMDs such as those of SFH1 and SFH4. star tests. In the following we briefly summarize common- alities and differences of the two methods. In both cases the synthetic CMDs have been built using the results of the artificial star tests men- tioned above to assign photometric errors and in-

4 Fig. 2.— CMDs of the SFH1 (left panel) and SFH4 (right panel) field observed with the ACS/WFC. The solid and the dashed red lines indicate the 50% and 25% levels of completeness, respectively. Formal errors on the estimated photometry are shown on the left side of each CMD (see text for details).

5 completeness to the synthetic stars. tometry of one of 18 age bins and one of five metallicities (Z = 0.008, Z = 0.004, Z = 0.002, 3.1. Bologna procedure Z = 0.001, Z = 0.0004). In order to reduce the Poisson fluctuations, all partial CMDs are gener- The Bologna approach involves the compari- ated with more than 106 solar masses. son of the observed CMD with a library of syn- thetic CMDs computed with different values of Within the framework of the adopted stellar metallicity, initial mass function (IMF), binary tracks and atmosphere models, the most likely so- fraction, distance modulus and reddening. Mod- lution to the underlying SFH is the one which min- els have been calculated with the latest Padova imizes the differences between data and synthetic stellar models (Bertelli et al. 2008, Bertelli et al. star counts over strategic regions of the CMD. 2 2009) for masses between the hydrogen burn- The degree of likelihood is assessed through a χ ing limit and 20 M⊙. Theoretical temperature function of the residuals. Our experience suggests and luminosity are transformed to the observa- that the performance of such minimization is very tional plane by means of the relations obtained by sensitive to the adopted CMD binning scheme. Origlia & Leitherer (2000) for the HST Vegamag Both fine and coarse grids offer advantages, as well photometric system. The synthetic CMD is cre- as disadvantages (see e.g. Cignoni et al. 2006). ated following a classical Monte Carlo procedure: Along the MS, a fine grid is mandatory to study 1) stellar masses and ages are randomly extracted the old stellar generations, since they are tightly from a time independent IMF and a star formation packed together in the CMD. However, such a so- law; 2) the stellar tracks are interpolated deriving lution would pay the penalty of under-weighting the absolute photometry for the synthetic pop- the star counts along the upper MS which is Pois- ulation; 3) corrections for the distance modulus son dominated. Vice versa, a coarse grid would and the foreground reddening are applied. Then be more adequate to map the recent activity, but the synthetic CMD is degraded to match both it would allow a worse resolution at early times. the completeness profile and the photometric er- The situation is even worse with the evolved stars, ror properties as derived from extensive artificial where the theoretical uncertainties are typically star tests (Sabbi et al. 2009). To be conserva- larger than for MS stars: the CMD position of tive, we have limited our fitting procedure to stars horizontal branch (HB) and RC stars is affected brighter than V = 23, whose completeness is bet- by a complex interplay of age, metallicity, mass ter than 40% and photometric errors smaller than loss and helium content (see e.g. Castellani et al. 0.05 mag. 2000); not only the color but also the shape of the RGB are strongly affected by the color trans- The full SFH of a galaxy can be a complex func- formations from the theoretical to the observa- tion of time. To make the problem manageable, we tional plane and by the mixing length parameter were forced to limit the range of parameter space (in a way which depends on the used wavelength); that our models cover on the basis of previous re- the BL morphology is affected by the He burning sults and indications from our data. In order to cross sections (especially the 12C(α,γ)16O) as well reduce the computational time, the star formation as by the efficiencies of external convection, over- rate is parametrized as a set of constant values shooting and mass loss. Nevertheless, the study of over adjacent temporal steps: a generic CMD is a evolved stars can prevent the ambiguities in mod- linear combination of chosen basis CMDs, where els restricted to MS star counts, for example, when each basis is a Monte Carlo extraction from a step photometric errors and incompleteness hinder the star formation episode. The duration of each step possibility to trace the old activity by means of low is chosen in relation to the evolutionary timescale mass MS stars. Given these issues, the problem of of the average stellar mass of the step, ranging the grid choice can be solved only by using a vari- from 50 Myr (approximately the MS timescale for able grid spacing. Several tests were conducted for an 8 M⊙ star) to 1 Gyr (as representative of the finding the optimal configuration. Figure 3 shows theoretical precision at 10 Gyr, i.e. 10%) for ages our best scheme. To capture all the information above 3 Gyr. The star formation rate is constant contained along the MS, the grid spacing shrinks within each step. In this work we deal with 90 with luminosity, balancing the longer evolutionary basis CMDs, each representing the synthetic pho-

6 Fig. 3.— CMD grid used to derive the SFH via the Bologna procedure.

times of lower mass stars and, at the same time, on the fit quality). To reduce the parameter space, providing sufficiently good statistics in the upper only three metallicities, Z =0.004, Z =0.002 and MS. Towards the red part of the CMD the average Z =0.001, have been adopted. The distance mod- grid spacing is coarser: this helps to take the un- ulus and the foreground reddening were allowed to certainties into account, while preserving the num- vary freely. The solutions with the lowest reduced ber of stars in specific phases. χ2 are listed as SFH1-A and SFH4-A for the re- Finally, the χ2 is minimized by means of a spective fields. downhill simplex algorithm. In order to avoid local In order to test how much the details of the minima, the simplex is re-started from thousands solution depend on the allowed metallicity range, of initial random positions. A bootstrap method we calculated a second solution where the Padova is used to assess the effect of random errors. The tracks with Z = 0.0004 and Z = 0.008 were also in- search of the best SFH is repeated for each boot- cluded. These solutions are presented as SFH1-B strapped data set, producing a distribution of best and SFH4-B. These are generally consistent with solutions. The errorbars on the final SFH repre- the solutions obtained with the restricted metal- sent one standard deviation using 100 bootstraps. licity range but provide a somewhat better match With this procedure we have analyzed the CMD to the observed CMD as described in §4.1 below. of SFH1 & SFH4. As a first step we have de- rived the SFH assuming a two-exponent power law 3.2. Simulated annealing procedure (Cole’s IMF (see Kroupa 2001) with exponent s=2.3, close method) 2 to Salpeter’s 1955 for masses above 0.6 M⊙ and The code developed by Cole has many features s =1.3 below, and a binary fraction (with primary in common with the Bologna code, but also incor- and secondary mass extracted from the same IMF) porates some differences in the treatment of the of 30% (this parameter has no substantial impact Hess diagram data and in the calculation of the merit function that measures the relative likeli- 2Salpeter’s IMF has exponent s=2.35 and the form +∞ −s hood of various solutions. Because the errors in R0 mφ(m) dm = 1 and φ(m) = m .

7 extracting detailed information from a CMD are a metallicity of Z = 0.004 based on HII region dominated by systematic effects, the application and Cepheid metallicities for the youngest stellar of multiple fit procedures to the same data can populations in the SMC Bar. For stars with ages provide important insight into which features of on the order of ≈1–7×108 yr, this metallicity also the SFH are robust and which may be artifacts. gives a good match between the colors of the blue Cole’s SFH-fitting code begins with a set of the- supergiant stars in SFH1 and the Padova models. oretical isochrones interpolated to a fine grid of For older ages, ranges of metallicity at each age age and metallicity in order to create a synthetic are allowed. CMD with no gaps. The most recent isochrones A synthetic CMD is constructed by convolving from the Padova models are used; all calculations the weighted, color-magnitude binned isochrones are made in the HST Vegamag system to minimize with the color and magnitude errors and incom- the possibility of introducing errors in transforma- pleteness functions determined by artificial star tion equations. The synthetic CMDs are binned tests. In the fit process, linear combinations of the in age and metallicity (Z) to increase the speed of individual synthetic CMDs are added to find the computing and to avoid “overfitting” noise in the composite CMD that best matches the data. Be- CMDs. We begin with bins that are evenly spaced cause many cells in the Hess diagram are empty by ≈0.10 in ∆log(age) and 0.20 in ∆log(Z), and or contain few stars, a maximum likelihood test merge adjacent bins when the noise level of the based on the Poisson distribution must be used in- solutions indicates there is insufficient information stead of a χ2 statistic (Cash 1979). A direct search content in the CMD to resolve the fine age bins. of parameter space is infeasible because there may The CMD is divided into a regular grid of color- be dozens of age and metallicity bins in the so- magnitude cells, and the expectation value of the lution; additionally, effects such as age-metallicity number of stars in each cell for a SFR of 1 M⊙ degeneracy can easily produce a large number of yr−1 is calculated from the isochrones. There are false, local maxima in the likelihood space. Be- several parameters that are taken as fixed con- cause of this we use a simulated annealing ap- stants during the solution. These include the proach, in which a simplex of initial guesses at IMF, fraction of binaries and mass ratio distri- the SFH is randomly perturbed and the changes bution function, and the distance modulus and are rejected with some probability if they worsen reddening. The adopted IMF is from Chabrier the fit. The star formation rates are transformed (2003), and the binary fraction and mass ratios according to an arcsin function prior to fitting in are parametrized based on Duquennoy & Mayor order to prevent negative star formation rates from (1991) and Mazeh et al. (1992). In this prescrip- being considered. tion, 35% of stars are single and the rest are bi- Errorbars on the SFH are calculated by test- nary. The binaries are divided into “wide” and ing each age-metallicity bin of the best-fit solu- “close” binaries in a 3:1 ratio; the secondaries in tion in turn. The SFR in the bin under considera- the wide systems are drawn from the same IMF tion is forced away from the optimal solution and as the primaries, but in the close systems the sec- the fitting procedure is redone, with the tested ondary masses are drawn from a flat IMF. The bin held fixed. The 1σ error on the SFR of the distances and reddenings are initially constrained bin is taken to be the value beyond which no fit to the values given in Sabbi et al. (2009), but are can be found that is within 1σ of the global best varied if the resulting synthetic CMDs are mis- fit. Because the total number of stars is fixed, matched to the data. a deficiency of stars in one age bin can often be No age-metallicity relation is explicitly as- partially compensated by increasing the SFR in sumed, but a range of metallicities at each age nearby bins; this means that errorbars are fairly is allowed, constrained by the color range of the large and the SFR in adjacent bins can be strongly data. In some cases, there is little leverage in anticorrelated. Because there are always unmod- the CMD to constrain the metallicity, so out- elled populations (including Galactic foregrounds, side information is used to choose the metallicity. background galaxies, and simple bad data) and For example, the V−I color of the upper MS is there may be poorly-fit stellar sequences (e.g., the not strongly metallicity-dependent, so we adopt colors of RGB stars), the overall fit quality as mea-

8 sured by the equivalent of a χ2 statistic is usually the values adopted here for SFH1. Zaritsky et al. found to be quite poor; only the relative likeli- (2002) discuss the physical reality of their derived hoods are of any meaning. The most likely SFH reddening distributions and conclude that small to match the observed CMD is driven by the most scale reddening variations may reasonably be at- populous cells in the CMD; because these are fre- tributed to residual gas and dust in star-forming quently near the faint end of the data where in- regions. The vast majority of old, cool stars will completeness is high, it is absolutely essential to not be seen through star-forming regions, and so have a very accurate model of the incompleteness the small highly-reddened tail of the old popula- to obtain robust results. tion will not strongly influence the mean redden- The ACS data are of high quality, but incom- ing. Because the amount of differential reddening pleteness sets in at a relatively bright level. The has high spatial frequency, the numerical agree- large number of stars and low average error leads ment (as in SFH1) or lack thereof (as in SFH4) us to use magnitude bins 0.08 mag high and color with the conclusions of Zaritsky et al. (2002) must bins 0.04 mag wide. In order to avoid systemat- be considered fortuitous. ics due to the incompleteness, we restricted the Cole’s best-fit solutions are presented here as comparison to magnitudes 15 ≤ I ≤ 22. Because SFH1-C and SFH4-C respectively for the two there are significant uncertainties in modeling the fields. colors of RGB and RC stars, we restricted our fits to the MS and SGB. We considered a set of 4. Results: SFH of two fields in the SMC isochrones binned by 0.10 in log(age) from 4 Myr Bar ≤ age ≤ 13.5 Gyr (6.60 ≤ log(age) ≤ 10.13); below log(age) = 8.60 the decreasing number of bright The properties of the three SFH solutions for tracer stars led us to double the size of the age each field are summarized in Table 1. While bins. We considered metallicities of Z = 0.00015, the two procedures use the same photometry and 0.0004, 0.0006, 0.001, 0.0015, 0.0024, and 0.004. isochrones, they differ in a number of important respects. The application of two completely in- We begin with the distance moduli and red- dependent codes to derive the SFH allows us to denings from Sabbi et al. (2009), but found that identify the most model-independent features of these needed to be modified in order to obtain the SFH, gives insight into the areas of the CMD the best match to the data with the Padova that are driving the fitting procedures to the re- isochrones. For SFH1 we found the best distance sults they give, and gives an estimate of the sys- modulus (m−M) = 18.83 and reddening E(B−V) 0 tematic uncertainties in the fits. = 0.08, while in SFH4 we used 18.85 and 0.12. In both fields, the values are within 1σ of the initial 4.1. SFH1 field guesses. In SFH1 we found that the model upper MS was significantly narrower and bluer than the The SFH1-A (restricted metallicity) solution data with the adopted reddening value, so we were which gives the lowest (2.6) reduced χ2 is shown forced to introduce an extra reddening component. in Fig. 4 (error bars represent 1-σ uncertainty), We adopted the simple expedient that all popula- while the corresponding synthetic CMD is shown tions younger than log(age) = 7.60 were subject in the middle panel of Fig. 5. The left panel to double the mean reddening of the field, while of the same figure shows the observational CMD stars with 7.60 ≤ log(age) < 8.00 were assigned a for a direct comparison. Our best distance mod- reddening halfway between the younger and older ulus and reddening are (m − M)0 = 18.77 and stars. Differential reddening in the central SMC E(B − V )=0.11, respectively. According to this has previously been observed by Zaritsky et al. solution, the SFR has been extremely low during (2002), who used multicolor photometry to derive the earliest 6-7 Gyr, with a strong and rapid in- line of sight extinction corrections to individual crease around 5 Gyr ago. From then on, the star bright SMC stars. They found a mean AV = 0.18 formation activity has been typically gasping, with mag for stars with 5500 K ≤ Teff ≤ 6500 K, and ups and downs with respect to an average almost AV = 0.46 mag for stars with Teff ≥ 12000 K; constant rate. these numbers are in reasonable agreement with In the middle row of Figure 6 we plot the lu-

9 Table 1: Summary of SFH Solution Parametersa Solution Method (m−M)0 E(B−V) IMF CMDbinning Metallicities 3 (mag) (mag) (color×mag) (Z×10 ) Field SFH1: b SFH1-A Bologna 18.77 0.11 Kroupa (2001) variable 1, 2, 4 SFH1-B Bologna 0.4, 1, 2, 4, 8 c SFH1-C Cole 18.83 0.08 Chabrier (2003) 0.04×0.08 0.15, 0.4, 0.6, 1.0, 1.5, 2.4, 4.0 Field SFH4: b SFH4-A Bologna 18.80 0.11 Kroupa (2001) variable 1, 2, 4 SFH4-B Bologna 0.4, 1, 2, 4, 8 SFH4-C Cole 18.85 0.12 Chabrier (2003) 0.04×0.08 0.15, 0.4, 0.6, 1.0, 1.5, 2.4, 4.0

a All models based on the Padova isochrone set; see text for details. b See Figure 3. c Differential reddening assumed for stars younger than 100 Myr.

minosity function (LF) of blue and red stars (first IMF could mitigate such discrepancy (see the top and second column from the left, respectively) in panel of the first column in Figure 6), but the cor- the observed CMD and in the one resulting from responding model would underestimate the num- SFH1-A. The top and bottom panels show for ber of low mass stars (see the discrepancy between comparison the result when the SFH is searched V = 24 and V = 26 in the top panel of the second using IMF exponents s = 2.0 and 2.5, respectively, column in Fig. 6). for stars above 0.6 M⊙. As far as the evolved stars are concerned, the Despite a somewhat high reduced χ2 2.6, this synthetic CMD matches well both the SGB magni- model is successful in describing the main features tude spread (that is a signature of the age spread) of the data. and the RGB color dispersion below the RC lumi- However when the results are examined in nosity. However, there are several differences as detail a few noticeable issues can be identified. well: 1) the predicted RC morphology is irregu- First of all the synthetic upper MS is sharper larly shaped, while the observational RC is smooth than the observational one. As already noted and rather elliptical. This difference may be par- by Sabbi et al. (2009), such a broadening may tially explained by the coarse metallicity resolu- indicate an additional absorption. One explana- tion of our model, but also by a small amount of tion could be that the SFH1 field is in the main differential reddening; 2) the synthetic CMD over- body of the SMC and the reddening material may estimates the number of RC stars by about a fac- be patchy and cover not uniformly all lines of tor of two; a steeper IMF could slightly mitigate sight. Indeed, this is what is observed in the this mismatch (see the bottom panel of the second extensive analysis of Haschke et al. (2011), who column in Fig. 6) but at the expenses of the up- suggests for SFH1 a broader reddening distribu- per MS fit; 3) the RGB stars brighter than the RC 3 tion (σE(V −I) ≈ 0.12 mag) with respect to SFH4 are too blue (even if the predicted counts match (σE(V −I) ≈ 0.09 mag). Yet other possibilities are exactly the observed ones). This suggests that our that a fraction of massive MS stars is affected by models are too metal poor or, more likely, that the rotation, which can lead to widened upper MSs in adopted color transformations systematically fail young clusters (see e.g. Grebel et al. 1996), or a near the RGB tip. different SMC depth along SFH1 and SFH4’s line We investigate the effect of metallicity on the of sights. SFH found via the Bologna procedure by adding Our model slightly but systematically underes- the Padova tracks with Z =0.008 and Z =0.0004 timates the number of blue massive stars brighter (solution SFH1-B). The top-left panel of Fig. 7 than V ≈ 19 by 20−30 percent (see the middle shows the resulting SFH, while the other panels panel of the first column in Fig. 6). A flatter of the same figure present the mass fraction con- tributed by each metallicity. 3Reddening values from the red clump method are available Overall, the solutions SFH1-A and SFH1-B are through the German Astrophysical Virtual Observatory in- rather similar (see Fig. 8), both being character- terface at http://dc.zah.uni-heidelberg.de/mcx.

10 Fig. 4.— SFH for the SFH1 field, obtained with a restricted metallicity range (solution SFH1-A). The −1 −1 −2 star-formation rate is given both in M⊙ yr (left ordinate) and per unit area (M⊙ yr kpc , in the right ordinate). ized by a long quiet period at the oldest epochs. SFH1-A (compare the middle panel of the third However, solution SFH1-B shows some different column with the middle panel of the first column features. First, the early activity is slightly higher, in Fig.6), but some problems still affect the sim- a direct consequence of the lower initial metallic- ulation: 1) the synthetic CMD shows a hint of a ity; then, the activity between 0.5 and 3 Gyr ago red HB, which is not observed in the data; 2) the is smoother than in SFH1-A. number of BL stars is still under-predicted; 3) the We find with additional test models that the RGB above the RC is too blue; 4) the predicted epoch of the first peak of star formation activity RC is still overpopulated and irregularly shaped. is progressively earlier for decreasing values of the Given these results, especially the improvement adopted lowest metallicity. We are thus confident along the upper MS, model SFH1-B is globally that our conclusion of a long almost quiescent pe- better than SFH1-A, even if neither of them is riod earlier than 5 Gyr ago is robust, although we fully satisfactory. do see and predict some stars as old as 10-12 Gyr. The result of applying Cole’s CMD-fitting pro- Concerning the recovered chemical history, we cedure to SFH1 is shown in Figure 9. The find that the largest stellar mass fraction is pro- left panel shows the data, binned 0.04×0.08 in duced at Z = 0.001, while the contribution from (mV −mI , mI ), while the middle and the right Z = 0.004 is only a few percent. Generally, the panels show the best-fit synthetic CMD (SFH1-C) metallicity increases with time. and map of residuals, respectively. The fit proce- From the point of view of the fit quality, the dure has matched the LF and mean color of the new solution improves the CMD match, yielding a stellar sequences well in general. The data and reduced χ2 of about 2.3. The differences are vis- the model can be readily distinguished from one ible both in the synthetic CMD of Fig. 5 (right another because the data is just one instance of panel) and in the LF (third and fourth columns) a random draw from the parent population and of Fig. 6. Model SFH1-B reproduces the upper contains unmodelled noise, while the model repre- MS star counts and morphology better than model sents the best guess at the pure parent population

11 Fig. 5.— Synthetic CMD for the solution SFH1-A (middle panel) and SFH1-B (right panel) compared with the observational CMD (left panel).

12 Fig. 6.— SFH1: Predicted (red dot-dashed histograms) and observational (black solid histograms) LFs for stars bluer (first and third columns) and redder (second and fourth columns) than V − I =0.6 for the model solutions SFH1-A (first two columns) and SFH1-B (third and fourth colums). Models in the top, middle and bottom panels are computed with IMF exponents above 0.6 M⊙ s = 2.0, 2.3, 2.5, respectively. and is therefore more smoothly distributed. and/or to gaps in modeling the physics of HB en- Among the well-matched features are the color velopes. and vertical extent of the RC and its upward ex- The SFH over the period from 0.5 to 13.5 Gyr tension to I ≈17, the enhancement in the SGB ago for SFH1-C is given in Figure 10. The SFH is at I ≈20.8, and the nearly-vertical finger of stars characterized by a very low star formation rate for at V−I ≈0.5 corresponding to an intermediate-age several billion years before a rapid increase about MSTO. Notable areas of mismatch include the fail- 5 Gyr ago. This is the event that produces the ure to reproduce the width of the MS (most ob- features in the MS and SGB of the SFH1 CMD. vious for I.21), and the factor of 2 overproduc- The SFH remains at a similar elevated level after tion of RC stars (I ≈18.5, V−I ≈1). There is a the 5 Gyr event, with some fluctuations. sparsely populated red HB in the model that is Overall, the SFH of the SFH1 field, based on not apparent in the data; this may be related to the synthetic CMDs from the Bologna and Cole the general factor of 2 overproduction of low-mass procedures, is characterized by the following fea- core helium-burning stars, to poor constraints on tures: the detailed element abundances at ancient times,

13 Fig. 7.— SFH for the SFH1 field resulting from models which include also the metallicities Z = 0.008 and Z = 0.0004 (solution SFH1-B). The top-left panel shows the total SFH while the others display the contribution from each of the labeled metallicities.

• The first 6–8 Gyr were rather quiescent and • About 5–6 Gyr ago SFH1 experienced a re- only a small fraction of the stars in this markable enhancement in the stellar produc- field was formed during these old epochs. tion: over ≈1 Gyr the activity ramped up to −2 −1 −2 The inferred average rate of star forma- about 2.2× 10 M⊙ yr kpc . The age tion for ages older than ≈5 Gyr is only ≈3 and magnitude of this enhancement is ro- −3 −1 −2 ×10 M⊙ yr kpc . bust to choices of IMF, CMD gridding, red- dening, assumed metallicity, and details of

14 Fig. 8.— Recovered SFHs for SFH1 using the Bologna procedure: SFH1-A vs. SFH1-B. The latter is a better match to the data.

the fitting procedure. 4.2. SFH4 field • The average SFH has not dropped signif- The history of the SFH4 field was determined icantly from its peak at 5 Gyr ago, but in similar fashion to the SFH1 field, using three the degree of burstiness in the solutions is sets of simulations (see Table 1). The recovered model-dependent. SFH1-B shows a rela- SFH for case SFH4-A is shown in Fig. 12. tively smooth recent history, while SFH1-C As in SFH1, the star formation proceeded at shows a factor of two drop in SFR between a low level until 5–6 Gyr ago, when it rose to al- 3–4 Gyr ago with a subsequent recovery and most the same amplitude of the first peak in SFH1. SFH1-A shows repeated gasps from 0–2 Gyr Afterward, the actvity started a slow but steady ago. The reasons for this range of behavior decline until now and only in the last 50 Myr it are considered in section 4.3 below. reached higher levels (mimicking in this respect the behavior of SFH1-A). Concerning distance and It is worthwhile to quantify the stellar mass reddening, the recovered values ((m−M)0 = 18.80 produced in each age interval (see the cumulative and E(B − V )=0.11) are not significantly dif- mass function in Fig. 11). According to solution ferent from those obtained for SFH1. The cor- SFH1-B, the SFH1 field assembled a small fraction responding reduced χ2, about 1.8, is significantly (11%) of its stellar mass before 6 Gyr ago, while lower than the best value obtained for SFH1. In- a significant fraction (about 40%) was assembled deed, a visual inspection of the synthetic CMD only in the last 2 Gyr. For comparison, SFH1-C (middle panel of Fig. 13) confirms a very good predicts that 17% of the stars were formed prior agreement with the observational counterpart: our to 6 Gyr, and 33% over the past 2 Gyr. best model can reproduce the position and the morphology of the MS, the SGB and the RGB (above and below the RC). The only minor mis- matches are in the BL region, which is underpop-

15 Fig. 9.— Left panel: observed Hess diagram for the SFH1 field of the SMC Bar, created by binning the data by 0.04 in color and 0.08 in magnitude; middle panel: synthetic CMD from the SFH1-C model (see text for details); right panel: difference Hess diagram between the data and SFH1-C solution.

ulated in our model, and in the SGB luminosity narrower in color, and an overproduction of RC distribution, which is slightly more discontinuous stars. The observed Hess diagram, synthetic di- than in the data. agram from solution SFH4-C, and the difference As for SFH1, we expanded the metallicity range between the two are shown in Figure 17. to include the values Z = 0.008 and Z = 0.0004 SFH4-C is broadly similar to SFH1-C, but the and re-derived the SFH. The essential features of mean SFR is lower at all ages, commensurate the new solution (SFH4-B; see Fig. 14) are not with the smaller number of stars in the field (Fig- changed (see Figure 15 for a comparison SFH4- ure 18). It is notable that the onset of significant A vs SFH4-B): there is still a long quiet initial star formation occurs simultaneously in the two period, followed by a rapid enhancement of the fields, to within the ±1 Gyr precision of our data. activity and a subsequent decline. According to However, the SFH4 field shows a decrease in SFR the solution SFH4-B, most of the stars ever formed relative to SFH1 over the last ≈1.5 Gyr (which is in SFH4 have metallicities around Z =0.001. obvious from a comparison of the RC morphology In terms of fit quality, solution SFH4-B im- and upper MS LF of the two fields). proves the match to the data significantly, leading to a reduced χ2 of about 1.3. The corresponding According to the solution SFH4-B (see the cu- synthetic CMD is excellent (Fig.13, right panel) mulative mass function of Figure 19) about 60% of and the only detectable discrepancy is around the the stellar mass has been produced between 6 Gyr BL phase, still slightly underpopulated by the and 2 Gyr ago, which is higher than the fraction of model. Unlike in SFH1, the number of RC stars is mass (48%) produced in SFH1 in the same period. well matched (see the LF in the right panel of Fig. On the other hand, 16% of the total mass was al- 16). The fit to SFH4 using Cole’s method shows ready in place before 6 Gyr ago, which is compara- similar results, with a better fit to the upper MS ble with the production in SFH1 (11%). The for- than in SFH1 because the observed sequence is mer difference is compensated in the last two Gyr,

16 Fig. 10.— Long term history of the SFH1 field from solution SFH1-C. Significant star formation commenced ≈5 Gyr ago and the rate has remained high, with some fluctuation, ever since. when SFH1 has turned into stars about twice the between 4 and 5 Gyr ago, followed by a slow de- mass fraction of SFH4. While in the SFH1 field cline and a very recent burst. The major difference the Bologna solution produced a younger result here is the rate of decline in SFR after the 4–5 Gyr than the Cole solution, in SFH4 the situation is enhancement, with Cole’s solution showing a con- reversed. The solution SFH4-C predicts that only stant activity for about 2.5 Gyr and Bologna’s so- 7% of the stellar mass was in place before 6 Gyr lution showing a faster decline followed by a mild ago, while 29% was astrated in the past 2 Gyr. enhancement around 1.5 Gyr ago. In both fields, the only ages in which the star 4.3. Comparison of the two methods formation rates differ by significantly more than We have derived the SFH of the field stars in the formal errorbars on the solutions are in the the SMC Bar using two different synthetic CMD period 1.5–3 Gyr ago, with Cole’s solution show- techniques. Figure 20 shows the results plotted ing consistently higher SFR. Part of this differ- together on the same scale to illustrate the simi- ence stems from the slightly higher SFR derived larities and differences. For SFH1 both SFHs indi- by Cole’s method at all ages, owing to the different cate a major activity in the last 5 Gyr, an overall mean stellar masses resulting from the different plateau (on average) since then on, and a spike assumed IMFs and binary fractions. The effect in the last 50 Myr. The major differences con- may be exacerbated in the 1.5–3 Gyr age range cern: 1) The exact epochs of the peaks between because of the interplay between age and the dif- 250 Myr and 4 Gyr ago (around 1.5 Gyr ago in ferent metallicity values considered, and the fact Bologna’s solution and 2.5 Gyr ago in Cole’s so- that Cole’s procedure does not consider the RC lution); 2) The rate in the period 8-13 Gyr ago, in calculating the best-fit SFH. Consideration of which is slightly stronger in Cole’s solution. a very large number of metallicity values tends to produce smoother SFHs than those derived using The agreement for the field SFH4 is good as only a few metallicities (e.g., Cole et al. 2009). well: both models predict a star formation onset

17 Fig. 11.— Cumulative mass function in SFH1 according to the solutions SFH1-B and SFH1-C.

Concerning the distance modulus, Cole’s best 10 and 17 kpc according to Glatt et al. 2008b) values (18.83 for SFH1 and 18.85 for SFH4) are which may depend on the direction. slightly higher than Bologna’s (18.77 and 18.80). The age-metallicity relations inferred by both In SFH1 these offsets may be due to the higher methods are consistent, implying a metallicity reddening required by the Bologna method, but slowly increasing with time and the bulk of the in SFH4 the converse is true. Yet the distance stars with metallicities between 0.001 . Z . 0.002. modulus offsets are similar. The offsets are quite The models derived with the Bologna procedure small and can be considered to be within the consistently find better fits when a metallicity of noise, but a possible reason could be the fact Z = 0.008 is included for the younger ages, while that Cole’s method ignores the RGB and RC, such a population is not required in Cole’s models so the distance is essentially a MS-fitting dis- (see Fig. 21). However, Cole adopted a higher dif- tance, while the Bologna method includes RGB ferential reddening in his modeling of SFH1, which and RC information. More interestingly, both dis- has the higher fraction of young stars; see Table 1. tance estimates are shorter than recent determi- nations based on RR-Lyrae ((m − M)0 = 18.90, 4.4. Comparison Between the Two Fields Kapakos et al. 2011) and eclipsing binaries ((m − 4.4.1. Spatially Resolved SFH of the SMC Bar M)0 = 19.11, North et al. 2010), but still com- patible with the average distance of star clusters Direct comparison of the two fields shows that (around 18.87 for Glatt et al. 2008b and between colors and magnitudes of the prominent stellar se- 18.71 and 18.82 for Crowl et al. 2001). This vari- quences are similar, and the main differences con- ance may indicate a different sensitivity to redden- cern the younger stellar populations. In SFH1 the ing, which is highly variable across the SMC (see upper MS, BL, and red supergiant region are more e.g. Zaritsky et al. 2002; Haschke et al. 2011), or a populous, indicating a higher SFR at ages of a few line of sight depth effect (up to 4.9 kpc according hundred Myr or less. A challenge for both SFH- to Subramanian & Subramaniam 2009, between fitting methods was to reproduce the breadth of

18 Fig. 12.— Best recovered SFH with restricted metallicity range (solution SFH4-A) for the field SFH4. the upper MS in SFH1, as well as its RC, which (2011) presented orbital calculations for the sys- was overproduced by both procedures. Integrat- tem MW/LMC/SMC, finding evidence for a close ing up the predicted stellar masses from the SFH- encounter about 1.5-2.0 Gyr ago. Analyzing SMC 5 fitting procedures we find that ≈4×10 M⊙ stars clusters, Piatti (2011) also suggested two enhanced 5 formed over the lifetime of SFH1, and 2×10 M⊙ formation processes that peaked at 2 and 5-6 Gyr in SFH4. ago. A comparison with our solutions suggests The general similarity of the two CMDs indi- that the former episode might be associated with cates that their features are representative of the the Bologna (1.5 Gyr) and Cole (2.5 Gyr) inter- general properties of the central SMC over most mediate age peak, while the latter coincides with of its history. Small localized bursts/gasps of star our strong rise 5 Gyr ago. formation may be smoothed away by the drift- The inset panels in Fig. 20 show the recent ing of stars away from their birthplace except for 1 Gyr history of both fields according to Bologna the youngest ages. The dramatic increase in SFR and Cole methods. In the SFH1 field the mean at 5 Gyr is a strong feature in both fields and specific SFR over the past 1 Gyr is nearly as high suggests a global change in environment to one as the Gyr-averaged peak rate during the 5 Gyr that strongly favored star formation throughout event, but in SFH4 the recent SFR is reduced by a the SMC at this time. Backward integration of factor of ≈3 from the intermediate-age peak. How- the Magellanic Clouds orbital paths around each ever, in both fields the recent average specific SFR other and the Milky Way over this length of time is still several times higher than the rate prior to are quite uncertain, and it is therefore difficult to the major episode 5 Gyr ago. The increase in the say if this increase is associated with a tidal inter- very youngest age bin may be an artifact of the fit- action or not. ting process, because the only stars younger than More recent dynamical signatures should be ≈50 Myr in our CMDs are still on the MS - all easier to locate in time, and thus to correlate of the evolved stars are above our bright limit- with global star formation events. Diaz & Bekki and so the code is free to vary the SFR at young

19 Fig. 13.— Synthetic CMD for the solution SFH4-A (middle panel) and SFH4-B (right panel) compared with the observational CMD (left panel).

ages somewhat arbitrarily to fit the residuals left 150 Myr) and the steepening of the MS LF at after older populations are constrained. However, around I ≈19.5–20, with the resulting large num- in literature there are also reports indicating very ber of red core helium-burning stars (vertical RC; recent bursts of activity; see e.g. the analysis by 400–800 Myr). In SFH1 the average recent SFR Indu & Subramaniam (2011) who found peaks at is comparable to the peak of the longterm average 0-10 Myr and 50-60 Myr ago. Stars formed in SFH since the 5 Gyr episode, but in SFH4 the rate these events, however, may not have had time to has dropped by a factor of ≈3 from the peak. If diffuse through the SMC, and thus may not be it is proposed that a tidal interaction between the present in our two small HST fields. SMC and LMC (e.g. Besla et al. 2012) or Milky There are also clear differences between the two Way resulted in a period of increased SFR at 100– fields, with SFH4 showing a lower overall SFR den- 200 Myr as seen in SFH1, then it remains to be −3 −1 −2 sity averaging ≈5×10 M⊙ yr kpc and SFH1 explained why this event is not seen in the SFH4 forming stars at the average rate of 2.8×10−2 field. −1 −2 M⊙ yr kpc with an apparently higher de- gree of burstiness. The features in the CMD that 4.4.2. Chemical Evolution indicate high SFR at the “burst” ages are the Because a range of metallicities is considered large number of MS stars at I ≈17–17.5 (100– at most ages, we can plot the resulting age-

20 Fig. 14.— Best SFH (solution SFH4-B) resulting from models which include also the metallicities Z =0.008 and Z =0.0004. The top-left panel shows the total SFH while the others display the contribution from each of the labeled metallicities.

21 Fig. 15.— Comparison between the solutions SFH4-A and SFH4-B.

metallicity relation for the best-fit solutions. This reddening between the two fits, it is not clear these is shown for the Cole solutions (which use a finer differences are significant. grid of metallicities, see Table 1) in Figure 21, The field star AMR shown in Fig. 21 is con- which shows that the two fields have experi- sistent with Piatti’s (2012) AMR, whose AMR is enced similar longterm chemical evolution his- obtained from CCD Washington CT1 photome- tories. This is not a surprising result because try. Moreover the metallicity values for ages older 9 over time periods longer than ∼10 yr the entire than 1 Gyr match well the spectroscopic field star galaxy should be mixed, but it is a confirmation measurements of red giants in the central SMC that the fit procedure is giving consistent results reported by Parisi et al. (2010). despite the differences between the two fields. For Compared to AMRs derived from clusters, our the youngest populations, only a single metal- result is in good agreement with Da Costa & Hatzidimitriou licity (Z=0.004) has been used, so the errorbars (1998) and Glatt et al. (2008b), while it dif- simply reflect the spacing of the metallicity grid. fers from the cluster AMR of Piatti (2011) and Where multiple metallicities were used, the mass- Parisi et al. (2009). Piatti’s (2011) metallicity dis- averaged mean metallicity is given, and the error persion is much higher than our formal errors at bar shows the range of metallicities within which any age, suggesting that SMC clusters may have 90% of the stars in each age bin are expected originated from a less mixed interstellar medium. to fall. The mean metallicity is quite similar to Parisi et al.’s (2009) AMR rises faster than ours that obtained from the Bologna method results, prior to 10 Gyr, suggesting that SMC clusters with the exceptions that Z=0.008 populations are may have experienced an independent chemical not present in large numbers in the Cole solution, evolution history. and the Bologna solution for SFH4 appears to in- It is worth to remind that we have derived the crease its metallicity more quickly for ages older best AMRs from two small regions of the SMC than ≈6 Gyr. Given the uncertainties in RGB and more fields are needed to confirm our results. model colors and the differences in distance and A detailed comparison of our six AMRs (sampling

22 Fig. 16.— Predicted LF (red histograms) for the solution SFH4-B vs the observed (black histograms) LF. The left-hand and the right-hand histograms are for stars bluer and redder than V − I =0.6 respectively. regions of the Bar, Wing and Halo) with spectro- reach the oldest MSTO. This last point is partic- scopic metallicity measurements will be presented ularly relevant, since studies that do reach it indi- in a future paper. cate, instead, that, although present, stars older than 8 Gyr do not dominate the SMC popula- 5. Comparison with other studies tion (Dolphin et al. 2001; McCumber et al. 2005; No¨el et al. 2007, 2009; Chiosi & Vallenari 2007). Harris & Zaritsky (2004) (hereafter HZ04) were Further support to this is also provided by the rel- the first to apply the synthetic CMD method to atively low number of RR Lyrae stars detected in the derivation of the SMC SFH. They used wide- ◦ ◦ the SMC compared to the LMC (Soszy˜nski et al. field UBVI photometry of an area of 4 × 4.5 2010). (MCPS catalog, Zaritsky et al. 2002), suggesting Dolphin et al. (2001) analyzed an external re- that 50% of SMC stars formed earlier than 8.4 gion close to the globular cluster NGC 121, and Gyr ago, and that the period between 3 and 8.4 concluded that stars in the outskirts of the SMC Gyr ago experienced very little, if any, star for- formed during a broadly peaked episode of star mation activity. Figure 22 shows HZ04’s results formation, with the largest rate between 5 and 8 for their regions including our SFH1 and SFH4 Gyr ago. As expected, this outer field is signifi- compared with our derivations. Both HZ04’s so- cantly older than ours, although the contribution lutions show a significantly higher stellar produc- from the age bins older than 10 Gyr is again a tion prior to 8.4 Gyr ago. At intermediate ages minor fraction of the total mass produced. HZ04’s solution for SFH1 is characterized by a clear and long lull between 2 Gyr and 8.4 Gyr McCumber et al. (2005) studied an external ago, which is in striking contrast with our SFHs, field in the SMC Wing, and found an increas- whereas HZ04’s SFH for SFH4 is stronger and in ing rate from 12 to 4 Gyr ago, and then over the satisfactory agreement with our predictions. In past 1.5 Gyr, with a significantly quieter phase the last 2 Gyr both HZ04’s star formation rates between 4 and 1.7 Gyr ago. This analysis was are much stronger than our findings. not conducted using a statistical significance test- ing and their favored solution was the one (among One should note that the HZ04 SFHs are de- three ansatz models) which best matched the star rived from a much larger field of view (about counts in five strategic CMD regions. Although 16 times wider) and their photometry did not this makes the comparison with our finding diffi-

23 Fig. 17.— Left panel: observed Hess diagram for the SFH4 field; middle panel: synthetic Hess diagram from solution SFH4-C; right panel: difference between the two. cult, two considerations can be made: 1) also in typical SFR density in their fields seems compara- this case the stellar production prior to 10 Gyr ago ble with ours, suggesting that the SMC periphery is low, in agreement with our result; 2) apart from and Bar share a common old population but that a recent burst, the average activity of their SFH is the initiation of significant star formation at 5–6 earlier than in our fields, as expected for a Wing Gyr affected the Bar far more than the external region, which represents a population with inter- regions of the SMC, marking the birth of the Bar mediate age between the Bar and the outskirts. as a distinct feature. No¨el et al. (2007) and No¨el et al. (2009) pro- A more direct comparison can be made with duced one of the most extensive and accurate Chiosi & Vallenari (2007), whose analysis was photometric campaign in the SMC outskirts (see focused on three deep HST/ACS fields located also Nidever et al. 2011). Their 12 fields are suf- in the Bar. In this case, our solutions are in ficiently deep to reach the old MSTO and none good quantitative and qualitative agreement. All of them showed a clear extended HB, representa- Chiosi & Vallenari’s (2007) SFHs show both an tive of a very old and metal poor stellar popu- unambiguous rise between 7 and 5 Gyr ago (with lation. No¨el et al. (2009) used a population syn- the precise epoch varying depending on the spe- thesis technique and recovered a great variety of cific direction) and a negligible earlier activity. SFHs, with only the outer fields characterized by However, despite the low rate, their and our early a strong old and intermediate age activity. Al- activities are not zero, offering a natural expla- though the spatial distribution of these fields, nation to the relatively low number of RR-Lyrae mainly enclosed in the periphery of the SMC, stars found in the SMC. Moreover, no two solu- makes the comparison with our results less direct, tions are alike, suggesting that the star formation we note that among No¨el et al.’s (2009) solutions, rate at young and intermediate ages was strongly their innermost field (qj0112) most closely resem- variable across the SMC Bar. bles our history of the Bar. Prior to 8 Gyr ago, the

24 Fig. 18.— Long term history of the SFH4 field from solution SFH4-C. Similar to SFH1, the field “switched on” ≈5 Gyr ago. The SFR has been declining over the past ≈1.5 Gyr in this locale.

Our results are rather similar to the SFHs of In spite of these similarities, the SMC and LMC other two Local Group irregulars such as IC 1613 AMRs (see e.g. Carrera et al. 2008, Harris & Zaritsky (Skillman et al. 2003) and Leo A (Cole et al. 2009) remain rather different. While the LMC 2007): both galaxies show an initial quiet pe- metallicity has increased faster both prior to riod followed by a prompt rise in the SF activity 10 Gyr ago and over the last 4 Gyrs (reaching about 5-7 Gyr ago, with a peak between 3 and 6 [Fe/H] ∼−0.7 and [Fe/H]= −0.2, respectively), Gyr and between 2 and 3 Gyr ago, respectively. between these periods it has progressed at much In both cases the recent activity (< 2 Gyr ago) slower rate. is more similar to that in SFH4 than in SFH1, Moreover, it is well known that the LMC clus- presumably because the IC 1613 explored field is ter distribution shows a long gap between 3 and in the outskirt and the Leo A field is large enough 12 Gyr. If this gap is a consequence of the qui- to include central as well as peripheral regions. escent period between 6 and 10 Gyr ago (see e.g. Finally, comparing our SFH with the solutions Harris & Zaritsky 2009), then it is intriguing that found in the literature for the LMC (see e.g. a similarly long quiescent period has no counter- Harris & Zaritsky 2009 and Smecker-Hane et al. part in the SMC cluster distribution (which has 2002) we find an overall similarity. Apart from been steadily increasing with time). the very early SF enhancement predicted by both Harris et al. and Smecker-Hane et al. 6. Summary and conclusions that has no counterpart in our SMC solutions, our results for SFH1 share key properties with This paper is the first of a series devoted to Smecker-Hane et al.’s (2002) findings for the LMC quantitative reconstruction of the SMC SFH from Bar: 1) a lull over the period 6-10 Gyr ago; 2) a the deep HST/ACS observations presented in steep increase at 5-6 Gyr; 3) a gasping regime over Sabbi et al. (2009). Here we have explored the the last 5 Gyr. directions SFH1 and SFH4, both located in the SMC’s Bar, by comparing the observational CMDs

25 Fig. 19.— Cumulative mass function in SFH4 according to the solutions SFH4-B and SFH4-C. with a library of model CMDs incorporating pho- lowing picture. At early times, both fields experi- tometric uncertainties and incompleteness as esti- enced a long quiescent phase characterized by low mated by Sabbi et al. (2009). To provide a robust SFRs, followed by a rapid SF increase around 5-6 characterization of the SFH, the choice of the best Gyr ago. Since then, the mode of star formation model CMD was independently conducted with has been somewhat different in the two fields. In two objective statistical methods, namely Cole’s SFH1 the star formation was gasping and reason- (Cole et al. 2007) and Bologna’s (Cignoni & Tosi ably high up to today. In SFH4 it was smoother 2010) procedures. and slowly declining. Our best simulated CMDs exhibit an overall To account for these differences and similarities good agreement with observational CMDs. The possible explanations are: star counts along the MS and the SGB morphol- Recent burstiness: The different level of bursti- ogy are generally well reproduced, indicating that ness is not surprising because these fields are sep- our recovered SFHs and metallicity distributions arated by a distance (850 pc) larger than most of are reasonable. However, while SFH4 CMD is well the star forming complexes discovered in the SMC fitted, there are some difficulties to reproduce the (see Livanou et al. 2007), thus allowing the recent exact morphology of SFH1’s CMD, especially the activity to fluctuate independently; upper-MS spread and RC counts, which are under- Recent systematic behavior: The stellar density estimated and overestimated respectively. Con- in the SFH4 region is lower than in SFH1. Hence cerning the resulting SFHs, a good consistency the systematic decrease of SF activity in SFH4 is found between the two methods in both fields. may be connected with a minor amount of fuel The only significant difference is the stronger rate available to support it up to today; suggested by Cole’s SFH4 solution between 1.5 Early quiescence and prompt rise: Are the qui- and 3 Gyr ago. escent period and rapid SF increase 5-6 Gyr ago The combination of synthetic CMDs which a global property of the Bar? Recent observations most resembles the observations suggests the fol- of RGB stars have revealed that older stellar com-

26 Fig. 20.— Bologna SFHs (black line) compared to Cole’s (red line). The left panel and the right panel refer to the fields SFH1 and SFH4, respectively. The inset panels zoom in on the last 1 Gyr.

ponents of the SMC have a velocity dispersion of rate of mutual interactions started around 5 Gyr about 27.5 km s−1 (Harris & Zaritsky 2006), high ago. enough to distribute the stars over a large distance Finally, Rafelski & Zaritsky (2005) noted that (of the order of few kpc) from their birth places the spatial distribution of the younger and older within few Gyr. Hence, the low early activity and clusters in the SMC is statistically different, lead- the prompt rise are not peculiarities of our fields ing to the inference that a significant accretion or but global features of the SMC Bar. Moreover, merger event may have taken place around 3−5 Subramanian & Subramaniam (2012) find no ev- Gyr ago. idence for a Bar in older stars, which is consis- The study reported in this paper was the first tent with our low early activity. Further support step in a wider research activity aimed to charac- is also provided by Chiosi & Vallenari (2007), who terize the SMC SFH through deep HST/ACS ob- found similar star formation trends for three other servations. A forthcoming paper will be dedicated Bar fields located around the SMC clusters K 29, to the analysis of the Wing and Halo fields. This NGC 290, and NGC 265. From the theoretical will allow a comparative analysis to look for global point of view it is not clear what mechanism is physical characteristics in the SMC star formation responsible for the rapid rise of stellar produc- process. tion. Was it externally triggered by the MW or the LMC, or self-initiated? The striking similar- ity between the SMC and LMC SFH is suggestive MC and MT have been partially funded with of the former. Pointing in this direction are the contracts COFIS ASI-INAF I/016/07/0 and ASI- recent calculations by Diaz & Bekki (2011) who INAF I/009/10/0. EKG acknowledges partial presented evidence that around 5.5 Gyr ago the funding from Sonderforschungsbereich ”The Milky LMC and SMC were within 160 kpc of the MW Way System” (SFB 881) of the German Research and 200 kpc of each other, therefore arguing for an Foundation (DFG), especially via subproject A2. independent origin of the Clouds. In this scenario Partial support for JSG’s analysis of data from the transition between quiescent and active phases HST program GO-10396 was provided by NASA could be naturally explained in terms of growing through a grant from the Space Telescope Science Institute, which is operated by the Association of

27 Fig. 21.— Best-fit age-metallicity relation of the two Bar fields derived by Cole’s method. The two fields show similar chemical evolution; within the errorbars they are the same. Taken together they paint a consistent picture of average metallicity increasing steadily over time.

Universities for Research in Astronomy, Inc., un- Besla, G., Kallivayalil, N., Hernquist, L., et al. der NASA contract NAS 5-26555. 2012, MNRAS, 2457

REFERENCES Carlson, L. R., Sewilo, M., Meixner, M., et al. 2011, ApJ, 730, 78 Anderson, J., & King, I. R. 2006, PSFs, Pho- tometry, and Astronomy for the ACS/WFC Carrera, R., Gallart, C., Hardy, E., Aparicio, A., (Instrum. Sci. Rep. ACS 2006-01; Baltimore: & Zinn, R. 2008, AJ, 135, 836 STScI) Cash, W. 1979, ApJ, 228, 939 Anderson, J., Sarajedini, A., Bedin, L., King, I.R., Piotto, G., Reid, I.N., Siegel, M., Majewski, Castellani, V., Degl’Innocenti, S., Girardi, L., S.R., Paust, N.E.Q., Aparicio, A., Milone, A.P., Marconi, M., Prada Moroni, P. G., & Weiss, Chaboyer, B., & Rosenberg, A. 2008, AJ, 135, A. 2000, A&A, 354, 150 2055 Chabrier, G. 2003, PASP, 115, 763 Angeretti, L., Tosi, M., Greggio, L., Sabbi, E., Chiosi, E., & Vallenari, A. 2007, A&A, 466, 165 Aloisi, A., Leitherer, C. 2005, AJ, 129, 2203 Bertelli, G., Girardi, L., Marigo, P., & Nasi, E. Cignoni, M., Degl’Innocenti, S., Prada Moroni, P., 2008, A&A, 484, 815 Shore, S.N. 2006, A&A, 459, 783 Bertelli, G., Nasi, E., Girardi, L., & Marigo, P. Cignoni, M. 2006, Ph.D. Thesis, 2009, A&A, 508, 355

28 Fig. 22.— Comparison between our results (Bologna’s solution is in black, Cole’s is in red) and Harris & Zaritsky (2004) solutions (HZ) for SFH1 (left panel) and SFH4 (right panel).

Cignoni, M., Sabbi, E., Nota, A., et al. 2009, AJ, Dolphin, A.E., Walker, A.R., Hodge, P.W., Ma- 137, 3668 teo, M., Olszewski, E.W., Scommer, R.A., & Suntzeff, N.B. 2001, ApJ, 562, 303 Cignoni, M., Tosi, M., Sabbi, E., et al. 2010, ApJ, 712, L63 Dolphin, A. E., et al. 2003, AJ, 126, 187 Cignoni, M, Tosi, M.2010, in Dwarf Galaxies and Duquennoy, A., & Mayor, M. 1991, A&A, 248, 485 Cosmology, Advances in Astronomy (Hindawi), vol. 3 Glatt, K., Gallagher, J. S., III, Grebel, E. K., et al. 2008a, AJ, 135, 1106 Cignoni, M., Tosi, M., Sabbi, E., Nota, A., & Gal- lagher, J. S. 2011, AJ, 141, 31 Glatt, K., Grebel, E. K., Sabbi, E., et al. 2008b, AJ, 136, 1703 Cole, A.A., Skillman, E.D., Tolstoy, E., Gallagher, J.S.III, Aparicio, A. et al 2007, ApJ, 659, L17 Glatt, K., Grebel, E. K., Gallagher, J. S., III, et al. 2009, AJ, 138, 1403 Cole, A.A., Grocholski, A.J., Geisler, D., Sara- jedini, A., Smith, V.V., & Tolstoy, E. 2009, Glatt, K., Grebel, E. K., Jordi, K., et al. 2011, AJ, in “The Magellanic System: Stars, Gas, and 142, 36 Galaxies”, IAUS 256, ed. J.Th. van Loon & J.M. Oliveira, (Cambridge), p. 263 Grebel, E. K., Roberts, W. J., & Brandner, W. 1996, A&A, 311, 470 Crowl, H. H., Sarajedini, A., Piatti, A. E., et al. 2001, AJ, 122, 220 Greggio, L., Tosi, M., Clampin, M., de Marchi, G., Leitherer, C., Nota, A., & Sirianni, M. 1998, Da Costa, G.S., & Hatzidimitriou, D. 1998, AJ, AJ, 504, 725 115, 1934 Harris, J., & Zaritsky, D. 2004, AJ, 127, 1531 Diaz, J. D., & Bekki, K. 2011, arXiv:1112.6191 Harris, J., & Zaritsky, D. 2006, AJ, 131, 2514 Dohm-Palmer, R. C., Skillman, E. D., Mateo, M., et al. 2002, AJ, 123, 813 Harris, J., & Zaritsky, D. 2009, AJ, 138, 1243

29 Haschke, R., Grebel, E. K., & Duffau, S. 2011, AJ, Ripepi, V., et al. 2006, Mem.S.A.It, 9, 267 141, 158 Sabbi, E., Gallagher, J.S., Tosi, M., , Anderson, Haschke, R., Grebel, E. K., Duffau, S., & Jin, S. J., Nota, A., Grebel, E.K., Cignoni, M., Cole, 2012, AJ, 143, 48 A.A. et al . 2009, ApJ, 703, 721 Indu, G., & Subramaniam, A. 2011, A&A, 535, Salpeter, E. E. 1955, ApJ, 121, 161 A115 Sirianni, M., Jee, M. J., Ben´ıtez, N., et al. 2005, Kallivayalil, N., van der Marel, R. P., & Alcock, PASP, 117, 1049 C. 2006, ApJ, 652, 1213 Skillman, E.D. & Gallart, C. 2002, in Observed Kapakos, E., Hatzidimitriou, D., & Soszy´nski, I. HR Diagrams and stellar evolution, T. Lejeune 2011, MNRAS, 415, 1366 & J. Fernandes eds, ASP Conf Ser., 274, 535 Kroupa, P. 2001, MNRAS, 322, 231 Skillman, E.D., Tolstoy, E., Cole, A.A., Dolphin, A.E., Saha, A., et al. 2003, ApJ, 596, 253 Livanou, E., Gonidakis, I., Kontizas, E., et al. 2007, AJ, 133, 2179 Smecker-Hane, T. A., Cole, A. A., Gallagher, J. S., III, & Stetson, P. B. 2002, ApJ, 566, 239 Mazeh, T., Goldberg, D., Duquennoy, A., & Mayor, M. 1992, ApJ, 401, 265 Soszy˜nski, I., Udalski, A., Szyma˜nski, M. K., et al. 2010, Acta Astron., 60, 165 McCumber, M.P., Garnett, D.R., Dufour, R.J. 2005, AJ, 130, 1083 Stanimirovi´c, S., Staveley-Smith, L., Dickey, J.M., Sault, R.J., & Snowden, S.L. 1999, MNRAS, Monelli, M., et al. 2010, ApJ, 720, 1225 302, 417 Nidever, D. L., Majewski, S. R., Mu˜noz, R. R., et Subramanian, S., & Subramaniam, A. 2009, A&A, al. 2011, ApJ, 733, L10 496, 399 No¨el, N. E. D.;,Aparicio, A. , Gallart, C., Hidalgo, Subramanian, S., & Subramaniam, A. 2012, ApJ, S. L., Costa, E., Mndez, R. A. 2009, ApJ, 705, 744, 128 1260 Tolstoy, E., Hill, V., Tosi, M. 2009, ARAA 47, 371 No¨el, N.E.D., Gallart, C., Costa, E., & M´endez, R.A. 2007, AJ, 133, 2037 Tosi, M., Gallagher, J.S.III, Sabbi, E., Glatt, K., Grebel, E.K., Christian, C., Cignoni, M., North, P., Gauderon, R., Barblan, F., & Royer, F. Clementini, G., et al. 2008, in IAU Symp.255, 2010, A&A, 520, A74 L.K. Hunt, S. Madden, R. Schneider eds, Nota, A., Sirianni, M., Sabbi, E., et al. 2006, ApJ, (CUP) p.381 640, L29 Tosi, M., Greggio, L., Marconi, G., & Focardi, P. Origlia, L., & Leitherer, C. 2000, AJ, 119, 2018 1991, AJ, 102, 951 Parisi, M.C., Grocholski, A.J., Geisler, D., Sara- Udalski, A., Soszy´nski, I., Szyma´nski, M. K., et jedini, A., & Clar´ıa, J.J. 2009, AJ, 138, 517 al. 2008, Acta Astron., 58, 329 Parisi, M.C., Geisler, D., Grocholski, A.J., Clar´ıa, Zaritsky, D., Harris, J., Thompson, I. B., Grebel, J.J., & Sarajedini, A. 2010, AJ, 139, 1168 E. K., & Massey, P. 2002, AJ, 123, 855 Piatti, A. E. 2011, MNRAS, 418, L69 Piatti, A. E. 2012, MNRAS, 422, 1109

A Rafelski, M., & Zaritsky, D. 2005, AJ, 129, 2701 This 2-column preprint was prepared with the AAS L TEX macros v5.2.

30