Arxiv:1711.03735V2 [Astro-Ph.GA] 9 May 2018
Total Page:16
File Type:pdf, Size:1020Kb
ApJ in Press Preprint typeset using LATEX style emulateapj v. 04/17/13 SILVERRUSH. V. CENSUS OF LYα, [OIII]λ5007,Hα, AND [CII]158µm LINE EMISSION WITH ∼ 1000 LAES AT z = 4:9 − 7:0 REVEALED WITH SUBARU/HSC Yuichi Harikane1,2, Masami Ouchi1,3, Takatoshi Shibuya1,4, Takashi Kojima1,2, Haibin Zhang1,2, Ryohei Itoh1,2, Yoshiaki Ono1, Ryo Higuchi1,2, Akio K. Inoue5, Jacopo Chevallard6,18, Peter L. Capak7,8, Tohru Nagao9, Masato Onodera10, Andreas L. Faisst8, Crystal L. Martin11, Michael Rauch12, Gustavo A. Bruzual13, Stephane Charlot14, Iary Davidzon8, Seiji Fujimoto1,15, Miftahul Hilmi1,15, Olivier Ilbert16, Chien-Hsiu Lee10, Yoshiki Matsuoka9, John D. Silverman3, and Sune Toft17 ApJ in Press ABSTRACT We investigate Lyα, [Oiii]λ5007,Hα, and [Cii]158µm emission from 1; 124 galaxies at z = 4:9 − 7:0. Our sample is composed of 1; 092 Lyα emitters (LAEs) at z = 4:9, 5:7, 6:6, and 7:0 identified by Subaru/Hyper Suprime-Cam (HSC) narrowband surveys covered by Spitzer large area survey with Subaru/HSC (SPLASH) and 34 galaxies at z = 5:148−7:508 with deep ALMA [Cii]158µm data in the literature. Fluxes of strong rest-frame optical lines of [Oiii] and Hα (Hβ) are constrained by significant excesses found in the SPLASH 3:6 and 4:5 µm photometry. At z = 4:9, we find that the rest-frame Hα equivalent width and the Lyα escape fraction fLyα positively correlate with the rest-frame Lyα 0 0 equivalent width EWLyα. The fLyα − EWLyα correlation is similarly found at z ∼ 0 − 2, suggesting no evolution of the correlation over z ' 0 − 5. The typical ionizing photon production efficiency of −1 LAEs is logξion=[Hz erg ] ' 25:5 significantly (60-100%) higher than those of LBGs at a given UV magnitude. At z = 5:7 − 7:0, there exists an interesting turn-over trend that the [Oiii]=Hα flux ratio 0 ˚ 0 ˚ increases in EWLyα ' 0 − 30 A, and then decreases out to EWLyα ' 130 A. We also identify an 0 anti-correlation between a [Cii] luminosity to star-formation rate ratio (L[CII]=SF R) and EWLyα at the > 99% confidence level. We carefully investigate physical origins of the correlations with stellar- 0 synthesis and photoionization models, and find that a simple anti-correlation between EWLyα and metallicity explains self-consistently all of the correlations of Lyα,Hα,[Oiii]=Hα, and [Cii] identified 0 ˚ in our study, indicating detections of metal-poor (∼ 0:03Z ) galaxies with EWLyα ' 200A. Keywords: galaxies: formation | galaxies: evolution | galaxies: high-redshift [email protected] 1 Institute for Cosmic Ray Research, The University of Tokyo, 5-1-5 Kashiwanoha, Kashiwa, Chiba 277-8582, Japan 2 Department of Physics, Graduate School of Science, The University of Tokyo, 7-3-1 Hongo, Bunkyo, Tokyo, 113-0033, Japan 3 Kavli Institute for the Physics and Mathematics of the Uni- verse (Kavli IPMU, WPI), University of Tokyo, Kashiwa, Chiba 277-8583, Japan 4 Department of Computer Sciences, Kitami Institute of Tech- nology, 165 Koen-cho, Kitami, Hokkaido 090-8507, Japan 5 Department of Environmental Science and Technology, Fac- ulty of Design Technology, Osaka Sangyo University, 3-1-1, Na- gaito, Daito, Osaka 574-8530, Japan 6 Scientific Support Office, Directorate of Science and Robotic arXiv:1711.03735v2 [astro-ph.GA] 9 May 2018 Exploration, ESA/ESTEC, Keplerlaan 1, 2201 AZ Noordwijk, The Netherlands 7 California Institute of Technology, MC 105-24, 1200 East California Blvd., Pasadena, CA 91125, USA 8 Infrared Processing and Analysis Center, California In- stitute of Technology, MC 100-22, 770 South Wilson Ave., Pasadena, CA 91125, USA 9 Research Center for Space and Cosmic Evolution, Ehime University, Bunkyo-cho, Matsuyama, Ehime 790-8577, Japan 10 Subaru Telescope, National Astronomical Observatory of Japan, 650 North Aohoku Place, Hilo, HI 96720, USA dAstrophysique de Paris, F-75014, Paris, France 11 Department of Physics, University of California, Santa Bar- 15 Department of Astronomy, Graduate School of Science, bara, CA, 93106, USA The University of Tokyo, 7-3-1 Hongo, Bunkyo, Tokyo 113-0033, 12 Carnegie Observatories, 813 Santa Barbara Street, Japan Pasadena, CA 91101, USA 16 Aix Marseille Univ, CNRS, LAM, Laboratoire dAstro- 13 Centro de Radioastronomia y Astrofisica (CRyA), UNAM, physique de Marseille, Marseille, France Campus Morelia Apartado Postal 3-72, 58089 Morelia, Michoa- 17 Cosmic Dawn Center (DAWN), Niels Bohr Institute, Ju- can, Mexico liane Mariesvej 30, DK-2100 Copenhagen, Denmark 14 Sorbonne Universites, UPMC-CNRS, UMR7095, Institut 18 ESA Research Fellow 2 Harikane et al. 1. INTRODUCTION veal that the redshifted emission lines significantly af- Probing physical conditions of the inter-stellar medium fect infrared broad-band photometry (e.g., Stark et al. (ISM) is fundamental in understanding star formation 2013; Smit et al. 2014, 2015; Rasappu et al. 2016; Faisst and gas reprocessing in galaxies across cosmic time. Re- et al. 2016a; Roberts-Borsani et al. 2016; Castellano et al. cent ALMA observations are uncovering interesting fea- 2017). Thus, infrared broad-band photometry can be tures of the ISM in high-redshift galaxies. Early obser- useful to estimate the rest-frame optical emission line vations found surprisingly weak [Cii]158µm emission in fluxes which are not accessible with the ground-based Lyα emitters (LAEs) at z ∼ 6 − 7 ([Cii] deficit; e.g., telescopes before the JWST era. Ouchi et al. 2013; Ota et al. 2014; Schaerer et al. 2015; The Subaru/Hyper Suprime-Cam Subaru strategic Maiolino et al. 2015). On the other hand, recent studies program (HSC-SSP) survey started in early 2014, and its detected strong [Cii] emission in galaxies at z = 5 − 7, first data release took place in 2017 February (Miyazaki whose [Cii] luminosities are comparable to local star- et al. 2012 and Aihara et al. 2018a,b; see also Miyazaki forming galaxies (e.g., Capak et al. 2015; Pentericci et al. et al. 2017; Komiyama et al. 2017; Furusawa et al. 2017; 2016; Bradaˇcet al. 2017). A theoretical study discusses Kawanomoto et al. 2017). The HSC-SSP survey provides that the [Cii] deficit can be explained by very low metal- a large high-redshift galaxy sample, especially LAEs se- lected with the narrow-band (NB) filters. The NB816 licity (0:05 Z ) in the ISM (Vallini et al. 2015; Olsen et al. 2017). Thus estimating metallicities of the high- and NB921 imaging data are already taken in the HSC- redshift galaxies is crucial to our understanding of the SSP survey. In addition, the NB718 and NB973 data origin of the [Cii] deficit. are taken in the Cosmic HydrOgen Reionization Unveiled The ISM property is also important for cosmic reion- with Subaru (CHORUS) project (PI: A. K. Inoue; A. K. ization. Observations by the Planck satellite and high Inoue et al. in prep), which is an independent program redshift UV luminosity functions (LFs) suggest that faint of the HSC-SSP survey. Spitzer large area survey with Subaru/HSC (SPLASH; PI: P. Capak; P. Capak et al. in and abundant star-forming galaxies dominate the reion- 19 ization process (e.g., Robertson et al. 2015). Further- preparation.) has obtained the Spitzer/Infrared Array more, Ishigaki et al. (2018) claim that the ionizing pho- Camera (IRAC) images overlapped with these NB data, ton budget of star-forming galaxies is sufficient for reion- which allow us to conduct statistical studies of the rest- izing the universe with the escape fraction of ionizing frame optical emission in the high redshift LAEs. Fur- photons of f = 0:15+0:06 and the faint limit of the UV thermore, the number of galaxies observed with ALMA esc −0:02 is increasing, which will improve our understanding of LF of Mtrunc > −12:5 for an assumed constant ioniz- −1 the [Cii] deficit. Thus in this study, we investigate the ing photon production efficiency of logξion=[Hz erg ] = 25:34, which is the number of Lyman continuum pho- ISM properties of high redshift galaxies by measuring the Lyα,[ ]λ5007, Hα,Hβ, and [ emission line tons per UV (1500 A)˚ luminosity (see also Faisst 2016). Oiii Cii] strength (Figure 1). On the other hand, Giallongo et al. (2015) argue that This paper is one in a series of papers from twin faint AGNs are important contributors to the reioniza- programs devoted to scientific results on high redshift tion from their estimates of number densities and ioniz- galaxies based on the HSC-SSP survey data. One pro- ing emissivities (c.f., Madau & Haardt 2015; Parsa et al. gram is our LAE study with the large-area narrow- 2018). One caveat in these two contradictory results is band images complemented by spectroscopic observa- that properties of ionizing sources (i.e., f and ξ ) esc ion tions, named Systematic Identification of LAEs for Visi- are not guaranteed to be the same as the typically as- ble Exploration and Reionization Research Using Sub- sumed values. Various studies constrain ionizing pho- aru HSC (SILVERRUSH; Ouchi et al. 2018; Shibuya ton production efficiencies of star forming galaxies to et al. 2018a,b; Konno et al. 2018; Inoue et al. 2018; be logξ =[Hz erg−1] = 24:8 − 25:3 at z ∼ 0 − 2 (e.g., ion Higuchi et al. 2018). The other one is a luminous Ly- Matthee et al. 2017a; Izotov et al. 2017; Shivaei et al. man break galaxy (LBG) studies, named Great Optically 2017; see also Sobral et al. 2018). Recently, Bouwens −1 Luminous Dropout Research Using Subaru HSC (GOL- et al. (2016) report logξion=[Hz erg ] = 25:3 − 25:8 for DRUSH; Ono et al. 2018; Harikane et al. 2018; Toshikawa Lyman break galaxies (LBGs) at z ∼ 4 − 5, relatively et al.