• We report the result of optical identifications of FIRST radio sources with the Hyper Suprime-Cam Subaru Strategic Program survey (HSC-SSP). The positional cross-match within 1" between the FIRST and HSC-SSP catalogs (i ~< 26) produced more than 3600 optical counterparts in the 156 deg^2 of the HSC-SSP field. The matched counterparts account for more than 50 % of the FIRST sources in the search field, which substantially exceeds previously reported fractions of SDSS counterparts (i ~< 22) of ~ 30 %. Among the matched sample, 9 % are optically unresolved sources such as radio-loud quasars. The optically faint (i > 21) radio galaxies (RGs) show a flatter slope of a fitting linear function of the 1.4 GHz source counts than the bright RGs, while optically faint radio quasars show a steeper slope than the bright radio quasars. The optically faint RGs shows a flat slope in the i-band number counts down to 24 mag, implying either less-massive or distant radio-AGNs beyond 24 mag. The photometric redshift and the comparison of colors with the galaxy models show that most of the matched RGs are distributed at redshifts from 0 to 1.5. The optically faint sample includes the high radio-loudness sources that are not seen in the optically bright sample. Such sources are located at a redshift of larger than 1. This study provides a large number of radio-AGNs lying at the optically faint end and high redshift regime that are not probed by the previous searches.
  • We report optical properties of infrared (IR)-bright dust-obscured galaxies (DOGs) that are defined as (i-[22])_AB > 7.0. Since supermassive black holes (SMBHs) in IR-bright DOGs are expected to be rapidly growing in the major merger scenario, they provide useful clues to understand the co-evolution of SMBHs and their host galaxies. However the optical properties of IR-bright DOGs are still unclear, because the optical emission of DOGs is very faint. By combining ~ 53 deg^2 images of the optical, near-IR, and mid-IR data that are obtained with the Subaru Hyper Suprime Cam (HSC) survey, VISTA VIKING survey, and WISE all-sky survey, respectively, 427 IR-bright DOGs are selected. We find that IR-bright DOGs show a redder (g-z)_AB color than other populations of dusty galaxies such as ultra-luminous IR galaxies (ULIRGs) at a similar redshift, with a significantly large dispersion. Among the selected DOGs, star-formation (SF)-dominated DOGs show a relatively red color while active galactic nucleus (AGN)-dominated DOGs show a rather blue color in optical. This result is consistent with the idea that the relative AGN contribution in the optical emission becomes more significant at a late stage in the major merger scenario. Roughly 2% of AGN-dominated DOGs show a significant blue excess in blue HSC bands (BluDOGs), which can be interpreted as a leaked AGN emission as proposed for some blue-excess Hot DOGs in earlier works.
  • Narrow-line regions excited by active galactic nuclei (AGN) are important for studying AGN photoionization and feedback. Their strong [O III] lines can be detected with broadband images, allowing morphological studies of these systems with large-area imaging surveys. We develop a new technique to reconstruct the [O III] images using the Subaru Hyper Suprime-Cam (HSC) Survey aided with spectra from the Sloan Digital Sky Survey (SDSS). The technique involves a careful subtraction of the galactic continuum to isolate emission from the [O III]$\lambda$5007 and [O III]$\lambda$4959 lines. Compared to traditional targeted observations, this technique is more efficient at covering larger samples with less dedicated observational resources. We apply this technique to an SDSS spectroscopically selected sample of 300 obscured AGN at redshifts 0.1 - 0.7, uncovering extended emission-line region candidates with sizes up to tens of kpc. With the largest sample of uniformly derived narrow-line region sizes, we revisit the narrow-line region size-luminosity relation. The area and radii of the [O III] emission-line regions are strongly correlated with the AGN luminosity inferred from the mid-infrared (15 $\mu$m rest-frame) with a power-law slope of $0.62^{+0.05}_{-0.06} \pm 0.10$ (statistical and systemic errors), consistent with previous spectroscopic findings. We discuss the implications for the physics of AGN emission-line region and future applications of this technique, which should be useful for current and next-generation imaging surveys to study AGN photoionization and feedback with large statistical samples.
  • We present the discovery of an extremely-luminous dust-obscured galaxy (DOG) at $z_{\rm spec}$ = 3.703, WISE J101326.25+611220.1. This DOG is selected as a candidate of extremely-luminous infrared (IR) galaxies based on the photometry from the Sloan Digital Sky Survey and Wide-field Infrared Survey Explorer. In order to derive its accurate IR luminosity, we perform follow-up observations at 450 and 850 $\mu$m using the Submillimetre Common User Bolometer Array 2 on the James Clerk Maxwell Telescope, and at 870 and 1300 $\mu$m using the Submillimeter Array, which enable us to pin down its IR Spectral Energy Distribution (SED). We perform SED fitting using 14 photometric data (0.4 - 1300 $\mu$m) and estimate its IR luminosity, $L_{\rm IR}$ (8-1000 $\mu$m), to be $2.2^{+1.5}_{-1.0}$ $\times 10^{14}$ $L_{\odot}$, making it one of the most luminous IR galaxies in the Universe. The energy contribution from an active galactic nucleus (AGN) to the IR luminosity is $94^{+6}_{-20}$%, which indicates it is an AGN-dominated DOG. On the other hand, its stellar mass ($M_*$) and star formation rate (SFR) are $\log \,(M_\ast/M_{\odot})$ = $11.2^{+0.6}_{-0.2}$ and $\log \,({\rm SFR}/M_{\odot}\,{\rm yr}^{-1}$) = $3.1^{+0.2}_{-0.1}$, respectively, which means that this DOG can be considered as a starburst galaxy in $M_*$--SFR plane. This extremely-luminous DOG shows significant AGN and star forming activity that provides us an important laboratory to probe the maximum phase of the co-evolution of galaxies and supermassive black holes.
  • 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 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.
  • The SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES) is a three-year JCMT Large Program aiming at reaching the 450 $\mu$m confusion limit in the COSMOS-CANDELS region, to study a representative sample of the high-redshift far-infrared galaxy population that gives rise to the bulk of the far-infrared background. We present the first-year data from STUDIES. We have reached a 450 $\mu$m noise level of 0.91~mJy for point sources at the map center, covered an area of 151 arcmin$^2$, and detected 98 and 141 sources at 4.0 and 3.5 $\sigma$, respectively. Our derived counts are best constrained in the 3.5-25 mJy regime using directly detected sources. Below the detection limits, our fluctuation analysis further constrains the slope of the counts down to 1 mJy. The resulting counts at 1-25 mJy are consistent with a power law having a slope of $-2.59$ ($\pm0.10$ for 3.5-25 mJy, and $^{+0.4}_{-0.7}$ for 1-3.5 mJy). There is no evidence of a faint-end termination or turn-over of the counts in this flux density range. Our counts are also consistent with previous SCUBA-2 blank-field and lensing cluster surveys. The integrated surface brightness from our counts down to 1 mJy is $90.0\pm17.2$ Jy deg$^{-2}$, which can account for up to $83^{+15}_{-16}\%$ of the COBE 450 $\mu$m background. We show that Herschel counts at 350 and 500 $\mu$m are significantly higher than our 450 $\mu$m counts, likely caused by its large beam and source clustering. High-angular resolution instruments like SCUBA-2 at 450 $\mu$m are therefore highly beneficial for measuring the luminosity and spatial density of high-redshift dusty galaxies.
  • We present the ionized gas properties of infrared (IR)-bright dust-obscured galaxies (DOGs) that show an extreme optical/IR color, $(i - [22])_{\rm AB} > 7.0$, selected with the Sloan Digital Sky Survey (SDSS) and Wide-field Infrared Survey Explorer (WISE). For 36 IR-bright DOGs that show [OIII]$\lambda$5007 emission in the SDSS spectra, we performed a detailed spectral analysis to investigate their ionized gas properties. In particular, we measured the velocity offset (the velocity with respect to the systemic velocity measured from the stellar absorption lines) and the velocity dispersion of the [OIII] line. We found that the derived velocity offset and dispersion of most IR-bright DOGs are larger than those of Seyfert 2 galaxies (Sy2s) at $z < 0.3$, meaning that the IR-bright DOGs show relatively strong outflows compared to Sy2s. This can be explained by the difference of IR luminosity contributed from active galactic nucleus, $L_{\rm IR}$ (AGN), because we found that (i) $L_{\rm IR}$ (AGN) correlates with the velocity offset and dispersion of [OIII] and (ii) our IR-bright DOGs sample has larger $L_{\rm IR}$ (AGN) than Sy2s. Nevertheless, the fact that about 75% IR-bright DOGs have a large ($>$ 300 km s$^{-1}$) velocity dispersion, which is a larger fraction compared to other AGN populations, suggests that IR-bright DOGs are good laboratories to investigate AGN feedback. The velocity offset and dispersion of [OIII] and [NeIII]$\lambda$3869 are larger than those of [OII]$\lambda$3727, which indicates that the highly ionized gas tends to show more stronger outflows.
  • We construct a sample of X-ray bright optically faint active galactic nuclei by combining Subaru Hyper Suprime-Cam, XMM-Newton, and infrared source catalogs. 53 X-ray sources satisfying i band magnitude fainter than 23.5 mag and X-ray counts with EPIC-PN detector larger than 70 are selected from 9.1 deg^2, and their spectral energy distributions (SEDs) and X-ray spectra are analyzed. 44 objects with an X-ray to i-band flux ratio F_X/F_i>10 are classified as extreme X-ray-to-optical flux sources. SEDs of 48 among 53 are represented by templates of type 2 AGNs or starforming galaxies and show signature of stellar emission from host galaxies in the optical in the source rest frame. Infrared/optical SEDs indicate significant contribution of emission from dust to infrared fluxes and that the central AGN is dust obscured. Photometric redshifts determined from the SEDs are in the range of 0.6-2.5. X-ray spectra are fitted by an absorbed power law model, and the intrinsic absorption column densities are modest (best-fit log N_H = 20.5-23.5 cm^-2 in most cases). The absorption corrected X-ray luminosities are in the range of 6x10^42 - 2x10^45 erg s^-1. 20 objects are classified as type 2 quasars based on X-ray luminsosity and N_H. The optical faintness is explained by a combination of redshifts (mostly z>1.0), strong dust extinction, and in part a large ratio of dust/gas.
  • We have measured the clustering of galaxies around active galactic nuclei (AGN) for which single-epoch virial masses of the super-massive black hole (SMBH) are available to investigate the relation between the large scale environment of AGNs and the evolution of SMBHs. The AGN samples used in this work were derived from the Sloan Digital Sky Survey (SDSS) observations and the galaxy samples were from 240~deg$^{2}$ S15b data of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). The investigated redshift range is 0.6--3.0, and the masses of the SMBHs lie in the range $10^{7.5}$--$10^{10}$~M$_{\odot}$. The absolute magnitude of the galaxy samples reaches to $M_{\lambda 310}$ $\sim$ $-$18 at rest frame wavelength 310~nm for the low-redshift end of the samples. More than 70\% of the galaxies in the analysis are blue. We found a significant dependence of the cross-correlation length on redshift, which primarily reflects the brightness dependence of the galaxy clustering. At the lowest redshifts the cross-correlation length increases from 7~$h^{-1}$~Mpc around $M_{\lambda 310}$$ = $$-19$~mag to $>$10~$h^{-1}$~Mpc beyond $M_{\lambda 310}$$ = $$-20$~mag. No significant dependence of the cross-correlation length on BH mass was found for whole galaxy samples dominated by blue galaxies, while there was an indication of BH mass dependence in the cross-correlation with red galaxies. These results provides a picture of the environment of AGNs studied in this paper being enriched with blue starforming galaxies, and a fraction of the galaxies are being evolved to red galaxies along with the evolution of SMBHs in that system.
  • 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.
  • We present the environmental dependence of colour, stellar mass, and star formation (SF) activity in H-alpha-selected galaxies along the large-scale structure at z=0.4 hosting twin clusters in DEEP2-3 field, discovered by Subaru Strategic Programme of Hyper Suprime-Cam (HSC SSP). By combining photo-z selected galaxies and H-alpha emitters selected with broad-band and narrow-band (NB) data of the recent data release of HSC SSP (DR1), we confirm that galaxies in higher-density environments or galaxies in the cluster central regions show redder colours. We find that there still remains a possible colour-density and colour-radius correlation even if we restrict the sample to H-alpha-selected galaxies, likely due to the presence of massive H-alpha emitters in denser regions. We also find a hint of increased star formation rates (SFR) amongst H-alpha emitters towards the highest-density environment, again primarily driven by the excess of red/massive H-alpha emitters in high-density environment, while their specific SFR does not significantly change with environment. This work demonstrates the power of the HSC SSP NB data to study SF galaxies across environment in the distant universe.
  • 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/.
  • Low-metallicity active galactic nuclei (AGNs) are interesting to study the early phase of the AGN evolution. However most AGNs are chemically matured and accordingly low-metallicity AGNs are extremely rare. One approach to search for low-metallicity AGNs systematically is utilizing the so-called BPT diagram that consists of the [O III]$\lambda$5007/H$\beta$$\lambda4861$ and [N II]$\lambda6584$/H$\alpha$$\lambda6563$ flux ratios. Specifically, photoionization models predict that low-metallicity AGNs show a high [O III]$\lambda$5007/H$\beta$$\lambda$4861 ratio and a relatively low [N II]$\lambda$6584/H$\alpha$$\lambda$6563 ratio, that corresponds to the location between the sequence of star-forming galaxies and that of usual AGNs on the BPT diagram (hereafter "the BPT valley"). However, other populations of galaxies such as star-forming galaxies and AGNs with a high electron density or a high ionization parameter could be also located in the BPT valley, not only low-metallicity AGNs. In this paper, we examine whether most of emission-line galaxies at the BPT valley are low-metallicity AGNs or not. We select 70 BPT-valley objects from 212,866 emission line galaxies obtained by the Sloan Digital Sky Survey. Among the 70 BPT-valley objects, 43 objects show firm evidence of the AGN activity; i.e., the He II$\lambda$4686 emission and/or weak but significant broad H$\alpha$ emission. Our analysis shows that those 43 BPT-valley AGNs are not characterized by a very high gas density nor ionization parameter, inferring that at least 43 among 70 BPT-valley objects (i.e., $>60$ %) are low-metallicity AGNs. This suggests that the BPT diagram is an efficient tool to search for low-metallicity AGNs.
  • 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.
  • We examine the clustering of quasars over a wide luminosity range, by utilizing 901 quasars at $\overline{z}_{\rm phot}\sim3.8$ with $-24.73<M_{\rm 1450}<-22.23$ photometrically selected from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) S16A Wide2 date release and 342 more luminous quasars at $3.4<z_{\rm spec}<4.6$ having $-28.0<M_{\rm 1450}<-23.95$ from the Sloan Digital Sky Survey (SDSS) that fall in the HSC survey fields. We measure the bias factors of two quasar samples by evaluating the cross-correlation functions (CCFs) between the quasar samples and 25790 bright $z\sim4$ Lyman Break Galaxies (LBGs) in $M_{\rm 1450}<-21.25$ photometrically selected from the HSC dataset. Over an angular scale of \timeform{10.0"} to \timeform{1000.0"}, the bias factors are $5.93^{+1.34}_{-1.43}$ and $2.73^{+2.44}_{-2.55}$ for the low and high luminosity quasars, respectively, indicating no luminosity dependence of quasar clustering at $z\sim4$. It is noted that the bias factor of the luminous quasars estimated by the CCF is smaller than that estimated by the auto-correlation function (ACF) over a similar redshift range, especially on scales below \timeform{40.0"}. Moreover, the bias factor of the less-luminous quasars implies the minimal mass of their host dark matter halos (DMHs) is $0.3$-$2\times10^{12}h^{-1}M_{\odot}$, corresponding to a quasar duty cycle of $0.001$-$0.06$.
  • 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.
  • We investigate the star forming activity of a sample of infrared (IR)-bright dust-obscured galaxies (DOGs) that show an extreme red color in the optical and IR regime, $(i - [22])_{\rm AB} > 7.0$. Combining an IR-bright DOG sample with the flux at 22 $\mu$m $>$ 3.8 mJy discovered by Toba & Nagao (2016) with IRAS faint source catalog version 2 and AKARI far-IR (FIR) all-sky survey bright source catalog version 2, we selected 109 DOGs with FIR data. For a subsample of 7 IR-bright DOGs with spectroscopic redshift ($0.07 < z < 1.0$) that was obtained from literature, we estimated their IR luminosity, star formation rate (SFR), and stellar mass based on the spectral energy distribution fitting. We found that (i) WISE 22 $\mu$m luminosity at observed frame is a good indicator of IR luminosity for IR-bright DOGs and (ii) the contribution of active galactic nucleus (AGN) to IR luminosity increases with IR luminosity. By comparing the stellar mass and SFR relation for our DOG sample and literature, we found that most of IR-bright DOGs lie significantly above the main sequence of star-forming galaxies at similar redshift, indicating that the majority of IRAS- and/or AKARI-detected IR-bright DOGs are starburst galaxies.
  • 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 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 systematically investigate the near- (NIR) to far-infrared (FIR) photometric properties of a nearly complete sample of local active galactic nuclei (AGN) detected in the Swift/Burst Alert Telescope (BAT) all-sky ultra hard X-ray (14-195 keV) survey. Out of 606 non-blazar AGN in the Swift/BAT 70-month catalog at high galactic latitude of $|b|>10^{\circ}$, we obtain IR photometric data of 604 objects by cross-matching the AGN positions with catalogs from the WISE, AKARI, IRAS, and Herschel infrared observatories. We find a good correlation between the ultra-hard X-ray and mid-IR (MIR) luminosities over five orders of magnitude ($41 < \log (L_{14-195}/{\rm erg}~{\rm s}^{-1})< 46$). Informed by previous measures of the intrinsic spectral energy distribution of AGN, we find FIR pure-AGN candidates whose FIR emission is thought to be AGN-dominated with low starformation activity. We demonstrate that the dust covering factor decreases with the bolometric AGN luminosity, confirming the luminosity-dependent unified scheme. We also show that the completeness of the WISE color-color cut in selecting Swift/BAT AGN increases strongly with 14-195 keV luminosity.
  • One of the unsettled, important problems in active galactic nuclei (AGNs) is the major ionization mechanism of gas clouds in their narrow-line regions (NLRs). In order to investigate this issue, we present our $J$-band spectroscopic observations for a sample of 26 nearby Seyfert galaxies. In our study, we use the flux ratio of the following two forbidden emission lines, [Fe II]1.257$\mu$m and [P II]1.188$\mu$m because it is known that this ratio is sensitive to the ionization mechanism. We obtain the [Fe II]/[P II] flux ratio or its lower limit for 19 objects. In addition to our data, we compile this flux ratio (or its lower limit) of 23 nearby Seyfert galaxies from the literature. Based on the collected data, we find that three Seyfert galaxies show very large lower limits of [Fe II]/[P II] flux ratios ($\gtrsim$10); NGC 2782, 5005, and Mrk 463. It is thus suggested that the contribution of the fast shock in the gas excitation is significantly large for them. However, more than half of Seyfert galaxies in our sample show moderate [Fe II]/[P II] flux ratios ($\sim$2), that is consistent to the pure photoionization by power-law ionizing continuum emission. We also find that the [Fe II]/[P II] flux ratio shows no clear correlation with the radio loudness, suggesting that the radio jet is not the primary origin of shocks in NLRs of Seyfert galaxies.
  • Recent studies have derived quasar luminosity functions (QLFs) at various redshifts. However, the faint side of the QLF at high redshifts is still too uncertain. An accurate estimate of the survey completeness is essential to derive an accurate QLF for use in studying the luminosity-dependent density evolution of the quasar population. Here we investigate how the luminosity dependence of quasar spectra (the Baldwin effect) and the attenuation model for the inter-galactic medium (IGM) affect the completeness estimates. For this purpose, we revisit the completeness of quasar surveys specifically at $z\sim4-5$, using the COSMOS images observed with Subaru/Suprime-Cam. As the result, we find that the completeness estimates are sensitive to the luminosity dependence of the quasar spectrum and difference in the IGM attenuation models. At $z\sim4$, the number density of quasars when we adopt the latest IGM model and take the luminosity dependence of spectra into account are $(3.49\pm1.62)\times10^{-7}$ Mpc$^{-3}$ mag$^{-1}$ for $-24.09<M_{1450}<-23.09$ and $(5.24\pm2.13)\times10^{-7}$ Mpc$^{-3}$ mag$^{-1}$ for $-23.09<M_{1450}<-22.09$, respectively, which are $\sim$$24$$\%$ lower than that estimated by the conventional method. On the other hand, at $z\sim5$, a $1\sigma$ confidence upper limit of the number density at $-24.5<M_{1450}<-22.5$ in our new estimates is $\sim$$43$$\%$ higher than that estimated previously. The results suggest that the luminosity dependence of the quasar spectrum and the IGM model are important for deriving accurate number density of high-$z$ quasars. Even taking these effects into account, the inferred luminosity-dependent density evolution of quasars is consistent with the AGN down-sizing evolutionary picture.
  • 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.
  • Nearby active galactic nuclei were diagnosed in the X-ray and mid-to-far infrared wavelengths, with Monitor of All-sky X-ray Image (MAXI) and the Japanese infrared observatory AKARI, respectively. Among the X-ray sources listed in the second release of the MAXI all-sky X-ray source catalog, 100 ones are currently identified as a non-blazar-type active galactic nucleus. These include 95 Seyfert galaxies and 5 quasars, and they are composed of 73 type-1 and 27 type-2 objects. The AKARI all-sky survey point source catalog was searched for their mid- and far-infrared counterparts at 9, 18, and 90 $\mu$m. As a result, 69 Seyfert galaxies in the MAXI catalog (48 type-1 and 21 type-2 ones) were found to be detected with AKARI. The X-ray (3-4 keV and 4-10 keV) and infrared luminosities of these objects were investigated, together with their color information. Adopting the canonical photon index, $\Gamma = 1.9$, of the intrinsic X-ray spectrum of the Seyfert galaxies, the X-ray hardness ratio between the 3-4 and 4-10 keV ranges derived with MAXI was roughly converted into the absorption column density. After the X-ray luminosity was corrected for absorption from the estimated column density, the well-known X-ray-to-infrared luminosity correlation was confirmed at least in the Compton-thin regime. In contrast, NGC 1365, only one Compton-thick object in the MAXI catalog, was found to deviate from the correlation toward a significantly lower X-ray luminosity by nearly an order of magnitude. It was verified that the relation between the X-ray hardness below 10 keV and X-ray-to-infrared color acts as an effective tool to pick up Compton-thick objects. The difference in the infrared colors between the type-1 and type-2 Seyfert galaxies and its physical implication on the classification and unification of active galactic nuclei were briefly discussed.