Draft version May 19, 2021 Typeset using LATEX twocolumn style in AASTeX63 A recently quenched isolated dwarf galaxy outside of the Local Group environment Ava Polzin,1 Pieter van Dokkum,1 Shany Danieli,1, 2, 3, 4, ∗ Johnny P. Greco,5, y and Aaron J. Romanowsky6, 7 1Department of Astronomy, Yale University, New Haven, CT 06511, USA 2Department of Physics, Yale University, New Haven, CT 06520, USA 3Yale Center for Astronomy and Astrophysics, Yale University, New Haven, CT 06511, USA 4Institute for Advanced Study, 1 Einstein Drive, Princeton, NJ 08540, USA 5Center for Cosmology and AstroParticle Physics (CCAPP), The Ohio State University, Columbus, OH 43210, USA 6Department of Physics & Astronomy, One Washington Square, San Jos´eState University, San Jose, CA 95192, USA 7University of California Observatories, 1156 High Street, Santa Cruz, CA 95064, USA (Received March 22, 2021; Revised May 14, 2021; Accepted May 17, 2021) Submitted to ApJL ABSTRACT 6 We report the serendipitous identification of a low mass (M∗ ∼ 2×10 M ), isolated, likely quenched dwarf galaxy in the \foreground" of the COSMOS-CANDELS field. From deep Hubble Space Telescope (HST) imaging we infer a surface brightness fluctuation distance for COSMOS-dw1 of DSBF = 22 ± 3 Mpc, which is consistent with its radial velocity of cz = 1222 ± 64 km s−1 via Keck/LRIS. At this distance, the galaxy is 1.4 Mpc in projection from its nearest massive neighbor. We do not detect significant Hα emission (EW(Hα)= −0:4 ± 0:5 A),˚ suggesting that COSMOS-dw1 is likely quenched. Very little is currently known about isolated quenched galaxies in this mass regime. Such galaxies are thought to be rare, as there is no obvious mechanism to permanently stop star formation in them; to date there are only four examples of well-studied quenched field dwarfs, only two of which appear to have quenched in isolation. COSMOS-dw1 is the first example outside of the immediate vicinity of the Local Group. COSMOS-dw1 has a relatively weak Dn4000 break and the HST data show a clump of blue stars indicating that star formation ceased only recently. We speculate that COSMOS-dw1 was quenched due to internal feedback, which was able to temporarily suspend star formation. In this 6 7 scenario the expectation is that quenched isolated galaxies with masses M∗ = 10 − 10 M generally have luminosity-weighted ages . 1 Gyr. Keywords: Dwarf galaxies (416), Galaxy quenching (2040), Quenched galaxies (2016), Galaxy evolution (594) 1. INTRODUCTION can be used to determine how typical Local Group galax- Dwarf galaxies in the Local Group play a pivotal role ies are for the general population (Weisz et al. 2011a). in many areas of astrophysics, including star formation, Most general galaxy catalogs are biased against the arXiv:2105.08061v1 [astro-ph.GA] 17 May 2021 galaxy formation, and cosmology (see, e.g., Mateo 1998; lowest mass dwarfs due to incompleteness (e.g., Tollerud Simon 2019). Given their outsized importance and im- et al. 2008; Walsh et al. 2009), which is caused by their pact there is great interest in finding dwarf galaxies at low surface brightness (see Danieli et al. 2018). To over- larger distances, as they provide information on the envi- come this barrier, various telescopes and surveys have ronmental dependence of low mass galaxy formation and been designed to be sensitive to low surface brightness emission, making use of specialized algorithms and in- struments to detect/characterize low luminosity galaxies Corresponding author: Ava Polzin (e.g., Greco et al. 2018; van Dokkum et al. 2020). As a [email protected] result, wide field surveys are now providing comprehen- ∗ NASA Hubble Fellow sive catalogs of low luminosity galaxies. y NSF Astronomy & Astrophysics Postdoctoral Fellow 2 Polzin et al. Past studies of star formation in dwarf galaxies have red, likely post-main sequence, stars that are distributed been subject to these same limitations. For instance, more evenly. while Geha et al.(2012) were able to study complete 8 2.2. GALFIT samples of galaxies with M∗ & 10 M with SDSS spec- troscopy, they probe only a very small volume at lower We begin by measuring the global structural param- 7 masses. Galaxies with masses M∗ . 10 M have not eters of the galaxy, such as its apparent size, bright- yet been exhaustively analyzed. ness, and color, using GALFIT (Peng et al. 2010). We The interpretation of detected low surface brightness first ran GALFIT on the combined F814W+F606W im- objects generally requires ancillary data, such as spec- age, after smoothing it with a Gaussian kernel with troscopy or high resolution imaging (see, e.g., Greco σ = 2 pix (000: 1). The galaxy was modeled with a single et al. 2021). Archival data offers a shortcut: if dwarf S´ersicfit and the image was aggressively masked, using galaxies are sought and found in fields that already a SEP (Bertin & Arnouts 1996; Barbary 2016) segmen- have a suite of ancillary data, the task of determin- tation map and manually masking additional potential ing distances, structural parameters, and stellar popu- contaminants. The resulting S´ersicindex is n = 0:69, lation parameters is far more efficient. Here we present and the effective radius is 400: 20. Next, GALFIT was run the serendipitous discovery of an isolated and seemingly in each band separately, holding the reff , S´ersicindex, quenched dwarf galaxy in the well-studied COSMOS- position angle, ellipticity, and x; y position fixed and let- CANDELS field. ting only the brightness vary in the fit. The resultant parameters of our GALFIT runs are listed in Table1. Er- rors are determined in the following way. Eleven copies 2. A FAINT, EXTENDED OBJECT IN THE of the best-fitting GALFIT model were injected into un- COSMOS FIELD crowded areas in our images and then fitted in the same 2.1. HST observations way as the actual data. The rms variation in the result- The COSMOS-CANDELS (Koekemoer et al. 2011; ing parameters was taken as the uncertainty for each Grogin et al. 2011) field covers ∼ 200 square arcmin parameter. This method captures errors that are intro- and is one of the most observed regions in the sky, with duced by improper masking of background sources and data taken in all major wavelength regimes from X-ray noise. However, it assumes that the galaxy is smooth to radio. and has a perfect S´ersicprofile, and it does not take COSMOS-dw1 (α = 10h00m30:03s, δ = deviations from those assumptions into account. +02◦08059:4700) was identified in archival Hubble Space 3. VELOCITY AND DISTANCE Telescope (HST ) data of the COSMOS field as an ob- ject with a semi-resolved appearance, suggesting that it 3.1. Keck Spectroscopy is nearby. The galaxy had been previously overlooked, We observed COSMOS-dw1 with the Low Resolution probably because it is split in multiple very faint objects Imaging Spectrograph (LRIS, Oke et al. 1995) on Keck in standard catalogs (Lilly et al. 2007; Alam et al. 2015). I on 2018 November 5. The 100: 5 long slit was used, with We obtained the CTE-corrected individual flc files in the 300 l mm−1 grism blazed at 5000 A.˚ The total ex- the ACS/WFC F475W, F606W, and F814W bands from posure time was 4600 s in excellent conditions. The the HST Archive. We use DrizzlePac (STSCI Develop- data reduction followed standard procedures for long ment Team 2012) to align the images in each filter via slit data. The spectrum is shown in Fig.2. The most TweakReg and then combine them via AstroDrizzle. prominent features are strong Balmer absorption lines, 1 Total exposure times are 2028 s in F475W, 3328 s in indicating a dominant population of A stars and an age F606W, and 6864 s in F814W. of ∼ 1 Gyr. There are no clearly detected emission lines. An RGB image combining the HST data is shown We use a χ2 minimization scheme over the spectral in Figure1. Several notable features are apparent. range 3875 − 5200 A˚ to determine the radial veloc- COSMOS-dw1 is a low surface brightness, semi-resolved ity, with Flexible Stellar Population Synthesis object. Its appearance is somewhat asymmetric, as (FSPS; Conroy et al. 2009) template spectra smoothed it has a clump of blue stars off-center to the North. to the instrumental resolution. We measure a heliocen- The rest of the galaxy appears to be dominated by tric radial velocity of 1222 ± 64 km s−1; the uncertainty in this result includes fits for a range of ages (0.5{8 Gyr) 1 We did not use existing data products from CANDELS or 3D- and metallicities (−2 ≤ [Fe/H] ≤ −1). HST for consistency with our analysis of individual flc files in § This velocity is consistent with COSMOS-dw1 being 4. +0:8 at a distance of 20:0−0:7 Mpc via the Cosmicflows-3 COSMOS-dw1 3 Figure 1. Both panels show color composite image of COSMOS-dw1 using F475W, F606W, and F814W. The left panel is 6800 × 6800; the right panel is 1800 × 1800. −1 −1 Table 1. COSMOS-dw1 { Observed (Kourkchi et al. 2020; H0 = 75 km s Mpc ,ΩM = properties 0:27, and ΩΛ = 0:73) calculator. However, this uncer- tainty ignores the peculiar velocity of COSMOS-dw1; a parameter value peculiar velocity of 300 km s−1 would correspond to a mF475W 19:43 ± 0:07 distance uncertainty of 5 Mpc. mF606W 19:31 ± 0:04 mF814W 19:17 ± 0:03 3.2. Surface brightness fluctuations F475W−F814W 0:26 ± 0:08 A distance can also be obtained from the HST imag- F606W−F814W 0:14 ± 0:05 −2 ing. The galaxy is only semi-resolved, and we cannot µ0;g (mag arcsec ) 23:03 ± 0:08 −2 obtain a distance from the tip of the red giant branch µ0;V (mag arcsec ) 22:90 ± 0:06 (see Section 4.3).