Arxiv:2105.11497V2 [Astro-Ph.GA] 17 Aug 2021
Total Page:16
File Type:pdf, Size:1020Kb
Draft version August 18, 2021 Typeset using LATEX twocolumn style in AASTeX63 First HETDEX Spectroscopic Determinations of Lyα and UV Luminosity Functions at z = 2 − 3: Bridging a Gap Between Faint AGN and Bright Galaxies Yechi Zhang,1, 2 Masami Ouchi,3, 1, 4 Karl Gebhardt,5 Erin Mentuch Cooper,5 Chenxu Liu,5 Dustin Davis,5 Donghui Jeong,6, 7 Daniel J. Farrow,8, 9 Steven L. Finkelstein,5 Eric Gawiser,10 Gary J. Hill,11, 5 Yuichi Harikane,1, 12 Ryota Kakuma,1 Viviana Acquaviva,13, 14 Caitlin M. Casey,5 Maximilian Fabricius,8, 9 Ulrich Hopp,8, 9 Matt J. Jarvis,15, 16 Martin Landriau,17 Ken Mawatari,1, 3 Shiro Mukae,1 Yoshiaki Ono,1 Nao Sakai,1 and Donald P. Schneider6, 7 1Institute for Cosmic Ray Research, The University of Tokyo, 5-1-5 Kashiwanoha, Kashiwa, Chiba 277-8582, Japan 2Department of Astronomy, Graduate School of Science, the University of Tokyo, 7-3-1 Hongo, Bunkyo, Tokyo 113-0033, Japan 3National Astronomical Observatory of Japan, 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan 4Kavli Institute for the Physics and Mathematics of the Universe (Kavli IPMU, WPI), The University of Tokyo, 5-1-5 Kashiwanoha, Kashiwa, Chiba, 277-8583, Japan 5Department of Astronomy, The University of Texas at Austin, 2515 Speedway, Stop C1400, Austin, Texas 78712, USA 6Department of Astronomy and Astrophysics, The Pennsylvania State University, University Park, PA 16802, USA 7Institute for Gravitation and the Cosmos, The Pennsylvania State University, University Park, PA 16802, USA 8Max-Planck Institut f¨urextraterrestrische Physik, Giessenbachstrasse 1, 85748 Garching, Germany 9University Observatory, Fakult¨atf¨urPhysik, Ludwig-Maximilians University Munich, Scheiner Strasse 1, 81679 Munich, Germany 10Department of Physics and Astronomy, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA 11McDonald Observatory, University of Texas at Austin, 2515 Speedway, Stop C1402, Austin, TX 78712, USA 12Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, UK 13Physics Department, NYC College of Technology, 300 Jay Street, Brooklyn, NY 11201, USA 14Center for Computational Astrophysics, Flatiron Institute, New York, NY 10010, USA 15Astrophysics, Department of Physics, Keble Road, Oxford, OX1 3RH, UK 16Department of Physics & Astronomy, University of the WesternCape, Private Bag X17, Bellville, Cape Town, 7535, South Africa 17Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, CA 94720, USA ABSTRACT We present Lyα and ultraviolet-continuum (UV) luminosity functions (LFs) of galaxies and active galactic nuclei (AGN) at z = 2:0 − 3:5 determined by the un-targetted optical spectroscopic survey of the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). We combine deep Subaru imaging with HETDEX spectra resulting in 11:4 deg2 of fiber-spectra sky coverage, obtaining 18320 galax- ies spectroscopically identified with Lyα emission, 2126 of which host type 1 AGN showing broad (FWHM > 1000 km s−1) Lyα emission lines. We derive the Lyα (UV) LF over 2 orders of magnitude −1 covering bright galaxies and AGN in log LLyα=[erg s ] = 43:3 − 45:5 (−27 < MUV < −20) by the 1=Vmax estimator. Our results reveal the bright-end hump of the Lyα LF is composed of type 1 AGN. In conjunction with previous spectroscopic results at the faint end, we measure a slope of the best-fit +0:13 Schechter function to be αSch = −1:70−0:14, which indicates αSch steepens from z = 2−3 towards high redshift. Our UV LF agrees well with previous AGN UV LFs, and extends to faint-AGN and bright- galaxy regimes. The number fraction of Lyα-emitting objects (X ) increases from M ∗ ∼ −21 arXiv:2105.11497v2 [astro-ph.GA] 17 Aug 2021 LAE UV to bright magnitude due to the contribution of type 1 AGN, while previous studies claim that XLyα ∗ ∗ decreases from faint magnitude to MUV, suggesting a valley in the XLyα−magnitude relation at MUV. Comparing our UV LF of type 1 AGN at z = 2 − 3 with those at z = 0, we find that the number density of faint (MUV > −21) type 1 AGN increases from z ∼ 2 to z ∼ 0 as opposed to the evolution of bright (MUV < −21) type 1 AGN, suggesting AGN downsizing in the rest-frame UV luminosity. Corresponding author: Yechi Zhang [email protected] ii Zhang et al. Keywords: galaxies: formation | galaxies: evolution | galaxies: high-redshift 1. INTRODUCTION suggesting no evolution of Lyα LF from z = 5:7 − 6:5 At high redshift, Lyman α emitters (LAEs) are a (e.g., Malhotra & Rhoads 2004). widely studied population of objects that feature strong Despite the efforts of previous studies of the Lyα LFs, Lyman α (Lyα λ1216 A)˚ emission lines (e.g., Rhoads several open questions still remain. One open question et al. 2000; Gronwall et al. 2007; Pentericci et al. 2009; is the steepness of the faint-end slope of the Lyα LFs Ouchi et al. 2020). Typical LAEs are interpreted as that describes the fraction of faint galaxies relative to 10 brighter ones. Theoretical models of hierarchical struc- young (. 50 Myr), low-mass (. 10 M ) galaxies with −1 ture formation predict that low-mass galaxies are more high star formation rates (SFR) of ∼ 1 − 100 M yr (e.g., Nagao et al. 2005; Gawiser et al. 2006; Finkelstein dominant at higher redshift, which results in a steeper et al. 2007, 2008, 2009; Ono et al. 2010a,b; Kashikawa faint-end slope. Such a redshift evolution of the faint- et al. 2012; Harikane et al. 2018). Such properties make end slope has been identified in the UV continuum LF LAEs important tracers for galaxy formation at the low- (hereafter UV LF) of LBGs (e.g., Bouwens et al. 2015; mass end of the spectrum in the early universe, comple- Finkelstein et al. 2015). Since the dust attenuation of mentary to the continuum selected Lyman Break Galax- Lyα emission in the interstellar medium (ISM) becomes ies (LBGs) that are relatively massive. larger towards fainter UV luminosity (e.g., Ando et al. A key statistical property of LAEs is the luminosity 2006; Ouchi et al. 2008) and higher redshift (e.g.,Blanc function (LF), which is defined as the number density as et al. 2011; Hayes et al. 2011), the observed faint-end a function of luminosity. The LAE LFs and their evolu- slope of Lyα LF is expected to be steeper than that of tion can provide valuable insights into the evolution of UV LF towards higher redshift. However, the faint-end young, star forming (SF) galaxies over the cosmic time. slope of the Lyα LF is poorly constrained in the previous Over the past decades, LAEs have been identified and studies due to the large uncertainty of contamination in studied with LFs over the redshift range z ∼ 0 − 9 us- photometric LAE samples and the limited spectroscopic ing deep narrow band (e.g., Gronwall et al. 2007; Ouchi LAE samples. et al. 2008, 2010; Konno et al. 2016, 2018; Tilvi et al. Another open question in relation to the Lyα LF is 2020) and spectroscopic surveys (e.g., Blanc et al. 2011; the shape of the bright end. Several studies have at- Cassata et al. 2011; Zheng et al. 2013; Drake et al. tempted to statistically characterize the bright LAEs −1 2017; Herenz et al. 2019). These studies have found with log LLyα=[erg s ] & 43:5 based on photometrically 43 −1 selected samples, reaching various conclusions. For ex- that at Lyα luminosity LLyα . 10 erg s , the Lyα LF of LAEs can be described by the Schechter function ample, Konno et al.(2016) identified an excess in the (Schechter 1976), number density with respect to the Schechter function at the bright end of their Lyα LFs at z = 2:2. Such an excess was also found in Lyα LFs over other red- φ (L)dlogL = Sch shift ranges (e.g., Wold et al. 2017; Zheng et al. 2017). αSch+1 ∗ L L Sobral et al.(2017, 2018a), and Matthee et al.(2017) ln 10 φSch ∗ exp − ∗ dlogL: (1) LSch LSch have shown a similar but less significant bright-end ex- cess at z = 2 − 3 and demonstrated that such an ex- ∗ ∗ The values of LSch, αSch, and φSch represent the char- cess can be fitted by a power law. At similar red- acteristic luminosity, the faint end slope, and the char- shifts, Spinoso et al.(2020) found ∼ 14; 500 LAEs with acteristic number density, respectively. Comparing Lyα −1 log LLyα=[erg s ] > 43:3. The Lyα LFs of their bright LFs at different epochs also shows the redshift evolution LAEs are described by Schechter exponential decays of LAEs. From z ∼ 0 to z ∼ 3, the number density of ∗ −1 with log LSch=[erg s ] ∼ 44:5 − 44:8. Despite these LAEs increases rapidly (e.g., Deharveng et al. 2008). At efforts, the precise shape of the Lyα LF, especially at z ∼ 3 − 6, there is little evolution in the number density the bright end, is in need of further constraints from of LAEs (e.g., Dawson et al. 2004; Ouchi et al. 2008). spectroscopic studies. Beyond z ∼ 6, the observed number density of LAEs Along with the shape of the Lyα LF at the bright begins to decrease due to the resonant scattering of Lyα end, it is also important to understand the nature of photon by the increasing neutral hydrogen (HI) fraction the extreme objects that cause the bright-end excess. in the intergalactic medium (IGM) towards the epoch Although the differential attenuation of Lyα photons in of reionization (e.g., Kashikawa et al. 2006; Hu et al. the clumsy ISM may cause the Lyα LFs to have a non- 2010; Itoh et al.