We present the stellar population and ionized-gas outflow properties of ultra-luminous IR galaxies (ULIRGs) at $z=$ 0.1-1.0, which are selected from AKARI FIR all-sky survey. We construct a catalog of 1077 ULIRGs to examine feedback effect after major mergers. 202 out of the 1077 ULIRGs are spectroscopically identified by SDSS and Subaru/FOCAS observations. Thanks to deeper depth and higher resolution of AKARI compared to the previous IRAS survey, and reliable identification from WISE MIR pointing, the sample is unique in identifying optically-faint (i$\sim$20) IR-bright galaxies, which could be missed in previous surveys. A self-consistent spectrum-SED decomposition method, which constrains stellar population properties in SED modeling based on spectral fitting results, has been employed for 149 ULIRGs whose optical continua are dominated by host galaxies. They are massive galaxies ($M_{\rm star}\sim10^{11}$-$10^{12}$ M$_{\odot}$), associated with intense star formation activities (SFR $\sim$ 200-2000 M$_{\odot}$ yr$^{-1}$). The sample covers a range of AGN bolometric luminosity of $10^{10}$-$10^{13}$ L$_{\odot}$, and the outflow velocity measured from [OIII] 5007A line shows a correlation with AGN luminosity. Eight galaxies show extremely fast outflows with velocity up to 1500-2000 km s$^{-1}$. However, the co-existence of vigorous starbursts and strong outflows suggests the star formation has not been quenched during the ULIRG phase. By deriving stellar mass and mass fraction of young stellar population, we find no significant discrepancies between stellar properties of ULIRGs with weak and powerful AGNs. The results are not consistent with the merger-induced evolutionary scenario, which predicts that SF-dominated ULIRGs show smaller stellar mass and younger stellar populations compared to AGN-dominated ULIRGs.

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 SDSS quasars, we identify those with two (or more) distinct optical point sources in HSC images covering 796 deg^2. 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/LRIS and Gemini-N/NIFS, 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 [OIII]5007 emission lines have mild velocity offsets, thus the joint [OIII] 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.

We present the physical properties of AKARI sources without optical counterparts in optical images from the Hyper Suprime-Cam (HSC) on the Subaru telescope. Using the AKARI infrared (IR) source catalog and HSC optical catalog, we select 583 objects that do not have HSC counterparts in the AKARI North Ecliptic Pole (NEP) wide survey field ($\sim 5$ deg$^{2}$). Because the HSC limiting magnitude is deep ($g_{\rm AB}$ $\sim 28.6$), these are good candidates for extremely red star-forming galaxies (SFGs) and/or active galactic nuclei (AGNs), possibly at high redshifts. We compile multi-wavelength data out to 500 $\mu$m and use it for Spectral Energy Distribution (SED) fitting with CIGALE to investigate the physical properties of AKARI galaxies without optical counterparts. We also compare their physical quantities with AKARI mid-IR selected galaxies with HSC counterparts. The estimated redshifts of AKARI objects without HSC counterparts range up to $z\sim 4$, significantly higher than that of AKARI objects with HSC counterparts. We find that: (i) 3.6 $-$ 4.5 $\mu$m color, (ii) AGN luminosity, (iii) stellar mass, (iv) star formation rate, and (v) $V$-band dust attenuation in the interstellar medium of AKARI objects without HSC counterparts are systematically larger than those of AKARI objects with counterparts. These results suggest that our sample includes luminous, heavily dust-obscured SFGs/AGNs at $z\sim 1-4$ that are missed by previous optical surveys, providing very interesting targets for the coming James Webb Space Telescope era.

We report a discovery of $z = 4.72$ radio galaxy, HSC J083913.17+011308.1, by using the Lyman break technique with the Hyper Suprime-Cam Subaru Strategic Survey (HSC-SSP) catalog for VLA FIRST radio sources. The number of known high-$z$ radio galaxies (HzRGs) at $z > 3$ is quite small to constrain the evolution of HzRGs so far. The deep and wide-area optical survey by HSC-SSP enables us to apply the Lyman break technique to a large search for HzRGs. For an HzRG candidate among pre-selected $r$-band dropouts with a radio detection, a follow-up optical spectroscopy with GMOS/Gemini has been performed. The obtained spectrum presents a clear Ly$\alpha$ emission line redshifted to $z=4.72$. The SED fitting analysis with the rest-frame UV and optical photometries suggests the massive nature of this HzRG with $\log{M_*/M_{\odot}} = 11.4$. The small equivalent width of Ly$\alpha$ and the moderately red UV colors indicate its dusty host galaxy, implying a chemically evolved and dusty system. The radio spectral index does not meet a criterion for an ultra-steep spectrum: $\alpha^{325}_{1400}$ of $-1.1$ and $\alpha^{150}_{1400}$ of $-0.9$, demonstrating that the HSC-SSP survey compensates for a sub-population of HzRGs which are missed in surveys focusing on an ultra-steep spectral index.

We analyze an extremely deep 450-$\mu$m image ($1\sigma=0.56$\,mJy\,beam$^{-1}$) of a $\simeq 300$\,arcmin$^{2}$ area in the CANDELS/COSMOS field as part of the SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES). We select a robust (signal-to-noise ratio $\geqslant 4$) and flux-limited ($\geqslant 4$\,mJy) sample of 164 sub-millimeter galaxies (SMGs) at 450-$\mu$m that have $K$-band counterparts in the COSMOS2015 catalog identified from radio or mid-infrared imaging. Utilizing this SMG sample and the 4705 $K$-band-selected non-SMGs that reside within the noise level $\leqslant 1$\,mJy\,beam$^{-1}$ region of the 450-$\mu$m image as a training set, we develop a machine-learning classifier using $K$-band magnitude and color-color pairs based on the thirteen-band photometry available in this field. We apply the trained machine-learning classifier to the wider COSMOS field (1.6\,deg$^{2}$) using the same COSMOS2015 catalog and identify a sample of 6182 450-$\mu$m SMG candidates with similar colors. The number density, radio and/or mid-infrared detection rates, redshift and stellar mass distributions, and the stacked 450-$\mu$m fluxes of these SMG candidates, from the S2COSMOS observations of the wide field, agree with the measurements made in the much smaller CANDELS field, supporting the effectiveness of the classifier. Using this 450-$\mu$m SMG candidate sample, we measure the two-point autocorrelation functions from $z=3$ down to $z=0.5$. We find that the 450-$\mu$m SMG candidates reside in halos with masses of $\simeq (2.0\pm0.5) \times10^{13}\,h^{-1}\,\rm M_{\odot}$ across this redshift range. We do not find evidence of downsizing that has been suggested by other recent observational studies.

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$ $<$ 20$^{\circ}$). The intrinsic 2-10 keV luminosities are $3.9^{+2.2}_{-1.7}$ $\times$ $10^{43}$ erg s$^{-1}$ (UGC 2608) and $2.0^{+3.3}_{-1.0}$ $\times$ $10^{43}$ erg s$^{-1}$ (NGC 5135). The [O IV]-to-nuclear-12 $\mu$m luminosity ratios are larger than those of typical Seyferts, which are consistent with the torus covering factors ($C_{\rm T} \lesssim$ 0.7) estimated from the torus angular widths and column densities by X-ray spectroscopy. The torus covering factors and Eddington ratios ($\lambda_{\rm Edd} \sim$ 0.1) follow the relation found by Ricci et al. (2017c) for local AGNs, implying that their tori become geometrically thin due to significant radiation pressure of the AGN that blows out some part of the tori. These results indicate that the CT AGNs in these "non-merger" LIRGs are just a normal AGN population seen edge-on through a large line-of-sight column density. They are in contrast to the buried CT AGNs in late-stage mergers that have large torus covering factors even at large Eddington ratios.

We investigated the properties of AGN environments, particularly environments where the association of luminous galaxies (LGs) is found within 4~Mpc from AGNs with redshifts of 0.8 -- 1.1. For comparison, three additional AGN environments, (namely, AGNs of all types, type~1 AGNs with X-ray and/or radio detection, and type~2 AGNs) and an environment of blue $M_{*}$, characteristic luminosity of the Schechter function, galaxies were investigated. The cross-correlation function with the surrounding galaxies was measured and compared between the AGN and blue galaxy samples. We also compared the distributions of color, absolute magnitude, and stellar mass of the galaxies around such target objects. The properties of clusters detected using surrounding galaxies selected based on a photometric redshift were examined and compared for different samples. The target AGNs were drawn from the Million Quasars (MILLIQUAS) catalog, and the blue galaxies were drawn from six redshift survey catalogs (SDSS, WiggleZ, DEEP2, VVDS, VIPERS, and PRIMUS). The galaxies used as a measure of the environment around the targets are drawn from S18a internal data released by the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). We found that, among the five AGN and blue galaxy samples considered, the environment of AGN-LG pairs is the most enriched with luminous galaxies. We also found an enhancement in the number of mass-selected clusters in the AGN-LG pair sample against those in the other samples. The results obtained in this study indicate that existence of multiple clusters is the major driver in the association of AGNs and LGs, rather than a single large-mass dark matter halo hosting the AGN.

We construct a SCUBA-2 450-$\mu$m map in the COSMOS field that covers an area of 300 arcmin$^{2}$ and reaches a 1$\sigma$ noise level of 0.65 mJy in the deepest region. We extract 256 sources detected at 450 $\mu$m with signal-to-noise ratio $>$ 4.0 and analyze the physical properties of their multi-wavelength counterparts. We find that most of the sources are at $z\lesssim3$, with a median of $z = 1.79^{+0.03}_{-0.15}$. About $35^{+32}_{-25}$% of our sources are classified as starburst galaxies based on their total star-formation rates (SFRs) and stellar masses ($M_{\ast}$). By fitting the far-infrared spectral energy distributions, we find that our 450-$\mu$m-selected sample has a wide range of dust temperatures (20 K $ \lesssim T_{\rm d} \lesssim$ 60 K), with a median of ${T}_{\rm d} = 38.3^{+0.4}_{-0.9}$ K. We do not find a redshift evolution in dust temperature for sources with $L_{\rm IR}$ > $10^{12}$ $\rm L_\odot$ at $z<3$. However, we find a moderate correlation where dust temperature increases with the deviation from the SFR-$M_{\ast}$ relation. The increase in dust temperature also correlates with optical morphology, which is consistent with merger-triggered starbursts in sub-millimeter galaxies. Our galaxies do not show the tight IRX-$\beta_{\rm UV}$ correlation that has been observed in the local Universe. We construct the infrared luminosity functions of our 450-$\mu$m sources and measure their comoving SFR densities. The contribution of the $L_{\rm IR}$ > $10^{12}$ $\rm L_\odot$ population to the SFR density rises dramatically from $z$ = 0 to 2 ($\propto$ ($1+z$)$^{3.9\pm1.1}$) and dominates the total SFR density at $z \gtrsim 2$.

We investigate the evolution of the gas mass fraction for galaxies in the COSMOS field using submillimetre emission from dust at 850$\mu$m. We use stacking methodologies on the 850$\mu$m S2COSMOS map to derive the gas mass fraction of galaxies out to high redshifts, 0 <= $z$ <= 5, for galaxies with stellar masses of $10^{9.5} < M_* (\rm M_{\odot}) < 10^{11.75}$. In comparison to previous literature studies we extend to higher redshifts, include more normal star-forming galaxies (on the main sequence), and also investigate the evolution of the gas mass fraction split by star-forming and passive galaxy populations. We find our stacking results broadly agree with scaling relations in the literature. We find tentative evidence for a peak in the gas mass fraction of galaxies at around $z$ ~ 2.5-3, just before the peak of the star formation history of the Universe. We find that passive galaxies are particularly devoid of gas, compared to the star-forming population. We find that even at high redshifts, high stellar mass galaxies still contain significant amounts of gas.

We present far-infrared (FIR) properties of an extremely luminous infrared galaxy (ELIRG) at $z_{\rm spec}$ = 3.703, WISE J101326.25+611220.1 (WISE1013+6112). This ELIRG is selected as an IR-bright dust-obscured galaxy (DOG) based on the photometry from the Sloan digital sky survey (SDSS) and wide-field infrared survey explorer (WISE). In order to derive its accurate IR luminosity, we perform follow-up observations at 89 and 154 $\mu$m using the high-resolution airborne wideband camera-plus (HAWC+) on board the 2.7-m stratospheric observatory for infrared astronomy (SOFIA) telescope. We conduct spectral energy distribution (SED) fitting with CIGALE using 15 photometric data (0.4-1300 $\mu$m). We successfully pin down FIR SED of WISE1013+6112 and its IR luminosity is estimated to be $L_{\rm IR}$ = (1.62 $\pm$ 0.08) $\times 10^{14}$ $L_{\odot}$, making it one of the most luminous IR galaxies in the universe. We determine the dust temperature of WISE1013+6112 is $T_{\rm dust}$ = 89 $\pm$ 3 K, which is significantly higher than that of other populations such as SMGs and FIR-selected galaxies at similar IR luminosities. The resultant dust mass is $M_{\rm dust}$ = (2.2 $\pm$ 0.1) $\times 10^{8}$ $M_{\odot}$. This indicates that WISE1013+6112 has a significant active galactic nucleus (AGN) and star-forming activity behind a large amount of dust.

We report the discovery of a Compton-thick (CT) dust-obscured galaxy (DOG) at $z$ = 0.89, WISE J082501.48+300257.2 (WISE0825+3002), observed by Nuclear Spectroscopic Telescope Array (NuSTAR). X-ray analysis with the XCLUMPY model revealed that hard X-ray luminosity in the rest-frame 2-10 keV band of WISE0825+3002 is $L_{\rm X}$ (2-10 keV) = $4.2^{+2.8}_{-1.6} \times 10^{44}$ erg s$^{-1}$ while its hydrogen column density is $N_{\rm H}$ = $1.0^{+0.8}_{-0.4} \times 10^{24}$ cm$^{-2}$, indicating that WISE0825+3002 is a mildly CT active galactic nucleus (AGN). We performed the spectral energy distribution (SED) fitting with CIGALE to derive its stellar mass, star formation rate, and infrared luminosity. The estimated Eddington ratio based on stellar mass and integration of the best-fit SED of AGN component is $\lambda_{\rm Edd}$ = 0.70, which suggests that WISE0825+3002 harbors an actively growing black hole behind a large amount of gas and dust. We found that the relationship between luminosity ratio of X-ray and 6 $\mu$m, and Eddington ratio follows an empirical relation for AGNs reported by Toba et al. (2019a).

The properties of low-mass galaxies hosting central black holes provide clues about the formation and evolution of the progenitors of supermassive black holes. In this letter, we present HSC-XD 52, a spectroscopically confirmed low-mass active galactic nucleus (AGN) at an intermediate redshift of $z\sim0.56$. We detect this object as a very luminous X-ray source coincident with a galaxy observed by the Hyper Suprime-Cam (HSC) as part of a broader search for low-mass AGN. We constrain its stellar mass through spectral energy distribution modeling to be LMC-like at $M_\star \approx 3 \times 10^9 M_\odot$, placing it in the dwarf regime. We estimate a central black hole mass of $M_\mathrm{BH} \sim 10^{6} M_\odot$. With an average X-ray luminosity of $L_X \approx 3.5 \times 10^{43}~\mathrm{erg}~\mathrm{s}^{-1}$, HSC-XD 52 is among the most luminous X-ray selected AGN in dwarf galaxies. The spectroscopic and photometric properties of HSC-XD 52 indicate that it is an intermediate redshift counterpart to local low-mass AGN.

We report on the torus constraints of the Compton-thick AGN with double-peaked optical narrow line region (NLR) emission lines, ANEPD-CXO245, at z=0.449 in the AKARI NEP Deep Field. The unique infrared data on this field, including those from the nine-band photometry over 2-24 $\mu$m with the AKARI Infrared Camera (IRC), and the X-ray spectrum from Chandra allow us to constrain torus parameters such as the torus optical depth, X-ray absorbing column, torus angular width ($\sigma$) and viewing angle ($i$). We analyze the X-ray spectrum as well as the UV-optical-infrared spectral energy distribution (UOI-SED) with clumpy torus models in X-ray (XCLUMPY; Tanimoto et al. 2019) and infrared (CLUMPY; Nenkova et al. 2008) respectively. From our current data, the constraints on $\sigma$--$i$ from both X-rays and UOI show that the line of sight crosses the torus as expected for a type 2 AGN. We obtain a small X-ray scattering fraction (<0.1%), which suggests narrow torus openings, giving preference to the bi-polar outflow picture of the double-peaked profile. Comparing the optical depth of the torus from the UOI-SED and the absorbing column density $N_{\rm H}$ from the X-ray spectrum, we find that the gas-to-dust ratio is $\gtrsim 4$ times larger than the Galactic value.

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$\sigma$ 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.

Recent active galactic nucleus (AGN) and quasar surveys have revealed a population showing rapid AGN luminosity variability by a factor of $\sim10$. Here we present the most drastic AGN luminosity decline by a factor of $\gtrsim 10^{3}$ constrained by a NuSTAR X-ray observation of the nearby galaxy Arp 187, which is a promising "dead" quasar whose current activity seems quiet but whose past activity of $L_\mathrm{bol} \sim 10^{46}$ erg s$^{-1}$ is still observable at a large scale by its light echo. The obtained upper bound of the X-ray luminosity is $\log (L_{\rm 2-10 keV}/{\rm erg} {\rm s}^{-1}) < 41.2$, corresponding to $\log (L_\mathrm{bol}/{\rm erg} {\rm s}^{-1}) < 42.5$, indicating an inactive central engine. Even if a putative torus model with $N_\mathrm{H} \sim 1.5 \times 10^{24}$ cm$^{-2}$ is assumed, the strong upper-bound still holds with $\log (L_{\rm 2-10 keV}/{\rm erg} {\rm s}^{-1}) < 41.8$ or $\log (L_\mathrm{bol}/{\rm erg} {\rm s}^{-1}) < 43.1$. Given the expected size of the narrow line region, this luminosity decrease by a factor of $\gtrsim 10^3$ must have occurred within $\lesssim 10^4$ yr. This extremely rapid luminosity/accretion shutdown is puzzling and it requires one burst-like accretion mechanism producing a clear outer boundary for an accretion disk. We raise two possible scenarios realizing such an accretion mechanism: a mass accretion 1) by the tidal disruption of a molecular cloud and/or 2) by the gas depletion as a result of vigorous nuclear starformation after rapid mass inflow to the central engine.

This paper presents the second data release of the Hyper Suprime-Cam Subaru Strategic Program, a wide-field optical imaging survey on the 8.2 meter Subaru Telescope. The release includes data from 174 nights of observation through January 2018. The Wide layer data cover about 300 deg^2 in all five broadband filters (grizy) to the nominal survey exposure (10min in gr and 20min in izy). Partially observed areas are also included in the release; about 1100 deg^2 is observed in at least one filter and one exposure. The median seeing in the i-band is 0.6 arcsec, demonstrating the superb image quality of the survey. The Deep (26 deg^2) and UltraDeep (4 deg^2) data are jointly processed and the UltraDeep-COSMOS field reaches an unprecedented depth of i~28 at 5 sigma for point sources. In addition to the broad-bands, narrow-band data are also available in the Deep and UltraDeep fields. This release includes a major update to the processing pipeline, including improved sky subtraction, PSF modeling, object detection, and artifact rejection. The overall data quality has been improved, but this release is not without problems; there is a persistent deblender problem as well as new issues with masks around bright stars. The user is encouraged to review the issue list before utilizing the data for scientific explorations. All the image products as well as catalog products are available for download. The catalogs are also loaded to a database, which provides an easy interface for users to retrieve data for objects of interest. In addition to these main data products, detailed galaxy shape measurements withheld from the Public Data Release 1 (PDR1) are now available to the community. The shape catalog is drawn from the S16A internal release, which has a larger area than PDR1 (160 deg^2). All products are available at the data release site, https://hsc-release.mtk.nao.ac.jp/.

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.

We present deep near-infrared spectroscopy of six quasars at 6.1<z<6.7 with VLT/X-Shooter and Gemini-N/GNIRS. Our objects, originally discovered through a wide-field optical survey with the Hyper Suprime-Cam (HSC) Subaru Strategic Program (HSC-SSP), have the lowest luminosities (-25.5< M1450<-23.1 mag) of the z>5.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.

Context. Although heavily obscured active galactic nuclei (AGNs) have been found by many observational studies, the properties of the surrounding dust are poorly understood. Using AKARI/IRC spectroscopy, we discover a new sample of a heavily obscured AGN in LEDA 1712304 which shows a deep spectral absorption feature due to silicate dust. Aims. We study the infrared (IR) spectral properties of circumnuclear silicate dust in LEDA 1712304. Methods. We perform IR spectral fitting, considering silicate dust properties such as composition, porosity, size and crystallinity. Spectral energy distribution (SED) fitting is also performed to the flux densities in the UV to sub-millimeter range to investigate the global spectral properties. Results. The best-fit model indicates 0.1 $\mu$m-sized porous amorphous olivine (${\rm Mg_{2x}Fe_{2-2x}SiO_4}$; $x=0.4$) with $4\%$ crystalline pyroxene. The optical depth is $\tau_{\rm sil}{\sim}2.3$, while the total IR luminosity and stellar mass are estimated to be $L_{\rm IR}=(5\pm1){\times}10^{10}\,L_{\odot}$ and $M_{\rm star}=(2.7\pm0.8){\times}10^{9}\,M_{\odot}$, respectively. In such low $L_{\rm IR}$ and $M_{\rm star}$ ranges, there are few galaxies which show that large ${\tau}_{\rm sil}$. Conclusions. The silicate dust in the AGN torus of LEDA 1712304 has properties notably similar to those in other AGNs as a whole, but slightly different in the wing shape of the absorption profile. The porosity of the silicate dust suggests dust coagulation or processing in the circumnuclear environments, while the crystallinity suggests that the silicate dust is relatively fresh.

We present physical properties of radio galaxies (RGs) with $f_{\rm 1.4 GHz} >$ 1 mJy discovered by Subaru Hyper Supreme-Cam (HSC) and VLA Faint Images of the Radio Sky at Twenty-Centimeters (FIRST) survey. For 1056 FIRST RGs at $0 < z \leq 1.7$ with HSC counterparts in about 100 deg$^2$, we compiled multi-wavelength data of optical, near-infrared (IR), mid-IR, far-IR, and radio (150 MHz). We derived their color excess ($E (B-V)_{*}$), stellar mass, star formation rate (SFR), IR luminosity, the ratio of IR and radio luminosity ($q_{\rm IR}$), and radio spectral index ($\alpha_{\rm radio}$) that are derived from the SED fitting with CIGALE. We also estimated Eddington ratio based on stellar mass and integration of the best-fit SEDs of AGN component. We found that $E (B-V)_{*}$, SFR, and IR luminosity clearly depend on redshift while stellar mass, $q_{\rm IR}$, and $\alpha_{\rm radio}$ do not significantly depend on redshift. Since optically-faint ($i_{\rm AB} \geq 21.3$) RGs that are newly discovered by our RG survey tend to be high redshift, they tend to not only have a large dust extinction and low stellar mass but also have high SFR and AGN luminosity, high IR luminosity, and high Eddington ratio compared to optically-bright ones. The physical properties of a fraction of RGs in our sample seem to differ from a classical view of RGs with massive stellar mass, low SFR, and low Eddington ratio, demonstrating that our RG survey with HSC and FIRST provides us curious RGs among entire RG population.

We report on the optical properties of infrared (IR)-bright dust-obscured galaxies (DOGs) that are defined as (i-[22])_AB > 7.0. Because supermassive black holes (SMBHs) in IR-bright DOGs are expected to be rapidly growing in the major merger scenario, they provide useful clues for understanding the co-evolution of SMBHs and their host galaxies. However, the optical properties of IR-bright DOGs remain unclear because the optical emission of a DOG is very faint. By combining ~ 105 deg^2 images of the optical, near-IR, and mid-IR data obtained from the Subaru Hyper Suprime-Cam (HSC) survey, the VISTA VIKING survey, and the WISE all-sky survey, respectively, 571 IR-bright DOGs were selected. We found 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 later stage in the major merger scenario. We discovered eight IR-bright DOGs showing a significant blue excess in blue HSC bands (BluDOGs). This blue excess can be interpreted as a leaked AGN emission that is either a directly leaking or a scattered AGN emission, as proposed for some blue-excess Hot DOGs in earlier studies.

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 ($<$0.5%) tend to show weaker [O IV]-to-12 $\mu$m luminosity ratios than the average of Seyfert 2 galaxies. This suggests that this ratio is a good indicator for identifying buried AGNs. Then, we apply this criterion to 23 local ultra/luminous infrared galaxies (U/LIRGs) in various merger stages hosting AGNs. We find that AGNs in most of mid- to late-stage mergers are buried, while those in earlier stage ones (including non-merger) are not. This result suggests that the fraction of buried AGNs in U/LIRGs increases as the galaxy-galaxy interaction becomes more significant.

Much of the star formation is obscured by dust. For the complete understanding of the cosmic star formation history (CSFH), infrared (IR) census is indispensable. AKARI carried out deep mid-infrared observations using its continuous 9-band filters in the North Ecliptic Pole (NEP) field (5.4 deg$^2$). This took significant amount of satellite's lifetime, $\sim$10\% of the entire pointed observations. By combining archival Spitzer (5 bands) and WISE (4 bands) mid-IR photometry, we have, in total, 18 band mid-IR photometry, which is the most comprehensive photometric coverage in mid-IR for thousands of galaxies. However previously, we only had shallow optical imaging ($\sim$25.9ABmag) in a small area of 1.0 deg$^2$. As a result, there remained thousands of AKARI's infrared sources undetected in optical. Using the new Hyper Suprime-Cam on Subaru telescope, we obtained deep enough optical images of the entire AKARI NEP field in 5 broad bands ($g\sim$27.5mag). These provided photometric redshift, and thereby IR luminosity for the previously undetected faint AKARI IR sources. Combined with the accurate mid-IR luminosity measurement, we constructed mid-IR LFs, and thereby performed a census of dust-obscured CSFH in the entire AKARI NEP field. We have measured restframe 8$\mu$m, 12$\mu$m luminosity functions (LFs), and estimated total infrared LFs at 0.35$<$z$<$2.2. Our results are consistent with our previous work, but with much reduced statistical errors thanks to the large area coverage of the new data. We have possibly witnessed the turnover of CSFH at $z\sim$2.

We report the discovery of a quasar at z = 7.07, which was selected from the deep multi-band imaging data collected by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. This quasar, HSC J124353.93+010038.5, has an order of magnitude lower luminosity than do the other known quasars at z > 7. The rest-frame ultraviolet absolute magnitude is M1450 = -24.13 +/- 0.08 mag and the bolometric luminosity is Lbol = (1.4 +/- 0.1) x 10^{46} erg/s. Its spectrum in the optical to near-infrared shows strong emission lines, and shows evidence for a fast gas outflow, as the C IV line is blueshifted and there is indication of broad absorption lines. The Mg II-based black hole mass is Mbh = (3.3 +/- 2.0) x 10^8 Msun, thus indicating a moderate mass accretion rate with an Eddington ratio 0.34 +/- 0.20. It is the first z > 7 quasar with sub-Eddington accretion, besides being the third most distant quasar, known to date. The luminosity and black hole mass are comparable to, or even lower than, those measured for the majority of low-z quasars discovered by the Sloan Digital Sky Survey, and thus this quasar likely represents a z > 7 counterpart to quasars commonly observed in the low-z universe.

In order to construct a sample of ultra-luminous infrared galaxies (ULIRGs, with infrared luminosity, $L_{\rm IR} > 10^{12}$ L$_{\odot}$) at 0.5 < z < 1, we are conducting an optical follow-up program for bright 90-$\mu$m FIR sources with a faint optical (i < 20) counterpart selected in the AKARI Far-Infrared Surveyor (FIS) Bright Source catalog (Ver.2). AKARI-FIS-V2 J0916248+073034, identified as a ULIRG at z = 0.49 in the spectroscopic follow-up observation, indicates signatures of an extremely strong outflow in its emission line profiles. Its [OIII] 5007\AA\ emission line shows FWHM of 1830 km s$^{-1}$ and velocity shift of -770 km s$^{-1}$ in relative to the stellar absorption lines. Furthermore, low-ionization [OII] 3726\AA\ 3729\AA\ doublet also shows large FWHM of 910 km s$^{-1}$ and velocity shift of -380 km s$^{-1}$. After the removal of an unresolved nuclear component, the long-slit spectroscopy 2D image possibly shows that the outflow extends to radius of 4 kpc. The mass outflow and energy ejection rates are estimated to be 500 M$_{\odot}$ yr$^{-1}$ and $4\times10^{44}$ erg s$^{-1}$, respectively, which imply that the outflow is among the most powerful ones observed in ULIRGs and QSOs at 0.3 < z < 1.6. The co-existence of the strong outflow and intense star formation (star formation rate of 990 M$_{\odot}$ yr$^{-1}$) indicates that the feedback of the strong outflow has not severely affect the star-forming region of the galaxy.

We revisit the correlation between the mid-infrared (6 $\mu$m) and hard X-ray (2--10 keV) luminosities of active galactic nuclei (AGNs) to understand the physics behind it. We construct an X-ray flux-limited sample of 571 type 1 AGNs with $f_{0.5-2.0 \,{\rm keV}} > 2.4 \times 10^{-12}$ erg cm$^{-2}$ s$^{-1}$, drawn from the ROSAT Bright Survey catalog. Cross-matching the sample with infrared data taken from Wide-field Infrared Survey Explorer, we investigate the relation between the rest-frame 6 $\mu$m luminosity ($L_{\rm 6}$) and the rest-frame 2--10 keV luminosity ($L_{\rm X}$), where $L_{\rm 6}$ is corrected for the contamination of host galaxies by using the spectral energy distribution fitting technique. We confirm that $L_{\rm 6}$ and $L_{\rm X}$ are correlated over four orders of magnitude, in the range of $L_{\rm X} = 10^{42-46}$ erg s$^{-1}$. We investigate what kinds of physical parameters regulate this correlation. We find that $L_{\rm X}$/$L_{\rm 6}$ clearly depends on the Eddington ratio ($\lambda_{\rm Edd}$) as $\log \lambda_{\rm Edd} = -(0.56 \pm 0.10) \log \, (L_{\rm X}/L_{\rm 6}) - (1.07 \pm 0.05)$, even taking into account quasars that are undetected by ROSAT as well as those detected by XMM-Newton in the literature. We also add hyper-luminous quasars with $L_{\rm 6}$ $>$ 10$^{46}$ erg s$^{-1}$ in the literature and perform a correlation analysis. The resultant correlation coefficient is $-0.41 \pm 0.07$, indicating a moderately tight correlation between $L_{\rm X}$/$L_{\rm 6}$ and $\lambda_{\rm Edd}$. This means that AGNs with high Eddington ratios tend to have lower X-ray luminosities with respect to the mid-infrared luminosities. This dependence can be interpreted as a change in the structure of the accretion flow.

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.

Arp 187 is one of the fading active galactic nuclei (AGN), whose AGN activity is currently decreasing in luminosity. We investigate the observational signatures of AGN in Arp 187, which trace various physical scales from less than 0.1 pc to the nearly 10 kpc, to estimate the longterm luminosity change over $10^{4}$ years. The VLA 5 GHz, 8 GHz, and the ALMA 133 GHz images reveal bimodal jet lobes with $\sim$5 kpc size and the absence of the central radio-core. The 6dF optical spectrum shows that Arp 187 hosts narrow line region with the estimated size of $\sim$1 kpc, and the line strengths give the AGN luminosity of $L_{\rm bol}=1.5 \times 10^{46}$ erg s$^{-1}$. On the other hand, the current AGN activity estimated from the AGN torus emission gives the upper bound of $L_{\rm bol} < 2.2 \times 10^{43}$ erg s$^{-1}$. The absence of the radio-core gives the more strict upper bound of the current AGN luminosity of $L_{\rm bol} < 8.0 \times 10^{40}$ erg s$^{-1}$, suggesting that the central engine is already quenched. These multi-wavelength signatures indicate that Arp 187 hosts a "dying" AGN: the central engine is already dead, but the large scale AGN indicators are still observable as the remnant of the past AGN activity. The central engine has experienced the drastic luminosity decline by a factor of $\sim10^{3-5}$ fainter within $\sim10^{4}$ years, which is roughly consistent with the viscous timescale of the inner part of the accretion disk within $\sim$500 years.

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.

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 exceed 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 that the fitting linear function of the 1.4 GHz source counts has a slope that is flatter than that of the bright RGs, while optically faint radio quasars show a slope steeper than that of bright radio quasars. The optically faint RGs show a flat slope in the i-band number counts down to 24 mag, implying either less massive or distant radio-active galactic nuclei (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 redshift z > 1. This study gives ~1500 radio AGNs lying at the optically faint end and high-redshift regime not probed by previous searches.

We present structural parameters and morphological properties of faint 450-um selected submillimeter galaxies (SMGs) from the JCMT Large Program, STUDIES, in the COSMOS-CANDELS region. Their properties are compared to an 850um selected and a matched star-forming samples. We investigate stellar structures of 169 faint 450-um sources (S450=2.8-29.6mJy; S/N>4) at z<3 using HST near-infrared observations. Based on our spectral energy distribution fitting, half of such faint SMGs (LIR=10^11.65+-0.98Lsun) lie above the star-formation rate (SFR)/stellar mass plane. The size-mass relation shows that these SMGs are generally similar to less-luminous star-forming galaxies selected by NUV-r vs. r-J colors. Because of the intrinsic luminosity of the sample, their rest-frame optical emission is less extended than the 850um sources (S850>2mJy), and more extended than the star-forming galaxies in the same redshift range. For the stellar mass and SFR matched sample at z~=1 and z~=2, the size differences are marginal between faint SMGs and the matched galaxies. Moreover, faint SMGs have similar Sersic indices and projected axis ratios as star-forming galaxies with the same stellar mass and SFR. Both SMGs and the matched galaxies show high fractions (~70%) of disturbed features at z~=2, and the fractions depend on the SFRs. These suggest that their star formation activity is related to galaxy merging, and the stellar structures of SMGs are similar to those of star-forming galaxies. We show that the depths of submillimeter surveys are approaching the lower luminosity end of star-forming galaxies, allowing us to detect galaxies on the main sequence.

Deep MIR surveys have revealed numerous strongly star-forming galaxies at redshift z~<2. Their MIR fluxes are produced by a combination of continuum and PAH emission features. The PAH features can dominate the total MIR flux, but are difficult to measure without spectroscopy. We aim to study star-forming galaxies by using a blind spectroscopic survey at MIR wavelengths to understand evolution of their star formation rate (SFR) and specific SFR up to z~=0.5, by paying particular attention to their PAH properties. We conducted a low-resolution (R~=50) slitless spectroscopic survey at 5-13um of 9um flux-selected sources (>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 z<0.5. The rest-frame optical-MIR spectral energy distributions (SEDs) based on CFHT and AKARI/IRC imaging covering 0.37-18um were produced, and analysed in conjunction with the PAH spectroscopy. We defined the PAH enhancement by using the luminosity ratio of the 7.7um PAH feature over the 3.5um stellar component of the SEDs. The rest-frame SEDs of all PAH galaxies have a universal shape with stellar and 7.7um bumps, except that the PAH enhancement significantly varies as a function of the PAH luminosities. We identified a PAH-enhanced population at z~>0.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.

We report discovery of 41 new high-z quasars and luminous galaxies, which were spectroscopically identified at 5.7 < z < 6.9. This is the fourth in a series of papers from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, based on the deep multi-band imaging data collected by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. We selected the photometric candidates by a Bayesian probabilistic algorithm, and then carried out follow-up spectroscopy with the Gran Telescopio Canarias and the Subaru Telescope. Combined with the sample presented in the previous papers, we have now spectroscopically identified 137 extremely-red HSC sources over about 650 deg2, which include 64 high-z quasars, 24 high-z luminous galaxies, 6 [O III] emitters at z ~ 0.8, and 43 Galactic cool dwarfs (low-mass stars and brown dwarfs). The new quasars span the luminosity range from M1450 ~ -26 to -22 mag, and continue to populate a few magnitude lower luminosities than have been probed by previous wide-field surveys. In a companion paper, we derive the quasar luminosity function at z ~ 6 over an unprecedentedly wide range of M1450 ~ -28 to -21 mag, exploiting the SHELLQs and other survey outcomes.

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.

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.

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.

Understanding infrared (IR) luminosity is fundamental to understanding the cosmic star formation history and AGN evolution. Japanese infrared satellite, AKARI, provided unique data sets to probe this both at low and high redshift; the AKARI all sky survey in 6 bands (9-160 $\mu$m), and the AKARI NEP survey in 9 bands (2-24$\mu$m). The AKARI performed all sky survey in 6 IR bands (9, 18, 65, 90, 140, and 160 $\mu$m) with 3-10 times better sensitivity than IRAS, covering the crucial far-IR wavelengths across the peak of the dust emission. Combined with a better spatial resolution, we measure the total infrared luminosity ($L_{TIR}$) of individual galaxies, and thus, the total infrared luminosity density of the local Universe much more precisely than previous work. In the AKARI NEP wide field, AKARI has obtained deep images in the mid-infrared (IR), covering 5.4 deg$^2$. However, our previous work was limited to the central area of 0.25 deg$^2$ due to the lack of deep optical coverage. To rectify the situation, we used the newly advent Subaru telescope's Hyper Suprime-Cam to obtain deep optical images over the entire 5.4 deg$^2$ of the AKARI NEP wide field. With this deep and wide optical data, we, for the first time, can use the entire AKARI NEP wide data to construct restframe 8$\mu$m, 12$\mu$m, and total infrared (TIR) luminosity functions (LFs) at 0.15$<z<$2.2. A continuous 9-band filter coverage in the mid-IR wavelength (2.4, 3.2, 4.1, 7, 9, 11, 15, 18, and 24$\mu$m) by the AKARI satellite allowed us to estimate restframe 8$\mu$m and 12$\mu$m luminosities without using a large extrapolation based on a SED fit, which was the largest uncertainty in previous work. By combining these two results, we reveal dust-hidden cosmic star formation history and AGN evolution from z=0 to z=2.2, all probed by the AKARI satellite.

The AKARI North Ecliptic Pole (NEP) survey consists of two survey projects: NEP-Deep (0.5 sq.deg) and NEP-Wide (5.4 sq.deg), providing with tens of thousands of galaxies. A continuous filter coverage in the mid-infrared wavelengths (7, 9, 11, 15, 18 and 24 $\mu$m) is unique to diagnose the contributions from dusty star-formation activity and AGNs. Here we present current status focused on the newly obtained optical images and near-future prospects with a new X-ray telescope. Hyper Suprime-Cam on Subaru telescope is a gigantic optical camera with huge Field of View (FoV). Thanks to the wide FoV, we successfully obtained deep optical images at g, r, i, z and Y-bands covering most of the NEP-Wide field. Using the deep optical images, we identified over 5000 optical counterparts of the mid-IR sources, presumably deeply obscured galaxies in NEP-Wide field. We also investigated properties of these infrared sources with SED-fitting. eROSITA, to be launched early 2018, is a new all-sky X-ray survey telescope, and expected to conduct ultra deep 2-10 keV imaging toward NEP. We expect unprecedentedly numerous Compton-thick AGN candidates when combined with the multi-wavelength data in NEP region.

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 present a study of the connection between black hole accretion, star formation, and galaxy morphology at z<=2.5. We focus on active galactic nuclei (AGNs) selected by their mid-IR power-law emission. By fitting optical to far-IR photometry with state-of-the-art spectral energy distribution (SED) techniques, we derive stellar masses, star formation rates, dust properties, and AGN contributions in galaxies over the whole COSMOS field. We find that obscured AGNs lie within or slightly above the star-forming sequence. We confirm our previous finding about compact host galaxies of obscured AGNs at z~1, and find that galaxies with 20-50% AGN contributions tend to have smaller sizes, by ~25-50%, compared to galaxies without AGNs. Furthermore, we find that a high merger fraction of up to 0.5 is appropriate for the most luminous (log (LIR/Lsun) ~ 12.5) AGN hosts and non-AGN galaxies, but not for the whole obscured AGN sample. Moreover, merger fraction depends on the total and star-forming infrared luminosity, rather than the decomposed AGN infrared luminosity. Our results suggest that major mergers are not the main driver of AGN activity, and therefore obscured AGNs might be triggered by internal mechanisms, such as secular processes, disk instabilities, and compaction in a particular evolutionary stage. We make the SED modeling results publicly available.

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.

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

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.

Aims. We report our finding of a high excitation emission line nebula associated with an Ultra Luminous X-ray source (ULX) at $z=$ 0.027, which we found in our Chandra observation of the AKARI North Ecliptic Pole (NEP) Deep Field. Methods. We present a Chandra X-ray and Gran Telescopio Canarias (GTC) optical spectral analysis of the ULX blob. We investigate the nature of the emission line nebula by using line ratio diagnostic diagrams, and its physical properties estimated from Spectral Energy Distribution (SED) fitting. Results. The optical spectrum of this ULX blob shows emission line ratios that are located on the borderlines between star-formation and Seyfert regimes in [OIII]/H$\beta$-[OI]/H$\alpha$, [OIII]/H$\beta$-[SII]/H$\alpha$ and [OIII]/H$\beta$-[OIII]/[OII] diagnostic diagrams. These are in contrast with those of a nearby blob observed with the same slit, which occupy the HII regimes. This result suggests that the ionization of the emission line nebula associated with the ULX is significantly contributed by energy input from the accretion power of the ULX, in addition to the star formation activity in the blob, suggesting the existence of an accretion disk in the ULX emitting UV radiation, or exerting shock waves.

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.

The identification of high-redshift massive galaxies with old stellar populations may pose challenges to some models of galaxy formation. However, to securely classify a galaxy as quiescent, it is necessary to exclude significant ongoing star formation, something that can be challenging to achieve at high redshift. In this letter, we analyse deep ALMA/870um and SCUBA-2/450um imaging of the claimed "post-starburst" galaxy ZF-20115 at z=3.717 that exhibits a strong Balmer break and absorption lines. The far-infrared imaging reveals a luminous starburst located 0.4+/-0.1 arcsec (~3kpc in projection) from the position of the rest-frame ultra-violet/optical emission, with an obscured star-formation rate of 100 Mo/yr. This star-forming component is undetected in the rest-frame ultraviolet but contributes significantly to the lower angular resolution photometry at restframe wavelengths >3500A, significantly complicating the determination of a reliable stellar mass. Importantly, in the presence of dust obscuration, strong Balmer features are not a unique signature of a post-starburst galaxy and are indeed frequently observed in infrared-luminous galaxies. We conclude that the ZF20015 system does not pose a challenge to current models of galaxy formation and that deep sub-/millimeter observations are a prerequisite for any claims of quiescence. The multi-wavelength observations of ZF20115 unveil a complex system with an intricate and spatially-varying star-formation history. ZF20115 demonstrates that understanding high-redshift obscured starbursts will only be possible with multi-wavelength studies that include high-resolution observations, available with the JWST, at mid-infrared wavelengths.

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.

We aim to search for hyperliminous infrared (IR) galaxies (HyLIRGs) with IR luminosity $L_{{\rm IR}}$ $>$ 10$^{13}$ $L_{\odot}$ by applying the selection method of Dust Obscured Galaxies (DOGs). They are spatially rare but could correspond to a maximum phase of cosmic star formation and/or active galactic nucleus (AGN) activity, hence they are a crucial population for understanding the star formation and mass assembly history of galaxies. Combining the optical and IR catalogs obtained from Sloan Digital Sky Survey (SDSS) and Wide-field Infrared Survey Explorer (WISE), we performed the extensive HyLIRGs survey; we selected 5,311 IR-bright DOGs with $i$ -- [22] $>$ 7.0 and flux at 22 $\mu$m $>$ 3.8 mJy in 14,555 deg$^2$, where $i$ and [22] are $i$-band and 22 $\mu$m AB magnitudes, respectively. Among them, 67 DOGs have reliable spectroscopic redshifts that enable us to estimate their total IR luminosity based on the SED fitting. Consequently, we successfully discovered 24 HyLIRGs among the 67 spectroscopically-confirmed DOGs. We found that (i) $i$ - [22] color of IR-bright DOGs correlates with the total IR luminosity and (ii) surface number density of HyLIRGs is $>$ 0.17 deg$^{-2}$. A high fraction ($\sim$ 73%) of IR-bright DOGs with $i$ - [22] $>$ 7.5 shows $L_{{\rm IR}}$ $>$ 10$^{13}$ $L_{\odot}$, and the DOGs criterion we adopted could be independently-effective against "W1W2-dropout method" based on four WISE bands, for searching hyper IR luminous populations of galaxies.

We report a new sample of obscured active galactic nuclei (AGNs) selected from the XMM serendipitous source and AKARI point-source catalogs. We match X-ray sources with infrared (18 and 90 micron) sources located at |b|>10 deg to create a sample consisting of 173 objects. Their optical classifications and absorption column densities measured by X-ray spectra are compiled and study efficient selection criteria to find obscured AGNs. We apply the criteria (1) X-ray hardness ratio defined by using the 2-4.5 keV and 4.5-12 keV bands >-0.1 and (2) EPIC-PN count rate (CR) in the 0.2-12 keV to infrared flux ratio CR/F90<0.1 or CR/F18<1, where F18 and F90 are infrared fluxes at 18 and 90 micron in Jy, respectively, to search for obscured AGNs. X-ray spectra of 48 candidates, for which no X-ray results have been published, are analyzed and X-ray evidence for the presence of obscured AGNs such as a convex shape X-ray spectrum indicative of absorption of NH~10^{22-24} cm^-2, a very flat continuum, or a strong Fe-K emission line with an equivalent width of >700 eV is found in 26 objects. Six among them are classified as Compton-thick AGNs, and four are represented by either Compton-thin or Compton-thick spectral models. The success rate of finding obscured AGNs combining our analysis and the literature is 92% if the 18 micron condition is used. Of the 26 objects, 4 are optically classified as an HII nucleus and are new "elusive AGNs" in which star formation activity likely overwhelms AGN emission in the optical and infrared bands.

We present infrared galaxy luminosity functions (LFs) in the AKARI North Ecliptic Pole (NEP) deep field using recently-obtained, wider CFHT optical/near-IR images. AKARI has obtained deep images in the mid-infrared (IR), covering 0.6 deg$^2$ of the NEP deep field. However, our previous work was limited to the central area of 0.25 deg$^2$ due to the lack of optical coverage of the full AKARI NEP survey. To rectify the situation, we recently obtained CFHT optical and near-IR images over the entire AKARI NEP deep field. These new CFHT images are used to derive accurate photometric redshifts, allowing us to fully exploit the whole AKARI NEP deep field. AKARI's deep, continuous filter coverage in the mid-IR wavelengths (2.4, 3.2, 4.1, 7, 9, 11, 15, 18, and 24$\mu$m) exists nowhere else, due to filter gaps of other space telescopes. It allows us to estimate restframe 8$\mu$m and 12$\mu$m luminosities without using a large extrapolation based on spectral energy distribution (SED) fitting, which was the largest uncertainty in previous studies. Total infrared luminosity (TIR) is also obtained more reliably due to the superior filter coverage. The resulting restframe 8$\mu$m, 12$\mu$m, and TIR LFs at $0.15<z<2.2$ are consistent with previous works, but with reduced uncertainties, especially at the high luminosity-end, thanks to the wide field coverage. In terms of cosmic infrared luminosity density ($\Omega_{\mathrm{IR}}$), we found that the $\Omega_{\mathrm{IR}}$ evolves as $\propto (1+z)^{4.2\pm 0.4}$.

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.

In this work, we investigate the dependence of the covering factor (CF) of active galactic nuclei (AGNs) (i) on the mid-infrared (MIR) luminosity and (ii) on the redshift. We constructed 12- and 22-micron luminosity functions (LFs) at 0.006 < z < 0.3 using the Wide-field Infrared Survey Explorer} (WISE) data. Combining the WISE catalog with the Sloan Digital Sky Survey (SDSS) spectroscopic data, we selected 223,982 galaxies at 12 micron and 25,721 galaxies at 22 micron for spectroscopic classification. We then identified 16,355 AGNs at 12 micron and 4,683 AGNs at 22 micron by their optical emission lines and cataloged classifications in the SDSS. Following that, we estimated the CF as the fraction of type 2 AGN in all AGNs whose MIR emissions are dominated by the active nucleus (not their host galaxies) based on their MIR colors. We found that (i) the CF decreased with increasing MIR luminosity, regardless of the choice of type 2 AGN classification criteria, and (ii) the CF did not change significantly with the redshift for z < 0.2. Furthermore, we carried out various tests to determine the influence of selection bias and confirmed similar dependences exist even when taking these uncertainties into account. The luminosity dependence of the CF can be explained by the receding torus model, but the "modified" receding torus model gives a slightly better fit, as suggested by Simpson.

We present the 9 and 18 micron luminosity functions (LFs) of galaxies at 0.006 < z < 0.8 (with an average redshift of ~ 0.04) using the AKARI mid-infrared all-sky survey catalog. We selected 243 galaxies at 9 micron and 255 galaxies at 18 micron from the Sloan Digital Sky Survey (SDSS) spectroscopy region. These galaxies were then classified by their optical emission lines, such as the line width of H_alpha or by their emission line ratios of [OIII]/H_beta and [NII]/H_alpha into five types: Type 1 active galactic nuclei (AGN) (Type 1); Type 2 AGN (Type 2); low-ionization narrow emission line galaxies (LINER); galaxies with both star formation and narrow-line AGN activity (composite galaxies); and star-forming galaxies (SF). We found that (i) the number density ratio of Type 2 to Type 1 AGNs is 1.73 +/- 0.36, which is larger than a result obtained from the optical LF and (ii) this ratio decreases with increasing 18 micron luminosity.

We present a new sample of active galactic nuclei (AGNs) identified using the catalog of the AKARI Mid-infrared(MIR) All-Sky Survey. Our MIR search has an advantage in detecting AGNs that are obscured at optical wavelengths due to extinction. We first selected AKARI 9micron excess sources with F(9micron)/F(K_S)>2 where K_S magnitudes were taken from the Two Micron All Sky Survey. We then obtained follow-up near-infrared spectroscopy with the AKARI/IRC, to confirm that the excess is caused by hot dust. We also obtained optical spectroscopy with the Kast Double Spectrograph on the Shane 3-m telescope at Lick Observatory. On the basis of on these observations, we detected hot dust with a characteristic temperature of ~500K in two luminous infrared galaxies. The hot dust is suspected to be associated with AGNs that exhibit their nonstellar activity not in the optical, but in the near- and mid-infrared bands, i.e., they harbor buried AGNs. The host galaxy stellar masses of 4-6 x 10^9 M_sun are small compared with the hosts in optically-selected AGN populations. These objects were missed by previous surveys, demonstrating the power of the AKARI MIR All-Sky Survey to widen AGN searches to include more heavily obscured objects. The existence of multiple dusty star clusters with massive stars cannot be completely ruled out with our current data.

We present Subaru observations of the newly discovered luminous quasar AKARI J1757+5907, which shows an absorption outflow in its spectrum. The absorption consists of 9 distinct troughs, and our analysis focuses on the troughs at ~ -1000$ km s^{-1} for which we can measure accurate column densities of He I*, Fe II and Mg II. We use photoionization models to constrain the ionization parameter, total hydrogen column density, and the number density of the outflowing gas. These constraints yield lower limits for the distance, mass flow rate and kinetic luminosity for the outflow of 3.7 kpc, 70 M_{sun} yr^{-1}, and 2.0 x 10^{43} ergs s^{-1}, respectively. Such mass flow rate value can contribute significantly to the metal enrichment of the intra-cluster medium. We find that this moderate velocity outflow is similar to those recently discovered in massive post-starburst galaxies. Finally, we describe the scientific potential of future observations targeting this object.