Showing 1–50 of 164 results for author: Toba, Y

Ongoing and fossil large-scale outflows detected in a high-redshift radio galaxy: [C II] observations of TN J0924$-$2201 at $z=5.174$

Authors: Kianhong Lee, Masayuki Akiyama, Kotaro Kohno, Daisuke Iono, Masatoshi Imanishi, Bunyo Hatsukade, Hideki Umehata, Tohru Nagao, Yoshiki Toba, Xiaoyang Chen, Fumi Egusa, Kohei Ichikawa, Takuma Izumi, Naoki Matsumoto, Malte Schramm, Kenta Matsuoka

Abstract: We present Atacama Large Millimeter/submillimeter Array observations of the [C II] 158 $μ$m line and the underlying continuum emission of TN J0924$-$2201, which is one of the most distant known radio galaxies at $z>5$. The [C II] line and 1-mm continuum emission are detected at the host galaxy. The systemic redshift derived from the [C II] line is $z_{\rm [C II]}=5.1736\pm0.0002$, indicating that… ▽ More We present Atacama Large Millimeter/submillimeter Array observations of the [C II] 158 $μ$m line and the underlying continuum emission of TN J0924$-$2201, which is one of the most distant known radio galaxies at $z>5$. The [C II] line and 1-mm continuum emission are detected at the host galaxy. The systemic redshift derived from the [C II] line is $z_{\rm [C II]}=5.1736\pm0.0002$, indicating that the Ly$α$ line is redshifted by a velocity of $1035\pm10$ km s$^{-1}$, marking the largest velocity offset between the [C II] and Ly$α$ lines recorded at $z>5$ to date. In the central region of the host galaxy, we identified a redshifted substructure of [C II] with a velocity of $702\pm17$ km s$^{-1}$, which is close to the CIV line with a velocity of $500\pm10$ km s$^{-1}$. The position and the velocity offsets align with a model of an outflowing shell structure, consistent with the large velocity offset of Ly$α$. The non-detection of [C II] and dust emission from the three CO(1--0)-detected companions indicates their different nature compared to dwarf galaxies based on the photodissociation region model. Given their large velocity of $\sim1500$ km s$^{-1}$, outflowing molecular clouds induced by the AGN is the most plausible interpretation, and they may exceed the escape velocity of a $10^{13}\,M_{\odot}$ halo. These results suggest that TN J0924$-$2201, with the ongoing and fossil large-scale outflows, is in a distinctive phase of removing molecular gas from a central massive galaxy in an overdense region in the early universe. A dusty HI absorber at the host galaxy is an alternative interpretation. △ Less

Submitted 1 July, 2024; originally announced July 2024.

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

SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES). V. Confusion-limited Submillimeter Galaxy Number Counts at 450 $μ$m and Data Release for the COSMOS Field

Authors: Zhen-Kai Gao, Chen-Fatt Lim, Wei-Hao Wang, Chian-Chou Chen, Ian Smail, Scott C. Chapman, Xian Zhong Zheng, Hyunjin Shim, Tadayuki Kodama, Yiping Ao, Siou-Yu Chang, David L. Clements, James S. Dunlop, Luis C. Ho, Yun-Hsin Hsu, Chorng-Yuan Hwang, Ho Seong Hwang, M. P. Koprowski, Douglas Scott, Stephen Serjeant, Yoshiki Toba, Sheona A. Urquhart

Abstract: We present confusion-limited SCUBA-2 450-$μ$m observations in the COSMOS-CANDELS region as part of the JCMT Large Program, SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES). Our maps at 450 and 850 $μ$m cover an area of 450 arcmin$^2$. We achieved instrumental noise levels of $σ_{\mathrm{450}}=$ 0.59 mJy beam$^{-1}$ and $σ_{\mathrm{850}}=$ 0.09 mJy beam$^{-1}$ in the deepest area of each map. The co… ▽ More We present confusion-limited SCUBA-2 450-$μ$m observations in the COSMOS-CANDELS region as part of the JCMT Large Program, SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES). Our maps at 450 and 850 $μ$m cover an area of 450 arcmin$^2$. We achieved instrumental noise levels of $σ_{\mathrm{450}}=$ 0.59 mJy beam$^{-1}$ and $σ_{\mathrm{850}}=$ 0.09 mJy beam$^{-1}$ in the deepest area of each map. The corresponding confusion noise levels are estimated to be 0.65 and 0.36 mJy beam$^{-1}$. Above the 4 (3.5) $σ$ threshold, we detected 360 (479) sources at 450 $μ$m and 237 (314) sources at 850 $μ$m. We derive the deepest blank-field number counts at 450 $μ$m, covering the flux-density range of 2 to 43 mJy. These are in agreement with other SCUBA-2 blank-field and lensing-cluster observations, but are lower than various model counts. We compare the counts with those in other fields and find that the field-to-field variance observed at 450 $μ$m at the $R=6^\prime$ scale is consistent with Poisson noise, so there is no evidence of strong 2-D clustering at this scale. Additionally, we derive the integrated surface brightness at 450 $μ$m down to 2.1 mJy to be $57.3^{+1.0}_{-6.2}$~Jy deg$^{-2}$, contributing to (41$\pm$4)\% of the 450-$μ$m extragalactic background light (EBL) measured by COBE and Planck. Our results suggest that the 450-$μ$m EBL may be fully resolved at $0.08^{+0.09}_{-0.08}$~mJy, which extremely deep lensing-cluster observations and next-generation submillimeter instruments with large aperture sizes may be able to achieve. △ Less

Submitted 31 May, 2024; originally announced May 2024.

Comments: 29 pages, 14 figures, accepted for publication in 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)

Galaxy Morphologies Revealed with Subaru HSC and Super-Resolution Techniques II: Environmental Dependence of Galaxy Mergers at z~2-5

Authors: Takatoshi Shibuya, Yohito Ito, Kenta Asai, Takanobu Kirihara, Seiji Fujimoto, Yoshiki Toba, Noriaki Miura, Takuya Umayahara, Kenji Iwadate, Sadman S. Ali, Tadayuki Kodama

Abstract: We super-resolve the seeing-limited Subaru Hyper Suprime-Cam (HSC) images for 32,187 galaxies at z~2-5 in three techniques, namely, the classical Richardson-Lucy (RL) point spread function (PSF) deconvolution, sparse modeling, and generative adversarial networks to investigate the environmental dependence of galaxy mergers. These three techniques generate overall similar high spatial resolution im… ▽ More We super-resolve the seeing-limited Subaru Hyper Suprime-Cam (HSC) images for 32,187 galaxies at z~2-5 in three techniques, namely, the classical Richardson-Lucy (RL) point spread function (PSF) deconvolution, sparse modeling, and generative adversarial networks to investigate the environmental dependence of galaxy mergers. These three techniques generate overall similar high spatial resolution images but with some slight differences in galaxy structures, for example, more residual noises are seen in the classical RL PSF deconvolution. To alleviate disadvantages of each technique, we create combined images by averaging over the three types of super-resolution images, which result in galaxy sub-structures resembling those seen in the Hubble Space Telescope images. Using the combined super-resolution images, we measure the relative galaxy major merger fraction corrected for the chance projection effect, f_merg, for galaxies in the ~300 deg^2-area data of the HSC Strategic Survey Program and the CFHT Large Area U-band Survey. Our f_merg measurements at z~3 validate previous findings showing that f_merg is higher in regions with a higher galaxy overdensity delta at z~2-3. Thanks to the large galaxy sample, we identify a nearly linear increase in f_merg with increasing delta at z~4-5, providing the highest-z observational evidence that galaxy mergers are related to delta. In addition to our f_merg measurements, we find that the galaxy merger fractions in the literature also broadly align with the linear f_merg-delta relation across a wide redshift range of z~2-5. This alignment suggests that the linear f_merg-delta relation can serve as a valuable tool for quantitatively estimating the contributions of galaxy mergers to various environmental dependences. This super-resolution analysis can be readily applied to datasets from wide field-of-view space telescopes such as Euclid and Roman. △ Less

Submitted 11 March, 2024; originally announced March 2024.

Comments: 29 pages, 13 figures, 6 tables. Submitted to PASJ. Comments welcome

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

Steep-spectrum AGN in eROSITA Final Equatorial-Depth Survey (eFEDS): Their host galaxies and multi-wavelength properties

Authors: K. Iwasawa, T. Liu, Th. Boller, J. Buchner, J. Li, T. Kawaguchi, T. Nagao, Y. Terashima, Y. Toba, J. D. Silverman, R. Arcodia, Th. Dauser, M. Krumpe, K. Nandra, J. Wilms

Abstract: We selected sources with a steep soft-X-ray-band spectrum with a photon index larger than 2.5 -- measured by eROSITA on board the Spectrum-Roentgen-Gamma (SRG) -- from the eFEDS AGN catalogue as candidates of highly accreting supermassive black holes, and investigated their multi-wavelength properties. Among 601 bright AGN with 0.2-5 keV counts of greater than 100, 83 sources (~14%) are classified… ▽ More We selected sources with a steep soft-X-ray-band spectrum with a photon index larger than 2.5 -- measured by eROSITA on board the Spectrum-Roentgen-Gamma (SRG) -- from the eFEDS AGN catalogue as candidates of highly accreting supermassive black holes, and investigated their multi-wavelength properties. Among 601 bright AGN with 0.2-5 keV counts of greater than 100, 83 sources (~14%) are classified as steep-spectrum sources. These sources have typical 0.5-2 keV luminosities of L(SX) ~ 1e44 erg/s and the majority of them are found at redshifts below z=1. In comparison with sources with flatter spectra, these sources have, on average, a UV (or optical) to 2 keV luminosity ratio that is larger by ~0.3 dex and bluer optical-to-UV continuum emission. They also appear to be radio quiet based on the detection rate in the FIRST and VLASS surveys. Their host galaxies -- at least in the redshift range of z=0.2-0.8, where the AGN-galaxy decomposition results from the Subaru Hyper Suprime-Cam imaging are available -- tend to be late-type and have smaller stellar masses than those of sources with flatter spectra. These properties are similar to those found in nearby narrow-line Seyfert 1 galaxies, in agreement with the picture that they are AGN with elevated accretion rates and are in the early growth phase of black hole and galaxy co-evolution. However, the steep-spectrum sources are not exclusively narrow-line Seyfert 1 galaxies; indeed many are broad-line Seyfert 1 galaxies, as found by a catalogue search. This suggests that these steep-spectrum sources may be black holes generally with high accretion rates but of a wide mass range, including a few objects emitting at L(SX)>1e45 erg/s, of which black hole masses can be close to 10^9 M_sun. △ Less

Submitted 21 February, 2024; originally announced February 2024.

Comments: 8 pages, accepted for publication in A&A

AGN properties of ~1 million member galaxies of galaxy groups and clusters at z < 1.4 based on the Subaru Hyper Suprime-Cam survey

Authors: Yoshiki Toba, Aoi Hashiguchi, Naomi Ota, Masamune Oguri, Nobuhiro Okabe, Yoshihiro Ueda, Masatoshi Imanishi, Atsushi J. Nishizawa, Tomotsugu Goto, Bau-Ching Hsieh, Marie Kondo, Shuhei Koyama, Kianhong Lee, Ikuyuki Mitsuishi, Tohru Nagao, Taira Oogi, Koki Sakuta, Malte Schramm, Anri Yanagawa, Anje Yoshimoto

Abstract: Herein, we present the statistical properties of active galactic nuclei (AGNs) for approximately 1 million member galaxies of galaxy groups and clusters, with 0.1 $<$ cluster redshift ($z_{\rm cl}$) $<$ 1.4, selected using Subaru Hyper Suprime-Cam, the so-called CAMIRA clusters. In this research, we focused on the AGN power fraction ($f_{\rm AGN}$), which is defined as the proportion of the contri… ▽ More Herein, we present the statistical properties of active galactic nuclei (AGNs) for approximately 1 million member galaxies of galaxy groups and clusters, with 0.1 $<$ cluster redshift ($z_{\rm cl}$) $<$ 1.4, selected using Subaru Hyper Suprime-Cam, the so-called CAMIRA clusters. In this research, we focused on the AGN power fraction ($f_{\rm AGN}$), which is defined as the proportion of the contribution of AGNs to the total infrared (IR) luminosity, $L_{\rm IR}$ (AGN)/$L_{\rm IR}$, and examined how $f_{\rm AGN}$ depends on (i) $z_{\rm cl}$ and (ii) the distance from the cluster center. We compiled multiwavelength data using the ultraviolet--mid-IR range. Moreover, we performed spectral energy distribution fits to determine $f_{\rm AGN}$ using the CIGALE code with the SKIRTOR AGN model. We found that (i) the value of $f_{\rm AGN}$ in the CAMIRA clusters is positively correlated with $z_{\rm cl}$, with the correlation slope being steeper than that for field galaxies, and (ii) $f_{\rm AGN}$ exhibits a high value at the cluster outskirts. These results indicate that the emergence of AGN population depends on the redshift and environment and that galaxy groups and clusters at high redshifts are important in AGN evolution. Additionally, we demonstrated that cluster--cluster mergers may enhance AGN activity at the outskirts of particularly massive galaxy clusters. Our findings are consistent with a related study on the CAMIRA clusters that was based on the AGN number fraction. △ Less

Submitted 16 February, 2024; originally announced February 2024.

Comments: 25 pages, 24 figures, and 3 tables, accepted for publication in ApJ. A value-added CAMIRA member galaxy catalog and the best-fit SED for each member galaxy will be available as FITS or machine-readable tables

Euclid preparation. Measuring detailed galaxy morphologies for Euclid with Machine Learning

Authors: Euclid Collaboration, B. Aussel, S. Kruk, M. Walmsley, M. Huertas-Company, M. Castellano, C. J. Conselice, M. Delli Veneri, H. Domínguez Sánchez, P. -A. Duc, U. Kuchner, A. La Marca, B. Margalef-Bentabol, F. R. Marleau, G. Stevens, Y. Toba, C. Tortora, L. Wang, N. Aghanim, B. Altieri, A. Amara, S. Andreon, N. Auricchio, M. Baldi, S. Bardelli , et al. (233 additional authors not shown)

Abstract: The Euclid mission is expected to image millions of galaxies with high resolution, providing an extensive dataset to study galaxy evolution. We investigate the application of deep learning to predict the detailed morphologies of galaxies in Euclid using Zoobot a convolutional neural network pretrained with 450000 galaxies from the Galaxy Zoo project. We adapted Zoobot for emulated Euclid images, g… ▽ More The Euclid mission is expected to image millions of galaxies with high resolution, providing an extensive dataset to study galaxy evolution. We investigate the application of deep learning to predict the detailed morphologies of galaxies in Euclid using Zoobot a convolutional neural network pretrained with 450000 galaxies from the Galaxy Zoo project. We adapted Zoobot for emulated Euclid images, generated based on Hubble Space Telescope COSMOS images, and with labels provided by volunteers in the Galaxy Zoo: Hubble project. We demonstrate that the trained Zoobot model successfully measures detailed morphology for emulated Euclid images. It effectively predicts whether a galaxy has features and identifies and characterises various features such as spiral arms, clumps, bars, disks, and central bulges. When compared to volunteer classifications Zoobot achieves mean vote fraction deviations of less than 12% and an accuracy above 91% for the confident volunteer classifications across most morphology types. However, the performance varies depending on the specific morphological class. For the global classes such as disk or smooth galaxies, the mean deviations are less than 10%, with only 1000 training galaxies necessary to reach this performance. For more detailed structures and complex tasks like detecting and counting spiral arms or clumps, the deviations are slightly higher, around 12% with 60000 galaxies used for training. In order to enhance the performance on complex morphologies, we anticipate that a larger pool of labelled galaxies is needed, which could be obtained using crowdsourcing. Finally, our findings imply that the model can be effectively adapted to new morphological labels. We demonstrate this adaptability by applying Zoobot to peculiar galaxies. In summary, our trained Zoobot CNN can readily predict morphological catalogues for Euclid images. △ Less

Submitted 15 February, 2024; originally announced February 2024.

Comments: 27 pages, 26 figures, 5 tables, submitted to A&A

[O IV] and [Ne V]-weak AGNs Hidden by Compton-thick Material in Late Mergers

Authors: Satoshi Yamada, Yoshihiro Ueda, Taiki Kawamuro, Claudio Ricci, Yoshiki Toba, Masatoshi Imanishi, Takamitsu Miyaji, Atsushi Tanimoto, Kohei Ichikawa, Martin Herrera-Endoqui, Shoji Ogawa, Ryosuke Uematsu, Keiichi Wada

Abstract: We study "buried" active galactic nuclei (AGNs) almost fully covered by circumnuclear material in ultra-/luminous infrared galaxies (U/LIRGs), which show weak ionized lines from narrow line regions. Employing an indicator of [O IV] 25.89-um or [Ne V] 14.32-um line to 12-um AGN luminosity ratio, we find 17 buried AGN candidates that are [O IV]-weak ($L_{\rm [O\,IV]}$/$L_{\rm 12,AGN} \leq -$3.0) or… ▽ More We study "buried" active galactic nuclei (AGNs) almost fully covered by circumnuclear material in ultra-/luminous infrared galaxies (U/LIRGs), which show weak ionized lines from narrow line regions. Employing an indicator of [O IV] 25.89-um or [Ne V] 14.32-um line to 12-um AGN luminosity ratio, we find 17 buried AGN candidates that are [O IV]-weak ($L_{\rm [O\,IV]}$/$L_{\rm 12,AGN} \leq -$3.0) or [Ne V]-weak ($L_{\rm [Ne\,V]}$/$L_{\rm 12,AGN} \leq -$3.4) among 30 AGNs in local U/LIRGs. For the [O IV]-weak AGNs, we estimate their covering fractions of Compton-thick (CT; $N_{\rm H} \geq 10^{24}$ cm$^{-2}$) material with an X-ray clumpy torus model to be $f^{\rm (spec)}_{\rm CT} = 0.55\pm0.19$ on average. This value is consistent with the fraction of CT AGNs ($f^{\rm (stat)}_{\rm CT} = 53\pm12$%) among the [O IV]-weak AGNs in U/LIRGs and much larger than that in Swift/BAT AGNs ($23\pm6$%). The fraction of [O IV]-weak AGNs increases from $27^{+13}_{-10}$% (early) to $66^{+10}_{-12}$% (late mergers). Similar results are obtained with the [Ne V] line. The [O IV] or [Ne V]-weak AGNs in late mergers show larger $N_{\rm H}$ and Eddington ratios ($λ_{\rm Edd}$) than those of the Swift/BAT AGNs, and the largest $N_{\rm H}$ is $\gtrsim$10$^{25}$ cm$^{-2}$ at ${\log}λ_{\rm Edd} \sim -$1, close to the effective Eddington limit for CT material. These suggest that (1) the circumnuclear material in buried AGNs is regulated by the radiation force from high-$λ_{\rm Edd}$ AGNs on the CT obscurers, and (2) their dense material with large $f^{\rm (spec)}_{\rm CT}$ ($\sim 0.5 \pm 0.1$) in U/LIRGs is a likely cause of a unique structure of buried AGNs, whose amount of material may be maintained through merger-induced supply from their host galaxies. △ Less

Submitted 8 February, 2024; originally announced February 2024.

Comments: 13 pages, 10 figures, and 1 table, accepted for publication in ApJ

ALMA Lensing Cluster Survey: Full SED Analysis of z~0.5-6 Lensed Galaxies Detected with Millimeter Observations

Authors: Ryosuke Uematsu, Yoshihiro Ueda, Kotaro Kohno, Yoshiki Toba, Satoshi Yamada, Ian Smail, Hideki Umehata, Seiji Fujimoto, Bunyo Hatsukade, Yiping Ao, Franz Erik Bauer, Gabriel Brammer, Miroslava Dessauges-Zavadsky, Daniel Espada, Jean-Baptiste Jolly, Anton M. Koekemoer, Vasily Kokorev, Georgios E. Magdis, Masamune Oguri, Fengwu Sun

Abstract: Sub/millimeter galaxies are a key population for the study of galaxy evolution because the majority of star formation at high redshifts occurred in galaxies deeply embedded in dust. To search for this population, we have performed an extensive survey with ALMA, called the ALMA Lensing Cluster Survey (ALCS). This survey covers 133 arcmin^2 area and securely detects 180 sources at z~0.5-6 with a flu… ▽ More Sub/millimeter galaxies are a key population for the study of galaxy evolution because the majority of star formation at high redshifts occurred in galaxies deeply embedded in dust. To search for this population, we have performed an extensive survey with ALMA, called the ALMA Lensing Cluster Survey (ALCS). This survey covers 133 arcmin^2 area and securely detects 180 sources at z~0.5-6 with a flux limit of ~0.2 mJy at 1.2 mm (Fujimoto et al. 2023). Here we report the results of multi-wavelength spectral energy distribution (SED) analysis of the whole ALCS sample, utilizing the observed-frame UV to millimeter photometry. We find that the majority of the ALCS sources lie on the star-forming main sequence, with a smaller fraction showing intense starburst activities. The ALCS sample contains high infrared-excess sources IRX=log(Ldust/LUV)>1), including two extremely dust-obscured galaxies (IRX>5). We also confirm that the ALCS sample probes a broader range in lower dust mass than conventional SMG samples in the same redshift range. We identify six heavily obscured AGN candidates that are not detected in the archival Chandra data in addition to the three X-ray AGNs reported by Uematsu et al. (2023). The inferred AGN luminosity density shows a possible excess at z=2-3 compared with that determined from X-ray surveys below 10 keV. △ Less

Submitted 13 February, 2024; v1 submitted 8 February, 2024; originally announced February 2024.

Comments: 27 pages, 20 figures

Ionised AGN outflows in the Goldfish galaxy -- The illuminating and interacting red quasar eFEDSJ091157.4+014327 at z $\sim$ 0.6

Authors: Blessing Musiimenta, Giovanna Speranza, Tanya Urrutia, Marcella Brusa, Cristina Ramos Almeida, Michele Perna, Ivàn Ezequiel López, David M. Alexander, Brivael Laloux, Francesco Shankar, Andrea Lapi, Mara Salvato, Yoshiki Toba, Carolina Andonie, Ivàn Munoz Rodríguez

Abstract: Evolutionary models suggest that the initial growth phases of active galactic nuclei (AGN) are dust-enshrouded, and characterised by jet/wind outflows that should gradually clear the interstellar medium (ISM) in the host by heating and/or expelling the surrounding gas. eFEDSJ091157.4$+$014327 (z$\sim$0.6) was selected from X-ray samples for its characteristics that are similar to sources with outf… ▽ More Evolutionary models suggest that the initial growth phases of active galactic nuclei (AGN) are dust-enshrouded, and characterised by jet/wind outflows that should gradually clear the interstellar medium (ISM) in the host by heating and/or expelling the surrounding gas. eFEDSJ091157.4$+$014327 (z$\sim$0.6) was selected from X-ray samples for its characteristics that are similar to sources with outflows which include red, obscured and X-ray luminous. We aim to explore the environment around this red quasar and characterise kinematics within the system. We used spatially resolved spectroscopic data from Multi Unit Spectroscopic Explorer (MUSE) with an average seeing of 0.6" to construct flux, velocity and velocity dispersions maps. We found that the quasar is embedded in an interacting and merging system with three other galaxies $\sim$ 50 kpc from its nucleus. Spatially resolved kinematics reveal that the quasar has extended ionised outflows of up to 9.2 kpc with positive and negative velocities up to 1000 km s$^{-1}$ and -1200 km s$^{-1}$, respectively. The velocity dispersion (W$_{80}$) ranges from 600-1800 km s$^{-1}$. We associate the presence of turbulent and high-velocity components with the outflow. The total mass outflow rate is estimated to be $\sim$ 10 M$_{\odot}$ yr$^{-1}$ and kinetic power of 2x10$^{42}$ erg s$^{-1}$. The kinetic coupling efficiencies range from 0.01%-0.03% and the momentum boosts are $\sim$ 0.2. These low values indicate that the ionised outflow is not energetically relevant. These values don't align with the theoretical predictions of both radiation-pressure-driven outflows and energy-conserving mechanisms. However, note that our results are based only on the ionised phase while theoretical predictions are multiphase. Moreover, the mass loading factor of $\sim$ 5 is an indication that these outflows are more likely AGN-driven than star formation-driven. △ Less

Submitted 6 May, 2024; v1 submitted 30 January, 2024; originally announced January 2024.

Comments: Version after addressing referee comments, accepted in Astronomy and Astrophysics(A&A), 17 pages, 12 figures, 1 table

Journal ref: A&A 687, A111 (2024)

AGN X-ray luminosity function and absorption function in the Early Universe ($3\leq z \leq 6$)

Authors: E. Pouliasis, A. Ruiz, I. Georgantopoulos, F. Vito, R. Gilli, C. Vignali, Y. Ueda, E. Koulouridis, M. Akiyama, S. Marchesi, B. Laloux, T. Nagao, S. Paltani, M. Pierre, Y. Toba, M. Habouzit, B. Vijarnwannaluk, C. Garrel

Abstract: The XLF of AGN offers a robust tool to study the evolution and the growth of SMBHs over cosmic time. Owing to the limited area probed by X-ray surveys, optical surveys are routinely used to probe the accretion in the high redshift Universe $z\geq 3$. However, optical surveys may be incomplete because they are strongly affected by dust redenning. In this work, we derive the XLF and its evolution at… ▽ More The XLF of AGN offers a robust tool to study the evolution and the growth of SMBHs over cosmic time. Owing to the limited area probed by X-ray surveys, optical surveys are routinely used to probe the accretion in the high redshift Universe $z\geq 3$. However, optical surveys may be incomplete because they are strongly affected by dust redenning. In this work, we derive the XLF and its evolution at high redshifts using a large sample of AGNs selected in different fields with various areas and depths covering a wide range of luminosities. Additionally, we put the tightest yet constraints on the absorption function in this redshift regime. In particular, we use more than 600 soft X-ray selected high-z sources in the Chandra Deep fields, the Chandra COSMOS Legacy survey and the XMM-XXL northern field. We derive the X-ray spectral properties for all sources via spectral fitting, using a consistent technique and model. For modeling the parametric form of the XLF and the absorption function, we use a Bayesian methodology allowing us to correctly propagate the uncertainties for the observed X-ray properties of our sources and also the absorption effects. The evolution of XLF is in agreement with a pure density evolution model similar to what is witnessed at optical wavelengths, although a luminosity dependent density evolution model cannot be securely ruled out. A large fraction ($60\%)$ of our sources are absorbed by column densities of $\rm N_H \geq 10^{23} cm^{-2} $, while $17$\% of the sources are CTK. Our results favor a scenario where both the ISM of the host and the AGN torus contribute to the obscuration. The derived BHAD is in agreement with the simulations, if one takes into account that the X-ray AGN are hosted by massive galaxies, while it differs from the one derived using JWST data. The latter could be due to the differences in the AGN and host-galaxy properties. △ Less

Submitted 24 January, 2024; originally announced January 2024.

Comments: Accepted for publication in A&A. 25 pages, 15 figures (+4 in Appendix), 2 tables (+1 in Appendix)

Journal ref: A&A 685, A97 (2024)

A population of Optically Quiescent Quasars from WISE and SDSS

Authors: Claire Greenwell, Poshak Gandhi, Daniel Stern, George Lansbury, Vincenzo Mainieri, Peter Boorman, Yoshiki Toba

Abstract: The growth of active galactic nuclei (AGN) occurs under some form of obscuration in a large fraction of the population. The difficulty in constraining this population leads to high uncertainties in cosmic X-ray background and galaxy evolution models. Using an SDSS-WISE cross-match, we target infrared luminous AGN ($W1-W2$ > 0.8, and monochromatic rest-frame luminosity above $λL_λ$(12$μ$m)… ▽ More The growth of active galactic nuclei (AGN) occurs under some form of obscuration in a large fraction of the population. The difficulty in constraining this population leads to high uncertainties in cosmic X-ray background and galaxy evolution models. Using an SDSS-WISE cross-match, we target infrared luminous AGN ($W1-W2$ > 0.8, and monochromatic rest-frame luminosity above $λL_λ$(12$μ$m) $\approx$ 3 $\times$ 10$^{44}$ erg s$^{-1}$), but with passive galaxy-like optical spectra (Optically Quiescent Quasars; OQQs). We find 47 objects that show no significant [O III]$λ$5007 emission, a typically strong AGN optical emission line. As a comparison sample, we examine SDSS-selected Type 2 quasars (QSO2s), which show a significant [O III]$λ$5007 line by definition. We find a 1:16 ratio of OQQs compared to QSO2s, suggesting that the OQQ duty cycle is likely much shorter than that of QSO2s (though selection biases are not fully quantified). We consider observed properties in comparison with other galaxy types, and examine them for consistency with theories on their intrinsic nature: chiefly (a) a high covering factor for surrounding obscuring matter, preventing the detection of high-ionisation emission lines - `cocooned AGN'; or (b) ionised gas being absent on the kpc scales of the Narrow Line Region (NLR), perhaps due to a `switching on' or `young' AGN. OQQs do not obviously fit the standard paradigm for merger-driven AGN and host galaxy evolution, implying we may be missing part of the flow of AGN evolution. △ Less

Submitted 3 January, 2024; originally announced January 2024.

Comments: Accepted for publication in MNRAS (20 pages, 18 figures)

Multiwavelength SED Analysis of X-Ray Selected AGNs at $z=0.2-0.8$ in Stripe 82 Region

Authors: Kenta Setoguchi, Yoshihiro Ueda, Yoshiki Toba, Junyao Li, John Silverman, Ryosuke Uematsu

Abstract: We perform a systematic, multiwavelength spectral energy distribution (SED) analysis of X-ray detected Active Galactic Nuclei (AGNs) at $z=0.2-0.8$ with SDSS counterparts in the Stripe 82 region, consisting of 60 type-1 and 137 type-2 AGNs covering a 2--10 keV luminosity range of $41.6 < {\rm log}\ L_{\rm x} < 44.7$. The latest CIGALE code, where dusty polar components are included, is employed. T… ▽ More We perform a systematic, multiwavelength spectral energy distribution (SED) analysis of X-ray detected Active Galactic Nuclei (AGNs) at $z=0.2-0.8$ with SDSS counterparts in the Stripe 82 region, consisting of 60 type-1 and 137 type-2 AGNs covering a 2--10 keV luminosity range of $41.6 < {\rm log}\ L_{\rm x} < 44.7$. The latest CIGALE code, where dusty polar components are included, is employed. To obtain reliable host and AGN parameters in type-1 AGNs, we utilize the image decomposed optical SED of host galaxies by Li et al. (2021) based on the Subaru Hyper Suprime-Cam (HSC) images. The mean ratio of black hole masses ($M_{\rm BH}$) and stellar masses ($M_{\rm stellar}$) of our X-ray detected type-1 AGN sample, $\log (M_{\rm BH}/M_{\rm stellar}) = -2.7\pm0.5$, is close to the local relation between black hole and stellar masses, as reported by Li et al. (2021) for SDSS quasars. This ratio is slightly lower than that found for more luminous ($\log L_{\rm bol} > 45$) type-1 AGNs at $z\sim1.5$. This can be explained by the AGN-luminosity dependence of $\log (M_{\rm BH}/M_{\rm stellar})$, which little evolves with redshift. We confirm the trend that the UV-to-X-ray slope ($α_{\rm OX}$) or X-ray-to-bolometric correction factor ($κ_{2-10}$) increases with AGN luminosity or Eddington ratio. We find that type-1 and type-2 AGNs with the same luminosity ranges share similar host stellar-mass distributions, whereas type-2s tend to show smaller AGN luminosities than type-1s. This supports the luminosity (or Eddington ratio) dependent unified scheme. △ Less

Submitted 6 December, 2023; originally announced December 2023.

Comments: 22 pages, 9 figures, 7 tables, accepted for ApJ publication

Apparent effect of dust extinction on the observed outflow velocity of ionized gas in galaxy mergers

Authors: Naomichi Yutani, Yoshiki Toba, Keiichi Wada

Abstract: In this study, we examine photoionization outflows during the late stages of galaxy mergers, with a specific focus on the relation between observed velocity of outflowing gas and the apparent effects of dust extinction. We used the N-body/smoothed particle hydrodynamics (SPH) code ASURA for galaxy merger simulations. These simulations concentrated on identical galaxy mergers featuring supermassive… ▽ More In this study, we examine photoionization outflows during the late stages of galaxy mergers, with a specific focus on the relation between observed velocity of outflowing gas and the apparent effects of dust extinction. We used the N-body/smoothed particle hydrodynamics (SPH) code ASURA for galaxy merger simulations. These simulations concentrated on identical galaxy mergers featuring supermassive black holes (SMBHs) of 10$^8$ M$_\odot$ and gas fractions of 30\% and 10 \%. From the simulation data, we derived velocity and velocity dispersion diagrams for the AGN-driven ionized outflowing gas. Our findings show that high-velocity outflows with velocity dispersions of 500 km s$^{-1}$ or greater can be observed in the late stages of galactic mergers. Particularly, in buried AGNs, both the luminosity-weighted outflow velocity and velocity dispersion increase owing to the apparent effects of dust extinction. Owing to these effects, the velocity--velocity dispersion diagrams display a noticeable blue-shifted tilt in models with higher gas fractions. Crucially, this tilt is not influenced by the AGN luminosity but emerges from the observational impacts of dust extinction. Our results imply that the observed high-velocity \OIII outflow exceeding 1000 km s$^{-1}$ in buried AGNs may be linked to the dust extinction that occurs during the late stages of gas-rich galaxy mergers. △ Less

Submitted 19 November, 2023; originally announced November 2023.

Comments: 14 pages, accepted for publication in ApJ

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

BASS XLII: The relation between the covering factor of dusty gas and the Eddington ratio in nearby active galactic nuclei

Authors: C. Ricci, K. Ichikawa, M. Stalevski, T. Kawamuro, S. Yamada, Y. Ueda, R. Mushotzky, G. C. Privon, M. J. Koss, B. Trakhtenbrot, A. C. Fabian, L. C. Ho, D. Asmus, F. E. Bauer, C. S. Chang, K. K. Gupta, K. Oh, M. Powell, R. W. Pfeifle, A. Rojas, F. Ricci, M. J. Temple, Y. Toba, A. Tortosa, E. Treister , et al. (3 additional authors not shown)

Abstract: Accreting supermassive black holes (SMBHs) located at the center of galaxies are typically surrounded by large quantities of gas and dust. The structure and evolution of this circumnuclear material can be studied at different wavelengths, from the submillimeter to the X-rays. Recent X-ray studies have shown that the covering factor of the obscuring material tends to decrease with increasing Edding… ▽ More Accreting supermassive black holes (SMBHs) located at the center of galaxies are typically surrounded by large quantities of gas and dust. The structure and evolution of this circumnuclear material can be studied at different wavelengths, from the submillimeter to the X-rays. Recent X-ray studies have shown that the covering factor of the obscuring material tends to decrease with increasing Eddington ratio, likely due to radiative feedback on dusty gas. Here we study a sample of 549 nearby (z<0.1) hard X-ray (14-195 keV) selected non-blazar active galactic nuclei (AGN), and use the ratio between the AGN infrared and bolometric luminosity as a proxy of the covering factor. We find that, in agreement with what has been found by X-ray studies of the same sample, the covering factor decreases with increasing Eddington ratio. We also confirm previous findings which showed that obscured AGN typically have larger covering factors than unobscured sources. Finally, we find that the median covering factors of AGN located in different regions of the column density-Eddington ratio diagram are in good agreement with what would be expected from a radiation-regulated growth of SMBHs. △ Less

Submitted 11 November, 2023; v1 submitted 2 November, 2023; originally announced November 2023.

Comments: ApJ in press

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

New technique to select recent fast-quenching galaxies at $z\sim2$ using the optical colors

Authors: Mariko Kubo, Tohru Nagao, Hisakazu Uchiyama, Takuji Yamashita, Yoshiki Toba, Masaru Kajisawa, Yuta Yamamoto

Abstract: Many massive quiescent galaxies have been discovered at $z>2$ thanks to multi-wavelength deep and wide surveys, however, substantial deep near-infrared spectroscopic observations are needed to constrain their star-formation histories statistically. Here, we present a new technique to select quiescent galaxies with a short quenching timescale ($\leq0.1$ Gyr) at $z\sim2$ photometrically. We focus on… ▽ More Many massive quiescent galaxies have been discovered at $z>2$ thanks to multi-wavelength deep and wide surveys, however, substantial deep near-infrared spectroscopic observations are needed to constrain their star-formation histories statistically. Here, we present a new technique to select quiescent galaxies with a short quenching timescale ($\leq0.1$ Gyr) at $z\sim2$ photometrically. We focus on a spectral break at $\sim1600$ Å~that appears for such fast-quenching galaxies $\sim1$ Gyr after quenching when early A-type stars go out, but late A-type stars still live. This spectral break at $z\sim2$ is similar to a Lyman break at $z\sim4$. We construct a set of color criteria for $z\sim2$ fast-quenching galaxies on $g-r$ vs. $r-i$ and $i-J$ vs. $J-H$ or $\rm i-[3.6]$ vs. $\rm [3.6]-[4.5]$ color diagrams, which are available with the existing and/or future wide imaging surveys, by simulating various model galaxy spectra and test their robustnesses using the COSMOS2020 catalog. Galaxies with photometric and/or spectroscopic redshifts $z\sim2$ and low specific star formation rates are successfully selected using these colors. The number density of these fast-quenching galaxy candidates at $z\sim2$ suggests that massive galaxies not so far above the star-formation main sequence at $z=3-4$ should be their progenitors. △ Less

Submitted 14 October, 2023; originally announced October 2023.

Comments: 12 pages, 11 figures, accepted for publication in MNRAS

A new discovery space opened by eROSITA: Ionised AGN outflows from X-ray selected samples

Authors: Blessing Musiimenta, Marcella Brusa, Teng Liu, Mara Salvato, Johannes Buchner, Zsofi Igo, Sophia G. H. Waddell, Yoshiki Toba, Riccardo Arcodia, Johan Comparat, David M. Alexander, Francesco Shankar, Andrea Lapi, Cristina Ramos Almeida, Antonis Georgakakis, Andrea Merloni, Tanya Urrutia, Junyao Li, Yuichi Terashima, Yue Shen, Qiaoya Wu, Tom Dwelly, Kirpal Nandra, Julien Wolf

Abstract: In the context of an evolutionary model, the outflow phase of an Active Galactic Nuclei (AGN) occurs at the peak of its activity, once the central SMBH is massive enough to generate sufficient power to counterbalance the potential well of the host galaxy. This phase plays a vital role in galaxy evolution. We aim to apply various selection methods to isolate powerful AGNs in the feedback phase, tra… ▽ More In the context of an evolutionary model, the outflow phase of an Active Galactic Nuclei (AGN) occurs at the peak of its activity, once the central SMBH is massive enough to generate sufficient power to counterbalance the potential well of the host galaxy. This phase plays a vital role in galaxy evolution. We aim to apply various selection methods to isolate powerful AGNs in the feedback phase, trace and characterise their outflows, and explore the link between AGN luminosity and outflow properties. We applied a combination of methods to the eROSITA Final Equatorial Depth survey (eFEDS) catalogue and isolated ~1400 candidates at z>0.5 out of ~11750 AGNs (~12\%). We tested the robustness of our selection on the small subsample of 50 sources with available good quality SDSS spectra at 0.5<z<1, for which we fitted the [OIII] emission line complex and searched for the presence of ionised gas outflows. We identified 23 quasars (~45\%) with evidence of ionised outflows based on the presence of significant broad and shifted components in the [OIII] line. They are on average more luminous and more obscured than the parent sample, although this may be ascribed to selection effects affecting the good quality SDSS spectra sample. By adding 118 outflowing quasars at 0.5<z<3.5 from the literature, we find a weak correlation between the maximum outflow velocity and AGN bolometric luminosity. On the contrary, we find strong correlations between mass outflow rate and outflow kinetic power with the AGN bolometric luminosity. About 30\% of our sample have kinetic coupling efficiencies >1\%. We find that the majority of the outflows have momentum flux ratios lower than 20 which rules out an energy-conserving nature. Our present work points to the unequivocal existence of a rather short AGN outflow phase, paving the way towards a new avenue to dissect AGN outflows in large samples within eROSITA and beyond. △ Less

Submitted 28 September, 2023; originally announced September 2023.

Comments: Accepted for publication in Astronomy and Astrophysics(A&A), 24 pages, 18 figures, 5 tables

Journal ref: A&A 679, A84 (2023)

Galaxy mergers in Subaru HSC-SSP: a deep representation learning approach for identification and the role of environment on merger incidence

Authors: Kiyoaki Christopher Omori, Connor Bottrell, Mike Walmsley, Hassen M. Yesuf, Andy D. Goulding, Xuheng Ding, Gergö Popping, John D. Silverman, Tsutomu T. Takeuchi, Yoshiki Toba

Abstract: We take a deep learning-based approach for galaxy merger identification in Subaru HSC-SSP, specifically through the use of deep representation learning and fine-tuning, with the aim of creating a pure and complete merger sample within the HSC-SSP survey. We can use this merger sample to conduct studies on how mergers affect galaxy evolution. We use Zoobot, a deep learning representation learning m… ▽ More We take a deep learning-based approach for galaxy merger identification in Subaru HSC-SSP, specifically through the use of deep representation learning and fine-tuning, with the aim of creating a pure and complete merger sample within the HSC-SSP survey. We can use this merger sample to conduct studies on how mergers affect galaxy evolution. We use Zoobot, a deep learning representation learning model pre-trained on citizen science votes on Galaxy Zoo DeCALS images. We fine-tune Zoobot for the purpose of merger classification of images of SDSS and GAMA galaxies in HSC-SSP PDR 3. Fine-tuning is done using 1200 synthetic HSC-SSP images of galaxies from the TNG simulation. We then find merger probabilities on observed HSC images using the fine-tuned model. Using our merger probabilities, we examine the relationship between merger activity and environment. We find that our fine-tuned model returns an accuracy on the synthetic validation data of 76%. This number is comparable to those of previous studies where convolutional neural networks were trained with simulation images, but with our work requiring a far smaller number of training samples. For our synthetic data, our model is able to achieve completeness and precision values of 80%. In addition, our model is able to correctly classify both mergers and non-mergers of diverse morphologies and structures, including those at various stages and mass ratios, while distinguishing between projections and merger pairs. For the relation between galaxy mergers and environment, we find two distinct trends. Using stellar mass overdensity estimates for TNG simulations and observations using SDSS and GAMA, we find that galaxies with higher merger scores favor lower density environments on scales of 0.5 to 8 h^-1 Mpc. However, below these scales in the simulations, we find that galaxies with higher merger scores favor higher density environments. △ Less

Submitted 27 September, 2023; originally announced September 2023.

Comments: 36 pages, 15 figures, accepted to Astronomy and Astrophysics

AGN number fraction in galaxy groups and clusters at z < 1.4 from the Subaru Hyper Suprime-Cam survey

Authors: Aoi Hashiguchi, Yoshiki Toba, Naomi Ota, Masamune Oguri, Nobuhiro Okabe, Yoshihiro Ueda, Masatoshi Imanishi, Satoshi Yamada, Tomotsugu Goto, Shuhei Koyama, Kianhong Lee, Ikuyuki Mitsuishi, Tohru Nagao, Atsushi J. Nishizawa, Akatoki Noboriguchi, Taira Oogi, Koki Sakuta, Malte Schramm, Mio Shibata, Yuichi Terashima, Takuji Yamashita, Anri Yanagawa, Anje Yoshimoto

Abstract: One of the key questions on active galactic nuclei (AGN) in galaxy clusters is how AGN could affect the formation and evolution of member galaxies and galaxy clusters in the history of the Universe. To address this issue, we investigate the dependence of AGN number fraction ($f_{\rm AGN}$) on cluster redshift ($z_{\rm cl}$) and distance from the cluster center ($R/R_{\rm 200}$). We focus on more t… ▽ More One of the key questions on active galactic nuclei (AGN) in galaxy clusters is how AGN could affect the formation and evolution of member galaxies and galaxy clusters in the history of the Universe. To address this issue, we investigate the dependence of AGN number fraction ($f_{\rm AGN}$) on cluster redshift ($z_{\rm cl}$) and distance from the cluster center ($R/R_{\rm 200}$). We focus on more than 27,000 galaxy groups and clusters at $0.1 < z_{\rm cl} < 1.4$ with more than 1 million member galaxies selected from the Subaru Hyper Suprime-Cam. By combining various AGN selection methods based on infrared (IR), radio, and X-ray data, we identify 2,688 AGN. We find that (i) $f_{\rm AGN}$ increases with $z_{\rm cl}$ and (ii) $f_{\rm AGN}$ decreases with $R/R_{\rm 200}$. The main contributors to the rapid increase of $f_{\rm AGN}$ towards high-$z$ and cluster center are IR- and radio-selected AGN, respectively. Those results indicate that the emergence of the AGN population depends on the environment and redshift, and galaxy groups and clusters at high-$z$ play an important role in AGN evolution. We also find that cluster-cluster mergers may not drive AGN activity in at least the cluster center, while we have tentative evidence that cluster-cluster mergers would enhance AGN activity in the outskirts of (particularly massive) galaxy clusters. △ Less

Submitted 4 September, 2023; originally announced September 2023.

Comments: 16 pages, 21 figures, and 2 tables, accepted for publication in PASJ

Similarity between compact extremely red objects discovered with JWST in cosmic dawn and blue-excess dust-obscured galaxies known in cosmic noon

Authors: Akatoki Noboriguchi, Akio K. Inoue, Tohru Nagao, Yoshiki Toba, Toru Misawa

Abstract: Spatially compact objects with extremely red color in the rest-frame optical to near-infrared (0.4--1 ${\rm μm}$) and blue color in the rest-frame ultraviolet (UV; 0.2--0.4 ${\rm μm}$) have been discovered at $5 < z < 9$ using the James Webb Space Telescope (JWST). These extremely red objects (JWST-EROs) exhibit spectral energy distributions (SEDs) that are difficult to explain using a single comp… ▽ More Spatially compact objects with extremely red color in the rest-frame optical to near-infrared (0.4--1 ${\rm μm}$) and blue color in the rest-frame ultraviolet (UV; 0.2--0.4 ${\rm μm}$) have been discovered at $5 < z < 9$ using the James Webb Space Telescope (JWST). These extremely red objects (JWST-EROs) exhibit spectral energy distributions (SEDs) that are difficult to explain using a single component of either star-forming galaxies or quasars, leading to two-component models in which the blue UV and extremely red optical are explained using less-dusty and dusty spectra of galaxies or quasars, respectively. Here, we report the remarkable similarity in SEDs between JWST-EROs and blue-excess dust-obscured galaxies (BluDOGs) identified at $2 < z < 3$. BluDOGs are a population of active galactic nuclei (AGNs) with blackhole masses of $\sim10^{8-9}$ M$_\odot$, which are one order of magnitude larger than those in some JWST-EROs. The Eddington ratios of BluDOGs are one or higher, whereas those of JWST-EROs are in the range of 0.1--1. Therefore, JWST-EROs are less massive, less active, and more common counterparts in higher-$z$ of BluDOGs in cosmic noon. Conversely, JWST-EROs have a significantly higher fraction of those with blue-excess than DOGs. We present the average UV spectra of BluDOGs as a comparison to JWST-EROs and discuss a coherent evolutionary scenario for dusty AGN populations. △ Less

Submitted 5 December, 2023; v1 submitted 2 September, 2023; originally announced September 2023.

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

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

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

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

Submitted 5 July, 2023; originally announced July 2023.

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

Quasar Luminosity Function at z = 7

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

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

Submitted 18 May, 2023; originally announced May 2023.

Comments: The Astrophysical Journal Letters, in press

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

Morphological asymmetries of quasar host galaxies with Subaru Hyper Suprime-Cam

Authors: Shenli Tang, John D. Silverman, Hassen M. Yesuf, Xuheng Ding, Junyao Li, Connor Bottrell, Andy Goulding, Kiyoaki Christopher Omori, Yoshiki Toba, Toshihiro Kawaguchi

Abstract: How does the host galaxy morphology influence a central quasar or vice versa? We address this question by measuring the asymmetries of 2424 SDSS quasar hosts at $0.2<z<0.8$ using broad-band ($grizy$) images from the Hyper Suprime-Cam Subaru Strategic Program. Control galaxies (without quasars) are selected by matching the redshifts and stellar masses of the quasar hosts. A two-step pipeline is run… ▽ More How does the host galaxy morphology influence a central quasar or vice versa? We address this question by measuring the asymmetries of 2424 SDSS quasar hosts at $0.2<z<0.8$ using broad-band ($grizy$) images from the Hyper Suprime-Cam Subaru Strategic Program. Control galaxies (without quasars) are selected by matching the redshifts and stellar masses of the quasar hosts. A two-step pipeline is run to decompose the PSF and \sersic\ components, and then measure asymmetry indices ($A_{\rm CAS}$, $A_{\rm outer}$, and $A_{\rm shape}$) of each quasar host and control galaxy. We find a mild correlation between host asymmetry and AGN bolometric luminosity ($L_{\rm bol}$) for the full sample (spearman correlation of 0.37) while a stronger trend is evident at the highest luminosities ($L_{\rm bol}>45$). This then manifests itself into quasar hosts being more asymmetric, on average, when they harbor a more massive and highly accreting black hole. The merger fraction also positively correlates with $L_{\rm bol}$ and reaches up to 35\% for the most luminous. Compared to control galaxies, quasar hosts are marginally more asymmetric (excess of 0.017 in median at 9.4$σ$ level) and the merger fractions are similar ($\sim 16.5\%$). We quantify the dependence of asymmetry on optical band which demonstrates that mergers are more likely to be identified with the bluer bands and the correlation between $L_{\rm bol}$ and asymmetry is also stronger in such bands. We stress that the band dependence, indicative of a changing stellar population, is an important factor in considering the influence of mergers on AGN activity. △ Less

Submitted 12 April, 2023; originally announced April 2023.

Comments: 27 pages, 28 figures

Journal ref: MNRAS, volume 521, number 4, pages 5272--5297 (2023)

J1406+0102: Dust Obscured Galaxy Hiding Super Eddington Accretion System with Bright Radio Emission

Authors: Hikaru Fukuchi, Kohei Ichikawa, Masayuki Akiyama, Shigeo Kimura, Yoshiki Toba, Kohei Inayoshi, Akatoki Noboriguchi, Toshihiro Kawaguchi, Xiaoyang Chen, Itsna K. Fitriana

Abstract: Recent high-$z$ quasar observations strongly indicate that super-Eddington accretion is a crucial phase to describe the existence of supermassive black holes (SMBHs) with $M_\mathrm{BH} \gtrsim 10^9 M_\odot$ at $z \gtrsim 7$. Motivated by the theoretical suggestion that the super-Eddington phase efficiently produces outflows and jets bright in radio bands, we search and find a super-Eddington radi… ▽ More Recent high-$z$ quasar observations strongly indicate that super-Eddington accretion is a crucial phase to describe the existence of supermassive black holes (SMBHs) with $M_\mathrm{BH} \gtrsim 10^9 M_\odot$ at $z \gtrsim 7$. Motivated by the theoretical suggestion that the super-Eddington phase efficiently produces outflows and jets bright in radio bands, we search and find a super-Eddington radio-loud dust-obscured galaxy (DOG) J1406+0102 at $z=0.236$, through cross-matching of the infrared-bright DOGs of Noboriguchi et al. (2019) with the VLA/FIRST 1.4 GHz radio and the SDSS optical spectral catalog. DOG J1406+0102 shows broad components in the Balmer lines. Assuming those lines are from the broad line region, it gives BH mass estimation of $\log\ (M_\mathrm{BH}/M_\odot)=7.30 \pm 0.25$, and AGN luminosity of $\log (L_\mathrm{bol,[OIII]}/\mathrm{erg}~\mathrm{s}^{-1}) = 45.91\pm0.38$ estimated from the intrinsic [OIII] luminosity, resulting in super-Eddington accretion of $λ_\mathrm{Edd}\simeq 3$. We show that 1) DOG J1406+0102 is operating strong AGN feedback: the [OIII] outflow velocity exceeds the escape velocity of the host galaxy halo and the kinetic efficiency is obtained as $\approx$ 8% that can be sufficient to quench the host galaxy, 2) the expected future growth pathway of DOG J1406+0102 would join an over-massive BH trajectory and 3) radio-loud DOGs can provide a significant contribution to the high-energy ($\gtrsim$ 100 TeV) cosmic neutrino background if we assume DOG J1406+0102 as a representative of radio-loud DOGs. △ Less

Submitted 9 March, 2023; originally announced March 2023.

Comments: 10 pages, 5 figures, submitted to ApJL

The eROSITA Final Equatorial-Depth Survey (eFEDS): Host-galaxy Demographics of X-ray AGNs with Subaru Hyper Suprime-Cam

Authors: Junyao Li, John D. Silverman, Andrea Merloni, Mara Salvato, Johannes Buchner, Andy Goulding, Teng Liu, Riccardo Arcodia, Johan Comparat, Xuheng Ding, Kohei Ichikawa, Masatoshi Imanishi, Toshihiro Kawaguchi, Lalitwadee Kawinwanichakij, Yoshiki Toba

Abstract: We investigate the physical properties, such as star-forming activity, disk vs. bulge nature, galaxy size, and obscuration of 3796 X-ray selected AGNs at $0.2<z<0.8$ in the eFEDS field. Using Subaru Hyper Suprime-Cam imaging data in the $grizy$ bands for SRG/eROSITA-detected AGNs, we measure the structural parameters for AGN host galaxies by performing a 2D AGN-host image decomposition. We then co… ▽ More We investigate the physical properties, such as star-forming activity, disk vs. bulge nature, galaxy size, and obscuration of 3796 X-ray selected AGNs at $0.2<z<0.8$ in the eFEDS field. Using Subaru Hyper Suprime-Cam imaging data in the $grizy$ bands for SRG/eROSITA-detected AGNs, we measure the structural parameters for AGN host galaxies by performing a 2D AGN-host image decomposition. We then conduct spectral energy distribution fitting to derive stellar mass and rest-frame colors for AGN hosts. We find that (1) AGNs can contribute significantly to the total optical light down to ${\rm log}\,L_{\rm X}\sim 42.5\ \rm erg\,s^{-1}$, thus ignoring the AGN component can significantly bias the structural measurements; (2) AGN hosts are predominately star-forming galaxies at ${\rm log}\,\mathcal{M}_\star \lesssim 11.3\ M_\odot$; (3) the bulk of AGNs (64%) reside in galaxies with significant stellar disks, while their host galaxies become increasingly bulge dominated and quiescent at ${\rm log}\,\mathcal{M}_\star \gtrsim 11.0\ M_\odot$; (4) the size-stellar mass relation of AGN hosts tends to lie between that of inactive star-forming and quiescent galaxies, suggesting that the physical mechanism responsible for building the central stellar density also efficiently fuel the black hole growth; (5) the hosts of X-ray unobscured AGNs are biased towards face-on systems and the average $E(B-V)/N_{\rm H}$ is similar to the galactic dust-to-gas ratio, suggesting that some of the obscuration of the nuclei could come from galaxy-scale gas and dust, which may partly account for (up to 30%) the deficiency of star-forming disks as host galaxies for the most massive AGNs. These results are consistent with a scenario in which the black hole and galaxy grow in mass while transform in structure and star-forming activity, as desired to establish the local scaling relations. △ Less

Submitted 23 February, 2023; originally announced February 2023.

Comments: 15 pages, 13 figures, submitted to MNRAS

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)

ALMA Lensing Cluster Survey: Properties of Millimeter Galaxies Hosting X-ray Detected Active Galactic Nuclei

Authors: Ryosuke Uematsu, Yoshihiro Ueda, Kotaro Kohno, Satoshi Yamada, Yoshiki Toba, Seiji Fujimoto, Bunyo Hatsukade, Hideki Umehata, Daniel Espada, Fengwu Sun, Georgios E. Magdis, Vasily Kokorev, Yiping Ao

Abstract: We report the multi-wavelength properties of millimeter galaxies hosting X-ray detected active galactic nuclei (AGNs) from the ALMA Lensing Cluster Survey (ALCS). ALCS is an extensive survey of well-studied lensing clusters with ALMA, covering an area of 133 arcmin$^2$ over 33 clusters with a 1.2 mm flux-density limit of ${\sim}$60 $\mathrm{μJy}$ ($1σ$). Utilizing the archival data of Chandra, we… ▽ More We report the multi-wavelength properties of millimeter galaxies hosting X-ray detected active galactic nuclei (AGNs) from the ALMA Lensing Cluster Survey (ALCS). ALCS is an extensive survey of well-studied lensing clusters with ALMA, covering an area of 133 arcmin$^2$ over 33 clusters with a 1.2 mm flux-density limit of ${\sim}$60 $\mathrm{μJy}$ ($1σ$). Utilizing the archival data of Chandra, we identify three AGNs at $z=$1.06, 2.09, and 2.84 among the 180 millimeter sources securely detected in the ALCS (of which 155 are inside the coverage of Chandra). The X-ray spectral analysis shows that two AGNs are not significantly absorbed ($\log N_{\mathrm{H}}/\mathrm{cm}^{-2} < 23$), while the other shows signs of moderate absorption ($\log N_{\mathrm{H}}/\mathrm{cm}^{-2}\sim 23.5$). We also perform spectral energy distribution (SED) modelling of X-ray to millimeter photometry. We find that our X-ray AGN sample shows both high mass accretion rates (intrinsic 0.5--8 keV X-ray luminosities of ${\sim}10^{\text{44--45}}\,\mathrm{erg\ s^{-1}}$) and star-formation rates (${\gtrsim}100\,M_{\odot}\,\mathrm{yr}^{-1}$). This demonstrates that a wide-area survey with ALMA and Chandra can selectively detect intense growth of both galaxies and supermassive black holes (SMBHs) in the high-redshift universe. △ Less

Submitted 23 January, 2023; originally announced January 2023.

Comments: 14 pages, 3 figures, 2 tables

Hard X-Ray to Radio Multiwavelength SED Analysis of Local U/LIRGs in GOALS Sample with Self-consistent AGN Model Including Polar-dust Component

Authors: Satoshi Yamada, Yoshihiro Ueda, Martín Herrera-Endoqui, Yoshiki Toba, Takamitsu Miyaji, Shoji Ogawa, Ryosuke Uematsu, Atsushi Tanimoto, Masatoshi Imanishi, Claudio Ricci

Abstract: We conduct a hard X-ray to radio multiwavelength spectral energy distribution (SED) decomposition for 57 local luminous and ultraluminous infrared galaxies (U/LIRGs) observed with Nuclear Spectroscopic Telescope Array and/or Swift/Burst Alert Telescope in GOALS (Armus et al. 2009) sample. We modify the latest SED-fitting code X-CIGALE by implementing the infrared (IR) CLUMPY model, allowing the mu… ▽ More We conduct a hard X-ray to radio multiwavelength spectral energy distribution (SED) decomposition for 57 local luminous and ultraluminous infrared galaxies (U/LIRGs) observed with Nuclear Spectroscopic Telescope Array and/or Swift/Burst Alert Telescope in GOALS (Armus et al. 2009) sample. We modify the latest SED-fitting code X-CIGALE by implementing the infrared (IR) CLUMPY model, allowing the multiwavelength study with the X-ray torus model (XCLUMPY) self-consistently. Adopting the torus parameters obtained by the X-ray fitting (Yamada et al. 2021), we estimate the properties of host galaxies, active galactic nucleus (AGN) tori, and polar dust. The star formation rates (SFRs) become larger with merger stage and most of them are above the main sequence. The SFRs are correlated with radio luminosity, indicating starburst emission is dominant in the radio band. Although polar-dust extinction is much smaller than torus extinction, the UV-to-IR (mainly IR) polar dust luminosities are $\sim$2 times larger than the torus ones. The polar-dust temperature decreases while the physical size, estimated by the temperature and dust sublimation radius, increases with AGN luminosity from a few tens of parsec (early mergers) to kiloparsec scales (late mergers), where the polar dust is likely the expanding (i.e., evolving) dusty outflows. The comparison between SFRs and intrinsic AGN luminosities suggests that the starbursts occur first and AGNs arise later, and overall their growth rates follow the simultaneous coevolution relation that can establish the local galaxy-SMBH mass relation. We confirm the coexistence of intense starbursts, AGNs, and large-scale outflows in late mergers, supporting a standard AGN feedback scenario. △ Less

Submitted 9 January, 2023; originally announced January 2023.

Comments: 84 pages (5 tables/29 figures in the main text and 8 tables/18 figures in the Appendix), accepted in ApJS

Massive molecular gas companions uncovered by VLA CO(1-0) observations of the $z$ = 5.2 radio galaxy TN J0924$-$2201

Authors: Kianhong Lee, Kotaro Kohno, Bunyo Hatsukade, Fumi Egusa, Takuji Yamashita, Malte Schramm, Kohei Ichikawa, Masatoshi Imanishi, Takuma Izumi, Tohru Nagao, Yoshiki Toba, Hideki Umehata

Abstract: We present Karl G. Jansky Very Large Array (VLA) K-band (19 GHz) observations of the redshifted CO(1-0) line emission toward the radio galaxy TN J0924$-$2201 at $z=5.2$, which is one of the most distant CO-detected radio galaxies. With the angular resolution of $\sim2''$, the CO(1-0) line emission is resolved into three clumps, within $\pm500$ km\,s$^{-1}$ relative to its redshift, where is determ… ▽ More We present Karl G. Jansky Very Large Array (VLA) K-band (19 GHz) observations of the redshifted CO(1-0) line emission toward the radio galaxy TN J0924$-$2201 at $z=5.2$, which is one of the most distant CO-detected radio galaxies. With the angular resolution of $\sim2''$, the CO(1-0) line emission is resolved into three clumps, within $\pm500$ km\,s$^{-1}$ relative to its redshift, where is determined by Ly$α$. We find that they locate off-center and 12-33 kpc away from the center of the host galaxy, which has counterparts in $HST$ $i$-band, $Spitzer$/IRAC and ALMA Band-6 (230 GHz; 1.3 mm). With the ALMA detection, we estimate $L_{\rm IR}$ and SFR of the host galaxy to be $(9.3\pm1.7)\times10^{11} L_{\odot}$ and $110\pm20$ $M_{\odot}\,\rm yr^{-1}$, respectively. We also derive the $3σ$ upper limit of $M_{\rm H_{2}}<1.3\times10^{10}$ $M_{\odot}$ at the host galaxy. The detected CO(1-0) line luminosities of three clumps, $L'_{\rm CO(1-0)}$ = (3.2-4.7)$\times10^{10}$ $\rm\,K\,km\,s^{-1}pc^{2}$, indicate the presence of three massive molecular gas reservoirs with $M_{\rm H_{2}}$ = (2.5-3.7)$\times10^{10}$ $M_{\odot}$, by assuming the CO-to-H$_{2}$ conversion factor $α_{\rm CO} = 0.8$ $M_{\rm \odot}\rm\,(K\,km\,s^{-1}pc^{2})^{-1}$, although the star formation rate (SFR) is not elevated because of the non-detection of ALMA 1.3 mm continuum (SFR $<$ 40 $M_\odot$ yr$^{-1}$). From the host galaxy, the nearest molecular gas clump labeled as clump A, is apparently aligning with the radio jet axis, showing the radio-CO alignment. The possible origin of these three clumps around TN J0924$-$2201 can be interpreted as merger, jet-induced metal enrichment and outflow. △ Less

Submitted 7 January, 2023; originally announced January 2023.

Comments: 14 pages, 9 figures. Accepted for publication in ApJ

The cosmic web of X-ray active galactic nuclei seen through the eROSITA Final Equatorial Depth Survey (eFEDS)

Authors: Johan Comparat, Wentao Luo, Andrea Merloni, Surhud More, Mara Salvato, Mirko Krumpe, Takamitsu Miyaji, William Brandt, Antonis Georgakakis, Masayuki Akiyama, Johannes Buchner, Tom Dwelly, Toshihiro Kawaguchi, Teng Liu, Tohru Nagao, Kirpal Nandra, John Silverman, Yoshiki Toba, Scott F. Anderson, Juna Kollmeier

Abstract: Which galaxies in the general population turn into active galactic nuclei (AGNs) is a keystone of galaxy formation and evolution. Thanks to SRG/eROSITA's contiguous 140 square degree pilot survey field, we constructed a large, complete, and unbiased soft X-ray flux-limited ($F_X>6.5\times 10^{-15}$ erg s$^{-1}$ cm$^{-2}$) AGN sample at low redshift, $0.05<z<0.55$. Two summary statistics, the clust… ▽ More Which galaxies in the general population turn into active galactic nuclei (AGNs) is a keystone of galaxy formation and evolution. Thanks to SRG/eROSITA's contiguous 140 square degree pilot survey field, we constructed a large, complete, and unbiased soft X-ray flux-limited ($F_X>6.5\times 10^{-15}$ erg s$^{-1}$ cm$^{-2}$) AGN sample at low redshift, $0.05<z<0.55$. Two summary statistics, the clustering using spectra from SDSS-V and galaxy-galaxy lensing with imaging from HSC, are measured and interpreted with halo occupation distribution and abundance matching models. Both models successfully account for the observations. We obtain an exceptionally complete view of the AGN halo occupation distribution. The population of AGNs is broadly distributed among halos with a mean mass of $3.9 _{- 2.4 }^{+ 2.0 }\times10^{12}M_\odot$. This corresponds to a large-scale halo bias of $b(z=0.34)= 0.99 ^{+0.08}_{-0.10}$. The central occupation has a large transition parameter, $σ_{\log_{10}(M)}=1.28\pm0.2$. The satellite occupation distribution is characterized by a shallow slope, $α_{\rm sat}=0.73\pm0.38$. We find that AGNs in satellites are rare, with $f_{\rm sat}<20\%$. Most soft X-ray-selected AGNs are hosted by central galaxies in their dark matter halo. A weak correlation between soft X-ray luminosity and large-scale halo bias is confirmed (3.3$σ$). We discuss the implications of environmental-dependent AGN triggering. This study paves the way toward fully charting, in the coming decade, the coevolution of X-ray AGNs, their host galaxies, and dark matter halos by combining eROSITA with SDSS-V, 4MOST, DESI, LSST, and \textit{Euclid} data. △ Less

Submitted 14 April, 2023; v1 submitted 3 January, 2023; originally announced January 2023.

Comments: 17 pages, 9 figures, 2 tables. A&A: Received: 19 December 2022 / Accepted: 21 March 2023

Journal ref: A&A 673, A122 (2023)

Metallicity-PAH Relation of MIR-selected Star-forming Galaxies in AKARI North Ecliptic Pole-wide Survey

Authors: Hyunjin Shim, Ho Seong Hwang, Woong-Seob Jeong, Yoshiki Toba, Minjin Kim, Dohyeong Kim, Hyunmi Song, Tetsuya Hashimoto, Takao Nakagawa, Ambra Nanni, William J. Pearson, Toshinobu Takagi

Abstract: We investigate the variation in the mid-infrared spectral energy distributions of 373 low-redshift ($z<0.4$) star-forming galaxies, which reflects a variety of polycyclic aromatic hydrocarbon (PAH) emission features. The relative strength of PAH emission is parameterized as $q_\mathrm{PAH}$, which is defined as the mass fraction of PAH particles in the total dust mass. With the aid of continuous m… ▽ More We investigate the variation in the mid-infrared spectral energy distributions of 373 low-redshift ($z<0.4$) star-forming galaxies, which reflects a variety of polycyclic aromatic hydrocarbon (PAH) emission features. The relative strength of PAH emission is parameterized as $q_\mathrm{PAH}$, which is defined as the mass fraction of PAH particles in the total dust mass. With the aid of continuous mid-infrared photometric data points covering 7-24$μ$m and far-infrared flux densities, $q_\mathrm{PAH}$ values are derived through spectral energy distribution fitting. The correlation between $q_\mathrm{PAH}$ and other physical properties of galaxies, i.e., gas-phase metallicity ($12+\mathrm{log(O/H)}$), stellar mass, and specific star-formation rate (sSFR) are explored. As in previous studies, $q_\mathrm{PAH}$ values of galaxies with high metallicity are found to be higher than those with low metallicity. The strength of PAH emission is also positively correlated with the stellar mass and negatively correlated with the sSFR. The correlation between $q_\mathrm{PAH}$ and each parameter still exists even after the other two parameters are fixed. In addition to the PAH strength, the application of metallicity-dependent gas-to-dust mass ratio appears to work well to estimate gas mass that matches the observed relationship between molecular gas and physical parameters. The result obtained will be used to calibrate the observed PAH luminosity-total infrared luminosity relation, based on the variation of MIR-FIR SED, which is used in the estimation of hidden star formation. △ Less

Submitted 8 December, 2022; originally announced December 2022.

Comments: 19 pages, 9 figures, AJ, in press

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

Spitzer/IRS full spectral modeling to characterize mineralogical properties of silicate dust in heavily obscured AGNs

Authors: T. Tsuchikawa, H. Kaneda, S. Oyabu, T. Kokusho, H. Kobayashi, Y. Toba

Abstract: Mid-Infrared (IR) silicate dust bands observed in heavily obscured active galactic nuclei (AGNs) include information on the mineralogical properties of silicate dust. We aim to investigate the mineralogical picture of the circumnuclear region of heavily obscured AGNs to reveal obscured AGN activities through the picture. In our previous study Tsuchikawa et al. (2021), we investigated the propertie… ▽ More Mid-Infrared (IR) silicate dust bands observed in heavily obscured active galactic nuclei (AGNs) include information on the mineralogical properties of silicate dust. We aim to investigate the mineralogical picture of the circumnuclear region of heavily obscured AGNs to reveal obscured AGN activities through the picture. In our previous study Tsuchikawa et al. (2021), we investigated the properties of silicate dust in heavily obscured AGNs focusing on the mineralogical composition and the crystallinity with Spitzer/IRS 5.3-12 micron spectra. In this study, we model the full-range Spitzer/IRS 5-30 micron spectra of 98 heavily obscured AGNs using a one-dimensional radiative transfer calculation with four dust species in order to evaluate wider ranges of the properties of silicate dust more reliably. Comparing fitting results between four dust models with a different size and porosity, 95 out of the 98 galaxies prefer a porous silicate dust model without micron-sized large grains. The pyroxene mass fraction and the crystallinity are overall consistent with but significantly different from the previous results for the individual galaxies. The pyroxene-poor composition, small dust size and high porosity are similar to newly formed dust around mass-loss stars as seen in our Galaxy, which presumably originates from the recent circumnuclear starburst activity. The high crystallinity on average suggests dust processing induced by AGN activities. △ Less

Submitted 22 October, 2022; originally announced October 2022.

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

A machine learning approach to assessing the presence of substructure in quasar host galaxies using the Hyper Suprime-Cam Subaru Strategic Program

Authors: Chris Nagele, John D. Silverman, Tilman Hartwig, Junyao Li, Connor Bottrell, Xuheng Ding, Yoshiki Toba

Abstract: The conditions under which galactic nuclear regions become active are largely unknown, although it has been hypothesized that secular processes related to galaxy morphology could play a significant role. We investigate this question using optical i-band images of 3096 SDSS quasars and galaxies at 0.3<z<0.6 from the Hyper Suprime-Cam Subaru Strategic Program, which possess a unique combination of a… ▽ More The conditions under which galactic nuclear regions become active are largely unknown, although it has been hypothesized that secular processes related to galaxy morphology could play a significant role. We investigate this question using optical i-band images of 3096 SDSS quasars and galaxies at 0.3<z<0.6 from the Hyper Suprime-Cam Subaru Strategic Program, which possess a unique combination of area, depth and resolution, allowing the use of residual images, after removal of the quasar and smooth galaxy model, to investigate internal structural features. We employ a variational auto-encoder which is a generative model that acts as a form of dimensionality reduction. We analyze the lower dimensional latent space in search of features which correlate with nuclear activity. We find that the latent space does separate images based on the presence of nuclear activity which appears to be associated with more pronounced components (i.e., arcs, rings and bars) as compared to a matched control sample of inactive galaxies. These results suggest the importance of secular processes, and possibly mergers (by their remnant features) in activating or sustaining black hole growth. Our study highlights the breadth of information available in ground-based imaging taken under optimal seeing conditions and having accurate characterization of the point spread function (PSF) thus demonstrating future science to come from the Rubin Observatory. △ Less

Submitted 17 April, 2023; v1 submitted 21 October, 2022; originally announced October 2022.

Journal ref: 2023 ApJ 947 30

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

XXL-HSC: The link between AGN activity and star formation in the Early Universe ($z\geqslant3.5$)

Authors: E. Pouliasis, G. Mountrichas, I. Georgantopoulos, A. Ruiz, R. Gilli, E. Koulouridis, M. Akiyama, Y. Ueda, C. Garrel, T. Nagao, S. Paltani, M. Pierre, Y. Toba, C. Vignali

Abstract: In this work, we aimed at investigating the star formation rate of active galactic nuclei host galaxies in the early Universe. To this end, we constructed a sample of 149 luminous ($\rm L_{2-10keV} > 10^{44}\,erg\,s^{-1}$) X-ray AGNs at $\rm z \geq3.5$ selected in three fields with different depths and observed areas (Chandra COSMOS Legacy survey, XMM-XXL North and eFEDS). We built their spectral… ▽ More In this work, we aimed at investigating the star formation rate of active galactic nuclei host galaxies in the early Universe. To this end, we constructed a sample of 149 luminous ($\rm L_{2-10keV} > 10^{44}\,erg\,s^{-1}$) X-ray AGNs at $\rm z \geq3.5$ selected in three fields with different depths and observed areas (Chandra COSMOS Legacy survey, XMM-XXL North and eFEDS). We built their spectral energy distributions (SED) using available multi-wavelength photometry from X-rays up to far-IR. Then, we estimated the stellar mass, M$_{*}$, and the SFR of the AGNs using the X-CIGALE SED fitting algorithm. After applying several quality criteria, we ended up with 89 high-z sources. More than half (55\%) of the X-ray sample have spectroscopic redshifts. Based on our analysis, our high-z X-ray AGNs live in galaxies with median $\rm M_{*}=5.6 \times10^{10}~M_\odot$ and $\rm SFR_{*}\approx240\,M_\odot yr^{-1}$. The majority of the high-z sources ($\sim89$\%) were found inside or above the main sequence (MS) of star-forming galaxies. Estimation of the normalised SFR, $\rm SFR_{NORM}$, defined as the ratio of the SFR of AGNs to the SFR of MS galaxies, showed that the SFR of AGNs is enhanched by a factor of $\sim 1.8$ compared to non-AGN star-forming systems. Combining our results with previous studies at lower redshifts, we confirmed that $\rm SFR_{NORM}$ does not evolve with redshift. Using the specific BHAR (i.e., $\rm L_X$ divided by $\rm M_{*}$), $\rm λ_{BHAR}$, that can be used as a tracer of the Eddington ratio, we found that the bulk of AGNs that lie inside or above the MS have higher specific accretion rates compared to sources below the MS. Finally, we found indications that the SFR of the most massive AGN host galaxies ($\rm log\,(M_{*}/ M_\odot) >10^{11.5-12}$) remains roughly constant as a function of M$_*$, in agreement with the SFR of MS star-forming galaxies. △ Less

Submitted 28 June, 2022; originally announced June 2022.

Comments: Accepted for publication in A&A. 19 pages, 12 figures (+6 in Appendix), 4 tables

Journal ref: A&A 667, A56 (2022)

A Wide and Deep Exploration of Radio Galaxies with Subaru HSC (WERGS). VII. Redshift Evolution of Radio Galaxy Environments at z=0.3-1.4

Authors: Hisakazu Uchiyama, Takuji Yamashita, Tohru Nagao, Kohei Ichikawa, Yoshiki Toba, Shogo Ishikawa, Mariko Kubo, Masaru Kajisawa, Toshihiro Kawaguchi, Nozomu Kawakatu, Chien-Hsiu Lee, Akatoki Noboriguchi

Abstract: We examine the redshift evolution of density environments around 2,163 radio galaxies with the stellar masses of $\sim10^{9}-10^{12} M_\odot$ between redshifts of $z=0.3-1.4$, based on the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) and Faint Images of the Radio Sky at Twenty-cm (FIRST). We use the $k$-nearest neighbor method to measure the local galaxy number density around our radio gal… ▽ More We examine the redshift evolution of density environments around 2,163 radio galaxies with the stellar masses of $\sim10^{9}-10^{12} M_\odot$ between redshifts of $z=0.3-1.4$, based on the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) and Faint Images of the Radio Sky at Twenty-cm (FIRST). We use the $k$-nearest neighbor method to measure the local galaxy number density around our radio galaxy sample. We find that the overdensities of the radio galaxies are weakly but significantly anti-correlated with redshift. This is consistent with the known result that the relative abundance of less-massive radio galaxies increases with redshift, because less-massive radio galaxies reside in relatively low density regions. Massive radio galaxies with stellar mass of $M_* >10^{11}M_\odot$ are found in high density environments compared with the control sample galaxies with radio-non-detection and matched-stellar-mass. Less-massive radio galaxies with $M_* <10^{11}M_\odot$ reside in average density environments. The fraction of the radio galaxies associated with the neighbors within a typical major merger scale, $<70$ kpc, is higher than (comparable to) that of the control galaxies at $M_* >10^{11}M_\odot$ ($M_* <10^{11}M_\odot$). We also find that the local densities around the radio galaxies are anti-correlated with the radio luminosities and black hole mass accretion rates at fixed stellar mass. These findings suggest that massive radio galaxies have matured through galaxy mergers in the past, and have supermassive black holes whose mass accretion almost ceased at $z>1.4$, while less-massive radio galaxies undergo active accretion just at this epoch, as they have avoided such merger events. △ Less

Submitted 14 June, 2022; originally announced June 2022.

Comments: 14 pages, 8 figures, ApJ accepted

Origin and evolution of dust-obscured galaxies in galaxy mergers

Authors: Naomichi Yutani, Yoshiki Toba, Shunsuke Baba, Keiichi Wada

Abstract: Dust Obscured Galaxies (DOGs), which are observationally characterized as faint in the optical and bright in the infrared, are the final stage of galaxy mergers and are essential objects in the evolution of galaxies and active galactic nuclei (AGNs). However, the relationship between torus-scale gas dynamics around AGNs and DOGs lifetime remain unclear. We obtained evolution of the spectral energy… ▽ More Dust Obscured Galaxies (DOGs), which are observationally characterized as faint in the optical and bright in the infrared, are the final stage of galaxy mergers and are essential objects in the evolution of galaxies and active galactic nuclei (AGNs). However, the relationship between torus-scale gas dynamics around AGNs and DOGs lifetime remain unclear. We obtained evolution of the spectral energy distributions (SEDs) of a galaxy merger system with AGN feedback, from post-processed pseudo-observations based on an N-body/Smoothed Particle Hydrodynamics (SPH) simulation. We focused on a late stage merger of two identical galaxies with a supermassive black hole (SMBH) of 10$^8$ M$_\odot$. We found that the infrared luminosity of the system reaches ultra- and hyper-luminous infrared galaxy classes (10$^{12}$ and 10$^{13}$ L$_\odot$, respectively). The DOGs phase corresponds to a state in which the AGNs are buried in dense gas and dust, with the infrared luminosity exceeding 3.3 $\times$ 10$^{12}$ L$_\odot$. We also identified the sub-categories of DOGs, namely bump and power-law DOGs from the SEDs and their evolution. The bump DOGs tend to evolve to power-law DOGs on several Myrs. We found that contribution from the hot dust around the nucleus in the infrared radiation is essential for identifying the system as a power-law DOG; the gas and dust distribute non-spherically around the nucleus, therefore, the observed properties of DOGs depend on the viewing angle. In our model, the lifetime of merger-driven DOGs is less than 4 Myrs, suggesting that the observed DOGs phase is a brief aspect of galaxy mergers. △ Less

Submitted 6 August, 2022; v1 submitted 1 May, 2022; originally announced May 2022.

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

Multi-wavelength properties of 850-$μ$m selected sources from the North Ecliptic Pole SCUBA-2 survey

Authors: H. Shim, D. Lee, Y. Kim, D. Scott, S. Serjeant, Y. Ao, L. Barrufet, S. C. Chapman, D. Clements, C. J. Conselice, T. Goto, T. R. Greve, H. S. Hwang, M. Im, W. -S. Jeong, H. K. Kim, M. Kim, S. J. Kim, A. K. H. Kong, M. P. Koprowski, M. A. Malkan, M. Michalowski, C. Pearson, H. Seo, T. Takagi , et al. (3 additional authors not shown)

Abstract: We present the multi-wavelength counterparts of 850-$μ$m selected submillimetre sources over a 2-deg$^2$ field centred on the North Ecliptic Pole. In order to overcome the large beam size (15 arcsec) of the 850-$μ$m images, deep optical to near-infrared (NIR) photometric data and arcsecond-resolution 20-cm images are used to identify counterparts of submillimetre sources. Among 647 sources, we ide… ▽ More We present the multi-wavelength counterparts of 850-$μ$m selected submillimetre sources over a 2-deg$^2$ field centred on the North Ecliptic Pole. In order to overcome the large beam size (15 arcsec) of the 850-$μ$m images, deep optical to near-infrared (NIR) photometric data and arcsecond-resolution 20-cm images are used to identify counterparts of submillimetre sources. Among 647 sources, we identify 514 reliable counterparts for 449 sources (69 per cent in number), based either on probabilities of chance associations calculated from positional offsets or offsets combined with the optical-to-NIR colours. In the radio imaging, the fraction of 850-$μ$m sources having multiple counterparts is 7 per cent. The photometric redshift, infrared luminosity, stellar mass, star-formation rate (SFR), and the AGN contribution to the total infrared luminosity of the identified counterparts are investigated through spectral energy distribution fitting. The SMGs are infrared-luminous galaxies at an average $\langle z\rangle=2.5$ with $\mathrm{log}_{10} (L_\mathrm{IR}/\mathrm{L}_\odot)=11.5-13.5$, with a mean stellar mass of $\mathrm{log}_{10} (M_\mathrm{star}/\mathrm{M}_\odot)=10.90$ and SFR of $\mathrm{log}_{10} (\mathrm{SFR/M_\odot\,yr^{-1}})=2.34$. The SMGs show twice as large SFR as galaxies on the star-forming main sequence, and about 40 per cent of the SMGs are classified as objects with bursty star formation. At $z\ge4$, the contribution of AGN luminosity to total luminosity for most SMGs is larger than 30 per cent. The FIR-to-radio correlation coefficient of SMGs is consistent with that of main-sequence galaxies at $z\simeq2$. △ Less

Submitted 24 April, 2022; originally announced April 2022.

Comments: 24 pages, 16 figures. Accepted for publication in MNRAS

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