Search

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 4422$. Those obscured and radio AGN reside in $0.4, Comment: 20 pages, 11 figures, accepted for publication in A&A

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., 84 pages (5 tables/29 figures in the main text and 8 tables/18 figures in the Appendix), accepted in ApJS

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., 7 pages, 3 figures, accepted in PASJ

We present the results of our multifrequency observations for 10 ultraluminous infrared galaxies (ULIRGs) made by the Karl G. Jansky Very Large Array at 1.4, 5.5, 9.0, and 14.0 GHz. Our sample is selected from ULIRGs whose active galactic nuclei (AGNs) are not found at optical wavelengths (∼70 per cent of the total ULIRGs), but whose presence is suggested by mid-infrared or submillimetre observations (>50 per cent of the non-AGN ULIRGs at optical wavelengths). The statistical properties of the targets are similar to those of the entire ULIRG sample, which implies that ULIRGs have common radiative processes regardless of the presence of optical AGNs, and thus AGNs might equally contribute to the radio emission of every ULIRG. Although their spectra are mainly explained by starbursts and/or merger activity, some individual sources suggest contributions from AGNs. IRAS 00188−0856, whose optical morphology is not disturbed, shows a large non-thermal fraction and a spectral break at high frequency, which can be explained by synchrotron ageing of non-thermal plasma emitted from AGNs. In addition, we find 100-kpc-scale extended emission associated with IRAS 01004−2237. The two-sided morphology and absence of extended X-ray emission suggest that this system is not induced by a merger in a cluster but originates from AGN activity.

We investigate spatial distributions of iron Kα (Fe-Kα) lines in the cental 100 pc of the Seyfert 2 galaxy NGC 1068 observed with Chandra. The spatial distributions of Fe-Kα lines, neutral and highly ionized, around the center of the galactic nucleus are not isotropic, as consistently confirmed in both image and spectral analyses. The hydrogen number density of the gas clouds responsible for the neutral Fe-Kα line emission is estimated to be 102–103 cm−3 for the sampled regions near the galactic core. The photo-ionization model, where iron is assumed to be ionized by X-rays from the galactic nucleus, yields ionization parameters lower than 19 for these clouds. The range of this ionization parameter is two or three orders of magnitude lower than the theoretically expected value to produce the observed helium-like Fe-Kα line intensities. Therefore, the photo-ionization model is excluded from the explanation of the amount of highly ionized iron that is responsible for the observed Fe-Kα lines. Also, we find anti-correlation in the spatial distributions between the molecular cloud in the area observed with ALMA and that of the Fe-Kα lines, including that from neutral iron. We suggest that X-ray iron-line and radio molecular cloud observations are complementary to probe the distribution of matters in the central regions around the cores of active galactic nuclei.

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., 28 pages, 16 figures, accepted for publication in ApJ

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 $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 and a blazar., Comment: 50 pages, 16 Figures. ApJ accepted

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., 13 pages, 9 figures, accepted for publication in ApJ

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 $\Omega_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 $\alpha_{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., Comment: 14 pages, 11 figures, 3 tables, accepted for publication in the Astronomy & Astrophysics (A&A) journal

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$\alpha$ emission, the CO($J$=2-1) observations reveal that of cold molecular gas. At high resolutions $\lesssim$ 1 arcsec, we derive Fe-K$\alpha$ 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$\alpha$ emission for six AGNs above 2$\sigma$. 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$\alpha$, and CO($J$=2-1) images for three AGNs with the highest significant Fe-K$\alpha$ 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$\alpha$-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., Comment: 50 pages (27 pages for the main body), 26 figures, accepted for publication in ApJS

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., 85 pages, 88 figures, accepted in ApJS

Onoue, Masafusa, et al., We present deep Keck/MOSFIRE near-infrared spectroscopy of a strong Lyα 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 narrow (FWHM < 500 km s-1) and prominent (LLyα > 1044 erg s-1) Lyα emission lines at absolute 1450 Å continuum magnitudes of M1450∼-22 mag. Its rest-frame Lyα equivalent width (EW) is 370 ± 30 Å. In the 2 hr Keck/MOSFIRE spectrum in Y band, the high-ionization C IV λλ1548,1550 doublet emission line was clearly detected with = - FWHM 120+20 20 km s-1 and a total rest-frame EW of - 37+5 6 Å. We also report the detection of weak continuum emission, and the tentative detection of O III]λλ1661,1666 in the 4 hr J-band spectrum. Judging from the UV magnitude, line widths, luminosities, and EWs of Lyα and C IV, 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 C IV emitters seen in z = 6-7 galaxy samples. Further observations of other metal emission lines in the rest-frame UV or optical, along with X-ray follow-up observations of the z = 6-7 narrow-line quasars, are needed for more robust diagnostics and to determine their nature., Funding was contributed by the FIRST program from Japanese Cabinet Office, the Ministry of Education, Culture, Sports, Science and Technology (MEXT), the Japan Society for the Promotion of Science (JSPS), Japan Science and Technology Agency (JST), the Toray Science Foundation, NAOJ, Kavli IPMU, KEK, ASIAA, and Princeton University. This work was supported by JSPS KAKENHI grant Nos. 17H04830 (YM), 20H01949 (TN), 18J01050 (YT). M.O. was supported by the ERC Advanced Grant 740246 Cosmic gas. Y. M. was supported by the Mitsubishi Foundation grant No. 30140. T.H. was supported by Leading Initiative for Excellent Young Researchers, MEXT, Japan (HJH02007).

We present the physical properties of Lyα emitters (LAEs) in a “DLA-concentrated regions” where there are 3 or more DLA within (50 Mpc)3 cubic box. We observed LAEs in a DLA-concentrated region at z = 2.3, the J1230+34 field, with Subaru/Suprime-Cam. In the 50 Mpc scale, we found no deferences in properties of LAEs such as Lyα luminosity function in the DLA-concentrated region compared to other fields at similar redshift. On the other hand, we found a ∼10 Mpc scale LAE overdensity around a strong DLA with NHI = 1021.08 cm−2.

We present ALMA [CII] line and far-infrared (FIR) continuum observations of seven z > 6 low-luminosity quasars (M1450 > −25 mag) discovered by our on-going Subaru Hyper Suprime-Cam survey. The [CII] line was detected in all targets with luminosities of ∼(2−10) × 108 L⊙, about one order of magnitude smaller than optically luminous quasars. Also found was a wide scatter of FIR continuum luminosity, ranging from LFIR < 1011L⊙ to ∼2 × 1012L⊙. With the [CII]-based dynamical mass, we suggest that a significant fraction of low-luminosity quasars are located on or even below the local Magorrian relation, particularly at the massive end of the galaxy mass distribution. This is a clear contrast to the previous finding that luminous quasars tend to have overmassive black holes relative to the relation. Our result is expected to show a less-biased nature of the early co-evolution of black holes and their host galaxies.

Izumi, T., et al., We present Atacama Large Millimeter/submillimeter Array [C ii] 158 μm line and far-infrared (FIR) continuum emission observations toward HSC J120505.09-000027.9 (J1205-0000) at z = 6.72 with a beam size of ∼0.″8 × 0.″5 (or 4.1 kpc × 2.6 kpc), the most distant red quasar known to date. Red quasars are modestly reddened by dust and are thought to be in rapid transition from an obscured starburst to an unobscured normal quasar, driven by powerful active galactic nucleus (AGN) feedback that blows out a cocoon of interstellar medium. The FIR continuum of J1205-0000 is bright, with an estimated luminosity of L FIR ∼ 3 × 1012 L o˙. The [C ii] line emission is extended on scales of r ∼ 5 kpc, greater than that of the FIR continuum. The line profiles at the extended regions are complex and broad (FWHM ∼ 630-780 km s-1). Although it is not practical to identify the nature of this extended structure, possible explanations include (i) companion/merging galaxies and (ii) massive AGN-driven outflows. For the case of (i), the companions are modestly star-forming (∼10 M o˙ yr-1) but are not detected by our Subaru optical observations (y AB,5σ = 24.4 mag). For the case of (ii), our lower limit to the cold neutral outflow rate is ∼100 M o˙ yr-1. The outflow kinetic energy and momentum are both much lower than predicted in energy-conserving wind models, suggesting that the AGN feedback in this quasar is not capable of completely suppressing its star formation., Funding was contributed by the FIRST program of the Japanese Cabinet Office, the Ministry of Education, Culture, Sports, Science, and Technology, the Japan Society for the Promotion of Science, the Japan Science and Technology Agency, the Toray Science Foundation, NAOJ, Kavli IPMU, KEK, ASIAA, and Princeton University. T.I. is supported by an ALMA Japan Research Grant of the NAOJ ALMA Project, NAOJ-ALMA-249.

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., 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

Full list of authors: Falstad, N.; Aalto, S.; König, S.; Onishi, K.; Muller, S.; Gorski, M.; Sato, M.; Stanley, F.; Combes, F.; González-Alfonso, E.; Mangum, J. G.; Evans, A. S.; Barcos-Muñoz, L.; Privon, G. C.; Linden, S. T.; Díaz-Santos, T.; Martín, S.; Sakamoto, K.; Harada, N.; Fuller, G. A.; Gallagher, J. S.; van der Werf, P. P.; Viti, S.; Greve, T. R.; García-Burillo, S.; Henkel, C.; Imanishi, M.; Izumi, T.; Nishimura, Y.; Ricci, C.; Mühle, S., Context. Some luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs) host extremely compact (r < 100 pc) and dusty nuclei. The high extinction associated with large column densities of gas and dust toward these objects render them hard to detect at many wavelengths. The intense infrared radiation arising from warm dust in these sources can provide a significant fraction of the bolometric luminosity of the galaxy and is prone to excite vibrational levels of molecules such as HCN. This results in emission from the rotational transitions of vibrationally excited HCN (HCN-vib); the brightest emission is found in compact obscured nuclei (CONs; ςHCN-vib > 1 L⊙ pc-2 in the J = 3-2 transition). However, there have been no systematic searches for CONs, and it is unknown how common they are. Aims. We aim to establish how common CONs are in the local Universe (z < 0.08), and whether their prevalence depends on the luminosity or other properties of the host galaxy. Methods. We conducted an Atacama Large Millimeter/submillimeter Array survey of the rotational J = 3-2 transition of HCN-vib in a volume-limited sample of 46 far-infrared luminous galaxies. Results. Compact obscured nuclei are identified in 38-13+18% of the ULIRGs, 21-6+12% of the LIRGs, and 0-0+9% of the lower luminosity galaxies. We find no dependence on the inclination of the host galaxy, but strong evidence of lower IRAS 25 μm to 60 μm flux density ratios (f25/f60) in CONs (with the exception of one galaxy, NGC 4418) compared to the rest of the sample. Furthermore, we find that CONs have stronger silicate features (s9.7 μm), but similar polycyclic aromatic hydrocarbon equivalent widths (EQW6.2 μm) compared to other galaxies. Along with signatures of molecular inflows seen in the far-infrared in most CONs, submillimeter observations also reveal compact, often collimated, outflows. Conclusions. In the local Universe, CONs are primarily found in (U)LIRGs, in which they are remarkably common. As such systems are often highly disturbed, inclinations are difficult to estimate, and high-resolution continuum observations of the individual nuclei are required to determine if the CON phenomenon is related to the inclinations of the nuclear disks. Further studies of the in- A nd outflow properties of CONs should also be conducted to investigate how these are connected to each other and to the CON phenomenon. The lower f25/f60 ratios in CONs as well as the results for the mid-infrared diagnostics investigated (EQW6.2 μm and s9.7 μm) are consistent with the notion that large dust columns gradually shift the radiation from the hot nucleus to longer wavelengths, making the mid- A nd far-infrared "photospheres"significantly cooler than the interior regions. Finally, to assess the importance of CONs in the context of galaxy evolution, it is necessary to extend this study to higher redshifts where (U)LIRGs are more common. © 2021 ESO., This paper makes use of the following ALMA data: ADS/JAO.ALMA#2018.1.01344.S., ADS/JAO.ALMA#2015.1.00404.S, ADS/ JAO.ALMA#2017.1.00082.S, ADS/JAO.ALMA#2017.1.00057.S, ADS/JAO. ALMA#2013.1.00032.S, ADS/JAO.ALMA#2017.1.00598.S, ADS/JAO.ALMA#2015.1.00412.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), MOST and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Repub lic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. We acknowledge support from the Nordic ALMA Regional Centre (ARC) node based at Onsala Space Observatory. The Nordic ARC node is funded through Swedish Research Council grant No 2017-00648. S.A. gratefully acknowledges support from an ERC AdvancedGrant 789410 and from the Swedish Research Council. K.S. S.G-B. acknowledges support from the research projects PGC2018-094671-B-I00 (MCIU/AEI/FEDER, UE) and PID2019-106027GA C44 from the Spanish Ministerio de Ciencia e Innovación T.D-S. acknowledges support from the CASSACA and CONICYT fund CAS-CONICYT Call 2018. G.A.F acknowledges financial support from the State Agency for Research of the Spanish MCIU through the AYA2017-84390-C2-1-R grant (co-funded by FEDER) and through the “Center of Excellence Severo Ochoa” award for the Instituto de Astrofísica de Andalucia (SEV-2017-0709). T.R.G. acknowledges the Cosmic Dawn Center of Excellence funded by the Danish National Research Foundation under grant No. 140. E.G-A. is a Research Associate at the Harvard Smithsonian Center for Astrophysics, and thanks the Spanish Ministerio de Economía y Competitividad for support under projects ESP2017-86582-C4-1-R and PID2019-105552RB-C41. This research has made use of NASA’s Astrophysics Data System. This research has made use of the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion Labora tory, California Institute of Technology, under contract with the National Aeronautics and Space Administration.

Izumi, T., et al., We present Atacama Large Millimeter/submillimeter Array [C ii] 158 μm line and underlying far-infrared (FIR) continuum emission observations (0.″70 × 0.″56 resolution) toward HSC J124353.93+010038.5 (J1243+0100) at z = 7.07, the only low-luminosity (M 1450 > -25 mag) quasar currently known at z > 7. The FIR continuum is bright (1.52 mJy) and resolved with a total luminosity of L FIR = 3.5 × 1012 L o˙. The spatially extended component is responsible for ∼40% of the emission. The area-integrated [C ii] spectrum shows a broad wing (FWHM = 997 km s-1, L [C ii] = 1.2 × 109 L o˙), as well as a bright core (FWHM = 235 km s-1, L [C ii] = 1.9 × 109 L o˙). This wing is the first detection of a galactic-scale quasar-driven outflow (atomic outflow rate >447 M o˙ yr-1) at z > 7. The estimated large mass-loading factor of the total outflow (e.g., ⪆9 relative to the [C ii]-based star formation rate) suggests that this outflow will soon quench the star formation of the host. The core gas dynamics are governed by rotation, with a rotation curve suggestive of a compact bulge (∼3.3 × 1010 M o˙), although it is not yet spatially resolved. Finally, we found that J1243+0100 has a black hole mass-to-dynamical mass (and -to-bulge mass) ratio of ∼0.4% (∼1%), consistent with the local value within the uncertainties. Our results therefore suggest that the black hole-host coevolution relation is already in place at z ∼ 7 for this object., Funding was contributed by the FIRST program from the Japanese Cabinet Office, the Ministry of Education, Culture, Sports, Science and Technology (MEXT), the Japan Society for the Promotion of Science (JSPS), the Japan Science and Technology Agency (JST), the Toray Science Foundation, NAOJ, Kavli IPMU, KEK, ASIAA, and Princeton University. T.H. was supported by th

We present the results of ALMA band-5 (~170 GHz) observations of the merging ultraluminous infrared galaxy, the "Superantennae" (IRAS 19254-7245) at z=0.0617, which has been diagnosed as containing a luminous obscured active galactic nucleus (AGN). In addition to dense molecular line emission (HCN, HCO+, and HNC J = 2-1), we detect a highly luminous (~6e4Lsun) 183 GHz H2O 3(1,3)-2(2,0) emission line. We interpret the strong H2O emission as largely originating in maser amplification in AGN-illuminated dense and warm molecular gas, based on (1) the spatially compact (500 pc) dense molecular emission, and (2) a strikingly different velocity profile from, and (3) significantly elevated flux ratio relative to, dense molecular emission lines. H2O maser emission, other than the widely studied 22 GHz 6(1,6)-5(2,3) line, has been expected to provide important information on the physical properties of gas in the vicinity of a central mass-accreting supermassive black hole (SMBH), because of different excitation energy. We here demonstrate that with highly sensitive ALMA, millimetre 183 GHz H2O maser detection is feasible out to >270 Mpc, opening a new window to scrutinize molecular gas properties around a mass-accreting SMBH far beyond the immediately local universe., 6 pages, 2 figures. MNRAS Letters accepted

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., Accepted for publication in The Astrophysical Journal Supplement Series

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., 41 pages, new version revised and accepted by A&A on the 22nd of June 2021

We present a supermassive black hole (SMBH) mass measurement in the Seyfert 1 galaxy NGC7469 using Atacama Large Millimeter/submillimeter Array (ALMA) observations of the atomic-${\rm [CI]}$(1-0) and molecular-$^{12}$CO(1-0) emission lines at the spatial resolution of $\approx0.3$" (or $\approx$ 100 pc). These emissions reveal that NGC7469 hosts a circumnuclear gas disc (CND) with a ring-like structure and a two-arm/bi-symmetric spiral pattern within it, surrounded by a starbursting ring. The CND has a relatively low $\sigma/V\approx0.35$ ($r\sim0.5$") and $\sim0.19$ ($r>0.5"$), suggesting that the gas is dynamically settled and suitable for dynamically deriving the mass of its central source. As is expected from X-ray dominated region (XDR) effects that dramatically increase an atomic carbon abundance by dissociating CO molecules, we suggest that the atomic [CI](1-0) emission is a better probe of SMBH masses than CO emission in AGNs. Our dynamical model using the ${\rm [CI]}$(1-0) kinematics yields a $M_{\rm BH}=1.78^{+2.69}_{-1.10}\times10^7$M$_\odot$ and $M/L_{\rm F547M}=2.25^{+0.40}_{-0.43}$ (M$_\odot$/L$_\odot$). The model using the CO(1-0) kinematics also gives a consistent $M_{\rm BH}$ with a larger uncertainty, up to an order of magnitude, i.e.\ $M_{\rm BH}=1.60^{+11.52}_{-1.45}\times10^7$M$_\odot$. This newly dynamical $M_{\rm BH}$ is $\approx$ 2 times higher than the mass determined from the reverberation mapped (RM) method using emissions arising in the unresolved broad-line region (BLR). Given this new $M_{\rm BH}$, we are able to constrain the specific RM dimensionless scaling factor of $f=7.2^{+4.2}_{-3.4}$ for the AGN BLR in NGC7469. The gas within the unresolved BLR thus has a Keplerian virial velocity component and the inclination of $i\approx11.0^\circ$$_{-2.5}^{+2.2}$, confirming its face-on orientation in a Seyfert 1 AGN by assuming a geometrically thin BLR model., Comment: 22 pages, 16 figures, 7 tables. Accepted for publication on MNRAS

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\sigma$ 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\'{e}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$., Comment: 27 pages, 18 Figures, 8 Tables. Accepted to MNRAS

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$., Comment: 22 pages, 8 figures, accepted in ApJ

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

arXiv:2007.08685v2, We present the first discovery of dust-reddened quasars (red quasars) in the high-z universe (z > 5.6). This is a result from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which is based on the sensitive multi-band optical imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. We identified four red quasar candidates from the 93 spectroscopically confirmed high-z quasars in the SHELLQs sample, based on detections in the Wide-field Infrared Survey Explorer (WISE) data at 3.4 and 4.6 μm (rest-frame ∼5000–6500 Å). The amount of dust reddening was estimated with spectral energy distribution (SED) fits over optical and mid-infrared wavelengths. Two of the four candidates were found to be red quasars with dust reddening of E(B − V) > 0.1. The remaining SHELLQs quasars without individual WISE detections are significantly fainter in the WISE bands and bluer than the red quasars, although we did detect them in the W1 band in a stacked image. We also conducted the same SED fits for high-z optically-luminous quasars, but no red quasar was found. This demonstrates the power of Subaru HSC to discover high-z red quasars, which are fainter than the limiting magnitudes of past surveys in the rest-frame ultraviolet, due to dust extinction., We are grateful to the referee for his/her useful comments to improve this paper. YM was supported by the Japan Society for the Promotion of Science (JSPS) KAKENHI grant No. JP17H04830 and the Mitsubishi Foundation grant No. 30140. TI acknowledges supports from the JSPS grant No. JP17K14247. The Hyper Suprime-Cam (HSC) collaboration includes the astronomical communities of Japan and Taiwan, and Princeton University. The HSC instrumentation and software were developed by the National Astronomical Observatory of Japan (NAOJ), the Kavli Institute for the Physics and Mathematics of the Universe (Kavli IPMU), the University of Tokyo, the High Energy Accelerator Research Organization (KEK), the Academia Sinica Institute for Astronomy and Astrophysics in Taiwan (ASIAA), and Princeton University. Funding was contributed by the FIRST program from Japanese Cabinet Office, the Ministry of Education, Culture, Sports, Science and Technology (MEXT), the Japan Society for the Promotion of Science (JSPS), Japan Science and Technology Agency (JST), the Toray Science Foundation, NAOJ, Kavli IPMU, KEK, ASIAA, and Princeton University. The Pan-STARRS1 Surveys (PS1) have been made possible through contributions of the Institute for Astronomy, the University of Hawaii, the Pan-STARRS Project Office, the Max-Planck Society and its participating institutes, the Max Planck Institute for Astronomy, Heidelberg and the Max Planck Institute for Extraterrestrial Physics, Garching, The Johns Hopkins University, Durham University, the University of Edinburgh, Queen’s University Belfast, the Harvard-Smithsonian Center for Astrophysics, the Las Cumbres Observatory Global Telescope Network Incorporated, the National Central University of Taiwan, the Space Telescope Science Institute, the National Aeronautics and Space Administration under Grant No. NNX08AR22G issued through the Planetary Science Division of the NASA Science Mission Directorate, the National Science Foundation under Grant No. AST-1238877, the University of Maryland, and Eotvos Lorand University (ELTE). This paper makes use of software developed for the Large Synoptic Survey Telescope. We thank the LSST Project for making their code available as free software at 〈http://dm.lsst.org. This publication makes use of data products from the Wide-field Infrared Survey Explorer, which is a joint project of the University of California, Los Angeles, and the Jet Propulsion Laboratory/California Institute of Technology, funded by the National Aeronautics and Space Administration.