Showing 1–50 of 203 results for author: Imanishi, M

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

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

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

Submitted 1 July, 2024; originally announced July 2024.

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

GATOS: missing molecular gas in the outflow of NGC5728 revealed by JWST

Authors: R. Davies, T. Shimizu, M. Pereira-Santaella, A. Alonso-Herrero, A. Audibert, E. Bellocchi, P. Boorman, S. Campbell, Y. Cao, F. Combes, D. Delaney, T. Diaz-Santos, F. Eisenhauer, D. Esparza Arredondo, H. Feuchtgruber, N. M. Forster Schreiber, L. Fuller, P. Gandhi, I. Garcia-Bernete, S. Garcia-Burillo, B. Garcia-Lorenzo, R. Genzel, S. Gillessen, O. Gonzalez Martin, H. Haidar , et al. (27 additional authors not shown)

Abstract: The ionisation cones of NGC5728 have a deficit of molecular gas based on millimetre observations of CO(2-1) emission. Although photoionisation from the active nucleus may lead to suppression of this transition, warm molecular gas can still be present. We report the detection of eight mid-infrared rotational H$_2$ lines throughout the central kiloparsec, including the ionisation cones, using integr… ▽ More The ionisation cones of NGC5728 have a deficit of molecular gas based on millimetre observations of CO(2-1) emission. Although photoionisation from the active nucleus may lead to suppression of this transition, warm molecular gas can still be present. We report the detection of eight mid-infrared rotational H$_2$ lines throughout the central kiloparsec, including the ionisation cones, using integral field spectroscopic observations with JWST/MIRI MRS. The H$_2$ line ratios, characteristic of a power-law temperature distribution, indicate that the gas is warmest where it enters the ionisation cone through disk rotation, suggestive of shock excitation. In the nucleus, where the data can be combined with an additional seven ro-vibrational H$_2$ transitions, we find that moderate velocity (30 km s$^{-1}$) shocks in dense ($10^5$ cm$^{-3}$) gas, irradiated by an external UV field ($G_0 = 10^3$), do provide a good match to the full set. The warm molecular gas in the ionisation cone that is traced by the H$_2$ rotational lines has been heated to temperatures $>200$ K. Outside of the ionisation cone the molecular gas kinematics are undisturbed. However, within the ionisation cone, the kinematics are substantially perturbed, indicative of a radial flow, but one that is quantitatively different from the ionised lines. We argue that this outflow is in the plane of the disk, implying a short 50 pc acceleration zone up to speeds of about 400 km s$^{-1}$ followed by an extended deceleration over $\sim$700 pc where it terminates. The deceleration is due to both the radially increasing galaxy mass, and mass-loading as ambient gas in the disk is swept up. △ Less

Submitted 24 June, 2024; originally announced June 2024.

Comments: A&A accepted; 16 pages

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)

A spectacular galactic scale magnetohydrodynamic powered wind in ESO 320-G030

Authors: M. D. Gorski, S. Aalto, S. König, C. F. Wethers, C. Yang, S. Muller, K. Onishi, M. Sato, N. Falstad, Jeffrey G. Mangum, S. T. Linden, F. Combes, S. Martín, M. Imanishi, Keiichi Wada, L. Barcos-Muñoz, F. Stanley, S. García-Burillo, P. P. van der Werf, A. S. Evans, C. Henkel, S. Viti, N. Harada, T. Díaz-Santos, J. S. Gallagher , et al. (1 additional authors not shown)

Abstract: How galaxies regulate nuclear growth through gas accretion by supermassive black holes (SMBHs) is one of the most fundamental questions in galaxy evolution. One potential way to regulate nuclear growth is through a galactic wind that removes gas from the nucleus. It is unclear whether galactic winds are powered by jets, mechanical winds, radiation, or via magnetohydrodynamic (MHD) processes. Compa… ▽ More How galaxies regulate nuclear growth through gas accretion by supermassive black holes (SMBHs) is one of the most fundamental questions in galaxy evolution. One potential way to regulate nuclear growth is through a galactic wind that removes gas from the nucleus. It is unclear whether galactic winds are powered by jets, mechanical winds, radiation, or via magnetohydrodynamic (MHD) processes. Compact obscured nuclei (CONs) represent a significant phase of galactic nuclear growth. These galaxies hide growing SMBHs or unusual starbursts in their very opaque, extremely compact (r $<$ 100 pc) centres. They are found in approximately 30 % of the luminous and ultra-luminous infrared galaxy (LIRG and ULIRG) population. Here, we present high-resolution ALMA observations ($\sim$30 mas, $\sim$5 pc) of ground-state and vibrationally excited HCN towards ESO 320-G030 (IRAS 11506-3851). ESO 320-G030 is an isolated luminous infrared galaxy known to host a compact obscured nucleus and a kiloparsec-scale molecular wind. Our analysis of these high-resolution observations excludes the possibility of a starburst-driven wind, a mechanically or energy driven active galactic nucleus (AGN) wind, and exposes a molecular MDH wind. These results imply that the nuclear evolution of galaxies and the growth of SMBHs are similar to the growth of hot cores or protostars where gravitational collapse of the nuclear torus drives a MHD wind. These results mean galaxies are capable, in part, of regulating the evolution of their nuclei without feedback. △ Less

Submitted 25 March, 2024; originally announced March 2024.

Comments: 9 pages, 10 figures

Journal ref: A&A 684, L11 (2024)

Spectroscopic Constraints on the Mid-Infrared Attenuation Curve: I -- Attenuation Model using PAH Emission

Authors: Thomas S. -Y. Lai, J. D. T. Smith, Els Peeters, Henrik W. W. Spoon, Shunsuke Baba, Masatoshi Imanishi, Takao Nakagawa

Abstract: We introduce a novel model to spectroscopically constrain the mid-infrared (MIR) extinction/attenuation curve from 3--17 um, using Polycyclic Aromatic Hydrocarbon (PAH) emission drawn from an AKARI-Spitzer extragalactic cross-archival dataset. Currently proposed MIR extinction curves vary significantly in their slopes toward the near-infrared, and the variation of the strengths and shapes of the 9… ▽ More We introduce a novel model to spectroscopically constrain the mid-infrared (MIR) extinction/attenuation curve from 3--17 um, using Polycyclic Aromatic Hydrocarbon (PAH) emission drawn from an AKARI-Spitzer extragalactic cross-archival dataset. Currently proposed MIR extinction curves vary significantly in their slopes toward the near-infrared, and the variation of the strengths and shapes of the 9.7 um and 18 um silicate absorption features make MIR spectral modeling and interpretation challenging, particularly for heavily obscured galaxies. By adopting the basic premise that PAH bands have relatively consistent intrinsic ratios within dusty starbursting galaxies, we can, for the first time, empirically determine the overall shape of the MIR attenuation curve by measuring the differential attenuation at specific PAH wavelengths. Our attenuation model shows PAH emission in most (U)LIRGs is unambiguously subjected to attenuation, and we find strong evidence that PAH bands undergo differential attenuation as obscuration increases. Compared to pre-existing results, the MIR attenuation curve derived from the model favors relatively gray continuum absorption from 3-8 $μ$m and silicate features with intermediate strength at 9.7 um but with stronger than typical 18 um opacity. △ Less

Submitted 4 March, 2024; originally announced March 2024.

Comments: 15 pages, 5 figures. Paper re-submitted to ApJ after addressing the first round of comments by the referee

Double, double, toil, and trouble: The tails, bubbles, and knots of the local compact obscured nucleus galaxy NGC4418

Authors: C. F. Wethers, S. Aalto, G. C. Privon, F. Stanley, J. Gallagher, M. Gorski, S. König, K. Onishi, M. Sato, C. Yang, R. Beswick, L. Barcos-Munoz F. Combes, T. Diaz-Santos, A. S. Evans, I. Garcia-Bernete, C. Henkel, M. Imanishi, S. Martín, S. Muller, Y. Nishimura, C. Ricci, D. Rigopoulou, S. Viti

Abstract: Compact obscured nuclei (CONs) are an extremely obscured (N$_{H2}$ >10$^{25}$ cm$^{-2}$) class of galaxy nuclei thought to exist in 20-40 per cent of nearby (ultra-)luminous infrared galaxies. While they have been proposed to represent a key phase of the active galactic nucleus (AGN) feedback cycle, the nature of these CONs - what powers them, their dynamics, and their impact on the host galaxy -… ▽ More Compact obscured nuclei (CONs) are an extremely obscured (N$_{H2}$ >10$^{25}$ cm$^{-2}$) class of galaxy nuclei thought to exist in 20-40 per cent of nearby (ultra-)luminous infrared galaxies. While they have been proposed to represent a key phase of the active galactic nucleus (AGN) feedback cycle, the nature of these CONs - what powers them, their dynamics, and their impact on the host galaxy - remains unknown. This work analyses the large-scale optical properties of the local CON, NGC4418 (z=0.00727). We present new, targeted integral field unit observations of the galaxy with the Multi-Unit Spectroscopic Explorer (MUSE). For the first time, we map the ionised and neutral gas components of the galaxy, along with their dynamical structure, to reveal several previously unknown features of the galaxy. We confirm the presence of a previously postulated blueshifted outflow along the minor axis of NGC4418. We find this outflow to be decelerating and, for the first time, show it to extend bilaterally from the nucleus. We report the discovery of two further outflow structures: a redshifted southern outflow connected to a tail of ionised gas surrounding the galaxy and a blueshifted bubble to the north. In addition to these features, we find the [OIII] emission reveals the presence of knots across the galaxy, which are consistent with regions of the galaxy that have been photoionised by an AGN. Based on the properties of these features, we conclude that the CON in NGC4418 is most likely powered by AGN activity. △ Less

Submitted 27 February, 2024; originally announced February 2024.

Comments: 17 pages, 15 figures, accepted to A&A

CON-quest II. Spatially and spectrally resolved HCN/HCO+ line ratios in local luminous and ultraluminous infrared galaxies

Authors: Y. Nishimura, S. Aalto, M. D. Gorski, S. König, K. Onishi, C. Wethers, C. Yang, L. Barcos-Muñoz, F. Combes, T. Díaz-Santos, J. S. Gallagher, S. García-Burillo, E. González-Alfonso, T. R. Greve, N. Harada, C. Henkel, M. Imanishi, K. Kohno, S. T. Linden, J. G. Mangum, S. Martín, S. Muller, G. C. Privon, C. Ricci, F. Stanley , et al. (2 additional authors not shown)

Abstract: Nuclear regions of ultraluminous and luminous infrared galaxies (U/LIRGs) are powered by starbursts and/or active galactic nuclei (AGNs). These regions are often obscured by extremely high columns of gas and dust. Molecular lines in the submillimeter windows have the potential to determine the physical conditions of these compact obscured nuclei (CONs). We aim to reveal the distributions of HCN an… ▽ More Nuclear regions of ultraluminous and luminous infrared galaxies (U/LIRGs) are powered by starbursts and/or active galactic nuclei (AGNs). These regions are often obscured by extremely high columns of gas and dust. Molecular lines in the submillimeter windows have the potential to determine the physical conditions of these compact obscured nuclei (CONs). We aim to reveal the distributions of HCN and HCO$^+$ emission in local U/LIRGs and investigate whether and how they are related to galaxy properties. Using ALMA, we have conducted sensitive observations of the HCN J=3--2 and HCO$^+$ J=3--2 lines toward 23 U/LIRGs in the local Universe (z < 0.07) with a spatial resolution of ~0.3" (~50--400 pc). We detected both HCN and HCO$^+$ in 21 galaxies, only HCN in one galaxy, and neither in one galaxy. The global HCN/HCO$^+$ line ratios, averaged over scales of ~0.5--4 kpc, range from 0.4 to 2.3, with an unweighted mean of 1.1. These line ratios appear to have no systematic trend with bolometric AGN luminosity or star formation rate. The line ratio varies with position and velocity within each galaxy, with an average interquartile range of 0.38 on a spaxel-by-spaxel basis. In eight out of ten galaxies known to have outflows and/or inflows, we found spatially and kinematically symmetric structures of high line ratios. These structures appear as a collimated bicone in two galaxies and as a thin spherical shell in six galaxies. Non-LTE analysis suggests that the high HCN/HCO$^+$ line ratio in outflows is predominantly influenced by the abundance ratio. Chemical model calculations indicate that the enhancement of HCN abundance in outflows is likely due to high-temperature chemistry triggered by shock heating. These results imply that the HCN/HCO$^+$ line ratio can aid in identifying the outflow geometry when the shock velocity of the outflows is sufficiently high to heat the gas. △ Less

Submitted 25 April, 2024; v1 submitted 23 February, 2024; originally announced February 2024.

Comments: 52 pages, 35 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 686, A48 (2024)

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

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

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

Submitted 16 February, 2024; originally announced February 2024.

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

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

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

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

Submitted 8 February, 2024; originally announced February 2024.

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

VLBI Detection of an Active Radio Source Potentially Driving 100-kpc Scale Emission in the Ultraluminous Infrared Galaxy IRAS F01004$-$2237

Authors: Takayuki J. Hayashi, Yoshiaki Hagiwara, Masatoshi Imanishi

Abstract: The nearby ultraluminous infrared galaxy (ULIRG) IRAS F01004$-$2237 exhibits 100-kpc scale continuum emission at radio wavelengths. The absence of extended X-ray emission in IRAS F01004$-$2237 has suggested an active galactic nucleus (AGN) origin for the extended radio emission, whose properties and role in merging systems still need to be better understood. We present the results of multi-frequen… ▽ More The nearby ultraluminous infrared galaxy (ULIRG) IRAS F01004$-$2237 exhibits 100-kpc scale continuum emission at radio wavelengths. The absence of extended X-ray emission in IRAS F01004$-$2237 has suggested an active galactic nucleus (AGN) origin for the extended radio emission, whose properties and role in merging systems still need to be better understood. We present the results of multi-frequency observations of IRAS F01004$-$2237 conducted by the Very Long Baseline Array at 2.3 and 8.4 GHz. Compact 8.4-GHz continuum emission was detected on a 1-pc scale in the nuclear region with an intrinsic brightness temperature of $10^{8.1}$ K suggesting that the radio source is originated from an AGN, potentially driving the extended emission. In contrast, no significant emission was observed at 2.3 GHz, indicating the presence of low-frequency absorption. This absorption cannot be attributed solely to synchrotron self-absorption; alternatively, free-free absorption due to thermal plasma is mainly at work in the spectrum. From combined perspectives, including mid-infrared and X-ray data, the AGN is obscured in a dense environment. The kinetic power of the nonthermal jet, as inferred from the extended emission, can play a more important role in dispersing the surrounding medium than the thermal outflow in IRAS F01004$-$2237. These findings hint that jet activities in ULIRGs may contribute to AGN feedback during galaxy evolution induced by merger events. △ Less

Submitted 21 May, 2024; v1 submitted 15 January, 2024; originally announced January 2024.

Comments: 14 pages, 4 figures, 3 tables, accepted for publication in ApJ

Deconvolution of JWST/MIRI Images: Applications to an AGN Model and GATOS Observations of NGC 5728

Authors: M. T. Leist, C. Packham, D. J. V. Rosario, D. A. Hope, A. Alonso-Herrero, E. K. S. Hicks, S. Hönig, L. Zhang, R. Davies, T. Díaz-Santos, O. Ganzález-Martín, E. Bellocchi, P. G. Boorman, F. Combes, I. García-Bernete, S. García-Burillo, B. García-Lorenzo, H. Haidar, K. Ichikawa, M. Imanishi, S. M. Jefferies, Á. Labiano, N. A. Levenson, R. Nikutta, M. Pereira-Santaella , et al. (10 additional authors not shown)

Abstract: The superb image quality, stability and sensitivity of the JWST permit deconvolution techniques to be pursued with a fidelity unavailable to ground-based observations. We present an assessment of several deconvolution approaches to improve image quality and mitigate effects of the complex JWST point spread function (PSF). The optimal deconvolution method is determined by using WebbPSF to simulate… ▽ More The superb image quality, stability and sensitivity of the JWST permit deconvolution techniques to be pursued with a fidelity unavailable to ground-based observations. We present an assessment of several deconvolution approaches to improve image quality and mitigate effects of the complex JWST point spread function (PSF). The optimal deconvolution method is determined by using WebbPSF to simulate JWST's complex PSF and MIRISim to simulate multi-band JWST/Mid-Infrared Imager Module (MIRIM) observations of a toy model of an active galactic nucleus (AGN). Five different deconvolution algorithms are tested: (1) Kraken deconvolution, (2) Richardson-Lucy, (3) Adaptive Imaging Deconvolution Algorithm, (4) Sparse regularization with the Condat-Vũ algorithm, and (5) Iterative Wiener Filtering and Thresholding. We find that Kraken affords the greatest FWHM reduction of the nuclear source of our MIRISim observations for the toy AGN model while retaining good photometric integrity across all simulated wavebands. Applying Kraken to Galactic Activity, Torus, and Outflow Survey (GATOS) multi-band JWST/MIRIM observations of the Seyfert 2 galaxy NGC 5728, we find that the algorithm reduces the FWHM of the nuclear source by a factor of 1.6-2.2 across all five filters. Kraken images facilitate detection of a SE to NW $\thicksim$2".5 ($\thicksim$470 pc, PA $\simeq$115°) extended nuclear emission, especially in the longest wavelengths. We demonstrate that Kraken is a powerful tool to enhance faint features otherwise hidden in the complex JWST PSF. △ Less

Submitted 14 February, 2024; v1 submitted 21 December, 2023; originally announced December 2023.

Comments: 32 pages, 23 figures, published in AJ 2024 February 7

SERENADE II: An ALMA Multi-Band Dust-Continuum Analysis of 28 Galaxies at $5<z<8$ and the Physical Origin of the Dust Temperature Evolution

Authors: Ikki Mitsuhashi, Yuichi Harikane, Franz E. Bauer, Tom Bakx, Andrea Ferrara, Seiji Fujimoto, Takuya Hashimoto, Akio K. Inoue, Kazushi Iwasawa, Yuri Nishimura, Masatoshi Imanishi, Yoshiaki Ono, Toshiki Saito, Yuma Sugahara, Hideki Umehata, Livia Vallini, Tao Wang

Abstract: We present an analysis of ALMA multi-band dust-continuum observations for 28 spectroscopically-confirmed bright Lyman-break galaxies at $5<z<8$. Our sample consists of 11 galaxies at $z\sim6$ newly observed in our ALMA program, which substantially increases the number of $5<z<8$ galaxies with both rest-frame 88 and 158 $μ{\rm m}$ continuum observations, allowing us to simultaneously measure the IR… ▽ More We present an analysis of ALMA multi-band dust-continuum observations for 28 spectroscopically-confirmed bright Lyman-break galaxies at $5<z<8$. Our sample consists of 11 galaxies at $z\sim6$ newly observed in our ALMA program, which substantially increases the number of $5<z<8$ galaxies with both rest-frame 88 and 158 $μ{\rm m}$ continuum observations, allowing us to simultaneously measure the IR luminosity and dust temperature for a statistical sample of $z\gtrsim5$ galaxies for the first time. We derive the relationship between the UV slope ($β_{\rm UV}$) and infrared excess (IRX) for the $z\sim6$ galaxies, and find a shallower IRX-$β_{\rm UV}$ relation compared to the previous results at $z\sim2$--4. Based on the IRX-$β_{\rm UV}$ relation consistent with our results and the $β_{\rm UV}$-$M_{\rm UV}$ relation including fainter galaxies in the literature, we find a limited contribution of the dust-obscured star formation to the total SFR density, $\sim30\%$ at $z\sim6$. Our measurements of the dust temperature at $z\sim6-7$, $T_{\rm dust}=40.9_{-9.1}^{+10.0}\,{\rm K}$ on average, supports a gentle increase of $T_{\rm dust}$ from $z=0$ to $z\sim6$--7. Using an analytic model with parameters consistent with recent {\it{JWST}} results, we discuss that the observed redshift evolution of the dust temperature can be reproduced by an $\sim0.6\,{\rm dex}$ increase in the gas depletion timescale and $\sim0.4\,{\rm dex}$ decrease of the metallicity. The variety of $T_{\rm dust}$ observed at high redshifts can also be naturally explained by scatters around the star-formation main sequence and average mass-metallicity relation, including an extremely high dust temperature of $T_{\rm dust}>80\,{\rm K}$ observed in a galaxy at $z=8.3$. △ Less

Submitted 28 November, 2023; originally announced November 2023.

Comments: Submitted to ApJ

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

BASS XXXIV: A Catalog of the Nuclear Mm-wave Continuum Emission Properties of AGNs Constrained on Scales $\lesssim$ 100--200 pc

Authors: Taiki Kawamuro, Claudio Ricci, Richard F. Mushotzky, Masatoshi Imanishi, Franz E. Bauer, Federica Ricci, Michael J. Koss, George C. Privon, Benny Trakhtenbrot, Takuma Izumi, Kohei Ichikawa, Alejandra F. Rojas, Krista Lynne Smith, Taro Shimizu, Kyuseok Oh, Jakob S. den Brok, Shunsuke Baba, Mislav Balokovic, Chin-Shin Chang, Darshan Kakkad, Ryan W. Pfeifle, Matthew J. Temple, Yoshihiro Ueda, Fiona Harrison, Meredith C. Powell , et al. (3 additional authors not shown)

Abstract: We present a catalog of the millimeter-wave (mm-wave) continuum properties of 98 nearby ($z <$ 0.05) active galactic nuclei (AGNs) selected from the 70-month Swift/BAT hard X-ray catalog that have precisely determined X-ray spectral properties and subarcsec-resolution ALMA Band-6 (211--275 GHz) observations as of 2021 April. Due to the hard-X-ray ($>$ 10 keV) selection, the sample is nearly unbias… ▽ More We present a catalog of the millimeter-wave (mm-wave) continuum properties of 98 nearby ($z <$ 0.05) active galactic nuclei (AGNs) selected from the 70-month Swift/BAT hard X-ray catalog that have precisely determined X-ray spectral properties and subarcsec-resolution ALMA Band-6 (211--275 GHz) observations as of 2021 April. Due to the hard-X-ray ($>$ 10 keV) selection, the sample is nearly unbiased for obscured systems at least up to Compton-thick-level obscuration, and provides the largest number of AGNs with high physical resolution mm-wave data ($\lesssim$ 100--200 pc). Our catalog reports emission peak coordinates, spectral indices, and peak fluxes and luminosities at 1.3 mm (230 GHz). Additionally, high-resolution mm-wave images are provided. Using the images and creating radial surface brightness profiles of mm-wave emission, we identify emission extending from the central source and isolated blob-like emission. Flags indicating the presence of these emission features are tabulated. Among 90 AGNs with significant detections of nuclear emission, 37 AGNs ($\approx$ 41%) appear to have both or one of extended or blob-like components. We, in particular, investigate AGNs that show well-resolved mm-wave components and find that these seem to have a variety of origins (i.e., a jet, radio lobes, a secondary AGN, stellar clusters, a narrow line region, galaxy disk, active star-formation regions, and AGN-driven outflows), and some components have currently unclear origins. △ Less

Submitted 6 September, 2023; originally announced September 2023.

Comments: 49 pages, 7 figures, 3 tables, accepted for publication in ApJS

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

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

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

Submitted 4 September, 2023; originally announced September 2023.

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

ALMA 0.5 kpc Resolution Spatially Resolved Investigations of Nuclear Dense Molecular Gas Properties in Nearby Ultraluminous Infrared Galaxies Based on HCN and HCO$^{+}$ Three Transition Line Data

Authors: Masatoshi Imanishi, Shunsuke Baba, Kouichiro Nakanishi, Takuma Izumi

Abstract: We present the results of our ALMA $\lesssim$0.5 kpc-resolution dense molecular line (HCN and HCO$^{+}$ J=2-1, J=3-2, and J=4-3) observations of 12 nearby (ultra)luminous infrared galaxies ([U]LIRGs). After matching beam sizes of all molecular line data to the same values in all (U)LIRGs, we derive molecular line flux ratios, by extracting spectra in the central 0.5, 1, 2 kpc circular regions, and… ▽ More We present the results of our ALMA $\lesssim$0.5 kpc-resolution dense molecular line (HCN and HCO$^{+}$ J=2-1, J=3-2, and J=4-3) observations of 12 nearby (ultra)luminous infrared galaxies ([U]LIRGs). After matching beam sizes of all molecular line data to the same values in all (U)LIRGs, we derive molecular line flux ratios, by extracting spectra in the central 0.5, 1, 2 kpc circular regions, and 0.5-1 and 1-2 kpc annular regions. Based on non-LTE model calculations, we quantitatively confirm that the innermost ($\lesssim$0.5 kpc) molecular gas is very dense ($\gtrsim$10$^{5}$ cm$^{-3}$) and warm ($\gtrsim$300 K) in ULIRGs, and that in one LIRG is also modestly dense (10$^{4-5}$ cm$^{-3}$) and warm ($\sim$100 K). We then investigate the spatial variation of the HCN-to-HCO$^{+}$ flux ratios and high-J to low-J flux ratios of HCN and HCO$^{+}$. A subtle sign of decreasing trend of these ratios from the innermost ($\lesssim$0.5 kpc) to outer nuclear (0.5-2 kpc) region is discernible in a significant fraction of the observed ULIRGs. For two AGN-hosting ULIRGs which display the trend most clearly, we find based on a Bayesian approach that the HCN-to-HCO$^{+}$ abundance ratio and gas kinetic temperature systematically increase from the outer nuclear to the innermost region. We suggest that this trend comes from potential AGN effects, because no such spatial variation is found in a starburst-dominated LIRG. △ Less

Submitted 27 July, 2023; originally announced July 2023.

Comments: 62 pages, 29 figures. ApJ accepted. Higher resolution version is available at http://www2.nao.ac.jp/~masaimanishi/Aph/ImanishiULIRG500pcH.pdf

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

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

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

Submitted 5 July, 2023; originally announced July 2023.

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

Molecular Abundance of the Circumnuclear Region Surrounding an Active Galactic Nucleus in NGC 1068 based on Imaging Line Survey in the 3-mm Band with ALMA

Authors: Taku Nakajima, Shuro Takano, Tomoka Tosaki, Akio Taniguchi, Nanase Harada, Toshiki Saito, Masatoshi Imanishi, Yuri Nishimura, Takuma Izumi, Yoichi Tamura, Kotaro Kohno, Eric Herbst

Abstract: We present an imaging molecular line survey in the 3-mm band (85-114 GHz) focused on one of the nearest galaxies with an active galactic nucleus (AGN), NGC 1068, based on observations taken with the Atacama Large Millimeter/submillimeter Array (ALMA). Distributions of 23 molecular transitions are obtained in the central ~3 kpc region, including both the circumnuclear disk (CND) and starburst ring… ▽ More We present an imaging molecular line survey in the 3-mm band (85-114 GHz) focused on one of the nearest galaxies with an active galactic nucleus (AGN), NGC 1068, based on observations taken with the Atacama Large Millimeter/submillimeter Array (ALMA). Distributions of 23 molecular transitions are obtained in the central ~3 kpc region, including both the circumnuclear disk (CND) and starburst ring (SBR) with 60 and 350 pc resolution. The column densities and relative abundances of all the detected molecules are estimated under the assumption of local thermodynamic equilibrium in the CND and SBR. Then, we discuss the physical and chemical effects of the AGN on molecular abundance corresponding to the observation scale. We found that H13CN, SiO, HCN, and H13CO+ are abundant in the CND relative to the SBR. In contrast, 13CO is more abundant in the SBR. Based on the calculated column density ratios of N(HCN)/N(HCO+), N(HCN)/N(CN), and other molecular distributions, we conclude that the enhancement of HCN in the CND may be due to high-temperature environments resulting from strong shocks, which are traced by the SiO emission. Moreover, the abundance of CN in the CND is significantly lower than the expected value of the model calculations in the region affected by strong radiation. The expected strong X-ray irradiation from the AGN has a relatively lower impact on the molecular abundance in the CND than mechanical feedback. △ Less

Submitted 5 July, 2023; originally announced July 2023.

Comments: 33 pages, 20 figures, 4 tables, accepted for publication in ApJ

Off-nuclear H2O maser and dense molecular gas in NGC1068

Authors: Yoshiaki Hagiwara, Willem A. Baan, Masatoshi Imanishi, Philip Diamond

Abstract: The results of high-resolution spectral-line observations of dense molecular gas are presented towards the nuclear region of the type 2 Seyfert galaxy NGC1068. MERLIN observations of the 22 GHz H2O maser were made for imaging the known off-nuclear maser emission at radio jet component located about 0.3" north-east of the radio nucleus in the galaxy. High angular resolution ALMA observations have s… ▽ More The results of high-resolution spectral-line observations of dense molecular gas are presented towards the nuclear region of the type 2 Seyfert galaxy NGC1068. MERLIN observations of the 22 GHz H2O maser were made for imaging the known off-nuclear maser emission at radio jet component located about 0.3" north-east of the radio nucleus in the galaxy. High angular resolution ALMA observations have spatially resolved the molecular gas emissions of HCN and HCO$^{+}$ in this region. The off-nuclear maser spots are found to nearly overlap with a ring-like molecular gas structure and are tracing an evolving shock-like structure, which appears to be energized by interaction between the radio jet and circumnuclear medium. A dynamic jet-ISM interaction is further supported by a systematic shift of the centroid velocities of the off-nuclear maser features over a period of 35 years. The integrated flux ratios of the HCO$^{+}$ line emission features at component C suggest a kinetic temperature T$_{k}$ $\gtrsim$ 300K and an H$_2$ density of $\gtrsim$ 10$^6$ cm$^{-3}$, which are conditions where water masers may be formed. The diagnostics of the masering action in this jet-ISM interaction region is exemplary for galaxies hosting off-nuclear H2O maser emission. △ Less

Submitted 2 January, 2024; v1 submitted 2 June, 2023; originally announced June 2023.

Comments: 12 pages, 10 figures, accepted (2023 December 29) for publication in the Monthly Notices of the Royal Astronomical Society

AGN feedback in action in the molecular gas ring of the Seyfert galaxy NGC7172

Authors: A. Alonso-Herrero, S. Garcia-Burillo, M. Pereira-Santaella, T. Shimizu, F. Combes, E. K. S. Hicks, R. Davies, C. Ramos Almeida, I. Garcia-Bernete, S. F. Hoenig, N. A. Levenson, C. Packham, E. Bellocchi, L. K. Hunt, M. Imanishi, C. Ricci, P. Roche

Abstract: We present new ALMA observations of the CO(3-2) transition and 854micron continuum at 0.06-0.3" resolution, together with new VLT/SINFONI observations of NGC7172. This is a luminous (bolometric luminosity of ~10^44 erg/s) Seyfert galaxy that belongs to the Galaxy Activity, Torus, and Outflow Survey (GATOS). The CO(3-2) observations reveal the presence of a highly inclined cold molecular gas ring w… ▽ More We present new ALMA observations of the CO(3-2) transition and 854micron continuum at 0.06-0.3" resolution, together with new VLT/SINFONI observations of NGC7172. This is a luminous (bolometric luminosity of ~10^44 erg/s) Seyfert galaxy that belongs to the Galaxy Activity, Torus, and Outflow Survey (GATOS). The CO(3-2) observations reveal the presence of a highly inclined cold molecular gas ring with an approximate radius of 3-4"~540-720 pc, which is likely associated with an inner Lindblad resonance of a putative stellar bar. There are noncircular motions in the VLT/SINFONI [SiVI]1.96micron and H2 at 2.12micron, and ALMA CO(3-2) velocity fields. After subtracting the stellar velocity field, we detected [SiVI] blueshifted velocities of a few hundred km/s to the south of the AGN. They trace outflowing ionized gas outside the plane of the galaxy and out to projected distances of ~200 pc. The CO(3-2) position-velocity diagram along the kinematic minor axis displays noncircular motions with observed velocities of up to ~150 km/s. Assuming that these are taking place in the disk of the galaxy, the observed velocity signs imply that the molecular gas ring is not only rotating but also outflowing. We derived an integrated cold molecular gas mass outflow rate of ~40 Msun/yr for the ring. Using the 854micron map, we resolved a 32 pc radius torus with a gas mass of 8x10^5 Msun. These torus properties are similar to other Seyfert galaxies in the GATOS sample. We measured a decreased cold molecular gas concentration in the nuclear-torus region relative to the circumnuclear region when compared to other less luminous Seyfert galaxies. We conclude that the effects of AGN feedback in NGC7172, which are likely caused by the AGN wind and/or the moderate luminosity radio jet, are seen as a large-scale outflowing molecular gas ring and accompanying redistribution of molecular gas in the nuclear regions. △ Less

Submitted 24 May, 2023; originally announced May 2023.

Comments: Accepted for publication to A&A

Quasar Luminosity Function at z = 7

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

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

Submitted 18 May, 2023; originally announced May 2023.

Comments: The Astrophysical Journal Letters, in press

Supermassive black hole feeding and feedback observed on sub-parsec scales

Authors: Takuma Izumi, Keiichi Wada, Masatoshi Imanishi, Kouichiro Nakanishi, Kotaro Kohno, Yuki Kudoh, Taiki Kawamuro, Shunsuke Baba, Naoki Matsumoto, Yutaka Fujita, Konrad R. W. Tristram

Abstract: Active galaxies contain a supermassive black hole at their center, which grows by accreting matter from the surrounding galaxy. The accretion process in the central ~10 parsecs has not been directly resolved in previous observations, due to the small apparent angular sizes involved. We observed the active nucleus of the Circinus Galaxy using sub-millimeter interferometry. A dense inflow of molecul… ▽ More Active galaxies contain a supermassive black hole at their center, which grows by accreting matter from the surrounding galaxy. The accretion process in the central ~10 parsecs has not been directly resolved in previous observations, due to the small apparent angular sizes involved. We observed the active nucleus of the Circinus Galaxy using sub-millimeter interferometry. A dense inflow of molecular gas is evident on sub-parsec scales. We calculate that less than 3% of this inflow is accreted by the black hole, with the rest being ejected by multiphase outflows, providing feedback to the host galaxy. The observations also reveal a dense gas disk surrounding the inflow; the disk is gravitationally unstable which drives the accretion into the central ~1 parsec. △ Less

Submitted 13 November, 2023; v1 submitted 6 May, 2023; originally announced May 2023.

Comments: First release on Nov 3, 2023 in Science. 32 pages in one column = Main (13 pages, 4 figures) + Supplement (19 pages, 9 figures + 2 tables). This is the accepted version after peer review

Journal ref: Science, Volume 382, Year 2023, Issue 6670, Pages 554-559

Dense Molecular Gas Properties of the Central Kpc of Nearby Ultraluminous Infrared Galaxies Constrained by ALMA Three Transition-line Observations

Authors: Masatoshi Imanishi, Shunsuke Baba, Kouichiro Nakanishi, Takuma Izumi

Abstract: We report the results of ALMA 1-2 kpc-resolution, three rotational transition line (J=2-1, J=3-2, and J=4-3) observations of multiple dense molecular gas tracers (HCN, HCO$^{+}$, and HNC) for ten nearby (ultra)luminous infrared galaxies ([U]LIRGs). Following the matching of beam sizes to 1-2 kpc for each (U)LIRG, the high-J to low-J transition-line flux ratios of each molecule and the emission lin… ▽ More We report the results of ALMA 1-2 kpc-resolution, three rotational transition line (J=2-1, J=3-2, and J=4-3) observations of multiple dense molecular gas tracers (HCN, HCO$^{+}$, and HNC) for ten nearby (ultra)luminous infrared galaxies ([U]LIRGs). Following the matching of beam sizes to 1-2 kpc for each (U)LIRG, the high-J to low-J transition-line flux ratios of each molecule and the emission line flux ratios of different molecules at each J transition are derived. We conduct RADEX non-LTE model calculations and find that, under a wide range of gas density and kinetic temperature, the observed HCN-to-HCO$^{+}$ flux ratios in the overall (U)LIRGs are naturally reproduced with enhanced HCN abundance compared to HCO$^{+}$. Thereafter, molecular gas properties are constrained primarily through the use of HCN and HCO$^{+}$ data and the adoption of fiducial values for the HCO$^{+}$ column density and HCN-to-HCO$^{+}$ abundance ratio. We quantitatively confirm the following: (1) Molecular gas at the (U)LIRGs' nuclei is dense ($\gtrsim$10$^{3-4}$ cm$^{-3}$) and warm ($\gtrsim$100 K). (2) Molecular gas density and temperature in nine ULIRGs' nuclei are significantly higher than that of one LIRG's nucleus. (3) Molecular gas in starburst-dominated sources tends to be less dense and cooler than ULIRGs with luminous AGN signatures. For six selected sources, we also apply a Bayesian approach by freeing all parameters and support the above main results. Our ALMA 1-2 kpc resolution, multiple transition-line data of multiple molecules are a very powerful tool for scrutinizing the properties of molecular gas concentrated around luminous energy sources in nearby (U)LIRGs' nuclei. △ Less

Submitted 14 March, 2023; originally announced March 2023.

Comments: 42 pages, 13 figures. ApJ accepted

UV & Ly$α$ halos of Ly$α$ emitters across environments at z=2.84

Authors: Satoshi Kikuta, Yuichi Matsuda, Shigeki Inoue, Charles C. Steidel, Renyue Cen, Zheng Zheng, Hidenobu Yajima, Rieko Momose, Masatoshi Imanishi, Yutaka Komiyama

Abstract: We present UV and Ly$α$ radial surface brightness (SB) profiles of Ly$α$ emitters (LAEs) at $z=2.84$ detected with the Hyper Suprime-Cam (HSC) on the Subaru Telescope. The depth of our data, together with the wide field coverage including a protocluster, enable us to study the dependence of Ly$α$ halos (LAHs) on various galaxy properties, including Mpc-scale environments. UV and Ly$α$ images of 34… ▽ More We present UV and Ly$α$ radial surface brightness (SB) profiles of Ly$α$ emitters (LAEs) at $z=2.84$ detected with the Hyper Suprime-Cam (HSC) on the Subaru Telescope. The depth of our data, together with the wide field coverage including a protocluster, enable us to study the dependence of Ly$α$ halos (LAHs) on various galaxy properties, including Mpc-scale environments. UV and Ly$α$ images of 3490 LAEs are extracted, and stacking the images yields SB sensitivity of $\sim1\times10^{-20}\mathrm{~erg~s^{-1}~cm^{-2}~arcsec^{-2}}$ in Ly$α$, reaching the expected level of optically thick gas illuminated by the UV background at $z\sim3$. Fitting of the two-component exponential function gives the scale-lengths of $1.56\pm0.01$ and $10.4\pm0.3$ pkpc. Dividing the sample according to their photometric properties, we find that while the dependence of halo scale-length on environment outside of the protocluster core is not clear, LAEs in the central regions of protoclusters appear to have very large LAHs which could be caused by combined effects of source overlapping and diffuse Ly$α$ emission from cool intergalactic gas permeating the forming protocluster core irradiated by active members. For the first time, we identify ``UV halos'' around bright LAEs which are probably due to a few lower-mass satellite galaxies. Through comparison with recent numerical simulations, we conclude that, while scattered Ly$α$ photons from the host galaxies are dominant, star formation in satellites evidently contributes to LAHs, and that fluorescent Ly$α$ emission may be boosted within protocluster cores at cosmic noon and/or near bright QSOs. △ Less

Submitted 24 February, 2023; originally announced February 2023.

Comments: 29 pages, 20 Figures, Accepted for publication in ApJ

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)

Resurrection of type IIL supernova 2018ivc: Implications for a binary evolution sequence connecting hydrogen-rich and -poor progenitors

Authors: Keiichi Maeda, Tomonari Michiyama, Poonam Chandra, Stuart Ryder, Hanindyo Kuncarayakti, Daichi Hiramatsu, Masatoshi Imanishi

Abstract: Long-term observations of synchrotron emission from supernovae (SNe), covering more than a year after the explosion, provide a unique opportunity to study the poorly-understood evolution of massive stars in the final millennium of their lives via changes in the mass-loss rate. Here, we present a result of our long-term monitoring of a peculiar type IIL SN 2018ivc, using the Atacama Large Millimete… ▽ More Long-term observations of synchrotron emission from supernovae (SNe), covering more than a year after the explosion, provide a unique opportunity to study the poorly-understood evolution of massive stars in the final millennium of their lives via changes in the mass-loss rate. Here, we present a result of our long-term monitoring of a peculiar type IIL SN 2018ivc, using the Atacama Large Millimeter/submillimeter Array (ALMA). Following the initial decay, it showed unprecedented rebrightening starting at ~ a year after the explosion. This is one of the rare examples showing such rebrightening in the synchrotron emission, and the first case at millimeter wavelengths. We find it to be in the optically-thin regime unlike the optically-thick centimeter emission. As such, we can robustly reconstruct the distribution of the circumstellar matter (CSM) and thus the mass-loss history in the final ~1,000 years. We find that the progenitor of SN 2018ivc had experienced a very high mass-loss rate >~10^{-3} Msun/yr ~1,500 years before the explosion, which was followed by a moderately high mass-loss rate (>~10^{-4} Msun/yr) up until the explosion. From this behavior, we suggest SN 2018ivc represents an extreme version of a binary evolution toward SNe IIb, which bridges the hydrogen-poor SNe (toward SNe Ib/c, without a hydrogen envelope) and hydrogen-rich SNe (SNe IIP, with a massive envelope). △ Less

Submitted 18 January, 2023; originally announced January 2023.

Comments: 13 pages, 4 figures, 1 table. Accepted for publication in ApJL

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

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

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

Submitted 9 January, 2023; originally announced January 2023.

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

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

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

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

Submitted 7 January, 2023; originally announced January 2023.

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

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

ALMA Observations for CO Emission from Luminous Lyman-break Galaxies at $z=6.0293$-$6.2037$

Authors: Yoshiaki Ono, Seiji Fujimoto, Yuichi Harikane, Masami Ouchi, Livia Vallini, Andrea Ferrara, Takatoshi Shibuya, Andrea Pallottini, Akio K. Inoue, Masatoshi Imanishi, Kazuhiro Shimasaku, Takuya Hashimoto, Chien-Hsiu Lee, Yuma Sugahara, Yoichi Tamura, Kotaro Kohno, Malte Schramm

Abstract: We present our new Atacama Large Millimeter/submillimeter Array (ALMA) observations targeting CO(6-5) emission from three luminous Lyman break galaxies (LBGs) at $z_{\rm spec} = 6.0293$-$6.2037$ found in the Subaru/Hyper Suprime-Cam survey, whose [OIII]$88μ$m and [CII]$158μ$m emission have been detected with ALMA. We find a marginal detection of the CO(6-5) line from one of our LBGs, J0235-0532, a… ▽ More We present our new Atacama Large Millimeter/submillimeter Array (ALMA) observations targeting CO(6-5) emission from three luminous Lyman break galaxies (LBGs) at $z_{\rm spec} = 6.0293$-$6.2037$ found in the Subaru/Hyper Suprime-Cam survey, whose [OIII]$88μ$m and [CII]$158μ$m emission have been detected with ALMA. We find a marginal detection of the CO(6-5) line from one of our LBGs, J0235-0532, at the $\simeq 4 σ$ significance level and obtain upper limits for the other two LBGs, J1211-0118 and J0217-0208. Our $z=6$ luminous LBGs are consistent with the previously found correlation between the CO luminosity and the infrared luminosity. The unique ensemble of the multiple far-infrared emission lines and underlying continuum fed to a photodissociation region model reveal that J0235-0532 has a relatively high hydrogen nucleus density that is comparable to those of low-$z$ (U)LIRGs, quasars, and Galactic star-forming regions with high $n_{\rm H}$ values, while the other two LBGs have lower $n_{\rm H}$ consistent with local star-forming galaxies. By carefully taking account of various uncertainties, we obtain total gas mass and gas surface density constraints from their CO luminosity measurements. We find that J0235-0532 locates below the Kennicutt-Schmidt (KS) relation, comparable to the previously CO(2-1) detected $z=5.7$ LBG, HZ10. Combined with previous results for dusty starbursts at similar redshifts, the KS relation at $z=5$-$6$ is on average consistent with the local one. △ Less

Submitted 1 November, 2022; originally announced November 2022.

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

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

BASS XXXII: Studying the Nuclear Mm-wave Continuum Emission of AGNs with ALMA at Scales $\lesssim$ 100-200 pc

Authors: Taiki Kawamuro, Claudio Ricci, Masatoshi Imanishi, Richard F. Mushotzky, Takuma Izumi, Federica Ricci, Franz E. Bauer, Michael J. Koss, Benny Trakhtenbrot, Kohei Ichikawa, Alejandra F. Rojas, Krista Lynne Smith, Taro Shimizu, Kyuseok Oh, Jakob S. den Brok, Shunsuke Baba, Mislav Baloković, Chin-Shin Chang, Darshan Kakkad, Ryan W. Pfeifle, George C. Privon, Matthew J. Temple, Yoshihiro Ueda, Fiona Harrison, Meredith C. Powell , et al. (3 additional authors not shown)

Abstract: To understand the origin of nuclear ($\lesssim$ 100 pc) millimeter-wave (mm-wave) continuum emission in active galactic nuclei (AGNs), we systematically analyzed sub-arcsec resolution Band-6 (211-275 GHz) ALMA data of 98 nearby AGNs ($z <$ 0.05) from the 70-month Swift/BAT catalog. The sample, almost unbiased for obscured systems, provides the largest number of AGNs to date with high mm-wave spati… ▽ More To understand the origin of nuclear ($\lesssim$ 100 pc) millimeter-wave (mm-wave) continuum emission in active galactic nuclei (AGNs), we systematically analyzed sub-arcsec resolution Band-6 (211-275 GHz) ALMA data of 98 nearby AGNs ($z <$ 0.05) from the 70-month Swift/BAT catalog. The sample, almost unbiased for obscured systems, provides the largest number of AGNs to date with high mm-wave spatial resolution sampling ($\sim$ 1-200 pc), and spans broad ranges of 14-150 keV luminosity {$40 < \log[L_{\rm 14-150}/({\rm erg\,s^{-1}})] < 45$}, black hole mass [$5 < \log(M_{\rm BH}/M_\odot) < 10$], and Eddington ratio ($-4 < \log λ_{\rm Edd} < 2$). We find a significant correlation between 1.3 mm (230 GHz) and 14-150 keV luminosities. Its scatter is $\approx$ 0.36 dex, and the mm-wave emission may serve as a good proxy of the AGN luminosity, free of dust extinction up to $N_{\rm H} \sim 10^{26}$ cm$^{-2}$. While the mm-wave emission could be self-absorbed synchrotron radiation around the X-ray corona according to past works, we also discuss different possible origins of the mm-wave emission; AGN-related dust emission, outflow-driven shocks, and a small-scale ($<$ 200 pc) jet. The dust emission is unlikely to be dominant, as the mm-wave slope is generally flatter than expected. Also, due to no increase in the mm-wave luminosity with the Eddington ratio, a radiation-driven outflow model is possibly not the common mechanism. Furthermore, we find independence of the mm-wave luminosity on indicators of the inclination angle from the polar axis of the nuclear structure, which is inconsistent with a jet model whose luminosity depends only on the angle. △ Less

Submitted 10 August, 2022; v1 submitted 7 August, 2022; originally announced August 2022.

Comments: 44 pages, 35 figures, 3 tables, accepted for publication in ApJ; modified format, added supplementary figure (Fig. 32)

From the Circumnuclear Disk in the Galactic Center to thick, obscuring tori of AGNs -- Modelling the molecular emission of a parsec-scale torus as found in NGC1068

Authors: B. Vollmer, R. I. Davies, P. Gratier, Th. Lizee, M. Imanishi, J. F. Gallimore, C. M. V. Impellizzeri, S. Garcia-Burillo, F. Le Petit

Abstract: The accretion rates needed to fuel the central black hole in a galaxy can be achieved via viscous torques in thick disks and rings, which can be resolved by millimetre interferometry within the inner ~20pc of the active galaxy NGC1068 at comparable scales and sensitivity to single dish observations of the Circumnuclear Disk (CND) in the Galactic Center. To interpret observations of these regions a… ▽ More The accretion rates needed to fuel the central black hole in a galaxy can be achieved via viscous torques in thick disks and rings, which can be resolved by millimetre interferometry within the inner ~20pc of the active galaxy NGC1068 at comparable scales and sensitivity to single dish observations of the Circumnuclear Disk (CND) in the Galactic Center. To interpret observations of these regions and determine the physical properties of their gas distribution, we present a modelling effort that includes (i) a simple dynamical simulations involving partially inelastic collisions between disk gas clouds, (ii) an analytical model of a turbulent clumpy gas disk calibrated by the dynamical model and observations, (iii) local turbulent and cosmic ray gas heating and cooling via H2O, H2, and CO emission, and (iv) determination of the molecular abundances. We also consider photodissociation regions (PDR) where gas is directly illuminated by the central engine. We compare the resulting model datacubes of the CO, HCN, HCO+, and CS brightness temperatures to available observations. In both cases the kinematics can be explained by one or two clouds colliding with a pre-existing ring, in a prograde sense for the CND and retrograde for NGC1068. And, with only dense disk clouds, the line fluxes can be reproduced to within a factor of about two. To avoid self-absorption of the intercloud medium, turbulent heating at the largest scales, comparable to the disk height, has to be decreased by a factor of 50-200. Our models indicate that turbulent mechanical energy input is the dominant gas heating mechanism within the thick gas disks. In N1068, while the bulk of the AGN X-ray radiation is absorbed in a layer of Compton-thick gas inside the dust sublimation radius, the optical/UV radiation may enhance the molecular line emission from photodissociation regions by ~50% at the inner edge of the gas ring. △ Less

Submitted 29 June, 2022; originally announced June 2022.

Comments: 36 pages, 32 figures, accepted for publication in A&A

Journal ref: A&A 665, A102 (2022)

Extremely Buried Nucleus of IRAS 17208$-$0014 Observed at Sub-Millimeter and Near-Infrared Wavelengths

Authors: Shunsuke Baba, Masatoshi Imanishi, Takuma Izumi, Taiki Kawamuro, Dieu D. Nguyen, Takao Nakagawa, Naoki Isobe, Shusuke Onishi, Kosei Matsumoto

Abstract: The ultraluminous infrared galaxy IRAS 17208$-$0014 is a late-stage merger that hosts a buried active galactic nucleus (AGN). To investigate its nuclear structure, we performed high spatial resolution ($\sim0.\!\!^{\prime\prime}04\sim32\,\mathrm{pc}$) Atacama Large Millimeter/submillimeter Array (ALMA) observations in Band 9 ($\sim$450\,\micron\ or $\sim$660\,GHz), along with near-infrared AKARI s… ▽ More The ultraluminous infrared galaxy IRAS 17208$-$0014 is a late-stage merger that hosts a buried active galactic nucleus (AGN). To investigate its nuclear structure, we performed high spatial resolution ($\sim0.\!\!^{\prime\prime}04\sim32\,\mathrm{pc}$) Atacama Large Millimeter/submillimeter Array (ALMA) observations in Band 9 ($\sim$450\,\micron\ or $\sim$660\,GHz), along with near-infrared AKARI spectroscopy in 2.5--5.0\,\micron. The Band 9 dust continuum peaks at the AGN location, and toward this position CO($J$=6--5) and CS($J$=14--13) are detected in absorption. Comparison with non-local thermal equilibrium calculations indicates that, within the central beam ($r\sim20\,\mathrm{pc}$), there exists a concentrated component that is dense ($10^7\,\mathrm{cm}^{-2}$) and warm ($>$200\,K) and has a large column density ($N_\mathrm{H_2}>10^{23}\,\mathrm{cm}^{-2}$). The AKARI spectrum shows deep and broad CO ro-vibrational absorption at 4.67\,\micron. Its band profile is well reproduced with a similarly dense and large column but hotter ($\sim$1000\,K) gas. The region observed through absorption in the near-infrared is highly likely in the nuclear direction, as in the sub-millimeter, but with a narrower beam including a region closer to the nucleus. The central component is considered to possess a hot structure where vibrationally excited HCN emission originates. The most plausible heating source for the gas is X-rays from the AGN. The AKARI spectrum does not show other AGN signs in 2.5--4\,\micron, but this absence may be usual for AGNs buried in a hot mid-infrared core. Besides, based on our ALMA observations, we relate various nuclear structures of IRAS 17208$-$0014 that have been proposed in the literature. △ Less

Submitted 22 February, 2022; originally announced February 2022.

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

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

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

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

Submitted 2 December, 2021; originally announced December 2021.

Comments: 22 pages, 8 figures, accepted in ApJ

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

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

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

Submitted 24 November, 2021; originally announced November 2021.

Comments: Accepted for publication in The Astrophysical Journal Supplement Series

ALMA Sub-arcsec-resolution 183 GHz H2O and Dense Molecular Line Observations of Nearby Ultraluminous Infrared Galaxies

Authors: Masatoshi Imanishi, Kouichiro Nakanishi, Takuma Izumi, Shunsuke Baba

Abstract: We present the results of ALMA $\sim$2 mm, $\lesssim$1$''$-resolution observations of ten (ultra)luminous infrared galaxies ([U]LIRGs; infrared luminosity $\gtrsim$ 10$^{11.7}$L$_{\odot}$) at $z <$ 0.15, targeting dense ($>$10$^{4}$ cm$^{-3}$) molecular (HCN, HCO$^{+}$, and HNC J=2-1) and 183 GHz H$_{2}$O 3$_{1,3}$-2$_{2,0}$ emission lines. Higher HCN to HCO$^{+}$ J=2-1 flux ratios are observed in… ▽ More We present the results of ALMA $\sim$2 mm, $\lesssim$1$''$-resolution observations of ten (ultra)luminous infrared galaxies ([U]LIRGs; infrared luminosity $\gtrsim$ 10$^{11.7}$L$_{\odot}$) at $z <$ 0.15, targeting dense ($>$10$^{4}$ cm$^{-3}$) molecular (HCN, HCO$^{+}$, and HNC J=2-1) and 183 GHz H$_{2}$O 3$_{1,3}$-2$_{2,0}$ emission lines. Higher HCN to HCO$^{+}$ J=2-1 flux ratios are observed in some, but not all, AGN-important ULIRGs than in starburst-classified sources. We detect 183 GHz H$_{2}$O emission in almost all AGN-important ULIRGs, and elevated H$_{2}$O emission is found in two sources with elevated HCN J=2-1 emission, relative to HCO$^{+}$ J=2-1. Except one ULIRG (the Superantennae), the H$_{2}$O emission largely comes from the entire nuclear regions ($\sim$1 kpc), rather than AGN-origin megamaser at the very center ($<<$1 kpc). Nuclear ($\sim$1 kpc) dense molecular gas mass derived from HCO$^{+}$ J=2-1 luminosity is $\gtrsim$a few $\times$ 10$^{8}$M$_{\odot}$, and its depletion time is estimated to be $\gtrsim$10$^{6}$ yr in all sources. Vibrationally excited J=2-1 emission lines of HCN and HNC are detected in a few (U)LIRGs, but those of HCO$^{+}$ are not. It is suggested that in mid-infrared-radiation-exposed innermost regions around energy sources, HCO$^{+}$ and HNC are substantially less abundant than HCN. In our ALMA $\sim$2 mm data of ten (U)LIRGs, two continuum sources are serendipitously detected within $\sim$10$''$, which are likely to be an infrared luminous dusty galaxy at $z >$ 1 and a blazar. △ Less

Submitted 17 November, 2021; originally announced November 2021.

Comments: 50 pages, 16 Figures. ApJ accepted

Second Epoch ALMA Observations of 321 GHz Water Maser Emission in NGC 4945 and the Circinus Galaxy

Authors: Yoshiaki Hagiwara, Shinji Horiuchi, Masatoshi Imanishi, Philip G. Edwards

Abstract: We present the results of second epoch ALMA observations of 321 GHz H$_2$O emission toward two nearby active galactic nuclei, NGC 4945 and the Circinus galaxy, together with Tidbinbilla 70-m monitoring of their 22 GHz H$_2$O masers. The two epoch ALMA observations show that the strengths of the 321 GHz emission are variable by a factor of at least a few, confirming a maser origin. In the second ep… ▽ More We present the results of second epoch ALMA observations of 321 GHz H$_2$O emission toward two nearby active galactic nuclei, NGC 4945 and the Circinus galaxy, together with Tidbinbilla 70-m monitoring of their 22 GHz H$_2$O masers. The two epoch ALMA observations show that the strengths of the 321 GHz emission are variable by a factor of at least a few, confirming a maser origin. In the second epoch, 321 GHz maser emission from NGC 4945 was not detected, while for the Circinus galaxy the flux density significantly increased and the velocity gradient and dispersion have been measured. With the velocity gradient spanning $\sim$110 km/s, we calculate the disk radius to be $\sim$28 pc, assuming disk rotation around the nucleus. We also estimate the dynamical mass within the central 28 pc to be 4.3 $\times$ 10$^8$ M$_{\odot}$, which is significantly larger than the larger scale dynamical mass, suggesting the velocity gradient does not trace circular motions on that scale. The overall direction of the velocity gradient and velocity range of the blueshifted features are largely consistent with those of the 22 GHz maser emission in a thin disk with smaller radii of 0.1--0.4 pc and molecular outflows within $\sim$1 pc from the central engine of the galaxy, implying that the 321 GHz masers could trace part of the circumnuclear disk or the nuclear outflows. △ Less

Submitted 1 November, 2021; originally announced November 2021.

Comments: 13 pages, 9 figures, accepted for publication in ApJ

The MBHBM$^{\star}$ Project -- II. Molecular Gas Kinematics in the Lenticular Galaxy NGC 3593 Reveal a Supermassive Black Hole

Authors: Dieu D. Nguyen, Martin Bureau, Sabine Thater, Kristina Nyland, Mark den Brok, Michelle Cappellari, Timothy A. Davis, Jenny E. Greene, Nadine Neumayer, Masatoshi Imanishi, Takuma Izumi, Taiki Kawamuro, Shunsuke Baba, Phuong M. Nguyen, Satoru Iguchi, Takafumi Tsukui, Lam N. T., Than Ho

Abstract: As part of the Measuring Black Holes in Below Milky Way-mass (M$^\star$) galaxies (MBHBM$^\star$) Project, we present a dynamical measurement of the supermassive black hole (SMBH) mass in the nearby lenticular galaxy NGC 3593, using cold molecular gas $^{12}$CO(2-1) emission observed at an angular resolution of $\approx0''.3$ ($\approx10$ pc) with the Atacama Large Millimeter/submillimeter Array (… ▽ More As part of the Measuring Black Holes in Below Milky Way-mass (M$^\star$) galaxies (MBHBM$^\star$) Project, we present a dynamical measurement of the supermassive black hole (SMBH) mass in the nearby lenticular galaxy NGC 3593, using cold molecular gas $^{12}$CO(2-1) emission observed at an angular resolution of $\approx0''.3$ ($\approx10$ pc) with the Atacama Large Millimeter/submillimeter Array (ALMA). Our ALMA observations reveal a circumnuclear molecular gas disc (CND) elongated along the galaxy major axis and rotating around the SMBH. Using dynamical modelling, the molecular gas kinematics allow us to infer a SMBH mass $M_{\rm BH}=2.40_{-1.05}^{+1.87}\times10^6$ M$_\odot$ (only statistical uncertainties at the $3σ$ level). We also detect a massive core of cold molecular gas (CMC) of mass $M_{\rm CMC}=(5.4\pm1.2)\times10^6$ M$_\odot$ and effective (half-mass) radius $r_{\rm CMC,e}=11.2\pm2.8$ pc, co-spatial with a nuclear star cluster (NSC) of mass $M_{\rm NSC}=(1.67\pm0.48)\times10^7$ M$_\odot$ and effective radius $r_{\rm NSC,e}=5.0\pm1.0$~pc (or $0''.15\pm0''.03$). The mass profiles of the CMC and NSC are well described by Sérsic functions with indices $1-1.4$. Our $M_{\rm BH}$ and $M_{\rm NSC}$ estimates for NGC 3593 agree well with the recently compiled $M_{\rm BH}$-$M_{\rm NSC}$ scaling relation. Although the $M_{\rm NSC}$ uncertainty is twice the inferred $M_{\rm BH}$, the rapid central rise of the rotation velocities of the CND (as the radius decreases) clearly suggests a SMBH. Indeed, our dynamical models show that even if $M_{\rm NSC}$ is at the upper end of its allowed range, the evidence for a black hole does not vanish, but remains with a lower limit of $M_{\rm BH}>3\times10^5$ M$_\odot$. △ Less

Submitted 16 October, 2021; originally announced October 2021.

Comments: 27 pages, 18 Figures, 8 Tables. Accepted to MNRAS

Hard X-ray Irradiation Potentially Drives Negative AGN Feedback by Altering Molecular Gas Properties

Authors: Taiki Kawamuro, Claudio Ricci, Takuma Izumi, Masatoshi Imanishi, Shunsuke Baba, Dieu D. Nguyen, Kyoko Onishi

Abstract: To investigate the role of active galactic nucleus (AGN) X-ray irradiation on the interstellar medium (ISM), we systematically analyzed Chandra and ALMA CO($J$=2-1) data for 26 ultra-hard X-ray ($>$ 10 keV) selected AGNs at redshifts below 0.05. While Chandra unveils the distribution of X-ray-irradiated gas via Fe-K$α$ emission, the CO($J$=2-1) observations reveal that of cold molecular gas. At hi… ▽ More To investigate the role of active galactic nucleus (AGN) X-ray irradiation on the interstellar medium (ISM), we systematically analyzed Chandra and ALMA CO($J$=2-1) data for 26 ultra-hard X-ray ($>$ 10 keV) selected AGNs at redshifts below 0.05. While Chandra unveils the distribution of X-ray-irradiated gas via Fe-K$α$ emission, the CO($J$=2-1) observations reveal that of cold molecular gas. At high resolutions $\lesssim$ 1 arcsec, we derive Fe-K$α$ and CO($J$=2-1) maps for the nuclear 2 arcsec region, and for the external annular region of 2 arcsec-4 arcsec, where 2 arcsec is $\sim$ 100-600 pc for most of our AGNs. First, focusing on the external regions, we find the Fe-K$α$ emission for six AGNs above 2$σ$. Their large equivalent widths ($\gtrsim$ 1 keV) suggest a fluorescent process as their origin. Moreover, by comparing 6-7 keV/3-6 keV ratio, as a proxy of Fe-K$α$, and CO($J$=2-1) images for three AGNs with the highest significant Fe-K$α$ detections, we find a possible spatial separation. These suggest the presence of X-ray-irradiated ISM and the change in the ISM properties. Next, examining the nuclear regions, we find that (1) The 20-50 keV luminosity increases with the CO($J$=2-1) luminosity. (2) The ratio of CO($J$=2-1)-to-HCN($J$=1-0) luminosities increases with 20-50 keV luminosity, suggesting a decrease in the dense gas fraction with X-ray luminosity. (3) The Fe-K$α$-to-X-ray continuum luminosity ratio decreases with the molecular gas mass. This may be explained by a negative AGN feedback scenario: the mass accretion rate increases with gas mass, and simultaneously, the AGN evaporates a portion of the gas, which possibly affects star formation. △ Less

Submitted 20 September, 2021; originally announced September 2021.

Comments: 50 pages (27 pages for the main body), 26 figures, accepted for publication in ApJS

Synchronized Co-evolution between Supermassive Black Holes and Galaxies Over the Last Seven Billion Years as Revealed by the Hyper Suprime-Cam

Authors: Junyao Li, John D. Silverman, Xuheng Ding, Michael A. Strauss, Andy Goulding, Malte Schramm, Hassen M. Yesuf, Mouyuan Sun, Yongquan Xue, Simon Birrer, Jingjing Shi, Yoshiki Toba, Tohru Nagao, Masatoshi Imanishi

Abstract: We measure the evolution of the $M_{\rm BH}-M_*$ relation using 584 uniformly-selected SDSS quasars at $0.2<z<0.8$. The black-hole masses ($M_{\rm BH}$) are derived from the single-epoch virial mass estimator using the H$β$ emission line, and span the range $7.0<{\rm log}\,M_{\rm BH}/M_\odot<9.5$. The host-galaxy stellar masses ($M_*$), which cover the interval… ▽ More We measure the evolution of the $M_{\rm BH}-M_*$ relation using 584 uniformly-selected SDSS quasars at $0.2<z<0.8$. The black-hole masses ($M_{\rm BH}$) are derived from the single-epoch virial mass estimator using the H$β$ emission line, and span the range $7.0<{\rm log}\,M_{\rm BH}/M_\odot<9.5$. The host-galaxy stellar masses ($M_*$), which cover the interval $10.0<{\rm log}\,M_* / M_\odot <11.5$, are determined by performing two-dimensional quasar-host decomposition of the Hyper Suprime-Cam images and spectral energy distribution fitting. To quantify sample-selection biases and measurement uncertainties on the mass terms, a mock quasar sample is constructed to jointly constrain the redshift evolution of the $M_{\rm BH}-M_*$ relation and its intrinsic scatter ($σ_μ$) through forward modeling. We find that the level of evolution is degenerate with $σ_μ$, such that both positive mild evolution (i.e, $M_{\rm BH}/M_*$ increases with redshift) with a small $σ_μ$, and negative mild evolution with a larger $σ_μ$ are consistent with our data. The posterior distribution of $σ_μ$ enables us to put a strong constraint on the intrinsic scatter of the $M_{\rm BH}-M_*$ relation, which has a best inference of $0.25_{-0.04}^{+0.03}$ dex, consistent with the local value. The redshift evolution of the $M_{\rm BH}-M_*$ relation relative to the local relation is constrained to be $(1+z)^{0.12_{-0.27}^{+0.28}}$, in agreement with no significant evolution since $z\sim0.8$. The tight and non-evolving $M_{\rm BH}-M_*$ relation is suggestive of a coupling through AGN feedback or/and a common gas supply at work, thus restricting the mass ratio of galaxies and their black holes to a limited range. Given the considerable stellar disk component, the $M_{\rm BH}-M_{\rm bulge}$ relation may evolve as previously seen at higher redshifts. △ Less

Submitted 6 September, 2021; originally announced September 2021.

Comments: 17 pages, 6 figures, accepted by ApJ

Comprehensive Broadband X-ray and Multiwavelength Study of Active Galactic Nuclei in Local 57 Ultra/luminous Infrared Galaxies Observed with NuSTAR and/or Swift/BAT

Authors: Satoshi Yamada, Yoshihiro Ueda, Atsushi Tanimoto, Masatoshi Imanishi, Yoshiki Toba, Claudio Ricci, George C. Privon

Abstract: We perform a systematic X-ray spectroscopic analysis of 57 local ultra/luminous infrared galaxy systems (containing 84 individual galaxies) observed with Nuclear Spectroscopic Telescope Array and/or Swift/BAT. Combining soft X-ray data obtained with Chandra, XMM-Newton, Suzaku and/or Swift/XRT, we identify 40 hard ($>$10 keV) X-ray detected active galactic nuclei (AGNs) and constrain their torus p… ▽ More We perform a systematic X-ray spectroscopic analysis of 57 local ultra/luminous infrared galaxy systems (containing 84 individual galaxies) observed with Nuclear Spectroscopic Telescope Array and/or Swift/BAT. Combining soft X-ray data obtained with Chandra, XMM-Newton, Suzaku and/or Swift/XRT, we identify 40 hard ($>$10 keV) X-ray detected active galactic nuclei (AGNs) and constrain their torus parameters with the X-ray clumpy torus model XCLUMPY (Tanimoto et al. 2019). Among the AGNs at $z < 0.03$, for which sample biases are minimized, the fraction of Compton-thick ($N_{\rm H} \geq 10^{24}$ cm$^{-2}$) AGNs reaches 64$^{+14}_{-15}$% (6/9 sources) in late mergers, while 24$^{+12}_{-10}$% (3/14 sources) in early mergers, consistent with the tendency reported by Ricci et al. (2017). We find that the bolometric AGN luminosities derived from the infrared data increase, but the X-ray to bolometric luminosity ratios decrease, with merger stage. The X-ray weak AGNs in late mergers ubiquitously show massive outflows at sub-pc to kpc scales. Among them, the most luminous AGNs ($L_{\rm bol,AGN} \sim 10^{46}$ erg s$^{-1}$) have relatively small column densities of $\lesssim$10$^{23}$ cm$^{-2}$ and almost super-Eddington ratios ($λ_{\rm Edd} \sim$ 1.0). Their torus covering factors ($C_{\rm T}^{\rm (22)} \sim 0.6$) are larger than those of Swift/BAT selected AGNs with similarly high Eddington ratios. These results suggest a scenario that, in the final stage of mergers, multiphase strong outflows are produced due to chaotic quasi-spherical inflows and the AGN becomes extremely X-ray weak and deeply buried due to obscuration by inflowing and/or outflowing material. △ Less

Submitted 22 July, 2021; originally announced July 2021.

Comments: 85 pages, 88 figures, accepted in ApJS

The Galaxy Activity, Torus and Outflow Survey (GATOS): II. Torus and polar dust emission in nearby Seyfert galaxies

Authors: A. Alonso-Herrero, S. García-Burillo, S. F. Hoenig, I. García-Bernete, C. Ramos Almeida, O. González-Martín, E. López-Rodríguez, P. G. Boorman, A. J. Bunker, L. Burtscher, F. Combes, R. Davies, T. Díaz-Santos, P. Gandhi, B. García-Lorenzo, E. K. S. Hicks, L. K. Hunt, K. Ichikawa, M. Imanishi, T. Izumi, A. Labiano, N. A. Levenson, C. Packham, M. Pereira-Santaella, C. Ricci , et al. (8 additional authors not shown)

Abstract: We compare mid-IR and ALMA far-IR images of 12 nearby Seyferts selected from GATOS. The mid-IR unresolved emission contributes more than 60% of the nuclear emission in most galaxies. By contrast, the ALMA 870micron continuum emission is mostly resolved and typically along the torus equatorial direction (Paper I, Garcia-Burillo et al. 2021). The Eddington ratios and nuclear hydrogen column densitie… ▽ More We compare mid-IR and ALMA far-IR images of 12 nearby Seyferts selected from GATOS. The mid-IR unresolved emission contributes more than 60% of the nuclear emission in most galaxies. By contrast, the ALMA 870micron continuum emission is mostly resolved and typically along the torus equatorial direction (Paper I, Garcia-Burillo et al. 2021). The Eddington ratios and nuclear hydrogen column densities NH of half the sample are favorable to launching polar and/or equatorial dusty winds, according to simulations. Six show mid-IR extended emission in the polar direction as traced by the NLR and perpendicular to the ALMA emission. In a few, the nuclear NH might be too high to uplift large quantities of dusty material along the polar direction. Five galaxies have low NH and/or Eddington ratios and thus polar dusty winds are not likely. We generate new CAT3D-WIND disk-wind model images. At low wind-to-disk cloud ratios the far-IR model images have disk- and ring-like morphologies. The X-shape associated with dusty winds is seen better in the far-IR at intermediate-high inclinations for the extended-wind configurations. In most models, the mid-IR emission comes from the inner part of the disk/cone. Extended bi-conical and one-sided polar mid-IR emission is seen in extended-wind configurations and high wind-to-disk cloud ratios. When convolved to our resolution, the model images reproduce qualitative aspects of the observed morphologies. Low-intermediate wind-to-disk ratios are required to account for the large fractions of unresolved mid-IR emission. This work and Paper I provide observational support for the torus+wind scenario. The wind component is more relevant at high Eddington ratios and/or AGN luminosities, and polar dust emission is predicted at NH of up to $10^{24}$cm$^{-2}$. The torus/disk component, on the other hand, prevails at low luminosities and/or Eddington ratios. (Abridged) △ Less

Submitted 2 July, 2021; v1 submitted 1 July, 2021; originally announced July 2021.

Comments: 27 pages, 27 figures. Accepted for publication in A&A

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)

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

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

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

Submitted 25 June, 2021; originally announced June 2021.

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

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

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

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

Submitted 21 May, 2021; originally announced May 2021.

Comments: 36 pages, 21 figures, 4 tables

Molecular gas kinematics in the nuclear region of nearby Seyfert galaxies with ALMA

Authors: A. Bewketu Belete, P. Andreani, J. A. Fernández-Ontiveros, E. Hatziminaoglou, F. Combes, M. Sirressi, R. Slater, C. Ricci, K. Dasyra, C. Cicone, S. Aalto, L. Spinoglio, M. Imanishi, J. R. De Medeiros

Abstract: We present the analysis of the molecular gas in the nuclear regions of NGC 4968, NGC 4845, and MCG-06-30-15, with the help of ALMA observations of the CO(2-1) emission line. The aim is to determine the kinematics of the gas in the central (~ 1 kpc) region. We use the 3D-Based Analysis of Rotating Object via Line Observations ($^{3D}$BAROLO) and DiskFit softwares. Circular motions dominate the kine… ▽ More We present the analysis of the molecular gas in the nuclear regions of NGC 4968, NGC 4845, and MCG-06-30-15, with the help of ALMA observations of the CO(2-1) emission line. The aim is to determine the kinematics of the gas in the central (~ 1 kpc) region. We use the 3D-Based Analysis of Rotating Object via Line Observations ($^{3D}$BAROLO) and DiskFit softwares. Circular motions dominate the kinematics of the gas in the central discs, mainly in NGC 4845 and MCG-06-30-15, however there is a clear evidence of non-circular motions in the central ($\sim$ 1 kpc) region of NGC 4845 and NGC 4968. The strongest non-circular motion is detected in the inner disc of NGC 4968 with velocity $\sim 115\, \rm{km\,s^{-1}}$. The bisymmetric model is found to give the best-fit for NGC 4968 and NGC 4845. If the dynamics of NGC 4968 is modeled as a corotation pattern just outside of the bar, the bar pattern speed turns out to be at $Ω_b$ = $52\, \rm{km\,s^{-1}\,kpc^{-1}}$ the corotation is set at 3.5 kpc and the inner Lindblad resonance (ILR) ring at R = 300pc corresponding to the CO emission ring. The 1.2 mm ALMA continuum is peaked and compact in NGC 4968 and MCG-06-30-15, but their CO(2-1) has an extended distribution. Allowing the CO-to-H$_{2}$ conversion factor $α_{CO}$ between 0.8 and 3.2, typical of nearby galaxies of the same type, the molecular mass M(H$_{2}$) is estimated to be $\sim 3-12\times 10^{7} ~{\rm M_\odot}$ (NGC 4968), $\sim 9-36\times 10^{7}~ {\rm M_\odot}$ (NGC 4845), and $\sim 1-4\times 10^{7}~ {\rm M_\odot}$ (MCG-06-30-15). We conclude that the observed non-circular motions in the disc of NGC 4968 and likely that seen in NGC 4845 is due to the presence of the bar in the nuclear region. At the current spectral and spatial resolution and sensitivity we cannot claim any strong evidence in these sources of the long sought feedback/feeding effect due to the AGN presence. △ Less

Submitted 14 May, 2021; originally announced May 2021.

Comments: 14 pages, 11 figures, 3 tables, accepted for publication in the Astronomy & Astrophysics (A&A) journal

Journal ref: A&A 654, A24 (2021)

The Galaxy Activity, Torus and Outflow Survey (GATOS) I. ALMA images of dusty molecular tori in Seyfert galaxies

Authors: S. García-Burillo, A. Alonso-Herrero, C. Ramos Almeida, O. González-Martín, F. Combes, A. Usero, S. Hönig, M. Querejeta, E. K. S. Hicks, L. K. Hunt, D. Rosario, R. Davies, P. G. Boorman, A. J. Bunker, L. Burtscher, L. Colina, T. Díaz-Santos, P. Gandhi, I. García-Bernete, B. García-Lorenzo, K. Ichikawa, M. Imanishi, T. Izumi, A. Labiano, N. A. Levenson , et al. (10 additional authors not shown)

Abstract: We present the first results of the Galaxy Activity, Torus and Outflow Survey (GATOS), a project aimed at understanding the properties of the dusty molecular tori and their connection to the host galaxy in nearby Seyfert galaxies. Our project expands the range of AGN luminosities and Eddington ratios covered by previous surveys of Seyferts conducted by ALMA and allows us to study the gas feeding a… ▽ More We present the first results of the Galaxy Activity, Torus and Outflow Survey (GATOS), a project aimed at understanding the properties of the dusty molecular tori and their connection to the host galaxy in nearby Seyfert galaxies. Our project expands the range of AGN luminosities and Eddington ratios covered by previous surveys of Seyferts conducted by ALMA and allows us to study the gas feeding and feedback cycle in a combined sample of 19 Seyferts. We used ALMA to obtain new images of the emission of molecular gas and dust using the CO(3-2) and HCO+(4-3) lines as well as their underlying continuum emission at 870 microns with high spatial resolutions (0.1'' ~ 7 - 13 pc) in the CND of 10 nearby (D < 28 Mpc) Seyfert galaxies. Our new ALMA observations detect 870 micron continuum and CO line emission from spatially resolved disks located around the AGN in all the sources. The bulk of the continuum flux can be accounted for by thermal emission from dust in the majority of the targets. For most of the sources the disks show a preponderant orientation perpendicular to the AGN wind axes, as expected for dusty molecular tori. The median diameters and molecular gas masses of the tori are ~ 42 pc, and ~ 6 x 10**5 Msun, respectively. We find a positive correlation between the line-of-sight gas column densities responsible for the absorption of X-rays and the molecular gas column densities derived from CO towards the AGN in our sources. The radial distributions of molecular gas in the CND of our combined sample show signs of nuclear-scale molecular gas deficits. We also detect molecular outflows in the sources that show the most extreme nuclear-scale gas deficits in our sample. These observations find for the first time supporting evidence that the imprint of AGN feedback is more extreme in higher luminosity and/or higher Eddington ratio Seyfert galaxies. △ Less

Submitted 22 June, 2021; v1 submitted 20 April, 2021; originally announced April 2021.

Comments: 41 pages, new version revised and accepted by A&A on the 22nd of June 2021

Journal ref: A&A 652, A98 (2021)