• We present a broadband (0.4-70 keV) X-ray spectral analysis of the luminous infrared galaxy (LIRG) system Mrk 463 observed with Nuclear Spectroscopic Telescope Array (NuSTAR), Chandra, and XMM-Newton, which contains double active galactic nuclei (AGNs) (Mrk 463E and Mrk 463W) with a separation of $\sim$ 3.8 kpc. Detecting their transmitted hard X-ray $>10$ keV continua with NuSTAR, we confirm that Mrk 463E and Mrk 463W have AGNs with intrinsic luminosities of (1.6-2.2) $\times 10^{43}$ and (0.5-0.6) $\times 10^{43}$ erg s$^{-1}$ (2-10 keV) obscured by hydrogen column densities of 8 $\times 10^{23}$ and 3 $\times 10^{23}$ cm$^{-2}$, respectively. Both nuclei show strong reflection components from cold matter. The luminosity ratio between X-ray (2-10 keV) and [O IV] 25.89 $\mu$m of Mrk 463E is $\sim$5 times smaller than those of normal Seyfert galaxies, suggesting that the intrinsic SED is X-ray weak relative to the UV luminosity. In fact, the bolometric AGN luminosity of Mrk 463E estimated from $L^{\prime}$-band (3.8 $\mu$m), [O IV] 25.89 $\mu$m, and [Ne V] 14.32 $\mu$m lines indicate a large bolometric-to-X-ray luminosity ratio, $\kappa_{\rm 2-10 keV} \approx$ 110-410, and a high Eddington ratio, $\lambda_{\rm Edd} \sim$ 0.4-0.8. We suggest that the merger triggered a rapid growth of the black hole in Mrk 463E, which is not yet deeply "buried" by circumnuclear dust. By contrast, the $L^{\prime}$-band luminosity of Mrk 463W is unusually small relative to the X-ray luminosity, suggesting that the Eddington ratio is low ($<10^{-3}$) and it might be still in an early phase of merger-driven AGN activity.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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<M_r<-12. Most of our photometry and spectrophotometry is measured from the same region of each galaxy, and our spectrophotometry has been validated with SDSS photometry, so our random and systematic errors are small relative to the intrinsic scatter seen in star formation rate indicator calibrations. We find WISE W4 (22.8 micron), Spitzer 24 micron and 1.4 GHz have tight correlations with Balmer decrement corrected H-alpha luminosity, with scatter of only 0.2 dex. Our calibrations are comparable to those from the prior literature for L* galaxies, but for dwarf galaxies our calibrations can give star formation rates that are far greater than those derived from much of the prior literature.
  • 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 i<25. The projected separations of the two pairs are R_perp=1.75 and 1.04 proper Mpc, and their velocity offsets are Delta V=692 and 1448 km s^{-1}, respectively. This result is in clear contrast to the average z~4 quasar environments as discussed in Uchiyama et al. (2017) and implies that the quasar activities of the pair members are triggered via major mergers in proto-clusters, unlike the vast majority of isolated quasars in general fields that may turn on via non-merger events such as bar and disk instabilities. At z~1, we find 37 pairs with R_perp<2 pMpc and Delta V<2300 km s^{-1} in the current HSC-Wide coverage, including four from Hennawi et al. (2006). The distribution of the peak overdensity significance within two arcminutes around the pairs has a long tail toward high density (>4sigma) 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<1 pMpc) pairs, two of which are found to reside in cluster fields. Our results demonstrate that <2 pMpc-scale quasar pairs at both redshift range tend to occur in massive haloes, although perhaps not the most massive ones, and that they are useful to search for rare density peaks.
  • Hiroaki Aihara, Robert Armstrong, Steven Bickerton, James Bosch, Jean Coupon, Hisanori Furusawa, Yusuke Hayashi, Hiroyuki Ikeda, Yukiko Kamata, Hiroshi Karoji, Satoshi Kawanomoto, Michitaro Koike, Yutaka Komiyama, Robert H. Lupton, Sogo Mineo, Hironao Miyatake, Satoshi Miyazaki, Tomoki Morokuma, Yoshiyuki Obuchi, Yukie Oishi, Yuki Okura, Paul A. Price, Tadafumi Takata, Manobu M. Tanaka, Masayuki Tanaka, Yoko Tanaka, Tomohisa Uchida, Fumihiro Uraguchi, Yousuke Utsumi, Shiang-Yu Wang, Yoshihiko Yamada, Hitomi Yamanoi, Naoki Yasuda, Nobuo Arimoto, Masashi Chiba, Francois Finet, Hiroki Fujimori, Seiji Fujimoto, Junko Furusawa, Tomotsugu Goto, Andy Goulding, James E. Gunn, Yuichi Harikane, Takashi Hattori, Masao Hayashi, Krzysztof G. Helminiak, Ryo Higuchi, Chiaki Hikage, Paul T.P. Ho, Bau-Ching Hsieh, Kuiyun Huang, Song Huang, Masatoshi Imanishi, Ikuru Iwata, Anton T. Jaelani, Hung-Yu Jian, Nobunari Kashikawa, Nobuhiko Katayama, Takashi Kojima, Akira Konno, Shintaro Koshida, Haruka Kusakabe, Alexie Leauthaud, C.-H. Lee, Lihwai Lin, Yen-Ting Lin, Rachel Mandelbaum, Yoshiki Matsuoka, Elinor Medezinski, Shoken Miyama, Rieko Momose, Anupreeta More, Surhud More, Shiro Mukae, Ryoma Murata, Hitoshi Murayama, Tohru Nagao, Fumiaki Nakata, Hiroko Niikura, Atsushi J. Nishizawa, Masamune Oguri, Nobuhiro Okabe, Yoshiaki Ono, Masato Onodera, Masafusa Onoue, Masami Ouchi, Tae-Soo Pyo, Takatoshi Shibuya, Kazuhiro Shimasaku, Melanie Simet, Joshua Speagle, David N. Spergel, Michael A. Strauss, Yuma Sugahara, Naoshi Sugiyama, Yasushi Suto, Nao Suzuki, Philip J. Tait, Masahiro Takada, Tsuyoshi Terai, Yoshiki Toba, Edwin L. Turner, Hisakazu Uchiyama, Keiichi Umetsu, Yuji Urata, Tomonori Usuda, Sherry Yeh, Suraphong Yuma
    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/.
  • 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.
  • Active galactic nuclei (AGNs) in the high-redshift Universe are thought to reside in overdense environments. However, recent works provide controversial results partly due to the use of different techniques and possible suppression of nearby galaxy formation by AGN feedback. We conducted deep and wide-field imaging observations with the Suprime-Cam on the Subaru telescope and searched for Lyman-alpha emitters (LAEs) around two QSOs (quasi-stellar objects) at $z\sim4.9$ and a radio galaxy at $z\sim4.5$ by using narrow-band filters to address these issues more robustly. In the QSO fields, we obtained additional broad-band images to select Lyman-break galaxies (LBGs) at $z\sim5$ for comparison. We constructed a photometric sample of 301 LAEs and 170 LBGs in total. A wide field of view (34arcmin$\times$27arcmin, corresponding to 80$\times$60 comoving Mpc$^2$) of the Suprime-Cam enabled us to probe galaxies in the immediate vicinities of the AGNs and in the blank fields simultaneously and compare various properties of them in a consistent manner. The two QSOs are located near local density peaks ($<2\sigma$) and one of the QSOs has a close companion LAE with projected separation of 80 physical kpc. The radio galaxy is found to be near a void of LAEs. The number densities of LAEs/LGBs in a larger spatial scale around the AGNs are not significantly different from those in blank fields. No sign of feedback is found down to $L_{Ly\alpha}\sim10^{41.8}\mathrm{~erg~s^{-1}}$. Our results suggest that high-redshift AGNs are not associated with extreme galaxy overdensity and that this cannot be attributed to the effect of AGN feedback.
  • 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.
  • In order to investigate the physical relationship between strong Ly-alpha absorbers (logN_HI>20.0 cm^-2) such as damped Ly-alpha absorption systems (DLAs) and young star-forming galaxies at high redshift, we have conducted narrow-band observations of Ly-alpha emitters (LAEs) in a concentrated region of strong Ly-alpha absorbers at z= 2.3, the J1230+34 field. Using a catalog of Ly-alpha absorbers with logN_HI>20.0 cm^-2 based on the baryon oscillation spectroscopic survey (BOSS), we found 6 fields where 3 or more absorbers are concentrated within a (50 Mpc)^3 cubic box in the comoving scale. Among them, we focus on the J1230+34 field, where 2 DLAs and 2 sub-DLAs present. Our narrow-band imaging observations with Subaru/Suprime-Cam using a custom-made filter, NB400 (lambda_c=4003 A and FWHM=92 A) yield a sample of 149 LAEs in this field. In the large scale (~50 Mpc), we have found no differences between the obtained Ly-alpha luminosity function and those in the blank fields at similar redshifts. We also compare the frequency distribution of the Ly-alpha rest-frame equivalent width (EW_0) in the target field and other fields including both overdensity region and blank field, but find no differences. On the other hand, in the small scale (~10 Mpc), we have found a possible overdensity of LAEs around a DLA with the highest HI column density (N_HI=21.08 cm^-2) in the target field while there are no density excess around the other absorbers with a lower N_HI.
  • 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 $\mu$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.
  • We present measurements of the clustering properties of a sample of infrared (IR) bright dust-obscured galaxies (DOGs). Combining 125 deg$^2$ of wide and deep optical images obtained with the Hyper Suprime-Cam on the Subaru Telescope and all-sky mid-IR (MIR) images taken with Wide-Field Infrared Survey Explorer, we have discovered 4,367 IR-bright DOGs with $(i - [22])_{\rm AB}$ $>$ 7.0 and flux density at 22 $\mu$m $>$ 1.0 mJy. We calculate the angular autocorrelation function (ACF) for a uniform subsample of 1411 DOGs with 3.0 mJy $<$ flux (22 $mu$m) $<$ 5.0 mJy and $i_{\rm AB}$ $<$ 24.0. The ACF of our DOG subsample is well-fit with a single power-law, $\omega (\theta)$ = (0.010 $\pm$ 0.003) $\theta^{-0.9}$, where $\theta$ in degrees. The correlation amplitude of IR-bright DOGs is larger than that of IR-faint DOGs, which reflects a flux-dependence of the DOG clustering, as suggested by Brodwin et al. (2008). We assume that the redshift distribution for our DOG sample is Gaussian, and consider 2 cases: (1) the redshift distribution is the same as IR-faint DOGs with flux at 22 $\mu$m $<$ 1.0 mJy, mean and sigma $z$ = 1.99 $\pm$ 0.45, and (2) $z$ = 1.19 $\pm$ 0.30, as inferred from their photometric redshifts. The inferred correlation length of IR-bright DOGs is $r_0$ = 12.0 $\pm$ 2.0 and 10.3 $\pm$ 1.7 $h^{-1}$ Mpc, respectively. IR-bright DOGs reside in massive dark matter halos with a mass of $\log [\langle M_{\mathrm{h}} \rangle / (h^{-1} M_{\odot})]$ = 13.57$_{-0.55}^{+0.50}$ and 13.65$_{-0.52}^{+0.45}$ in the two cases, respectively.
  • We report the detection of two CH$_3$OH lines (J$_K$ = 2$_K$-1$_K$ and 3$_K$-2$_K$) between the progenitor's disks ("Overlap") of the mid-stage merging galaxy VV 114 obtained using the Atacama Large Millimeter/submillimeter Array (ALMA) Band 3 and Band 4. The detected CH$_3$OH emission show an extended filamentary structure (~ 3 kpc) across the progenitor's disks with relatively large velocity width (FWZI ~ 150 km/s). The emission is only significant in the "overlap" and not detected in the two merging nuclei. Assuming optically-thin emission and local thermodynamic equilibrium (LTE), we found the CH$_3$OH column density relative to H$_2$ ($X_{\rm CH_3OH}$) peaks at the "Overlap" (~ 8 $\times$ 10$^{-9}$), which is almost an order of magnitude larger than that at the eastern nucleus. We suggest that kpc-scale shocks driven by galaxy-galaxy collision may play an important role to enhance the CH$_3$OH abundance at the "Overlap". This scenario is consistent with that shock-induced large velocity dispersion components of ionized gas have been detected in optical wavelength at the same region. Conversely, low $X_{\rm CH_3OH}$ at the nuclear regions might be attributed to the strong photodissociation by nuclear starbursts and/or putative active galactic nucleus (AGN), or inefficient production of CH$_3$OH on dust grains due to initial high temperature conditions (i.e., desorption of the precursor molecule, CO, into gas-phase before forming CH$_3$OH on dust grains). These ALMA observations demonstrate that CH$_3$OH is a unique tool to address kpc-scale shock-induced gas dynamics and star formation in merging galaxies.
  • We report the discovery of 15 quasars and bright galaxies at 5.7 < z < 6.9. This is the initial result from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the exquisite multiband imaging data produced by the Subaru Hyper Suprime-Cam (HSC) Strategic Program survey. The candidate selection is performed by combining several photometric approaches including a Bayesian probabilistic algorithm to reject stars and dwarfs. The spectroscopic identification was carried out with the Gran Telescopio Canarias and the Subaru Telescope for the first 80 deg2 of the survey footprint. The success rate of our photometric selection is quite high, approaching 100 % at the brighter magnitudes (zAB < 23.5 mag). Our selection also recovered all the known high-z quasars on the HSC images. Among the 15 discovered objects, six are likely quasars, while the other six with interstellar absorption lines and in some cases narrow emission lines are likely bright Lyman-break galaxies. The remaining three objects have weak continua and very strong and narrow Ly alpha lines, which may be excited by ultraviolet light from both young stars and quasars. These results indicate that we are starting to see the steep rise of the luminosity function of z > 6 galaxies, compared with that of quasars, at magnitudes fainter than M1450 ~ -22 mag or zAB ~24 mag. Follow-up studies of the discovered objects as well as further survey observations are ongoing.
  • 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.
  • Compiling data from literature and the ALMA archive, we show enhanced HCN(4-3)/HCO$^+$(4-3) and/or HCN(4-3)/CS(7-6) integrated intensity ratios in circumnuclear molecular gas around active galactic nuclei (AGNs) compared to those in starburst (SB) galaxies (submillimeter HCN-enhancement). The number of sample galaxies is significantly increased from our previous work. We expect this feature could potentially be an extinction-free energy diagnostic tool of nuclear regions of galaxies. Non-LTE radiative transfer modelings of the above molecular emission lines involving both collisional and radiative excitation, as well as a photon trapping effect were conducted to investigate the cause of the high line ratios in AGNs. As a result, we found that enhanced abundance ratios of HCN-to-HCO$^+$ and HCN-to-CS in AGNs as compared to SB galaxies by a factor of a few to even $>$ 10 is a plausible explanation for the submillimeter HCN-enhancement. However, a counter argument of a systematically higher gas density in AGNs than in SB galaxies can also be a plausible scenario. Although we could not fully discriminate these two scenarios at this moment due to insufficient amount of multi-transition, multi-species data, the former scenario equivalently claims for abnormal chemical composition in AGNs. Regarding the actual mechanism to realize the composition, we suggest it is difficult with conventional gas phase X-ray dominated region (XDR) ionization models to reproduce the observed high line ratios. We might have to take into account other mechanisms such as neutral-neutral reactions that are efficiently activated at high temperature environments and/or mechanically heated regions to further understand the high line ratios in AGNs.
  • We present ALMA Cycle 1 observations of the central kpc region of the luminous type-1 Seyfert galaxy NGC 7469 with unprecedented high resolution (0.5$"$ $\times$ 0.4$"$ = 165 pc $\times$ 132 pc) at submillimeter wavelengths. Utilizing the wide-bandwidth of ALMA, we simultaneously obtained HCN(4-3), HCO$^+$(4-3), CS(7-6), and partially CO(3-2) line maps, as well as the 860 $\mu$m continuum. The region consists of the central $\sim$ 1$"$ component and the surrounding starburst ring with a radius of $\sim$ 1.5$"$-2.5$"$. Several structures connect these components. Except for CO(3-2), these dense gas tracers are significantly concentrated towards the central $\sim$ 1$"$, suggesting their suitability to probe the nuclear regions of galaxies. Their spatial distribution resembles well those of centimeter and mid-infrared continuum emissions, but it is anti-correlated with the optical one, indicating the existence of dust obscured star formation. The integrated intensity ratios of HCN(4-3)/HCO$^+$(4-3) and HCN(4-3)/CS(7-6) are higher at the AGN position than at the starburst ring, which is consistent to our previous findings (submm-HCN enhancement). However, the HCN(4-3)/HCO$^+$(4-3) ratio at the AGN position of NGC 7469 (1.11$\pm$0.06) is almost half of the corresponding value of the low-luminosity type-1 Seyfert galaxy NGC 1097 (2.0$\pm$0.2), despite the more than two orders of magnitude higher X-ray luminosity of NGC 7469. But the ratio is comparable to that of the close vicinity of the AGN of NGC 1068 ($\sim$ 1.5). Based on these results, we speculate that some other heating mechanisms than X-ray (e.g., mechanical heating due to AGN jet) can contribute significantly for shaping the chemical composition in NGC 1097.
  • We present the results from a large near-infrared spectroscopic survey with Subaru/FMOS (\textit{FastSound}) consisting of $\sim$ 4,000 galaxies at $z\sim1.4$ with significant H$\alpha$ detection. We measure the gas-phase metallicity from the [N~{\sc ii}]$\lambda$6583/H$\alpha$ emission line ratio of the composite spectra in various stellar mass and star-formation rate bins. The resulting mass-metallicity relation generally agrees with previous studies obtained in a similar redshift range to that of our sample. No clear dependence of the mass-metallicity relation with star-formation rate is found. Our result at $z\sim1.4$ is roughly in agreement with the fundamental metallicity relation at $z\sim0.1$ with fiber aperture corrected star-formation rate. We detect significant [S~{\sc ii}]$\lambda\lambda$6716,6731 emission lines from the composite spectra. The electron density estimated from the [S~{\sc ii}]$\lambda\lambda$6716,6731 line ratio ranges from 10 -- 500 cm$^{-3}$, which generally agrees with that of local galaxies. On the other hand, the distribution of our sample on [N~{\sc ii}]$\lambda$6583/H$\alpha$ vs. [S~{\sc ii}]$\lambda\lambda$6716,6731/H$\alpha$ is different from that found locally. We estimate the nitrogen-to-oxygen abundance ratio (N/O) from the N2S2 index, and find that the N/O in galaxies at $z\sim1.4$ is significantly higher than the local values at a fixed metallicity and stellar mass. The metallicity at $z\sim1.4$ recalculated with this N/O enhancement taken into account decreases by 0.1 -- 0.2 dex. The resulting metallicity is lower than the local fundamental metallicity relation.
  • We present an unprecedented measurement of the disc stability and local instability scales in the luminous infrared Seyfert 1 host, NGC7469, based on ALMA observations of dense gas tracers and with a synthesized beam of 165 x 132 pc. While we confirm that non-circular motions are not significant in redistributing the dense interstellar gas in this galaxy, we find compelling evidence that the dense gas is a suitable tracer for studying the origin of its intensely high-mass star forming ring-like structure. Our derived disc stability parameter accounts for a thick disc structure and its value falls below unity at the radii in which intense star formation is found. Furthermore, we derive the characteristic instability scale and find a striking agreement between our measured scale of ~ 180 pc, and the typical sizes of individual complexes of young and massive star clusters seen in high-resolution images.
  • We present the photometric properties of a sample of infrared (IR) bright dust obscured galaxies (DOGs). Combining wide and deep optical images obtained with the Hyper Suprime-Cam (HSC) on the Subaru Telescope and all-sky mid-IR (MIR) images taken with Wide-Field Infrared Survey Explorer (WISE), we discovered 48 DOGs with $i - K_\mathrm{s} > 1.2$ and $i - [22] > 7.0$, where $i$, $K_\mathrm{s}$, and [22] represent AB magnitude in the $i$-band, $K_\mathrm{s}$-band, and 22 $\mu$m, respectively, in the GAMA 14hr field ($\sim$ 9 deg$^2$). Among these objects, 31 ($\sim$ 65 %) show power-law spectral energy distributions (SEDs) in the near-IR (NIR) and MIR regime, while the remainder show a NIR bump in their SEDs. Assuming that the redshift distribution for our DOGs sample is Gaussian, with mean and sigma $z$ = 1.99 $\pm$ 0.45, we calculated their total IR luminosity using an empirical relation between 22 $\mu$m luminosity and total IR luminosity. The average value of the total IR luminosity is (3.5 $\pm$ 1.1) $\times$ $10^{13}$ L$_{\odot}$, which classifies them as hyper-luminous infrared galaxies (HyLIRGs). We also derived the total IR luminosity function (LF) and IR luminosity density (LD) for a flux-limited subsample of 18 DOGs with 22 $\mu$m flux greater than 3.0 mJy and with $i$-band magnitude brighter than 24 AB magnitude. The derived space density for this subsample is log $\phi$ = -6.59 $\pm$ 0.11 [Mpc$^{-3}$]. The IR LF for DOGs including data obtained from the literature is well fitted by a double-power law. The derived lower limit for the IR LD for our sample is $\rho_{\mathrm{IR}}$ $\sim$ 3.8 $\times$ 10$^7$ [L$_{\odot}$ Mpc$^{-3}$] and its contributions to the total IR LD, IR LD of all ultra-luminous infrared galaxies (ULIRGs), and that of all DOGs are $>$ 3 %, $>$ 9 %, and $>$ 15 %, respectively.
  • We report the discovery of a new dwarf galaxy (NGC6503-d1) during the Subaru extended ultraviolet (XUV) disk survey. It is a likely companion of the spiral galaxy NGC6503. The resolved images, in B, V, R, i, and Halpha, show an irregular appearance due to bright stars with underlying, smooth and unresolved stellar emission. It is classified as the transition type (dIrr/dSph). Its structural properties are similar to those of the dwarfs in the Local Group, with a V absolute magnitude ~ -10.5, half-light radius ~400 pc, and central surface brightness ~25.2. Despite the low stellar surface brightness environment, one HII region was detected, though its Halpha luminosity is low, indicating an absence of any appreciable O-stars at the current epoch. The presence of multiple stellar populations is indicated by the color-magnitude diagram of ~300 bright resolved stars and the total colors of the dwarf, with the majority of its total stellar mass ~4x10^6 Msun in an old stellar population.
  • We present ALMA cycle 0 observations of the molecular gas and dust in the IR-bright mid-stage merger VV114 obtained at 160 - 800 pc resolution. The main aim of this study is to investigate the distribution and kinematics of the cold/warm gas and to quantify the spatial variation of the excitation conditions across the two merging disks. The data contain 10 molecular lines, including the first detection of extranuclear CH3OH emission in interacting galaxies, as well as continuum emission. We map the 12CO(3-2)/12CO(1-0) and the 12CO(1-0)/13CO(1-0) line ratio at 800 pc resolution (in the units of K km/s), and find that these ratios vary from 0.2 - 0.8 and 5 - 50, respectively. Conversely, the 200 pc resolution HCN(4-3)/HCO+(4-3) line ratio shows low values (< 0.5) at a filament across the disks except for the unresolved eastern nucleus which is three times higher (1.34 +/- 0.09). We conclude from our observations and a radiative transfer analysis that the molecular gas in the VV114 system consists of five components with different physical and chemical conditions; i.e., 1) dust-enshrouded nuclear starbursts and/or AGN, 2) wide-spread star forming dense gas, 3) merger-induced shocked gas, 4) quiescent tenuous gas arms without star formation, 5) H2 gas mass of (3.8 +/- 0.7) * 10^7 Msun (assuming a conversion factor of {\alpha}_CO = 0.8 Msun (K km s^-1 pc^2)^-1) at the tip of the southern tidal arm, as a potential site of tidal dwarf galaxy formation.
  • We present results from the fitting of infrared (IR) spectral energy distributions of 21 active galactic nuclei (AGN) with clumpy torus models. We compiled high spatial resolution ($\sim 0.3$--$0.7$ arcsec) mid-IR $N$-band spectroscopy, $Q$-band imaging and nuclear near- and mid-IR photometry from the literature. Combining these nuclear near- and mid-IR observations, far-IR photometry and clumpy torus models, enables us to put constraints on the torus properties and geometry. We divide the sample into three types according to the broad line region (BLR) properties; type-1s, type-2s with scattered or hidden broad line region (HBLR) previously observed, and type-2s without any published HBLR signature (NHBLR). Comparing the torus model parameters gives us the first quantitative torus geometrical view for each subgroup. We find that NHBLR AGN have smaller torus opening angles and larger covering factors than those of HBLR AGN. This suggests that the chance to observe scattered (polarized) flux from the BLR in NHBLR could be reduced by the dual effects of (a) less scattering medium due to the reduced scattering volume given the small torus opening angle and (b) the increased torus obscuration between the observer and the scattering region. These effects give a reasonable explanation for the lack of observed HBLR in some type-2 AGN.