Search

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.

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

To study the high-transition dense-gas tracers and their relationships to the star formation of the inner $\sim$ 2 kpc circumnuclear region of NGC253, we present HCN $J=4-3$ and HCO$^+ J=4-3$ maps obtained with the James Clerk Maxwell Telescope (JCMT). With the spatially resolved data, we compute the concentration indices $r_{90}/r_{50}$ for the different tracers. HCN and HCO$^+$ 4-3 emission features tend to be centrally concentrated, which is in contrast to the shallower distribution of CO 1-0 and the stellar component. The dense-gas fraction ($f_\text{dense}$, traced by the velocity-integrated-intensity ratios of HCN/CO and HCO$^+$/CO) and the ratio $R_\text{31}$ (CO 3-2/1-0) decline towards larger galactocentric distances, but increase with higher SFR surface density. The radial variation and the large scatter of $f_\text{dense}$ and $R_\text{31}$ imply distinct physical conditions in different regions of the galactic disc. The relationships of $f_\text{dense}$ versus $\Sigma_\text{stellar}$, and SFE$_\text{dense}$ versus $\Sigma_\text{stellar}$ are explored. SFE$_\text{dense}$ increases with higher $\Sigma_\text{stellar}$ in this galaxy, which is inconsistent with previous work that used HCN 1-0 data. This implies that existing stellar components might have different effects on the high-$J$ HCN and HCO$^+$ than their low-$J$ emission. We also find that SFE$_\text{dense}$ seems to be decreasing with higher $f_\text{dense}$, which is consistent with previous works, and it suggests that the ability of the dense gas to form stars diminishes when the average density of the gas increases. This is expected in a scenario where only the regions with high-density contrast collapse and form stars., Comment: accepted to MNRAS

We present the first measurement of the mass of a supermassive black hole (SMBH) in the nearby double-barred spiral galaxy NGC 3504 as part of the Measuring Black Holes Below the Milky Way ($M_{\star}$) mass galaxies (MBHBM$_{\star}$) Project. Our analysis is based on Atacama Large Millimeter/submillimeter Array (ALMA) Cycle-5 observations of the ${\rm ^{12}CO(2-1)}$ emission line. NGC 3504 has a circumnuclear gas disk (CND), which has a relatively high-velocity dispersion of 30 \kms. Our dynamical models of the CND yield a \Mbh~of $M_{\rm BH}=1.02^{+0.18}_{-0.15}\times10^7$\Msun and a mass-to-light ratio in $H$-band of \ml$_{\rm H}=0.66^{+1.44}_{-0.65}$ (\Msun/\Lsun). This black hole (BH) mass is consistent with BH--galaxy scaling relations. We also detect a central deficiency in the ${\rm ^{12}CO(2-1)}$ integrated intensity map with a diameter of 2.7 pc at the putative position of the SMBH. However, this hole is filled by a dense gas tracer ${\rm CS(5-4)}$ that peaks at the galaxy center found in one of the three low-velocity-resolution continuum spectral correlators. The ${\rm CS(5-4)}$ line has the same kinematics with the ${\rm ^{12}CO(2-1)}$ line within the CND, suggesting that it is also an alternative transition for measuring the central \Mbh~in NGC 3504 probably more accurately than the current commonly used of ${\rm ^{12}CO(2-1)}$ due to its centralization., Comment: 25 pages, 15 Figures, 6 Tables, Submitted to ApJ

AGN-driven outflows are believed to play an important role in regulating the growth of galaxies mostly via negative feedback. However, their effects on their hosts are far from clear, especially for low and moderate luminosity Seyferts. To investigate this issue, we have obtained cold molecular gas observations, traced by the CO(2-1) transition, using the NOEMA interferometer of five nearby (distances between 19 and 58 Mpc) Seyfert galaxies. The resolution of approx. 0.3-0.8 arcsec (approx. 30-100 pc) and field of view of NOEMA allowed us to study the CO(2-1) morphology and kinematics in the nuclear regions (approx. 100 pc) and up to radial distances of approx. 900 pc. We have detected CO(2-1) emission in all five galaxies with disky or circumnuclear ring like morphologies. We derived cold molecular gas masses on nuclear (approx. 100 pc) and circumnuclear (approx. 650 pc) scales in the range from $10^6$ to $10^7$M$_{\odot}$ and from $10^7$ to $10^8$ $M_{\odot}$, respectively. In all of our galaxies the bulk of this gas is rotating in the plane of the galaxy. However, non-circular motions are also present. In NGC 4253, NGC 4388 and NGC 7465, we can ascribe the streaming motions to the presence of a large-scale bar. In Mrk 1066 and NGC 4388, the non-circular motions in the nuclear regions are explained as outflowing material due to the interaction of the AGN wind with molecular gas in the galaxy disk. We conclude that for an unambiguous and precise interpretation of the kinematics of the cold molecular gas we need a detailed knowledge of the host galaxy (i.e., presence of bars, interactions, etc) as well as of the ionized gas kinematics and the ionization cone geometry., Accepted for publication in A&A, 23 pages, 17 figures

We present infrared K'- (2.1 micron) and L'-band (3.8 micron) high-spatial-resolution (2 micron are used that should be sensitive to buried AGNs due to small dust extinction effects. For ULIRGs with resolved multiple nuclear K'-band emission, we estimate the activation of supermassive black holes (SMBHs) at individual galaxy nuclei in the form of AGN luminosity normalized by SMBH mass inferred from host galaxy stellar luminosity. We confirm a trend that more massive SMBHs in K'-band brighter primary galaxy nuclei are generally more active, with higher SMBH-mass-normalized AGN luminosity than less massive SMBHs in K'-band fainter secondary galaxy nuclei, as predicted by numerical simulations of gas-rich major galaxy mergers. In two sources, the presence of even infrared elusive extremely deeply buried AGNs is indicated by comparisons with available (sub)millimeter data. Non-synchronous SMBH activation (i.e., less activation of less massive SMBHs) and the possible presence of such infrared elusive AGNs may be responsible for the small fraction of infrared-detected luminous dual AGNs in nearby merging ULIRGs., 29 pages, 6 Figures (preprint style). Accepted for publication in ApJ. Resolutions of some figures are degraded

Full author list: John D. Silverman, Shenli Tang, Khee-Gan Lee, Tilman Hartwig, Andy Goulding, Michael A. Strauss, Malte Schramm, Xuheng Ding, Rogemar A. Riffel, Seiji Fujimoto, Chiaki Hikage, Masatoshi Imanishi, Kazushi Iwasawa, Knud Jahnke, Issha Kayo, Nobunari Kashikawa, Toshihiro Kawaguchi, Kotaro Kohno, Wentao Luo, Yoshiki Matsuoka, Yuichi Matsuda, Tohru Nagao, Masamune Oguri, Yoshiaki Ono, Masafusa Onoue, Masami Ouchi, Kazuhiro Shimasaku, Hyewon Suh, Nao Suzuki, Yoshiaki Taniguchi, Yoshiki Toba, Yoshihiro Ueda, and Naoki Yasuda, The unique combination of superb spatial resolution, wide-area coverage, and deep depth of the optical imaging from the Hyper Suprime-Cam (HSC) Subaru Strategic Program is utilized to search for dual quasar candidates. Using an automated image analysis routine on 34,476 known Sloan Digital Sky Survey quasars, we identify those with two (or more) distinct optical point sources in HSC images covering 796 deg2. We find 421 candidates out to a redshift of 4.5 of which one hundred or so are more likely after filtering out contaminating stars. Angular separations of 0.″6-4.″0 correspond to projected separations of 3-30 kpc, a range relatively unexplored for population studies of luminous dual quasars. Using Keck I/Low Resolution Imaging Spectrometer and Gemini-N/Near-Infrared Integral Field Spectrometer, we spectroscopically confirm three dual quasar systems at z < 1, two of which are previously unknown out of eight observed, based on the presence of characteristic broad emission lines in each component, while highlighting that the continuum of one object in one of the pairs is reddened. In all cases, the [O iii]λ5007 emission lines have mild velocity offsets, thus the joint [O iii] line profile is not double-peaked. We find a dual quasar fraction of 0.26 ± 0.18% and no evidence for evolution. A comparison with the Horizon-AGN simulation seems to support the case of no evolution in the dual quasar fraction when broadly matching the quasar selection. These results may indicate a scenario in which the frequency of the simultaneous triggering of luminous quasars is not as sensitive as expected to the cosmic evolution of the merger rate or gas content of galaxies., R. Morgan thanks the LSSTC Data Science Fellowship Program, which is funded by LSSTC, NSF Cybertraining Grant #1829740, the Brinson Foundation, and the Moore Foundation; his participation in the program has benefited this work. F. O.E. acknowledges support from the FONDECYT grant No. 1201223. This material is based upon work supported by the National Science Foundation Graduate Research Fellowship Program under grant No. 1744555. Any opinions, findings, and conclusions or recommendations expressed in this material are those of the author(s) and do not necessarily reflect the views of the National Science Foundation. Funding for the DES Projects has been provided by the U.S. Department of Energy, the U.S. National Science Foundation, the Ministry of Science and Education of Spain, the Science and Technology Facilities Council of the United Kingdom, the Higher Education Funding Council for England, the National Center for Supercomputing Applications at the University of Illinois at Urbana-Champaign, the Kavli Institute of Cosmological Physics at the University of Chicago, the Center for Cosmology and AstroParticle Physics at The Ohio State University, the Mitchell Institute for Fundamental Physics and Astronomy at Texas A&M University, Financiadora de Estudos e Projetos, Fundação Carlos Chagas Filho de Amparo à Pesquisa do Estado do Rio de Janeiro, Conselho Nacional de Desenvolvimento Científico e Tecnológico and the Ministério da Ciência, Tecnologia e Inovação, the Deutsche Forschungsgemeinschaft, and the Collaborating Institutions in the Dark Energy Survey. The Collaborating Institutions are Argonne National Laboratory, the University of California at Santa Cruz, the University of Cambridge, Centro de Investigaciones Energéticas, Medioambientales y Tecnológicas-Madrid, the University of Chicago, University College London, the DES-Brazil Consortium, the University of Edinburgh, the Eidgenössische Technische Hochschule (ETH) Zürich, Fermi National Accelerator Laboratory, the University of Illinois at UrbanaChampaign, the Institut de Ciències de l’Espai (IEEC/CSIC), the Institut de Física d’Altes Energies, Lawrence Berkeley National Laboratory, the Ludwig-Maximilians Universität München and the associated Excellence Cluster Universe, the University of Michigan, NFS’s NOIRLab, the University of Nottingham, The Ohio State University, the University of Pennsylvania, the University of Portsmouth, SLAC National Accelerator Laboratory, Stanford University, the University of Sussex, Texas A&M University, and the OzDES Membership Consortium. The UCSC team is supported in part by NASA grant NNG17PX03C, NSF grant AST-1815935, the Gordon & Betty Moore Foundation, the Heising-Simons Foundation, and by fellowships from the David and Lucile Packard Foundation to R.J.F. Based in part on observations at Cerro Tololo InterAmerican Observatory at NSFʼs NOIRLab (NOIRLab Prop. ID 2012B-0001; PI: J. Frieman), which is managed by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation. Based on observations obtained at the Southern Astrophysical Research (SOAR) telescope, which is a joint project of the Ministério da Ciência, Tecnologia, Inovação e Comunicações (MCTIC) da República Federativa do Brasil, the U.S. NSF’s NOIRLab, the University of North Carolina at Chapel Hill (UNC), and Michigan State University (MSU). The DES data management system is supported by the National Science Foundation under grant Nos. AST-1138766 and AST-1536171. The DES participants from Spanish institutions are partially supported by MICINN under grants ESP2017-89838, PGC2018-094773, PGC2018-102021, SEV2016-0588, SEV-2016-0597, and MDM-2015-0509, some of which include ERDF funds from the European Union. IFAE is partially funded by the CERCA program of the Generalitat de Catalunya. Research leading to these results has received funding from the European Research Council under the European Union’s Seventh Framework Program (FP7/2007- 2013), including ERC grant agreements 240672, 291329, and 306478. We acknowledge support from the Brazilian Instituto Nacional de Ciência e Tecnologia (INCT) e-Universe (CNPq grant 465376/2014-2). This paper has gone through internal review by the DES collaboration. This manuscript has been authored by the Fermi Research Alliance, LLC under Contract No. DE-AC02- 07CH11359 with the U.S. Department of Energy, Office of Science, Office of High Energy Physics.

We present ALMA ~0.02"-resolution observations of the nucleus of the nearby (~14 Mpc) type-2 AGN NGC 1068 at HCN/HCO+/HNC J=3-2 lines, as well as at their 13C isotopologue and vibrationally excited lines, to scrutinize the morphological/dynamical/chemical/physical properties of dense molecular gas in the putative dusty molecular torus around a mass-accreting supermassive black hole. We confirm almost east-west-oriented dense molecular gas emission both morphologically and dynamically, which we regard as coming from the torus. Bright emission is compact (2 pc) part, in both the eastern and western torus. We prefer a scenario of counter-rotating dense molecular gas with innermost almost-Keplerian-rotation and outer slowly rotating (far below Keplerian) components. Our high-spatial-resolution dense molecular line data reveal that torus properties of NGC 1068 are much more complicated than the simple axi-symmetrically rotating torus picture in the classical AGN unification paradigm., 45 pages, 20 Figures. Accepted for publication in ApJ

We present a large sample of 2.5-38 $\mu m$ galaxy spectra drawn from a cross-archival comparison in the AKARI-Spitzer Extragalactic Spectral Survey (ASESS), and investigate a subset of 113 star-forming galaxies with prominent polycyclic aromatic hydrocarbon (PAH) emission spanning a wide range of star formation properties. With AKARI's extended 2.5-5 $\mu m$ wavelength coverage, we self-consistently model for the first time all PAH emission bands using a modified version of PAHFIT. We find $L_\mathrm{PAH\,3.3}/L_\mathrm{IR}$ $\sim$ 0.1% and the 3.3 $\mu m$ PAH feature contributes $\sim$1.5-3% to the total PAH power -- somewhat less than earlier dust models have assumed. We establish a calibration between 3.3 $\mu m$ PAH emission and star formation rate, but also find regimes where it loses reliability, including at high luminosity and low metallicity. The 3.4 $\mu m$ aliphatic emission and a broad plateau feature centered at 3.47 $\mu m$ are also modeled. As the shortest wavelength PAH feature, 3.3 $\mu m$ is susceptible to attenuation, leading to a factor of $\sim$ 3 differences in the inferred star formation rate at high obscuration with different assumed attenuation geometries. Surprisingly, $L_\mathrm{PAH\, 3.3} / L_\mathrm{\Sigma\, PAH}$ shows no sign of decline at high luminosities, and the low metallicity dwarf galaxy II Zw 40 exhibits an unusually strong 3.3 $\mu m$ band; both results suggest either that the smallest PAHs are better able to survive under intense radiation fields than presumed, or that PAH emission is shifted to shorter wavelengths in intense and high energy radiation environments. A photometric surrogate for 3.3 $\mu m$ PAH luminosity using JWST/NIRCam is provided and found to be highly reliable at low redshift., Comment: 29 pages, 18 figures, accepted in ApJ

A prominent jet-driven outflow of CO(2-1) molecular gas is found along the kinematic minor axis of the Seyfert 2 galaxy ESO 420-G13, at a distance of $340-600\, \rm{pc}$ from the nucleus. The wind morphology resembles a characteristic funnel shape, formed by a highly collimated filamentary emission at the base, likely tracing the jet propagation through a tenuous medium, until a bifurcation point at $440\, \rm{pc}$ where the jet hits a dense molecular core and shatters, dispersing the molecular gas into several clumps and filaments within the expansion cone. We also trace the jet in ionised gas within the inner $\lesssim 340\, \rm{pc}$ using the [NeII]$_{\rm 12.8 \mu m}$ line emission, where the molecular gas follows a circular rotation pattern. The wind outflow carries a mass of $\sim 8 \times 10^6\, \rm{M_\odot}$ at an average wind projected speed of $\sim 160\, \rm{km\,s^{-1}}$, which implies a mass outflow rate of $\sim 14\, \rm{M_\odot\,yr^{-1}}$. Based on the structure of the outflow and the budget of energy and momentum, we discard radiation pressure from the active nucleus, star formation, and supernovae as possible launching mechanisms. ESO 420-G13 is the second case after NGC 1377 where the presence of a previously unknown jet is revealed due to its interaction with the interstellar medium, suggesting that unknown jets in feeble radio nuclei might be more common than expected. Two possible jet-cloud configurations are discussed to explain the presence of an outflow at such distance from the AGN. The outflowing gas will likely not escape, thus a delay in the star formation rather than quenching is expected from this interaction, while the feedback effect would be confined within the central few hundred parsecs of the galaxy., Comment: 12 pages, 8 figures. Accepted for publication in A&A

The obscuring structure surrounding active galactic nuclei (AGN) can be explained as a dust and gas flow cycle that fundamentally connects the AGN with their host galaxies. This structure is believed to be associated with dusty winds driven by radiation pressure. However, the role of magnetic fields, which are invoked in almost all models for accretion onto a supermassive black hole and outflows, is not thoroughly studied. Here we report the first detection of polarized thermal emission by means of magnetically aligned dust grains in the dusty torus of NGC 1068 using ALMA Cycle 4 polarimetric dust continuum observations ($0.07"$, $4.2$ pc; 348.5 GHz, $860$ $\mu$m). The polarized torus has an asymmetric variation across the equatorial axis with a peak polarization of $3.7\pm0.5$\% and position angle of $109\pm2^{\circ}$ (B-vector) at $\sim8$ pc east from the core. We compute synthetic polarimetric observations of magnetically aligned dust grains assuming a toroidal magnetic field and homogeneous grain alignment. We conclude that the measured 860 $\mu$m continuum polarization arises from magnetically aligned dust grains in an optically thin region of the torus. The asymmetric polarization across the equatorial axis of the torus arises from 1) an inhomogeneous optical depth, and 2) a variation of the velocity dispersion, i.e. variation of the magnetic field turbulence at sub-pc scales, from the eastern to the western region of the torus. These observations and modeling constrain the torus properties beyond spectral energy distribution results. This study strongly supports that magnetic fields up to a few pc contribute to the accretion flow onto the active nuclei., Comment: 19 pages, 11 figures (Accepted for Publication to ApJ)

We present ALMA [CII] line and far-infrared (FIR) continuum observations of three $z > 6$ low-luminosity quasars ($M_{\rm 1450} > -25$) discovered by our Subaru Hyper Suprime-Cam (HSC) survey. The [CII] line was detected in all three targets with luminosities of $(2.4 - 9.5) \times 10^8~L_\odot$, about one order of magnitude smaller than optically luminous ($M_{\rm 1450} \lesssim -25$) quasars. The FIR continuum luminosities range from $< 9 \times 10^{10}~L_\odot$ (3$��$ limit) to $\sim 2 \times 10^{12}~L_\odot$, indicating a wide range in star formation rates in these galaxies. Most of the HSC quasars studied thus far show [CII]/FIR luminosity ratios similar to local star-forming galaxies. Using the [CII]-based dynamical mass ($M_{\rm dyn}$) as a surrogate for bulge stellar mass ($M_{\rm bulge}$), we find that a significant fraction of low-luminosity quasars are located on or even below the local $M_{\rm BH} - M_{\rm bulge}$ relation, particularly at the massive end of the galaxy mass distribution. In contrast, previous studies of optically luminous quasars have found that black holes are overmassive relative to the local relation. Given the low luminosities of our targets, we are exploring the nature of the early co-evolution of supermassive black holes and their hosts in a less biased way. Almost all of the quasars presented in this work are growing their black hole mass at much higher pace at $z \sim 6$ than the parallel growth model, in which supermassive black holes and their hosts grow simultaneously to match the local $M_{\rm BH} - M_{\rm bulge}$ relation at all redshifts. As the low-luminosity quasars appear to realize the local co-evolutionary relation even at $z \sim 6$, they should have experienced vigorous starbursts prior to the currently observed quasar phase to catch up with the relation., 19 pages, 11 figures, 4 tables. Accepted for publication in Publications of the Astronomical Society of Japan (PASJ)

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.

Tang, S., et al., 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, 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), Japan Science and Technology Agency (JST), the Toray Science Foundation, NAOJ, Kavli IPMU, KEK, ASIAA, and Princeton University

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), 27 pages, 27 figures. Accepted for publication in A&A

We explore X-ray evidence for the presence of active galactic nuclei (AGNs) in the two most actively star-forming Green Pea galaxies (GPs), SDSS J0749+3337 and SDSS J0822+2241, which have star-formation rates (SFRs) of $123~M_\odot$ yr$^{-1}$ and $78~M_\odot$ yr$^{-1}$, respectively. The GPs have red mid-infrared (MIR) spectral energy distributions and higher 22 $\mu$m luminosities than expected from a proxy of the SFR (H$\alpha$ luminosity), consistent with hosting AGNs with 2-10 keV luminosities of $\sim10^{44}$ erg s$^{-1}$. We thus obtain and analyze the first hard ($>$ 10 keV) X-ray data observed with NuSTAR and archival XMM-Newton data below 10 keV. From the NuSTAR $\approx$20 ksec data, however, we find no significant hard X-ray emission. By contrast, soft X-ray emission with 0.5--8 keV luminosities of $\approx10^{42}$ erg s$^{-1}$ is significantly detected in both targets, which can be explained only by star formation (SF). A possible reason for the lack of clear evidence is that a putative AGN torus absorbs most of the X-ray emission. Applying a smooth-density AGN torus model, we determine minimum hydrogen column densities along the equatorial plane ($N_{\rm H}^{\rm eq}$) consistent with the non-detection. The results indicate $N_{\rm H}^{\rm eq} \gtrsim 2\times10^{24}$ cm$^{-2}$ for SDSS J0749+3337 and $N_{\rm H}^{\rm eq} \gtrsim 5\times10^{24}$ cm$^{-2}$ for SDSS J0822+2241. Therefore, the GPs may host such heavily obscured AGNs. Otherwise, no AGN exists and the MIR emission is ascribed to SF. Active SF in low-mass galaxies is indeed suggested to reproduce red MIR colors. This would imply that diagnostics based on MIR photometry data alone may misidentify such galaxies as AGNs., Comment: 12 pages, 3 tables, 5 figures, accepted for publication in ApJ

We present ALMA Band~7 observations at $850��$m of 20 luminous ($\log\, L_{\rm bol}>46.9$ [erg s$^{-1}$]) unobscured quasars at $z\sim2$. We detect continuum emission for 19/20 quasars. After subtracting an AGN contribution, we measure the total far-IR luminosity for 18 quasars, assuming a modified blackbody model, and attribute the emission as indicative of the star formation rate (SFR). Our sample can be characterized with a log-normal SFR distribution having a mean of 140 $M_\odot$ yr$^{-1}$ and a dispersion of 0.5 dex. Based on an inference of their stellar masses, the SFRs are similar, in both the mean and dispersion, with star-forming main-sequence galaxies at the equivalent epoch. Thus, there is no evidence for a systematic enhancement or suppression (i.e., regulation or quenching) of star formation in the hosts of the most luminous quasars at $z\sim2$. These results are consistent with the Magneticum cosmological simulation, while in disagreement with a widely recognized phenomenological model that predicts higher SFRs than observed here based on the high bolometric luminosities of this sample. Furthermore, there is only a weak relation between SFR and accretion rate onto their supermassive black holes both for average and individual measurements. We interpret these results as indicative of star formation and quasar accretion being fed from the available gas reservoir(s) in their host with a disconnect due to their different physical sizes, temporal scales, and means of gas processing., 20 pages, 12 figures, re-submitted to MNRAS

We report on the results of deep and wide-field (1.1 deg$^2$) narrow-band observations with Subaru/Hyper Suprime-Cam (HSC) of a field around a hyperluminous QSO (HLQSO), HS1549+1919, residing in a protocluster at $z=2.84$, to map the large-scale structure of Ly$\alpha$ emitters (LAEs). One HSC pointing enables us to detect 3490 LAEs and 76 extended Ly$\alpha$ blobs (LABs), probing diverse environments from voids to protoclusters. The HLQSO is found to be near the center of the protocluster, which corresponds to the intersection of $\sim$100 cMpc-scale structures of LAEs. LABs are basically distributed along the large-scale structure, with larger ones particularly clustered around the HLQSO, confirming a previously noted tendency of LABs to prefer denser environments. Moreover, the shapes of LABs near the HLQSO appear to be aligned with the large-scale structure. Finally, a deep Ly$\alpha$ image reveals a diffuse Ly$\alpha$ nebula along a filamentary structure with no luminous UV/sub-mm counterpart. We suggest that the diffuse nebula is due to a cold filament with high clumping factor illuminated by the QSO, with a required high clumpiness provided by unresolved residing halos of mass $\leq 10^{9-10}M_\odot$., Comment: 6 pages, 3 figures, submitted to PASJ

We present deep near-infrared spectroscopy of six quasars at 6.15.8 quasars with measured black hole masses. From single-epoch mass measurements based on MgII2798, we find a wide range in black hole masses, from M_BH=10^7.6 to 10^9.3 Msun. The Eddington ratios L_bol/L_Edd range from 0.16 to 1.1, but the majority of the HSC quasars are powered by M_BH=10^9 Msun supermassive black holes (SMBHs) accreting at sub-Eddington rates. The Eddington ratio distribution of the HSC quasars is inclined to lower accretion rates than those of Willott et al. (2010a), who measured the black hole masses for similarly faint z=6 quasars. This suggests that the global Eddington ratio distribution is wider than has previously been thought. The presence of M_BH=10^9 Msun SMBHs at z=6 cannot be explained with constant sub-Eddington accretion from stellar remnant seed black holes. Therefore, we may be witnessing the first buildup of the most massive black holes in the first billion years of the universe, the accretion activity of which is transforming from active growth to a quiescent phase. Measurements of a larger complete sample of z>6 low-luminosity quasars, as well as deeper observations with future facilities will enable us to better understand the early SMBH growth in the reionization epoch., Comment: 20 pages, 10 figures, accepted for publication in ApJ

We report a study of X-ray-irradiated gas in the central ~100 pc of the Circinus galaxy, hosting a Compton-thick active galactic nucleus (AGN), at 10-pc resolution using Chandra and ALMA. Based on ~200 ksec Chandra/ACIS-S data, we created an image of the Fe Kalpha line at 6.4 keV, tracing X-ray-irradiated dense gas. The ALMA data in Bands 6 (~270 GHz) and 7 (~350 GHz) cover five molecular lines: CO(J=3--2), HCN(J=3--2), HCN(J=4--3), HCO^+(J=3--2), and HCO^+(J=4--3). The detailed spatial distribution of dense molecular gas was revealed, and compared to the iron line image. The molecular gas emission appeared faint in regions with bright iron emission. Motivated by this, we quantitatively discuss the possibility that the molecular gas is efficiently dissociated by AGN X-ray irradiation (i.e., creating an X-ray-dominated region). Based on a non-local thermodynamic equilibrium model, we constrained the molecular gas densities and found that they are as low as interpreted by X-ray dissociation. Furthermore, judging from inactive star formation (SF) reported in the literature, we suggest that the X-ray emission has potential to suppress SF, particularly in the proximity of the AGN., 20 pages, 47 figures, accepted for publication in PASJ

Understanding the nuclear growth and feedback processes in galaxies requires investigating their often obscured central regions. One way to do this is to use (sub)millimeter line emission from vibrationally excited HCN (HCN-vib). It has been suggested that the most intense HCN-vib emission from a galaxy is connected to a phase of nuclear growth that occurs before the nuclear feedback processes have been fully developed. We aim to investigate if there is a connection between the presence of strong HCN-vib emission and the development of feedback in (U)LIRGs. We collected literature and archival data to compare the luminosities of rotational lines of HCN-vib, normalized to the total infrared luminosity, to the median velocities of 119 {\mu}m OH absorption lines, potentially indicating outflows, in a total of 17 (U)LIRGs. The most HCN-vib luminous systems all lack signatures of significant molecular outflows in the far-infrared OH absorption lines. However, at least some of the systems with bright HCN-vib emission do have fast and collimated outflows that can be seen in spectral lines at longer wavelengths. We conclude that the galaxy nuclei with the highest L(HCN-vib)/L(IR) do not drive wide-angle outflows detectable using the median velocities of far-infrared OH absorption lines. It is possible that this is due to an orientation effect where sources which are oriented in such a way that their outflows are not along our line of sight also radiate a smaller proportion of their infrared luminosity in our direction. It could also be that massive wide-angle outflows destroy the deeply embedded regions responsible for bright HCN-vib emission, so that the two phenomena cannot coexist. This would strengthen the idea that vibrationally excited HCN traces a heavily obscured stage of evolution before nuclear feedback mechanisms are fully developed., Comment: 10 pages, 4 figures, Accepted for publication in Astronomy and Astrophysics

We propose new diagnostics that utilize the [O IV] 25.89 $\mu$m and nuclear (subarcsecond scale) 12 $\mu$m luminosity ratio for identifying whether an AGN is deeply `buried' in their surrounding material. Utilizing a sample of 16 absorbed AGNs at redshifts lower than 0.03 in the Swift/BAT catalog observed with Suzaku, we find that AGNs with small scattering fractions ($, Comment: 8 pages, 4 figures, (accepted in the Astrophysical Journal, ApJ)

Hydrogen recombination lines at the submillimetre band (submm-RLs) can serve as probes of ionized gas without dust extinction. One therefore expects to probe the broad line region (BLR) of an obscured (type 2) active galactic nucleus (AGN) with those lines. However, admitting the large uncertainty in the continuum level, here we report on the non-detection of both broad and narrow H26$\alpha$ emission line (rest frequency = 353.62 GHz) towards the prototypical type 2 AGN of NGC 1068 with the Atacama Large Millimeter/submillimeter Array (ALMA). We also investigate the nature of BLR clouds that can emit submm-RLs with model calculations. As a result, we suggest that clouds with an electron density ($N_e$) of $\sim$ 10$^9$ cm$^{-3}$ can mainly contribute to broad submm-RLs in terms of the line flux. On the other hand, line flux from other density clouds would be insignificant considering their too large or too small line optical depths. However, even for the case of $N_e \sim 10^9$ cm$^{-3}$ clouds, we also suggest that the expected line flux is extremely low, which is impractical to detect even with ALMA., Comment: 7 pages, 4 figures, and 1 table; To be published in MNRAS; This is the version 2 of the manuscript. One significant computational error that influences one of our previous conclusions (1602.07316) was modified accordingly. Our sincere apology for this issue

Angel Castro,Takamitsu Miyaji, Mai Shirahata, Kohei Ichikawa, Shinki Oyabu, David M. Clark, Masatoshi Imanishi, Takao Nakagawa, Yoshihiro Ueda 1Instituto de Astronomia, Universidad Nacional Autonoma de Mexico (UNAM), Ensenada, Baja California, Mexico (mailing address: PO Box 439027, San Diego, CA 92143-9027, USA) 2University of California, San Diego, Center for Astrophysics and Space Sciences, 9500 Gilman Drive, La Jolla, CA 92093-0424, USA 3National Astronomical Observatory of Japan (NAOJ), 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan 4Institute of Space and Astronautical Science (ISAS), Japan Aerospace Exploration Agency (JAXA), 3-1-1 Yoshino-dai, Chuo-ku, Sagamihara 252-5210 5Department of Astronomy, Kyoto University, Kitashirakawa-Oiwake-cho, Sakyo-ku, Kyoto 606-8502 6Graduate School of Science, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8602, Japan 7Subaru Telescope, 650 North Aohoku Place, Hilo, Hawaii, 96720, USA acastro@astrosen.unam.mx

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

We used the Atacama Large Millimeter/submillimeter Array (ALMA) to map $^{12}$CO($J$ = 1-0), $^{12}$CO($J$ = 2-1), $^{12}$CO($J$ = 3-2), $^{13}$CO($J$ = 2-1), and [CI]($^3P_1$-$^3P_0$) emission lines around the type 1 active galactic nucleus (AGN) of NGC 7469 ($z = 0.0164$) at $\sim 100$ pc resolutions. The CO lines are bright both in the circumnuclear disk (central $\sim 300$ pc) and the surrounding starburst (SB) ring ($\sim 1$ kpc diameter), with two bright peaks on either side of the AGN. By contrast, the [CI]($^3P_1$-$^3P_0$) line is strongly peaked on the AGN. Consequently, the brightness temperature ratio of [CI]($^3P_1$-$^3P_0$) to $^{13}$CO(2-1) is $\sim 20$ at the AGN, as compared to $\sim 2$ in the SB ring. Our local thermodynamic equilibrium (LTE) and non-LTE models indicate that the enhanced line ratios (or CI enhancement) are due to an elevated C$^0$/CO abundance ratio ($\sim 3-10$) and temperature ($\sim 100-500$ K) around the AGN as compared to the SB ring (abundance ratio $\sim 1$, temperature $\lesssim 100$ K), which accords with the picture of the X-ray-dominated Region (XDR). Based on dynamical modelings, we also provide CO(1-0)-to- and [CI]($^3P_1$-$^3P_0$)-to-molecular mass conversion factors at the central $\sim 100$ pc of this AGN as $\alpha_{\rm CO} = 4.1$ and $\alpha_{\rm CI} = 4.4~M_\odot$ (K km s$^{-1}$ pc$^2$)$^{-1}$, respectively. Our results suggest that the CI enhancement is potentially a good marker of AGNs that could be used in a new submillimeter diagnostic method toward dusty environments., Comment: 25 pages, 12 figures, 8 tables. Accepted for publication in ApJ. Minor updates in a replacement

We have analyzed the broadband X-ray spectra of active galactic nuclei (AGNs) in two "non-merging" luminous infrared galaxies (LIRGs) UGC 2608 and NGC 5135, utilizing the data of NuSTAR, Suzaku, XMM-Newton, and Chandra. Applying the X-ray clumpy-torus model (XCLUMPY: Tanimoto et al. 2019), we find that both sources have similar spectra characterized by Compton-thick (CT) absorption ($N_{\rm H} \sim$ 5-7 $\times$ $10^{24}$ cm$^{-2}$) and small torus angular width ($\sigma$ $, Comment: 15 pages, 12 figures, and 3 tables, accepted for publication in ApJ

We report spatial distributions of the Fe-K$\alpha$ line at 6.4 keV and the CO($J$ = 2--1) line at 230.538 GHz in NGC 2110, which are respectively revealed by $Chandra$ and ALMA at $\approx$ 0.5 arcsec. A $Chandra$ 6.2--6.5 keV-to-3.0--6.0 keV image suggests that the Fe-K$\alpha$ emission extends preferentially in a northwest-to-southeast direction out to $\sim$ 3 arcsec, or 500 pc, on each side. Spatially-resolved spectral analyses support this by finding significant Fe-K$\alpha$ emission lines only in northwest and southeast regions. Moreover, their equivalent widths are found $\sim$ 1.5 keV, indicative for the fluorescence by nuclear X-ray irradiation as the physical origin. By contrast, CO($J$ = 2--1) emission is weak therein. For quantitative discussion, we derive ionization parameters by following an X-ray dominated region (XDR) model. We then find them high enough to interpret the weakness as the result of X-ray dissociation of CO and/or H$_2$. Another possibility also remains that CO molecules follow a super-thermal distribution, resulting in brighter emission in higher-$J$ lines. Further follow-up observations are encouraged to draw a conclusion on what predominantly changes the inter-stellar matter properties, and whether the X-ray irradiation eventually affects the surrounding star formation as an AGN feedback., Comment: 15 pages, 11 figures, 3 tables. Accepted for publication in ApJ

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

With the imminent launch of the JWST, the field of thermal-infrared (TIR) astronomy will enjoy a revolution. It is easy to imagine that all areas of infrared (IR) astronomy will be greatly advanced, but perhaps impossible to conceive of the new vistas that will be opened. To allow both follow-up JWST observations and a continuance of work started on the ground-based 8m’s, we continue to plan the science cases and instrument design for a TIR imager and spectrometer for early operation on the TMT. We present the current status of our science cases and the instrumentation plans, harnessing expertise across the TMT partnership. This instrument will be proposed by the MICHI team as a second-generation instrument in any upcoming calls for proposals., Ground-based and Airborne Instrumentation for Astronomy VII, June 10-15, 2018, Austin, USA, Series: Proceedings of SPIE; no. 10702

We present HCN J=4-3 and HCO^+ J=4-3 maps of six nearby star-forming galaxies, NGC 253, NGC 1068, IC 342, M82, M83, and NGC 6946, obtained with the James Clerk Maxwell Telescope as part of the MALATANG survey. All galaxies were mapped in the central 2 arcmin $\times$ 2 arcmin region at 14 arcsec (FWHM) resolution (corresponding to linear scales of ~ 0.2-1.0 kpc). The L_IR-L'_dense relation, where the dense gas is traced by the HCN J=4-3 and the HCO^+ J=4-3 emission, measured in our sample of spatially-resolved galaxies is found to follow the linear correlation established globally in galaxies within the scatter. We find that the luminosity ratio, L_IR/L'_dense, shows systematic variations with L_IR within individual spatially resolved galaxies, whereas the galaxy-integrated ratios vary little. A rising trend is also found between L_IR/L'_dense ratio and the warm-dust temperature gauged by the 70 \mu m/100 \mu m flux ratio. We find the luminosity ratios of IR/HCN(4-3) and IR/HCO^+(4-3), which can be taken as a proxy for the efficiency of star formation in the dense molecular gas (SFE_dense), appears to be nearly independent of the dense-gas fraction (f_dense) for our sample of galaxies. The SFE of the total molecular gas (SFE_mol) is found to increase substantially with f_dense when combining our data with that on local (ultra)luminous infrared galaxies and high-z quasars. The mean L'_HCN(4-3)/L'_HCO^+(4-3) line ratio measured for the six targeted galaxies is 0.9+/-0.6. No significant correlation is found for the L'_HCN(4-3)/L'_HCO^+(4-3) ratio with the SFR as traced by L_IR, nor with the warm-dust temperature, for the different populations of galaxies., Comment: 29 pages, 12 figures, accepted for publication in ApJ

We present our ALMA Cycle 4 measurements of the [CII] emission line and the underlying far-infrared (FIR) continuum emission from four optically low-luminosity ($M_{\rm 1450} > -25$) quasars at $z \gtrsim 6$ discovered by the Subaru Hyper Suprime Cam (HSC) survey. The [CII] line and FIR continuum luminosities lie in the ranges $L_{\rm [CII]} = (3.8-10.2) \times 10^8~L_\odot$ and $L_{\rm FIR} = (1.2-2.0) \times 10^{11}~L_\odot$, which are at least one order of magnitude smaller than those of optically-luminous quasars at $z \gtrsim 6$. We estimate the star formation rates (SFR) of our targets as $\simeq 23-40~M_\odot ~{\rm yr}^{-1}$. Their line and continuum-emitting regions are marginally resolved, and found to be comparable in size to those of optically luminous quasars, indicating that their SFR or likely gas mass surface densities (key controlling parameter of mass accretion) are accordingly different. The $L_{\rm [CII]}/L_{\rm FIR}$ ratios of the hosts, $\simeq (2.2-8.7) \times 10^{-3}$, are fully consistent with local star-forming galaxies. Using the [CII] dynamics, we derived their dynamical masses within a radius of 1.5-2.5 kpc as $\simeq (1.4-8.2) \times 10^{10}~M_\odot$. By interpreting these masses as stellar ones, we suggest that these faint quasar hosts are on or even below the star-forming main sequence at $z \sim 6$, i.e., they appear to be transforming into quiescent galaxies. This is in contrast to the optically luminous quasars at those redshifts, which show starburst-like properties. Finally, we find that the ratios of black hole mass to host galaxy dynamical mass of the most of low-luminosity quasars including the HSC ones are consistent with the local value. The mass ratios of the HSC quasars can be reproduced by a semi-analytical model that assumes merger-induced black hole-host galaxy evolution., 20 pages, 9 figures, 4 tables. Accepted for publication in PASJ

We present the results of our ALMA Cycle 4 high-spatial-resolution (0.04-0.07") observations, at HCN J=3-2 and HCO+ J=3-2 lines, of the nucleus of NGC 1068, the nearby prototypical type 2 active galactic nucleus (AGN). Our previous ALMA observations identified the compact emission of these lines at the putative location of the torus around a mass-accreting supermassive black hole. We now report that we have detected the rotation of this compact emission, with the eastern and western sides being redshifted and blueshifted, respectively. Unlike the previously reported CO J=6-5 emission, both the morphological and dynamical alignments of the HCN J=3-2 and HCO+ J=3-2 emission are roughly aligned along the east-west direction (i.e., the expected torus direction), suggesting that these molecular lines are better probes of a rotating dense molecular gas component in the torus. The western part of the torus exhibits larger velocity dispersion and stronger emission in the HCN J=3-2 and HCO+ J=3-2 lines than the eastern part, revealing a highly inhomogeneous molecular torus. The dense molecular gas in the torus and that of the host galaxy at 0.5-2.0" from the AGN along the torus direction are found to be counter-rotating, suggesting an external process happened in the past at the NGC 1068 nucleus., Comment: 7 pages, 3 figures, accepted for publication in ApJL. Figures 1 and 2 are degraded. High resolution version is available at http://www2.nao.ac.jp/~masaimanishi/N1068/N1068bMemuH.pdf

Deep MIR surveys have revealed numerous strongly star-forming galaxies at redshift z~0.3 mJy) around the North Ecliptic Pole with the Infrared Camera (IRC) onboard AKARI. After removing 11 AGN candidates by using the IRC photometry, we identified 48 PAH galaxies with PAH 6.2, 7.7, and 8.6um features at z0.35, whose SEDs and luminosities are typical of luminous infrared galaxies. They show particularly larger PAH enhancement at high luminosity, implying that they are vigorous star-forming galaxies with elevated specific SFR. Our composite starburst model that combines a very young and optically very thick starburst with a very old population can successfully reproduce most of their SED characteristics, although we could not confirm this optically thick component from our spectral analysis., 61 pages, 24 figures, accepted for publication in Astronomy and Astrophysics

We present a near-infrared (NIR) spectrum of WISE J1042+1641, an extremely red dust-obscured galaxy (DOG), which has been observed with the LIRIS on the 4.2m William Hershel Telescope. This object was selected as a hyper-luminous DOG candidate at z ~ 2 by combining the optical and IR photometric data based on the SDSS and WISE, although its redshift had not yet been confirmed. Based on the LIRIS observation, we confirmed its redshift of 2.521 and total IR luminosity of log(L_IR/L_sun) = 14.57, which satisfies the criterion for an extremely luminous IR galaxy (ELIRG). Moreover, we indicate that this object seems to have an extremely massive black hole with M_BH = 10^10.92 M_sun based on the broad Halpha line: the host stellar mass is derived as M_star = 10^13.55 M_sun by a fit of the spectral energy distribution. Very recently, it has been reported that this object is an anomalous gravitationally lensed quasar based on near-IR high-resolution imaging data obtained with the Hubble Space Telescope. Its magnification factor has also been estimated with some uncertainty (i.e., mu = 53-122). We investigate the ratio of the black hole to galaxy mass, which is less strongly affected by a lensing magnification factor, instead of the absolute values of the luminosities and masses. We find that the M_BH/M_star ratio (i.e., 0.0140-0.0204) is significantly higher than the local relation, following a sequence of unobscured quasars instead of obscured objects (e.g., submillimeter galaxies) at the same redshift. Moreover, the LIRIS spectrum shows strongly blueshifted oxygen lines with an outflowing velocity of ~ 1100 km/s, and our Swift X-ray observation also supports that this source is an absorbed AGN with an intrinsic column density of N_H = 4.9 x 10^23 cm^-2. These results imply that WISE J1042+1641 is in a blow-out phase at the end of the buried rapid black hole growth., Comment: 6 pages, 4 figures, accepted for publication in A&A

We present the luminosity function of z=4 quasars based on the Hyper Suprime-Cam Subaru Strategic Program Wide layer imaging data in the g, r, i, z, and y bands covering 339.8 deg^2. From stellar objects, 1666 z~4 quasar candidates are selected by the g-dropout selection down to i=24.0 mag. Their photometric redshifts cover the redshift range between 3.6 and 4.3 with an average of 3.9. In combination with the quasar sample from the Sloan Digital Sky Survey in the same redshift range, the quasar luminosity function covering the wide luminosity range of M1450=-22 to -29 mag is constructed. It is well described by a double power-law model with a knee at M1450=-25.36+-0.13 mag and a flat faint-end slope with a power-law index of -1.30+-0.05. The knee and faint-end slope show no clear evidence of redshift evolution from those at z~2. The flat slope implies that the UV luminosity density of the quasar population is dominated by the quasars around the knee, and does not support the steeper faint-end slope at higher redshifts reported at z>5. If we convert the M1450 luminosity function to the hard X-ray 2-10keV luminosity function using the relation between UV and X-ray luminosity of quasars and its scatter, the number density of UV-selected quasars matches well with that of the X-ray-selected AGNs above the knee of the luminosity function. Below the knee, the UV-selected quasars show a deficiency compared to the hard X-ray luminosity function. The deficiency can be explained by the lack of obscured AGNs among the UV-selected quasars., 24 pages, 21 figures, submitted to PASJ

We used the Atacama Large Millimeter/submillimeter Array (ALMA) to map the CO(3-2) and the underlying continuum emissions around the type 1 low-luminosity active galactic nucleus (LLAGN; bolometric luminosity $\lesssim 10^{42}$ erg~s$^{-1}$) of NGC 1097 at $\sim 10$ pc resolution. These observations revealed a detailed cold gas distribution within a $\sim 100$ pc of this LLAGN. In contrast to the luminous Seyfert galaxy NGC 1068, where a $\sim 7$ pc cold molecular torus was recently revealed, a distinctively dense and compact torus is missing in our CO(3-2) integrated intensity map of NGC 1097. Based on the CO(3-2) flux, the gas mass of the torus of NGC 1097 would be a factor of $\gtrsim 2-3$ less than that found for NGC 1068 by using the same CO-to-H$_2$ conversion factor, which implies less active nuclear star formation and/or inflows in NGC 1097. Our dynamical modeling of the CO(3-2) velocity field implies that the cold molecular gas is concentrated in a thin layer as compared to the hot gas traced by the 2.12 $\mu$m H$_2$ emission in and around the torus. Furthermore, we suggest that NGC 1097 hosts a geometrically thinner torus than NGC 1068. Although the physical origin of the torus thickness remains unclear, our observations support a theoretical prediction that geometrically thick tori with high opacity will become deficient as AGNs evolve from luminous Seyferts to LLAGNs., Comment: 9 pages, 5 figures. Accepted for publication in ApJL

We present the cross-correlation between 151 luminous quasars ($M_{ \mathrm{UV}} < -26$) and 179 protocluster candidates at $z \sim 3.8$, extracted from the Wide imaging survey ($ \sim 121~ $deg$^2$) performed with a part of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). We find that only two out of 151 quasars reside in regions that are more overdense compared to the average field at $ > 4 ��$. The distributions of the distance between quasars and the nearest protoclusters and the significance of the overdensity at the position of quasars are statistically identical to those found for $g$-dropout galaxies, suggesting that quasars tend to reside in almost the same environment as star-forming galaxies at this redshift. Using stacking analysis, we find that the average density of $g$-dropout galaxies around quasars is slightly higher than that around $g$-dropout galaxies on $1.0 - 2.5$ pMpc scales, while at $ < 0.5$ pMpc that around quasars tends to be lower. We also find that quasars with higher UV-luminosity or with more massive black holes tend to avoid the most overdense regions, and that the quasar near zone sizes are anti-correlated with overdensity. These findings are consistent with a scenario in which the luminous quasar at $z \sim4 $ resides in structures that are less massive than those expected for the progenitors of today's rich clusters of galaxies, and possibly that luminous quasars may be suppressing star formation in their close vicinity., 15 pages, 11 figures, accepted for publication in PASJ special issue

We investigate the galaxy overdensity around proto-cluster scale quasar pairs at high (z>3) and low (z~1) redshift based on the unprecedentedly wide and deep optical survey of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). Using the first-year survey data covering effectively ~121 deg^2 with the 5sigma depth of i~26.4 and the SDSS DR12Q catalog, we find two luminous pairs at z~3.3 and 3.6 which reside in >5sigma overdense regions of g-dropout galaxies at i4sigma) regions. Thanks to the large sample size, we find a statistical evidence that this excess is unique to the pair environments when compared to single quasar and randomly selected galaxy environments at the same redshift range. Moreover, there are nine small-scale (R_perp, Comment: 22 pages, 9 figures, 7 tables, accepted for publication in PASJ special issue

Hyper Suprime-Cam (HSC) is a wide-field imaging camera on the prime focus of the 8.2m Subaru telescope on the summit of Maunakea in Hawaii. A team of scientists from Japan, Taiwan and Princeton University is using HSC to carry out a 300-night multi-band imaging survey of the high-latitude sky. The survey includes three layers: the Wide layer will cover 1400 deg$^2$ in five broad bands ($grizy$), with a $5\,\sigma$ point-source depth of $r \approx 26$. The Deep layer covers a total of 26~deg$^2$ in four fields, going roughly a magnitude fainter, while the UltraDeep layer goes almost a magnitude fainter still in two pointings of HSC (a total of 3.5 deg$^2$). Here we describe the instrument, the science goals of the survey, and the survey strategy and data processing. This paper serves as an introduction to a special issue of the Publications of the Astronomical Society of Japan, which includes a large number of technical and scientific papers describing results from the early phases of this survey., Comment: 14 pages, 7 figures, 5 tables. Corrected for a typo in the coordinates of HSC-Wide spring equatorial field in Table 5

We present spectroscopic identification of 32 new quasars and luminous galaxies discovered at 5.7 < z < 6.8. This is the second in a series of papers presenting the results of the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The photometric candidates were selected by a Bayesian probabilistic algorithm, and then observed with spectrographs on the Gran Telescopio Canarias and the Subaru Telescope. Combined with the sample presented in the previous paper, we have now identified 64 HSC sources over about 430 deg2, which include 33 high-z quasars, 14 high-z luminous galaxies, 2 [O III] emitters at z ~ 0.8, and 15 Galactic brown dwarfs. The new quasars have considerably lower luminosity (M1450 ~ -25 to -22 mag) than most of the previously known high-z quasars. Several of these quasars have luminous (> 10^(43) erg/s) and narrow (< 500 km/s) Ly alpha lines, and also a possible mini broad absorption line system of N V 1240 in the composite spectrum, which clearly separate them from typical quasars. On the other hand, the high-z galaxies have extremely high luminosity (M1450 ~ -24 to -22 mag) compared to other galaxies found at similar redshift. With the discovery of these new classes of objects, we are opening up new parameter spaces in the high-z Universe. Further survey observations and follow-up studies of the identified objects, including the construction of the quasar luminosity function at z ~ 6, are ongoing., Comment: Accepted for publication in PASJ (HSC special issue); minor revisions made since ver. 1 and 2

The Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) is a three-layered imaging survey aimed at addressing some of the most outstanding questions in astronomy today, including the nature of dark matter and dark energy. The survey has been awarded 300 nights of observing time at the Subaru Telescope and it started in March 2014. This paper presents the first public data release of HSC-SSP. This release includes data taken in the first 1.7 years of observations (61.5 nights) and each of the Wide, Deep, and UltraDeep layers covers about 108, 26, and 4 square degrees down to depths of i~26.4, ~26.5, and ~27.0 mag, respectively (5sigma for point sources). All the layers are observed in five broad bands (grizy), and the Deep and UltraDeep layers are observed in narrow bands as well. We achieve an impressive image quality of 0.6 arcsec in the i-band in the Wide layer. We show that we achieve 1-2 per cent PSF photometry (rms) both internally and externally (against Pan-STARRS1), and ~10 mas and 40 mas internal and external astrometric accuracy, respectively. Both the calibrated images and catalogs are made available to the community through dedicated user interfaces and database servers. In addition to the pipeline products, we also provide value-added products such as photometric redshifts and a collection of public spectroscopic redshifts. Detailed descriptions of all the data can be found online. The data release website is https://hsc-release.mtk.nao.ac.jp/., 34 pages, 20 figures, 7 tables, moderate revision, accepted for publication in PASJ

We report the discovery of an infrared (IR)-bright dust-obscured galaxy (DOG) that shows a strong ionized-gas outflow but no significant molecular gas outflow. Based on detail analysis of their optical spectra, we found some peculiar IR-bright DOGs that show strong ionized-gas outflow ([OIII]$\lambda$5007) from the central active galactic nucleus (AGN). For one of these DOGs (WISE J102905.90+050132.4) at $z_{\rm spec} = 0.493$, we performed follow-up observations using ALMA to investigate their CO molecular gas properties. As a result, we successfully detected $^{12}$CO($J$=2-1) and $^{12}$CO($J$=4-3) lines, and the continuum of this DOG. The intensity-weighted velocity map of both lines shows a gradient, and the line profile of those CO lines is well-fitted by a single narrow Gaussian, meaning that this DOG has no sign of strong molecular gas outflow. The IR luminosity of this object is $\log\,(L_{\rm IR}/L_{\odot})$ = 12.40 that is classified as ultraluminous IR galaxy (ULIRG). We found that (i) the stellar mass and star-formation rate relation and (ii) the CO luminosity and far-IR luminosity relation are consistent with those of typical ULIRGs at similar redshifts. These results indicate that the molecular gas properties of this DOG are normal despite that its optical spectrum showing a powerful AGN outflow. We conclude that a powerful ionized-gas outflow caused by the AGN does not necessarily affect the cold interstellar medium in the host galaxy at least for this DOG., Comment: 8 pages, 5 figures, and 1 table, accepted for publication in ApJ

We present calibrations for star formation rate indicators in the ultraviolet, mid-infrared and radio continuum bands, including one of the first direct calibrations of 150 MHz as a star formation rate indicator. Our calibrations utilize 66 nearby star forming galaxies with Balmer decrement corrected H-alpha luminosities, which span 5 orders of magnitude in star formation rate and have absolute magnitudes of -24, Comment: Accepted for publication in the Astrophysical Journal; 39 pages, 14 figures and 4 tables

We constrain the quasar contribution to cosmic reionization based on our deep optical survey of z~6 quasars down to z_R=24.15 using Subaru/Suprime-Cam in three UKIDSS-DXS fields covering 6.5 deg^2. In Kashikawa et al. (2015), we select 17 quasar candidates and report our initial discovery of two low-luminosity quasars (M_1450~ -23) from seven targets, one of which might be a Lyman alpha emitting galaxy. From an additional optical spectroscopy, none of the four candidates out of the remaining ten turn out to be genuine quasars. Moreover, the deeper optical photometry provided by the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) shows that, unlike the two already-known quasars, the i-z and z-y colors of the last six candidates are consistent with M- or L-type brown dwarfs. Therefore, the quasar luminosity function (QLF) in the previous paper is confirmed. Compiling QLF measurements from the literature over a wide magnitude range, including an extremely faint AGN candidate from Parsa et al. (2017}, to fit them with a double power-law, we find that the best-fit faint-end slope is alpha=-2.04^+0.33_-0.18 (-1.98^+0.48_-0.21) and characteristic magnitude is M_1450^*=-25.8^+1.1_-1.9 (-25.7^+1.0_-1.8) in the case of two (one) quasar detection. Our result suggests that, if the QLF is integrated down to M_1450=-18, quasars produce ~1-12% of the ionizing photons required to ionize the whole universe at z~6 with 2sigma confidence level, assuming that the escape fraction is f_esc=1 and the IGM clumpy factor is C=3. Even when the systematic uncertainties are taken into account, our result supports the scenario that quasars are the minor contributors of reionization., Comment: 8 pages, 3 figures, ApJL accepted

We present ALMA cycle 0 observations of the luminous merger VV 114. One of the main goals is to investigate mechanisms of molecular line ratio enhancement. Regions with the high 12CO (1–0)/13CO (1–0) and 12CO (3–2)/12CO (1–0) is located at a central filamentary structure (∼6 kpc) in VV 114. The filament consists of the eastern nucleus and the overlap region, where the galaxy disks are colliding. We also investigate these molecular line ratios on the Kennicutt-Schmidt law. VV 114 fills a gap between the “starburst” sequence and the “normal disk” sequence, and regions with the high ratios show the high ΣSFR and ΣH2. We suggest that the high ratios in VV 114 are due to star-forming activities in the both progenitor's nuclei and the merger-induced overlap region.

We report the broadband X-ray spectra of the ultra-luminous infrared galaxy (ULIRG) UGC 5101 in the 0.25-100 keV band observed with Swift/Burst Alert Telescope (BAT), NuSTAR, Suzaku, XMM-Newton, and Chandra. A Compton-thick AGN obscured with a hydrogen column density of $\approx 1.3\times10^{24}$ cm$^{-2}$ is detected above 10 keV. A spectral fit with a numerical torus model favors a large half opening angle of the torus, $>41$ degrees, suggesting that the covering fraction of material heavily obscuring the X-ray source is moderate. The intrinsic 2-10 keV luminosity is determined to be $\approx 1.4\times 10^{43}$ erg s$^{-1}$, which is $\approx$2.5 times larger than the previous estimate using only data below 10 keV with a simple spectral model. We find that UGC 5101 shows the ratio between the [O IV] 26 $��$m line and 2-10 keV luminosities similar to those of normal Seyfert galaxies, along with other ULIRGs observed with NuSTAR, indicating that a significant portion of local ULIRGs are not really "X-ray faint" with respect to the flux of forbidden lines originating from the narrow line region (NLR). We propose a possible scenario that (1) the AGN in UGC 5101 is surrounded not only by Compton-thick matter located close to the equatorial plane but also by Compton-thin ($N_\mathrm{H} \sim 10^{21}$ cm$^{-2}$) matter in the torus-hole region and (2) it is accreting at a high Eddington rate with a steep UV to X-ray spectral energy distribution. Nevertheless, we argue that AGNs in many ULIRGs do not look extraordinary (i.e., extremely X-ray faint), as suggested by recent works, compared with normal Seyferts., 11 pages, 7 figures, accepted for publication in ApJ

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