Showing 1–50 of 121 results for author: Matsuoka, Y

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

Euclid. I. Overview of the Euclid mission

Authors: Euclid Collaboration, Y. Mellier, Abdurro'uf, J. A. Acevedo Barroso, A. Achúcarro, J. Adamek, R. Adam, G. E. Addison, N. Aghanim, M. Aguena, V. Ajani, Y. Akrami, A. Al-Bahlawan, A. Alavi, I. S. Albuquerque, G. Alestas, G. Alguero, A. Allaoui, S. W. Allen, V. Allevato, A. V. Alonso-Tetilla, B. Altieri, A. Alvarez-Candal, A. Amara, L. Amendola , et al. (1086 additional authors not shown)

Abstract: The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14… ▽ More The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue`Euclid on Sky'

Merging gas-rich galaxies that harbor low-luminosity twin quasars at z = 6.05: a promising progenitor of the most luminous quasars

Authors: Takuma Izumi, Yoshiki Matsuoka, Masafusa Onoue, Michael A. Strauss, Hideki Umehata, John D. Silverman, Tohru Nagao, Masatoshi Imanishi, Kotaro Kohno, Yoshiki Toba, Kazushi Iwasawa, Kouichiro Nakanishi, Mahoshi Sawamura, Seiji Fujimoto, Satoshi Kikuta, Toshihiro Kawaguchi, Kentaro Aoki, Tomotsugu Goto

Abstract: We present ALMA [CII] 158 $μ$m line and underlying far-infrared continuum emission observations ($0''.57 \times 0''.46$ resolution) toward a quasar-quasar pair system recently discovered at $z = 6.05$ (Matsuoka et al. 2024). The quasar nuclei (C1 and C2) are very faint ($M_{\rm 1450} \gtrsim -23$ mag), but we detect very bright [CII] emission bridging the 12 kpc between the two objects and extendi… ▽ More We present ALMA [CII] 158 $μ$m line and underlying far-infrared continuum emission observations ($0''.57 \times 0''.46$ resolution) toward a quasar-quasar pair system recently discovered at $z = 6.05$ (Matsuoka et al. 2024). The quasar nuclei (C1 and C2) are very faint ($M_{\rm 1450} \gtrsim -23$ mag), but we detect very bright [CII] emission bridging the 12 kpc between the two objects and extending beyond them (total luminosity $L_{\rm [CII]} \simeq 6 \times 10^9~L_\odot$). The total star formation rate of the system is $\sim 100$ to 550 $M_\odot$ yr$^{-1}$, with a [CII]-based total gas mass of $\sim 10^{11}~M_\odot$. The dynamical masses of the two galaxies are large ($\sim 9 \times 10^{10}~M_\odot$ for C1 and $\sim 5 \times 10^{10}~M_\odot$ for C2). There is a smooth velocity gradient in [CII], indicating that these quasars are a tidally interacting system. We identified a dynamically distinct, fast [CII] component around C1. Detailed inspection of the line spectrum reveals the presence of a broad wing component, which we interpret as the indication of fast outflows with a velocity of $\sim 600$ km s$^{-1}$. The expected mass loading factor of the outflows, after accounting for multiphase gas, would be $\gtrsim 3-8$, as expected for an AGN-driven outflow. Hydrodynamic simulations in the literature predicted that this pair will evolve to a luminous ($M_{\rm 1450} \lesssim -26$ mag), starbursting ($\gtrsim 1000~M_\odot$ yr$^{-1}$) quasar after coalescence, one of the most extreme populations in the early universe. △ Less

Submitted 3 May, 2024; originally announced May 2024.

Comments: 14 pages incl. Appendix, 4 main figures, 1 table, 1 appendix figure. Submitted to AAS journal

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)

The relationship of SMBHs and host galaxies at $z<4$ in the deep optical variability-selected AGN sample in the COSMOS field

Authors: Atsushi Hoshi, Toru Yamada, Mitsuru Kokubo, Yoshiki Matsuoka, Tohru Nagao

Abstract: We present the study on the relationship between SMBHs and their host galaxies using our variability-selected AGN sample ($i_\mathrm{AB} \leq 25.9,\ z \leq 4.5$) constructed from the HSC-SSP Ultra-deep survey in the COSMOS field. We estimated the BH mass ($M_\mathrm{BH}=10^{5.5-10}\ M_{\odot}$) based on the single-epoch virial method and the total stellar mass (… ▽ More We present the study on the relationship between SMBHs and their host galaxies using our variability-selected AGN sample ($i_\mathrm{AB} \leq 25.9,\ z \leq 4.5$) constructed from the HSC-SSP Ultra-deep survey in the COSMOS field. We estimated the BH mass ($M_\mathrm{BH}=10^{5.5-10}\ M_{\odot}$) based on the single-epoch virial method and the total stellar mass ($M_\mathrm{star}=10^{10-12}\ M_{\odot}$) by separating the AGN component with SED fitting. We found that the redshift evolution of the BH-stellar mass ratio ($M_\mathrm{BH}/M_\mathrm{star}$) depends on the $M_\mathrm{BH}$ which is caused by the no significant correlation between $M_\mathrm{BH}$ and $M_\mathrm{star}$. Variable AGNs with massive SMBHs ($M_\mathrm{BH}>10^{9}\ M_{\odot}$) at $1.5<z<3$ show considerably higher BH-stellar mass ratios ($>\sim1\%$) than the BH-bulge ratios ($M_\mathrm{BH}/M_\mathrm{bulge}$) observed in the local universe for the same BH range. This implies that there is a typical growth path of massive SMBHs which is faster than the formation of the bulge component as final products seen in the present day. For the low-mass SMBHs ($M_\mathrm{BH}<10^{8}\ M_{\odot}$) at $0.5<z<3$, on the other hand, variable AGNs show the similar BH-stellar mass ratios with the local objects ($\sim 0.1\%$) but smaller than those observed at $z > 4$. We interpret that host galaxies harboring less massive SMBHs at intermediate redshift have already acquired sufficient stellar mass, although high-z galaxies are still in the early stage of galaxy formation relative to those at the intermediate/local universe. △ Less

Submitted 21 April, 2024; originally announced April 2024.

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

Journal ref: 2024, ApJ, 969, 11

Observational Properties of AGN Obscuration During the Peak of Accretion Growth

Authors: Bovornpratch Vijarnwannaluk, Masayuki Akiyama, Malte Schramm, Yoshihiro Ueda, Yoshiki Matsuoka, Yoshiki Toba, Naoki Matsumoto, Angel Ruiz, Ioannis Georgantopoulos, Ektoras Pouliasis, Elias Koulouridis, Kohei Ichikawa, Marcin Sawicki, Stephen Gwyn

Abstract: We investigated the gas obscuration and host galaxy properties of active galactic nuclei (AGN) during the peak of cosmic accretion growth of supermassive black holes (SMBHs) at redshift 0.8-1.8 using X-ray detected AGN with mid-infrared and far-infrared detection. The sample was classified as type-1 and type-2 AGN using optical spectral and morphological classification while the host galaxy proper… ▽ More We investigated the gas obscuration and host galaxy properties of active galactic nuclei (AGN) during the peak of cosmic accretion growth of supermassive black holes (SMBHs) at redshift 0.8-1.8 using X-ray detected AGN with mid-infrared and far-infrared detection. The sample was classified as type-1 and type-2 AGN using optical spectral and morphological classification while the host galaxy properties were estimated with multiwavelength SED fitting. For type-1 AGN, the black hole mass was determined from MgII emission lines while the black hole mass of type-2 AGN was inferred from the host galaxy's stellar mass. Based on the derived parameters, the distribution of the sample in the absorption hydrogen column density ($N_{\rm H}$) vs. Eddington ratio diagram is examined. Among the type-2 AGN, $28\pm5$\% are in the forbidden zone, where the obscuration by dust torus cannot be maintained due to radiation pressure on dusty material. The fraction is higher than that observed in the local universe from the BAT AGN Spectroscopic Survey (BASS) data release 2 ($11\pm3$\%). The higher fraction implies that the obscuration of the majority of AGN is consistent with the radiation pressure regulated unified model but with an increased incidence of interstellar matter (ISM) obscured AGN. We discuss the possibility of dust-free absorption in type-1 AGN and heavy ISM absorption in type-2 AGN. We also find no statistical difference in the star-formation activity between type-1 and type-2 AGN which may suggest that obscuration triggered by a gas-rich merging is not common among X-ray detected AGN in this epoch. △ Less

Submitted 10 March, 2024; originally announced March 2024.

Comments: Accepted for publication in 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

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

Reconstruction of Cosmic Black Hole Growth and Mass Distribution from Quasar Luminosity Functions at $z>4$: Implications for Faint and Low-mass Populations in JWST

Authors: Wenxiu Li, Kohei Inayoshi, Masafusa Onoue, Wanqiu He, Yoshiki Matsuoka, Zhiwei Pan, Masayuki Akiyama, Takuma Izumi, Tohru Nagao

Abstract: The evolution of the quasar luminosity function (QLF) is fundamental to understanding the cosmic evolution of black holes (BHs) through their accretion phases. In the era of the James Webb Space Telescope (JWST), Euclid, and Nancy Grace Roman Space Telescope, their unprecedented detection sensitivity and wide survey area can unveil the low-luminosity quasar and low-mass BH population, and provide… ▽ More The evolution of the quasar luminosity function (QLF) is fundamental to understanding the cosmic evolution of black holes (BHs) through their accretion phases. In the era of the James Webb Space Telescope (JWST), Euclid, and Nancy Grace Roman Space Telescope, their unprecedented detection sensitivity and wide survey area can unveil the low-luminosity quasar and low-mass BH population, and provide new insights into quasar host galaxies. We present a theoretical model describing BH growth from initial seeding at $z>20$ to $z\sim 4$,incorporating the duration of accretion episodes, the distribution of Eddington ratios, and the mass dependency of BH accretion rates. By constraining the model parameters with the observed QLFs at $4\leq z\leq6$ across a wide UV luminosity range ($-29<M_{\rm 1450}<-24$), we find that the high-redshift BH population grows rapidly at $z\gtrsim6$, and decelerates the pace in subsequent epochs. Toward lower redshifts ($z<6$), mass-dependent accretion inhibits the growth of high-mass BHs with $M_{\bullet}>10^8~M_\odot$, leading to mass saturation at $M_\bullet\gtrsim 10^{10}~M_\odot$. We predict the BH mass function down to $M_{\bullet}\sim 10^6~M_\odot$ for both unobscured and obscured quasar populations at $4\leq z \leq 11$, offering a benchmark for future observational tests. Our model accounts for the presence of both bright and faint quasars at $z>4$, including those discovered by JWST. Furthermore, our findings suggest two distinct pathways for the early assembly of the BH-galaxy mass correlation: the population with a BH-to-stellar mass ratio near the local value of $M_\bullet/M_{\star}\simeq5\times10^{-3}$ maintains a proximity to the relation through its evolution via moderate growth, while the population that begins to grow above the local relation accretes mass rapidly and becomes as overmassive as $M_\bullet/M_\star \sim 0.01-0.1$ by $z\sim 6$. △ Less

Submitted 30 April, 2024; v1 submitted 9 June, 2023; originally announced June 2023.

Comments: 16 pages, 8 figures; Accepted by 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

Supermassive primordial black holes: a view from clustering of quasars at $z \sim 6$

Authors: Takumi Shinohara, Wanqiu He, Yoshiki Matsuoka, Tohru Nagao, Teruaki Suyama, Tomo Takahashi

Abstract: We investigate a scenario where primordial black holes (PBHs) can be the progenitors of supermassive black holes (SMBHs) observed at $z\sim6$. To this end, we carried out clustering analysis using a sample of 81 quasars at $5.88 <z<6.49$, which is constructed in Subaru High-$z$ Exploration of Low-Luminosity Quasars (SHELLQs) project, and 11 quasars in the same redshift range selected from the lite… ▽ More We investigate a scenario where primordial black holes (PBHs) can be the progenitors of supermassive black holes (SMBHs) observed at $z\sim6$. To this end, we carried out clustering analysis using a sample of 81 quasars at $5.88 <z<6.49$, which is constructed in Subaru High-$z$ Exploration of Low-Luminosity Quasars (SHELLQs) project, and 11 quasars in the same redshift range selected from the literature. The resulting angular auto-correlation function (ACF) can be fitted to a power-law form of $ω_θ= 0.045^{+0.114}_{-0.106}~θ^{-0.8}$ over a scale of $0.2\!-\!10$ degrees. We compare the ACF of the quasars to that predicted for the PBH model at $z\sim 6$ and found that such a scenario is excluded for a broad range of parameter space, from which we can conclude that a scenario with PBHs as SMBHs is not viable. We also discuss a model in which SMBHs at $z \sim 6$ originate from the direct collapse of PBH clumps and argue that the observed ACF excludes such a scenario in the context of our PBH model. △ Less

Submitted 17 April, 2023; originally announced April 2023.

Comments: 23 pages, 7 figures

GREX-PLUS Science Book

Authors: GREX-PLUS Science Team, :, Akio K. Inoue, Yuichi Harikane, Takashi Moriya, Hideko Nomura, Shunsuke Baba, Yuka Fujii, Naoteru Gouda, Yasuhiro Hirahara, Yui Kawashima, Tadayuki Kodama, Yusei Koyama, Hiroyuki Kurokawa, Taro Matsuo, Yoshiki Matsuoka, Shuji Matsuura, Ken Mawatari, Toru Misawa, Kentaro Nagamine, Kimihiko Nakajima, Shota Notsu, Takafumi Ootsubo, Kazumasa Ohno, Hideo Sagawa , et al. (7 additional authors not shown)

Abstract: GREX-PLUS (Galaxy Reionization EXplorer and PLanetary Universe Spectrometer) is a mission candidate for a JAXA's strategic L-class mission to be launched in the 2030s. Its primary sciences are two-fold: galaxy formation and evolution and planetary system formation and evolution. The GREX-PLUS spacecraft will carry a 1.2 m primary mirror aperture telescope cooled down to 50 K. The two science instr… ▽ More GREX-PLUS (Galaxy Reionization EXplorer and PLanetary Universe Spectrometer) is a mission candidate for a JAXA's strategic L-class mission to be launched in the 2030s. Its primary sciences are two-fold: galaxy formation and evolution and planetary system formation and evolution. The GREX-PLUS spacecraft will carry a 1.2 m primary mirror aperture telescope cooled down to 50 K. The two science instruments will be onboard: a wide-field camera in the 2-8 $μ$m wavelength band and a high resolution spectrometer with a wavelength resolution of 30,000 in the 10-18 $μ$m band. The GREX-PLUS wide-field camera aims to detect the first generation of galaxies at redshift $z>15$. The GREX-PLUS high resolution spectrometer aims to identify the location of the water ``snow line'' in proto-planetary disks. Both instruments will provide unique data sets for a broad range of scientific topics including galaxy mass assembly, origin of supermassive blackholes, infrared background radiation, molecular spectroscopy in the interstellar medium, transit spectroscopy for exoplanet atmosphere, planetary atmosphere in the Solar system, and so on. △ Less

Submitted 30 May, 2023; v1 submitted 17 April, 2023; originally announced April 2023.

Comments: This document is the first version of a collection of scientific themes which can be achieved with GREX-PLUS. Each section in Chapters 2 and 3 is based on the presentation at the GREX-PLUS Science Workshop held on 24-25 March, 2022 at Waseda University

eROSITA Final Equatorial-Depth Survey (eFEDS): eFEDS X-ray view of WERGS Radio Galaxies selected by the Subaru/HSC and VLA/FIRST survey

Authors: Kohei Ichikawa, Takuji Yamashita, Andrea Merloni, Junyao Li, Teng Liu, Mara Salvato, Masayuki Akiyama, Riccardo Arcodia, Tom Dwelly, Xiaoyang Chen, Masatoshi Imanishi, Kohei Inayoshi, Toshihiro Kawaguchi, Taiki Kawamuro, Mitsuru Kokubo, Yoshiki Matsuoka, Tohru Nagao, Malte Schramm, Hyewon Suh, Masayuki Tanaka, Yoshiki Toba, Yoshihiro Ueda

Abstract: We construct the eROSITA X-ray catalog of radio galaxies discovered by the WERGS survey that is made by the cross-matching of the wide-area Subaru/HSC optical survey and VLA/FIRST 1.4 GHz radio survey. We find 393 eROSITA detected radio galaxies in the 0.5--2 keV band in the eFEDS field covering 140~deg$^2$. Thanks to the wide and medium depth eFEDS X-ray survey, the sample contains the rare and m… ▽ More We construct the eROSITA X-ray catalog of radio galaxies discovered by the WERGS survey that is made by the cross-matching of the wide-area Subaru/HSC optical survey and VLA/FIRST 1.4 GHz radio survey. We find 393 eROSITA detected radio galaxies in the 0.5--2 keV band in the eFEDS field covering 140~deg$^2$. Thanks to the wide and medium depth eFEDS X-ray survey, the sample contains the rare and most X-ray luminous radio galaxies above the knee of the X-ray luminosity function, spanning 44<log L(0.5-2keV,abs)<46.5 at $1<z<4$. Based on the X-ray properties obtained by the spectral fitting, 37 sources show obscured AGN signature with $\log (N_\mathrm{H}/\mathrm{cm}^{-2})>22$. Those obscured and radio AGN reside in $0.4<z<3.2$, indicating that they are obscured counterparts of the radio-loud quasar, which are missed in the previous optical quasar surveys. By combining radio and X-ray luminosities, the jet production efficiency $η_\mathrm{jet}$ is investigated, and we find 14 sources with extremely high jet production efficiency at $η_\mathrm{jet}\approx1$. This high $η_\mathrm{jet}$ value might be a result of 1) the decreased radiation efficiency of $η_\mathrm{rad}<0.1$ due to the low accretion rate for those sources and/or 2) the boosting due to the decline of $L_\mathrm{bol}$ by a factor of 10--100 by keeping $P_\mathrm{jet}$ constant in the previous Myr, indicating the experience of the AGN feedback. Finally, inferring the BH masses from the stellar-mass, we find that X-ray luminous sources show the excess of the radio emission with respect to the value estimated from the fundamental plane. Such radio emission excess cannot be explained by the Doppler booming alone, and therefore disk-jet connection of X-ray luminous eFEDS-WERGS is fundamentally different from the conventional fundamental plane which mainly covers low accretion regime. △ Less

Submitted 15 February, 2023; originally announced February 2023.

Comments: 20 pages, 11 figures, accepted for publication in A&A

Journal ref: A&A 672, A171 (2023)

Extreme nature of four blue-excess dust-obscured galaxies revealed by optical spectroscopy

Authors: Akatoki Noboriguchi, Tohru Nagao, Yoshiki Toba, Kohei Ichikawa, Masaru Kajisawa, Nanako Kato, Toshihiro Kawaguchi, Hideo Matsuhara, Yoshiki Matsuoka, Kyoko Onishi, Masafusa Onoue, Nozomu Tamada, Koki Terao, Yuichi Terashima, Yoshihiro Ueda, Takuji Yamashita

Abstract: We report optical spectroscopic observations of four blue-excess dust-obscured galaxies (BluDOGs) identified by Subaru Hyper Suprime-Cam. BluDOGs are a sub-class of dust-obscured galaxies (DOGs, defined with the extremely red color $(i-[22])_{\rm AB} \geq 7.0$; Toba et al. 2015), showing a significant flux excess in the optical $g$- and $r$-bands over the power-law fits to the fluxes at the longer… ▽ More We report optical spectroscopic observations of four blue-excess dust-obscured galaxies (BluDOGs) identified by Subaru Hyper Suprime-Cam. BluDOGs are a sub-class of dust-obscured galaxies (DOGs, defined with the extremely red color $(i-[22])_{\rm AB} \geq 7.0$; Toba et al. 2015), showing a significant flux excess in the optical $g$- and $r$-bands over the power-law fits to the fluxes at the longer wavelengths. Noboriguchi et al. (2019) has suggested that BluDOGs may correspond to the blowing-out phase involved in a gas-rich major merger scenario. However the detailed properties of BluDOGs are not understood because of the lack of spectroscopic information. In this work, we carry out deep optical spectroscopic observations of four BluDOGs using Subaru/FOCAS and VLT/FORS2. The obtained spectra show broad emission lines with extremely large equivalent widths, and a blue wing in the CIV line profile. The redshifts are between 2.2 and 3.3. The averaged rest-frame equivalent widths of the CIV lines are $160\pm33$ $\mathrm{\mathring{A}}$, $\sim$7 times higher than the average of a typical type-1 quasar. The FWHMs of their velocity profiles are between 1990 and 4470 ${\rm km\ s^{-1}}$, and their asymmetric parameters are 0.05 and 0.25. Such strong CIV lines significantly affect the broad-band magnitudes, which is partly the origin of the blue excess seen in the spectral energy distribution of BluDOGs. Their estimated supermassive black hole masses are $1.1\times10^8 < M_{\rm BH}/M_\odot < 5.5 \times 10^8$. The inferred Eddington ratios of the BluDOGs are higher than 1 ($1.1< λ_{\rm Edd} < 3.8$), suggesting that the BluDOGs are in a rapidly evolving phase of supermassive black holes. △ Less

Submitted 30 November, 2022; originally announced November 2022.

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

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 Obscured Fraction of Quasars at Cosmic Noon

Authors: Bovornpratch Vijarnwannaluk, Masayuki Akiyama, Malte Schramm, Yoshihiro Ueda, Yoshiki Matsuoka, Yoshiki Toba, Marcin Sawicki, Stephen Gwyn, Janek Pflugradt

Abstract: Statistical studies of X-ray selected Active Galactic Nuclei (AGN) indicate that the fraction of obscured AGN increases with increasing redshift, and the results suggest that a significant part of the accretion growth occurs behind obscuring material in the early universe. We investigate the obscured fraction of highly accreting X-ray AGN at around the peak epoch of supermassive black hole growth… ▽ More Statistical studies of X-ray selected Active Galactic Nuclei (AGN) indicate that the fraction of obscured AGN increases with increasing redshift, and the results suggest that a significant part of the accretion growth occurs behind obscuring material in the early universe. We investigate the obscured fraction of highly accreting X-ray AGN at around the peak epoch of supermassive black hole growth utilizing the wide and deep X-ray and optical/IR imaging datasets. A unique sample of luminous X-ray selected AGNs above $z>2$ was constructed by matching the XMM-SERVS X-ray point-source catalog with a PSF-convolved photometric catalog covering from $u^*$ to 4.5$μ\mathrm{m}$ bands. Photometric redshift, hydrogen column density, and 2-10 keV AGN luminosity of the X-ray selected AGN candidates were estimated. Using the sample of 306 2-10 keV detected AGN at above redshift 2, we estimate the fraction of AGN with $\log N_{\rm H}\ (\rm cm^{-2})>22$, assuming parametric X-ray luminosity and absorption functions. The results suggest that $76_{-3}^{+4}\%$ of luminous quasars ($\log L_X\ (\rm erg\ s^{-1}) >44.5$) above redshift 2 are obscured. The fraction indicates an increased contribution of obscured accretion at high redshift than that in the local universe. We discuss the implications of the increasing obscured fraction with increasing redshift based on the AGN obscuration scenarios, which describe obscuration properties in the local universe. Both the obscured and unobscured $z>2$ AGN show a broad range of SEDs and morphology, which may reflect the broad variety of host galaxy properties and physical processes associated with the obscuration. △ Less

Submitted 30 September, 2022; v1 submitted 16 September, 2022; originally announced September 2022.

Comments: 25 Pages, 19 Figure, ApJ accepted. Updated version contains corrections based on comments from the community

A Wide and Deep Exploration of Radio Galaxies with Subaru HSC (WERGS). IX. The Most Overdense Region at z~5 Inhabited by a Massive Radio Galaxy

Authors: Hisakazu Uchiyama, Takuji Yamashita, Tohru Nagao, Yoshiaki Ono, Jun Toshikawa, Kohei Ichikawa, Nozomu Kawakatu, Masaru Kajisawa, Yoshiki Toba, Yoshiki Matsuoka, Mariko Kubo, Masatoshi Imanishi, Kei Ito, Toshihiro Kawaguchi, Chien-Hsiu Lee, Tomoki Saito

Abstract: We report on the galaxy density environment around a high-z radio galaxy (HzRG) at z=4.72, HSC J083913.17+011308.1 (HSC J0839+0113), probed using an r-dropout Lyman break galaxy (LBG) sample from the Hyper Suprime-Cam Subaru Strategic Program data. We find that HSC J0839+0113 resides in the outskirt of an overdense region identified by the r-dropout galaxies at a 4.7 sigma significance level. The… ▽ More We report on the galaxy density environment around a high-z radio galaxy (HzRG) at z=4.72, HSC J083913.17+011308.1 (HSC J0839+0113), probed using an r-dropout Lyman break galaxy (LBG) sample from the Hyper Suprime-Cam Subaru Strategic Program data. We find that HSC J0839+0113 resides in the outskirt of an overdense region identified by the r-dropout galaxies at a 4.7 sigma significance level. The projected distance between HSC J0839+0113 and the peak position of the overdense region is 0.4 physical Mpc which is shorter than the typical protocluster radius in this epoch. According to the extended Press Schechter and the light cone models, the HSC J0839+0113-hosted overdense region is expected to evolve into a halo > 10^14 Msun at z=0 with a high probability of >80 %. These findings suggest that HSC J0839+0113 is associated with a protocluster. The HSC J0839+0113 rich-system is the most overdense region of LBGs among the known protoclusters with LBGs in the same cosmic epoch. △ Less

Submitted 31 August, 2022; originally announced August 2022.

Comments: 7 pages, 3 figures, accepted in PASJ

Optical IFU Observations of GOALS Sample with KOOLS-IFU on Seimei Telescope: Initial results of 9 U/LIRGs at $z <$ 0.04

Authors: Yoshiki Toba, Satoshi Yamada, Kazuya Matsubayashi, Koki Terao, Aoi Moriya, Yoshihiro Ueda, Kouji Ohta, Aoi Hashiguchi, Kazuharu G. Himoto, Hideyuki Izumiura, Kazuma Joh, Nanako Kato, Shuhei Koyama, Hiroyuki Maehara, Rana Misato, Akatoki Noboriguchi, Shoji Ogawa, Naomi Ota, Mio Shibata, Nozomu Tamada, Anri Yanagawa, Naoki Yonekura, Tohru Nagao, Masayuki Akiyama, Masaru Kajisawa , et al. (1 additional authors not shown)

Abstract: We present ionized gas properties of 9 local ultra/luminous infrared galaxies (U/LIRGs) at $z <$ 0.04 through IFU observations with KOOLS-IFU on Seimei Telescope. The observed targets are drawn from the Great Observatories All-sky LIRG Survey (GOALS), covering a wide range of merger stages. We successfully detect emission lines such as H$β$, [OIII]$λ$5007, H$α$, [NII]$λλ$6549,6583, and [SII]$λλ$67… ▽ More We present ionized gas properties of 9 local ultra/luminous infrared galaxies (U/LIRGs) at $z <$ 0.04 through IFU observations with KOOLS-IFU on Seimei Telescope. The observed targets are drawn from the Great Observatories All-sky LIRG Survey (GOALS), covering a wide range of merger stages. We successfully detect emission lines such as H$β$, [OIII]$λ$5007, H$α$, [NII]$λλ$6549,6583, and [SII]$λλ$6717,6731 with a spectral resolution of $R$ = 1500-2000, which provides (i) spatially-resolved ($\sim$200-700 pc) moment map of ionized gas and (ii) diagnostics for active galactic nucleus (AGN) within the central $\sim$3--11 kpc in diameter for our sample. We find that [OIII] outflow that is expected to be driven by AGN tends to be stronger (i) towards the galactic center and (ii) as a sequence of merger stage. In particular, the outflow strength in the late-stage (stage D) mergers is about 1.5 times stronger than that in the early-state (stage B) mergers, which indicates that galaxy mergers could induce AGN-driven outflow and play an important role in the co-evolution of galaxies and supermassive black holes. △ Less

Submitted 24 August, 2022; originally announced August 2022.

Comments: 12 pages, 8 figures, and 2 tables, accepted for publication in PASJ

Evolution of the mass relation between supermassive black holes and dark matter halos across the cosmic time

Authors: Aryan Bansal, Kiyotomo Ichiki, Hiroyuki Tashiro, Yoshiki Matsuoka

Abstract: A positive observational proof suggests that most galaxies contain a central supermassive black hole (SMBH) with mass in the range of $10^6M_\odot$ $-$ $10^{10}M_\odot$. It is suggested that the mass of SMBHs is proportionally related to that of the dark matter (DM) halo, even at $z= 6$. This implies that these SMBHs could coevolve with the host DM halos. Here, we investigate the mass evolution of… ▽ More A positive observational proof suggests that most galaxies contain a central supermassive black hole (SMBH) with mass in the range of $10^6M_\odot$ $-$ $10^{10}M_\odot$. It is suggested that the mass of SMBHs is proportionally related to that of the dark matter (DM) halo, even at $z= 6$. This implies that these SMBHs could coevolve with the host DM halos. Here, we investigate the mass evolution of SMBHs in a hierarchical structure formation by building halo merger trees using the extended Press-Schechter formalism. An SMBH with a mass that follows various power-law relations with the DM halo mass is assigned as an initial condition. Assuming that the mass growth of all black holes is due to mergers, we obtain the relation between SMBH and DM halo masses at the present epoch. By requiring that the mass of the SMBHs at $z=0$ should not be greater than the currently observed SMBH-DM halo mass relation, a lower bound on the mass of the DM halo that can contain a SMBH can be imposed at $z= 6$ as $M_{\rm lim} > 3.6 \times 10^{10} M_\odot\times (1.4-n)^{2.3}$, where $n$ is the power-law index of the SMBH-DM halo mass relation at $z=6$. Because we only consider mergers for the mass evolution of SMBHs, this model is simplistic and should underestimate the mass of SMBHs relative to the mass of the DM halo at the present epoch. Our work aims to constrain the SMBH-DM mass relation at z=6 and not build a new model to explain the observations △ Less

Submitted 25 May, 2023; v1 submitted 3 June, 2022; originally announced June 2022.

Comments: 10 pages, 7 figures

Detection of extended millimeter emission in the host galaxy of 3C273 and its implications for QSO feedback via high dynamic range ALMA imaging

Authors: Shinya Komugi, Yoshiki Toba, Yoshiki Matsuoka, Toshiki Saito, Takuji Yamashita

Abstract: We estimate the amount of negative feedback energy injected into the ISM of the host galaxy of 3C273, a prototypical radio loud quasar. We obtained 93, 233 and 343 GHz continuum images with the Atacama Large Millimeter/Sub-millimeter Array (ALMA). After self calibration and point source subtraction, we reach an image dynamic range of $\sim 85000$ at 93\ GHz, $\sim 39000$ at 233\ GHz and… ▽ More We estimate the amount of negative feedback energy injected into the ISM of the host galaxy of 3C273, a prototypical radio loud quasar. We obtained 93, 233 and 343 GHz continuum images with the Atacama Large Millimeter/Sub-millimeter Array (ALMA). After self calibration and point source subtraction, we reach an image dynamic range of $\sim 85000$ at 93\ GHz, $\sim 39000$ at 233\ GHz and $\sim 2500$ at 343\ GHz. These are currently the highest image dynamic range obtained using ALMA. We detect spatially extended millimeter emission associated with the host galaxy, cospatial with the Extended Emission Line Region (EELR) observed in the optical. The millimeter spectral energy distribution and comparison with centimeter data show that the extended emission cannot be explained by dust thermal emission, synchrotron or thermal bremsstrahlung arising from massive star formation. We interpret the extended millimeter emission as thermal bremsstrahlung from gas directly ionized by the central source. The extended flux indicates that at least $\sim 7\%$ of the bolometric flux of the nuclear source was used to ionize atomic hydrogen in the host galaxy. The ionized gas is estimated to be as massive as $10^{10}$ to $10^{11}\ \mathrm{M_\odot}$, but the molecular gas fraction with respect to the stellar mass is consistent with other ellipticals, suggesting that direct ionization ISM by the QSO may not be sufficient to suppress star formation, or we are witnessing a short timescale before negative feedback becomes observable. The discovery of a radio counterpart to EELRs provides a new pathway to studying the QSO-host ISM interaction. △ Less

Submitted 28 March, 2022; originally announced March 2022.

On the Connection between Supermassive Black Hole and Galaxy Growth in the Reionization Epoch

Authors: Junyao Li, John D. Silverman, Takuma Izumi, Wanqiu He, Masayuki Akiyama, Kohei Inayoshi, Yoshiki Matsuoka, Masafusa Onoue, Yoshiki Toba

Abstract: The correlation between the mass of supermassive black holes (SMBHs; $\mathcal{M}_{\rm BH}$) and their host galaxies ($\mathcal{M}_\star$) in the reionization epoch provides valuable constraints on their early growth. High-redshift quasars typically have a $\mathcal{M}_{\rm BH}$/$\mathcal{M}_\star$ ratio significantly elevated in comparison to the local value. However, the degree to which this app… ▽ More The correlation between the mass of supermassive black holes (SMBHs; $\mathcal{M}_{\rm BH}$) and their host galaxies ($\mathcal{M}_\star$) in the reionization epoch provides valuable constraints on their early growth. High-redshift quasars typically have a $\mathcal{M}_{\rm BH}$/$\mathcal{M}_\star$ ratio significantly elevated in comparison to the local value. However, the degree to which this apparent offset is driven by observational biases is unclear for the most distant quasars. To address this issue, we model the sample selection and measurement biases for a compilation of 20 quasars at $z\sim6$ with host properties based on ALMA observations. We find that the observed distribution of quasars in the $\mathcal{M}_{\rm BH} - \mathcal{M}_\star$ plane can be reproduced by assuming that the underlying SMBH population at $z\sim6$ follows the relationship in the local universe. However, a positive or even a negative evolution in $\mathcal{M}_{\rm BH}$/$\mathcal{M}_\star$ can also explain the data, depending on whether the intrinsic scatter evolves and the strength of various systematic uncertainties. To break these degeneracies, an improvement in the accuracy of mass measurements and an expansion of the current sample to lower $\mathcal{M}_{\rm BH}$ limits are needed. Furthermore, assuming a radiative efficiency of 0.1 and quasar duty cycles estimated from the active SMBH fraction, significant outliers in $\mathcal{M}_{\rm BH}$/$\mathcal{M}_\star$ tend to move toward the local relation given their instantaneous BH growth rate and star formation rate. This may provide evidence for a self-regulated SMBH-galaxy coevolution scenario that is in place at $z\sim6$, with AGN feedback being a possible driver. △ Less

Submitted 29 May, 2022; v1 submitted 20 March, 2022; originally announced March 2022.

Comments: 10 pages, 5 figures, published in ApJL

Multi-line assessment of narrow-line regions in $z \sim$ 3 radio galaxies

Authors: Koki Terao, Tohru Nagao, Kyoko Onishi, Kenta Matsuoka, Masayuki Akiyama, Yoshiki Matsuoka, Takuji Yamashita

Abstract: In this paper, we utilize high-quality rest-UV spectra of three radio galaxies at $z \sim$ 3 observed with VLT/FORS2 to measure the flux of several emission lines including relatively faint ones, such as NIV]$λ$1486, OIII]$λ$1663, and [NeIV]$λ$2424. Additionally, we collect fluxes of faint rest-UV emission lines in 12 $z \sim$ 3 radio galaxies from the literature. Previously, physical and chemical… ▽ More In this paper, we utilize high-quality rest-UV spectra of three radio galaxies at $z \sim$ 3 observed with VLT/FORS2 to measure the flux of several emission lines including relatively faint ones, such as NIV]$λ$1486, OIII]$λ$1663, and [NeIV]$λ$2424. Additionally, we collect fluxes of faint rest-UV emission lines in 12 $z \sim$ 3 radio galaxies from the literature. Previously, physical and chemical properties of narrow-line regions (NLRs) in high-$z$ active galactic nuclei (AGNs) have been investigated mostly by using only strong rest-UV emission-lines (e.g., NV$λ$1240, CIV$λ$1549, HeII$λ$1640, and CIII]$λ$1909). Such strong-line diagnostics are based on various assumptions due to the limitation in the number of available emission-line constraints. In this work, both physical and chemical properties of NLR clouds in each object are estimated by fitting detailed photoionization models to the measured emission-line fluxes. We confirm that the metallicity of NLRs in AGNs at $z \sim$ 3 is solar or super-solar, without assuming the gas density and ionization parameter thanks to the constrains from the faint emission lines. This result suggests that high-$z$ radio galaxies are already chemically matured at $z \sim$ 3. △ Less

Submitted 7 March, 2022; originally announced March 2022.

Comments: 17 pages, 10 figures, accepted for publication in the 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

A Wide and Deep Exploration of Radio Galaxies with Subaru HSC (WERGS). IV. Rapidly Growing (Super-)Massive Black Holes in Extremely Radio-Loud Galaxies

Authors: Kohei Ichikawa, Takuji Yamashita, Yoshiki Toba, Tohru Nagao, Kohei Inayoshi, Maria Charisi, Wanqiu He, Alexander Y. Wagner, Masayuki Akiyama, Bovornpratch Vijarnwannaluk, Xaioyang Chen, Masaru Kajisawa, Taiki Kawamuro, Chien-Hsiu Lee, Yoshiki Matsuoka, Malte Schramm, Hyewon Suh, Masayuki Tanaka, Hisakazu Uchiyama, Yoshihiro Ueda, Janek Pflugradt, Hikaru Fukuchi

Abstract: We present the optical and infrared properties of 39 extremely radio-loud galaxies discovered by cross-matching the Subaru/Hyper Suprime-Cam (HSC) deep optical imaging survey and VLA/FIRST 1.4 GHz radio survey. The recent Subaru/HSC strategic survey revealed optically-faint radio galaxies (RG) down to $g_\mathrm{AB} \sim 26$, opening a new parameter space of extremely radio-loud galaxies (ERGs) wi… ▽ More We present the optical and infrared properties of 39 extremely radio-loud galaxies discovered by cross-matching the Subaru/Hyper Suprime-Cam (HSC) deep optical imaging survey and VLA/FIRST 1.4 GHz radio survey. The recent Subaru/HSC strategic survey revealed optically-faint radio galaxies (RG) down to $g_\mathrm{AB} \sim 26$, opening a new parameter space of extremely radio-loud galaxies (ERGs) with radio-loudness parameter of $\log \mathcal{R}_\mathrm{rest} = \log (f_{1.4 \mathrm{GHz,rest}}/f_{g,\mathrm{rest}}) >4$. Because of their optical faintness and small number density of $\sim1~$deg$^{-2}$, such ERGs were difficult to find in the previous wide but shallow, or deep but small area optical surveys. ERGs show intriguing properties that are different from the conventional RGs: (1) most ERGs reside above or on the star-forming main-sequence, and some of them might be low-mass galaxies with $\log (M_\star/M_\odot) < 10$. (2) ERGs exhibit a high specific black hole accretion rate, reaching the order of the Eddington limit. The intrinsic radio-loudness ($\mathcal{R}_\mathrm{int}$), defined by the ratio of jet power over bolometric radiation luminosity, is one order of magnitude higher than that of radio quasars. This suggests that ERGs harbor a unique type of active galactic nuclei (AGN) that show both powerful radiations and jets. Therefore, ERGs are prominent candidates of very rapidly growing black holes reaching Eddington-limited accretion just before the onset of intensive AGN feedback. △ Less

Submitted 5 August, 2021; originally announced August 2021.

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

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)

The eROSITA Final Equatorial-Depth Survey (eFEDS): The first archetypal Quasar in the feedback phase discovered by eROSITA

Authors: M. Brusa, T. Urrutia, Y. Toba, J. Buchner, J. -Y. Li, T. Liu, M. Perna, M. Salvato, A. Merloni, B. Musiimenta, K. Nandra, J. Wolf, R. Arcodia, T. Dwelly, A. Georgakakis, A. Goulding, Y. Matsuoka, T. Nagao, M. Schramm, J. D. Silverman, Y. Terashima

Abstract: Theoretical models of galaxy-AGN co-evolution ascribe an important role for the feedback process to a short, luminous, obscured, and dust-enshrouded phase during which the accretion rate of the SMBH is expected to be at its maximum and the associated AGN-driven winds are also predicted to be maximally developed. To test this scenario, we have isolated a text-book candidate from the eROSITA Final E… ▽ More Theoretical models of galaxy-AGN co-evolution ascribe an important role for the feedback process to a short, luminous, obscured, and dust-enshrouded phase during which the accretion rate of the SMBH is expected to be at its maximum and the associated AGN-driven winds are also predicted to be maximally developed. To test this scenario, we have isolated a text-book candidate from the eROSITA Final Equatorial-Depth Survey (eFEDS) obtained within the Performance and Verification program of the eROSITA telescope on board Spectrum Roentgen Gamma. From an initial catalog of 246 hard X-ray selected sources matched with the photometric and spectroscopic information available within the eROSITA and Hyper Suprime-Cam consortia, three candidates Quasars in the feedback phase have been isolated applying the diagnostic proposed in Brusa et al. (2015). Only one source (eFEDSU J091157.5+014327) has a spectrum already available (from SDSS-DR16, z=0.603) and it unambiguously shows the presence of a broad component (FWHM~1650 km/s) in the [OIII]5007 line. The associated observed L_[OIII] is ~2.6x10^{42} erg/s, one to two orders of magnitude larger than that observed in local Seyferts and comparable to those observed in a sample of z~0.5 Type 1 Quasars. From the multiwavelength data available we derive an Eddington Ratio (L_bol/L_Edd) of ~0.25, and a bolometric correction in the hard X-ray of k_bol~10, lower than those observed for objects at similar bolometric luminosity. The presence of an outflow, the high X-ray luminosity and moderate X-ray obscuration (L_X~10^44.8 erg/s, N_H~2.7x10^22 cm^-2) and the red optical color, all match the prediction of quasars in the feedback phase from merger driven models. Forecasting to the full eROSITA all-sky survey with its spectroscopic follow-up, we predict that by the end of 2024 we will have a sample of few hundreds such objects at z=0.5-2. △ Less

Submitted 28 June, 2021; originally announced June 2021.

Comments: 11 pages, 6 figures, 2 tables. Resubmitted to A&A for the Special Issue: "The Early Data Release of eROSITA and Mikhail Pavlinsky ART-XC on the SRG Mission"

Journal ref: A&A 661, A9 (2022)

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

The Sizes of Quasar Host Galaxies with the Hyper Suprime-Cam Subaru Strategic Program

Authors: Junyao Li, John D. Silverman, Xuheng Ding, Michael A. Strauss, Andy Goulding, Simon Birrer, Hassen M. Yesuf, Yongquan Xue, Lalitwadee Kawinwanichakij, Yoshiki Matsuoka, Yoshiki Toba, Tohru Nagao, Malte Schramm, Kohei Inayoshi

Abstract: The relationship between quasars and their host galaxies provides clues on how supermassive black holes (SMBHs) and massive galaxies are jointly assembled. To elucidate this connection, we measure the structural and photometric properties of the host galaxies of ~5000 SDSS quasars at 0.2<z<1 using five-band (grizy) optical imaging from the Hyper Suprime-Cam Subaru Strategic Program. An automated a… ▽ More The relationship between quasars and their host galaxies provides clues on how supermassive black holes (SMBHs) and massive galaxies are jointly assembled. To elucidate this connection, we measure the structural and photometric properties of the host galaxies of ~5000 SDSS quasars at 0.2<z<1 using five-band (grizy) optical imaging from the Hyper Suprime-Cam Subaru Strategic Program. An automated analysis tool is used to forward-model the blended emission of the quasar as characterized by the point spread function and the underlying host galaxy as a two-dimensional Sersic profile. In agreement with previous studies, quasars are preferentially hosted by massive star-forming galaxies with disk-like light profiles. Furthermore, we find that the size distribution of quasar hosts is broad at a given stellar mass and the average values exhibit a size-stellar mass relation as seen with inactive galaxies. In contrast, the sizes of quasar hosts are more compact than inactive star-forming galaxies on average, but not as compact as quiescent galaxies of similar stellar masses. This is true irrespective of quasar properties including bolometric luminosity, Eddington ratio, and black hole mass. These results are consistent with a scenario in which galaxies are concurrently fueling a SMBH and building their stellar bulge from a centrally-concentrated gas reservoir. Alternatively, quasar hosts may be experiencing a compaction process in which stars from the disk and inflowing gas are responsible for growing the bulge. In addition, we confirm that the host galaxies of type-1 quasars have a bias of being closer towards face-on, suggesting that galactic-scale dust can contribute to obscuring the broad-line region. △ Less

Submitted 13 May, 2021; originally announced May 2021.

Comments: 34 pages, 22 figures, 3 tables, resubmitted to ApJ after addressing referee's comments; catalog will be available after the paper is published

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

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

Constraining the quasar radio-loud fraction at $z \sim 6$ with deep radio observations

Authors: Yuanqi Liu, Ran Wang, Emmanuel Momjian, Eduardo Banados, Greg Zeimann, Chris J. Willott, Yoshiki Matsuoka, Alain Omont, Yali Shao, Qiong Li, Jianan Li

Abstract: We carry out a series of deep Karl G. Jansky Very Large Array (VLA) S-band observations of a sample of 21 quasars at $z\sim6$. The new observations expand the searches of radio continuum emission to the optically faint quasar population at the highest redshift with rest-frame $4400 \rm Å$ luminosities down to $3 \times10^{11} \ L_{\odot}$. We report the detections of two new radio-loud quasars: CF… ▽ More We carry out a series of deep Karl G. Jansky Very Large Array (VLA) S-band observations of a sample of 21 quasars at $z\sim6$. The new observations expand the searches of radio continuum emission to the optically faint quasar population at the highest redshift with rest-frame $4400 \rm Å$ luminosities down to $3 \times10^{11} \ L_{\odot}$. We report the detections of two new radio-loud quasars: CFHQS J2242+0334 (hereafter J2242+0334) at $z=5.88$ and CFHQS J0227$-$0605 (hereafter J0227$-$0605) at $z=6.20$, detected with 3 GHz flux densities of $87.0 \pm 6.3 \ μ\rm Jy$ and $55.4 \pm 6.7 \ μ\rm Jy$, respectively. Their radio \replaced{loudness}{loudnesses} are estimated to be $54.9 \pm 4.7$ and $16.5 \pm 3.2$, respectively. To better constrain the radio-loud fraction (RLF), we combine the new measurements with the archival VLA L-band data as well as available data from the literature, considering the upper limits for non-detections and \deleted{and} possible selection effects. The final derived RLF is $9.4 \pm 5.7\%$ for the optically selected quasars at $z\sim6$. We also compare the RLF to that of the quasar samples at low redshift and check the RLF in different quasar luminosity bins. The RLF for the optically faint objects is still poorly constrained due to the limited sample size. Our \replaced{result}{results} show no evidence of significant quasar RLF evolution with redshift. There is also no clear trend of RLF evolution with quasar UV/optical luminosity due to the limited sample size of optically faint objects with deep radio observations. △ Less

Submitted 14 December, 2020; originally announced December 2020.

Comments: 14 pages, 5 figures, accepted for publication in the Astrophysical Journal

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

Follow-up Observations for IceCube-170922A: Detection of Rapid Near-Infrared Variability and Intensive Monitoring of TXS 0506+056

Authors: Tomoki Morokuma, Yousuke Utsumi, Kouji Ohta, Masayuki Yamanaka, Koji S. Kawabata, Yoshiyuki Inoue, Masaomi Tanaka, Michitoshi Yoshida, Ryosuke Itoh, Mahito Sasada, Nozomu Tominaga, Hiroki Mori, Miho Kawabata, Tatsuya Nakaoka, Maiko Chogi, Taisei Abe, Ruochen Huang, Naoki Kawahara, Hiroki Kimura, Hiroki Nagashima, Kengo Takagi, Yuina Yamazaki, Wei Liu, Ryou Ohsawa, Shigeyuki Sako , et al. (13 additional authors not shown)

Abstract: We present our follow-up observations to search for an electromagnetic counterpart of the IceCube high-energy neutrino, IceCube-170922A. Monitoring observations of a likely counterpart, TXS 0506+056, are also described. First, we quickly took optical and near-infrared images of 7 flat-spectrum radio sources within the IceCube error region right after the neutrino detection and found a rapid flux d… ▽ More We present our follow-up observations to search for an electromagnetic counterpart of the IceCube high-energy neutrino, IceCube-170922A. Monitoring observations of a likely counterpart, TXS 0506+056, are also described. First, we quickly took optical and near-infrared images of 7 flat-spectrum radio sources within the IceCube error region right after the neutrino detection and found a rapid flux decline of TXS 0506+056 in Kanata/HONIR J-band data. Motivated by this discovery, intensive follow-up observations of TXS 0506+056 are continuously done, including our monitoring imaging observations, spectroscopic observations, and polarimetric observations in optical and near-infrared wavelengths. TXS 0506+056 shows a large amplitude (~1.0 mag) variability in a time scale of several days or longer, although no significant variability is detected in a time scale of a day or shorter. TXS 0506+056 also shows a bluer-when-brighter trend in optical and near-infrared wavelengths. Structure functions of variabilities are examined and indicate that TXS 0506+056 is not a special blazar in terms of optical variability. Polarization measurement results of TXS 0506+056 are also discussed. △ Less

Submitted 10 November, 2020; originally announced November 2020.

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

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

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

Subaru Hyper Suprime-Cam View of Quasar Host Galaxies at z < 1

Authors: Toru Ishino, Yoshiki Matsuoka, Shuhei Koyama, Yuya Saeda, Michael A. Strauss, Andy D. Goulding, Masatoshi Imanishi, Toshihiro Kawaguchi, Takeo Minezaki, Tohru Nagao, Akatoki Noboriguchi, Malte Schramm, John D. Silverman, Yoshiaki Taniguchi, Yoshiki Toba

Abstract: Active galactic nuclei (AGNs) are key for understanding the coevolution of galaxies and supermassive black holes (SMBHs). AGN activity is thought to affect the properties of their host galaxies, via a process called "AGN feedback", which drives the co-evolution. From a parent sample of 1151 z < 1 type-1 quasars from the Sloan Digital Sky Survey quasar catalog, we detected host galaxies of 862 of t… ▽ More Active galactic nuclei (AGNs) are key for understanding the coevolution of galaxies and supermassive black holes (SMBHs). AGN activity is thought to affect the properties of their host galaxies, via a process called "AGN feedback", which drives the co-evolution. From a parent sample of 1151 z < 1 type-1 quasars from the Sloan Digital Sky Survey quasar catalog, we detected host galaxies of 862 of them in the high-quality grizy images of the Subaru Hyper Suprime-Cam (HSC) survey. The unprecedented combination of the survey area and depth allows us to perform a statistical analysis of the quasar host galaxies, with small sample variance. We fit the radial image profile of each quasar as a linear combination of the point spread function and the Sersic function, decomposing the images into the quasar nucleus and the host galaxy components. We found that the host galaxies are massive, with stellar mass Mstar > 10^(10) Msun, and are mainly located on the green valley. This trend is consistent with a scenario in which star formation of the host galaxies is suppressed by AGN feedback, that is, AGN activity may be responsible for the transition of these galaxies from the blue cloud to the red sequence. We also investigated the SMBH mass to stellar mass relation of the z < 1 quasars, and found a consistent slope with the local relation, while the SMBHs may be slightly undermassive. However, the above results are subject to our sample selection, which biases against host galaxies with low masses and/or large quasar-to-host flux ratios. △ Less

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

Wide and Deep Exploration of Radio Galaxies with Subaru HSC (WERGS). III. Discovery of a z = 4.72 Radio Galaxy with Lyman Break Technique

Authors: Takuji Yamashita, Tohru Nagao, Hiroyuki Ikeda, Yoshiki Toba, Masaru Kajisawa, Yoshiaki Ono, Masayuki Tanaka, Masayuki Akiyama, Yuichi Harikane, Kohei Ichikawa, Toshihiro Kawaguchi, Taiki Kawamuro, Kotaro Kohno, Chien-Hsiu Lee, Kianhong Lee, Yoshiki Matsuoka, Mana Niida, Kazuyuki Ogura, Masafusa Onoue, Hisakazu Uchiyama

Abstract: We report a discovery of $z = 4.72$ radio galaxy, HSC J083913.17+011308.1, by using the Lyman break technique with the Hyper Suprime-Cam Subaru Strategic Survey (HSC-SSP) catalog for VLA FIRST radio sources. The number of known high-$z$ radio galaxies (HzRGs) at $z > 3$ is quite small to constrain the evolution of HzRGs so far. The deep and wide-area optical survey by HSC-SSP enables us to apply t… ▽ More We report a discovery of $z = 4.72$ radio galaxy, HSC J083913.17+011308.1, by using the Lyman break technique with the Hyper Suprime-Cam Subaru Strategic Survey (HSC-SSP) catalog for VLA FIRST radio sources. The number of known high-$z$ radio galaxies (HzRGs) at $z > 3$ is quite small to constrain the evolution of HzRGs so far. The deep and wide-area optical survey by HSC-SSP enables us to apply the Lyman break technique to a large search for HzRGs. For an HzRG candidate among pre-selected $r$-band dropouts with a radio detection, a follow-up optical spectroscopy with GMOS/Gemini has been performed. The obtained spectrum presents a clear Ly$α$ emission line redshifted to $z=4.72$. The SED fitting analysis with the rest-frame UV and optical photometries suggests the massive nature of this HzRG with $\log{M_*/M_{\odot}} = 11.4$. The small equivalent width of Ly$α$ and the moderately red UV colors indicate its dusty host galaxy, implying a chemically evolved and dusty system. The radio spectral index does not meet a criterion for an ultra-steep spectrum: $α^{325}_{1400}$ of $-1.1$ and $α^{150}_{1400}$ of $-0.9$, demonstrating that the HSC-SSP survey compensates for a sub-population of HzRGs which are missed in surveys focusing on an ultra-steep spectral index. △ Less

Submitted 8 June, 2020; originally announced June 2020.

Comments: 10 pages, 5 figures, accepted for publication in AJ

The Mean Absorption Line Spectra of a Selection of Luminous z~6 Lyman Break Galaxies

Authors: Yuichi Harikane, Nicolas Laporte, Richard S. Ellis, Yoshiki Matsuoka

Abstract: We examine the absorption line spectra of a sample of 31 luminous (M_UV=-23) Lyman break galaxies at redshift z~6 using data taken with the FOCAS and OSIRIS spectrographs on the Subaru and GTC telescopes. For two of these sources we present longer exposure data taken at higher spectral resolution from ESO's X-shooter spectrograph. Using these data, we demonstrate the practicality of stacking our l… ▽ More We examine the absorption line spectra of a sample of 31 luminous (M_UV=-23) Lyman break galaxies at redshift z~6 using data taken with the FOCAS and OSIRIS spectrographs on the Subaru and GTC telescopes. For two of these sources we present longer exposure data taken at higher spectral resolution from ESO's X-shooter spectrograph. Using these data, we demonstrate the practicality of stacking our lower resolution data to measure the depth of various interstellar and stellar absorption lines to probe the covering fraction of low ionization gas and the gas-phase and stellar metallicities near the end of the era of cosmic reionization. From maximum absorption line depths of SiII1260 and CII1334, we infer a mean covering fraction of >0.85+/-0.16 for our sample. This is larger than that determined using similar methods for lower luminosity galaxies at slightly lower redshifts, suggesting that the most luminous galaxies appear to have a lower escape fraction than fainter galaxies, and therefore may not play a prominent role in concluding reionization. Using various interstellar absorption lines we deduce gas-phase metallicities close to solar indicative of substantial early enrichment. Using selected stellar absorption lines, we model our spectra with a range of metallicities using techniques successfully employed at lower redshift and deduce a stellar metallicity of 0.4 +0.3/-0.1 solar, consistent with the stellar mass - stellar metallicity relation recently found at z~3-5. We discuss the implications of these metallicity estimates for the typical ages of our luminous galaxies and conclude our results imply initial star formation at redshifts z~10, consistent with independent analyses of earlier objects. △ Less

Submitted 7 October, 2020; v1 submitted 22 May, 2020; originally announced May 2020.

Comments: 13 pages, 10 figures, 3 tables, Accepted for Publication in ApJ

Properties of AGN Multiband Optical Variability in the HSC SSP Transient Survey

Authors: Yuki Kimura, Toru Yamada, Mitsuru Kokubo, Naoki Yasuda, Tomoki Morokuma, Tohru Nagao, Yoshiki Matsuoka

Abstract: We study variability of active galactic nuclei (AGNs) by using the deep optical multiband photometry data obtained from the Hyper Suprime-Cam Subaru Strategic Program (HSC SSP) survey in the COSMOS field. The images analyzed here were taken with 8, 10, 13, and 15 epochs over three years in the $g$, $r$, $i$, and $z$ bands, respectively. We identified 491 robust variable AGN candidates, down to… ▽ More We study variability of active galactic nuclei (AGNs) by using the deep optical multiband photometry data obtained from the Hyper Suprime-Cam Subaru Strategic Program (HSC SSP) survey in the COSMOS field. The images analyzed here were taken with 8, 10, 13, and 15 epochs over three years in the $g$, $r$, $i$, and $z$ bands, respectively. We identified 491 robust variable AGN candidates, down to $i=25$ mag and with redshift up to $4.26$. Ninety percent of the variability-selected AGNs are individually identified with the X-ray sources detected in the Chandra COSMOS Legacy survey. We investigate their properties in variability by using structure function analysis and find that the structure function for low-luminosity AGNs ($L_{\mathrm{bol}}\lesssim10^{45}$ erg s$^{-1}$) shows a positive correlation with luminosity, which is the opposite trend for the luminous quasars. This trend is likely to be caused by larger contribution of the host galaxy light for lower-luminosity AGNs. Using the model templates of galaxy spectra, we evaluate the amount of host galaxy contribution to the structure function analysis and find that dominance of the young stellar population is needed to explain the observed luminosity dependence. This suggests that low-luminosity AGNs at $0.8\lesssim z\lesssim1.8$ are predominantly hosted in star-forming galaxies. The X-ray stacking analysis reveals the significant emission from the individually X-ray undetected AGNs in our variability-selected sample. The stacked samples show very large hardness ratios in their stacked X-ray spectrum, which suggests that these optically variable sources have large soft X-ray absorption by dust-free gas. △ Less

Submitted 19 April, 2020; originally announced April 2020.

Comments: 27 pages, 22 figures. Accepted for publication in ApJ

Large Population of ALMA Galaxies at z>6 with Very High [OIII]88um to [CII]158um Flux Ratios: Evidence of Extremely High Ionization Parameter or PDR Deficit?

Authors: Yuichi Harikane, Masami Ouchi, Akio K. Inoue, Yoshiki Matsuoka, Yoichi Tamura, Tom Bakx, Seiji Fujimoto, Kana Moriwaki, Yoshiaki Ono, Tohru Nagao, Ken-ichi Tadaki, Takashi Kojima, Takatoshi Shibuya, Eiichi Egami, Andrea Ferrara, Simona Gallerani, Takuya Hashimoto, Kotaro Kohno, Yuichi Matsuda, Hiroshi Matsuo, Andrea Pallottini, Yuma Sugahara, Livia Vallini

Abstract: We present our new ALMA observations targeting [OIII]88um, [CII]158um, [NII]122um, and dust continuum emission for three Lyman break galaxies at z=6.0293-6.2037 identified in the Subaru/Hyper Suprime-Cam survey. We clearly detect [OIII] and [CII] lines from all of the galaxies at 4.3-11.8sigma levels, and identify multi-band dust continuum emission in two of the three galaxies, allowing us to esti… ▽ More We present our new ALMA observations targeting [OIII]88um, [CII]158um, [NII]122um, and dust continuum emission for three Lyman break galaxies at z=6.0293-6.2037 identified in the Subaru/Hyper Suprime-Cam survey. We clearly detect [OIII] and [CII] lines from all of the galaxies at 4.3-11.8sigma levels, and identify multi-band dust continuum emission in two of the three galaxies, allowing us to estimate infrared luminosities and dust temperatures simultaneously. In conjunction with previous ALMA observations for six galaxies at z>6, we confirm that all the nine z=6-9 galaxies have high [OIII]/[CII] ratios of L[OIII]/L[CII]~3-20, ~10 times higher than z~0 galaxies. We also find a positive correlation between the [OIII]/[CII] ratio and the Lya equivalent width (EW) at the ~90% confidence level. We carefully investigate physical origins of the high [OIII]/[CII] ratios at z=6-9 using Cloudy, and find that high density of the interstellar medium, low C/O abundance ratio, and the cosmic microwave background attenuation are responsible to only a part of the z=6-9 galaxies. Instead, the observed high [OIII]/[CII] ratios are explained by 10-100 times higher ionization parameters or low photodissociation region (PDR) covering fractions of 0-10%, both of which are consistent with our [NII] observations. The latter scenario can be reproduced with a density bounded nebula with PDR deficit, which would enhance the Lya, Lyman continuum, and C+ ionizing photons escape from galaxies, consistent with the [OIII]/[CII]-Lya EW correlation we find. △ Less

Submitted 24 May, 2020; v1 submitted 24 October, 2019; originally announced October 2019.

Comments: 20 pages, 18 figures, Accepted for publication in ApJ

Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs). X. Discovery of 35 Quasars and Luminous Galaxies at 5.7 $\le$ z $\le$ 7.0

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

Abstract: We report the discovery of 28 quasars and 7 luminous galaxies at 5.7 $\le$ z $\le$ 7.0. This is the tenth in a series of papers from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The total number of spectroscopically identified objects in SHELLQs has… ▽ More We report the discovery of 28 quasars and 7 luminous galaxies at 5.7 $\le$ z $\le$ 7.0. This is the tenth in a series of papers from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The total number of spectroscopically identified objects in SHELLQs has now grown to 93 high-z quasars, 31 high-z luminous galaxies, 16 [O III] emitters at z ~ 0.8, and 65 Galactic cool dwarfs (low-mass stars and brown dwarfs). These objects were found over 900 deg2, surveyed by HSC between 2014 March and 2018 January. The full quasar sample includes 18 objects with very strong and narrow Ly alpha emission, whose stacked spectrum is clearly different from that of other quasars or galaxies. While the stacked spectrum shows N V 1240 emission and resembles that of lower-z narrow-line quasars, the small Ly alpha width may suggest a significant contribution from the host galaxies. Thus these objects may be composites of quasars and star-forming galaxies. △ Less

Submitted 19 August, 2019; originally announced August 2019.

Comments: ApJ, in press. arXiv admin note: text overlap with arXiv:1803.01861

Journal ref: ApJ 883:183, 2019

SILVERRUSH. IX. Lya Intensity Mapping with Star-Forming Galaxies at z=5.7 and 6.6: A Possible Detection of Extended Lya Emission at $\gtrsim$100 comoving kpc around and beyond the Virial-Radius Scale of Galaxy Dark Matter Halos

Authors: Ryota Kakuma, Masami Ouchi, Yuichi Harikane, Yoshiaki Ono, Akio K. Inoue, Yutaka Komiyama, Haruka Kusakabe, Chien-Hsiu Lee, Yuichi Matsuda, Yoshiki Matsuoka, Ken Mawatari, Rieko Momose, Takatoshi Shibuya, Yoshiaki Taniguchi

Abstract: We present results of the cross-correlation Ly$α$ intensity mapping with Subaru/Hyper Suprime-Cam (HSC) ultra-deep narrowband images and Ly$α$ emitters (LAEs) at $z=5.7$ and $6.6$ in a total area of $4$ deg$^2$. Although overwhelming amount of data quality controls have been performed for the narrowband images, we further conduct extensive analysis evaluating systematics of large-scale point-sprea… ▽ More We present results of the cross-correlation Ly$α$ intensity mapping with Subaru/Hyper Suprime-Cam (HSC) ultra-deep narrowband images and Ly$α$ emitters (LAEs) at $z=5.7$ and $6.6$ in a total area of $4$ deg$^2$. Although overwhelming amount of data quality controls have been performed for the narrowband images, we further conduct extensive analysis evaluating systematics of large-scale point-spread-function wings, sky subtractions, and unknown errors based on physically uncorrelated signals and sources found in real HSC images and object catalogs, respectively. Removing the systematics, we carefully calculate cross-correlations between Ly$α$ intensity of the narrowband images and the LAEs. We tentatively identify very diffuse Ly$α$ emission with the $\simeq 3σ$ ($\simeq 2σ$) significance at $\gtrsim$ 100 comoving kpc (ckpc) far from the LAEs at $z=5.7$ ($6.6$), around and probably even beyond a virial radius of star-forming galaxies with $M_\mathrm{h}\sim10^{11}M_\odot$. The diffuse Ly$α$ emission possibly extends up to $1$,$000$ ckpc with the surface brightness of $10^{-20}$-$10^{-19}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$. We confirm that the small-scale ($<150$ ckpc) Ly$α$ radial profiles of LAEs are consistent with those obtained by recent MUSE observations. Comparisons with numerical simulations suggest that the large-scale ($\sim150$-$1$,$000$ ckpc) Ly$α$ emission are not explained by unresolved faint neighboring galaxies including satellites, but by a combination of Ly$α$ photons emitted from the central LAE and other unknown sources, such as cold-gas streams and galactic outflow. We find no evolution in the Ly$α$ radial profiles of our LAEs from $z=5.7$ to $6.6$, where theoretical models predict a flattening of the profile slope made by cosmic reionization, albeit with our moderately large observational errors. △ Less

Submitted 8 June, 2021; v1 submitted 1 June, 2019; originally announced June 2019.

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

Second Data Release of the Hyper Suprime-Cam Subaru Strategic Program

Authors: Hiroaki Aihara, Yusra AlSayyad, Makoto Ando, Robert Armstrong, James Bosch, Eiichi Egami, Hisanori Furusawa, Junko Furusawa, Andy Goulding, Yuichi Harikane, Chiaki Hikage, Paul T. P. Ho, Bau-Ching Hsieh, Song Huang, Hiroyuki Ikeda, Masatoshi Imanishi, Kei Ito, Ikuru Iwata, Anton T. Jaelani, Ryota Kakuma, Kojiro Kawana, Satoshi Kikuta, Umi Kobayashi, Michitaro Koike, Yutaka Komiyama , et al. (40 additional authors not shown)

Abstract: This paper presents the second data release of the Hyper Suprime-Cam Subaru Strategic Program, a wide-field optical imaging survey on the 8.2 meter Subaru Telescope. The release includes data from 174 nights of observation through January 2018. The Wide layer data cover about 300 deg^2 in all five broadband filters (grizy) to the nominal survey exposure (10min in gr and 20min in izy). Partially ob… ▽ More This paper presents the second data release of the Hyper Suprime-Cam Subaru Strategic Program, a wide-field optical imaging survey on the 8.2 meter Subaru Telescope. The release includes data from 174 nights of observation through January 2018. The Wide layer data cover about 300 deg^2 in all five broadband filters (grizy) to the nominal survey exposure (10min in gr and 20min in izy). Partially observed areas are also included in the release; about 1100 deg^2 is observed in at least one filter and one exposure. The median seeing in the i-band is 0.6 arcsec, demonstrating the superb image quality of the survey. The Deep (26 deg^2) and UltraDeep (4 deg^2) data are jointly processed and the UltraDeep-COSMOS field reaches an unprecedented depth of i~28 at 5 sigma for point sources. In addition to the broad-bands, narrow-band data are also available in the Deep and UltraDeep fields. This release includes a major update to the processing pipeline, including improved sky subtraction, PSF modeling, object detection, and artifact rejection. The overall data quality has been improved, but this release is not without problems; there is a persistent deblender problem as well as new issues with masks around bright stars. The user is encouraged to review the issue list before utilizing the data for scientific explorations. All the image products as well as catalog products are available for download. The catalogs are also loaded to a database, which provides an easy interface for users to retrieve data for objects of interest. In addition to these main data products, detailed galaxy shape measurements withheld from the Public Data Release 1 (PDR1) are now available to the community. The shape catalog is drawn from the S16A internal release, which has a larger area than PDR1 (160 deg^2). All products are available at the data release site, https://hsc-release.mtk.nao.ac.jp/. △ Less

Submitted 22 August, 2019; v1 submitted 29 May, 2019; originally announced May 2019.

Comments: 26 pages, 23 figures, accepted for publication in PASJ. Data available at https://hsc-release.mtk.nao.ac.jp/

A Wide and Deep Exploration of Radio Galaxies with Subaru HSC (WERGS). II. Physical Properties derived from the SED Fitting with Optical, Infrared, and Radio Data

Authors: Yoshiki Toba, Takuji Yamashita, Tohru Nagao, Wei-Hao Wang, Yoshihiro Ueda, Kohei Ichikawa, Toshihiro Kawaguchi, Masayuki Akiyama, Bau-Ching Hsieh, Masaru Kajisawa, Chien-Hsiu Lee, Yoshiki Matsuoka, Akatoki Noboriguchi, Masafusa Onoue, Malte Schramm, Masayuki Tanaka, Yutaka Komiyama

Abstract: We present physical properties of radio galaxies (RGs) with $f_{\rm 1.4 GHz} >$ 1 mJy discovered by Subaru Hyper Supreme-Cam (HSC) and VLA Faint Images of the Radio Sky at Twenty-Centimeters (FIRST) survey. For 1056 FIRST RGs at $0 < z \leq 1.7$ with HSC counterparts in about 100 deg$^2$, we compiled multi-wavelength data of optical, near-infrared (IR), mid-IR, far-IR, and radio (150 MHz). We deri… ▽ More We present physical properties of radio galaxies (RGs) with $f_{\rm 1.4 GHz} >$ 1 mJy discovered by Subaru Hyper Supreme-Cam (HSC) and VLA Faint Images of the Radio Sky at Twenty-Centimeters (FIRST) survey. For 1056 FIRST RGs at $0 < z \leq 1.7$ with HSC counterparts in about 100 deg$^2$, we compiled multi-wavelength data of optical, near-infrared (IR), mid-IR, far-IR, and radio (150 MHz). We derived their color excess ($E (B-V)_{*}$), stellar mass, star formation rate (SFR), IR luminosity, the ratio of IR and radio luminosity ($q_{\rm IR}$), and radio spectral index ($α_{\rm radio}$) that are derived from the SED fitting with CIGALE. We also estimated Eddington ratio based on stellar mass and integration of the best-fit SEDs of AGN component. We found that $E (B-V)_{*}$, SFR, and IR luminosity clearly depend on redshift while stellar mass, $q_{\rm IR}$, and $α_{\rm radio}$ do not significantly depend on redshift. Since optically-faint ($i_{\rm AB} \geq 21.3$) RGs that are newly discovered by our RG survey tend to be high redshift, they tend to not only have a large dust extinction and low stellar mass but also have high SFR and AGN luminosity, high IR luminosity, and high Eddington ratio compared to optically-bright ones. The physical properties of a fraction of RGs in our sample seem to differ from a classical view of RGs with massive stellar mass, low SFR, and low Eddington ratio, demonstrating that our RG survey with HSC and FIRST provides us curious RGs among entire RG population. △ Less

Submitted 3 May, 2019; originally announced May 2019.

Comments: 30 pages, 20 figures, and 4 tables, accepted for publication in ApJS. The catalog and SED template of radio galaxies will be accessible through an online service