Draft version September 25, 2019 Typeset using LATEX twocolumn style in AASTeX63 Comparing the Quenching Times of Faint M31 and Milky Way Satellite Galaxies Daniel R. Weisz,1 Nicolas F. Martin,2, 3 Andrew E. Dolphin,4, 5 Saundra M. Albers,1 Michelle L. M. Collins,6 Annette M. N. Ferguson,7 Geraint F. Lewis,8 Dougal Mackey,9 Alan McConnachie,10 R. Michael Rich,11 and Evan D. Skillman12 1Department of Astronomy, University of California, Berkeley, CA 94720 USA 2Universit´ede Strasbourg, CNRS, Observatoire astronomique de Strasbourg, UMR 7550, F-67000 Strasbourg France 3Max-Planck-Institut f¨urAstronomie, K¨onigstuhl17, D-69117 Heidelberg Germany 4Raytheon, 1151 E. Hermans Road, Tucson, AZ 85756 USA 5Steward Observatory, University of Arizona, 933 North Cherry Avenue, Tucson, AZ 85721-0065 USA 6Department of Physics University of Surrey, Guildford, GU2 7XH, Surrey UK 7Institute for Astronomy, University of Edinburgh Royal Observatory, Blackford Hill, Edinburgh EH9 3HJ UK 8Sydney Institute for Astronomy, School of Physics, A28, The University of Sydney, NSW 2006 Australia 9Research School of Astronomy and Astrophysics, Australian National University, Canberra, ACT 2611 Australia 10National Research Council, Herzberg Institute of Astrophysics, 5071 West Saanich Road, Victoria BC V9E 2E7 Canada 11Department of Physics and Astronomy, University of California, Los Angeles, PAB 430 Portola Plaza, Los Angeles, CA 90095 USA 12Minnesota Institute for Astrophysics, University of Minnesota, Minneapolis, MN 55441 USA Submitted to APJL ABSTRACT We present the star formation histories (SFHs) of 20 faint M31 satellites (−12 . MV . −6) that were measured by modeling sub-horizontal branch (HB) depth color-magnitude diagrams constructed from Hubble Space Telescope (HST) imaging. Reinforcing previous results, we find that virtually all galaxies quenched between 3 and 9 Gyr ago, independent of luminosity, with a notable concentration 3 − 6 Gyr ago. This is in contrast to the Milky Way (MW) satellites, which are generally either faint with ancient quenching times or luminous with recent (<3 Gyr) quenching times. We suggest that systematic differences in the quenching times of M31 and MW satellites may be a reflection of the varying accretion histories of M31 and the MW. This result implies that the formation histories of low-mass satellites may not be broadly representative of low-mass galaxies in general. Among the M31 satellite population we identify two distinct groups based on their SFHs: one with exponentially declining SFHs (τ ∼ 2 Gyr) and one with rising SFHs with abrupt quenching. We speculate how these two groups could be related to scenarios for a recent major merger involving M31. The Cycle 27 HST Treasury survey of M31 satellites will provide well-constrained ancient SFHs to go along with the quenching times we measure here. The discovery and characterization of M31 satellites with MV & −6 would help quantify the relative contributions of reionization and environment to quenching of the lowest-mass satellites. Keywords: Andromeda Galaxy (39), Stellar populations (1622), Dwarf spheroidal galaxies (420), Local Group (929) arXiv:1909.11036v1 [astro-ph.GA] 24 Sep 2019 1. INTRODUCTION ter cosmology on small scales (e.g., Bullock & Boylan- Milky Way (MW) satellite galaxies have long an- Kolchin 2017; Simon 2019). Their star formation histo- chored our understanding of low-mass mass galaxy for- ries (SFHs) and chemical content provide insight into mation. Their number counts, spatial distributions, and cosmic reionization and the baryonic processes that structural properties are used to constrain dark mat- uniquely shape low-mass galaxy evolution (e.g., Kirby et al. 2011; Brown et al. 2014; Weisz et al. 2014). More recently, their orbital histories, as measured by the Hub- Corresponding author: Daniel R. Weisz ble Space Telescope (HST) and Gaia, reveal the complex [email protected] effects of central galaxies on the evolution of low-mass 2 Weisz et al. satellites (e.g., Sohn et al. 2012; Fritz et al. 2018; Simon HST-specific modules (Dolphin 2000). We adopted the 2018). DOLPHOT parameters recommended in Williams et al. At the same time, there is growing evidence that the (2014). The raw photometric catalogs were culled to MW satellites may not be representative of low-mass include only good stars as described in Martin et al. satellites in general. Compared to the MW, satellite (2017). We ran ∼ 105 artificial star tests (ASTs) for each systems throughout the local Universe show varying lu- galaxy to quantify completeness and photometric uncer- minosity functions, stellar populations, quenching prop- tainties. The 50% completeness limit for a typical galaxy erties, and spatial configurations, often in excess of cos- in our sample is F 606W ∼ 27:1 and F 814W ∼ 26:2. mic variance (e.g., McConnachie & Irwin 2006; Brasseur Figure1 illustrates the quality of our data by show- et al. 2011; Tollerud et al. 2011; Geha et al. 2017; ing a Hess diagram, i.e., a binned CMD, for And XXI M¨ulleret al. 2018; Pawlowski 2018; Smercina et al. 2018; (MV = −9:2). The CMD shows a clear red giant branch Pawlowski et al. 2019). Even in our nearest neighbor, (RGB), red clump (RC), and a predominantly red HB, M31, there are hints that the internal (e.g., kinemat- as described in Martin et al.(2017). The faint limit ics, stellar content; Da Costa et al. 1996; Collins et al. of the CMD is set to the 50% completeness limit of 2010, 2014; Martin et al. 2017) and global (e.g., \plane F 814W = 26:3, which is ∼1.5 magnitudes fainter than of satellites"; Ibata et al. 2013; Pawlowski 2018) prop- the HB. erties of M31 and MW satellites are different. Thus, it is unclear whether the fundamental insights established 3. METHODOLOGY in MW satellites are applicable to all low-mass systems We measure the SFHs of the 20 galaxies in our sample or stem from the specific accretion history of the MW. using MATCH (Dolphin 2002), a software package that In this Letter, we present the first uniform SFH mea- forward models the CMD of a galaxy in order to measure surements of many faint M31 satellites. We use sub- its SFH as described in Weisz et al.(2014). Here, we horizontal branch (HB) HST imaging, which has been provide a brief summary of the pertinent details. presented in two previous papers in this series (Martin For this analysis, we adopt a Kroupa IMF (Kroupa et al. 2017; Weisz et al. 2019), to measure their SFHs, 2001), a binary fraction of 0.35, HB-based distances for and we compare our measurements to literature SFHs of the ACS data (Weisz et al. 2019) and self-consistent MW satellites. Relative to the MW satellites, the SFHs TRGB-based distances for the WFPC2 data (Weisz we measure for M31 satellites have coarser age resolu- et al. 2014), and foreground extinction values from tion, because the M31 CMDs do not reach the oldest Schlafly & Finkbeiner(2011). main sequence turn off (MSTO). We will acquire oldest We fit each entire CMD with five different stellar evo- MSTO photometry for all known M31 satellites during lution libraries: Dartmouth (Dotter et al. 2008), Padova HST Cycle 27 as part of HST-GO-15902 (PI D. Weisz), (Girardi et al. 2010), PARSEC (Bressan et al. 2012), which will strengthen the results we preview in this Let- MIST (Choi et al. 2016), and BaSTI (Hidalgo et al. ter. 2018). We find that for all CMDs analyzed in this paper the BaSTI 2018 models provide the best overall fits in 2. THE DATA terms of visual inspection of the residuals and through The observations and data reduction used in this pro- comparison of likelihood ratios between models. Thus, gram are described in detail in Martin et al.(2017). we adopt the BaSTI models for this paper. Here, we provide a brief summary. We adopt a metallicity grid that ranges from −2:3 ≤ Through HST-GO-13699 (PI N. Martin), we observed [M/H] ≤ −0:5 with a resolution of 0.1 dex and an age 16 faint M31 satellites with the Advanced Camera for grid that ranges from 9:00 ≤ log(t=yr) ≤ 10:15 with a Surveys (ACS) resolution of log(t=yr) = 0:05 dex. We find that includ- that had no previous HST imaging. Each galaxy was ing ages younger than log(t=yr) = 9:0 did not change observed for a single orbit with equal integration times the SFHs (as these galaxies have no young populations) in the F606W and F814W filters. We added archival but increased the computational time. Thus, we sim- F606W and F814W ACS data of comparable depth for ply exclude ages with log(t=yr) < 9:0 from our CMD And XVIII (HST-SNAP-13442; PI B. Tully). We also modeling. added archival observations of And XI, And XII, and Finally, given that our CMDs do not reach the oldest And XIII (HST-GO-11084; PI D. Zucker) which were MSTO, we follow Weisz et al.(2011) and Weisz et al. taken in F606W and F814W with the Wide Field Plan- (2014) and adopt a prior on the age-metallicity relation- etary Camera 2 (WFPC2). ship that requires the metallicity to increase monoton- In total, our sample has 20 systems with −12 . MV . ically with time, with a modest dispersion allowed at −6, including 7 ultra-faint dwarf galaxies (UFDs, MV > each age. This choice helps to mitigate some of the −7:7; Simon 2019). age-metallicity degeneracy on the RGB and RC. We For each galaxy, we performed point spread function compute random (which account for the finite number photometry with DOLPHOT, a widely-used package for re- of stars on the a CMD) and systematic uncertainties ducing observations of resolved stellar populations with (which are a proxy for uncertainties in the physics of the Formation Histories of M31 Satellites 3 Figure 1.