Showing 1–50 of 97 results for author: Kawaguchi, T

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)

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

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

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

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

Submitted 15 November, 2023; originally announced November 2023.

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

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

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

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

Submitted 18 October, 2023; originally announced October 2023.

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

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

EMPRESS. XIV. Strong High Ionization Lines of Young Galaxies at $z=0-8$: Ionizing Spectra Consistent with the Intermediate Mass Black Holes with $M_{\rm BH}\sim 10^3-10^6\ M_\odot$

Authors: Shun Hatano, Masami Ouchi, Hiroya Umeda, Kimihiko Nakajima, Toshihiro Kawaguchi, Yuki Isobe, Shohei Aoyama, Kuria Watanabe, Yuichi Harikane, Haruka Kusakabe, Akinori Matsumoto, Takashi J. Moriya, Moka Nishigaki, Yoshiaki Ono, Masato Onodera, Yuma Sugahara, Akihiro Suzuki, Yi Xu, Yechi Zhang

Abstract: We present ionizing spectra estimated at 13.6--100 eV for ten dwarf galaxies with strong high ionization lines of He {\sc {ii}}$λ$4686 and [Ne {\sc{v}}]$λ$3426 ([Ne {\sc{iv}}]$λ$2424) at $z=0$ ($z=8$) that are identified in our Keck/LRIS spectroscopy and the literature (the JWST ERO program). With the flux ratios of these high ionization lines and $>10$ low-ionization lines of hydrogen, helium, ox… ▽ More We present ionizing spectra estimated at 13.6--100 eV for ten dwarf galaxies with strong high ionization lines of He {\sc {ii}}$λ$4686 and [Ne {\sc{v}}]$λ$3426 ([Ne {\sc{iv}}]$λ$2424) at $z=0$ ($z=8$) that are identified in our Keck/LRIS spectroscopy and the literature (the JWST ERO program). With the flux ratios of these high ionization lines and $>10$ low-ionization lines of hydrogen, helium, oxygen, neon, and sulfur, we determine ionizing spectra consisting of stellar and non-thermal power-law radiation by photoionization modeling with free parameters of nebular properties including metallicity and ionization parameter, cancelling out abundance ratio differences. We find that all of the observed flux ratios are well reproduced by the photoinization models with the power law index $α_{\rm EUV}$ of $α_{\rm EUV}\sim (-1)-0$ and the luminosity $L_{\rm EUV}$ of $L_{\rm EUV}\sim 10^{40}-10^{42}$ erg s$^{-1}$ at $\sim 55-100$ eV for six galaxies, while four galaxies include large systematics in $α_{\rm EUV}$ caused by stellar radiation contamination. We then compare $α_{\rm EUV}$ and $L_{\rm EUV}$ of these six galaxies with those predicted by the black hole (BH) accretion disk models, and find that $α_{\rm EUV}$ and $L_{\rm EUV}$ are similar to those of the intermediate mass black holes (IMBHs) in BH accretion disk models {albeit with possibilities of the other scenarios.} Confirming these results with a known IMBH having a mass $M_{\rm BH}$ of $M_{\rm BH}=10^{5.75} \ M_\odot$, we find that four local galaxies and one $z=7.665$ galaxy have ionizing spectra consistent with those of IMBHs with $M_{\rm BH} \sim 10^3-10^5 \ M_\odot$. △ Less

Submitted 10 March, 2024; v1 submitted 3 May, 2023; originally announced May 2023.

Comments: Accepted for publication in ApJ and 25 pages, 12 figures

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)

Active Massive Black Hole Found in the Young Star-Forming Dwarf Galaxy SBS 0335-052E

Authors: Shun Hatano, Masami Ouchi, Kimihiko Nakajima, Toshihiro Kawaguchi, Mitsuru Kokubo, Satoshi Kikuta, Nozomu Tominaga, Yi Xu, Kuria Watanabe, Yuichi Harikane, Yuki Isobe, Akinori Matsumoto, Moka Nishigaki, Yoshiaki Ono, Masato Onodera, Yuma Sugahara, Hiroya Umeda, Yechi Zhang

Abstract: Since the late 20th century, young star-forming dwarf galaxies with extremely low heavy-element abundances have been identified in the local universe\cite{1972ApJ...173...25S, 1990Natur.343..238I, 2022ApJS..262....3N}. Such a population of galaxies, represented by SBS 0335-052E, is intensively studied as a laboratory of galaxy formation, and interpreted as galaxies with hard emission and energetic… ▽ More Since the late 20th century, young star-forming dwarf galaxies with extremely low heavy-element abundances have been identified in the local universe\cite{1972ApJ...173...25S, 1990Natur.343..238I, 2022ApJS..262....3N}. Such a population of galaxies, represented by SBS 0335-052E, is intensively studied as a laboratory of galaxy formation, and interpreted as galaxies with hard emission and energetic outflows driven by massive stars\cite{1990Natur.343..238I, 2004ApJ...606..213T, 2009AA...503...61I}. Here we report the temporal flux variability of SBS 0335-052E in $3-4 \ μ$m bands on timescale of months to years with dimming and brightening up to 50\% over 12 years. This is a clear signature of dust torus emission of an active massive black hole in SBS 0335-052E. The deep optical spectrum reveals a very broad component with $1.24\ (\pm 0.01) \times 10^{4}$ km s$^{-1}$ in H$α$ emission, suggesting emission originated from broad line regions around the massive black hole. The black hole mass is estimated to be $\sim 100$ thousand solar masses from the spectral energy distribution. This is the first example of the massive black hole existence in a long-considered young star-forming dwarf galaxy. △ Less

Submitted 7 April, 2023; originally announced April 2023.

Comments: 32 pages, 6 figures, 3 tables

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)

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)

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

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

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

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

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

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

Submitted 24 November, 2021; originally announced November 2021.

Comments: Accepted for publication in The Astrophysical Journal Supplement Series

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

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

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

Submitted 21 September, 2021; originally announced September 2021.

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

The eROSITA Final Equatorial-Depth Survey (eFEDS): A multiwavelength view of WISE mid-infrared galaxies/active galactic nuclei

Authors: Yoshiki Toba, Teng Liu, Tanya Urrutia, Mara Salvato, Junyao Li, Yoshihiro Ueda, Marcella Brusa, Naomichi Yutani, Keiichi Wada, Atsushi J. Nishizawa, Johannes Buchner, Tohru Nagao, Andrea Merloni, Masayuki Akiyama, Riccardo Arcodia, Bau-Ching Hsieh, Kohei Ichikawa, Masatoshi Imanishi, Kaiki T. Inoue, Toshihiro Kawaguchi, Georg Lamer, Kirpal Nandra, John D. Silverman, Yuichi Terashima

Abstract: We investigate the physical properties--such as the stellar mass, SFR, IR luminosity, X-ray luminosity, and hydrogen column density--of MIR galaxies and AGN at $z < 4$ in the 140 deg$^2$ field observed by SRG/eROSITA through the eFEDS survey. By cross-matching the WISE 22 $μ$m (W4)-detected sample and the eFEDS X-ray point-source catalog, we find that 692 extragalactic objects are detected by eROS… ▽ More We investigate the physical properties--such as the stellar mass, SFR, IR luminosity, X-ray luminosity, and hydrogen column density--of MIR galaxies and AGN at $z < 4$ in the 140 deg$^2$ field observed by SRG/eROSITA through the eFEDS survey. By cross-matching the WISE 22 $μ$m (W4)-detected sample and the eFEDS X-ray point-source catalog, we find that 692 extragalactic objects are detected by eROSITA. We have compiled a multiwavelength dataset. We have also performed (i) an X-ray spectral analysis, (ii) SED fitting using X-CIGALE, (iii) 2D image-decomposition analysis using Subaru HSC images, and (iv) optical spectral fitting with QSFit to investigate the AGN and host-galaxy properties. For 7,088 WISE W4 objects that are undetected by eROSITA, we have performed an X-ray stacking analysis to examine the typical physical properties of these X-ray faint and/or probably obscured objects. We find that (i) 82% of the eFEDS-W4 sources are classified as X-ray AGN with $\log\,L_{\rm X} >$ 42 erg s$^{-1}$; (ii) 67% and 24% of the objects have $\log\,(L_{\rm IR}/L_{\odot}) > 12$ and 13, respectively; (iii) the relationship between $L_{\rm X}$ and the 6 $μ$m luminosity is consistent with that reported in previous works; and (iv) the relationship between the Eddington ratio and $N_{\rm H}$ for the eFEDS-W4 sample and a comparison with a model prediction from a galaxy-merger simulation indicates that approximately 5% of the eFEDS-W4 sources in our sample are likely to be in an AGN-feedback phase, in which strong radiation pressure from the AGN blows out the surrounding material from the nuclear region. Thanks to the wide area coverage of eFEDS, we have been able to constrain the ranges of the physical properties of the WISE W4 sample of AGNs at $z < 4$, providing a benchmark for forthcoming studies on a complete census of MIR galaxies selected from the full-depth eROSITA all-sky survey. △ Less

Submitted 24 August, 2021; v1 submitted 28 June, 2021; originally announced June 2021.

Comments: 18 pages, 19 figures, and 3 tables, accepted to appear on A&A, Special Issue: The Early Data Release of eROSITA and Mikhail Pavlinsky ART-XC on the SRG Mission

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

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

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

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

Submitted 21 May, 2021; originally announced May 2021.

Comments: 36 pages, 21 figures, 4 tables

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

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

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

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

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

The eROSITA Final Equatorial-Depth Survey (eFEDS): An X-ray bright, extremely luminous infrared galaxy at z = 1.87

Authors: Yoshiki Toba, Marcella Brusa, Teng Liu, Johannes Buchner, Yuichi Terashima, Tanya Urrutia, Mara Salvato, Masayuki Akiyama, Riccardo Arcodia, Andy D. Goulding, Yuichi Higuchi, Kaiki T. Inoue, Toshihiro Kawaguchi, Georg Lamer, Andrea Merloni, Tohru Nagao, Yoshihiro Ueda, Kirpal Nandra

Abstract: In this study, we investigate the X-ray properties of WISE J090924.01+000211.1 (WISEJ0909+0002), an extremely luminous infrared (IR) galaxy (ELIRG) at $z_{\rm spec}$= 1.871 in the eROSITA final equatorial depth survey (eFEDS). WISEJ0909+0002 is a WISE 22 $μ$m source, located in the GAMA-09 field, detected by eROSITA during the performance and verification phase. The corresponding optical spectrum… ▽ More In this study, we investigate the X-ray properties of WISE J090924.01+000211.1 (WISEJ0909+0002), an extremely luminous infrared (IR) galaxy (ELIRG) at $z_{\rm spec}$= 1.871 in the eROSITA final equatorial depth survey (eFEDS). WISEJ0909+0002 is a WISE 22 $μ$m source, located in the GAMA-09 field, detected by eROSITA during the performance and verification phase. The corresponding optical spectrum indicates that this object is a type-1 active galactic nucleus (AGN). Observations from eROSITA combined with Chandra and XMM-Newton archival data indicate a very luminous ($L$ (2--10 keV) = ($2.1 \pm 0.2) \times 10^{45}$ erg s$^{-1}$) unobscured AGN with a power-law photon index of $Γ$ = 1.73$_{-0.15}^{+0.16}$, and an absorption hydrogen column density of $\log\,(N_{\rm H}/{\rm cm}^{-2}) < 21.0$. The IR luminosity was estimated to be $L_{\rm IR}$ = (1.79 $\pm$ 0.09) $\times 10^{14}\, L_{\odot}$ from spectral energy distribution modeling based on 22 photometric data (X-ray to far-IR) with X-CIGALE, which confirmed that WISEJ0909+0002 is an ELIRG. A remarkably high $L_{\rm IR}$ despite very low $N_{\rm H}$ would indicate that we are witnessing a short-lived phase in which hydrogen gas along the line of sight is blown outwards, whereas warm and hot dust heated by AGNs still exist. As a consequence of eROSITA all-sky survey, $6.8_{-5.6}^{+16}\times 10^2$ such X-ray bright ELIRGs are expected to be discovered in the entire extragalactic sky ($|b| > 10^\circ$). This can potentially be the key population to constrain the bright-end of IR luminosity functions. △ Less

Submitted 6 April, 2021; originally announced April 2021.

Comments: 10 pages, 5 figures, and 3 tables, accepted for publication in A&A Letters (special Issue: First science highlights from SRG/eROSITA)

Journal ref: A&A 649, L11 (2021)

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

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

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

Submitted 4 January, 2021; originally announced January 2021.

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

Destruction of the central black hole gas reservoir through head-on galaxy collisions

Authors: Yohei Miki, Masao Mori, Toshihiro Kawaguchi

Abstract: A massive black hole exists in almost every galaxy. They occasionally radiate a vast amount of light by releasing gravitational energy of accreting gas, with a cumulative active period of only a few $10^8$ years, so-called the duty cycle of the Active Galactic Nuclei. Namely, many galaxies today host a starving massive black hole. Although galaxy collisions have been thought to enhance nucleus act… ▽ More A massive black hole exists in almost every galaxy. They occasionally radiate a vast amount of light by releasing gravitational energy of accreting gas, with a cumulative active period of only a few $10^8$ years, so-called the duty cycle of the Active Galactic Nuclei. Namely, many galaxies today host a starving massive black hole. Although galaxy collisions have been thought to enhance nucleus activity, the origin of the duty cycle, especially the shutdown process, is a still critical issue. Here we show that galaxy collisions are also capable of suppressing black hole fueling by using an analytic model and three-dimensional hydrodynamic simulations, applying the well-determined parameter sets for the galactic collision in the Andromeda galaxy. Our models demonstrate that a central collision of galaxies can strip the torus-shaped gas surrounding the massive black hole, the putative fueling source. The derived condition for switching-off the black hole fueling indicates that a significant fraction of currently bright nuclei can become inactive, reminiscent of fading/dying active nucleus phenomena associated with galaxy merging events. Galaxy collisions may therefore be responsible both for switching-off and turning-on the nucleus activity, depending on the collision orbit (head-on or far-off-centre). △ Less

Submitted 25 July, 2021; v1 submitted 10 December, 2020; originally announced December 2020.

Comments: 41 pages, 10 figures, published in Nature Astronomy (read-only access is freely available at https://rdcu.be/cebgO)

Journal ref: Nature Astronomy 5, 478-484 (2021)

NuSTAR Non-detection of a Faint Active Galactic Nucleus in an Ultraluminous IR Galaxy with Kpc-scale Fast Wind

Authors: Xiaoyang Chen, Kohei Ichikawa, Hirofumi Noda, Taiki Kawamuro, Toshihiro Kawaguchi, Yoshiki Toba, Masayuki Akiyama

Abstract: Large-scale outflows are generally considered as a possible evidence that active galactic nuclei (AGNs) can severely affect their host galaxies. Recently an ultraluminous IR galaxy (ULIRG) at $z=0.49$, AKARI J0916248+073034, was found to have a galaxy-scale [OIII] $λ$5007 outflow with one of the highest energy-ejection rates at $z<1.6$. However, the central AGN activity estimated from its torus mi… ▽ More Large-scale outflows are generally considered as a possible evidence that active galactic nuclei (AGNs) can severely affect their host galaxies. Recently an ultraluminous IR galaxy (ULIRG) at $z=0.49$, AKARI J0916248+073034, was found to have a galaxy-scale [OIII] $λ$5007 outflow with one of the highest energy-ejection rates at $z<1.6$. However, the central AGN activity estimated from its torus mid-IR (MIR) radiation is weak relative to the luminous [OIII] emission. In this work we report the first NuSTAR hard X-ray follow-up of this ULIRG to constrain its current AGN luminosity. The intrinsic 2-10 keV luminosity shows a 90% upper-limit of $3.0\times10^{43}$ erg s$^{-1}$ assuming Compton-thick obscuration ($N_{\rm H}=1.5\times10^{24}$ cm$^{-2}$), which is only 3.6% of the luminosity expected from the extinction corrected [OIII] luminosity. With the NuSTAR observation, we succeed to identify that this ULIRG has a most extreme case of X-ray deficit among local ULIRGs. A possible scenario to explain the drastic declining in both of the corona (X-ray) and torus (MIR) is that the primary radiation from the AGN accretion disk is currently in a fading status, as a consequence of a powerful nuclear wind suggested by its powerful ionized outflow in the galaxy scale. △ Less

Submitted 13 November, 2020; originally announced November 2020.

Comments: 8 pages, 2 figures, accepted for publication in The Astrophysical Journal Letters

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

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

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

Submitted 1 October, 2020; originally announced October 2020.

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

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

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

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

Submitted 14 September, 2020; originally announced September 2020.

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

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

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

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

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

Comments: Accepted for publication in PASJ

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

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

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

Submitted 16 July, 2020; originally announced July 2020.

Comments: Accepted for publication in PASJ

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

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

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

Submitted 10 July, 2020; originally announced July 2020.

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

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

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

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

Submitted 8 June, 2020; originally announced June 2020.

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

Revisiting Soltan's argument based on a semi-analytical model for galaxy and black hole evolution

Authors: Hikari Shirakata, Toshihiro Kawaguchi, Takashi Okamoto, Masahiro Nagashima, Taira Oogi

Abstract: We show the significance of the super-Eddington accretion for the cosmic growth of supermassive black holes (SMBHs) with a semi-analytical model for galaxy and black hole evolution. The model explains various observed properties of galaxies and active galactic nuclei at a wide redshift range. By tracing the growth history of individual SMBHs, we find that the fraction of the SMBH mass acquired dur… ▽ More We show the significance of the super-Eddington accretion for the cosmic growth of supermassive black holes (SMBHs) with a semi-analytical model for galaxy and black hole evolution. The model explains various observed properties of galaxies and active galactic nuclei at a wide redshift range. By tracing the growth history of individual SMBHs, we find that the fraction of the SMBH mass acquired during the super-Eddington accretion phases to the total SMBH mass becomes larger for less massive black holes and at higher redshift. Even at z = 0, SMBHs with > 1e+9 Msun have acquired more than 50% of their mass by super-Eddington accretions, which is apparently inconsistent with classical Soltan's argument. However, the mass-weighted radiation efficiency of SMBHs with > 1e+8 Msun obtained with our model, is about 0.08 at z = 0, which is consistent with Soltan's argument within the observational uncertainties. We, therefore, conclude that Soltan's argument cannot reject the possibility that SMBHs are grown mainly by super-Eddington accretions. △ Less

Submitted 3 June, 2020; originally announced June 2020.

Comments: 10 pages, 6 figures. Accepted for publication in ApJ

Application of X-Ray Clumpy Torus Model (XCLUMPY) to 10 Obscured Active Galactic Nuclei Observed with Suzaku and NuSTAR

Authors: Atsushi Tanimoto, Yoshihiro Ueda, Hirokazu Odaka, Shoji Ogawa, Satoshi Yamada, Toshihiro Kawaguchi, Kohei Ichikawa

Abstract: We apply XCLUMPY, an X-ray spectral model from a clumpy torus in an active galactic nucleus (AGN), to the broadband X-ray spectra of 10 obscured AGNs observed with both Suzaku and NuSTAR. The infrared spectra of these AGNs were analyzed with the CLUMPY code. Since XCLUMPY adopts the same clump distribution as that in the CLUMPY, we can directly compare the torus parameters obtained from the X-ray… ▽ More We apply XCLUMPY, an X-ray spectral model from a clumpy torus in an active galactic nucleus (AGN), to the broadband X-ray spectra of 10 obscured AGNs observed with both Suzaku and NuSTAR. The infrared spectra of these AGNs were analyzed with the CLUMPY code. Since XCLUMPY adopts the same clump distribution as that in the CLUMPY, we can directly compare the torus parameters obtained from the X-ray spectra and those from the infrared ones. The torus angular widths determined from the infrared spectra ($σ_{\mathrm{IR}}$) are systematically larger than those from the X-ray data ($σ_{\mathrm{X}}$); the difference ($σ_{\mathrm{IR}}-σ_{\mathrm{X}}$) correlates with the inclination angle determined from the X-ray spectrum. These results can be explained by the contribution from dusty polar outflows to the observed infrared flux, which becomes more significant at higher inclinations (more edge-on views). The ratio of the hydrogen column density and V-band extinction in the line of sight absorber shows large scatter ($\simeq$1 dex) around the Galactic value, suggesting that a significant fraction of AGNs have dust-rich circumnuclear environments. △ Less

Submitted 26 May, 2020; originally announced May 2020.

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

A Comparison of Properties of Quasars with and without Rapid Broad Absorption Line Variability

Authors: Takashi Horiuchi, Tomoki Morokuma, Toru Misawa, Hidekazu Hanayama, Toshihiro Kawaguchi

Abstract: We investigate the correlation between rest-frame UV flux variability of broad absorption line (BAL) quasars and their variability in BAL equivalent widths (EWs) in a various timescale from $<10$~days to a few years in the quasar rest-frame. We use the data sets of BAL EWs taken by the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project and photometric data taken by the intermediate P… ▽ More We investigate the correlation between rest-frame UV flux variability of broad absorption line (BAL) quasars and their variability in BAL equivalent widths (EWs) in a various timescale from $<10$~days to a few years in the quasar rest-frame. We use the data sets of BAL EWs taken by the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project and photometric data taken by the intermediate Palomar Transient Factory (iPTF) in $g$ and $R$-bands and the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS) in $grizy$ bands. Our results are summarized as below; (1) the distributions of flux variability versus BAL variability show weak, moderate, or a strong positive correlation, (2) there is no significant difference in flux variability amplitudes between BAL quasar with significant short timescale EW variability (called class S1) and without (class S2), (3) in all time scales considered in this paper, the class S1 quasars show systematically larger BAL variability amplitudes than those of the class S2 quasars, and (4) there are possible correlations between BAL variability and physical parameters of the quasars such as black hole masses (moderate positive), Eddington ratios, and accretion disk temperature (strong negative) in the class S2 quasars. These results indicate that the BAL variability requires changing in the ionizing continuum and an ancillary mechanism such as variability in X-ray shielding gas located at the innermost region of an accretion disk. △ Less

Submitted 4 April, 2020; originally announced April 2020.

Comments: 22 pages

Semi-analytic modelling of AGNs: auto-correlation function and halo occupation

Authors: Taira Oogi, Hikari Shirakata, Masahiro Nagashima, Takahiro Nishimichi, Toshihiro Kawaguchi, Takashi Okamoto, Tomoaki Ishiyama, Motohiro Enoki

Abstract: The spatial clustering of active galactic nuclei (AGNs) is considered to be one of the important diagnostics for the understanding of the underlying processes behind their activities complementary to measurements of the luminosity function (LF). We analyse the AGN clustering from a recent semi-analytic model performed on a large cosmological $N$-body simulation covering a cubic gigaparsec comoving… ▽ More The spatial clustering of active galactic nuclei (AGNs) is considered to be one of the important diagnostics for the understanding of the underlying processes behind their activities complementary to measurements of the luminosity function (LF). We analyse the AGN clustering from a recent semi-analytic model performed on a large cosmological $N$-body simulation covering a cubic gigaparsec comoving volume. We have introduced a new time-scale of gas accretion on to the supermassive black holes to account for the loss of the angular momentum on small scales, which is required to match the faint end of the observed X-ray LF. The large simulation box allows us accurate determination of the auto-correlation function of the AGNs. The model prediction indicates that this time-scale plays a significant role in allowing massive haloes to host relatively faint population of AGNs, leading to a higher bias factor for those AGNs. The model predictions are in agreement with observations of X-ray selected AGNs in the luminosity range $10^{41.5}~\mathrm{erg} \ \mathrm{s}^{-1} \leq L_{2-10\mathrm{keV}} \leq 10^{44.5}~\mathrm{erg} \ \mathrm{s}^{-1}$, with the typical host halo mass of $10^{12.5-13.5} h^{-1}\,{\rm M}_{\odot}$ at $z \lesssim 1$. This result shows that the observational clustering measurements impose an independent constraint on the accretion time-scale complementary to the LF measurements. Moreover, we find that not only the effective halo mass corresponding to the overall bias factor, but the extended shape of the predicted AGN correlation function shows remarkable agreement with those from observations. Further observational efforts towards the low luminosity end at $z \sim 1$ would give us stronger constraints on the triggering mechanisms of AGN activities through their clustering. △ Less

Submitted 4 July, 2020; v1 submitted 25 December, 2019; originally announced December 2019.

Comments: 20 pages, 13 figures, accepted in MNRAS. Error estimates of the 2PCFs and the fitting to them have been improved. Further discussions and references have been added. Main results remain unchanged

Report number: YITP-19-125

Mass accretion toward black holes in the final phase of galaxy mergers

Authors: Takeru Kawaguchi, Naomichi Yutani, Keiichi Wada

Abstract: We studied the final phases of galactic mergers, focusing on interactions between supermassive black holes (SMBHs) and the interstellar medium in a central sub-kpc region, using an N-body/hydrodynamics code. This numerical experiment aims to understand the fate of the gas supplied by mergers of two or more galaxies with SMBHs, whose masses are $10^7 M_\odot$. We observed that the mass accretion ra… ▽ More We studied the final phases of galactic mergers, focusing on interactions between supermassive black holes (SMBHs) and the interstellar medium in a central sub-kpc region, using an N-body/hydrodynamics code. This numerical experiment aims to understand the fate of the gas supplied by mergers of two or more galaxies with SMBHs, whose masses are $10^7 M_\odot$. We observed that the mass accretion rate to one SMBH exceeds the Eddington accretion rate when the distance between two black holes (BHs) rapidly decreases. However, this rapid accretion phase does not last for more than $10^7$ yrs, and it drops to $\sim$ 10% of the Eddington rate in the quasi-steady accretion phase. The rapid accretion is caused by the angular momentum transfer from the gas to the stellar component, and the moderate accretion in the quasi-steady phase is caused {by turbulent viscosity and gravitational torque in the disk. The second merger event enhances the mass accretion to the BHs; however, this phase takes place on a similar timescale to the first merger event. We also found that the AGN feedback and the mass accretion to BHs can coexist in the central region of merged galaxies, if the amount of feedback energy is given as $(2 \times 10^{-4} - 2 \times 10^{-3} )\dot{M} c^2$, where $\dot{M}$ is the accretion rate to $r= 1$ pc. The accretion rate is suppressed by $\sim$ 1/50 in the quasi-steady accretion phase for $0.02 \dot{M} c^2$. The fraction of the gas that finally falls to each BH is approximately 5-7% of the supplied total gas mass ($10^8 M_\odot$), and 15-20% of the gas forms a circumnuclear gas inside 100 pc. This remnant gas heavily obscures the luminous phase of the active galactic nuclei (AGN) during merger events, and the moderate AGN feedback does not alter this property. △ Less

Submitted 23 December, 2019; originally announced December 2019.

Comments: 22 pages, 16 figures, accepted by ApJ

Tracing the co-evolution path of super massive black holes and spheroids with AKARI-selected ultra-luminous IR galaxies at intermediate redshifts

Authors: Xiao-Yang Chen, Masayuki Akiyama, Kohei Ichikawa, Hirofumi Noda, Yoshiki Toba, Issei Yamamura, Toshihiro Kawaguchi, Abdurro'uf, Mitsuru Kokubo

Abstract: We present the stellar population and ionized-gas outflow properties of ultra-luminous IR galaxies (ULIRGs) at $z=$ 0.1-1.0, which are selected from AKARI FIR all-sky survey. We construct a catalog of 1077 ULIRGs to examine feedback effect after major mergers. 202 out of the 1077 ULIRGs are spectroscopically identified by SDSS and Subaru/FOCAS observations. Thanks to deeper depth and higher resolu… ▽ More We present the stellar population and ionized-gas outflow properties of ultra-luminous IR galaxies (ULIRGs) at $z=$ 0.1-1.0, which are selected from AKARI FIR all-sky survey. We construct a catalog of 1077 ULIRGs to examine feedback effect after major mergers. 202 out of the 1077 ULIRGs are spectroscopically identified by SDSS and Subaru/FOCAS observations. Thanks to deeper depth and higher resolution of AKARI compared to the previous IRAS survey, and reliable identification from WISE MIR pointing, the sample is unique in identifying optically-faint (i$\sim$20) IR-bright galaxies, which could be missed in previous surveys. A self-consistent spectrum-SED decomposition method, which constrains stellar population properties in SED modeling based on spectral fitting results, has been employed for 149 ULIRGs whose optical continua are dominated by host galaxies. They are massive galaxies ($M_{\rm star}\sim10^{11}$-$10^{12}$ M$_{\odot}$), associated with intense star formation activities (SFR $\sim$ 200-2000 M$_{\odot}$ yr$^{-1}$). The sample covers a range of AGN bolometric luminosity of $10^{10}$-$10^{13}$ L$_{\odot}$, and the outflow velocity measured from [OIII] 5007A line shows a correlation with AGN luminosity. Eight galaxies show extremely fast outflows with velocity up to 1500-2000 km s$^{-1}$. However, the co-existence of vigorous starbursts and strong outflows suggests the star formation has not been quenched during the ULIRG phase. By deriving stellar mass and mass fraction of young stellar population, we find no significant discrepancies between stellar properties of ULIRGs with weak and powerful AGNs. The results are not consistent with the merger-induced evolutionary scenario, which predicts that SF-dominated ULIRGs show smaller stellar mass and younger stellar populations compared to AGN-dominated ULIRGs. △ Less

Submitted 13 July, 2020; v1 submitted 11 November, 2019; originally announced November 2019.

Comments: 28 pages, 25 figures, 7 tables; accepted by ApJ

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

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

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

Submitted 19 August, 2019; originally announced August 2019.

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

Journal ref: ApJ 883:183, 2019

Slowing down of cosmic growth of supermassive black holes: Theoretical prediction of the Eddington ratio distribution

Authors: Hikari Shirakata, Toshihiro Kawaguchi, Taira Oogi, Takashi Okamoto, Masahiro Nagashima

Abstract: We show the Eddington ratio distributions of supermassive black holes at a wide redshift range (0 < z < 8) obtained with a semi-analytic model of galaxy formation. The distribution is broadly consistent with observational estimates at low redshift. We find that the growth rate of black holes at higher redshift is more likely to exceed the Eddington limit because the typical gas fraction of the hos… ▽ More We show the Eddington ratio distributions of supermassive black holes at a wide redshift range (0 < z < 8) obtained with a semi-analytic model of galaxy formation. The distribution is broadly consistent with observational estimates at low redshift. We find that the growth rate of black holes at higher redshift is more likely to exceed the Eddington limit because the typical gas fraction of the host galaxies is higher at higher redshift. We also find that the super- Eddington growth is more common for less massive supermassive black holes, supporting an idea that supermassive black holes have been formed via super-Eddington accretion. These results indicate the "slowing down" of cosmic growth of supermassive black holes: the growth of supermassive black holes with a higher Eddington ratio peaks at higher redshift. We also show the effect of the sample selection on the shape of the Eddington ratio distribution functions and find that shallower observations will miss active galactic nuclei with not only the smaller but also higher Eddington ratios. △ Less

Submitted 7 May, 2019; originally announced May 2019.

Comments: 11 pages, 13 figures. Accepted for publication in MNRAS

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

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

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

Submitted 3 May, 2019; originally announced May 2019.

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

X-Ray Spectral Model from Clumpy Torus and Its Application to Circinus Galaxy

Authors: Atsushi Tanimoto, Yoshihiro Ueda, Hirokazu Odaka, Toshihiro Kawaguchi, Yasushi Fukazawa, Taiki Kawamuro

Abstract: We construct an X-ray spectral model from the clumpy torus in an active galactic nucleus (AGN), designated as "XCLUMPY", utilizing the Monte Carlo simulation for Astrophysics and Cosmology framework (MONACO: Odaka et al. 2011, 2016). The adopted geometry of the torus is the same as that in Nenkova et al. (2008), who assume a power law distribution of clumps in the radial direction and a normal dis… ▽ More We construct an X-ray spectral model from the clumpy torus in an active galactic nucleus (AGN), designated as "XCLUMPY", utilizing the Monte Carlo simulation for Astrophysics and Cosmology framework (MONACO: Odaka et al. 2011, 2016). The adopted geometry of the torus is the same as that in Nenkova et al. (2008), who assume a power law distribution of clumps in the radial direction and a normal distribution in the elevation direction. We investigate the dependence of the X-ray continuum and Fe K$α$ fluorescence line profile on the torus parameters. Our model is compared with other torus models: MYTorus model (Murphy & Yaqoob 2009), Ikeda model (Ikeda et al. 2009), and CTorus model (Liu & Li 2014). As an example, we also present the results applied to the broadband X-ray spectra of the Circinus galaxy observed with XMM-Newton, Suzaku, and NuSTAR. Our model can well reproduce the data, yielding a hydrogen column density along the equatorial plane $N_{\mathrm{H}}^{\mathrm{Equ}} = 9.08_{-0.08}^{+0.14} \times 10^{24}$ cm$^{-2}$, a torus angular width $σ= 14.7_{-0.39}^{+0.44}$ degree, and a 2--10 keV luminosity $\log L_{2-10}/\mathrm{erg \ s^{-1}} = 42.8$. These results are discussed in comparison with the observations in other wavelengths. △ Less

Submitted 18 April, 2019; originally announced April 2019.

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

Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs). VI. Black Hole Mass Measurements of Six Quasars at 6.1<z<6.7

Authors: Masafusa Onoue, Nobunari Kashikawa, Yoshiki Matsuoka, Nanako Kato, Takuma Izumi, Tohru Nagao, Michael A. Strauss, Yuichi Harikane, Masatoshi Imanishi, Kei Ito, Kazushi Iwasawa, Toshihiro Kawaguchi, Chien-Hsiu Lee, Akatoki Noboriguchi, Hyewon Suh, Masayuki Tanaka, Yoshiki Toba

Abstract: We present deep near-infrared spectroscopy of six quasars at 6.1<z<6.7 with VLT/X-Shooter and Gemini-N/GNIRS. Our objects, originally discovered through a wide-field optical survey with the Hyper Suprime-Cam (HSC) Subaru Strategic Program (HSC-SSP), have the lowest luminosities (-25.5< M1450<-23.1 mag) of the z>5.8 quasars with measured black hole masses. From single-epoch mass measurements based… ▽ More We present deep near-infrared spectroscopy of six quasars at 6.1<z<6.7 with VLT/X-Shooter and Gemini-N/GNIRS. Our objects, originally discovered through a wide-field optical survey with the Hyper Suprime-Cam (HSC) Subaru Strategic Program (HSC-SSP), have the lowest luminosities (-25.5< M1450<-23.1 mag) of the z>5.8 quasars with measured black hole masses. From single-epoch mass measurements based on MgII2798, we find a wide range in black hole masses, from M_BH=10^7.6 to 10^9.3 Msun. The Eddington ratios L_bol/L_Edd range from 0.16 to 1.1, but the majority of the HSC quasars are powered by M_BH=10^9 Msun supermassive black holes (SMBHs) accreting at sub-Eddington rates. The Eddington ratio distribution of the HSC quasars is inclined to lower accretion rates than those of Willott et al. (2010a), who measured the black hole masses for similarly faint z=6 quasars. This suggests that the global Eddington ratio distribution is wider than has previously been thought. The presence of M_BH=10^9 Msun SMBHs at z=6 cannot be explained with constant sub-Eddington accretion from stellar remnant seed black holes. Therefore, we may be witnessing the first buildup of the most massive black holes in the first billion years of the universe, the accretion activity of which is transforming from active growth to a quiescent phase. Measurements of a larger complete sample of z>6 low-luminosity quasars, as well as deeper observations with future facilities will enable us to better understand the early SMBH growth in the reionization epoch. △ Less

Submitted 12 July, 2019; v1 submitted 15 April, 2019; originally announced April 2019.

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

Application of Clumpy Torus Model to Broadband X-ray Spectra of Two Seyfert 1 Galaxies: IC 4329A and NGC 7469

Authors: Shoji Ogawa, Yoshihiro Ueda, Satoshi Yamada, Atsushi Tanimoto, Toshihiro Kawaguchi

Abstract: We apply a new X-ray clumpy torus model called XCLUMPY (Tanimoto et al. 2019), where the clump distribution is assumed to be the same as in the infrared clumpy torus model (CLUMPY) by Nenkova et al. 2008, to the broadband X-ray spectra of type-1 active galactic nuclei (AGNs) for the first time. We analyze the archival data of IC 4329A and NGC 7469 observed with NuSTAR/Suzaku and NuSTAR/XMM-Newton,… ▽ More We apply a new X-ray clumpy torus model called XCLUMPY (Tanimoto et al. 2019), where the clump distribution is assumed to be the same as in the infrared clumpy torus model (CLUMPY) by Nenkova et al. 2008, to the broadband X-ray spectra of type-1 active galactic nuclei (AGNs) for the first time. We analyze the archival data of IC 4329A and NGC 7469 observed with NuSTAR/Suzaku and NuSTAR/XMM-Newton, respectively, whose infrared spectra were studied with CLUMPY by Ichikawa et al. 2015 and optical extinctions ($A_{\rm V}$) of the tori were estimated. We consider two models, invoking (Model 1) a relativistic reflection component from the accretion disk and (Model 2) a partial absorber. Assuming that the narrow Fe K$α$ emission line at 6.4 keV originates from the torus, we separate the contribution of the torus reflection components in the total spectra. Our models yield equatorial hydrogen column densities of the tori to be $N^{\rm Equ}_{\rm H} =$ (0.53--1.43) $\times 10^{23}~\rm cm^{-2}$ and $N^{\rm Equ}_{\rm H} =$ (0.84--1.43) $\times 10^{24}~\rm cm^{-2}$, for IC 4329A and NGC 7469, respectively. We find that the $N_{\rm H}/A_{\rm V}$ ratios in the tori are by factors of 25--68 (IC 4329A) and 2.4--3.9 (NGC 7469) smaller than that in the Galactic interstellar medium (ISM). These results suggest that a non-negligible fraction of AGNs are 'dust-rich' compared with the Galactic ISM, as opposite to the general trend previously reported in many obscured AGNs. △ Less

Submitted 27 March, 2019; originally announced March 2019.

Comments: 13 pages, 11 figures, (accepted in the Astrophysical Journal, ApJ)

SILVERRUSH. VIII. Spectroscopic Identifications of Early Large Scale Structures with Protoclusters Over 200 Mpc at z~6-7: Strong Associations of Dusty Star-Forming Galaxies

Authors: Yuichi Harikane, Masami Ouchi, Yoshiaki Ono, Seiji Fujimoto, Darko Donevski, Takatoshi Shibuya, Andreas L. Faisst, Tomotsugu Goto, Bunyo Hatsukade, Nobunari Kashikawa, Kotaro Kohno, Takuya Hashimoto, Ryo Higuchi, Akio K. Inoue, Yen-Ting Lin, Crystal L. Martin, Roderik Overzier, Ian Smail, Jun Toshikawa, Hideki Umehata, Yiping Ao, Scott Chapman, David L. Clements, Myungshin Im, Yipeng Jing , et al. (10 additional authors not shown)

Abstract: We have obtained three-dimensional maps of the universe in $\sim200\times200\times80$ comoving Mpc$^3$ (cMpc$^3$) volumes each at $z=5.7$ and $6.6$ based on a spectroscopic sample of 179 galaxies that achieves $\gtrsim80$\% completeness down to the Ly$α$ luminosity of $\log(L_{\rm Lyα}/[\mathrm{erg\ s^{-1}}])=43.0$, based on our Keck and Gemini observations and the literature. The maps reveal fila… ▽ More We have obtained three-dimensional maps of the universe in $\sim200\times200\times80$ comoving Mpc$^3$ (cMpc$^3$) volumes each at $z=5.7$ and $6.6$ based on a spectroscopic sample of 179 galaxies that achieves $\gtrsim80$\% completeness down to the Ly$α$ luminosity of $\log(L_{\rm Lyα}/[\mathrm{erg\ s^{-1}}])=43.0$, based on our Keck and Gemini observations and the literature. The maps reveal filamentary large-scale structures and two remarkable overdensities made out of at least 44 and 12 galaxies at $z=5.692$ (z57OD) and $z=6.585$ (z66OD), respectively, making z66OD the most distant overdensity spectroscopically confirmed to date with $>10$ spectroscopically confirmed galaxies. We compare spatial distributions of submillimeter galaxies at $z\simeq 4-6$ with our $z=5.7$ galaxies forming the large-scale structures, and detect a $99.97\%$ signal of cross correlation, indicative of a clear coincidence of dusty star-forming galaxy and dust unobscured galaxy formation at this early epoch. The galaxies in z57OD and z66OD are actively forming stars with star formation rates (SFRs) $\gtrsim5$ times higher than the main sequence, and particularly the SFR density in z57OD is 10 times higher than the cosmic average at the redshift (a.k.a. the Madau-Lilly plot). Comparisons with numerical simulations suggest that z57OD and z66OD are protoclusters that are progenitors of the present-day clusters with halo masses of $\sim10^{14}\ \mathrm{M_\odot}$. △ Less

Submitted 24 June, 2019; v1 submitted 25 February, 2019; originally announced February 2019.

Comments: 17 pages, 11 figures, Accepted for Publication in ApJ

Discovery of the First Low-Luminosity Quasar at z > 7

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

Abstract: We report the discovery of a quasar at z = 7.07, which was selected from the deep multi-band imaging data collected by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. This quasar, HSC J124353.93+010038.5, has an order of magnitude lower luminosity than do the other known quasars at z > 7. The rest-frame ultraviolet absolute magnitude is M1450 = -24.13 +/- 0.08 mag and the bolometric l… ▽ More We report the discovery of a quasar at z = 7.07, which was selected from the deep multi-band imaging data collected by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. This quasar, HSC J124353.93+010038.5, has an order of magnitude lower luminosity than do the other known quasars at z > 7. The rest-frame ultraviolet absolute magnitude is M1450 = -24.13 +/- 0.08 mag and the bolometric luminosity is Lbol = (1.4 +/- 0.1) x 10^{46} erg/s. Its spectrum in the optical to near-infrared shows strong emission lines, and shows evidence for a fast gas outflow, as the C IV line is blueshifted and there is indication of broad absorption lines. The Mg II-based black hole mass is Mbh = (3.3 +/- 2.0) x 10^8 Msun, thus indicating a moderate mass accretion rate with an Eddington ratio 0.34 +/- 0.20. It is the first z > 7 quasar with sub-Eddington accretion, besides being the third most distant quasar, known to date. The luminosity and black hole mass are comparable to, or even lower than, those measured for the majority of low-z quasars discovered by the Sloan Digital Sky Survey, and thus this quasar likely represents a z > 7 counterpart to quasars commonly observed in the low-z universe. △ Less

Submitted 29 January, 2019; originally announced January 2019.

Comments: Accepted for publication in ApJ Letters

Journal ref: ApJL 872:L2, 2019

Discovery of a strong ionized-gas outflow in an AKARI-selected Ultra-luminous Infrared Galaxy at z = 0.5

Authors: Xiaoyang Chen, Masayuki Akiyama, Hirofumi Noda, Abdurrouf, Yoshiki Toba, Issei Yamamura, Toshihiro Kawaguchi, Mitsuru Kokubo, Kohei Ichikawa

Abstract: In order to construct a sample of ultra-luminous infrared galaxies (ULIRGs, with infrared luminosity, $L_{\rm IR} > 10^{12}$ L$_{\odot}$) at 0.5 < z < 1, we are conducting an optical follow-up program for bright 90-$μ$m FIR sources with a faint optical (i < 20) counterpart selected in the AKARI Far-Infrared Surveyor (FIS) Bright Source catalog (Ver.2). AKARI-FIS-V2 J0916248+073034, identified as a… ▽ More In order to construct a sample of ultra-luminous infrared galaxies (ULIRGs, with infrared luminosity, $L_{\rm IR} > 10^{12}$ L$_{\odot}$) at 0.5 < z < 1, we are conducting an optical follow-up program for bright 90-$μ$m FIR sources with a faint optical (i < 20) counterpart selected in the AKARI Far-Infrared Surveyor (FIS) Bright Source catalog (Ver.2). AKARI-FIS-V2 J0916248+073034, identified as a ULIRG at z = 0.49 in the spectroscopic follow-up observation, indicates signatures of an extremely strong outflow in its emission line profiles. Its [OIII] 5007Å emission line shows FWHM of 1830 km s$^{-1}$ and velocity shift of -770 km s$^{-1}$ in relative to the stellar absorption lines. Furthermore, low-ionization [OII] 3726Å 3729Å doublet also shows large FWHM of 910 km s$^{-1}$ and velocity shift of -380 km s$^{-1}$. After the removal of an unresolved nuclear component, the long-slit spectroscopy 2D image possibly shows that the outflow extends to radius of 4 kpc. The mass outflow and energy ejection rates are estimated to be 500 M$_{\odot}$ yr$^{-1}$ and $4\times10^{44}$ erg s$^{-1}$, respectively, which imply that the outflow is among the most powerful ones observed in ULIRGs and QSOs at 0.3 < z < 1.6. The co-existence of the strong outflow and intense star formation (star formation rate of 990 M$_{\odot}$ yr$^{-1}$) indicates that the feedback of the strong outflow has not severely affect the star-forming region of the galaxy. △ Less

Submitted 16 January, 2019; originally announced January 2019.

Comments: 20 pages, 18 figures, accepted to be published in PASJ

Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs). V. Quasar Luminosity Function and Contribution to Cosmic Reionization at z = 6

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

Abstract: We present new measurements of the quasar luminosity function (LF) at $z \sim 6$, over an unprecedentedly wide range of the rest-frame ultraviolet luminosity $M_{1450}$ from $-30$ to $-22$ mag. This is the fifth in a series of publications from the Subaru High-$z$ Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-… ▽ More We present new measurements of the quasar luminosity function (LF) at $z \sim 6$, over an unprecedentedly wide range of the rest-frame ultraviolet luminosity $M_{1450}$ from $-30$ to $-22$ mag. This is the fifth in a series of publications from the Subaru High-$z$ Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The LF was calculated with a complete sample of 110 quasars at $5.7 \le z \le 6.5$, which includes 48 SHELLQs quasars discovered over 650 deg$^2$, and 63 brighter quasars discovered by the Sloan Digital Sky Survey and the Canada-France-Hawaii Quasar Survey (including one overlapping object). This is the largest sample of $z \sim 6$ quasars with a well-defined selection function constructed to date, and has allowed us to detect significant flattening of the LF at its faint end. A double power-law function fit to the sample yields a faint-end slope $α= -1.23^{+0.44}_{-0.34}$, a bright-end slope $β= -2.73^{+0.23}_{-0.31}$, a break magnitude $M_{1450}^* = -24.90^{+0.75}_{-0.90}$, and a characteristic space density $Φ^* = 10.9^{+10.0}_{-6.8}$ Gpc$^{-3}$ mag$^{-1}$. Integrating this best-fit model over the range $-18 < M_{1450} < -30$ mag, quasars emit ionizing photons at the rate of $\dot{n}_{\rm ion} = 10^{48.8 \pm 0.1}$ s$^{-1}$ Mpc$^{-3}$ at $z = 6.0$. This is less than 10 % of the critical rate necessary to keep the intergalactic medium ionized, which indicates that quasars are not a major contributor to cosmic reionization. △ Less

Submitted 5 November, 2018; originally announced November 2018.

Comments: Accepted for publication in ApJ

Journal ref: ApJ 869:150, 2018

Compton Shoulder Diagnostics in Active Galactic Nuclei for Probing the Metalicity of the Obscuring Compton-Thick Tori

Authors: Masaya Hikitani, Masanori Ohno, Yasushi Fukazawa, Toshihiro Kawaguchi, Hirokazu Odaka

Abstract: We analyzed the spectral shape of the Compton shoulder around the neutral Fe-K$_α$ line of the Compton-thick type II Seyfert nucleus of the Circinus galaxy. The characteristics of this Compton shoulder with respect to the reflected continuum and Fe-K$_α$ line core intensity are a powerful diagnostics tool for analyzing the structure of the molecular tori, which obscure the central engine. We appli… ▽ More We analyzed the spectral shape of the Compton shoulder around the neutral Fe-K$_α$ line of the Compton-thick type II Seyfert nucleus of the Circinus galaxy. The characteristics of this Compton shoulder with respect to the reflected continuum and Fe-K$_α$ line core intensity are a powerful diagnostics tool for analyzing the structure of the molecular tori, which obscure the central engine. We applied our Monte-Carlo-based X-ray reflection spectral model to the Chandra High Energy Transmission Grating data and successfully constrained the various spectral parameters independently, using only the spectral data only around the Fe-K$_α$ emission line. The obtained column density and inclination angle are consistent with the previous observations and the Compton-thick type II Seyfert picture. In addition, we determined the metal abundance of the molecular torus for the case of the smooth and clumpy torus to be 1.75$^{+0.19}_{-0.17}$ and 1.74$\pm$0.16 solar abundance, respectively. Such slightly over-solar abundance can be useful information for discussing the star formation rate in the molecular tori of active galactic nuclei. △ Less

Submitted 2 October, 2018; v1 submitted 2 October, 2018; originally announced October 2018.

Comments: 8 pages, 3 figures, 1 table, accepted for publication in The Astrophysical Journal

A 100-pc Scale, Fast and Dense Outflow in Narrow-Line Seyfert 1 Galaxy IRAS04576+0912

Authors: Toshihiro Kawaguchi, Shinobu Ozaki, Hajime Sugai, Kazuya Matsubayashi, Takashi Hattori, Atsushi Shimono, Kentaro Aoki, Yutaka Hayano, Yosuke Minowa, Kazuma Mitsuda, Yasuhito Hashiba

Abstract: We report the initial result of an adaptive-optics assisted, optical integral-field-unit observation on IRAS04576+0912, the nearest (z=0.039) active galactic nucleus with a prominent blueshift/tail in [O III] emission from a sample of such objects that we have collected from the literature. We aim at addressing the putative quasar-mode feedback process with Subaru/Kyoto 3D II+AO188. The optical wa… ▽ More We report the initial result of an adaptive-optics assisted, optical integral-field-unit observation on IRAS04576+0912, the nearest (z=0.039) active galactic nucleus with a prominent blueshift/tail in [O III] emission from a sample of such objects that we have collected from the literature. We aim at addressing the putative quasar-mode feedback process with Subaru/Kyoto 3D II+AO188. The optical waveband (6400--7500 AA) enables us to measure the gas density via the [S II] doublets, in contrast to earlier Near-IR studies. Since the fast [O III] outflow happens only around rapidly growing central black holes, this object is suitable for investigating the black hole-galaxy coevolution. The obtained data cube exhibits blue tail in the [S II] emission at many lenslets. By fitting the spectrum with the high excess flux at the [S II] blue tail, we find the fast (~ 860 km/s), dense (>3000/cc), wide-angle and offset outflow in central 100-pc scales. Although the large opening angle and the high gas outflow-to-accretion ratio may favour the feedback hypothesis, the inferred kinetic power injection rate of this ionized gas outflow seems insufficient to influence the whole host galaxy. A conventional assumption of a low density must have overestimated the feedback process. △ Less

Submitted 19 July, 2018; v1 submitted 17 July, 2018; originally announced July 2018.

Comments: 7 pages, 2 figures, Accepted for publication in PASJ