Showing 1–50 of 243 results for author: Kashikawa, N

JWST, ALMA, and Keck Spectroscopic Constraints on the UV Luminosity Functions at z~7-14: Clumpiness and Compactness of the Brightest Galaxies in the Early Universe

Authors: Yuichi Harikane, Akio K. Inoue, Richard S. Ellis, Masami Ouchi, Yurina Nakazato, Naoki Yoshida, Yoshiaki Ono, Fengwu Sun, Riku A. Sato, Seiji Fujimoto, Nobunari Kashikawa, Derek J. McLeod, Pablo G. Perez-Gonzalez, Marcin Sawicki, Yuma Sugahara, Yi Xu, Satoshi Yamanaka, Adam C. Carnall, Fergus Cullen, James S. Dunlop, Eiichi Egami, Norman Grogin, Yuki Isobe, Anton M. Koekemoer, Nicolas Laporte , et al. (10 additional authors not shown)

Abstract: We present the number densities and physical properties of the bright galaxies spectroscopically confirmed at $z\sim7-14$. Our sample is composed of 53 galaxies at $z_\mathrm{spec}\sim7-14$, including recently-confirmed galaxies at $z_\mathrm{spec}=12.34-14.32$ with JWST, as well as new confirmations at $z_\mathrm{spec}=6.583-7.643$ with $-24< M_\mathrm{UV}< -21$ mag using ALMA and Keck. Our JWST/… ▽ More We present the number densities and physical properties of the bright galaxies spectroscopically confirmed at $z\sim7-14$. Our sample is composed of 53 galaxies at $z_\mathrm{spec}\sim7-14$, including recently-confirmed galaxies at $z_\mathrm{spec}=12.34-14.32$ with JWST, as well as new confirmations at $z_\mathrm{spec}=6.583-7.643$ with $-24< M_\mathrm{UV}< -21$ mag using ALMA and Keck. Our JWST/NIRSpec observations have also revealed that very bright galaxy candidates at $z\sim10-13$ identified from ground-based telescope images before JWST are passive galaxies at $z\sim3-4$, emphasizing the necessity of strict screening and spectroscopy in the selection of the brightest galaxies at $z>10$. The UV luminosity functions derived from these spectroscopic results are consistent with a double power-law function, showing tensions with theoretical models at the bright end. To understand the origin of the overabundance of bright galaxies, we investigate their morphologies using JWST/NIRCam high-resolution images obtained in various surveys including PRIMER and COSMOS-Web. We find that $\sim70\%$ of the bright galaxies at $z\sim7$ exhibit clumpy morphologies with multiple sub-components, suggesting merger-induced starburst activity, which is consistent with SED fitting results showing bursty star formation histories. At $z\gtrsim10$, bright galaxies are classified into two types of galaxies; extended ones with weak high-ionization emission lines, and compact ones with strong high-ionization lines including NIV]$λ$1486, indicating that at least two different processes (e.g., merger-induced starburst and compact star formation/AGN) are shaping the physical properties of the brightest galaxies at $z\gtrsim10$ and are responsible for their overabundance. △ Less

Submitted 26 June, 2024; originally announced June 2024.

Comments: 32 pages, 22 figures, submitted to ApJ

MAMMOTH-Subaru. II. Diverse Populations of Circumgalactic Ly$α$ Nebulae at Cosmic Noon

Authors: Mingyu Li, Haibin Zhang, Zheng Cai, Yongming Liang, Nobunari Kashikawa, Ke Ma, Xiaohui Fan, J. Xavier Prochaska, Bjorn H. C. Emonts, Xin Wang, Yunjing Wu, Shiwu Zhang, Qiong Li, Sean D. Johnson, Minghao Yue, Fabrizio Arrigoni Battaia, Sebastiano Cantalupo, Joseph F. Hennawi, Satoshi Kikuta, Yuanhang Ning, Masami Ouchi, Rhythm Shimakawa, Ben Wang, Weichen Wang, Zheng Zheng , et al. (1 additional authors not shown)

Abstract: Circumgalactic Lyman-alpha (Ly$α$) nebulae are gaseous halos around galaxies exhibiting luminous extended Ly$α$ emission. This work investigates Ly$α$ nebulae from deep imaging of $\sim12~\mathrm{deg}^2$ sky, targeted by the MAMMOTH-Subaru survey. Utilizing the wide-field capability of Hyper Suprime-Cam (HSC), we present one of the largest blind Ly$α$ nebula selections, including QSO nebulae, Ly… ▽ More Circumgalactic Lyman-alpha (Ly$α$) nebulae are gaseous halos around galaxies exhibiting luminous extended Ly$α$ emission. This work investigates Ly$α$ nebulae from deep imaging of $\sim12~\mathrm{deg}^2$ sky, targeted by the MAMMOTH-Subaru survey. Utilizing the wide-field capability of Hyper Suprime-Cam (HSC), we present one of the largest blind Ly$α$ nebula selections, including QSO nebulae, Ly$α$ blobs, and radio galaxy nebulae down to typical $2σ$ Ly$α$ surface brightness of $(5-10)\times10^{-18}\mathrm{~erg~s^{-1}~cm^{-2}~arcsec^{-2}}$. The sample contains 117 nebulae with Ly$α$ sizes of 40 - 400 kpc, and the most gigantic one spans about 365 kpc, referred to as the Ivory Nebula. Combining with multiwavelength data, we investigate diverse nebula populations and associated galaxies. We find a small fraction of Ly$α$ nebulae have QSOs ($\sim7\%$), luminous infrared galaxies ($\sim1\%$), and radio galaxies ($\sim 2\%$). Remarkably, among the 28 enormous Ly$α$ nebulae (ELANe) exceeding 100 kpc, about $80\%$ are associated with UV-faint galaxies ($M_\mathrm{UV} > -22$), categorized as Type II ELANe. We underscore that Type II ELANe constitute the majority but remain largely hidden in current galaxy and QSO surveys. Dusty starburst and obscured AGN activity are proposed to explain the nature of Type II ELANe. The SED of stacking all Ly$α$ nebulae also reveals signs of massive dusty star-forming galaxies with obscured AGNs. We propose a model to explain the dusty nature where the diverse populations of Ly$α$ nebula capture massive galaxies at different evolutionary stages undergoing violent assembling. Ly$α$ nebulae provide critical insights into the formation and evolution of today's massive cluster galaxies at cosmic noon. △ Less

Submitted 21 May, 2024; originally announced May 2024.

Comments: 26 pages, 10 figures, 3 tables, submitted to ApJS, comments welcome

Discovery of merging twin quasars at z = 6.05

Authors: Yoshiki Matsuoka, Takuma Izumi, Masafusa Onoue, Michael A. Strauss, Kazushi Iwasawa, Nobunari Kashikawa, Masayuki Akiyama, Kentaro Aoki, Junya Arita, Masatoshi Imanishi, Rikako Ishimoto, Toshihiro Kawaguchi, Kotaro Kohno, Chien-Hsiu Lee, Tohru Nagao, John D. Silverman, Yoshiki Toba

Abstract: We report the discovery of two quasars at a redshift of $z$ = 6.05, in the process of merging. They were serendipitously discovered from the deep multi-band imaging data collected by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The quasars, HSC $J$121503.42$-$014858.7 (C1) and HSC $J$121503.55$-$014859.3 (C2), both have luminous ($>$10$^{43}$ erg s$^{-1}$) Ly$α$ emission with a cle… ▽ More We report the discovery of two quasars at a redshift of $z$ = 6.05, in the process of merging. They were serendipitously discovered from the deep multi-band imaging data collected by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The quasars, HSC $J$121503.42$-$014858.7 (C1) and HSC $J$121503.55$-$014859.3 (C2), both have luminous ($>$10$^{43}$ erg s$^{-1}$) Ly$α$ emission with a clear broad component (full width at half maximum $>$1000 km s$^{-1}$). The rest-frame ultraviolet (UV) absolute magnitudes are $M_{1450} = -23.106 \pm 0.017$ (C1) and $-22.662 \pm 0.024$ (C2). Our crude estimates of the black hole masses provide $\log (M_{\rm BH}/M_\odot) = 8.1 \pm 0.3$ in both sources. The two quasars are separated by 12 kpc in projected proper distance, bridged by a structure in the rest-UV light suggesting that they are undergoing a merger. This pair is one of the most distant merging quasars reported to date, providing crucial insight into galaxy and black hole build-up in the hierarchical structure formation scenario. A companion paper will present the gas and dust properties captured by Atacama Large Millimeter/submillimeter Array observations, which provide additional evidence for and detailed measurements of the merger and also demonstrate that the two sources are not gravitationally-lensed images of a single quasar. △ Less

Submitted 3 May, 2024; originally announced May 2024.

Comments: Published from ApJL. A companion paper (T. Izumi et al.) presents the associated ALMA observations

Journal ref: The Astrophysical Journal Letters, Volume 965, Issue 1, id.L4, 8 pp (2024)

Cosmic Himalayas: The Highest Quasar Density Peak Identified in a 10,000 deg$^2$ Sky with Spatial Discrepancies between Galaxies, Quasars, and IGM HI

Authors: Yongming Liang, Masami Ouchi, Dongsheng Sun, Nobunari Kashikawa, Zheng Cai, Sebastiano Cantalupo, Kentaro Nagamine, Hidenobu Yajima, Takanobu Kirihara, Haibin Zhang, Mingyu Li, Rhythm Shimakawa, Xiaohui Fan, Kei Ito, Masayuki Tanaka, Yuichi Harikane, J. Xavier Prochaska, Andrea Travascio, Weichen Wang, Martin Elvis, Giuseppina Fabbiano, Junya Arita, Masafusa Onoue, John D. Silverman, Dongdong Shi , et al. (5 additional authors not shown)

Abstract: We report the identification of a quasar overdensity in the BOSSJ0210 field, dubbed Cosmic Himalayas, consisting of 11 quasars at $z=2.16-2.20$, the densest overdensity of quasars ($17σ$) in the $\sim$10,000 deg$^2$ of the Sloan Digital Sky Survey. We present the spatial distributions of galaxies and quasars and an HI absorption map of the intergalactic medium (IGM). On the map of 465 galaxies sel… ▽ More We report the identification of a quasar overdensity in the BOSSJ0210 field, dubbed Cosmic Himalayas, consisting of 11 quasars at $z=2.16-2.20$, the densest overdensity of quasars ($17σ$) in the $\sim$10,000 deg$^2$ of the Sloan Digital Sky Survey. We present the spatial distributions of galaxies and quasars and an HI absorption map of the intergalactic medium (IGM). On the map of 465 galaxies selected from the MAMMOTH-Subaru survey, we find two galaxy density peaks that do not fall on the quasar overdensity but instead exist at the northwest and southeast sides, approximately 25 $h^{-1}$ comoving-Mpc apart from the quasar overdensity. With a spatial resolution of 15 $h^{-1}$ comoving Mpc in projection, we produce a three-dimensional HI tomography map by the IGM Ly$α$ forest in the spectra of 23 SDSS/eBOSS quasars behind the quasar overdensity. Surprisingly, the quasar overdensity coincides with neither an absorption peak nor a transmission peak of IGM HI but lies near the border separating opaque and transparent volumes, with the more luminous quasars located in an environment with lesser IGM HI. Hence remarkably, the overdensity region traced by the 11 quasars, albeit all in coherently active states, has no clear coincidence with peaks of galaxies or HI absorption densities. Current physical scenarios with mixtures of HI overdensities and quasar photoionization cannot fully interpret the emergence of Cosmic Himalayas, suggesting this peculiar structure is an excellent laboratory to unveil the interplay between galaxies, quasars, and the IGM. △ Less

Submitted 24 April, 2024; originally announced April 2024.

Comments: 19 pages, 11 figures, submitted to ApJ, comments are welcome

Galaxy properties from the outskirts to the core of a protocluster at z=3.699

Authors: Jun Toshikawa, Stijn Wuyts, Nobunari Kashikawa, Hisakazu Uchiyama, Malcolm Bremer, Marcin Sawicki, Yoshiaki Ono, Mariko Kubo, Kei Ito

Abstract: We present follow-up spectroscopy on a protocluster candidate selected from the wide-field imaging of the Hyper SuprimeCam Subaru Strategic Programme. The target protocluster candidate was identified as a $4.5σ$ overdense region of $g$-dropout galaxies, and the redshifts of $g$-dropout galaxies are determined by detecting their Ly$α$ emission. Thirteen galaxies, at least, are found to be clusterin… ▽ More We present follow-up spectroscopy on a protocluster candidate selected from the wide-field imaging of the Hyper SuprimeCam Subaru Strategic Programme. The target protocluster candidate was identified as a $4.5σ$ overdense region of $g$-dropout galaxies, and the redshifts of $g$-dropout galaxies are determined by detecting their Ly$α$ emission. Thirteen galaxies, at least, are found to be clustering in the narrow redshift range of $Δz<0.05$ at $z=3.699$. This is clear evidence of the presence of a protocluster in the target region. Following the discovery of the protocluster at $z=3.699$, the physical properties and three-dimensional distribution of its member galaxies are investigated. Based on spectroscopically-confirmed $g$-dropout galaxies, we find an overabundance of rest-frame ultraviolet (UV) bright galaxies in the protocluster. The UV brightest protocluster member turns out to be an active galactic nucleus, and the other UV brighter members tend to show smaller Ly$α$ equivalent widths than field counterparts. The member galaxies tend to densely populate near the centre of the protocluster, but the separation from the nearest neighbour rather than the distance from the centre of the protocluster is more tightly correlated to galaxy properties, implying that the protocluster is still in an early phase of cluster formation and only close neighbours have a significant impact on the physical properties of protocluster members. The number density of massive galaxies, selected from an archival photometric-redshift catalogue, is higher near the centre of the protocluster, while dusty starburst galaxies are distributed on the outskirts. These observational results suggest that the protocluster consists of multiple galaxy populations, whose spatial distributions may hint at the developmental phase of the galaxy cluster. △ Less

Submitted 24 April, 2024; originally announced April 2024.

Comments: 11 pages, 8 figures, submitted to MNRAS

Black hole mass and Eddington ratio distributions of less-luminous quasars at $z\sim4$ in the Subaru Hyper Suprime-Cam Wide field

Authors: Wanqiu He, Masayuki Akiyama, Motohiro Enoki, Kohei Ichikawa, Kohei Inayoshi, Nobunari Kashikawa, Toshihiro Kawaguchi, Yoshiki Matsuoka, Tohru Nagao, Masafusa Onoue, Taira Oogi, Andreas Schulze, Yoshiki Toba, Yoshihiro Ueda

Abstract: We investigate the black hole mass function (BHMF) and Eddington ratio distribution function (ERDF) of broad-line AGNs at z=4, based on a sample of 52 quasars with i<23.2 at 3.50 < z < 4.25 from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) S16A-Wide2 dataset, and 1,462 quasars with i<20.2 in the same redshift range from the Sloan Digital Sky Survey (SDSS) DR7 quasar catalog. Virial BH… ▽ More We investigate the black hole mass function (BHMF) and Eddington ratio distribution function (ERDF) of broad-line AGNs at z=4, based on a sample of 52 quasars with i<23.2 at 3.50 < z < 4.25 from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) S16A-Wide2 dataset, and 1,462 quasars with i<20.2 in the same redshift range from the Sloan Digital Sky Survey (SDSS) DR7 quasar catalog. Virial BH masses of quasars are estimated using the width of the CIV 1549Å line and the continuum luminosity at 1350Å. To obtain the intrinsic broad-line AGN BHMF and ERDF, we correct for the incompleteness in the low-mass and/or low-Eddington-ratio ranges caused by the flux-limited selection. The resulting BHMF is constrained down to $\log M_{\rm BH}/M_{\odot}\sim7.5$. In comparison with broad-line AGN BHMFs at z=2 in the literature, we find that the number density of massive SMBHs peaks at higher redshifts, consistent with the "down-sizing" evolutionary scenario. Additionally, the resulting ERDF shows a negative dependence on BH mass, suggesting more massive SMBHs tend to accrete at lower Eddington ratios at z=4. With the derived intrinsic broad-line AGN BHMF, we also evaluate the active fraction of broad-line AGNs among the entire SMBH population at z=4. The resulting active fraction may suggest a positive dependence on BH mass. Finally, we examine the time evolution of broad-line AGN BHMF between z=4 and 6 through solving the continuity equation. The results suggest that the broad-line AGN BHMFs at z=4-6 only show evolution in their normalization, but with no significant changes in their shape. △ Less

Submitted 15 November, 2023; originally announced November 2023.

Comments: 38 pages, 26 figures. Submitted to ApJ. Data is available in https://github.com/wanqqq31/BHMF_hsc_z4

Subaru High-$z$ Exploration of Low-luminosity Quasars (SHELLQs). XVII. Black Hole Mass Distribution at $z \sim 6$ Estimated via Spectral Comparison with Low-$z$ Quasars

Authors: Ayumi Takahashi, Yoshiki Matsuoka, Masafusa Onoue, Michael A. Strauss, Nobunari Kashikawa, Yoshiki Toba, Kazushi Iwasawa, Masatoshi Imanishi, Masayuki Akiyama, Toshihiro Kawaguchi, Akatoki Noboriguchi, Chien-Hsiu Lee

Abstract: We report the distribution of black hole (BH) masses and Eddingont ratios estimated for a sample of 131 low luminosity quasars in the early cosmic epoch ($5.6 < z < 7.0$). Our work is based on Subaru High-$z$ Exploration of Low-Luminosity Quasars (SHELLQs) project, which has constructed a low luminosity quasar sample down to $M_{1450} \sim - 21$ mag, exploiting the survey data of Hyper Suprime-Cam… ▽ More We report the distribution of black hole (BH) masses and Eddingont ratios estimated for a sample of 131 low luminosity quasars in the early cosmic epoch ($5.6 < z < 7.0$). Our work is based on Subaru High-$z$ Exploration of Low-Luminosity Quasars (SHELLQs) project, which has constructed a low luminosity quasar sample down to $M_{1450} \sim - 21$ mag, exploiting the survey data of Hyper Suprime-Cam installed on Subaru Telescope. The discovery spectra of these quasars are limited to the rest-frame wavelengths of $\sim$ 1200 -- 1400 Å, which contains no emission lines that can be used as BH mass estimators. In order to overcome this problem, we made use of low-$z$ counterpart spectra from the Sloan Digital Sky Survey, which are spectrally matched to the high-$z$ spectra in overlapping wavelengths. We then combined the C~{\sc iv} emission line widths of the counterparts with the continuum luminosity from the SHELLQs data to estimate BH masses. The resulting BH mass distribution has a range of $\sim 10^{7-10} M_{\odot}$, with most of the quasars having BH masses $\sim 10^{8} M_{\odot}$ with sub-Eddington accretion. The present study provides not only a new insight into normal quasars in the reionization epoch, but also a new promising way to estimate BH masses of high-$z$ quasars without near-infrared spectroscopy. △ Less

Submitted 18 October, 2023; originally announced October 2023.

Comments: 18 pages, 17 figures, accepted for publication in ApJ

An enhanced abundance of bright galaxies in protocluster candidates at z~3-5

Authors: Jun Toshikawa, Stijn Wuyts, Nobunari Kashikawa, Chengze Liu, Marcin Sawicki, Roderik Overzier, Mariko Kubo, Hisakazu Uchiyama, Kei Ito, Malcolm Bremer, Yoshiaki Ono, Tadayuki Kodama, Yen-Ting Lin, Tomoki Saito

Abstract: We present a protocluster search covering $z\sim3$ to $z\sim5$ based on the combination of the Hyper SuprimeCam Subaru Strategic Programme and the CFHT Large Area $U$-band Deep Survey. We identify about 30 protocluster candidates per unit redshift over the $\sim25\,\mathrm{deg^2}$ area of the Deep/Ultra-Deep layer. Protocluster candidates are selected as regions with a significantly enhanced surfa… ▽ More We present a protocluster search covering $z\sim3$ to $z\sim5$ based on the combination of the Hyper SuprimeCam Subaru Strategic Programme and the CFHT Large Area $U$-band Deep Survey. We identify about 30 protocluster candidates per unit redshift over the $\sim25\,\mathrm{deg^2}$ area of the Deep/Ultra-Deep layer. Protocluster candidates are selected as regions with a significantly enhanced surface density of dropout galaxies. With this large sample, we characterise the properties of their individual member galaxies. We compare the number counts of dropout galaxies in protocluster candidates with that of coeval field galaxies. Rest-frame UV bright galaxies are over-abundant in protocluster candidates, a trend seen across the full redshift range studied. We do not find evidence for their spatial distribution within protocluster candidates to be distinct from their fainter counterparts, nor for their UV colour to be different from that of field galaxies with the same brightness. Cosmological simulations predict this bright-end excess, with the main cause being a richer population of massive galaxies, with only a minor contribution from an enhancement in star formation activity (and therefore UV emission) at fixed mass. $U$-to-$K$ SED modelling of our observed samples supports this interpretation. This environmental differentiation in number counts is already in place at $z\sim5$, with no significant redshift dependence over the range in lookback times probed. These observational results and model predictions suggest that the cosmic clock is ahead in high-density environments. △ Less

Submitted 12 October, 2023; originally announced October 2023.

Comments: 17 pages, 13 figures, accepted for publication in MNRAS

Subaru High-$z$ Exploration of Low-Luminosity Quasars (SHELLQs). XVIII. The Dark Matter Halo Mass of Quasars at $z\sim6$

Authors: Junya Arita, Nobunari Kashikawa, Yoshiki Matsuoka, Wanqiu He, Kei Ito, Yongming Liang, Rikako Ishimoto, Takehiro Yoshioka, Yoshihiro Takeda, Kazushi Iwasawa, Masafusa Onoue, Yoshiki Toba, Masatoshi Imanishi

Abstract: We present, for the first time, dark matter halo (DMH) mass measurement of quasars at $z\sim6$ based on a clustering analysis of 107 quasars. Spectroscopically identified quasars are homogeneously extracted from the HSC-SSP wide layer over $891\,\mathrm{deg^2}$. We evaluate the clustering strength by three different auto-correlation functions: projected correlation function, angular correlation fu… ▽ More We present, for the first time, dark matter halo (DMH) mass measurement of quasars at $z\sim6$ based on a clustering analysis of 107 quasars. Spectroscopically identified quasars are homogeneously extracted from the HSC-SSP wide layer over $891\,\mathrm{deg^2}$. We evaluate the clustering strength by three different auto-correlation functions: projected correlation function, angular correlation function, and redshift-space correlation function. The DMH mass of quasars at $z\sim6$ is evaluated as $5.0_{-4.0}^{+7.4}\times10^{12}\,h^{-1}M_\odot$ with the bias parameter $b=20.8\pm8.7$ by the projected correlation function. The other two estimators agree with these values, though each uncertainty is large. The DMH mass of quasars is found to be nearly constant $\sim10^{12.5}\,h^{-1}M_\odot$ throughout cosmic time, suggesting that there is a characteristic DMH mass where quasars are always activated. As a result, quasars appear in the most massive halos at $z \sim 6$, but in less extreme halos thereafter. The DMH mass does not appear to exceed the upper limit of $10^{13}\,h^{-1}M_\odot$, which suggests that most quasars reside in DMHs with $M_\mathrm{halo}<10^{13}\,h^{-1}M_\odot$ across most of the cosmic time. Our results supporting a significant increasing bias with redshift are consistent with the bias evolution model with inefficient AGN feedback at $z\sim6$. The duty cycle ($f_\mathrm{duty}$) is estimated as $0.019\pm0.008$ by assuming that DMHs in some mass interval can host a quasar. The average stellar mass is evaluated from stellar-to-halo mass ratio as $M_*=6.5_{-5.2}^{+9.6}\times10^{10}\,h^{-1}M_\odot$, which is found to be consistent with [C II] observational results. △ Less

Submitted 5 July, 2023; originally announced July 2023.

Comments: 22 pages, 8 figures, accepted for publication in ApJ

Quasar Luminosity Function at z = 7

Authors: Yoshiki Matsuoka, Masafusa Onoue, Kazushi Iwasawa, Michael A. Strauss, Nobunari Kashikawa, Takuma Izumi, Tohru Nagao, Masatoshi Imanishi, Masayuki Akiyama, John D. Silverman, Naoko Asami, James Bosch, Hisanori Furusawa, Tomotsugu Goto, James E. Gunn, Yuichi Harikane, Hiroyuki Ikeda, Kohei Inayoshi, Rikako Ishimoto, Toshihiro Kawaguchi, Satoshi Kikuta, Kotaro Kohno, Yutaka Komiyama, Chien-Hsiu Lee, Robert H. Lupton , et al. (19 additional authors not shown)

Abstract: We present the quasar luminosity function (LF) at $z = 7$, measured with 35 spectroscopically confirmed quasars at $6.55 < z < 7.15$. The sample of 22 quasars from the Subaru High-$z$ Exploration of Low-Luminosity Quasars (SHELLQs) project, combined with 13 brighter quasars in the literature, covers an unprecedentedly wide range of rest-frame ultraviolet magnitudes over $-28 < M_{1450} < -23$. We… ▽ More We present the quasar luminosity function (LF) at $z = 7$, measured with 35 spectroscopically confirmed quasars at $6.55 < z < 7.15$. The sample of 22 quasars from the Subaru High-$z$ Exploration of Low-Luminosity Quasars (SHELLQs) project, combined with 13 brighter quasars in the literature, covers an unprecedentedly wide range of rest-frame ultraviolet magnitudes over $-28 < M_{1450} < -23$. We found that the binned LF flattens significantly toward the faint end populated by the SHELLQs quasars. A maximum likelihood fit to a double power-law model has a break magnitude $M^*_{1450} = -25.60^{+0.40}_{-0.30}$, a characteristic density $Φ^* = 1.35^{+0.47}_{-0.30}$ Gpc$^{-3}$ mag$^{-1}$, and a bright-end slope $β= -3.34^{+0.49}_{-0.57}$, when the faint-end slope is fixed to $α= -1.2$ as observed at $z \le 6$. The overall LF shape remains remarkably similar from $z = 4$ to $7$, while the amplitude decreases substantially toward higher redshifts, with a clear indication of an accelerating decline at $z \ge 6$. The estimated ionizing photon density, $10^{48.2 \pm 0.1}$ s$^{-1}$ Mpc$^{-3}$, is less than 1 % of the critical rate to keep the intergalactic medium ionized at $z = 7$, and thus indicates that quasars are not a major contributor to cosmic reionization. △ Less

Submitted 18 May, 2023; originally announced May 2023.

Comments: The Astrophysical Journal Letters, in press

SILVERRUSH. XIII. A Catalog of 20,567 Ly$α$ Emitters at $z=2-7$ Identified in the Full-depth Data of the Subaru/HSC-SSP and CHORUS Surveys

Authors: Satoshi Kikuta, Masami Ouchi, Takatoshi Shibuya, Yongming Liang, Hiroya Umeda, Akinori Matsumoto, Kazuhiro Shimasaku, Yuichi Harikane, Yoshiaki Ono, Akio K. Inoue, Satoshi Yamanaka, Haruka Kusakabe, Rieko Momose, Nobunari Kashikawa, Yuichi Matsuda, Chien-Hsiu Lee

Abstract: We present 20,567 Ly$α$ emitters (LAEs) at $z=2.2-7.3$ that are photometrically identified by the SILVERRUSH program in a large survey area up to 25 deg$^2$ with deep images of five broadband filters (grizy) and seven narrowband filters targeting Ly$α$ lines at $z=2.2$, $3.3$, $4.9$, $5.7$, $6.6$, $7.0$, and $7.3$ taken by the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) and the Cosmic Hyd… ▽ More We present 20,567 Ly$α$ emitters (LAEs) at $z=2.2-7.3$ that are photometrically identified by the SILVERRUSH program in a large survey area up to 25 deg$^2$ with deep images of five broadband filters (grizy) and seven narrowband filters targeting Ly$α$ lines at $z=2.2$, $3.3$, $4.9$, $5.7$, $6.6$, $7.0$, and $7.3$ taken by the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) and the Cosmic HydrOgen Reionization Unveiled with Subaru (CHORUS) survey. We select secure $>5σ$ sources showing narrowband color excesses via Ly$α$ break screening, taking into account the spatial inhomogeneity of limiting magnitudes. After removing spurious sources by careful masking and visual inspection of coadded and multi-epoch images obtained over the 7 yr of the surveys, we construct LAE samples consisting of 6995, 4641, 726, 6124, 2058, 18, and 5 LAEs at $z=2.2$, 3.3, 4.9, 5.7, 6.6, 7.0, and 7.3, respectively, although the $z=7.3$ candidates are tentative. Our LAE catalogs contain 241 spectroscopically confirmed LAEs at the expected redshifts from previous work. We demonstrate that the number counts of our LAEs are consistent with previous studies with similar LAE selection criteria. The LAE catalogs will be made public on our project webpage with detailed descriptions of the content and ancillary information about the masks and limiting magnitudes. △ Less

Submitted 1 August, 2023; v1 submitted 15 May, 2023; originally announced May 2023.

Comments: 38 pages, 19 Figures, 5 Tables. Accepted for publication in ApJS

Updated measurements of [O III] 88 $μ$m, [C II] 158 $μ$m, and Dust Continuum Emission from a z=7.2 Galaxy

Authors: Yi W. Ren, Yoshinobu Fudamoto, Akio K. Inoue, Yuma Sugahara, Tsuyoshi Tokuoka, Yoichi Tamura, Hiroshi Matsuo, Kotaro Kohno, Hideki Umehata, Takuya Hashimoto, Rychard J. Bouwens, Renske Smit, Nobunari Kashikawa, Takashi Okamoto, Takatoshi Shibuya, Ikkoh Shimizu

Abstract: We present updated measurements of the [O III] 88 $μ$m, [C II] 158 $μ$m, and dust continuum emission from a star-forming galaxy at $z=7.212$, SXDF-NB1006-2, by utilizing Atacama Large Millimeter/submillimeter Array (ALMA) archival data sets analysed in previous studies and data sets that have not been analysed before. The follow-up ALMA observations with higher angular resolution and sensitivity r… ▽ More We present updated measurements of the [O III] 88 $μ$m, [C II] 158 $μ$m, and dust continuum emission from a star-forming galaxy at $z=7.212$, SXDF-NB1006-2, by utilizing Atacama Large Millimeter/submillimeter Array (ALMA) archival data sets analysed in previous studies and data sets that have not been analysed before. The follow-up ALMA observations with higher angular resolution and sensitivity reveal a clumpy structure of the [O III] emission on a scale of $0.32-0.85\,\rm{kpc}$. We also combined all the ALMA [O III] ([C II]) data sets and updated the [O III] ([C II]) detection to $5.9σ$ ($3.6σ-4.5σ$). The non-detection of [C II] with data from the REBELS large program implies the incompleteness of spectral-scan surveys using [C II] to detect galaxies with high star formation rates (SFRs) but marginal [C II] emission at high-$z$. The dust continuum at 90 $μ$m and 160 $μ$m remains undetected, indicating little dust content of $<3.9\times10^{6}\,M_\odot\,(3σ)$, and we obtained a more stringent constraint on the total infrared luminosity. We updated the [O III]/[C II] luminosity ratios to $10.2\pm4.7~(6.1\pm3.5$) and $20\pm12~(9.6\pm6.1$) for $4.5σ$ and $3.6σ$ [C II] detections, respectively, where the ratios in the parentheses are corrected for the surface brightness dimming effect on the extended [C II] emission. We also found a strong [C II] deficit ($0.6-1.3$ dex) between SXDF-NB1006-2 and the mean $L_{\rm{[CII]}}-\rm{SFR}$ relation of galaxies at $0<z<9$. △ Less

Submitted 9 March, 2023; v1 submitted 5 February, 2023; originally announced February 2023.

Comments: 13 pages, 10 figures, 8 tables, published in ApJ

COSMOS2020: Discovery of a protocluster of massive quiescent galaxies at $z=2.77$

Authors: Kei Ito, Masayuki Tanaka, Francesco Valentino, Sune Toft, Gabriel Brammer, Katriona M. L. Gould, Olivier Ilbert, Nobunari Kashikawa, Mariko Kubo, Yongming Liang, Henry J. McCracken, John R. Weaver

Abstract: Protoclusters of galaxies have been found in the last quarter century. However, most of them have been found through the overdensity of star-forming galaxies, and there had been no known structures identified by multiple spectroscopically confirmed quiescent galaxies at $z>2.5$. In this letter, we report the discovery of an overdense structure of massive quiescent galaxies with the spectroscopic r… ▽ More Protoclusters of galaxies have been found in the last quarter century. However, most of them have been found through the overdensity of star-forming galaxies, and there had been no known structures identified by multiple spectroscopically confirmed quiescent galaxies at $z>2.5$. In this letter, we report the discovery of an overdense structure of massive quiescent galaxies with the spectroscopic redshift $z=2.77$ in the COSMOS field, QO-1000. We first photometrically identify this structure as a $4.2σ$ overdensity with 14 quiescent galaxies in $7\times4\ {\rm pMpc^2}$ from the COSMOS2020 catalog. We then securely confirm the spectroscopic redshifts of 4 quiescent galaxies by detecting multiple Balmer absorption lines with Keck/MOSFIRE. All the spectroscopically confirmed members are massive ($\log{(M_\star/M_\odot)}>11.0$) and located in a narrow redshift range ($2.76<z<2.79$). Moreover, three of them are in the $1\times1\ {\rm pMpc^2}$ in the transverse direction at the same redshift ($z=2.760-2.763$). Such a concentration of four spectroscopically confirmed quiescent galaxies implies that QO-1000 is $>68$ times denser than in the general field. In addition, we confirm that they form a red sequence in the $J-K_s$ color. This structure's halo mass is estimated as $\log{(M_{\rm halo}/M_\odot)}>13.2$ from their stellar mass. Similar structures found in the IllustrisTNG simulation are expected to evolve into massive galaxy clusters with $\log{(M_{\rm halo}/M_\odot)}\geq14.8$ at $z=0$. These results suggest that QO-1000 is a more mature protocluster than the other known protoclusters. It is likely in a transition phase between the star-forming protoclusters and the quenched galaxy clusters. △ Less

Submitted 20 January, 2023; originally announced January 2023.

Comments: 11 pages, 4 figures, 1 table, accepted for publication in ApJL

MAMMOTH-Subaru V. Effects of Cosmic Variance on Ly$α$ Luminosity Functions at $z=2.2-2.3$

Authors: Ke Ma, Haibin Zhang, Zheng Cai, Yongming Liang, Nobunari Kashikawa, Mingyu Li, Yunjing Wu, Qiong Li, Xiaohui Fan, Sean D. Johnson, Masami Ouchi

Abstract: Cosmic variance introduces significant uncertainties into galaxy number density properties when surveying the high-z Universe with a small volume, such uncertainties produce the field-to-field variance of galaxy number $σ_{g}$ in observational astronomy. This uncertainty significantly affects the Luminosity Functions (LF) measurement of Lya Emitters (LAEs). For most previous Lya LF studies,… ▽ More Cosmic variance introduces significant uncertainties into galaxy number density properties when surveying the high-z Universe with a small volume, such uncertainties produce the field-to-field variance of galaxy number $σ_{g}$ in observational astronomy. This uncertainty significantly affects the Luminosity Functions (LF) measurement of Lya Emitters (LAEs). For most previous Lya LF studies, $σ_{g}$ is often adopted from predictions by cosmological simulations, but barely confirmed by observations. Measuring cosmic variance requires a huge sample over a large volume, exceeding the capabilities of most astronomical instruments. In this study, we demonstrate an observational approach for measuring the cosmic variance contribution for $z\approx2.2$ Lya LFs. The LAE candidates are observed using narrowband and broadband of the Subaru/Hyper Suprime-Cam (HSC), with 8 independent fields, making the total survey area $\simeq11.62$deg$^2$ and a comoving volume of $\simeq8.71\times10^6$Mpc$^3$. These eight fields are selected using the project of MAMMOTH. We report a best-fit Schechter function with parameters $α=-1.75$ (fixed), $L_{Lyα}^{*}=5.18_{-0.40}^{+0.43} \times 10^{42}$erg s$^{-1}$ and $φ_{Lya}^{*}=4.87_{-0.55}^{+0.54}\times10^{-4}$Mpc$^{-3}$ for the overall Lya LFs. After clipping out the regions that can bias the cosmic variance measurements, we calculate $σ_{g}$, by sampling LAEs within multiple pointings assigned on the field image. We investigate the relation between $σ_{g}$ and survey volume $V$, and fit a simple power law: $σ_g=k\times(\frac{V_{\rm eff}}{10^5 {\rm Mpc}^3})^β$. We find best-fit values of $-1.209_{-0.106}^{+0.106}$ for $β$ and $0.986_{-0.100}^{+0.108}$ for k. We compare our measurements with predictions from simulations and find that the cosmic variance of LAEs might be larger than that of general star-forming galaxies. △ Less

Submitted 18 January, 2023; originally announced January 2023.

Comments: 14 pages, 7 figures

MAMMOTH-Subaru IV. Large Scale Structure and Clustering Analysis of Ly$α$ Emitters and Ly$α$ Blobs at $z=2.2-2.3$

Authors: Haibin Zhang, Zheng Cai, Mingyu Li, Yongming Liang, Nobunari Kashikawa, Ke Ma, Yunjing Wu, Qiong Li, Sean D. Johnson, Satoshi Kikuta, Masami Ouchi, Xiaohui Fan

Abstract: We report the large scale structure and clustering analysis of Ly$α$ emitters (LAEs) and Ly$α$ blobs (LABs) at $z=2.2-2.3$. Using 3,341 LAEs, 117 LABs, and 58 bright (Ly$α$ luminosity $L_{\rm Lyα}>10^{43.4}$ erg s$^{-1}$) LABs at $z=2.2-2.3$ selected with Subaru/Hyper Suprime-Cam (HSC), we calculate the LAE overdensity to investigate the large scale structure at $z=2$. We show that 74% LABs and 78… ▽ More We report the large scale structure and clustering analysis of Ly$α$ emitters (LAEs) and Ly$α$ blobs (LABs) at $z=2.2-2.3$. Using 3,341 LAEs, 117 LABs, and 58 bright (Ly$α$ luminosity $L_{\rm Lyα}>10^{43.4}$ erg s$^{-1}$) LABs at $z=2.2-2.3$ selected with Subaru/Hyper Suprime-Cam (HSC), we calculate the LAE overdensity to investigate the large scale structure at $z=2$. We show that 74% LABs and 78% bright LABs locate in overdense regions, which is consistent with the trend found by previous studies that LABs generally locate in overdense regions. We find that one of our 8 fields dubbed J1349 contains $39/117\approx33\%$ of our LABs and $22/58\approx38\%$ of our bright LABs. A unique and overdense $24'\times12'$ ($\approx 40\times20$ comoving Mpc$^2$) region in J1349 has 12 LABs (8 bright LABs). By comparing to SSA22 that is one of the most overdense LAB regions found by previous studies, we show that the J1349 overdense region contains $\geq 2$ times more bright LABs than the SSA22 overdense region. We calculate the angular correlation functions (ACFs) of LAEs and LABs in the unique J1349 field and fit the ACFs with a power-law function to measure the slopes. The slopes of LAEs and LABs are similar, while the bright LABs show a $\approx 2$ times larger slope suggesting that bright LABs are more clustered than faint LABs and LAEs. We show that the amplitudes of ACFs of LABs are higher than LAEs, which suggests that LABs have a $\approx 10$ times larger galaxy bias and field-to-field variance than LAEs. The strong field-to-field variance is consistent with the large differences of LAB numbers in our 8 fields. △ Less

Submitted 18 January, 2023; originally announced January 2023.

Comments: 8 pages, 5 figures, and 1 table

MAMMOTH-Subaru III. Ly$α$ Halo Extended to $\sim200$ kpc Identified by Stacking $\sim 3300$ Ly$α$ Emitters at $z=2.2-2.3$

Authors: Haibin Zhang, Zheng Cai, Yongming Liang, Ke Ma, Nobunari Kashikawa, Mingyu Li, Yunjing Wu, Qiong Li, Sean D. Johnson, Masami Ouchi, Xiaohui Fan

Abstract: In this paper, we present a Ly$α$ halo extended to $\sim200$ kpc identified by stacking $\sim 3300$ Ly$α$ emitters at $z=2.2-2.3$. We carry out imaging observations and data reduction with Subaru/Hyper Suprime-Cam (HSC). Our total survey area is $\sim12$ deg$^2$ and imaging depths are $25.5-27.0$ mag. Using the imaging data, we select 1,240 and 2,101 LAE candidates at $z=2.2$ and 2.3, respectively… ▽ More In this paper, we present a Ly$α$ halo extended to $\sim200$ kpc identified by stacking $\sim 3300$ Ly$α$ emitters at $z=2.2-2.3$. We carry out imaging observations and data reduction with Subaru/Hyper Suprime-Cam (HSC). Our total survey area is $\sim12$ deg$^2$ and imaging depths are $25.5-27.0$ mag. Using the imaging data, we select 1,240 and 2,101 LAE candidates at $z=2.2$ and 2.3, respectively. We carry out spectroscopic observations of our LAE candidates and data reduction with Magellan/IMACS to estimate the contamination rate of our LAE candidates. We find that the contamination rate of our sample is low (8%). We stack our LAE candidates with a median stacking method to identify the Ly$α$ halo at $z=2$. We show that the Ly$α$ halo is extended to $\sim200$ kpc at a surface brightness level of $10^{-20}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$. Comparing to previous studies, our Ly$α$ halo is more extended at radii of $\sim25-100$ kpc, which is not likely caused by the contamination in our sample but by different redshifts and fields instead. To investigate how central galaxies affect surrounding LAHs, we divide our LAEs into subsamples based on the Ly$α$ luminosity ($L_{\rm Lyα}$), rest-frame Ly$α$ equivalent width (EW$_0$), and UV magnitude (M$_{\rm uv}$). We stack the subsamples and find that higher $L_{\rm Lyα}$, lower EW$_0$, and brighter M$_{\rm uv}$ cause more extended halos. Our results suggest that more massive LAEs generally have more extended Ly$α$ halos. △ Less

Submitted 18 January, 2023; originally announced January 2023.

Comments: 9 pages, 7 figures, and 2 tables

Detection of stellar light from quasar host galaxies at redshifts above 6

Authors: Xuheng Ding, Masafusa Onoue, John D. Silverman, Yoshiki Matsuoka, Takuma Izumi, Michael A. Strauss, Knud Jahnke, Camryn L. Phillips, Junyao Li, Marta Volonteri, Zoltan Haiman, Irham Taufik Andika, Kentaro Aoki, Shunsuke Baba, Rebekka Bieri, Sarah E. I. Bosman, Connor Bottrell, Anna-Christina Eilers, Seiji Fujimoto, Melanie Habouzit, Masatoshi Imanishi, Kohei Inayoshi, Kazushi Iwasawa, Nobunari Kashikawa, Toshihiro Kawaguchi , et al. (19 additional authors not shown)

Abstract: The detection of starlight from the host galaxies of quasars during the reionization epoch ($z>6$) has been elusive, even with deep HST observations. The current highest redshift quasar host detected, at $z=4.5$, required the magnifying effect of a foreground lensing galaxy. Low-luminosity quasars from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) mitigate the challenge of detecting the… ▽ More The detection of starlight from the host galaxies of quasars during the reionization epoch ($z>6$) has been elusive, even with deep HST observations. The current highest redshift quasar host detected, at $z=4.5$, required the magnifying effect of a foreground lensing galaxy. Low-luminosity quasars from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) mitigate the challenge of detecting their underlying, previously-undetected host galaxies. Here we report rest-frame optical images and spectroscopy of two HSC-SSP quasars at $z>6$ with JWST. Using NIRCam imaging at 3.6$μ$m and 1.5$μ$m and subtracting the light from the unresolved quasars, we find that the host galaxies are massive (stellar masses of $13\times$ and $3.4\times$ $10^{10}$ M$_{\odot}$, respectively), compact, and disk-like. NIRSpec medium-resolution spectroscopy shows stellar absorption lines in the more massive quasar, confirming the detection of the host. Velocity-broadened gas in the vicinity of these quasars enables measurements of their black hole masses ($1.4\times 10^9$ and $2.0\times$ $10^{8}$ M$_{\odot}$, respectively). Their location in the black hole mass - stellar mass plane is consistent with the distribution at low redshift, suggesting that the relation between black holes and their host galaxies was already in place less than a billion years after the Big Bang. △ Less

Submitted 23 June, 2023; v1 submitted 25 November, 2022; originally announced November 2022.

Comments: Matched to the published Nature version of the article. 27 pages, 4 main figures, 1 table, 6 supplementary figures, 2 supplementary table

X-ray emission from a rapidly accreting narrow-line Seyfert 1 galaxy at z=6.56

Authors: Julien Wolf, Kirpal Nandra, Mara Salvato, Johannes Buchner, Masafusa Onoue, Teng Liu, Riccardo Arcodia, Andrea Merloni, Stefano Ciroi, Francesco Di Mille, Vadim Burwitz, Marcella Brusa, Rikako Ishimoto, Nobunari Kashikawa, Yoshiki Matsuoka, Tanya Urrutia, Sophia Waddell

Abstract: This study aims at identifying luminous quasars at $z>5.7$ among X-ray-selected sources in the eROSITA Final Equatorial-Depth Survey (eFEDS) in order to place a lower limit on black hole accretion well into the epoch of re-ionisation. We confirm the low significance detection with eROSITA of a previously known, optically faint $z=6.56$ quasar from the Subaru High-z Exploration of Low-luminosity Qu… ▽ More This study aims at identifying luminous quasars at $z>5.7$ among X-ray-selected sources in the eROSITA Final Equatorial-Depth Survey (eFEDS) in order to place a lower limit on black hole accretion well into the epoch of re-ionisation. We confirm the low significance detection with eROSITA of a previously known, optically faint $z=6.56$ quasar from the Subaru High-z Exploration of Low-luminosity Quasars (SHELLQs) survey. We obtained a pointed follow-up observation of the source with the Chandra X-ray telescope in order to confirm the eROSITA detection. Using new near-infrared spectroscopy, we derived the physical properties of the super-massive black hole. Finally, we used this detection to infer a lower limit on the black hole accretion density rate at $z>6$. The Chandra observation confirms the eFEDS source as the most distant blind X-ray detection to date. The derived X-ray luminosity is high with respect to the rest-frame optical emission of the quasar. With a narrow MgII line, low derived black hole mass, and high Eddington ratio, as well as its steep photon index, the source shows properties that are similar to local narrow-line Seyfert 1 galaxies, which are thought to be powered by young super-massive black holes. In combination with a previous high-redshift quasar detection in the field, we show that quasars with $L_{2-10 \, \mathrm{keV}} >10^{45} \, \mathrm{erg \, s^{-1}}$ dominate accretion onto super-massive black holes at $z\sim 6$. △ Less

Submitted 24 November, 2022; originally announced November 2022.

Comments: Accepted for publication in A&A

Journal ref: A&A 669, A127 (2023)

The physical origin for spatially large scatter of IGM opacity at the end of reionization: the IGM Ly$α$ opacity-galaxy density relation

Authors: Rikako Ishimoto, Nobunari Kashikawa, Daichi Kashino, Kei Ito, Yongming Liang, Zheng Cai, Takehiro Yoshioka, Katsuya Okoshi, Toru Misawa, Masafusa Onoue, Yoshihiro Takeda, Hisakazu Uchiyama

Abstract: The large opacity fluctuations in the $z > 5.5$ Ly$α$ forest may indicate inhomogeneous progress of reionization. To explain the observed large scatter of the effective Ly$α$ optical depth ($τ_{\rm eff}$) of the intergalactic medium (IGM), fluctuation of UV background ($Γ$ model) or the IGM gas temperature ($T$ model) have been proposed, which predict opposite correlations between $τ_{\rm eff}$ an… ▽ More The large opacity fluctuations in the $z > 5.5$ Ly$α$ forest may indicate inhomogeneous progress of reionization. To explain the observed large scatter of the effective Ly$α$ optical depth ($τ_{\rm eff}$) of the intergalactic medium (IGM), fluctuation of UV background ($Γ$ model) or the IGM gas temperature ($T$ model) have been proposed, which predict opposite correlations between $τ_{\rm eff}$ and galaxy density. In order to address which model can explain the large scatter of $τ_{\rm eff}$, we search for Ly$α$ emitters (LAEs) around two (J1137+3549 and J1602+4228) quasar sightlines with $τ_{\rm eff}\sim3$ and J1630+4012 sightline with $τ_{\rm eff}\sim5.5$. Using a narrowband imaging with Subaru/Hyper Suprime-Cam, we draw LAE density maps to explore their spatial distributions. Overdensities are found within 20 $h^{-1}$Mpc of the quasar sightlines in the low $τ_{\rm eff}$ regions, while a deficit of LAEs is found in the high $τ_{\rm eff}$ region. Although the $τ_{\rm eff}$ of the three quasar sightlines are neither high nor low enough to clearly distinguish the two models, these observed $τ_{\rm eff}$-galaxy density relations all consistently support the $Γ$ model rather than the $T$ model in the three fields, along with the previous studies. The observed overdensities near the low $τ_{\rm eff}$ sightlines may suggest that the relic temperature fluctuation does not affect reionization that much. Otherwise, these overdensities could be attributed to other factors besides the reionization process, such as the nature of LAEs as poor tracers of underlying large-scale structures. △ Less

Submitted 11 July, 2022; originally announced July 2022.

Comments: 13 pages, 14 figures, accepted for publication in MNRAS

Exploratory Study of Transverse Proximity Effect around BAL Quasars

Authors: Toru Misawa, Rikako Ishimoto, Satoshi Kobu, Nobunari Kashikawa, Katsuya Okoshi, Akatoki Noboriguchi, Malte Schramm, Qiang Liu

Abstract: We aim to find out the reason why there exists an anisotropic HI absorption around quasars; i.e., the environments around quasars are highly biased toward producing strong HI absorption in transverse direction while there exists a significant deficit of HI absorption within a few Mpc of quasars along line-of-sight. The most plausible explanation for this opposite trend is that the transverse direc… ▽ More We aim to find out the reason why there exists an anisotropic HI absorption around quasars; i.e., the environments around quasars are highly biased toward producing strong HI absorption in transverse direction while there exists a significant deficit of HI absorption within a few Mpc of quasars along line-of-sight. The most plausible explanation for this opposite trend is that the transverse direction is shadowed from the quasar UV radiation due to dust torus. However, a critical weakness of this idea is that we have no information on inclination angle of our sightline relative to the torus. In this study, we examine environments of quasars with broad absorption troughs in their spectra (i.e., BAL quasars) because it is widely believed that BAL troughs are observed if the central continuum is viewed from the side through their powerful outflows near the dust torus. With closely separated 12 projected quasar pairs at different redshift with separation angle of $θ$$<$120$^{\prime\prime}$, we examine HI absorption at foreground BAL quasars in spectra of background quasars. We confirm there exist optically thick gas around two of 12 BAL quasars, and that the mean HI absorption strength is EW$_{\rm rest}$$\sim$1A. These are consistent to the past results around non-BAL quasars, although not statistically significant. However, the origins of optically thick HI absorbers around BAL and non-BAL quasars could be different since their column densities are different by $\sim$3 orders of magnitude. The larger sample would be required for narrowing down possible scenarios for the anisotropic HI absorption around quasars. △ Less

Submitted 22 June, 2022; originally announced June 2022.

Comments: 17 pages, 10 figures, accepted for publication in ApJ

CHORUS IV: Mapping the Spatially Inhomogeneous Cosmic Reionization with Subaru HSC

Authors: Takehiro Yoshioka, Nobunari Kashikawa, Akio K. Inoue, Satoshi Yamanaka, Kazuhiro Shimasaku, Yuichi Harikane, Takatoshi Shibuya, Rieko Momose, Kei Ito, Yongming Liang, Rikako Ishimoto, Yoshihiro Takeda, Masami Ouchi, Chien-Hsiu Lee

Abstract: The spatial inhomogeneity is one of the important features for understanding the reionization process; however, it has not yet been fully quantified. To map this inhomogeneous distribution, we simultaneously detect Ly$α$ emitters (LAEs) and Lyman break galaxies (LBGs) at $z \sim 6.6$ from the Subaru/Hyper Suprime-Cam (HSC) large-area ($\sim1.5\,\mathrm{ deg}^2 = 34000\,\mathrm{cMpc}^2$) deep surve… ▽ More The spatial inhomogeneity is one of the important features for understanding the reionization process; however, it has not yet been fully quantified. To map this inhomogeneous distribution, we simultaneously detect Ly$α$ emitters (LAEs) and Lyman break galaxies (LBGs) at $z \sim 6.6$ from the Subaru/Hyper Suprime-Cam (HSC) large-area ($\sim1.5\,\mathrm{ deg}^2 = 34000\,\mathrm{cMpc}^2$) deep survey. We estimate the neutral fraction, $x_\mathrm{HI}$, from the observed number density ratio of LAEs to LBGs, $n(\mathrm{LAE})/n(\mathrm{LBG})$ based on numerical radiative transfer simulation, in which model galaxies are selected to satisfy the observed selection function. While the average $x_\mathrm{HI}$ within the field of view is found to be $x_\mathrm{HI} < 0.4$, which is consistent with previous studies, the variation of $n(\mathrm{LAE})/n(\mathrm{LBG})$ within the field of view for each $140\,\mathrm{pMpc}^2$ is found to be as large as a factor of three. This may suggest a spatially inhomogeneous topology of reionization, but it also leaves open the possibility that the variation is based on the inherent large-scale structure of the galaxy distribution. Based on the simulations, it may be difficult to distinguish between the two from the current survey. We also find that LAEs in the high LAE density region are more populate high $\mathrm{EW}_0$, supporting that the observed $n(\mathrm{LAE})/n(\mathrm{LBG})$ is more or less driven by the neutral fraction, though the statistical significance is not high. △ Less

Submitted 18 January, 2022; originally announced January 2022.

Comments: 20 pages, 22 figures, accepted for publication in ApJ

Reply to: Three papers regarding the origin of GN-z11-flash

Authors: Linhua Jiang, Shu Wang, Bing Zhang, Nobunari Kashikawa, Luis C. Ho, Zheng Cai, Eiichi Egami, Gregory Walth, Yi-Si Yang, Bin-Bin Zhang, Hai-Bin Zhao

Abstract: In Jiang et al., we detected a bright flash (hereafter GN-z11-flash) that appeared as compact continuum emission during our Keck MOSFIRE observations of the galaxy GN-z11 at z ~ 11. We performed a comprehensive analysis of the origin of the flash using all available information and our current understanding of known man-made objects or moving objects in the solar system. We found that GN-z11-flash… ▽ More In Jiang et al., we detected a bright flash (hereafter GN-z11-flash) that appeared as compact continuum emission during our Keck MOSFIRE observations of the galaxy GN-z11 at z ~ 11. We performed a comprehensive analysis of the origin of the flash using all available information and our current understanding of known man-made objects or moving objects in the solar system. We found that GN-z11-flash was likely a rest-frame UV flash associated with a long gamma-ray burst (GRB) from GN-z11. Recently, Steinhardt et al., Michalowski et al., and Nir et al. reported that GN-z11-flash was more likely from a satellite. While one cannot completely rule out the possibility of unknown satellites (or debris), we find that either the chance probabilities of being a satellite estimated by these authors have been largely overestimated or their identified satellites have been ruled out in our original analysis. Our new calculations show that the probability of GN-z11-flash being a satellite is still lower than that of it being a signal originated from GN-z11. △ Less

Submitted 14 December, 2021; originally announced December 2021.

Comments: Published in Nature Astronomy Matters Arising 5, 998-1000, 2021; Authors' version; Some content slightly overlaps with our previous reply in arXiv:2102.01239

A Wide and Deep Exploration of Radio Galaxies with Subaru HSC (WERGS). VI. Distant Filamentary Structures Pointed by High-z Radio Galaxies at z~4

Authors: Hisakazu Uchiyama, Takuji Yamashita, Jun Toshikawa, Nobunari Kashikawa, Kohei Ichikawa, Mariko Kubo, Kei Ito, Nozomu Kawakatu, Tohru Nagao, Yoshiki Toba, Yoshiaki Ono, Yuichi Harikane, Masatoshi Imanishi, Masaru Kajisawa, Chien-Hsiu Lee, Yongming Liang

Abstract: We present the environmental properties around high-$z$ radio galaxies (HzRGs) at $z\sim4$, which have been poorly investigated because of their rarity. We use the largest samples of HzRGs and $g$-dropout galaxy overdense regions at $z\sim4$, which were constructed from Hyper Suprime-Cam Subaru Strategic Program, to characterize the HzRG environments statistically. We measure the $g$-dropout galax… ▽ More We present the environmental properties around high-$z$ radio galaxies (HzRGs) at $z\sim4$, which have been poorly investigated because of their rarity. We use the largest samples of HzRGs and $g$-dropout galaxy overdense regions at $z\sim4$, which were constructed from Hyper Suprime-Cam Subaru Strategic Program, to characterize the HzRG environments statistically. We measure the $g$-dropout galaxy overdensities around 21 HzRGs whose rest-frame 1.4 GHz radio luminosities ($L_{1.4\mathrm{GHz}}$) are $10^{26-27}$ W Hz$^{-1}$. We find that the overdensities around the faint HzRGs with $L_{1.4\mathrm{GHz}}\sim10^{26.0-26.5}$ W Hz$^{-1}$ tend to be higher than that of the $g$-dropout galaxies. On the other hand, no significant difference of density environments is found between the luminous HzRGs with $ L_{1.4\mathrm{GHz}}\sim10^{26.5-27.0} $ W Hz$^{-1}$ and the $g$-dropout galaxies. The HzRGs are found to occupy more massive halos than $g$-dropout galaxies through a cross-correlation between the HzRGs and $g$-dropout galaxies. This trend is more pronounced in the faint HzRGs. These results are consistent with a scenario where HzRGs get older and more massive as the radio-luminosity decreases. The HzRGs are expected to trace the progenitors of local cluster halos from their calculated halo mass. In addition, we find that surrounding galaxies tend to distribute along the radio-jet major axis of the HzRGs at angular distances less than $\lesssim500$ physical kpc. Our findings imply the onset of the filamentary structures around the HzRGs at $z\sim4$. △ Less

Submitted 2 December, 2021; originally announced December 2021.

Comments: 22 pages, 8 figures, accepted in ApJ

Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs). XVI. 69 New Quasars at 5.8 < z < 7.0

Authors: Yoshiki Matsuoka, Kazushi Iwasawa, Masafusa Onoue, Takuma Izumi, Nobunari Kashikawa, Michael A. Strauss, Masatoshi Imanishi, Tohru Nagao, Masayuki Akiyama, John D. Silverman, Naoko Asami, James Bosch, Hisanori Furusawa, Tomotsugu Goto, James E. Gunn, Yuichi Harikane, Hiroyuki Ikeda, Rikako Ishimoto, Toshihiro Kawaguchi, Nanako Kato, Satoshi Kikuta, Kotaro Kohno, Yutaka Komiyama, Chien-Hsiu Lee, Robert H. Lupton , et al. (19 additional authors not shown)

Abstract: We present the spectroscopic discovery of 69 quasars at 5.8 < z < 7.0, drawn from the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP) imaging survey data. This is the 16th publication from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, and completes identification of all but the faintest candidates (i.e., i-band dropouts with zAB < 24 and y-band detections, and z… ▽ More We present the spectroscopic discovery of 69 quasars at 5.8 < z < 7.0, drawn from the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP) imaging survey data. This is the 16th publication from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, and completes identification of all but the faintest candidates (i.e., i-band dropouts with zAB < 24 and y-band detections, and z-band dropouts with yAB < 24) with Bayesian quasar probability Pq > 0.1 in the HSC-SSP third public data release (PDR3). The sample reported here also includes three quasars with Pq < 0.1 at z ~ 6.6, which we selected in an effort to completely cover the reddest point sources with simple color cuts. The number of high-z quasars discovered in SHELLQs has now grown to 162, including 23 type-II quasar candidates. This paper also presents identification of seven galaxies at 5.6 < z < 6.7, an [O III] emitter at z = 0.954, and 31 Galactic cool stars and brown dwarfs. High-z quasars and galaxies comprise 75 % and 16 % respectively of all the spectroscopic SHELLQs objects that pass our latest selection algorithm with the PDR3 photometry. That is, a total of 91 % of the objects lie at z > 5.6. This demonstrates that the algorithm has very high efficiency, even though we are probing an unprecedentedly low-luminosity population down to M1450 ~ -21 mag. △ Less

Submitted 24 November, 2021; originally announced November 2021.

Comments: Accepted for publication in The Astrophysical Journal Supplement Series

The SSA22 HI Tomography Survey (SSA22-HIT). I. Data Set and Compiled Redshift Catalog

Authors: Ken Mawatari, Akio K. Inoue, Toru Yamada, Tomoki Hayashino, J. Xavier Prochaska, Khee-Gan Lee, Nicolas Tejos, Nobunari Kashikawa, Takuya Otsuka, Satoshi Yamanaka, David J. Schlegel, Yuichi Matsuda, Joseph F. Hennawi, Ikuru Iwata, Hideki Umehata, Shiro Mukae, Masami Ouchi, Yuma Sugahara, Yoichi Tamura

Abstract: We conducted a deep spectroscopic survey, named SSA22-HIT, in the SSA22 field with the DEep Imaging MultiObject Spectrograph (DEIMOS) on the Keck telescope, designed to tomographically map high-z HI gas through analysis of Lya absorption in background galaxies' spectra. In total, 198 galaxies were spectroscopically confirmed at 2.5 < z < 6 with a few low-z exceptions in the 26 x 15 arcmin^2 area,… ▽ More We conducted a deep spectroscopic survey, named SSA22-HIT, in the SSA22 field with the DEep Imaging MultiObject Spectrograph (DEIMOS) on the Keck telescope, designed to tomographically map high-z HI gas through analysis of Lya absorption in background galaxies' spectra. In total, 198 galaxies were spectroscopically confirmed at 2.5 < z < 6 with a few low-z exceptions in the 26 x 15 arcmin^2 area, of which 148 were newly determined in this study. Our redshift measurements were merged with previously confirmed redshifts available in the 34 x 27 arcmin^2 area of the SSA22 field. This compiled catalog containing 730 galaxies of various types at z > 2 is useful for various applications, and it is made publicly available. Our SSA22-HIT survey has increased by approximately twice the number of spectroscopic redshifts of sources at z > 3.2 in the observed field. From a comparison with publicly available redshift catalogs, we show that our compiled redshift catalog in the SSA22 field is comparable to those among major extragalactic survey fields in terms of a combination of wide area and high surface number density of objects at z > 2. About 40 % of the spectroscopically confirmed objects in SSA22-HIT show reasonable quality of spectra in the wavelengths shorter than Lya when a sufficient amount of smoothing is adopted. Our data set enables us to make the HI tomographic map at z > 3, which we present in a parallel study. △ Less

Submitted 9 March, 2023; v1 submitted 22 October, 2021; originally announced October 2021.

Comments: Accepted for publication in Astronomical Journal, 29 pages, 17 figures, 7 tables

Where's Swimmy?: Mining unique color features buried in galaxies by deep anomaly detection using Subaru Hyper Suprime-Cam data

Authors: Takumi S. Tanaka, Rhythm Shimakawa, Kazuhiro Shimasaku, Yoshiki Toba, Nobunari Kashikawa, Masayuki Tanaka, Akio K. Inoue

Abstract: We present the Swimmy (Subaru WIde-field Machine-learning anoMalY) survey program, a deep-learning-based search for unique sources using multicolored ($grizy$) imaging data from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). This program aims to detect unexpected, novel, and rare populations and phenomena, by utilizing the deep imaging data acquired from the wide-field coverage of the H… ▽ More We present the Swimmy (Subaru WIde-field Machine-learning anoMalY) survey program, a deep-learning-based search for unique sources using multicolored ($grizy$) imaging data from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). This program aims to detect unexpected, novel, and rare populations and phenomena, by utilizing the deep imaging data acquired from the wide-field coverage of the HSC-SSP. This article, as the first paper in the Swimmy series, describes an anomaly detection technique to select unique populations as "outliers" from the data-set. The model was tested with known extreme emission-line galaxies (XELGs) and quasars, which consequently confirmed that the proposed method successfully selected 60-70% of the quasars and 60% of the XELGs without labeled training data. In reference to the spectral information of local galaxies at $z=$0.05-0.2 obtained from the Sloan Digital Sky Survey, we investigated the physical properties of the selected anomalies and compared them based on the significance of their outlier values. The results revealed that XELGs constitute notable fractions of the most anomalous galaxies, and certain galaxies manifest unique morphological features. In summary, a deep anomaly detection is an effective tool that can search rare objects, and ultimately, unknown unknowns with large data-sets. Further development of the proposed model and selection process can promote the practical applications required to achieve specific scientific goals. △ Less

Submitted 15 February, 2022; v1 submitted 11 October, 2021; originally announced October 2021.

Comments: 24 pages, 19 figures. Accepted for publication in Publications of the Astronomical Society of Japan (PASJ)

Journal ref: PASJ 74 (2022) 1-23

Multiple MgII Absorption Systems in the Lines of Sight to Quadruply Lensed Quasar H1413+1143

Authors: Katsuya Okoshi, Yosuke Minowa, Nobunari Kashikawa, Toru Misawa, Daichi Kashino, Hajime Sugai, Kazuya Matsubayashi, Atsushi Shimono, Shinobu Ozaki

Abstract: We find multiple Mg II absorption systems at redshift z=1.66, 2.069, and 2.097 in the spatially resolved spectra of the quadruply gravitationally lensed quasar H1413+1143 utilizing the Kyoto tridimensional spectrograph II (Kyoto 3DII) spectrograph on board the Subaru telescope. Here we present the first measurement of differences in Mg II absorption strength of the multiple intervening absorbers,… ▽ More We find multiple Mg II absorption systems at redshift z=1.66, 2.069, and 2.097 in the spatially resolved spectra of the quadruply gravitationally lensed quasar H1413+1143 utilizing the Kyoto tridimensional spectrograph II (Kyoto 3DII) spectrograph on board the Subaru telescope. Here we present the first measurement of differences in Mg II absorption strength of the multiple intervening absorbers, which include ones identified as damped Lyman alpha (DLA) absorption systems. Our detection of the significant Mg II absorptions in the spatially resolved spectra reveals the inhomogeneous chemical enrichment on scales of about 12 kpc within the separation of the four sightlines. For the DLA system at z=1.66, the rest equivalent widths of the Mg II absorption lines between the four spatially resolved lines of sight change by factors of up to about 6, which trace the variations in the H I absorption strength. This suggests that inhomogeneous cold absorbers that give rise to the strong H I/Mg II absorptions dwell on a scale of about 6-12 kpc between the four lines of sight. We also investigate the degree of variation in the equivalent width of the absorption lines between the lines of sight. We find that the systems giving rise to strong absorptions in the spectra of the quadruply lensed quasars tend to have a high degree of variation in absorption strength between the lines of sight toward the lensed quasars. △ Less

Submitted 4 October, 2021; originally announced October 2021.

Comments: 21 pages, 13 figures, and 3 tables Accepted for publication in the Astronomical Journal

Ionizing radiation from AGNs at z>3.3 with the Subaru Hyper Suprime-Cam Survey and the CFHT Large Area U-band Deep Survey (CLAUDS)

Authors: Ikuru Iwata, Marcin Sawicki, Akio K. Inoue, Masayuki Akiyama, Genoveva Micheva, Toshihiro Kawaguchi, Nobunari Kashikawa, Stephen Gwyn, Stephane Arnouts, Jean Coupon, Guillaume Desprez

Abstract: We use deep and wide imaging data from the CFHT Large Area U-band Deep Survey (CLAUDS) and the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) to constrain the ionizing radiation (Lyman Continuum; LyC) escape fraction from AGNs at $z \sim 3 - 4$. For 94 AGNs with spectroscopic redshifts at $3.3 < z < 4.0$, we use their U-band / i-band flux ratios to estimate LyC transmission of individual AGN… ▽ More We use deep and wide imaging data from the CFHT Large Area U-band Deep Survey (CLAUDS) and the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) to constrain the ionizing radiation (Lyman Continuum; LyC) escape fraction from AGNs at $z \sim 3 - 4$. For 94 AGNs with spectroscopic redshifts at $3.3 < z < 4.0$, we use their U-band / i-band flux ratios to estimate LyC transmission of individual AGNs. The distribution of their LyC transmission shows values lower than the range of LyC transmission values for IGM of the same redshift range, which suggests that LyC escape fraction of AGNs at $z>3.3$ is considerably lower than unity in most cases. We do not find any trend in LyC transmission values depending on their UV luminosities. Based on the photometry of stacked images we find the average flux ratio of LyC and non-ionizing UV photons escaping from the objects $(f_{LyC}/f_{UV})^{out} = 0.182 \pm 0.043$ for AGNs at $3.3<z<3.6$, which corresponds to LyC escape fraction $f_{esc} = 0.303 \pm 0.072$ if we assume a fiducial intrinsic SED of AGN. Based on the estimated LyC escape fraction and the UV luminosity function of AGNs, we argue that UV-selected AGNs' contribution to the LyC emissivity at the epoch is minor, although the size of their contribution largely depends on the shape of the UV luminosity function. △ Less

Submitted 21 September, 2021; originally announced September 2021.

Comments: 19 pages, 12 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society

SILVERRUSH. XI. Intensity Mapping for Lya Emission Extending over $100-1000$ comoving kpc around $z\sim2-7$ LAEs with Subaru HSC-SSP and CHORUS Data

Authors: Shotaro Kikuchihara, Yuichi Harikane, Masami Ouchi, Yoshiaki Ono, Takatoshi Shibuya, Ryohei Itoh, Ryota Kakuma, Akio K. Inoue, Haruka Kusakabe, Kazuhiro Shimasaku, Rieko Momose, Yuma Sugahara, Satoshi Kikuta, Shun Saito, Nobunari Kashikawa, Haibin Zhang, Chien-Hsiu Lee

Abstract: We conduct intensity mapping to probe for extended diffuse Ly$α$ emission around Ly$α$ emitters (LAEs) at $z\sim2-7$, exploiting very deep ($\sim26$ mag at $5σ$) and large-area ($\sim4.5$ deg$^2$) Subaru/Hyper Suprime-Cam narrow-band (NB) images and large LAE catalogs consisting of a total of 1781 LAEs at $z=2.2$, $3.3$, $5.7$, and $6.6$ obtained by the HSC-SSP SILVERRUSH and CHORUS projects. We c… ▽ More We conduct intensity mapping to probe for extended diffuse Ly$α$ emission around Ly$α$ emitters (LAEs) at $z\sim2-7$, exploiting very deep ($\sim26$ mag at $5σ$) and large-area ($\sim4.5$ deg$^2$) Subaru/Hyper Suprime-Cam narrow-band (NB) images and large LAE catalogs consisting of a total of 1781 LAEs at $z=2.2$, $3.3$, $5.7$, and $6.6$ obtained by the HSC-SSP SILVERRUSH and CHORUS projects. We calculate the spatial correlations of these LAEs with $\sim1-2$ billion pixel flux values of the NB images, deriving the average Ly$α$ surface brightness (${\rm SB_{Lyα}}$) radial profiles around the LAEs. By carefully estimating systematics such as fluctuations of sky background and point spread functions, we detect diffuse Ly$α$ emission ($\sim10^{-20}-10^{-19}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$) at $100-1000$ comoving kpc around $z=3.3$ LAEs at the $4.1σ$ level and tentatively ($\sim2σ$) at the other redshifts, beyond the virial radius of a dark-matter halo with a mass of $10^{11}\ M_\odot$. While the observed ${\rm SB_{Lyα}}$ profiles have similar amplitudes at $z=2.2-6.6$ within the uncertainties, the intrinsic ${\rm SB_{Lyα}}$ profiles (corrected for the cosmological dimming effect) increase toward high redshifts. This trend may be explained by increasing hydrogen gas density due to the evolution of the cosmic volume. Comparisons with theoretical models suggest that extended Ly$α$ emission around a LAE is powered by resonantly scattered Ly$α$ photons in the CGM and IGM that originates from the inner part of the LAE, and/or neighboring galaxies around the LAE. △ Less

Submitted 20 August, 2021; originally announced August 2021.

GOLDRUSH. IV. Luminosity Functions and Clustering Revealed with ~4,000,000 Galaxies at z~2-7: Galaxy-AGN Transition, Star Formation Efficiency, and Implication for Evolution at z>10

Authors: Yuichi Harikane, Yoshiaki Ono, Masami Ouchi, Chengze Liu, Marcin Sawicki, Takatoshi Shibuya, Peter S. Behroozi, Wanqiu He, Kazuhiro Shimasaku, Stephane Arnouts, Jean Coupon, Seiji Fujimoto, Stephen Gwyn, Jiasheng Huang, Akio K. Inoue, Nobunari Kashikawa, Yutaka Komiyama, Yoshiki Matsuoka, Chris J. Willott

Abstract: We present new measurements of rest-UV luminosity functions and angular correlation functions from 4,100,221 galaxies at z~2-7 identified in the Subaru/Hyper Suprime-Cam survey and CFHT Large-Area U-band Survey. The obtained luminosity functions at z~4-7 cover a very wide UV luminosity range of ~0.002-2000L*uv combined with previous studies, revealing that the dropout luminosity function is a supe… ▽ More We present new measurements of rest-UV luminosity functions and angular correlation functions from 4,100,221 galaxies at z~2-7 identified in the Subaru/Hyper Suprime-Cam survey and CFHT Large-Area U-band Survey. The obtained luminosity functions at z~4-7 cover a very wide UV luminosity range of ~0.002-2000L*uv combined with previous studies, revealing that the dropout luminosity function is a superposition of the AGN luminosity function dominant at Muv<-24 mag and the galaxy luminosity function dominant at Muv>-22 mag, consistent with galaxy fractions based on 1037 spectroscopically-identified sources. Galaxy luminosity functions estimated from the spectroscopic galaxy fractions show the bright end excess beyond the Schechter function at >2sigma levels, which is possibly made by inefficient mass quenching, low dust obscuration, and/or hidden AGN activity. By analyzing the correlation functions at z~2-6 with halo occupation distribution models, we find a weak redshift evolution (within 0.3 dex) of the ratio of the star formation rate (SFR) to the dark matter accretion rate, SFR/(dMh/dt), indicating the almost constant star formation efficiency at z~2-6, as suggested by our earlier work at z~4-7. Meanwhile, the ratio gradually increases with decreasing redshift at z<5 within 0.3 dex, which quantitatively reproduces the redshift evolution of the cosmic SFR density, suggesting that the evolution is primarily driven by the increase of the halo number density due to the structure formation, and the decrease of the accretion rate due to the cosmic expansion. Extrapolating this calculation to higher redshifts assuming the constant efficiency suggests a rapid decrease of the SFR density at z>10 with $\propto10^{-0.5(1+z)}$, which will be directly tested with JWST. △ Less

Submitted 23 November, 2021; v1 submitted 2 August, 2021; originally announced August 2021.

Comments: 53 pages, 26 figures, 10 tables, accepted for publication in ApJS. Dropout galaxy catalogs are available on our project webpage (http://cos.icrr.u-tokyo.ac.jp/rush.html)

Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs). XIV. A Candidate Type-II Quasar at z=6.1292

Authors: Masafusa Onoue, Yoshiki Matsuoka, Nobunari Kashikawa, Michael A. Strauss, Kazushi Iwasawa, Takuma Izumi, Tohru Nagao, Naoko Asami, Seiji Fujimoto, Yuichi Harikane, Takuya Hashimoto, Masatoshi Imanishi, Chien-Hsiu Lee, Takatoshi, Shibuya, Yoshiki Toba

Abstract: We present deep Keck/MOSFIRE near-infrared spectroscopy of a strong Lyman alpha emitting source at z=6.1292, HSC J142331.71-001809.1, which was discovered by the SHELLQS program from imaging data of the Subaru Hyper Suprime-Cam (HSC) survey. This source is one of five objects that show unresolved (<230 km s-1) and prominent (>10^44 erg s-1) Lyman alpha emission lines at absolute 1450 angstrom cont… ▽ More We present deep Keck/MOSFIRE near-infrared spectroscopy of a strong Lyman alpha emitting source at z=6.1292, HSC J142331.71-001809.1, which was discovered by the SHELLQS program from imaging data of the Subaru Hyper Suprime-Cam (HSC) survey. This source is one of five objects that show unresolved (<230 km s-1) and prominent (>10^44 erg s-1) Lyman alpha emission lines at absolute 1450 angstrom continuum magnitudes of M1450~-22 mag. Its rest-frame Lyman alpha equivalent width (EW) is 370+/-30 angstrom. In the 2 hour Keck/MOSFIRE spectrum in Y band, the high-ionization CIV 1548,1550 doublet emission line was clearly detected with FWHM =120+/-20 km s-1 and a total rest-frame EW of 37-5+6 angstrom. We also report the detection of weak continuum emission, and the tentative detection of OIII] 1661,1666 in the 4 hour J band spectrum. Judging from the UV magnitude, line widths, luminosities, and EWs of Lyman alpha and CIV, we suggest that this source is a reionization-era analog of classical type-II AGNs, although there is a possibility that it represents a new population of AGN/galaxy composite objects in the early universe. We compare the properties of J1423-0018 to intermediate-redshift type-II AGNs and CIV emitters seen in z=6-7 galaxy samples. Further observations of other metal emission lines in the rest-frame UV or optical, and X-ray follow-up observations of the z=6-7 narrow-line quasars are needed for more robust diagnostics and to determine their nature. △ Less

Submitted 25 June, 2021; originally announced June 2021.

Comments: 15 pages, 8 figures, accepted for publication in the Astrophysical Journal

Optical Spectroscopy of Dual Quasar Candidates from the Subaru HSC-SSP program

Authors: Shenli Tang, John D. Silverman, Xuheng Ding, Junyao Li, Khee-Gan Lee, Michael A. Strauss, Andy Goulding, Malte Schramm, Lalitwadee Kawinwanichakij, J. Xavier Prochaska, Joseph F. Hennawi, Masatoshi Imanishi, Kazushi Iwasawa, Yoshiki Toba, Issha Kayo, Masamune Oguri, Yoshiki Matsuoka, Kohei Ichikawa, Tilman Hartwig, Nobunari Kashikawa, Toshihiro Kawaguchi, Kotaro Kohno, Yuichi Matsuda, Tohru Nagao, Yoshiaki Ono , et al. (8 additional authors not shown)

Abstract: We report on a spectroscopic program to search for dual quasars using Subaru Hyper Suprime-Cam (HSC) images of SDSS quasars which represent an important stage during galaxy mergers. Using Subaru/FOCAS and Gemini-N/GMOS, we identify three new physically associated quasar pairs having projected separations less than 20 kpc, out of 26 observed candidates. These include the discovery of the highest re… ▽ More We report on a spectroscopic program to search for dual quasars using Subaru Hyper Suprime-Cam (HSC) images of SDSS quasars which represent an important stage during galaxy mergers. Using Subaru/FOCAS and Gemini-N/GMOS, we identify three new physically associated quasar pairs having projected separations less than 20 kpc, out of 26 observed candidates. These include the discovery of the highest redshift ($z=3.1$) quasar pair with a separation $<$ 10 kpc. Based on the sample acquired to date, the success rate of identifying physically associated dual quasars is $19\%$ when excluding stars based on their HSC colors. Using the full sample of six spectroscopically confirmed dual quasars, we find that the black holes in these systems have black hole masses ($M_{BH} \sim 10^{8-9}M_{\odot}$) similar to single SDSS quasars as well as their bolometric luminosities and Eddington ratios. We measure the stellar mass of their host galaxies based on 2D image decomposition of the five-band ($grizy$) optical emission and assess the mass relation between supermassive black holes (SMBHs) and their hosts. Dual SMBHs appear to have elevated masses relative to their host galaxies. Thus mergers may not necessarily align such systems onto the local mass relation, as suggested by the Horizon-AGN simulation. This study suggests that dual luminous quasars are triggered prior to the final coalescence of the two SMBHs, resulting in early mass growth of the black holes relative to their host galaxies. △ Less

Submitted 21 May, 2021; originally announced May 2021.

Comments: 36 pages, 21 figures, 4 tables

Interrelation of the environment of Ly$α$ emitters and massive galaxies at $2<z<4.5$

Authors: Kei Ito, Nobunari Kashikawa, Masayuki Tanaka, Mariko Kubo, Yongming Liang, Jun Toshikawa, Hisakazu Uchiyama, Rikako Ishimoto, Takehiro Yoshioka, Yoshihiro Takeda

Abstract: We present a comparison of the spatial distributions of Ly$α$ emitters (LAEs) and massive star-forming and quiescent galaxies (SFGs and QGs) at $2<z<4.5$. We use the photometric redshift catalog to select SFGs and QGs and a LAE catalog from intermediate/narrow bands obtained from the Subaru Telescope and Isaac-Newton Telescope in Cosmic Evolution Survey (COSMOS). We derive the auto-/cross- correla… ▽ More We present a comparison of the spatial distributions of Ly$α$ emitters (LAEs) and massive star-forming and quiescent galaxies (SFGs and QGs) at $2<z<4.5$. We use the photometric redshift catalog to select SFGs and QGs and a LAE catalog from intermediate/narrow bands obtained from the Subaru Telescope and Isaac-Newton Telescope in Cosmic Evolution Survey (COSMOS). We derive the auto-/cross- correlation signals of SFGs, QGs, and LAEs, and the galaxy overdensity distributions at the position of them. Whereas the cross-correlation signals of SFGs and QGs are explained solely by their halo mass differences, those of SFGs and LAEs are significantly lower than those expected from their auto-correlation signals, suggesting that some additional physical processes are segregating these two populations. Such segregation of SFGs and LAEs becomes stronger for rest-frame ultraviolet faint LAEs ($M_{\rm UV}>-20$). From the overdensity distributions, LAEs are located in less dense regions than SFGs and QGs, whereas SFGs and QGs tend to be in the same overdensity distributions. The different spatial distributions of LAEs compared to those of massive galaxies may be attributed to assembly bias or large amounts of neutral hydrogen gas associated with massive halos. These results reinforce the importance of exploring multiple galaxy populations in quantifying the intrinsic galaxy environment of the high-$z$ universe. △ Less

Submitted 4 May, 2021; originally announced May 2021.

Comments: 21 pages, 14 figures, 3 tables, Accepted for publication in ApJ

Catch Me if You Can: Biased Distribution of Ly$α$-emitting Galaxies according to the Viewing Direction

Authors: Rieko Momose, Kazuhiro Shimasaku, Kentaro Nagamine, Ikkoh Shimizu, Nobunari Kashikawa, Makoto Ando, Haruka Kusakabe

Abstract: We report that Ly$α$-emitting galaxies (LAEs) may not faithfully trace the cosmic web of neutral hydrogen (HI), but their distribution is likely biased depending on the viewing direction. We calculate the cross-correlation (CCF) between galaxies and Ly$α$ forest transmission fluctuations on the near and far sides of the galaxies separately, for three galaxy samples at $z\sim2$: LAEs, [OIII] emitte… ▽ More We report that Ly$α$-emitting galaxies (LAEs) may not faithfully trace the cosmic web of neutral hydrogen (HI), but their distribution is likely biased depending on the viewing direction. We calculate the cross-correlation (CCF) between galaxies and Ly$α$ forest transmission fluctuations on the near and far sides of the galaxies separately, for three galaxy samples at $z\sim2$: LAEs, [OIII] emitters (O3Es), and continuum-selected galaxies. We find that only LAEs have anisotropic CCFs, with the near side one showing lower signals up to $r=3-4~h^{-1}$ comoving Mpc. This means that the average HI density on the near side of LAEs is lower than that on the far-side by a factor of $2.1$ under the Fluctuating Gunn-Peterson Approximation. Mock LAEs created by assigning Ly$α$ equivalent width ($EW_\text{Ly$α$}^\text{obs}$) values to O3Es with an empirical relation also show similar, anisotropic CCFs if we use only objects with higher $EW_\text{Ly$α$}^\text{obs}$ than a certain threshold. These results indicate that galaxies on the far side of a dense region are more difficult to be detected ("hidden") in Ly$α$ because Ly$α$ emission toward us is absorbed by dense neutral hydrogen. If the same region is viewed from a different direction, a different set of LAEs will be selected as if galaxies are playing hide-and-seek using HI gas. Care is needed when using LAEs to search for overdensities. △ Less

Submitted 21 April, 2021; originally announced April 2021.

Comments: 10 pages, 5 figures, Accepted for publication on ApJL

Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs). XIII. Large-scale Feedback and Star Formation in a Low-Luminosity Quasar at z = 7.07 on the Local Black Hole to Host Mass Relation

Authors: Takuma Izumi, Yoshiki Matsuoka, Seiji Fujimoto, Masafusa Onoue, Michael A. Strauss, Hideki Umehata, Masatoshi Imanishi, Kotaro Kohno, Toshihiro Kawaguchi, Taiki Kawamuro, Shunsuke Baba, Tohru Nagao, Yoshiki Toba, Kohei Inayoshi, John D. Silverman, Akio K. Inoue, Soh Ikarashi, Kazushi Iwasawa, Nobunari Kashikawa, Takuya Hashimoto, Kouichiro Nakanishi, Yoshihiro Ueda, Malte Schramm, Chien-Hsiu Lee, Hyewon Suh

Abstract: We present ALMA [CII] 158 $μ$m line and underlying far-infrared (FIR) continuum emission observations ($0''.70 \times 0''.56$ resolution) toward HSC J124353.93$+$010038.5 (J1243$+$0100) at $z = 7.07$, the only low-luminosity ($M_{\rm 1450} > -25$ mag) quasar currently known at $z > 7$. The FIR continuum is bright (1.52 mJy) and resolved with a total luminosity of… ▽ More We present ALMA [CII] 158 $μ$m line and underlying far-infrared (FIR) continuum emission observations ($0''.70 \times 0''.56$ resolution) toward HSC J124353.93$+$010038.5 (J1243$+$0100) at $z = 7.07$, the only low-luminosity ($M_{\rm 1450} > -25$ mag) quasar currently known at $z > 7$. The FIR continuum is bright (1.52 mJy) and resolved with a total luminosity of $L_{\rm FIR} = 3.5 \times 10^{12}~L_\odot$. The spatially extended component is responsible for $\sim 40\%$ of the emission. The area-integrated [CII] spectrum shows a broad wing (${\rm FWHM} = 997$ km s$^{-1}$, $L_{\rm [CII]} = 1.2 \times 10^9~L_\odot$) as well as a bright core (${\rm FWHM} = 235$ km s$^{-1}$, $L_{\rm [CII]} = 1.9 \times 10^9~L_\odot$). This wing is the first detection of a galactic-scale quasar-driven outflow (atomic outflow rate $> 447~M_\odot$ yr$^{-1}$) at $z > 7$. The estimated large mass loading factor of the total outflow (e.g., $\gtrsim 9$ relative to the [CII]-based SFR) suggests that this outflow will soon quench the star-formation of the host. The core gas dynamics are governed by rotation, with a rotation curve suggestive of a compact bulge ($\sim 3.3 \times 10^{10}~M_\odot$), although it is not yet spatially resolved. Finally, we found that J1243$+$0100 has a black hole mass-to-dynamical mass ratio (and -to-bulge mass ratio) of $\sim 0.4\%$ ($\sim 1\%$), consistent with the local value within uncertainties. Our results therefore suggest that the black hole-host co-evolution relation is already in place at $z \sim 7$ for this object. △ Less

Submitted 20 May, 2021; v1 submitted 12 April, 2021; originally announced April 2021.

Comments: Slightly revised version after the proof correction. 22 pages, 13 figures, 6 tables. Accepted for publication in ApJ

SILVERRUSH X: Machine Learning-Aided Selection of $9318$ LAEs at $z=2.2$, $3.3$, $4.9$, $5.7$, $6.6$, and $7.0$ from the HSC SSP and CHORUS Survey Data

Authors: Yoshiaki Ono, Ryohei Itoh, Takatoshi Shibuya, Masami Ouchi, Yuichi Harikane, Satoshi Yamanaka, Akio K. Inoue, Toshiyuki Amagasa, Daichi Miura, Maiki Okura, Kazuhiro Shimasaku, Ikuru Iwata, Yoshiaki Taniguchi, Seiji Fujimoto, Masanori Iye, Anton T. Jaelani, Nobunari Kashikawa, Shotaro Kikuchihara, Satoshi Kikuta, Masakazu A. R. Kobayashi, Haruka Kusakabe, Chien-Hsiu Lee, Yongming Liang, Yoshiki Matsuoka, Rieko Momose , et al. (3 additional authors not shown)

Abstract: We present a new catalog of $9318$ Ly$α$ emitter (LAE) candidates at $z = 2.2$, $3.3$, $4.9$, $5.7$, $6.6$, and $7.0$ that are photometrically selected by the SILVERRUSH program with a machine learning technique from large area (up to $25.0$ deg$^2$) imaging data with six narrowband filters taken by the Subaru Strategic Program with Hyper Suprime-Cam (HSC SSP) and a Subaru intensive program, Cosmi… ▽ More We present a new catalog of $9318$ Ly$α$ emitter (LAE) candidates at $z = 2.2$, $3.3$, $4.9$, $5.7$, $6.6$, and $7.0$ that are photometrically selected by the SILVERRUSH program with a machine learning technique from large area (up to $25.0$ deg$^2$) imaging data with six narrowband filters taken by the Subaru Strategic Program with Hyper Suprime-Cam (HSC SSP) and a Subaru intensive program, Cosmic HydrOgen Reionization Unveiled with Subaru (CHORUS). We construct a convolutional neural network that distinguishes between real LAEs and contaminants with a completeness of $94$% and a contamination rate of $1$%, enabling us to efficiently remove contaminants from the photometrically selected LAE candidates. We confirm that our LAE catalogs include $177$ LAEs that have been spectroscopically identified in our SILVERRUSH programs and previous studies, ensuring the validity of our machine learning selection. In addition, we find that the object-matching rates between our LAE catalogs and our previous results are $\simeq 80$-$100$% at bright NB magnitudes of $\lesssim 24$ mag. We also confirm that the surface number densities of our LAE candidates are consistent with previous results. Our LAE catalogs will be made public on our project webpage. △ Less

Submitted 5 April, 2021; originally announced April 2021.

Comments: 19 pages, 10 figures, accepted for publication in ApJ. Our LAE catalogs will become available at http://cos.icrr.u-tokyo.ac.jp/rush.html

"A more probable explanation" is still impossible to explain GN-z11-flash: in response to Steinhardt et al. (arXiv:2101.12738)

Authors: Linhua Jiang, Shu Wang, Bing Zhang, Nobunari Kashikawa, Luis C. Ho, Zheng Cai, Eiichi Egami, Gregory Walth, Yi-Si Yang, Bin-Bin Zhang, Hai-Bin Zhao

Abstract: In Jiang et al. (2020), we reported a possible bright flash (hereafter GN-z11-flash) from a galaxy GN-z11 at z ~ 11. Recently, Steinhardt et al. (2021; arXiv:2101.12738) found 27 images with transient signals in Keck MOSFIRE archival data and claimed that GN-z11-flash was more likely from a moving object in our Solar system. We show that the Steinhardt et al.'s definition of the chance probability… ▽ More In Jiang et al. (2020), we reported a possible bright flash (hereafter GN-z11-flash) from a galaxy GN-z11 at z ~ 11. Recently, Steinhardt et al. (2021; arXiv:2101.12738) found 27 images with transient signals in Keck MOSFIRE archival data and claimed that GN-z11-flash was more likely from a moving object in our Solar system. We show that the Steinhardt et al.'s definition of the chance probability and their methodology of finding GN-z11-flash-like transients are problematic in several aspects. In particular, none of their transients is analogous to GN-z11-flash, and none of them is positionally coincident with a known object in their imaging data. In Jiang et al., we performed a comprehensive analysis of the origin of GN-z11-flash and ruled out, to the best of our knowledge, the possibility of known man-made objects or moving objects in the Solar system, based on all available information and our current understanding of these objects. Steinhardt et al. did not use such information and did not analyse the GN-z11-flash event itself. The majority of their transients are apparently low-Earth orbit satellites or aircrafts. Therefore, their analysis can neither prove nor disprove our results. Finally, we present a method to estimate the chance probability of finding GN-z11-flash-like transients in archival data. Based on this method and the archival data used by Steinhardt et al., we obtain a loose upper limit of the probability that actually support the original results of Jiang et al. (2020). △ Less

Submitted 1 February, 2021; originally announced February 2021.

Comments: 5 pages, 1 figure, in response to Steinhardt et al. (arXiv:2101.12738)

Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) XII. Extended [C II] Structure (Merger or Outflow) in a z = 6.72 Red Quasar

Authors: Takuma Izumi, Masafusa Onoue, Yoshiki Matsuoka, Michael A. Strauss, Seiji Fujimoto, Hideki Umehata, Masatoshi Imanishi, Taiki Kawamuro, Tohru Nagao, Yoshiki Toba, Kotaro Kohno, Nobunari Kashikawa, Kohei Inayoshi, Toshihiro Kawaguchi, Kazushi Iwasawa, Akio K. Inoue, Tomotsugu Goto, Shunsuke Baba, Malte Schram, Hyewon Suh, Yuichi Harikane, Yoshihiro Ueda, John D. Silverman, Takuya Hashimoto, Yasuhiro Hashimoto , et al. (9 additional authors not shown)

Abstract: We present ALMA [C II] 158 $μ$m line and far-infrared (FIR) continuum emission observations toward HSC J120505.09$-$000027.9 (J1205$-$0000) at $z = 6.72$ with the beam size of $\sim 0''.8 \times 0''.5$ (or 4.1 kpc $\times$ 2.6 kpc), the most distant red quasar known to date. Red quasars are modestly reddened by dust, and are thought to be in rapid transition from an obscured starburst to an unobsc… ▽ More We present ALMA [C II] 158 $μ$m line and far-infrared (FIR) continuum emission observations toward HSC J120505.09$-$000027.9 (J1205$-$0000) at $z = 6.72$ with the beam size of $\sim 0''.8 \times 0''.5$ (or 4.1 kpc $\times$ 2.6 kpc), the most distant red quasar known to date. Red quasars are modestly reddened by dust, and are thought to be in rapid transition from an obscured starburst to an unobscured normal quasar, driven by powerful active galactic nucleus (AGN) feedback which blows out a cocoon of interstellar medium (ISM). The FIR continuum of J1205$-$0000 is bright, with an estimated luminosity of $L_{\rm FIR} \sim 3 \times 10^{12}~L_\odot$. The [C II] line emission is extended on scales of $r \sim 5$ kpc, greater than the FIR continuum. The line profiles at the extended regions are complex and broad (FWHM $\sim 630-780$ km s$^{-1}$). Although it is not practical to identify the nature of this extended structure, possible explanations include (i) companion/merging galaxies and (ii) massive AGN-driven outflows. For the case of (i), the companions are modestly star-forming ($\sim 10~M_\odot$ yr$^{-1}$), but are not detected by our Subaru optical observations ($y_{\rm AB,5σ} = 24.4$ mag). For the case of (ii), our lower-limit to the cold neutral outflow rate is $\sim 100~M_\odot$ yr$^{-1}$. The outflow kinetic energy and momentum are both much smaller than what predicted in energy-conserving wind models, suggesting that the AGN feedback in this quasar is not capable of completely suppressing its star formation. △ Less

Submitted 4 January, 2021; originally announced January 2021.

Comments: 12 pages, 4 figures, 2 tables, accepted for publication in ApJ

A possible bright ultraviolet flash from a galaxy at redshift z ~ 11

Authors: Linhua Jiang, Shu Wang, Bing Zhang, Nobunari Kashikawa, Luis C. Ho, Zheng Cai, Eiichi Egami, Gregory Walth, Yi-Si Yang, Bin-Bin Zhang, Hai-Bin Zhao

Abstract: In the optical sky, minutes-duration transients from cosmological distances are rare. Known objects that give rise to such transients include gamma-ray bursts (GRBs), the most luminous explosions in the universe that have been detected at redshift as high as z ~ 9.4. These high-redshift GRBs and their associated emission can be used to probe the star formation and reionization history in the era o… ▽ More In the optical sky, minutes-duration transients from cosmological distances are rare. Known objects that give rise to such transients include gamma-ray bursts (GRBs), the most luminous explosions in the universe that have been detected at redshift as high as z ~ 9.4. These high-redshift GRBs and their associated emission can be used to probe the star formation and reionization history in the era of cosmic dawn. Here we report a near-infrared transient with an observed duration shorter than 245 s coincident with the luminous star-forming galaxy GN-z11 at z ~ 11. The telluric absorption shown in the near-infrared spectrum indicates its origin from above the atmosphere. We can rule out the possibility of known man-made objects or moving objects in the Solar system based on the observational information and our current understanding of the properties of these objects. Since some long-duration GRBs are associated with a bright ultraviolet (UV) or optical flash, we investigate the possibility that the detected signal arose from a rest-frame UV flash associated with a long GRB from GN-z11. Despite the very low probability of being a GRB, we find that the spectrum, brightness, and duration of the transient are consistent with such an interpretation. Our result may suggest that long GRBs can be produced as early as 420 million years after the Big Bang. △ Less

Submitted 12 December, 2020; originally announced December 2020.

Comments: Published in Nature Astronomy on Dec 14, 2020; 19 pages; authors' version

Evidence for GN-z11 as a luminous galaxy at redshift 10.957

Authors: Linhua Jiang, Nobunari Kashikawa, Shu Wang, Gregory Walth, Luis C. Ho, Zheng Cai, Eiichi Egami, Xiaohui Fan, Kei Ito, Yongming Liang, Daniel Schaerer, Daniel P. Stark

Abstract: GN-z11 was photometrically selected as a luminous star-forming galaxy candidate at redshift z > 10 based on Hubble Space Telescope (HST) imaging data. Follow-up HST near-infrared grism observations detected a continuum break that was explained as the Ly-alpha break corresponding to z = 11.09 (+0.08-0.12). However, its accurate redshift remained unclear. Here we report a probable detection of three… ▽ More GN-z11 was photometrically selected as a luminous star-forming galaxy candidate at redshift z > 10 based on Hubble Space Telescope (HST) imaging data. Follow-up HST near-infrared grism observations detected a continuum break that was explained as the Ly-alpha break corresponding to z = 11.09 (+0.08-0.12). However, its accurate redshift remained unclear. Here we report a probable detection of three ultraviolet (UV) emission lines from GN-z11, which can be interpreted as the [C III] 1907, C III] 1909 doublet and O III] 1666 at z = 10.957+/-0.001 (when the Universe was only ~420 Myr old, or ~3% of its current age). This is consistent with the redshift of the previous grism observations, supporting GN-z11 as the most distant galaxy known to date. Its UV lines likely originate from dense ionized gas that is rarely seen at low redshifts, and its strong [C III] and C III] emission is partly due to an active galactic nucleus (AGN) or enhanced carbon abundance. GN-z11 is luminous and young, yet moderately massive, implying a rapid build-up of stellar mass in the past. Future facilities will be able to find the progenitors of such galaxies at higher redshift and probe the cosmic epoch in the beginning of re-ionization. △ Less

Submitted 12 December, 2020; originally announced December 2020.

Comments: Published in Nature Astronomy on Dec 14, 2020; 21 pages; authors' version

FOREVER22: galaxy formation in protocluster regions

Authors: Hidenobu Yajima, Makito Abe, Sadegh Khochfar, Kentaro Nagamine, Akio K. Inoue, Tadayuki Kodama, Shohei Arata, Claudio Dalla-Vecchia, Hajime Fukushima, Takuya Hashimoto, Nobunari Kashikawa, Mariko Kubo, Yuexing Li, Yuichi Matsuda, Ken Mawatari, Masami Ouchi, Hideki Umehata

Abstract: We present results from a new cosmological hydrodynamics simulation campaign of protocluster (PC) regions, FOREVER22: FORmation and EVolution of galaxies in Extremely-overdense Regions motivated by SSA22. The simulations cover a wide range of cosmological scales using three different zoom set-ups in a parent volume of $(714.2~\rm cMpc)^{3}$: PCR (Proto-Cluster Region; $V= (28.6~{\rm cMpc})^{3} $,… ▽ More We present results from a new cosmological hydrodynamics simulation campaign of protocluster (PC) regions, FOREVER22: FORmation and EVolution of galaxies in Extremely-overdense Regions motivated by SSA22. The simulations cover a wide range of cosmological scales using three different zoom set-ups in a parent volume of $(714.2~\rm cMpc)^{3}$: PCR (Proto-Cluster Region; $V= (28.6~{\rm cMpc})^{3} $, SPH particle mass, $m_{\rm{SPH}} = 4.1 \times 10^{6}~\rm M_{\odot}$ and final redshift, $z_{\rm end}=2.0$), BCG (Brightest proto-Cluster Galaxy; $V \sim (10~{\rm cMpc})^{3} $, $m_{\rm SPH} = 5.0\times10^{5}~\rm M_{\odot}$ and $z_{\rm end}=4.0$ ), and First ( $V \sim (3~{\rm cMpc})^{3} $, $m_{\rm SPH} = 7.9 \times 10^{3}~\rm M_{\odot}$ and $z_{\rm end}=9.5$) runs, that allow to focus on different aspects of galaxy formation. In the PCR runs, we follow 10 PCs, each harbouring 1 - 4 SMBHs with $M_{\rm BH} \ge 10^{9}~\rm M_{\odot}$. One of the PC cores shows a spatially close arrangement of seven starburst galaxies with ${\rm SFR} \gtrsim 100~\rm M_{\odot}~yr^{-1}$ each, that are dust-obscured and would appear as submillimeter galaxies with flux $\gtrsim 1~$ mJy at $1.1~ \rm mm$ in observations. The BCG runs show that the total SFRs of haloes hosting BCGs are affected by AGN feedback, but exceed $1000~\rm M_{\odot}~yr^{-1}$ at $z \lesssim 6$. The First runs resolve mini-haloes hosting population (Pop) III stars and we show that, in PC regions, the dominant stellar population changes from Pop III to Pop II at $z \gtrsim 20$, and the first galaxies with ${\rm SFR} \gtrsim 18~\rm M_{\odot}~yr^{-1}$ form at $z \sim 10$. These can be prime targets for future observations with the James Webb Space Telescope. Our simulations successfully reproduce the global star formation activities in observed PCs and suggest that PCs can kickstart cosmic reionization. △ Less

Submitted 7 January, 2022; v1 submitted 23 November, 2020; originally announced November 2020.

Comments: 24 pages, 19 figures, accepted for publication in MNRAS

Faint Quasars Live in the Same Number Density Environments as Lyman Break Galaxies at z ~ 4

Authors: Hisakazu Uchiyama, Masayuki Akiyama, Jun Toshikawa, Nobunari Kashikawa, Roderik Overzier, Tohru Nagao, Kohei Ichikawa, Murilo Marinello, Masatoshi Imanishi, Masayuki Tanaka, Yoshiki Matsuoka, Yutaka Komiyama, Shogo Ishikawa, Masafusa Onoue, Mariko Kubo, Yuichi Harikane, Kei Ito, Shigeru Namiki, Yongming Liang

Abstract: Characterizing high-z quasar environments is key to understanding the co-evolution of quasars and the surrounding galaxies. To restrict their global picture, we statistically examine the g-dropout galaxy overdensity distribution around 570 faint quasar candidates at z ~ 4, based on the Hyper Suprime-Cam Subaru Strategic Program survey. We compare the overdensity significances of g-dropout galaxies… ▽ More Characterizing high-z quasar environments is key to understanding the co-evolution of quasars and the surrounding galaxies. To restrict their global picture, we statistically examine the g-dropout galaxy overdensity distribution around 570 faint quasar candidates at z ~ 4, based on the Hyper Suprime-Cam Subaru Strategic Program survey. We compare the overdensity significances of g-dropout galaxies around the quasars with those around g-dropout galaxies, and find no significant difference between their distributions. A total of 4 (22) out of the 570 faint quasars, 0.7_{-0.4}^{+0.4} (3.9_{-0.8}^{+0.8}) %, are found to be associated with the > 4 sigma overdense regions within an angular separation of 1.8 (3.0) arcmin, which is the typical size of protoclusters at this epoch. This is similar to the fraction of g-dropout galaxies associated with the > 4 sigma overdense regions. This result is consistent with our previous work that 1.3_{-0.9}^{+0.9} % and 2.0_{-1.1}^{+1.1} % of luminous quasars detected in the Sloan Digital Sky Survey exist in the > 4 sigma overdense regions within 1.8 and 3.0 arcmin separations, respectively. Therefore, we suggest that the galaxy number densities around quasars are independent of their luminosity, and most quasars do not preferentially appear in the richest protocluster regions at z ~ 4. The lack of an apparent positive correlation between the quasars and the protoclusters implies that: i) the gas-rich major merger rate is relatively low in the protocluster regions, ii) most high-z quasars may appear through secular processes, or iii) some dust-obscured quasars exist in the protocluster regions. △ Less

Submitted 19 November, 2020; originally announced November 2020.

Comments: 12 pages, 4 figures

CHORUS. I. Cosmic HydrOgen Reionization Unveiled with Subaru: Overview

Authors: Akio K. Inoue, Satoshi Yamanaka, Masami Ouchi, Ikuru Iwata, Kazuhiro Shimasaku, Yoshiaki Taniguchi, Tohru Nagao, Nobunari Kashikawa, Yoshiaki Ono, Ken Mawatari, Takatoshi Shibuya, Masao Hayashi, Hiroyuki Ikeda, Haibin Zhang, Yongming Liang, C. -H. Lee, Miftahul Hilmi, Satoshi Kikuta, Haruka Kusakabe, Hisanori Furusawa, Tomoki Hayashino, Masaru Kajisawa, Yuichi Matsuda, Kimihiko Nakajima, Rieko Momose , et al. (6 additional authors not shown)

Abstract: To determine the dominant sources for cosmic reionization, the evolution history of the global ionizing fraction, and the topology of the ionized regions, we have conducted a deep imaging survey using four narrow-band (NB) and one intermediate-band (IB) filters on the Subaru/Hyper Suprime-Cam (HSC), called Cosmic HydrOgen Reionization Unveiled with Subaru (CHORUS). The central wavelengths and full… ▽ More To determine the dominant sources for cosmic reionization, the evolution history of the global ionizing fraction, and the topology of the ionized regions, we have conducted a deep imaging survey using four narrow-band (NB) and one intermediate-band (IB) filters on the Subaru/Hyper Suprime-Cam (HSC), called Cosmic HydrOgen Reionization Unveiled with Subaru (CHORUS). The central wavelengths and full-widths-at-half-maximum of the CHORUS filters are, respectively, 386.2 nm and 5.5 nm for NB387, 526.0 nm and 7.9 nm for NB527, 717.1 nm and 11.1 nm for NB718, 946.2 nm and 33.0 nm for IB945, and 971.2 nm and 11.2 nm for NB973. This combination, including NB921 (921.5 nm and 13.5 nm) from the Subaru Strategic Program with HSC (HSC SSP), are carefully designed, as if they were playing a chorus, to observe multiple spectral features simultaneously, such as Lyman continuum, Ly$α$, C~{\sc iv}, and He~{\sc ii} for $z=2$--$7$. The observing field is the same as that of the deepest footprint of the HSC SSP in the COSMOS field and its effective area is about 1.6 deg$^2$. Here, we present an overview of the CHORUS project, which includes descriptions of the filter design philosophy, observations and data reduction, multiband photometric catalogs, assessments of the imaging quality, measurements of the number counts, and example use cases of the data. All the imaging data, photometric catalogs, masked pixel images, data of limiting magnitudes and point spread functions, results of completeness simulations, and source number counts are publicly available through the HSC SSP database. △ Less

Submitted 13 November, 2020; originally announced November 2020.

Comments: 17 pages, 6 figures, PASJ open access (https://academic.oup.com/pasj/advance-article-abstract/doi/10.1093/pasj/psaa100/5959796), all the data products are available at https://hsc-release.mtk.nao.ac.jp/doc/index.php/chorus

Journal ref: PASJ, psaa100, 2020

The Faint End of the Quasar Luminosity Function at $z \sim 5$ from the Subaru Hyper Suprime-Cam Survey

Authors: Mana Niida, Tohru Nagao, Hiroyuki Ikeda, Masayuki Akiyama, Yoshiki Matsuoka, Wanqiu He, Kenta Matsuoka, Yoshiki Toba, Masafusa Onoue, Masakazu A. R. Kobayashi, Yoshiaki Taniguchi, Hisanori Furusawa, Yuichi Harikane, Masatoshi Imanishi, Nobunari Kashikawa, Toshihiro Kawaguchi, Yutaka Komiyama, Hikari Shirakata, Yuichi Terashima, Yoshihiro Ueda

Abstract: We present the quasar luminosity function at $z \sim 5$ derived from the optical wide-field survey data obtained as a part of the Subaru strategic program (SSP) with Hyper Suprime-Cam (HSC). From $\sim$81.8 deg$^2$ area in the Wide layer of the HSC-SSP survey, we selected 224 candidates of low-luminosity quasars at $z \sim 5$ by adopting the Lyman-break method down to $i = 24.1$ mag. Based on our… ▽ More We present the quasar luminosity function at $z \sim 5$ derived from the optical wide-field survey data obtained as a part of the Subaru strategic program (SSP) with Hyper Suprime-Cam (HSC). From $\sim$81.8 deg$^2$ area in the Wide layer of the HSC-SSP survey, we selected 224 candidates of low-luminosity quasars at $z \sim 5$ by adopting the Lyman-break method down to $i = 24.1$ mag. Based on our candidates and spectroscopically-confirmed quasars from the Sloan Digital Sky Survey (SDSS), we derived the quasar luminosity function at $z \sim 5$ covering a wide luminosity range of $-28.76 < M_{\rm 1450} < -22.32$ mag. We found that the quasar luminosity function is fitted by a double power-law model with a break magnitude of $M^{*}_{1450} = -25.05^{+0.10}_{-0.24}$ mag. The inferred number density of low-luminosity quasars is lower, and the derived faint-end slope, $-1.22^{+0.03}_{-0.10}$, is flatter than those of previous studies at $z \sim 5$. A compilation of the quasar luminosity function at $4 \leq z \leq 6$ from the HSC-SSP suggests that there is little redshift evolution in the break magnitude and in the faint-end slope within this redshift range, although previous studies suggest that the faint-end slope becomes steeper at higher redshifts. The number density of low-luminosity quasars decreases more rapidly from $z \sim 5$ to $z \sim 6$ than from $z \sim 4$ to $z \sim 5$. △ Less

Submitted 1 October, 2020; originally announced October 2020.

Comments: 18 pages, 11 figures, accepted for publication in ApJ

Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs). XI. Proximity Zone Analysis for Faint Quasar Spectra at z~6

Authors: Rikako Ishimoto, Nobunari Kashikawa, Masafusa Onoue, Yoshiki Matsuoka, Takuma Izumi, Michael A. Strauss, Seiji Fujimoto, Masatoshi Imanishi, Kei Ito, Kazushi Iwasawa, Toshihiro Kawaguchi, Chien-Hsiu Lee, Yongming Liang, Ting-Yi Lu, Rieko Momose, Yoshiki Toba, Hisakazu Uchiyama

Abstract: We present measurements of the size of the quasar proximity zone ($R_p$) for eleven low-luminosity ($-26.16\leq M_{1450}\leq-22.83$) quasars at $z\sim6$, discovered by the Subaru High-$z$ Exploration of Low-Luminosity Quasars project (SHELLQs). Our faint quasar sample expands the $R_p$ measurement down to $M_{1450}=-22.83$ mag, where more common quasar populations dominate at the epoch. We restric… ▽ More We present measurements of the size of the quasar proximity zone ($R_p$) for eleven low-luminosity ($-26.16\leq M_{1450}\leq-22.83$) quasars at $z\sim6$, discovered by the Subaru High-$z$ Exploration of Low-Luminosity Quasars project (SHELLQs). Our faint quasar sample expands the $R_p$ measurement down to $M_{1450}=-22.83$ mag, where more common quasar populations dominate at the epoch. We restrict the sample to quasars whose systemic redshifts have been precisely measured by [CII] 158 $μ$m or MgII $λ$2798 emission lines. We also update the $R_p$ measurements for 26 luminous quasars presented in Eilers et al. (2017)(arXiv:1703.02539) by using the latest systemic redshift results. The luminosity dependence on $R_p$ is found to be consistent with the theoretical prediction assuming highly ionized intergalactic medium. We find a shallow redshift evolution of the luminosity-corrected $R_p$, $R_{p,{\rm corr}}^{-25}$ ($R_{p, \rm corr}^{-25}\propto(1+z)^{-3.79\pm1.72}$) over $5.8\lesssim z \lesssim6.6$. This trend is steeper than that of Eilers et al. (2017) but significantly shallower than those of the earlier studies. Our results suggest that $R_{p,\rm corr}$ is insensitive to the neutral fraction of the universe at $z\sim6$. Four quasars show exceptionally small $R_{p,\rm corr}^{-25}$ ($ \lesssim0.90$ proper Mpc), which could be the result of their young age ($<10^4$ yr) in the reionization epoch, though statistics is still small. △ Less

Submitted 14 September, 2020; originally announced September 2020.

Comments: 14 pages, 8 figures, Accepted for publication in ApJ

Statistical correlation between the distribution of Ly$α$ emitters and IGM HI at $z\sim2.2$ mapped by Subaru/Hyper Suprime-Cam

Authors: Yongming Liang, Nobunari Kashikawa, Zheng Cai, Xiaohui Fan, J. Xavier Prochaska, Kazuhiro Shimasaku, Masayuki Tanaka, Hisakazu Uchiyama, Kei Ito, Rhythm Shimakawa, Kentaro Nagamine, Ikkoh Shimizu, Masafusa Onoue, Jun Toshikawa

Abstract: The correlation between neutral Hydrogen (HI) in the intergalactic medium (IGM) and galaxies now attracts great interests. We select four fields which include several coherently strong Ly$α$ absorption systems at $z\sim2.2$ detected by using background quasars from the whole SDSS/(e)BOSS database. Deep narrow-band and $g$-band imaging are performed using the Hyper Suprime-Cam on the Subaru Telesco… ▽ More The correlation between neutral Hydrogen (HI) in the intergalactic medium (IGM) and galaxies now attracts great interests. We select four fields which include several coherently strong Ly$α$ absorption systems at $z\sim2.2$ detected by using background quasars from the whole SDSS/(e)BOSS database. Deep narrow-band and $g$-band imaging are performed using the Hyper Suprime-Cam on the Subaru Telescope. We select out 2,642 Ly$α$ emitter (LAE) candidates at $z=2.177\pm0.023$ down to the Ly$α$ luminosity of $L_{\text{Ly}α}\approx 2 \times 10^{42} {\rm erg~s}^{-1}$ to construct the galaxy overdensity maps, covering an effective area of 5.39 deg$^2$. Combining the sample with the Ly$α$ absorption estimated from 64 (e)BOSS quasar spectra, we find a moderate to strong correlation between the LAE overdensity $δ_{\rm LAE}$ and the effective optical depth $τ_{\rm LoS}$ in line-of-sights, with $P$-value$=0.09\%$ ($<0.01\%$) when the field that contains a significant quasar overdensity is in(ex)cluded. The cross-correlation analysis also clearly suggests that up to $4\pm1$ pMpc, LAEs tend to cluster in the regions rich in HI gas, indicated by the high $τ_{\rm LoS}$, and avoid the low $τ_{\rm LoS}$ region where the HI gas is deficient. By averaging the $τ_{\rm LoS}$ as a function of the projected distance ($d$) to LAEs, we find a $30\%$ excess signal at $2σ$ level at $d<200$ pkpc, indicating the dense HI in circumgalactic medium, and a tentative excess at $400<d<600$ pkpc in IGM regime, corroborating the cross-correlation signal detected at about $0.5$ pMpc. These statistical analyses indicate that galaxy$-$IGM HI correlations exist on scales ranging from several hundred pkpc to several pMpc at $z\sim2.2$. △ Less

Submitted 4 August, 2020; originally announced August 2020.

Comments: 27 pages, 15 figures; Resubmitted to ApJ after the first referee's report. Comments are welcome

Probing Feedback via IGM tomography and Ly$α$ forest with Subaru PFS, TMT/ELT, and JWST

Authors: Kentaro Nagamine, Ikkoh Shimizu, Katsumi Fujita, Nao Suzuki, Khee-Gan Lee, Rieko Momose, Shiro Mukae, Yongming Liang, Nobunari Kashikawa, Masami Ouchi, John Silverman

Abstract: In preparation for the IGM tomography study by Subaru Prime Focus Spectrograph (PFS) survey and other large future telescopes such as TMT/ELT/GMT, we present the results of our pilot study on Ly$α$ forest and IGM tomography statistics using the GADGET3-OSAKA cosmological smoothed particle hydrodynamical simulation. Our simulation includes models for star formation and supernova feedback, which ena… ▽ More In preparation for the IGM tomography study by Subaru Prime Focus Spectrograph (PFS) survey and other large future telescopes such as TMT/ELT/GMT, we present the results of our pilot study on Ly$α$ forest and IGM tomography statistics using the GADGET3-OSAKA cosmological smoothed particle hydrodynamical simulation. Our simulation includes models for star formation and supernova feedback, which enables more realistic cross-correlation studies between galaxies, neutral hydrogen (HI) and metals in circumgalactic and intergalactic medium. We create a light-cone data set at $z=2-3$ from our simulations and generate mock Ly$α$ forest data. As a first step, in this paper, we focus on the distribution of HI and galaxies, and present statistical results on 1D flux PDF, 1D power spectrum, flux contrast vs. impact parameter, HI$-$galaxy cross-correlations. Our results show overall agreement with current observational data, with some interesting discrepancies on small scales that are due to either feedback effects or varying observational conditions. Our simulation shows stronger HI absorption with decreasing transverse distance from galaxies. We find that the massive galaxies with $M_\star \ge 10^{10}\,M_\odot$ contribute strongly to the flux contrast signal, and that the lower-mass galaxies with $M_\star \sim 10^8-10^{10}\,M_\odot$ tend to dilute the flux contrast signal from massive galaxies. On large scales, the average flux contrast smoothly connects to the IGM level, supporting the concordance $Λ$ cold dark matter model. We also find an increase in the HI absorption toward the center of a protocluster. △ Less

Submitted 7 July, 2021; v1 submitted 28 July, 2020; originally announced July 2020.

Comments: 26 pages, 16 figures. Published in ApJ (matched to the published version)

Journal ref: ApJ, 914, 66 (2021)

Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs). IX. Identification of Two Red Quasars at z > 5.6

Authors: Nanako Kato, Yoshiki Matsuoka, Masafusa Onoue, Shuhei Koyama, Yoshiki Toba, Masayuki Akiyama, Seiji Fujimoto, Masatoshi Imanishi, Kazushi Iwasawa, Takuma Izumi, Nobunari Kashikawa, Toshihiro Kawaguchi, Chien-Hsiu Lee, Takeo Minezaki, Tohru Nagao, Akatoki Noboriguchi, Michael A. Strauss

Abstract: We present the first discovery of dust-reddened quasars (red quasars) in the high-z universe (z >5.6). This is a result from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which is based on the sensitive multi-band optical imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. We identified four red quasar candidates from the spectroscopic… ▽ More We present the first discovery of dust-reddened quasars (red quasars) in the high-z universe (z >5.6). This is a result from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which is based on the sensitive multi-band optical imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. We identified four red quasar candidates from the spectroscopically confirmed 93 high-z quasars in the SHELLQs sample, based on detections in the Wide-field Infrared Survey Explorer (WISE) data at 3.4 and 4.6 um (rest-frame ~5000-6500 A). The amount of dust reddening was estimated with spectral energy distribution (SED) fits over optical and mid-infrared wavelengths. Two of the four candidates were found to be red quasars with dust reddening of E(B-V) > 0.1. The remaining SHELLQs quasars without individual WISE detections are significantly fainter in the WISE bands and bluer than the red quasars, although we did detect them in the W1 band in a stacked image. We also conducted the same SED fits for high-z optically-luminous quasars, but no red quasar was found. This demonstrates the power of Subaru HSC to discover high-z red quasars, which are fainter than the limiting magnitudes of past surveys in the rest-frame ultraviolet, due to dust extinction. △ Less

Submitted 28 July, 2020; v1 submitted 16 July, 2020; originally announced July 2020.

Comments: Accepted for publication in PASJ

Dual supermassive black holes at close separation revealed by the Hyper Suprime-Cam Subaru Strategic Program

Authors: John D. Silverman, Shenli Tang, Khee-Gan Lee, Tilman Hartwig, Andy Goulding, Michael A. Strauss, Malte Schramm, Xuheng Ding, Rogemar Riffel, Seiji Fujimoto, Chiaki Hikage, Masatoshi Imanishi, Kazushi Iwasawa, Knud Jahnke, Issha Kayo, Nobunari Kashikawa, Toshihiro Kawaguchi, Kotaro Kohno, Wentao Luo, Yoshiki Matsuoka, Yuichi Matsuda, Tohru Nagao, Masamune Oguri, Yoshiaki Ono, Masafusa Onoue , et al. (8 additional authors not shown)

Abstract: The unique combination of superb spatial resolution, wide-area coverage, and deep depth of the optical imaging from the Hyper Suprime-Cam (HSC) Subaru Strategic Program is utilized to search for dual quasar candidates. Using an automated image analysis routine on 34,476 known SDSS quasars, we identify those with two (or more) distinct optical point sources in HSC images covering 796 deg^2. We find… ▽ More The unique combination of superb spatial resolution, wide-area coverage, and deep depth of the optical imaging from the Hyper Suprime-Cam (HSC) Subaru Strategic Program is utilized to search for dual quasar candidates. Using an automated image analysis routine on 34,476 known SDSS quasars, we identify those with two (or more) distinct optical point sources in HSC images covering 796 deg^2. We find 421 candidates out to a redshift of 4.5 of which one hundred or so are more likely after filtering out contaminating stars. Angular separations of 0.6 - 4.0" correspond to projected separations of 3 - 30 kpc, a range relatively unexplored for population studies of luminous dual quasars. Using Keck-I/LRIS and Gemini-N/NIFS, we spectroscopically confirm three dual quasar systems at z < 1, two of which are previously unknown out of eight observed, based on the presence of characteristic broad emission lines in each component, while highlighting that the continuum of one object in one of the pairs is reddened. In all cases, the [OIII]5007 emission lines have mild velocity offsets, thus the joint [OIII] line profile is not double-peaked. We find a dual quasar fraction of 0.26+/-0.18% and no evidence for evolution. A comparison with the Horizon-AGN simulation seems to support the case of no evolution in the dual quasar fraction when broadly matching the quasar selection. These results may indicate a scenario in which the frequency of the simultaneous triggering of luminous quasars is not as sensitive as expected to the cosmic evolution of the merger rate or gas content of galaxies. △ Less

Submitted 10 July, 2020; originally announced July 2020.

Comments: 11 pages; 12 figures; Accepted for publication in The Astrophysical Journal

The UV Luminosity Function of Protocluster Galaxies at $z\sim4$: the Bright-end Excess and the Enhanced Star Formation Rate Density

Authors: Kei Ito, Nobunari Kashikawa, Jun Toshikawa, Roderik Overzier, Mariko Kubo, Hisakazu Uchiyama, Yongming Liang, Masafusa Onoue, Masayuki Tanaka, Yutaka Komiyama, Chien-Hsiu Lee, Yen-Ting Lin, Murilo Marinello, Crystal L. Martin, Takatoshi Shibuya

Abstract: We report the rest-frame ultraviolet luminosity function of $g$-dropout galaxies in 177 protocluster candidates (PC UVLF) at $z\sim4$ selected in the Hyper Suprime-Cam Subaru Strategic Program. Comparing with the UVLF of field galaxies at the same redshift, we find that the PC UVLF shows a significant excess towards the bright-end. This excess can not be explained by the contribution of only activ… ▽ More We report the rest-frame ultraviolet luminosity function of $g$-dropout galaxies in 177 protocluster candidates (PC UVLF) at $z\sim4$ selected in the Hyper Suprime-Cam Subaru Strategic Program. Comparing with the UVLF of field galaxies at the same redshift, we find that the PC UVLF shows a significant excess towards the bright-end. This excess can not be explained by the contribution of only active galactic nuclei, and we also find that this is more significant in higher dense regions. Assuming that all protocluster members are located on the star formation main sequence, the PC UVLF can be converted into a stellar mass function. Consequently, our protocluster members are inferred to have a 2.8 times more massive characteristic stellar mass than that of the field Lyman break galaxies at the same redshift. This study, for the first time, clearly shows that the enhancement in star formation or stellar mass in overdense regions can generally be seen as early as at $z\sim4$. We also estimate the star formation rate density (SFRD) in protocluster regions as $\simeq 6-20\%$ of the cosmic SFRD, based on the measured PC UVLF after correcting for the selection incompleteness in our protocluster sample. This high value suggests that protoclusters make a non-negligible contribution to the cosmic SFRD at $z\sim4$, as previously suggested by simulations. Our results suggest that protoclusters are essential components for the galaxy evolution at $z\sim4$. △ Less

Submitted 6 July, 2020; originally announced July 2020.

Comments: 22 pages, 16 figures, Accepted for publication in ApJ