Draft version July 5, 2021 Typeset using LATEX twocolumn style in AASTeX63 The mass-metallicity relation at z ∼ 1 − 2 and its dependence on star formation rate Alaina Henry,1, 2 Marc Rafelski,1, 2 Ben Sunnquist,1 Norbert Pirzkal,1 Camilla Pacifici,1 Hakim Atek,3 Micaela Bagley,4 Ivano Baronchelli,5 Guillermo Barro,6 Andrew J. Bunker,7, 8 James Colbert,9 Y.Sophia Dai,10 Bruce G. Elmegreen,11 Debra Meloy Elmegreen,12 Steven Finkelstein,4 Dale Kocevski,13 Anton Koekemoer,1 Matthew Malkan,14 Crystal L. Martin,15 Vihang Mehta,16 Anthony Pahl,14 Casey Papovich,17 Michael Rutkowski,18 Jorge Sanchez´ Almeida,19 Claudia Scarlata,16 Gregory Snyder,1 and Harry Teplitz9 1Space Telescope Science Institute; 3700 San Martin Drive, Baltimore, MD, 21218, USA 2Department of Physics & Astronomy, Johns Hopkins University, Baltimore, MD 21218, USA 3Institut d'astrophysique de Paris, CNRS UMR7095, Sorbonne Universit´e,98bis Boulevard Arago, F-75014 Paris, France 4Department of Astronomy, The University of Texas at Austin, Austin, TX, 78712, USA 5Dipartmento di Fisica e Astronomia, Universit´adi Padova, Vicolo dell'Osservatorio, I-3 31522, Padova, Italy 6Department of Physics, University of the Pacific, 3601 Pacific Ave, Stockton, CA, 95211, USA 7Sub-department of Astrophysics, Department of Physics, University of Oxford, Denys Wilkinson Building, Keble Road, Oxford OX1 3RH, UK 8Kavli Institute for the Physics and Mathematics of the Universe (WPI), The University of Tokyo Institutes for Advanced Study, The University of Tokyo, Kashiwa, Chiba 277-8583, Japan 9Infrared Processing and analysis Center, California Institute of Technology, Pasadena, CA 91125, USA 10Chinese Academy of Sciences South America Center for Astronomy (CASSACA), NAOC, 20A Datun Road, Beijing, 100101, China 11IBM Research Division, T.J. Watson Research Center, 1101 Kitchawan Road, Yorktown Heights, NY10598, USA 12Department of Physics & Astronomy, Vassar College, Poughkeepsie, NY 12604, USA 13Department of Physics and Astronomy, Colby College, Waterville, ME 04901, USA 14Department of Physics and Astronomy, UCLA, Los Angeles, CA 90095-1547, USA 15Department of Physics, University of California, Santa Barbara, CA 93106, USA 16Minnesota Institute for Astrophysics, University of Minnesota, 116 Church Street SE, Minneapolis, MN 55455, USA 17Department of Physics and Astronomy, Texas A&M University, College Station, TX 77843-4343, USA 18Department of Physics and Astronomy, Minnesota State University, Mankato, MN, 56001, USA 19Instituto de Astrof´ısica de Canarias, La Laguna, Tenerife 38200, Spain ABSTRACT We present a new measurement of the gas-phase mass-metallicity relation (MZR), and its dependence on star formation rates (SFRs) at 1:3 < z < 2:3. Our sample comprises 1056 galaxies with a mean redshift of z = 1:9, identified from the Hubble Space Telescope Wide Field Camera 3 (WFC3) grism spectroscopy in the Cosmic Assembly Near-Infrared Deep Extragalactic Survey (CANDELS) and the WFC3 Infrared Spectroscopic Parallel Survey (WISP). This sample is four times larger than previous 8 metallicity surveys at z ∼ 2, and reaches an order of magnitude lower in stellar mass (10 M ). Using stacked spectra, we find that the MZR evolves by 0.3 dex relative to z ∼ 0:1. Additionally, we identify a subset of 49 galaxies with high signal-to-noise (SNR) spectra and redshifts between 1:3 < z < 1:5, where Hα emission is observed along with [O III] and [O II]. With accurate measurements of SFR in these objects, we confirm the existence of a mass-metallicity-SFR (M-Z-SFR) relation at arXiv:2107.00672v1 [astro-ph.GA] 1 Jul 2021 high redshifts. These galaxies show systematic differences from the local M-Z-SFR relation, which vary depending on the adopted measurement of the local relation. However, it remains difficult to ascertain whether these differences could be due to redshift evolution, as the local M-Z-SFR relation is poorly constrained at the masses and SFRs of our sample. Lastly, we reproduced our sample selection in the IllustrisTNG hydrodynamical simulation, demonstrating that our line flux limit lowers the normalization of the simulated MZR by 0.2 dex. We show that the M-Z-SFR relation in IllustrisTNG has an SFR dependence that is too steep by a factor of around three. 1. INTRODUCTION evolution. The growth of galaxies{ and the proper- The role of the baryon cycle in galaxy evolution is one ties that we observe today{ are determined by gas ac- of the most critical questions facing studies of galaxy cretion from the intergalactic medium (IGM) and cir- 2 cumgalactic medium (CGM), as well as star-formation lower metallicities, while those with lower SFRs have feedback that heats and removes gas from the interstel- more enriched ISM. Mannucci et al.(2010) dubbed this lar medium (ISM). Theoretical studies of galaxy forma- relation the \fundamental metallicity relation" (FMR). tion, both from hydrodynamical simulations and semi- Using the limited data available at the time, they argued analytic models, must balance these processes to pro- that the FMR did not evolve out to z ∼ 2:5; galaxies duce realistic properties of galaxies over cosmic history merely evolve within it. From a theoretical perspective, (Somerville & Dav´e 2015). the interpretation that the FMR is produced by varia- One of the key constraints on galaxy formation mod- tions in accretion, SFR, gas fractions, and feedback is els is the correlation between stellar mass and galaxy supported by both semi-analytic models and numerical metallicities (both stellar and gas-phase1), known as the simulations (Yates et al. 2012; Dayal et al. 2013; De mass-metallicity relation (MZR). Since metals are pro- Rossi et al. 2017; Torrey et al. 2018, 2019; Dav´eet al. duced by star formation, their production is also regu- 2017, 2019; De Lucia et al. 2020). lated by feedback. A large number of models, from an- In the low redshift universe (z . 0:3), the MZR alytical approaches (Dalcanton 2007; Erb 2008; Peeples and FMR have been well-characterized down to ap- 8 & Shankar 2011; Dav´eet al. 2012; Dayal et al. 2013; proximately 10 M , using the large sample of galaxies Lilly et al. 2013; Yabe et al. 2015b), to semi-analytical covered by the spectroscopic Sloan Digital Sky Survey models (Lu et al. 2014; Croton et al. 2016) and hydro- (SDSS; Tremonti et al. 2004; Andrews & Martini 2013; dynamical simulations (Finlator & Dav´e 2008; Ma et al. Brown et al. 2016; Curti et al. 2020). Additionally, ex- 2016; De Rossi et al. 2017; Dav´eet al. 2017, 2019; Tor- tensions of the MZR to low masses have been made us- rey et al. 2018, 2019), have all modeled the MZR. What ing nearby galaxies (Lee et al. 2006; Berg et al. 2012; is clear is that a large number of physical processes may Izotov et al. 2012; Ly et al. 2016). Moreover, at inter- contribute to the shape, normalization, and evolution of mediate redshifts, obtaining statistical samples is still the MZR. Supernova-driven outflows are usually consid- practical. The MZR derived from optical multi-object 8 9 ered a key ingredient, setting the slope of the MZR at spectroscopic surveys reaches M∗ ∼ 10 − 10 M for low masses (Tremonti et al. 2004; Henry et al. 2013a,b; samples of around 1000 galaxies at z ∼ 0:5 − 1:0 (Zahid De Rossi et al. 2017). These outflows may be mixed et al. 2011; Maier et al. 2015; Guo et al. 2016a). How- with the ISM, or metal-enriched from the sites of star ever, at z > 1, observational constraints on the evolution formation and supernovae (e.g. Chisholm et al. 2018). of the MZR require near infrared spectroscopy, and have Other studies have argued that the evolution of the MZR therefore been slower to develop (but see early single is largely tied to increasing gas fractions in lower mass object spectroscopy studies; Erb et al. 2006; Maiolino et and higher redshift systems (Ma et al. 2016; Torrey et al. 2008; Wuyts et al. 2012; Belli et al. 2013; Masters et al. 2019). Feedback from active galactic nuclei (AGN) in al. 2014). Nevertheless, in recent years, infrared grism simulations has also been shown to play an important spectroscopy with the Wide Field Camera 3 on the Hub- role in shaping the MZR at high masses (De Rossi et ble Space Telescope (MacKenty et al. 2010), along with al. 2017). Complicating all of these effects, gas is likely ground-based multi-object infrared spectrometers have recycled through the CGM; ejected gas may not leave opened up a new window on the evolution of galaxy the halo and re-accreted gas may be metal-deficient or metallicities, extending measurements of statistical sam- metal-enriched (S´anchez Almeida et al. 2014). ples to z ∼ 2 (Henry et al. 2013b; Zahid et al. 2014; Nonetheless, evidence of an interplay between inflows, Sanders et al. 2015; Kacprzak et al. 2016; Grasshorn star formation, metal production, and feedback is ob- Gebhardt et al. 2016; Wuyts et al. 2016; Salim et al. served. The scatter in the MZR is shown to correlate 2015; Yabe et al. 2012, 2015a; Kashino et al. 2017; with star-formation rates (SFRs) and (in some cases) Sanders et al. 2018; Papovich et al., in prep). No- gas fractions in galaxies out to z ∼ 2 (Ellison et al. tably, numerous studies show that a mass-metallicity- 2008; Mannucci et al. 2010, 2011; Cresci et al. 2012; SFR (M-Z-SFR) relation, similar to the local FMR, is Lara-L´opez et al. 2010, 2013; Bothwell et al. 2013; An- also present out to z ∼ 2. drews & Martini 2013; Henry et al. 2013a; Salim et al. One of the challenges faced by high redshift (z ∼ 1−2) 2014, 2015; Sanders et al. 2018; Gillman et al.