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.

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 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 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 an observed relationship between Eddington ratio ({\lambda}Edd) and optical narrow-emission-line ratio ([N II]{\lambda}6583/H{\alpha}) of X-ray-selected broad-line active galactic nuclei (AGN) at 0.6 < z < 1.7. We use 27 near-infrared spectra from the Fiber Multi-Object Spectrograph along with 26 sources from the literature. We show that the {\lambda}Edd and [N II]{\lambda}6583/H{\alpha} ratio at 0.6 < z < 1.7 exhibits a similar anti-correlation distribution of {\lambda}Edd-[N II]{\lambda}6583/H{\alpha} as has been found for local (<z> = 0.036), ultra-hard X-ray selected AGN. The observed distribution implies that there is a consistent relationship from local to z \sim 1.7 which corresponds from the present time to 4 Gyr old. Further study of high redshift low Eddington ratio AGN (log{\lambda}Edd < -2) is necessary to determine fully whether the {\lambda}Edd-[N II]{\lambda}6583/H{\alpha} anti-correlation still holds in high-redshift AGN at low Eddington ratios.

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.

Despite decreasing cosmic star formation rate density over the last 10 Gyr, the stellar mass ($M_{*}$) buildups in galaxies were still progressing during this epoch. About 50\% of the current $M_{*}$ density in the universe was built over the last $\sim 8.7$ Gyr. In this research, we investigated the stellar mass buildup and quenching of spatially resolved regions within massive disk galaxies over the last 10 Gyr. We apply the spectral energy distribution (SED) fitting method to SEDs of sub-galactic regions in galaxies to derive the spatially resolved distributions of SFR and $M_{*}$ in the galaxies. This namely \textit{pixel-to-pixel SED fitting} method is applied to massive disk galaxies at $0.01<z<0.02$ and $0.8<z<1.8$. We found that massive disk galaxies tend to build their $M_{*}$ and quench their star formation progressively from the central region to the outskirts, i.e. \textit{inside-out} stellar mass buildup and quenching.

In the SSA22 field which exhibits a large-scale proto-cluster at $z=3.1$, we carried out a spectroscopic survey for Lyman Break Galaxies (LBGs) with the VLT/VIMOS and identified 78 confident LBGs at $z=2.5$--4. We stacked their spectra in the observer's frame by using a sophisticated method. Analyzing the composite spectrum, we have revealed that the large-scale proto-cluster at $z=3.1$ has a strong HI absorption dip of rest-frame equivalent width of $-1.7$ A. Another strong absorption dip found at $z=3.28$ is associated with a modestly high-density LBG peak, similar to that at $z=3.1$. We have also detected an HI transparency peak at $z=2.98$ in the composite spectrum, coincident with a void in the LBG distribution. In this paper, we also investigated the relation between LBGs, HI gas and AGNs at $z=3$--4 in the SSA22 field. Two AGNs at $z=3.353$ and 3.801 are, respectively, associated with the LBG concentration of an overdensity factor $\delta_{\rm LBG}\simeq2$ in the present statistics. Another structure at $z=3.453$ is remarkable: 20 comoving Mpc scale dense HI gas which is not associated with any apparent LBG overdensity but involving a pair of AGNs. Such structure may be a new type of the AGN-matter correlation. If the inhomogeneous structures over a comoving Gpc scale found in this paper are confirmed with sufficient statistics in the future, the SSA22 field will become a key region to test the standard cold dark matter structure formation scenario.

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 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 near-IR spectroscopy in J- and H-band for a large sample of 243 X-ray selected moderate-luminosity type-1 AGN in the COSMOS, SXDS and E-CDF-S survey fields using the multi-object spectrograph Subaru/FMOS. Our sample covers the redshift range 0.5<z<3.0 and an X-ray luminosity range of $10^{43}\lesssim L_X \lesssim 10^{45}$~erg s$^{-1}$. We provide emission-line properties and derived virial black hole mass estimates, bolometric luminosities and Eddington ratios, based on H$\alpha$ (211), H$\beta$ (63) and MgII (4). We compare line widths, luminosities and black hole mass estimates from H$\alpha$ and H$\beta$ and augment these with commensurate measurements of MgII and CIV detected in optical spectra. We demonstrate the robustness of using H$\alpha$, H$\beta$ and MgII as reliable black hole mass estimators for high-z moderate-luminosity AGN, while the use of CIV is prone to large uncertainties (>0.4 dex). We extend a recently proposed correction based on the CIV blueshift to lower luminosities and black hole masses. While our sample shows an improvement in their CIV black hole mass estimates, the deficit of high blueshift sources reduces its overall importance for moderate-luminosity AGN, compared to the most luminous quasars. In addition, we revisit luminosity correlations between $L_{\rm{bol}}$, $L_X$, $L_{\rm{[OIII]}}$, $L_{5100}$ and $L_{\rm{H}\alpha}$ and find them to be consistent with a simple empirical model, based on a small number of well-established scaling relations. We finally highlight our highest redshift AGN, CID 781, at z=4.6 which shows the lowest black hole mass ($\sim10^8$M$_\odot$) among current near-IR samples at this redshift, and is in a state of fast growth.

We investigate a relation between surface densities of star formation rate (SFR) and stellar mass ($M_{*}$) at a $\sim 1$ kpc scale namely spatially resolved star formation main sequence (SFMS) in massive ($\log(M_{*}/M_{\odot})>10.5$) face-on disc galaxies at $0.01<z<0.02$ and $0.8<z<1.8$ and examine evolution of the relation. The spatially resolved SFMS of $z\sim 0$ galaxies is discussed in a companion paper. For $z\sim 1$ sample, we use 8 bands imaging dataset from CANDELS and 3D-HST and perform a pixel-to-pixel SED fitting to derive the spatially resolved SFR and $M_{*}$. We find a linear spatially resolved SFMS in the $z\sim 1$ galaxies that lie on the global SFMS, while a 'flattening' at high $\Sigma_{*}$ end is found in that relation for the galaxies that lie below the global SFMS. Comparison with the spatially resolved SFMS of the $z\sim 0$ galaxies shows smaller difference in the specific SFR (sSFR) at low $\Sigma_{*}$ than that at high $\Sigma_{*}$. This trend is consistent with the evolution of the sSFR$(r)$ radial profile, which shows a faster decrease in the central region than in the outskirt, agrees with the inside-out quenching scenario. We then derive an empirical model for the evolution of the $\Sigma_{*}(r)$, $\Sigma_{\rm SFR}(r)$ and sSFR$(r)$ radial profiles. Based on the empirical model, we estimate the radial profile of the quenching timescale and reproduce the observed spatially resolved SFMS at $z\sim 1$ and $z\sim 0$.

The Gemini Infrared Multi-Object Spectrograph (GIRMOS) is a powerful new instrument being built to facility-class standards for the Gemini telescope. It takes advantage of the latest developments in adaptive optics and integral field spectrographs. GIRMOS will carry out simultaneous high-angular-resolution, spatially-resolved infrared ($1-2.4$ $\mu$m) spectroscopy of four objects within a two-arcminute field-of-regard by taking advantage of multi-object adaptive optics. This capability does not currently exist anywhere in the world and therefore offers significant scientific gains over a very broad range of topics in astronomical research. For example, current programs for high redshift galaxies are pushing the limits of what is possible with infrared spectroscopy at $8-10$-meter class facilities by requiring up to several nights of observing time per target. Therefore, the observation of multiple objects simultaneously with adaptive optics is absolutely necessary to make effective use of telescope time and obtain statistically significant samples for high redshift science. With an expected commissioning date of 2023, GIRMOS's capabilities will also make it a key followup instrument for the James Webb Space Telescope when it is launched in 2021, as well as a true scientific and technical pathfinder for future Thirty Meter Telescope (TMT) multi-object spectroscopic instrumentation. In this paper, we will present an overview of this instrument's capabilities and overall architecture. We also highlight how this instrument lays the ground work for a future TMT early-light instrument.

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.

We investigate the X-ray active galactic nucleus (AGN) properties of millimeter galaxies in the Great Observatories Origins Deep Survey South (GOODS-S) field detected with the Atacama Large Millimeter/submillimeter Array (ALMA), by utilizing the Chandra 7-Ms data, the deepest X-ray survey to date. Our millimeter galaxy sample comes from the ASAGAO survey covering 26 arcmin$^2$ (12 sources at a 1.2-mm flux-density limit of $\approx$0.6 mJy), supplemented by the deeper but narrower 1.3-mm survey of a part of the ASAGAO field by Dunlop et al.\ (2017). Fourteen out of the total 25 millimeter galaxies have Chandra counterparts. The observed AGN fractions at $z=1.5-3$ is found to be 90$^{+8}_{-19}$\% and $57^{+23}_{-25}$\% for the ultra/luminous infrared galaxies with $\log L_{\rm IR}/L_{\odot} = 12-12.8$ and $\log L_{\rm IR}/L_{\odot} = 11.5-12$, respectively. The majority ($\sim$2/3) of the ALMA and/or Herschel detected X-ray AGNs at $z=1.5-3$ appear to be star-formation dominant populations, having $L_{\rm X}$/ $L_{\rm IR}$ ratios smaller than the "simultaneous evolution" value expected from the local black-hole mass to stellar mass ($M_{\rm BH}$-$M_*$) relation. On the basis of the $L_{\rm X}$ and stellar mass relation, we infer that a large fraction of star-forming galaxies at $z=1.5-3$ have black hole masses smaller than those expected from the local $M_{\rm BH}$-$M_*$ relation. This is opposite to previous reports on luminous AGN at same redshifts detected in wider and shallower surveys, which are subject to selection biases against lower luminosity AGN. Our results are consistent with an evolutionary scenario that star formation occurs first, and an AGN-dominant phase follows later, in objects finally evolving into galaxies with classical bulges.

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

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 conduct a systematic search for galaxy protoclusters at $z\sim3.8$ based on the latest internal data release (S16A) of the Hyper SuprimeCam Subaru strategic program (HSC-SSP). In the Wide layer of the HSC-SSP, we investigate the large-scale projected sky distribution of $g$-dropout galaxies over an area of $121\,\mathrm{deg^2}$, and identify 216 large-scale overdense regions ($>4\sigma$ overdensity significance) that are good protocluster candidates. Of these, 37 are located within $8\,\mathrm{arcmin}$ ($3.4\,\mathrm{physicalMpc}$) from other protocluster candidates of higher overdensity, and are expected to merge into a single massive structure by $z=0$. Therefore, we find 179 unique protocluster candidates in our survey. A cosmological simulation that includes projection effects predicts that more than 76\% of these candidates will evolve into galaxy clusters with halo masses of at least $10^{14}\,M_{\odot}$ by $z=0$. The unprecedented size of our protocluster candidate catalog allowed us to perform, for the first time, an angular clustering analysis of the systematic sample of protocluster candidates. We find a correlation length of $35.0\,h^{-1}\,\mathrm{Mpc}$. The relation between correlation length and number density of $z\sim3.8$ protocluster candidates is consistent with the prediction of the $\Lambda$CDM model, and the correlation length is similar to that of rich clusters in the local universe. This result suggests that our protocluster candidates are tracing similar spatial structures as those expected of the progenitors of rich clusters and enhances the confidence that our method to identify protoclusters at high redshifts is robust. In the coming years, our protocluster search will be extended to the entire HSC-SSP Wide sky coverage of $\sim1400\,\mathrm{deg^2}$ to probe cluster formation over a wide redshift range of $z\sim2\mathrm{-}6$.

We study the UV luminosity functions (LFs) at $z\sim 4$, $5$, $6,$ and $7$ based on the deep large-area optical images taken by the Hyper Suprime-Cam (HSC) Subaru strategic program (SSP). On the 100 deg$^2$ sky of the HSC SSP data available to date, we make enormous samples consisting of a total of 579,565 dropout candidates at $z\sim 4-7$ by the standard color selection technique, 358 out of which are spectroscopically confirmed by our follow-up spectroscopy and other studies. We obtain UV LFs at $z \sim 4-7$ that span a very wide UV luminosity range of $\sim 0.002 - 100 \, L_{\rm UV}^\ast$ ($-26 < M_{\rm UV} < -14$ mag) by combining LFs from our program and the ultra-deep Hubble Space Telescope legacy surveys. We derive three parameters of the best-fit Schechter function, $\phi^\ast$, $M_{\rm UV}^\ast$, and $\alpha$, of the UV LFs in the magnitude range where the AGN contribution is negligible, and find that $\alpha$ and $\phi^\ast$ decrease from $z\sim 4$ to $7$ with no significant evolution of $M_{\rm UV}^\ast$. Because our HSC SSP data bridge the LFs of galaxies and AGNs with great statistical accuracy, we carefully investigate the bright end of the galaxy UV LFs that are estimated by the subtraction of the AGN contribution either aided with spectroscopy or the best-fit AGN UV LFs. We find that the bright end of the galaxy UV LFs cannot be explained by the Schechter function fits at $> 2 \sigma$ significance, and require either double power-law functions or modified Schechter functions that consider a magnification bias due to gravitational lensing.

We present clustering properties from 579,492 Lyman break galaxies (LBGs) at z~4-6 over the 100 deg^2 sky (corresponding to a 1.4 Gpc^3 volume) identified in early data of the Hyper Suprime-Cam (HSC) Subaru strategic program survey. We derive angular correlation functions (ACFs) of the HSC LBGs with unprecedentedly high statistical accuracies at z~4-6, and compare them with the halo occupation distribution (HOD) models. We clearly identify significant ACF excesses in 10"<$\theta$<90", the transition scale between 1- and 2-halo terms, suggestive of the existence of the non-linear halo bias effect. Combining the HOD models and previous clustering measurements of faint LBGs at z~4-7, we investigate dark-matter halo mass (Mh) of the z~4-7 LBGs and its correlation with various physical properties including the star-formation rate (SFR), the stellar-to-halo mass ratio (SHMR), and the dark matter accretion rate (dotMh) over a wide-mass range of Mh/M$_\odot$=4x10^10-4x10^12. We find that the SHMR increases from z~4 to 7 by a factor of ~4 at Mh~1x10^11 M$_\odot$, while the SHMR shows no strong evolution in the similar redshift range at Mh~1x10^12 M$_\odot$. Interestingly, we identify a tight relation of SFR/dotMh-Mh showing no significant evolution beyond 0.15 dex in this wide-mass range over z~4-7. This weak evolution suggests that the SFR/dotMh-Mh relation is a fundamental relation in high-redshift galaxy formation whose star formation activities are regulated by the dark matter mass assembly. Assuming this fundamental relation, we calculate the cosmic SFR densities (SFRDs) over z=0-10 (a.k.a. Madau-Lilly plot). The cosmic SFRD evolution based on the fundamental relation agrees with the one obtained by observations, suggesting that the cosmic SFRD increase from z~10 to 4-2 (decrease from z~4-2 to 0) is mainly driven by the increase of the halo abundance (the decrease of the accretion rate).

We have measured the clustering of galaxies around active galactic nuclei (AGN) for which single-epoch virial masses of the super-massive black hole (SMBH) are available to investigate the relation between the large scale environment of AGNs and the evolution of SMBHs. The AGN samples used in this work were derived from the Sloan Digital Sky Survey (SDSS) observations and the galaxy samples were from 240~deg$^{2}$ S15b data of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). The investigated redshift range is 0.6--3.0, and the masses of the SMBHs lie in the range $10^{7.5}$--$10^{10}$~M$_{\odot}$. The absolute magnitude of the galaxy samples reaches to $M_{\lambda 310}$ $\sim$ $-$18 at rest frame wavelength 310~nm for the low-redshift end of the samples. More than 70\% of the galaxies in the analysis are blue. We found a significant dependence of the cross-correlation length on redshift, which primarily reflects the brightness dependence of the galaxy clustering. At the lowest redshifts the cross-correlation length increases from 7~$h^{-1}$~Mpc around $M_{\lambda 310}$$ = $$-19$~mag to $>$10~$h^{-1}$~Mpc beyond $M_{\lambda 310}$$ = $$-20$~mag. No significant dependence of the cross-correlation length on BH mass was found for whole galaxy samples dominated by blue galaxies, while there was an indication of BH mass dependence in the cross-correlation with red galaxies. These results provides a picture of the environment of AGNs studied in this paper being enriched with blue starforming galaxies, and a fraction of the galaxies are being evolved to red galaxies along with the evolution of SMBHs in that system.

We investigate the galaxy overdensity around proto-cluster scale quasar pairs at high (z>3) and low (z~1) redshift based on the unprecedentedly wide and deep optical survey of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). Using the first-year survey data covering effectively ~121 deg^2 with the 5sigma depth of i~26.4 and the SDSS DR12Q catalog, we find two luminous pairs at z~3.3 and 3.6 which reside in >5sigma overdense regions of g-dropout galaxies at i<25. The projected separations of the two pairs are R_perp=1.75 and 1.04 proper Mpc, and their velocity offsets are Delta V=692 and 1448 km s^{-1}, respectively. This result is in clear contrast to the average z~4 quasar environments as discussed in Uchiyama et al. (2017) and implies that the quasar activities of the pair members are triggered via major mergers in proto-clusters, unlike the vast majority of isolated quasars in general fields that may turn on via non-merger events such as bar and disk instabilities. At z~1, we find 37 pairs with R_perp<2 pMpc and Delta V<2300 km s^{-1} in the current HSC-Wide coverage, including four from Hennawi et al. (2006). The distribution of the peak overdensity significance within two arcminutes around the pairs has a long tail toward high density (>4sigma) regions. Thanks to the large sample size, we find a statistical evidence that this excess is unique to the pair environments when compared to single quasar and randomly selected galaxy environments at the same redshift range. Moreover, there are nine small-scale (R_perp<1 pMpc) pairs, two of which are found to reside in cluster fields. Our results demonstrate that <2 pMpc-scale quasar pairs at both redshift range tend to occur in massive haloes, although perhaps not the most massive ones, and that they are useful to search for rare density peaks.

We present the result of our spectroscopic follow-up observation for faint quasar candidates at z~5 in a part of the Canada-France-Hawaii Telescope Legacy Survey wide field. We select nine photometric candidates and identify three z~5 faint quasars, one z~4 faint quasar, and a late-type star. Since two faint quasar spectra show Civ emission line without suffering from a heavy atmospheric absorption, we estimate the black hole mass (M$_{BH}$) and Eddington ratio (L/L$_{Edd}$) of them. The inferred log M$_{BH}$ are 9.04+/-0.14 and 8.53+/-0.20, respectively. In addition, the inferred log (L/L$_{Edd}$) are -1.00+/-0.15 and -0.42+/-0.22, respectively. If we adopt that L/L$_{Edd}$= constant or $\propto$ (1+z)^2, the seed black hole masses (M$_{seed}$) of our z~5 faint quasars are expected to be >10^5 M$_\odot$ in most cases. We also compare the observational results with a mass accretion model where angular momentum is lost due to supernova explosions (Kawakatu & Wada 2008). Accordingly, M$_{BH}$ of the z~5 faint quasars in our sample can be explained even if M$_{seed}$ is ~10^3M$_\odot$. Since z~6 luminous qusars and our z~5 faint quasars are not on the same evolutionary track, z~6 luminous quasars and our z~5 quasars are not the same populations but different populations, due to the difference of a period of the mass supply from host galaxies. Furthermore, we confirm that one can explain M$_{BH}$ of z~6 luminous quasars and our z~5 faint quasars even if their seed black holes of them are formed at z~7.

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

We examine the clustering of quasars over a wide luminosity range, by utilizing 901 quasars at $\overline{z}_{\rm phot}\sim3.8$ with $-24.73<M_{\rm 1450}<-22.23$ photometrically selected from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) S16A Wide2 date release and 342 more luminous quasars at $3.4<z_{\rm spec}<4.6$ having $-28.0<M_{\rm 1450}<-23.95$ from the Sloan Digital Sky Survey (SDSS) that fall in the HSC survey fields. We measure the bias factors of two quasar samples by evaluating the cross-correlation functions (CCFs) between the quasar samples and 25790 bright $z\sim4$ Lyman Break Galaxies (LBGs) in $M_{\rm 1450}<-21.25$ photometrically selected from the HSC dataset. Over an angular scale of \timeform{10.0"} to \timeform{1000.0"}, the bias factors are $5.93^{+1.34}_{-1.43}$ and $2.73^{+2.44}_{-2.55}$ for the low and high luminosity quasars, respectively, indicating no luminosity dependence of quasar clustering at $z\sim4$. It is noted that the bias factor of the luminous quasars estimated by the CCF is smaller than that estimated by the auto-correlation function (ACF) over a similar redshift range, especially on scales below \timeform{40.0"}. Moreover, the bias factor of the less-luminous quasars implies the minimal mass of their host dark matter halos (DMHs) is $0.3$-$2\times10^{12}h^{-1}M_{\odot}$, corresponding to a quasar duty cycle of $0.001$-$0.06$.

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

We investigate the relation between star formation rate (SFR) and stellar mass ($M_*$) in sub-galactic ($\sim 1$kpc) scale of 93 local ($0.01<z<0.02$) massive ($M_*>10^{10.5}M_{\odot}$) spiral galaxies. To derive spatially-resolved SFR and stellar mass, we perform so-called pixel-to-pixel SED fitting, which fits an observed spatially-resolved multiband SED with a library of model SEDs using Bayesian statistics approach. We use 2 bands (FUV and NUV) and 5 bands ($u$, $g$, $r$, $i$, and $z$) imaging data from Galaxy Evolution Explorer (GALEX) and Sloan Digital Sky Survey (SDSS), respectively. We find a tight nearly linear relation between the local surface density of SFR ($\Sigma_{\rm{SFR}}$) and stellar mass ($\Sigma_{*}$) which has flattening in high $\Sigma_{*}$. The near linear relation between $\Sigma_{*}$ and $\Sigma_{\rm SFR}$ suggests constant sSFR throughout the galaxies, and the scatter of the relation is directly related to that of sSFR. Therefore, we analyse the variation of sSFR in various scales. More massive galaxies on average have lower sSFR throughout them than less massive galaxies. We also find that barred galaxies have lower sSFR in a core region than non-barred galaxies. However, in the outside region, sSFR of barred and non-barred galaxies are similar and lead to the similar total sSFR.

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.

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.

The chemical abundances for five metal-poor stars in and towards the Galactic bulge have been determined from H-band infrared spectroscopy taken with the RAVEN multi-object adaptive optics science demonstrator and the IRCS spectrograph at the Subaru 8.2-m telescope. Three of these stars are in the Galactic bulge and have metallicities between -2.1 < [Fe/H] < -1.5, and high [alpha/Fe] ~+0.3, typical of Galactic disk and bulge stars in this metallicity range; [Al/Fe] and [N/Fe] are also high, whereas [C/Fe] < +0.3. An examination of their orbits suggests that two of these stars may be confined to the Galactic bulge and one is a halo trespasser, though proper motion values used to calculate orbits are quite uncertain. An additional two stars in the globular cluster M22 show [Fe/H] values consistent to within 1 sigma, although one of these two stars has [Fe/H] = -2.01 +/- 0.09, which is on the low end for this cluster. The [alpha/Fe] and [Ni/Fe] values differ by 2 sigma, with the most metal-poor star showing significantly higher values for these elements. M22 is known to show element abundance variations, consistent with a multi-population scenario (i.e. Marino et al. 2009, 2011; Alves-Brito et al. 2012) though our results cannot discriminate this clearly given our abundance uncertainties. This is the first science demonstration of multi-object adaptive optics with high resolution infrared spectroscopy, and we also discuss the feasibility of this technique for use in the upcoming era of 30-m class telescope facilities.

We discuss images of the star clusters GLIMPSE C01 (GC01) and GLIMPSE C02 (GC02) that were recorded with the Subaru IRCS. Distortions in the wavefront were corrected with the RAVEN adaptive optics (AO) science demonstrator, allowing individual stars in the central regions of both clusters -- where the fractional contamination from non-cluster objects is lowest -- to be imaged. In addition to J, H, and K' images, both clusters were observed through a narrow-band filter centered near 3.05um; GC01 was also observed through two other narrow-band filters that sample longer wavelengths. Stars in the narrow-band images have a FWHM that is close to the telescope diffraction limit, demonstrating that open loop AO systems like RAVEN can deliver exceptional image quality. The near-infrared color magnitude diagram of GC01 is smeared by non-uniform extinction with a dispersion +/- 0.13 magnitudes in A_K. The Red Clump is identified in the K luminosity function (LF) of GC01, and a distance modulus of 13.6 is found. The K LF of GC01 is consistent with a system that is dominated by stars with an age > 1 Gyr. As for GC02, the K LF is flat for K > 16, and the absence of a sub-giant branch argues against an old age if the cluster is at a distance of ~ 7 kpc. Archival SPITZER [3.6] and [4.5] images of the clusters are also examined, and the red giant branch-tip is identified.

Prior statistical knowledge of the atmospheric turbulence is essential for designing, optimizing and evaluating tomographic adaptive optics systems. We present the statistics of the vertical profiles of $C_N^2$ and the outer scale at Maunakea estimated using a Slope Detection And Ranging (SLODAR) method from on-sky telemetry taken by RAVEN, which is a MOAO demonstrator in the Subaru telescope. In our SLODAR method, the profiles are estimated by a fit of the theoretical auto- and cross-correlation of measurements from multiple Shack-Haltmann wavefront sensors to the observed correlations via the non-linear Levenberg-Marquardt Algorithm (LMA), and the analytic derivatives of the spatial phase structure function with respect to its parameters for the LMA are also developed. The estimated profile has the median total seeing of 0.460$^{\prime\prime}$ and large $C_N^2$ fraction of the ground layer of 54.3%. The $C_N^2$ profile has a good agreement with the result from literatures, except for the ground layer. The median value of the outer scale is 25.5m and the outer scale is larger at higher altitudes, and these trends of the outer scale are consistent with findings in literatures.

We constrain the rate of gas inflow into and outflow from a main-sequence star-forming galaxy at z~1.4 by fitting a simple analytic model for the chemical evolution in a galaxy to the observational data of the stellar mass, metallicity, and molecular gas mass fraction. The molecular gas mass is derived from CO observations with a metallicity-dependent CO-to-H2 conversion factor, and the gas metallicity is derived from the H{\alpha} and [NII]{\lambda} 6584 emission line ratio. Using a stacking analysis of CO integrated intensity maps and the emission lines of H{\alpha} and [NII], the relation between stellar mass, metallicity, and gas mass fraction is derived. We constrain the inflow and outflow rates with least-chi-square fitting of a simple analytic chemical evolution model to the observational data. The best-fit inflow and outflow rates are ~1.7 and ~0.4 in units of star-formation rate, respectively. The inflow rate is roughly comparable to the sum of the star-formation rate and outflow rate, which supports the equilibrium model for galaxy evolution; i.e., all inflow gas is consumed by star formation and outflow.

We report the discovery of a new ultra-faint dwarf satellite companion of the Milky Way based on the early survey data from the Hyper Suprime-Cam Subaru Strategic Program. This new satellite, Virgo I, which is located in the constellation of Virgo, has been identified as a statistically significant (5.5 sigma) spatial overdensity of star-like objects with a well-defined main sequence and red giant branch in their color-magnitude diagram. The significance of this overdensity increases to 10.8 sigma when the relevant isochrone filter is adopted for the search. Based on the distribution of the stars around the likely main sequence turn-off at r ~ 24 mag, the distance to Virgo I is estimated as 87 kpc, and its most likely absolute magnitude calculated from a Monte Carlo analysis is M_V = -0.8 +/- 0.9 mag. This stellar system has an extended spatial distribution with a half-light radius of 38 +12/-11 pc, which clearly distinguishes it from a globular cluster with comparable luminosity. Thus, Virgo I is one of the faintest dwarf satellites known and is located beyond the reach of the Sloan Digital Sky Survey. This demonstrates the power of this survey program to identify very faint dwarf satellites. This discovery of VirgoI is based only on about 100 square degrees of data, thus a large number of faint dwarf satellites are likely to exist in the outer halo of the Milky Way.

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.

In tomographic adaptive-optics (AO) systems, errors due to tomographic wave-front reconstruction limit the performance and angular size of the scientific field of view (FoV), where AO correction is effective. We propose a multi time-step tomographic wave-front reconstruction method to reduce the tomographic error by using the measurements from both the current and the previous time-steps simultaneously. We further outline the method to feed the reconstructor with both wind speed and direction of each turbulence layer. An end-to-end numerical simulation, assuming a multi-object AO (MOAO) system on a 30 m aperture telescope, shows that the multi time-step reconstruction increases the Strehl ratio (SR) over a scientific FoV of 10 arcminutes in diameter by a factor of 1.5--1.8 when compared to the classical tomographic reconstructor, depending on the guide star asterism and with perfect knowledge of wind speeds and directions. We also evaluate the multi time-step reconstruction method and the wind estimation method on the RAVEN demonstrator under laboratory setting conditions. The wind speeds and directions at multiple atmospheric layers are measured successfully in the laboratory experiment by our wind estimation method with errors below 2 \ms. With these wind estimates, the multi time-step reconstructor increases the SR value by a factor of 1.2--1.5, which is consistent with a prediction from end-to-end numerical simulation.

We present basic properties of $\sim$3,300 emission line galaxies detected by the FastSound survey, which are mostly H$\alpha$ emitters at $z \sim$ 1.2-1.5 in the total area of about 20 deg$^2$, with the H$\alpha$ flux sensitivity limit of $\sim 1.6 \times 10^{-16} \rm erg \ cm^{-2} s^{-1}$ at 4.5 sigma. This paper presents the catalogs of the FastSound emission lines and galaxies, which will be open to the public in the near future. We also present basic properties of typical FastSound H$\alpha$ emitters, which have H$\alpha$ luminosities of $10^{41.8}$-$10^{43.3}$ erg/s, SFRs of 20--500 $M_\odot$/yr, and stellar masses of $10^{10.0}$--$10^{11.3}$ $M_\odot$. The 3D distribution maps for the four fields of CFHTLS W1--4 are presented, clearly showing large scale clustering of galaxies at the scale of $\sim$ 100--600 comoving Mpc. Based on 1,105 galaxies with detections of multiple emission lines, we estimate that contamination of non-H$\alpha$ lines is about 4% in the single-line emission galaxies, which are mostly [OIII]$\lambda$5007. This contamination fraction is also confirmed by the stacked spectrum of all the FastSound spectra, in which H$\alpha$, [NII]$\lambda \lambda$6548,6583, [SII]$\lambda \lambda$6717, 6731, and [OI]$\lambda \lambda$6300,6364 are seen.

We measure the redshift-space correlation function from a spectroscopic sample of 2783 emission line galaxies from the FastSound survey. The survey, which uses the Subaru Telescope and covers the redshift ranges of $1.19<z<1.55$, is the first cosmological study at such high redshifts. We detect clear anisotropy due to redshift-space distortions (RSD) both in the correlation function as a function of separations parallel and perpendicular to the line of sight and its quadrupole moment. RSD has been extensively used to test general relativity on cosmological scales at $z<1$. Adopting a LCDM cosmology with the fixed expansion history and no velocity dispersion $\sigma_{\rm v}=0$, and using the RSD measurements on scales above 8Mpc/h, we obtain the first constraint on the growth rate at the redshift, $f(z)\sigma_8(z)=0.482\pm 0.116$ at $z\sim 1.4$ after marginalizing over the galaxy bias parameter $b(z)\sigma_8(z)$. This corresponds to $4.2\sigma$ detection of RSD. Our constraint is consistent with the prediction of general relativity $f\sigma_8\sim 0.392$ within the $1-\sigma$ confidence level. When we allow $\sigma_{\rm v}$ to vary and marginalize it over, the growth rate constraint becomes $f\sigma_8=0.494^{+0.126}_{-0.120}$. We also demonstrate that by combining with the low-z constraints on $f\sigma_8$, high-z galaxy surveys like the FastSound can be useful to distinguish modified gravity models without relying on CMB anisotropy experiments.

We conducted observations of 12CO(J=5-4) and dust thermal continuum emission toward twenty star-forming galaxies on the main sequence at z~1.4 using ALMA to investigate the properties of the interstellar medium. The sample galaxies are chosen to trace the distributions of star-forming galaxies in diagrams of stellar mass-star formation rate and stellar mass-metallicity. We detected CO emission lines from eleven galaxies. The molecular gas mass is derived by adopting a metallicity-dependent CO-to-H2 conversion factor and assuming a CO(5-4)/CO(1-0) luminosity ratio of 0.23. Molecular gas masses and its fractions (molecular gas mass/(molecular gas mass + stellar mass)) for the detected galaxies are in the ranges of (3.9-12) x 10^{10} Msun and 0.25-0.94, respectively; these values are significantly larger than those in local spiral galaxies. The molecular gas mass fraction decreases with increasing stellar mass; the relation holds for four times lower stellar mass than that covered in previous studies, and that the molecular gas mass fraction decreases with increasing metallicity. Stacking analyses also show the same trends. The dust thermal emissions were clearly detected from two galaxies and marginally detected from five galaxies. Dust masses of the detected galaxies are (3.9-38) x 10^{7} Msun. We derived gas-to-dust ratios and found they are 3-4 times larger than those in local galaxies. The depletion times of molecular gas for the detected galaxies are (1.4-36) x 10^{8} yr while the results of the stacking analysis show ~3 x 10^{8} yr. The depletion time tends to decrease with increasing stellar mass and metallicity though the trend is not so significant, which contrasts with the trends in local galaxies.

Aims. We present a spectroscopic study of the properties of 64 Balmer break galaxies that show signs of star formation. The studied sample of star-forming galaxies spans a redshift range from 0.094 to 1.475 with stellar masses in the range 10$^{8}-$10$^{12}$ $M_{\odot}$. The sample also includes eight broad emission line galaxies with redshifts between 1.5 $<z<$ 3.0. Methods. We derived star formation rates (SFRs) from emission line luminosities and investigated the dependence of the SFR and specific SFR (SSFR) on the stellar mass and color. Furthermore, we investigated the evolution of these relations with the redshift. Results. We found that the SFR correlates with the stellar mass, our data is consistent with previous results from other authors in that there is a break in the correlation, which reveals the presence of massive galaxies with lower SFR values (i.e., decreasing star formation). We also note an anticorrelation for the SSFR with the stellar mass. Again in this case, our data is also consistent with a break in the correlation, revealing the presence of massive star-forming galaxies with lower SSFR values, thereby increasing the anticorrelation. These results might suggest a characteristic mass ($M_{0}$) at which the red sequence could mostly be assembled. In addition, at a given stellar mass, high-redshift galaxies have on average higher SFR and SSFR values than local galaxies. Finally, we explored whether a similar trend could be observed with redshift in the SSFR$-(u-B)$ color diagram, and we hypothesize that a possible $(u-B)_{0}$ break color may define a characteristic color for the formation of the red sequence.

We have conducted 1.1 mm ALMA observations of a contiguous $105'' \times 50''$ or 1.5 arcmin$^2$ window in the SXDF-UDS-CANDELS. We achieved a 5$\sigma$ sensitivity of 0.28 mJy, providing a flat sensus of dusty star-forming galaxies with $L_{\rm IR} \sim6\times10^{11}$ $L_\odot$ (for $T_{\rm dust}$ =40K) up to $z\sim10$ thanks to the negative K-correction at this wavelength. We detected 5 brightest sources (S/N$>$6) and 18 low-significance sources (5$>$S/N$>$4; these may contain spurious detections, though). One of the 5 brightest ALMA sources ($S_{\rm 1.1mm} = 0.84 \pm 0.09$ mJy) is extremely faint in the WFC3 and VLT/HAWK-I images, demonstrating that a contiguous ALMA imaging survey is able to uncover a faint dust-obscured population that is invisible in deep optical/near-infrared surveys. We found a possible [CII]-line emitter at $z=5.955$ or a low-$z$ CO emitting galaxy within the field, which may allow us to constrain the [CII] and/or the CO luminosity functions across the history of the universe.

Near-infrared (NIR) spectra that have an angular resolution of ~ 0.15 arcsec are used to examine the stellar content of the central regions of the nearby elliptical galaxy Maffei 1. The spectra were recorded at the Subaru Telescope, with wavefront distortions corrected by the RAVEN Multi-Object Adaptive Optics science demonstrator. The Ballick-Ramsey C_2 absorption bandhead near 1.76 microns is detected, and models in which 10 - 20% of the light near 1.8 microns originates from stars of spectral type C5 reproduce this feature. Archival NIR and mid-infrared images are also used to probe the structural and photometric properties of the galaxy. Comparisons with models suggest that an intermediate age population dominates the spectral energy distribution between 1 and 5 microns near the galaxy center. This is consistent not only with the presence of C stars, but also with the large HBeta index that has been measured previously for Maffei 1. The J-K color is more-or-less constant within 15 arcsec of the galaxy center, suggesting that the brightest red stars are well-mixed in this area.

In order to investigate the origin of quasars, we estimate the bias factor for low-luminosity quasars at high redshift for the first time. In this study, we use the two-point angular cross-correlation function (CCF) for both low-luminosity quasars at $-24<M_{\rm 1450}<-22$ and Lyman-break galaxies (LBGs). Our sample consists of both 25 low-luminosity quasars (16 objects are spectroscopically confirmed low-luminosity quasars) in the redshift range $3.1<z<4.5$ and 835 color-selected LBGs with $z^{\prime}_{\rm LBG}<25.0$ at $z\sim4$ in the COSMOS field. We have made our analysis for the following two quasar samples; (1) the spectroscopic sample (the 16 quasars confirmed by spectroscopy), and (2) the total sample (the 25 quasars including 9 quasars with photometric redshifts). The bias factor for low-luminosity quasars at $z\sim4$ is derived by utilizing the quasar-LBG CCF and the LBG auto-correlation function. We then obtain the $86\%$ upper limits of the bias factors for low-luminosity quasars, that are 5.63 and 10.50 for the total and the spectroscopic samples, respectively. These bias factors correspond to the typical dark matter halo masses, log $(M_{\rm DM}/(h^{-1}M_{\odot}))=$$12.7$ and $13.5$, respectively. This result is not inconsistent with the predicted bias for quasars which is estimated by the major merger models.

We present the results from a large near-infrared spectroscopic survey with Subaru/FMOS (\textit{FastSound}) consisting of $\sim$ 4,000 galaxies at $z\sim1.4$ with significant H$\alpha$ detection. We measure the gas-phase metallicity from the [N~{\sc ii}]$\lambda$6583/H$\alpha$ emission line ratio of the composite spectra in various stellar mass and star-formation rate bins. The resulting mass-metallicity relation generally agrees with previous studies obtained in a similar redshift range to that of our sample. No clear dependence of the mass-metallicity relation with star-formation rate is found. Our result at $z\sim1.4$ is roughly in agreement with the fundamental metallicity relation at $z\sim0.1$ with fiber aperture corrected star-formation rate. We detect significant [S~{\sc ii}]$\lambda\lambda$6716,6731 emission lines from the composite spectra. The electron density estimated from the [S~{\sc ii}]$\lambda\lambda$6716,6731 line ratio ranges from 10 -- 500 cm$^{-3}$, which generally agrees with that of local galaxies. On the other hand, the distribution of our sample on [N~{\sc ii}]$\lambda$6583/H$\alpha$ vs. [S~{\sc ii}]$\lambda\lambda$6716,6731/H$\alpha$ is different from that found locally. We estimate the nitrogen-to-oxygen abundance ratio (N/O) from the N2S2 index, and find that the N/O in galaxies at $z\sim1.4$ is significantly higher than the local values at a fixed metallicity and stellar mass. The metallicity at $z\sim1.4$ recalculated with this N/O enhancement taken into account decreases by 0.1 -- 0.2 dex. The resulting metallicity is lower than the local fundamental metallicity relation.

We report optical-infrared (IR) properties of faint 1.3 mm sources (S_1.3mm = 0.2-1.0 mJy) detected with the Atacama Large Millimeter/submillimeter Array (ALMA) in the Subaru/XMM-Newton Deep Survey (SXDS) field. We searched for optical/IR counterparts of 8 ALMA-detected sources (>=4.0 sigma, the sum of the probability of spurious source contamination is ~1) in a K-band source catalog. Four ALMA sources have K-band counterpart candidates within a 0.4" radius. Comparison between ALMA-detected and undetected K-band sources in the same observing fields shows that ALMA-detected sources tend to be brighter, more massive, and more actively forming stars. While many of the ALMA-identified submillimeter-bright galaxies (SMGs) in previous studies lie above the sequence of star-forming galaxies in stellar mass--star-formation rate plane, our ALMA sources are located in the sequence, suggesting that the ALMA-detected faint sources are more like `normal' star-forming galaxies rather than `classical' SMGs. We found a region where multiple ALMA sources and K-band sources reside in a narrow photometric redshift range (z ~ 1.3-1.6) within a radius of 5" (42 kpc if we assume z = 1.45). This is possibly a pre-merging system and we may be witnessing the early phase of formation of a massive elliptical galaxy.

The TMT Detailed Science Case describes the transformational science that the Thirty Meter Telescope will enable. Planned to begin science operations in 2024, TMT will open up opportunities for revolutionary discoveries in essentially every field of astronomy, astrophysics and cosmology, seeing much fainter objects much more clearly than existing telescopes. Per this capability, TMT's science agenda fills all of space and time, from nearby comets and asteroids, to exoplanets, to the most distant galaxies, and all the way back to the very first sources of light in the Universe. More than 150 astronomers from within the TMT partnership and beyond offered input in compiling the new 2015 Detailed Science Case. The contributing astronomers represent the entire TMT partnership, including the California Institute of Technology (Caltech), the Indian Institute of Astrophysics (IIA), the National Astronomical Observatories of the Chinese Academy of Sciences (NAOC), the National Astronomical Observatory of Japan (NAOJ), the University of California, the Association of Canadian Universities for Research in Astronomy (ACURA) and US associate partner, the Association of Universities for Research in Astronomy (AURA).

We report the multi-wavelength identification of the X-ray sources found in the Subaru-XMM-Newton Deep Survey (SXDS) using deep imaging data covering the wavelength range between the far-UV to the mid-IR. We select a primary counterpart of each X-ray source by applying the likelihood ratio method to R-band, 3.6micron, near-UV, and 24micron source catalogs as well as matching catalogs of AGN candidates selected in 1.4GHz radio and i'-band variability surveys. Once candidates of Galactic stars, ultra-luminous X-ray sources in a nearby galaxy, and clusters of galaxies are removed there are 896 AGN candidates in the sample. We conduct spectroscopic observations of the primary counterparts with multi-object spectrographs in the optical and NIR; 65\% of the X-ray AGN candidates are spectroscopically-identified. For the remaining X-ray AGN candidates, we evaluate their photometric redshift with photometric data in 15 bands. Utilising the multi-wavelength photometric data of the large sample of X-ray selected AGNs, we evaluate the stellar masses, M*, of the host galaxies of the narrow-line AGNs. The distribution of the stellar mass is remarkably constant from z=0.1 to 4.0. The relation between M* and 2--10 keV luminosity can be explained with strong cosmological evolution of the relationship between the black hole mass and M*. We also evaluate the scatter of the UV-MIR spectral energy distribution (SED) of the X-ray AGNs as a function of X-ray luminosity and absorption to the nucleus. The scatter is compared with galaxies which have redshift and stellar mass distribution matched with the X-ray AGN. The UV-NIR SEDs of obscured X-ray AGNs are similar to those of the galaxies in the matched sample. In the NIR-MIR range, the median SEDs of X-ray AGNs are redder, but the scatter of the SEDs of the X-ray AGN broadly overlaps that of the galaxies in the matched sample.

FastSound is a galaxy redshift survey using the near-infrared Fiber Multi-Object Spectrograph (FMOS) mounted on the Subaru Telescope, targeting H$\alpha$ emitters at $z \sim 1.18$--$1.54$ down to the sensitivity limit of H$\alpha$ flux $\sim 2 \times 10^{-16} \ \rm erg \ cm^{-2} s^{-1}$. The primary goal of the survey is to detect redshift space distortions (RSD), to test General Relativity by measuring the growth rate of large scale structure and to constrain modified gravity models for the origin of the accelerated expansion of the universe. The target galaxies were selected based on photometric redshifts and H$\alpha$ flux estimates calculated by fitting spectral energy distribution (SED) models to the five optical magnitudes of the Canada France Hawaii Telescope Legacy Survey (CFHTLS) Wide catalog. The survey started in March 2012, and all the observations were completed in July 2014. In total, we achieved $121$ pointings of FMOS (each pointing has a $30$ arcmin diameter circular footprint) covering $20.6$ deg$^2$ by tiling the four fields of the CFHTLS Wide in a hexagonal pattern. Emission lines were detected from $\sim 4,000$ star forming galaxies by an automatic line detection algorithm applied to 2D spectral images. This is the first in a series of papers based on FastSound data, and we describe the details of the survey design, target selection, observations, data reduction, and emission line detections.

Using the HST/WFC3 and ACS multi-band imaging data taken in CANDELS and 3D-HST, we study the general properties and the diversity of the progenitors of the Milky Way (MWs) and local massive galaxy (MGs) at 0.5 < z < 3.0, based on a constant cumulative number density analysis. After careful data reduction and stacking analysis, we conduct a radially resolved pixel SED fitting to obtain the radial distributions of the stellar mass and rest-frame colors. The stellar mass of MWs increases in self-similar way, irrespective of the radial distance, while that of MGs grows in inside-out way where they obtain ~ 75% of the total mass at outer (> 2.5 kpc) radius since z ~ 2. Although the radial mass profiles evolve in distinct ways, the formation and quenching of the central dense region (or bulge) ahead of the outer disk formation are found to be common for both systems. The sudden reddening of bulge at z ~ 1.6 and z ~ 2.4 for MWs and MGs, respectively, suggests the formation of bulge and would give a clue to the different gas accretion histories and quenching. A new approach to evaluate the morphological diversity is conducted by using the average surface density profile and its dispersion. The variety of the radial mass profiles for MGs peaks at higher redshift (z > 2.8), and then rapidly converges to more uniform shape at z < 1.5, while that for MWs remains in the outer region over the redshift. Compared with the observed star formation rates and color profiles, the evolution of variety is consistently explained by the star formation activities.

We construct the rest-frame 2--10 keV intrinsic X-ray luminosity function of Active Galactic Nuclei (AGNs) from a combination of X-ray surveys from the all-sky Swift BAT survey to the Chandra Deep Field-South. We use ~3200 AGNs in our analysis, which covers six orders of magnitude in flux. The inclusion of the XMM and Chandra COSMOS data has allowed us to investigate the detailed behavior of the XLF and evolution. In deriving our XLF, we take into account realistic AGN spectrum templates, absorption corrections, and probability density distributions in photometric redshift. We present an analytical expression for the overall behavior of the XLF in terms of the luminosity-dependent density evolution, smoothed two power-law expressions in 11 redshift shells, three-segment power-law expression of the number density evolution in four luminosity classes, and binned XLF. We observe a sudden flattening of the low luminosity end slope of the XLF slope at z>~0.6. Detailed structures of the AGN downsizing have been also revealed, where the number density curves have two clear breaks at all luminosity classes above log LX>43. The two break structure is suggestive of two-phase AGN evolution, consisting of major merger triggering and secular processes.

We try to constrain the gas inflow and outflow rate of star-forming galaxies at $z\sim1.4$ by employing a simple analytic model for the chemical evolution of galaxies. The sample is constructed based on a large near-infrared (NIR) spectroscopic sample observed with Subaru/FMOS. The gas-phase metallicity is measured from the [\ion{N}{2}]$\lambda$6584/H$\alpha$ emission line ratio and the gas mass is derived from the extinction corrected H$\alpha$ luminosity by assuming the Kennicutt-Schmidt law. We constrain the inflow and outflow rate from the least-$\chi^{2}$ fittings of the observed gas mass fraction, stellar mass, and metallicity with the analytic model. The joint $\chi^{2}$ fitting shows the best-fit inflow rate is $\sim1.8$ and the outflow rate is $\sim0.6$ in unit of star-formation rate (SFR). By applying the same analysis to the previous studies at $z\sim0$ and $z\sim2.2$, it is shown that the both inflow rate and outflow rate decrease with decreasing redshift, which implies the higher activity of gas flow process at higher redshift. The decreasing trend of the inflow rate from $z\sim2.2$ to $z\sim0$ agrees with that seen in the previous observational works with different methods, though the absolute value is generally larger than the previous works. The outflow rate and its evolution from $z\sim2.2$ to $z\sim0$ obtained in this work agree well with the independent estimations in the previous observational works.

We describe the development of automated emission line detection software for the Fiber Multi-Object Spectrograph (FMOS), which is a near-infrared spectrograph fed by $400$ fibers from the $0.2$ deg$^2$ prime focus field of view of the Subaru Telescope. The software, FIELD (FMOS software for Image-based Emission Line Detection), is developed and tested mainly for the FastSound survey, which is targeting H$\alpha$ emitting galaxies at $z \sim 1.3$ to measure the redshift space distortion as a test of general relativity beyond $z \sim 1$. The basic algorithm is to calculate the line signal-to-noise ratio ($S/N$) along the wavelength direction, given by a 2-D convolution of the spectral image and a detection kernel representing a typical emission line profile. A unique feature of FMOS is its use of OH airglow suppression masks, requiring the use of flat-field images to suppress noise around the mask regions. Bad pixels on the detectors and pixels affected by cosmic-rays are efficiently removed by using the information obtained from the FMOS analysis pipeline. We limit the range of acceptable line-shape parameters for the detected candidates to further improve the reliability of line detection. The final performance of line detection is tested using a subset of the FastSound data; the false detection rate of spurious objects is examined by using inverted frames obtained by exchanging object and sky frames. The false detection rate is $< 1$\% at $S/N > 5$, allowing an efficient and objective emission line search for FMOS data at the line flux level of $\gtrsim 1.0 \times 10^{-16}$[erg/cm$^2$/s].

We present the most up-to-date X-ray luminosity function (XLF) and absorption function of Active Galactic Nuclei (AGNs) over the redshift range from 0 to 5, utilizing the largest, highly complete sample ever available obtained from surveys performed with Swift/BAT, MAXI, ASCA, XMM-Newton, Chandra, and ROSAT. The combined sample, including that of the Subaru/XMM-Newton Deep Survey, consists of 4039 detections in the soft (0.5--2 keV) and/or hard ($>2$ keV) band. We utilize a maximum likelihood method to reproduce the count-rate versus redshift distribution for each survey, by taking into account the evolution of the absorbed fraction, the contribution from Compton-thick (CTK) AGNs, and broad band spectra of AGNs including reflection components from tori based on the luminosity and redshift dependent unified scheme. We find that the shape of the XLF at $z \sim 1-3$ is significantly different from that in the local universe, for which the luminosity dependent density evolution model gives much better description than the luminosity and density evolution model. These results establish the standard population synthesis model of the X-Ray Background (XRB), which well reproduces the source counts, the observed fractions of CTK AGNs, and the spectrum of the hard XRB. The number ratio of CTK AGNs to the absorbed Compton-thin (CTN) AGNs is constrained to be $\approx$0.5--1.6 to produce the 20--50 keV XRB intensity within present uncertainties, by assuming that they follow the same evolution as CTN AGNs. The growth history of supermassive black holes is discussed based on the new AGN bolometric luminosity function.

The efficient selection of high-redshift emission galaxies is important for future large galaxy redshift surveys for cosmology. Here we describe the target selection methods for the FastSound project, a redshift survey for H alpha emitting galaxies at z=1.2-1.5 using Subaru/FMOS to measure the linear growth rate f\sigma 8 via Redshift Space Distortion (RSD) and constrain the theory of gravity. To select ~400 target galaxies in the 0.2 deg^2 FMOS field-of-view from photometric data of CFHTLS-Wide (u*g'r'i'z'), we test several different methods based on color-color diagrams or photometric redshift estimates from spectral energy distribution (SED) fitting. We also test the improvement in selection efficiency that can be achieved by adding near-infrared data from the UKIDSS DXS (J). The success rates of H alpha detection with FMOS averaged over two observed fields using these methods are 11.3% (color-color, optical), 13.6% (color-color, optical+NIR), 17.3% (photo-z, optical), and 15.1% (photo-z, optical+NIR). Selection from photometric redshifts tends to give a better efficiency than color-based methods, although there is no significant improvement by adding J band data within the statistical scatter. We also investigate the main limiting factors for the success rate, by using the sample of the HiZELS H alpha emitters that were selected by narrow-band imaging. Although the number density of total H alpha emitters having higher H alpha fluxes than the FMOS sensitivity is comparable with the FMOS fiber density, the limited accuracy of photometric redshift and H alpha flux estimations have comparable effects on the success rate of <~20% obtained from SED fitting.

We present the first systematic search for extended metal-line [OII]{\lambda}{\lambda}3726,3729 nebulae, or [OII] blobs (OIIBs), at z=1.2 using deep narrowband imaging with a survey volume of 1.9x10^5 Mpc^3 on the 0.62 deg^2 sky of Subaru-XMM Deep Survey (SXDS) field. We discover a giant OIIB, dubbed 'OIIB 1', with a spatial extent over ~75 kpc at a spectroscopic redshift of z=1.18, and also identify a total of twelve OIIBs with a size of >30 kpc. Our optical spectrum of OIIB 1 presents [NeV]{\lambda}3426 line at the 6{\sigma} level, indicating that this object harbors an obscured type-2 AGN. The presence of gas outflows in this object is suggested by two marginal detections of FeII{\lambda}2587 absorption and FeII*{\lambda}2613 emission lines both of which are blueshifted at as large as 500-600 km/s, indicating that the heating source of OIIB 1 is AGN or associated shock excitation rather than supernovae produced by starbursts. The number density of OIIB 1-type giant blobs is estimated to be ~5x10^{-6} Mpc^{-3} at z~1.2, which is comparable with that of AGNs driving outflow at a similar redshift, suggesting that giant OIIBs are produced only by AGN activity. On the other hand, the number density of small OIIBs, 6x10^{-5} Mpc^{-3}, compared to that of z~1 galaxies in the blue cloud in the same M_B range, may imply that 3% of star-forming galaxies at z~1 are quenching star formation through outflows involving extended [OII] emission.

We present a stellar mass-metallicity relation at z~1.4 with an unprecedentedly large sample of ~340 star-forming galaxies obtained with FMOS on the Subaru Telescope. We observed K-band selected galaxies at 1.2 < z_{ph} < 1.6 in the SXDS/UDS fields with M_{*} > 10^{9.5} M_{\sun}, and expected F(H\alpha) > 5 \times 10^{-17} erg s^{-1} cm^{-2}. Among the observed ~1200 targets, 343 objects show significant H\alpha emission lines. The gas-phase metallicity is obtained from [NII]\lambda 6584/H\alpha line ratio, after excluding possible active galactic nuclei (AGNs). Due to the faintness of the [NII]\lambda 6584 lines, we apply the stacking analysis and derive the mass-metallicity relation at z~1.4. Our results are compared to past results at different redshifts in the literature. The mass-metallicity relation at z~1.4 is located between those at z~0.8 and z~2.2; it is found that the metallicity increases with decreasing redshift from z~3 to z~0 at fixed stellar mass. Thanks to the large size of the sample, we can study the dependence of the mass-metallicity relation on various galaxy physical properties. The average metallicity from the stacked spectra is close to the local FMR in the higher metallicity part but >0.1 dex higher in metallicity than the FMR in the lower metallicity part. We find that galaxies with larger E(B-V), B-R, and R-H colours tend to show higher metallicity by ~0.05 dex at fixed stellar mass. We also find relatively clearer size dependence that objects with smaller half light radius tend to show higher metallicity by ~0.1 dex at fixed stellar mass, especially in the low mass part.

The properties of K-band selected galaxies (K_AB<24) in the z = 3.09 SSA22 protocluster field are studied. 430 galaxies at 2.6 < z_phot < 3.6 are selected as potential protocluster members in a 112 arcmin^2 area based on their photometric redshifts. We find that \approx 20% of the massive galaxies with stellar masses >10^11 M_sun at z_phot \sim 3.1 have colors consistent with those of quiescent galaxies with ages > 0.5 Gyr. This fraction increases to \approx 50% after correcting for unrelated foreground/background objects. We also find that 30% of the massive galaxies are heavily reddened dusty star-forming galaxies. Few such quiescent galaxies at similar redshifts are seen in typical survey fields. An excess surface density of 24\mu m sources at z_phot \sim 3.1 is also observed, implying the presence of dusty star-formation activity in the protocluster. Cross-correlation with the X-ray data indicates that the fraction of K-band selected protocluster galaxies hosting active galactic nuclei (AGN) is also high compared with the field. The sky distribution of the quiescent galaxies, the 24\mu m sources, and the X-ray AGNs show clustering around a density peak of z=3.1 Ly\alpha emitters (LAEs). A significant fraction of the massive galaxies have already become quiescent, while the dusty star-formation is still active in the SSA22 protocluster. These findings indicate that we are witnessing the formation epoch of massive early-type galaxies at the center of predecessors to present-day rich galaxy clusters.

We compare a large sample of galaxies between 0.5<z<2.6 with theoretical predictions for how the optical diagnostic line ratios in galaxy ensembles change as a function of cosmic time. We show that star forming galaxies at high redshift (z>1.5) are consistent with a model in which the ISM conditions are more extreme at high redshift than seen in the global spectra of local galaxies. We speculate that global spectra of our high redshift galaxies may be dominated by HII regions similar to the extreme clumpy, dense star-forming complexes in the Antennae and M82. The transition to local-type conditions occurs between 0.8<z<1.5. We conclude that classification schemes developed for local samples should not be applied at high redshift (z > 1.5). We use our theoretical models to derive a new redshift-dependent classification line that utilizes the standard optical diagnostic line ratios [OIII]/H-beta and [NII]/H-alpha. Our new line can be used to separate star-forming galaxies from AGN between z=0 to z~3.5. We anticipate that our redshift-dependent optical classification line will be useful for future large surveys with near-infrared multi-object spectrographs. We apply our classification line to a sample of gravitationally lensed galaxies at z~2.5. Although limited by small numbers, we show that our classification line is consistent with the position of AGN that have been independently confirmed via other methods.

We present an investigation into how well the properties of the accretion flow onto a supermassive black hole may be coupled to those of the overlying hot corona. To do so, we specifically measure the characteristic spectral index, Gamma, of a power-law energy distribution, over an energy range of 2 to 10 keV, for X-ray selected, broad-lined radio-quiet AGN up to z~2 in COSMOS and E-CDF-S. We test the previously reported dependence between Gamma and black hole mass, FWHM and Eddington ratio using a sample of AGN covering a broad range in these parameters based on both the Mg ii and H-alpha emission lines with the later afforded by recent near infrared spectroscopic observations using Subaru/FMOS. We calculate the Eddington ratios, lambda_Edd, for sources where a bolometric luminosity (L_Bol) has been presented in the literature, based on SED fitting, or, for sources where these data do not exist, we calculate L_Bol using a bolometric correction to the LX, derived from a relationship between the bolometric correction, and LX/L3000. From a sample of 69 X-ray bright sources (>250 counts), where Gamma can be measured with greatest precision, with an estimate of L_Bol, we find a statistically significant correlation between Gamma and lambda_Edd, which is highly significant with a chance probability of 6.59x10^-8. A statistically significant correlation between Gamma and the FWHM of the optical lines is confirmed, but at lower significance than with lambda_Edd indicating that lambda_Edd is the key parameter driving conditions in the corona. Linear regression analysis reveals that Gamma=(0.32+/-0.05)log10 lambda_Edd+(2.27+/-0.06) and Gamma=(-0.69+/-0.11)log10(FWHM/km/s)+(4.44+/-0.42). Our results on Gamma-lambda_Edd are in very good agreement with previous results. (ABRIDGED)

We present the faint end of number counts at 1.3 mm (238 GHz) obtained with the Atacama Large Millimeter/submillimeter Array (ALMA). Band 6 observations were carried out targeting 20 star-forming galaxies at z ~ 1.4 in the Subaru/XMM-Newton Deep Survey field. In the observations, we serendipitously detect 15 sources (>=3.8 sigma, S(1.3 mm) = 0.15-0.61 mJy) other than the targeted sources. We create number counts by using these `sub-mJy sources', which probe the faintest flux range among surveys at millimeter wavelengths. The number counts are consistent with (flux-scaled) number counts at 850 um and 870 um obtained with gravitational lensing clusters. The ALMA number counts agree well with model predictions, which suggest that these sub-mJy populations are more like `normal' star-forming galaxies than `classical' SMGs with intense star-forming activity. In this flux range, ~80% of the extragalactic background light at 1.3 mm is resolved into individual sources.

We present an analysis of broad emission lines observed in moderate-luminosity active galactic nuclei (AGNs), typical of those found in X-ray surveys of deep fields, with the aim to test the validity of single-epoch virial black hole mass estimates. We have acquired near-infrared (NIR) spectra of AGNs up to z ~ 1.8 in the COSMOS and Extended Chandra Deep Field-South Survey, with the Fiber Multi-Object Spectrograph (FMOS) mounted on the Subaru Telescope. These low-resolution NIR spectra provide a significant detection of the broad Halpha line that has been shown to be a reliable probe of black hole mass at low redshift. Our sample has existing optical spectroscopy which provides a detection of MgII, a broad emission line typically used for black hole mass estimation at z > 1. We carry out a spectral-line fitting procedure using both Halpha and MgII to determine the virial velocity of gas in the broad line region, the monochromatic continuum luminosity at 3000 A, and the total Halpha line luminosity. With a sample of 43 AGNs spanning a range of two decades in luminosity (i.e., L ~ 10^44-46 ergs/s), we find a tight correlation between the continuum and line luminosity with a distribution characterized by <log(L_3000/L_Halpha)> = 1.52 and a dispersion sigma = 0.16. There is also a close one-to-one relationship between the FWHM of Halpha and of MgII up to 10000 km/s with a dispersion of 0.14 in the distribution of the logarithm of their ratios. Both of these then lead to there being very good agreement between Halpha- and MgII-based masses over a wide range in black hole mass (i.e., M_BH ~ 10^7-9 M_sun). We do find a small offset in MgII-based masses, relative to those based on Halpha, of +0.17 dex and a dispersion sigma = 0.32. In general, these results demonstrate that local scaling relations, using MgII or Halpha, are applicable for AGN at moderate luminosities and up to z ~ 2.

In order to investigate the growth of super-massive black holes (SMBHs), we construct the black hole mass function (BHMF) and Eddington ratio distribution function (ERDF) of X-ray-selected broad-line AGNs at z~1.4 in the Subaru XMM-Newton Deep Survey field. In this redshift range, a significant part of the accretion growth of SMBHs is thought to be taking place. Black hole masses of X-ray-selected broad-line AGNs are estimated using the width of the broad MgII line and the 3000A monochromatic luminosity. We supplement the MgII FWHM values with the Ha FWHM obtained from our NIR spectroscopic survey. Using the black hole masses of broad-line AGNs at redshifts between 1.18 and 1.68, the binned broad-line AGN BHMF and ERDF are calculated using the Vmax method. To properly account for selection effects that impact the binned estimates, we derive the corrected broad-line AGN BHMF and ERDF by applying the Maximum Likelihood method, assuming that the ERDF is constant regardless of the black hole mass. We do not correct for the non-negligible uncertainties in virial BH mass estimates. If we compare the corrected broad-line AGN BHMF with that in the local Universe, the corrected BHMF at z~1.4 has a higher number density above 10^8 Msolar but a lower number density below that mass range. The evolution may be indicative of a down-sizing trend of accretion activity among the SMBH population. The evolution of broad-line AGN ERDF from z=1.4 to 0 indicates that the fraction of broad-line AGNs with accretion rate close to the Eddington-limit is higher at higher redshifts.

We study the comoving space density of X-ray-selected luminous active galactic nuclei (AGNs) and the obscured AGN fraction at high redshifts ($3 < z < 5$) in the Subaru/{\it XMM-Newton} Deep Survey (SXDS) field. From an X-ray source catalog with high completeness of optical identification thanks to deep optical images, we select a sample of 30 AGNs at $z > 3$ with intrinsic (de-absorbed and rest-frame 2--10 keV) luminosities of $L_{\rm X} = 10^{44-45}$ erg s$^{-1}$ detected in the 0.5--2 keV band, consisting of 20 and 10 objects with spectroscopic and photometric redshifts, respectively. Utilizing the $1/V_{\rm max}$ method, we confirm that the comoving space density of luminous AGNs decreases with redshift above $z > 3$. When combined with the {\it Chandra}-COSMOS result of Civano et al.\ (2011), the density decline of AGNs with $L_{\rm X} = 10^{44-45}$ erg s$^{-1}$ is well represented by a power law of $(1 + z)^{-6.2 \pm 0.9}$. We also determine the fraction of X-ray obscured AGNs with $N_{\rm H} > 10^{22}$ cm$^{-2}$ in the Compton-thin population to be 0.54$^{+0.17}_{-0.19}$, by carefully taking into account observational biases including the effects of photon statistics for each source. This result is consistent with an independent determination of the type-2 AGN fraction based on optical properties, for which the fraction is found to be 0.59$\pm$0.09. Comparing our result with that obtained in the local Universe, we conclude that the obscured fraction of luminous AGNs increases significantly from $z=0$ to $z>3$ by a factor of 2.5$\pm$1.1.

We investigate the properties (e.g. star formation rate, dust attentuation, stellar mass and metallicity) of a sample of infrared luminous galaxies at z \sim 1 via near-IR spectroscopy with Subaru-FMOS. Our sample consists of Herschel SPIRE and Spitzer MIPS selected sources in the COSMOS field with photometric redshifts in the range 0.7 < z-phot < 1.8, which have been targeted in 2 pointings (0.5 sq. deg.) with FMOS. We find a modest success rate for emission line detections, with candidate H{\alpha} emission lines detected for 57 of 168 SPIRE sources (34 per cent). By stacking the near-IR spectra we directly measure the mean Balmer decrement for the H{\alpha} and H{\beta} lines, finding a value of <E(B-V)> = 0.51\pm0.27 for <LIR> = 10^12 Lsol sources at <z> = 1.36. By comparing star formation rates estimated from the IR and from the dust uncorrected H{\alpha} line we find a strong relationship between dust attenuation and star formation rate. This relation is broadly consistent with that previously seen in star-forming galaxies at z ~ 0.1. Finally, we investigate the metallicity via the N2 ratio, finding that z ~ 1 IR-selected sources are indistinguishable from the local mass-metallicity relation. We also find a strong correlation between dust attentuation and metallicity, with the most metal-rich IR-sources experiencing the largest levels of dust attenuation.

We present the results of wide-field deep JHK imaging of the SSA22 field using MOIRCS instrument equipped with Subaru telescope. The observed field is 112 arcmin^2 in area, which covers the z=3.1 protocluster characterized by the overdensities of Ly Alpha emitters (LAEs) and Ly Alpha Blobs (LABs). The 5 sigma limiting magnitude is K_{AB} = 24.3. We extract the potential protocluster members from the K-selected sample by using the multi-band photometric-redshift selection as well as the simple color cut for distant red galaxies (DRGs; J-K_{AB}>1.4). The surface number density of DRGs in our observed fields shows clear excess compared with those in the blank fields, and the location of the densest area whose projected overdensity is twice the average coincides with the large-scale density peak of LAEs. We also found that K-band counterparts with z_{phot} = 3.1 are detected for 75% (15/20) of the LABs within their Ly Alpha halo, and the 40 % (8/20) of LABs have multiple components, which gives a direct evidence of the hierarchical multiple merging in galaxy formation. The stellar mass ofLABs correlates with their luminosity, isophotal area, and the Ly Alpha velocity widths, implying that the physical scale and the dynamical motion of Ly Alpha emission are closely related to their previous star-formation activities. Highly dust-obscured galaxies such as hyper extremely red objects (HEROs; J-K_{AB}>2.1) and plausible K-band counterparts of submillimeter sources are also populated in the high density region.

We present spectroscopic and eleven-band photometric redshifts for galaxies in the 100-uJy Subaru/XMM-Newton Deep Field radio source sample. We find good agreement between our redshift distribution and that predicted by the SKA Simulated Skies project. We find no correlation between K-band magnitude and radio flux, but show that sources with 1.4-GHz flux densities below ~1mJy are fainter in the near-infrared than brighter radio sources at the same redshift, and we discuss the implications of this result for spectroscopically-incomplete samples where the K-z relation has been used to estimate redshifts. We use the infrared--radio correlation to separate our sample into radio-loud and radio-quiet objects and show that only radio-loud hosts have spectral energy distributions consistent with predominantly old stellar populations, although the fraction of objects displaying such properties is a decreasing function of radio luminosity. We calculate the 1.4-GHz radio luminosity function (RLF) in redshift bins to z=4 and find that the space density of radio sources increases with lookback time to z~2, with a more rapid increase for more powerful sources. We demonstrate that radio-loud and radio-quiet sources of the same radio luminosity evolve very differently. Radio-quiet sources display strong evolution to z~2 while radio-loud AGNs below the break in the radio luminosity function evolve more modestly and show hints of a decline in their space density at z>1, with this decline occurring later for lower-luminosity objects. If the radio luminosities of these sources are a function of their black hole spins then slowly-rotating black holes must have a plentiful fuel supply for longer, perhaps because they have yet to encounter the major merger that will spin them up and use the remaining gas in a major burst of star formation.

We present near-infrared spectroscopic observations of star-forming galaxies at z~1.4 with FMOS on the Subaru Telescope. We observed K-band selected galaxies in the SXDS/UDS fields with K<23.9 mag, 1.2<z_ph<1.6, M*>10^{9.5} Msun, and expected F(Halpha)>10^{-16} erg s^{-1} cm^{-2}. 71 objects in the sample have significant detections of Halpha. For these objects, excluding possible AGNs identified from the BPT diagram, gas-phase metallicities are obtained from [NII]/Halpha line ratio. The sample is split into three stellar mass bins, and the spectra are stacked in each stellar mass bin. The mass-metallicity relation obtained at z~1.4 is located between those at z~0.8 and z~2.2. We constrain an intrinsic scatter to be ~0.1 dex or larger in the mass-metallicity relation at z~1.4; the scatter may be larger at higher redshifts. We found trends that the deviation from the mass-metallicity relation depends on the SFR and the half light radius: Galaxies with higher SFR and larger half light radii show lower metallicities at a given stellar mass. One possible scenario for the trends is the infall of pristine gas accreted from IGM or through merger events. Our data points show larger scatter than the fundamental metallicity relation (FMR) at z~0.1 and the average metallicities slightly deviate from the FMR. The compilation of the mass-metallicity relations at z~3 to z~0.1 shows that they evolve smoothly from z~3 to z~0 without changing the shape so much except for the massive part at z~0.

The FIBRE-pac (FMOS image-based reduction package) is an IRAF-based reduction tool for the fiber multiple-object spectrograph (FMOS) of the Subaru telescope. To reduce FMOS images, a number of special techniques are necessary because each image contains about 200 separate spectra with airglow emission lines variable in spatial and time domains, and with complicated throughput patterns for the airglow masks. In spite of these features, almost all of the reduction processes except for a few steps are carried out automatically by scripts in text format making it easy to check the commands step by step. Wavelength- and flux-calibrated images together with their noise maps are obtained using this reduction package.

We study the evolution of quiescent galaxies at 0.5<z<2.5 as a function of stellar mass, using very deep NIR imaging data taken with the Multi-Object Infrared Camera and Spectrograph on the Subaru Telescope in the GOODS-North region. The deep NIR data allow us to construct a stellar mass-limited sample of quiescent galaxies down to ~10^{10} Msun even at z~2 for the first time. We selected quiescent galaxies with age/tau>6 by performing SED fitting of the multi broad-band photometry from the U to Spitzer 5.8um bands with the population synthesis model of Bruzual & Charlot (2003) where exponentially decaying star formation histories are assumed. The number density of quiescent galaxies increases by a factor of ~3 from 1.0<z<1.5 to 0.5<z<1.0, and by a factor of ~10 from 1.5<z<2.5 to 0.5<z<1.0, while that of star-forming galaxies with age/tau<4 increases only by factors of ~2 and ~3 in the same redshift ranges. At 0.5<z<2.5, the low-mass slope of the stellar mass function of quiescent galaxies is alpha ~ 0 -- 0.6, which is significantly flatter than those of star-forming galaxies (alpha ~ -1.3 -- -1.5). As a result, the fraction of quiescent galaxies in the overall galaxy population increases with stellar mass in the redshift range. The fraction of quiescent galaxies at 10^{11}-10^{11.5} Msun increases from ~20-30% at z~2 to ~40-60% at z~0.75, while that at 10^{10}-10^{10.5} Msun increases from <~ 5% to ~15% in the same redshift range. These results could suggest that the quenching of star formation had been more effective in more massive galaxies at 1<~z<~2. Such a mass-dependent quenching could explain the rapid increase of the number density of ~M* galaxies relative to lower-mass galaxies at z >~ 1-1.5.

We investigate mass-dependent galaxy evolution based on a large sample of (more than 50,000) K-band selected galaxies in a multi-wavelength catalog of the Subaru/XMM-Newton Deep Survey (SXDS) and the UKIRT Infrared Deep Sky Survey (UKIDSS)/Ultra Deep Survey (UDS). We employ the optical to near-infrared photometry to determine photometric redshifts of these galaxies. Then, we estimate the stellar mass of our sample galaxies using a standard fitting procedure. From the sample galaxies, we obtain the stellar mass function of galaxies and the cosmic stellar mass density up to z<4. Our results are consistent with previous studies and we find a considerable number of low-mass galaxies (M<10^{10.5}) at the redshift range 3<z<4. The stellar-mass dependent correlation functions of our sample galaxies show clear evolution and they connect to that in the local universe consistently. Also, the massive galaxies show strong clustering throughout our studied redshift range. The correlation length of massive galaxies rapidly decreases from z=4 to 2. We also find some high mass density regions of massive galaxies at 1.4<z<2.5 in our sample, which may be candidate progenitors of the present-day clusters of galaxies. At this redshift range, massive star-forming galaxies are the dominant population making up the structures and the passively evolving galaxies show stronger clustering and they may have formed earlier than those star-forming galaxies.

We present deep J-, H-, and Ks-band imaging data of the MOIRCS Deep Survey (MODS), which was carried out with Multi-Object Infrared Camera and Spectrograph (MOIRCS) mounted on the Subaru telescope in the GOODS-North region. The data reach 5sigma total limiting magnitudes for point sources of J=23.9, H=22.8, and Ks=22.8 (Vega magnitude) over 103 arcmin^2 (wide field). In 28 arcmin^2 of the survey area, which is ultra deep field of the MODS (deep field), the data reach the 5sigma depths of J=24.8, H=23.4, and Ks=23.8. The spatial resolutions of the combined images are FWHM ~ 0.6 arcsec and ~ 0.5 arcsec for the wide and deep fields in all bands, respectively. Combining the MODS data with the multi-wavelength public data taken with the HST, Spitzer, and other ground-based telescopes in the GOODS field, we construct a multi-wavelength photometric catalog of Ks-selected sources. Using the catalog, we present Ks-band number counts and near-infrared color distribution of the detected objects, and demonstrate some selection techniques with the NIR colors for high redshift galaxies. These data and catalog are publicly available via internet.

XMM and Chandra opened a new area for the study of clusters of galaxies. Not only for cluster physics but also, for the detection of faint and distant clusters that were inaccessible with previous missions. This article presents 66 spectroscopically confirmed clusters (0.05<z<1.5) within an area of 6 deg2 enclosed in the XMM-LSS survey. Almost two thirds have been confirmed with dedicated spectroscopy only and 10% have been confirmed with dedicated spectroscopy supplemented by literature redshifts. Sub-samples, or classes, of extended-sources are defined in a two-dimensional X-ray parameter space allowing for various degrees of completeness and contamination. We describe the procedure developed to assess the reality of these cluster candidates using the CFHTLS photometric data and spectroscopic information from our own follow-up campaigns. Most of these objects are low mass clusters, hence constituting a still poorly studied population. In a second step, we quantify correlations between the optical properties such as richness or velocity dispersion and the cluster X-ray luminosities. We examine the relation of the clusters to the cosmic web. Finally, we review peculiar structures in the surveyed area like very distant clusters and fossil groups.

(Abridged) Deep NB359 imaging with Subaru by Iwata et al. have detected surprisingly strong Lyman continuum (LyC; ~900A in the rest-frame) from some LAEs at z=3.1. However, the redshifts might be misidentified due to a narrow wavelength coverage in previous spectroscopy. We here present new deep spectroscopy covering the observed 4,000-7,000A with VLT/VIMOS and Subaru/FOCAS of 8 LAEs detected in NB359. All the 8 objects have only one detectable emission line around 4,970A which is most likely to be Ly-A at z=3.1, and thus, the objects are certainly LAEs at the redshift. However, 5 of them show a ~0.''8 spatial offset between the Ly-A emission and the source detected in NB359. No indications of the redshifts of the NB359 sources are found although it is statistically difficult that all the 5 LAEs have a foreground object accounting for the NB359 flux. The rest 3 LAEs show no significant offset from the NB359 position. Therefore, they are truly LyC emitting LAEs at z=3.1. We also examine the stellar population which simultaneously accounts for the strength of the LyC and the spectral slope of non-ionizing ultraviolet of the LAEs. We consider the latest statistics of Lyman limit systems to estimate the LyC optical depth in the IGM and an additional contribution of the bound-free LyC from photo-ionized nebulae to the LyC emissivity. As a result, we find that stellar populations with metallicity Z>=1/50Z_sun can explain the observed LyC strength only with a very top-heavy initial mass function (IMF; <m>~50 M_sun). However, the critical metallicity for such an IMF is expected to be much lower. A very young (~1 Myr) and massive (~100 M_sun) extremely metal-poor (Z<=5e-4Z_sun) or metal-free (so-called Population III) stellar population can reproduce the observed LyC strength. The required mass fraction of such `primordial' stellar population is ~1--10% in total stellar mass of the LAEs.

We investigate rest-frame near-infrared (NIR) morphologies of a sample of 139 galaxies with M_{s} >= 1 x 10^{10} M_{sun} at z=0.8-1.2 in the GOODS-North field using our deep NIR imaging data (MOIRCS Deep Survey, MODS). We focus on Luminous Infrared Galaxies (LIRGs), which dominate high star formation rate (SFR) density at z~1, in the sample identified by cross-correlating with the Spitzer/MIPS 24um source catalog. We perform two-dimensional light profile fitting of the z~1 galaxies in the Ks-band (rest-frame J-band) with a single component Sersic model. We find that at z~1, ~90% of LIRGs have low Sersic indices (n<2.5, similar to disk-like galaxies) in the Ks-band, and those disk-like LIRGs consist of ~60% of the whole disk-like sample above M_{s} >= 3 x 10^{10} M_{sun}. The z~1 disk-like LIRGs are comparable or ~20% small at a maximum in size compared to local disk-like galaxies in the same stellar mass range. If we examine rest-frame UV-optical morphologies using the HST/ACS images, the rest-frame B-band sizes of the z~1 disk-like galaxies are comparable to those of the local disk-like galaxies as reported by previous studies on size evolution of disk-like galaxies in the rest-frame optical band. Measuring color gradients (galaxy sizes as a function of wavelength) of the z~1 and local disk-like galaxies, we find that the z~1 disk-like galaxies have 3-5 times steeper color gradient than the local ones. Our results indicate that (i) more than a half of relatively massive disk-like galaxies at z~1 are in violent star formation epochs observed as LIRGs, and also (ii) most of those LIRGs are constructing their fundamental disk structure vigorously. The high SFR density in the universe at z~1 may be dominated by such star formation in disk region in massive galaxies.

We present Lya luminosity function (LF), clustering measurements, and Lya line profiles based on the largest sample, to date, of 207 Lya emitters (LAEs) at z=6.6 on the 1-deg^2 sky of Subaru/XMM-Newton Deep Survey (SXDS) field. Our z=6.6 Lya LF including cosmic variance estimates yields the best-fit Schechter parameters of phi*=8.5 +3.0/-2.2 x10^(-4) Mpc^(-3) and L*(Lya)=4.4 +/-0.6 x10^42 erg s^(-1) with a fixed alpha=-1.5, and indicates a decrease from z=5.7 at the >~90% confidence level. However, this decrease is not large, only =~30% in Lya luminosity, which is too small to be identified in the previous studies. A clustering signal of z=6.6 LAEs is detected for the first time. We obtain the correlation length of r_0=2-5 h^(-1) Mpc and bias of b=3-6, and find no significant boost of clustering amplitude by reionization at z=6.6. The average hosting dark halo mass inferred from clustering is 10^10-10^11 Mo, and duty cycle of LAE population is roughly ~1% albeit with large uncertainties. The average of our high-quality Keck/DEIMOS spectra shows an FWHM velocity width of 251 +/-16 km s^(-1). We find no large evolution of Lya line profile from z=5.7 to 6.6, and no anti-correlation between Lya luminosity and line width at z=6.6. The combination of various reionization models and our observational results about the LF, clustering, and line profile indicates that there would exist a small decrease of IGM's Lya transmission owing to reionization, but that the hydrogen IGM is not highly neutral at z=6.6. Our neutral-hydrogen fraction constraint implies that the major reionization process took place at z>~7.

We study the evolution of star formation activity of galaxies at 0.5<z<3.5 as a function of stellar mass, using very deep NIR data taken with Multi-Object Infrared Camera and Spectrograph (MOIRCS) on the Subaru telescope in the GOODS-North region. The NIR imaging data reach K ~ 23-24 Vega magnitude and they allow us to construct a nearly stellar mass-limited sample down to ~ 10^{9.5-10} Msun even at z~3. We estimated star formation rates (SFRs) of the sample with two indicators, namely, the Spitzer/MIPS 24um flux and the rest-frame 2800A luminosity. The SFR distribution at a fixed Mstar shifts to higher values with increasing redshift at 0.5<z<3.5. More massive galaxies show stronger evolution of SFR at z>~1. We found galaxies at 2.5<z<3.5 show a bimodality in their SSFR distribution, which can be divided into two populations by a constant SSFR of ~2 Gyr^{-1}. Galaxies in the low-SSFR group have SSFRs of ~ 0.5-1.0 Gyr^{-1}, while the high-SSFR population shows ~10 Gyr^{-1}. The cosmic SFRD is dominated by galaxies with Mstar = 10^{10-11} Msun at 0.5<z<3.5, while the contribution of massive galaxies with Mstar = 10^{11-11.5} Msun shows a strong evolution at z>1 and becomes significant at z~3, especially in the case with the SFR based on MIPS 24um. In galaxies with Mstar = 10^{10-11.5} Msun, those with a relatively narrow range of SSFR (<~1 dex) dominates the cosmic SFRD at 0.5<z<3.5. The SSFR of galaxies which dominate the SFRD systematically increases with redshift. At 2.5<z<3.5, the high-SSFR population, which is relatively small in number, dominates the SFRD. Major star formation in the universe at higher redshift seems to be associated with a more rapid growth of stellar mass of galaxies.

Fibre Multi-Object Spectrograph (FMOS) is the first near-infrared instrument with a wide field of view capable of acquiring spectra simultaneously from up to 400 objects. It has been developed as a common-use instrument for the F/2 prime-focus of the Subaru Telescope. The field coverage of 30' diameter is achieved using a new 3-element corrector optimized in the near-infrared (0.9-1.8um) wavelength range. Due to limited space at the prime-focus, we have had to develop a novel fibre positioner called "Echidna" together with two OH-airglow suppressed spectrographs. FMOS consists of three subsystems: the prime focus unit for IR, the fibre positioning system/connector units, and the two spectrographs. After full systems integration, FMOS was installed on the telescope in late 2007. Many aspects of performance were checked through various test and engineering observations. In this paper, we present the optical and mechanical components of FMOS and show the results of our on-sky engineering observations to date.

We present near-infrared (NIR) spectroscopy for 37 BzK-color-selected star-forming galaxies with MOIRCS on Subaru. The sample is drawn from the Ks-band selected catalog of the MOIRCS Deep Survey (MODS) in the GOODS-N region. About half of our samples are selected from the MIPS 24um-source catalog. Ha emission lines are detected from 23 galaxies, of which the median redshift is 2.12. We derived the star formation rates (SFRs) from extinction-corrected Ha luminosities. The extinction correction is estimated from the SED fitting of multi-band photometric data covering UV to NIR wavelengths. The Balmer decrement of the stacked emission lines shows that the amount of extinction for the ionized gas is larger than that for the stellar continuum. From a comparison of the extinction corrected Ha luminosity and other SFR indicators we found that the relation between the dust properties of stellar continuum and ionized gas is different depending on the intrinsic SFR. The comparison between the Ha SFR and stellar mass estimated from SED fitting shows no correlation between them. Some galaxies with stellar mass smaller than ~10^10 Msun show SFRs higher than ~100 Msun/yr. The specific SFRs (SSFRs) of these galaxies are remarkably high. The average stellar-population age of these high-SSFR galaxies is younger than 100 Myr, which is consistent with the implied high SSFR. The large SFR implies the possibility that the high SSFR galaxies significantly contribute to the cosmic SFR density of the universe at z~2. The total SFR density estimated from the Ha emission line galaxies is 0.089-0.136 Msun/yr/Mpc^3, which is consistent with the total SFR densities in the literature. The metallicity of the high-SSFR galaxies is larger than that expected from the mass-metallicity relation of UV-selected galaxies at z~2 by Erb et al. (2006a). [Abridged]

We present an 880 micron Submillimeter Array (SMA) detection of the submillimeter galaxy SXDF850.6. SXDF850.6 is a bright source (S(850 micron) = 8 mJy) detected in the SCUBA Half Degree Extragalactic Survey (SHADES), and has multiple possible radio counterparts in its deep radio image obtained at the VLA. Our new SMA detection finds that the submm emission coincides with the brightest radio emission that is found ~8" north of the coordinates determined from SCUBA. Despite the lack of detectable counterparts in deep UV/optical images, we find a source at the SMA position in near-infrared and longer wavelength images. We perform SED model fits to UV-optical-IR photometry (u, B, V, R, i', z', J, H, K, 3.6 micron, 4.5 micron, 5.8 micron, and 8.0 micron) and to submm-radio photometry (850 micron, 880 micron, 1100 micron, and 21 cm) independently, and we find both are well described by starburst templates at a redshift of z ~= 2.2 (+/- 0.3). The best-fit parameters from the UV-optical-IR SED fit are a redshift of z = 1.87 (+0.15/-0.07), a stellar mass of M_star = 2.5 +2.2/-0.3 x 10^11 M_sun, an extinction of A_V = 3.0 (+0.3/-1.0) mag, and an age of 720 (+1880/-210) Myr. The submm-radio SED fit provides a consistent redshift of z ~ 1.8-2.5, an IR luminosity of L_IR = (7-26) x 10^12 L_sun, and a star formation rate of 1300-4500 M_sun/yr. These results suggest that SXDF850.6 is a mature system already having a massive amount of old stellar population constructed before its submm bright phase and is experiencing a dusty starburst, possibly induced by major mergers.

We present a universal linear correlation between the stellar mass and surface brightness (SB) of galaxies at 0.3<z<3, using a deep K-band selected catalog in the GOODS-North region. The correlation has a nearly constant slope, independent of redshift and color of galaxies in the rest-z frame. Considering unresolved compact galaxies, the tight correlation gives a lower boundary of SB for a given stellar mass; lower SB galaxies are prohibited over the boundary. The universal slope suggests that the stellar mass in galaxies was build up over their cosmic histories in a similar manner irrelevant to galaxy mass, as oppose to the scenario that massive galaxies mainly accumulated their stellar mass by major merging. In contrast, SB shows a strong dependence on redshift for a given stellar mass. It evolves as (1+z)^(-2.0~-0.8), in addition to dimming as (1+z)^4 by the cosmological expansion effect. The brightening depends on galaxy color and stellar mass. The blue population (rest-frame U-V<0), which is dominated by young and star-forming galaxies, evolves as ~(1+z)^(-0.8 +-0.3) in the rest-V band. On the other hand, the red population (U-V>0) and the massive galaxies (M_*>10^(10)M_sun) shows stronger brightening, (1+z)^(-1.5+-0.1). Based on the comparison with galaxy evolution models, we find that the phenomena are well explained by the pure luminosity evolution of galaxies out to z~3.

We present the results of a search for galaxy clusters in Subaru-XMM Deep Field. We reach a depth for a total cluster flux in the 0.5-2 keV band of 2x10^{-15} ergs cm^{-2} s^{-1} over one of the widest XMM-Newton contiguous raster surveys, covering an area of 1.3 square degrees. Cluster candidates are identified through a wavelet detection of extended X-ray emission. The red sequence technique allows us to identify 57 cluster candidates. We report on the progress with the cluster spectroscopic follow-up and derive their properties based on the X-ray luminosity and cluster scaling relations. In addition, 3 sources are identified as X-ray counterparts of radio lobes, and in 3 further sources, X-ray counterpart of radio lobes provides a significant fraction of the total flux of the source. In the area covered by NIR data, our identification success rate achieves 86%. We detect a number of radio galaxies within our groups and for a luminosity-limited sample of radio galaxies we compute halo occupation statistics using a marked cluster mass function. We compare the cluster detection statistics in the SXDF with the predictions of concordance cosmology and current knowledge of the X-ray cluster properties, concluding that a reduction of concordance sigma_8 value by 5% is required in order to match the prediction of the model and the data. This conclusion still needs verification through the completion of cluster follow-up.

We investigate the stellar populations of Lyman alpha emitters (LAEs) at z=3.1 and 3.7 in 0.65 deg^2 of the Subaru/XMM-Newton Deep Field, based on rest-frame UV-to-optical photometry obtained from the Subaru/XMM-Newton Deep Survey, the UKIDSS/Ultra Deep Survey, and the Spitzer legacy survey of the UKIDSS/UDS. Among a total of 302 LAEs (224 for z=3.1 and 78 for z=3.7), only 11 are detected in the K band, i.e., brighter than K(3sigma)=24.1 mag. Eight of the 11 K-detected LAEs are spectroscopically confirmed. We find that the K-undetected objects, which should closely represent the LAE population as a whole, have low stellar masses of ~ 10^8 - 10^8.5 Msun, modest SFRs of 1 - 100 Msun yr^-1, and modest dust extinction of E(B-V) < 0.2. The K-detected objects are massive, Mstar ~ 10^9 - 10^10.5 Msun, and have significant dust extinction with a median of E(B-V) ~= 0.3. Four K-detected objects with the reddest spectral energy distributions, two of which are spectroscopically confirmed, are heavily obscured with E(B-V) ~ 0.65, and their continua resemble those of some local ULIRGs. Interestingly, they have large Lyman alpha equivalent widths ~= 70 - 250 A. If these four are excluded, our sample has a weak anti-correlation between EW(Lya) and Mstar. We compare the stellar masses and the specific star formation rates (sSFR) of LAEs with those of Lyman-break galaxies, distant red galaxies, submillimetre galaxies, and I- or K-selected galaxies with z_phot ~ 3. We find that the LAE population is the least massive among all the galaxy populations in question, but with relatively high sSFRs, while NIR-detected LAEs have Mstar and sSFR similar to LBGs. Our reddest four LAEs have very high sSFRs in spite of large Mstar, thus occupying a unique region in the Mstar versus sSFR space. (abridged)

We present results of optical spectroscopic observations of candidates of Lyman Break Galaxies (LBGs) at $z \sim 5$ in the region including the GOODS-N and the J0053+1234 region by using GMOS-N and GMOS-S, respectively. Among 25 candidates, five objects are identified to be at $z \sim 5$ (two of them were already identified by an earlier study) and one object very close to the color-selection window turned out to be a foreground galaxy. With this spectroscopically identified sample and those from previous studies, we derived the lower limits on the number density of bright ($M_{UV}<-22.0$ mag) LBGs at $z \sim 5$. These lower limits are comparable to or slightly smaller than the number densities of UV luminosity functions (UVLFs) that show the smaller number density among $z \sim 5$ UVLFs in literature. However, by considering that there remain many LBG candidates without spectroscopic observations, the number density of bright LBGs is expected to increase by a factor of two or more. The evidence for the deficiency of UV luminous LBGs with large Ly$\alpha$ equivalent widths was reinforced. We discuss possible causes for the deficiency and prefer the interpretation of dust absorption.

We use very deep near-infrared (NIR) imaging data obtained in MOIRCS Deep Survey (MODS) to investigate the evolution of the galaxy stellar mass function back to z~3. The MODS data reach J=24.2, H=23.1, K=23.1 (5sigma, Vega magnitude) over 103 arcmin^2 (wide) and J=25.1, H=23.7, K=24.1 over 28 arcmin^2 (deep) in the GOODS-North region. The wide and very deep NIR data allow us to measure the number density of galaxies down to low stellar mass (10^9-10^10 Msun) even at high redshift with high statistical accuracy. The normalization of the mass function decreases with redshift and the integrated stellar mass density becomes ~ 8-18% of the local value at z~2 and ~ 4-9% at z~3, which are consistent with results of previous studies in general fields. Furthermore, we found that the low-mass slope becomes steeper with redshift from alpha ~- 1.3 at z~1 to alpha ~- 1.6 at z~3, and that the evolution of the number density of low-mass (10^9-10^10 Msun) galaxies is weaker than that of M* (~10^11 Msun) galaxies. This indicates that the contribution of low-mass galaxies to the total stellar mass density has been significant at high redshift. The steepening of the low-mass slope with redshift is opposite trend expected from the stellar mass dependence of the specific star formation rate reported in previous studies. The present result suggests that the hierarchical merging process overwhelmed the effect of the stellar mass growth by star formation and was very important for the stellar mass assembly of these galaxies at 1<~z<~3.

We investigate the X-ray properties of the K-band-selected galaxies at redshift 2 < z < 4 by using our deep near-infrared images obtained in the MOIRCS Deep Survey project and the published Chandra X-ray source catalog. 61 X-ray sources with the 2-10 keV luminosity L_X = 10^{42}-10^{44} erg/s are identified with the K-selected galaxies and we found that they are exclusively (90%) associated with the massive objects with stellar mass larger than 10^{10.5} Msun. Our results are consistent with the idea that the M_BH/M_str ratio of the galaxies at z=2-4 is similar to the present-day value. On the other hand, the AGN detection rate among the very massive galaxies with the stellar mass larger than 10^{11} Msun is high, 33% (26/78). They are active objects in the sense that the black-hole mass accretion rate is ~ 1-50% of the Eddington limit if they indeed have similar M_BH/M_str ratio with those observed in the local universe. The active duration in the AGN duty cycle of the high-redshift massive galaxies seems large.

Knowing the amount of ionizing photons from young star-forming galaxies is of particular importance to understanding the reionization process. Here we report initial results of Subaru/Suprime-Cam deep imaging observation of the SSA22 proto-cluster region at z=3.09, using a special narrow-band filter to optimally trace Lyman continuum (LyC) from galaxies at z~3. The unique wide field-of-view of Suprime-Cam enabled us to search for ionizing photons from 198 galaxies (73 Lyman break galaxies (LBGs) and 125 Ly-alpha emitters (LAEs)) with spectroscopically measured redshifts z~3.1. We detected LyC from 7 LBGs, as well as from 10 LAE candidates. Some of the detected galaxies show significant spatial offsets of LyC from non-ionizing UV emission. For some LBGs the observed non-ionizing UV to LyC flux density ratios are smaller than values expected from population synthesis models with a standard Salpeter initial mass function (IMF) with moderate dust attenuation (which is suggested from the observed UV slopes), even if we assume very transparent IGM along the sightlines of these objects. This implies an intrinsically bluer spectral energy distribution, e.g, that produced by a top-heavy IMF, for these LBGs. The observed flux desity ratios of non-ionizing UV to LyC of 7 detected LBGs range from 2.4 to 23.8 and the median is 6.6. The observed flux density ratios of the detected LAEs are even smaller than LBGs, if they are truly at z~3.1. We find that the median value of the flux density ratio for the deteced LBGs suggest that their escape fractions is likely to be higher than 4%, if the LyC escape is isotropic. The results imply that some of the LBGs in the proto-cluster at z~3 have the escape fraction significantly higher than that of galaxies (in a general field) at z~1 studied previously.

We report the discovery of a giant Lya emitter (LAE) with a Spitzer/IRAC counterpart near the reionization epoch at z=6.595. The giant LAE is found from the extensive 1 deg^2 Subaru narrow-band survey for z=6.6 LAEs in the Subaru/XMM-Newton Deep Survey (SXDS) field, and subsequently identified by deep spectroscopy of Keck/DEIMOS and Magellan/IMACS. Among our 207 LAE candidates, this LAE is not only the brightest narrow-band object with L(Lya) = 3.9+/-0.2 x 10^43 erg/s in our survey volume of 10^6 Mpc^3, but also a spatially extended Lya nebula with the largest isophotal area whose major axis is at least ~3''. This object is more likely to be a large Lya nebula with a size of >~17 kpc than to be a strongly-lensed galaxy by a foreground object. Our Keck spectrum with medium-high spectral and spatial resolutions suggests that the velocity width is v(FWHM)=251+/-21 km/s, and that the line-center velocity changes by ~60 km/s in a 10-kpc range. The stellar mass and star-formation rate are estimated to be 0.9-5.0 x 10^10 Mo and >34 Mo/yr, respectively, from the combination of deep optical to infrared images of Subaru, UKIDSS-Ultra Deep Survey, and Spitzer/IRAC. Although the nature of this object is not yet clearly understood, this could be an important object for studying cooling clouds accreting onto a massive halo, or forming-massive galaxies with significant outflows contributing to cosmic reionization and metal enrichment of inter-galactic medium.

Efficient selection of emission line galaxies at z > 1 by photometric information in wide field surveys is one of the keys for future spectroscopic surveys to constrain dark energy using the baryon acoustic oscillation (BAO) signature. Here we estimate the H alpha and [O II] line luminosity functions of galaxies at z = 0.5-1.7 using a novel approach where multi-wavelength imaging data is used to jointly estimate both photometric redshifts and star-formation rates. These photometric estimates of line luminosities at high-redshift use the large data sets of the Subaru Deep Field and Subaru XMM-Newton Deep Field (covering \sim 1 deg^2) and are calibrated with the spectroscopic data of the local Sloan Digital Sky Survey galaxies. The derived luminosity functions (especially H alpha) are in reasonable agreement with the past estimates based on spectroscopic or narrow-band-filter surveys. This dataset is useful for examining the photometric selection of target galaxies for BAO surveys because of the large cosmological volume covered and the large number of galaxies with detailed photometric information. We use the sample to derive the photometric and physical properties of emission line galaxies to assist planning for future spectroscopic BAO surveys. We also show some examples of photometric selection procedures which can efficiently select these emission line galaxies.

We discuss the optical properties, X-ray detections, and Active Galactic Nucleus (AGN) populations of four clusters at z ~ 1 in the Subaru-XMM Deep Field (SXDF). The velocity distribution and plausible extended X-ray detections are examined, as well as the number of X-ray point sources and radio sources associated with the clusters. We find that the two clusters that appear virialised and have an extended X-ray detection contain few, if any, AGN, whereas the two pre-virialised clusters have a large AGN population. This constitutes evidence that the AGN fraction in clusters is linked to the clusters' evolutionary stage. The number of X-ray AGN in the pre-virialised clusters is consistent with an overdensity of factor ~ 200; the radio AGN appear to be clustered with a factor of three to six higher. The median K-band luminosities of L_K = 1.7 +/- 0.7 L* for the X-ray sources and L_K = 2.3 +/- 0.1 L* for the radio sources support the theory that these AGN are triggered by galaxy interaction and merging events in sub-groups with low internal velocity distributions, which make up the cluster environment in a pre-virialisation evolutionary stage.

We present the results of Spectral Energy Distribution(SED) fitting analysis for Lyman Break Galaxies(LBGs) at z~5 in the GOODS-N and its flanking fields (the GOODS-FF). With the publicly available IRAC images in the GOODS-N and IRAC data in the GOODS-FF, we constructed the rest-frame UV to optical SEDs for a large sample (~100) of UV-selected galaxies at z~5. Comparing the observed SEDs with model SEDs generated with a population synthesis code, we derived a best-fit set of parameters (stellar mass, age, color excess, and star formation rate) for each of sample LBGs. The derived stellar masses range from 10^8 to 10^11M_sun with a median value of 4.1x10^9M_sun. The comparison with z=2-3 LBGs shows that the stellar masses of z~5 LBGs are systematically smaller by a factor of 3-4 than those of z=2-3 LBGs in a similar rest-frame UV luminosity range. The star formation ages are relatively younger than those of the z=2-3 LBGs. We also compared the results for our sample with other studies for the z=5-6 galaxies. Although there seem to be similarities and differences in the properties, we could not conclude its significance. We also derived a stellar mass function of our sample by correcting for incompletenesses. Although the number densities in the massive end are comparable to the theoretical predictions from semi-analytic models, the number densities in the low-mass part are smaller than the model predictions. By integrating the stellar mass function down to 10^8 M_sun, the stellar mass density at z~5 is calculated to be (0.7-2.4)x10^7M_sun Mpc^-3. The stellar mass density at z~5 is dominated by massive part of the stellar mass function. Compared with other observational studies and the model predictions, the mass density of our sample is consistent with general trend of the increase of the stellar mass density with time.

We aim to understand the multi-wavelength properties of 2XMM J123204+215255, the source with the most extreme X-ray-to-optical flux ratio amongst a sample of bright X-ray selected EXOs drawn from a cross-correlation of the 2XMMp catalogue with the SDSS-DR5 catalogue. We use 2XMMp X-ray data, SDSS-DR5, NOT and UKIRT optical/NIR photometric data and Subaru MOIRCS IR spectroscopy to study the properties of 2XMM J123204+215255. We created a model SED including an obscured QSO and the host galaxy component to constrain the optical/IR extinction and the relative contribution of the AGN and the galaxy to the total emission. 2XMM J123204+215255 is a bright X-ray source with f_X~10^{-12} erg cm^{-2} s^{-1} (2-10 keV energy band) which has no detection down to a magnitude i' > 25.2. NIR imaging reveals a faint K-band counterpart and NIR spectroscopy shows a single broad (FWHM=5300 km/s) emission line, which is almost certainly H-alpha at z=1.87. The X-ray spectrum shows evidence of significant absorption (N_H > 10^{23} cm^{-2}), typical of type 2 AGN, but the broad H-alpha emission suggests a type 1 AGN classification. The very red optical/NIR colours (i'-K > 5.3) strongly suggest significant reddening however. We find that simple modelling can successfully reproduce the NIR continuum and strongly constrain the intrinsic nuclear optical/IR extinction to A_V~4, which turns out to be much smaller than the expected from the X-ray absorption (assuming Galactic gas-to-dust ratio).

We present the results of Spectral Energy Distribution (SED) fitting analysis for Lyman Break Galaxies (LBGs) at $z\sim5$ in the GOODS-N and its flanking fields. From the SED fitting for $\sim100$ objects, we found that the stellar masses range from $10^{8}$ to $10^{11}M_{\odot}$ with a median value of $4\times10^{9}M_{\odot}$. By using the large sample of galaxies at $z\sim5$, we construct the stellar mass function (SMF) with incompleteness corrections. By integrating down to $10^{8}M_{\odot}$, the cosmic stellar mass density at $z\sim5$ is calculated to be $7\times10^{6}M_{\odot}\textrm{Mpc}^{-3}$.

Properties of Lyman Break Galaxies at $z\sim5$ are presented and are compared with those at $z\sim 2-3$ LBGs to see evolution of LBGs from $z\sim5$ to $z \sim3-2$. The LBGs at $z\sim5$ are tend to have smaller stellar masses, but they may show large star formation rates and specific star formation rates.

We briefly summarize our findings from the unbiased surveys for $z$$\sim$5 LBGs based on Subaru/Suprime-Cam and follow-up optical spectroscopy.

We present the X-ray source catalog in the Subaru/XMM-Newton deep survey. A continuous area of 1.14 deg^2 centered at R.A. = 02h18m and Dec. = -05d is mapped by seven pointings with XMM-Newton covering the 0.2-10 keV band. From the combined images of the EPIC pn and MOS cameras, we detect 866, 1114, 645, and 136 sources with sensitivity limits of 6x10^{-16}, 8x10^{-16}, 3x10^{-15}, and 5x10^{-15} erg cm^{-2} s^{-1} in the 0.5-2, 0.5-4.5, 2-10, and 4.5-10 keV bands, respectively, with detection likelihood >= 7 (corresponding to a confidence level of 99.91%). The catalog consists of 1245 sources in total including 32 extended-source candidates. The averaged log N-log S relations are in good agreement with previous results, bridging the flux range between Chandra deep surveys and brighter surveys. The log N-log S relations show significant spatial variation among pointings on a scale of 0.2 deg^2. Analyzing the auto correlation function, we detect significant clustering signals from the 0.5-2 keV band sample, which can be fit with a power law form (\theta/\theta_c)^{-0.8} with a correlation length of \theta_c=5.9^{+1.0}_{-0.9} arcsec when the integral constraint term is included. In the 2-10 keV band, however, the clustering is not significant with a 90% upper limit of \theta_c < 1.5 arcsec.

We have undertaken a pilot survey for faint QSOs in the UKIDSS Ultra Deep Survey Field using the KX selection technique. These observations exploit the very deep near-infrared and optical imaging of this field from UKIRT and Subaru to select candidate QSOs based on their VJK colours and morphologies. We determined redshifts for 426 candidates using the AAOmega spectrograph on the AAT in service time. We identify 17 QSOs (M_B<= -23) in this pilot survey at z=1.57-3.29. We combine our sample with an X-ray selected sample of QSOs in the same field (a large fraction of which also comply with our KX selection) to constrain the surface density of QSOs with K<=20, deriving limits on the likely surface density of 85-150/deg^2. We use the good image quality available from our near-infrared imaging to detect a spatially extended component of the QSO light which probably represents the host galaxies. We also use our sample to investigate routes to improve the selection of KX QSOs at faint limits in the face of the significant contamination by compact, foreground galaxies. The brightest examples from our combined QSO sample will be used in conjunction with a large VLT VIMOS spectroscopic survey of high redshift galaxies in this region to study the structures inhabited by gas, galaxies and growing super-massive black holes at high redshifts in the UKIDSS UDS.

We present luminosity functions (LFs) and various properties of Lya emitters (LAEs) at z=3.1, 3.7, and 5.7, in a 1 deg^2 sky of the Subaru/XMM-Newton Deep Survey (SXDS) Field. We obtain a photometric sample of 858 LAE candidates based on deep Subaru/Suprime-Cam imaging data, and a spectroscopic sample of 84 confirmed LAEs from Subaru/FOCAS and VLT/VIMOS spectroscopy in a survey volume of ~10^6 Mpc^3 with a limiting Lya luminosity of ~3x10^42 erg/s. We derive the LFs of Lya and UV-continuum (~1500 \AA) for each redshift, taking into account the statistical error and the field-to-field variation. We find that the apparent Lya LF shows no significant evolution between z=3.1 and 5.7 within factors of 1.8 and 2.7 in L* and phi*, respectively. On the other hand, the UV LF of LAEs increases from z=3.1 to 5.7, indicating that galaxies with Lya emission are more common at earlier epochs. We identify six LAEs with AGN activities from our spectra combined with VLA, Spitzer, and XMM-Newton data. Among the photometrically selected LAEs at z=3.1 and 3.7, only ~1 % show AGN activities, while the brightest LAEs with logL(Lya) >~ 43.4-43.6 erg/s appear to always host AGNs. Our LAEs are bluer in UV-continuum color than dropout galaxies, suggesting lower extinction and/or younger stellar populations. Our stacking analyses provide upper limits to the radio luminosity and the f(HeII)/f(Lya) line fraction, and constrain the hidden star formation (+low-luminosity AGN) and the primordial population in LAEs.

We present the results of deep near-infrared imaging observations of the z=3.1 proto-cluster region in the SSA22a field taken by MOIRCS mounted on the Subaru Telescope. We observed a 21.7 arcmin^2 field to the depths of J=24.5, H=24.3, and K=23.9 (5 sigma). We examine the distribution of the K-selected galaxies at z~3 by using the simple color cut for distant red galaxies (DRGs) as well as the photometric-redshift selection technique. The marginal density excess of DRGs and the photo-z selected objects are found around the two most luminous Ly alpha blobs (LABs). We investigate the correlation between the K-selected objects and the LABs, and find that several galaxies with stellar mass M_* = 10^9-10^11 M_solar exist in vicinity of LABs, especially around the two most luminous ones. We also find that 7 of the 8 LABs in the field have plausible K_s-band counterparts and the sum of the stellar mass possibly associated with LABs correlates with the luminosity and surface brightness of them, which implies that the origin of Ly alpha emission may be closely correlated with their stellar mass or their previous star formation phenomena.

We present multi-waveband optical imaging data obtained from observations of the Subaru/XMM-Newton Deep Survey (SXDS). The survey field, centered at R.A.=02:18:00, decl.=-05:00:00, has been the focus of a wide range of multi-wavelength observing programs spanning from X-ray to radio wavelengths. A large part of the optical imaging observations are carried out with Suprime-Cam on Subaru Telescope at Mauna Kea in the course of Subaru Telescope Observatory Projects. This paper describes our optical observations, data reduction and analysis procedures employed, and the characteristics of the data products. A total area of 1.22 sqdeg is covered in five contiguous sub-fields, each of which corresponds to a single Suprime-Cam field of view (34'x27'), in five broad-band filters B, V, Rc, i', z' to the depths of B=28.4, V=27.8, Rc=27.7, i'=27.7 and z'=26.6 (AB, 3-sigma, 2-arcsec aperture). The data are reduced and compiled into five multi-waveband photometric catalogs, separately for each Suprime-Cam pointing. The i'-band catalogs contain about 900,000 objects, making the SXDS catalogs one of the largest multi-waveband catalogs in corresponding depth and area coverage. The SXDS catalogs can be used for an extensive range of astronomical applications such as the number density of the Galactic halo stars to the large scale structures at the distant universe. The number counts of galaxies are derived and compared with those of existing deep extragalactic surveys. The optical data, the source catalogs, and configuration files used to create the catalogs are publicly available via the SXDS web page (http://www.naoj.org/Science/SubaruProject/SXDS/index.html)

We present the properties of active galactic nuclei (AGN) selected by optical variability in the Subaru/XMM-Newton Deep Field (SXDF). Based on the locations of variable components and light curves, 211 optically variable AGN were reliably selected. We made three AGN samples; X-ray detected optically non-variable AGN (XA), X-ray detected optically variable AGN (XVA), and X-ray undetected optically variable AGN (VA). In the VA sample, we found a bimodal distribution of the ratio between the variable component flux and the host flux. One of these two components in the distribution, a class of AGN with a faint variable component $i'_{\rm{vari}}\sim25$ mag in bright host galaxies $i'\sim21$ mag, is not seen in the XVA sample. These AGN are expected to have low Eddington ratios if we naively consider a correlation between bulge luminosity and black hole mass. These galaxies have photometric redshifts $z_{\rm{photo}}\sim0.5$ and we infer that they are low-luminosity AGN with radiatively inefficient accretion flows (RIAFs). The properties of the XVA and VA objects and the differences from those of the XA objects can be explained within the unified scheme for AGN. Optical variability selection for AGN is an independent method and could provide a complementary AGN sample which even deep X-ray surveys have not found.

We present our survey for optically faint variable objects using multi-epoch (8-10 epochs over 2-4 years) $i'$-band imaging data obtained with Subaru Suprime-Cam over 0.918 deg$^2$ in the Subaru/XMM-Newton Deep Field (SXDF). We found 1040 optically variable objects by image subtraction for all the combinations of images at different epochs. This is the first statistical sample of variable objects at depths achieved with 8-10m class telescopes or HST. The detection limit for variable components is $i'_{\rm{vari}}\sim25.5$ mag. These variable objects were classified into variable stars, supernovae (SNe), and active galactic nuclei (AGN), based on the optical morphologies, magnitudes, colors, and optical-mid-infrared colors of the host objects, spatial offsets of variable components from the host objects, and light curves. Detection completeness was examined by simulating light curves for periodic and irregular variability. We detected optical variability for $36\pm2%$ ($51\pm3%$ for a bright sample with $i'<24.4$ mag) of X-ray sources in the field. Number densities of variable obejcts as functions of time intervals $\Delta{t}$ and variable component magnitudes $i'_{\rm{vari}}$ are obtained. Number densities of variable stars, SNe, and AGN are 120, 489, and 579 objects deg$^{-2}$, respectively. Bimodal distributions of variable stars in the color-magnitude diagrams indicate that the variable star sample consists of bright ($V\sim22$ mag) blue variable stars of the halo population and faint ($V\sim23.5$ mag) red variable stars of the disk population. There are a few candidates of RR Lyrae providing a possible number density of $\sim10^{-2}$ kpc$^{-3}$ at a distance of $>150$ kpc from the Galactic center.

In order to reveal the stellar mass distribution of z~3 galaxies, we are conducting deep imaging observations of U-dropout Lyman Break Galaxies (LBGs) with Adaptive Optics (AO) systems in K-band, which corresponds to rest-frame V-band of z~3 galaxies. The results of the Subaru intensive-program observations with AO36/NGS/IRCS indicate that 1) the K-band peaks of some of the LBGs brighter than K=22.0 mag show significant offset from those in the optical images, 2) the z~3 Mv* LBGs and serendipitously observed Distant Red Galaxies (DRGs) have flat profiles similar to disk galaxies in the local universe (i.e., Sersic with n<2), and 3) the surface stellar mass densities of the Mv* LBGs are 3-6 times larger than those of disk galaxies at z=0-1. Considering the lack of n>2 systems among the luminous z~3 LBGs and DRGs, and their strong spatial clustering, we infer that the dense n<2 disk-like structures evolve into the n>2 spheroids of nearby galaxies through relaxations due to major merger events.

In order to reveal the rest-frame V-band morphology of galaxies at z~3, we conducted AO-assisted K-band imaging observations of z~3 LBGs with Mv*-0.5 to Mv*+3.0 mag. LBGs brighter than Mv* have larger rHL (0.40") than the fainter LBGs (0.23") on average, and there is no bright LBGs with a small rHL. The LBGs brighter than Mv* have red rest-frame U-V colors (average of 0.2 mag) and most of the fainter LBGs show blue rest-frame U-V colors (average of -0.4 mag). The K-band peaks of some of the LBGs brighter than K=22.0 mag show significant shift from those in the optical images. The images of all but one of the LBGs with K<21.5 mag are fitted well with Sersic profile with n index less than 2, similar to disk galaxies in the local universe. Assuming that the LBGs have a disk-shape, we compared their size-luminosity and size-stellar mass relation with those of z=0 and z=1 disk galaxies. The LBGs are brighter than z=0 and z=1 disk galaxies at the same effective radius. The rest-frame V-band surface brightness of the LBGs are 2.2-2.9 mag and 1.2-1.9 mag brighter than the disk galaxies at z=0 and z=1, respectively. The size-stellar mass relation of the LBGs shows that the effective radii of the LBGs do not depend on their stellar mass. For the LBGs brighter than Mv*, the average surface stellar mass density is 3-6 times larger than those of the z=0 and z=1 disk galaxies. We also examine the profiles of the serendipitously observed DRGs. They are also fitted with the Sersic profiles with n<2. The implications of the dominance of n<2 population among galaxies at z~3 and the presence of the high surface stellar mass density disk systems are discussed.

Spatially extended Ly-alpha sources that are faint and/or compact in coninuum are candidates for extremely young (~< 10^7 yrs) galaxies at high redshifts. We present medium-resolution (R~2000) spectroscopy of such extended Ly-alpha sources found in our previous study at z~3-5, using VLT/VIMOS. The deep spectroscopy showed that all 18 objects we observed have large equivalent widths (EWs) exceeding 100 A. For about 30% of our sample (five objects), we identified conspicuous asymmetry on the profiles of the Ly-alpha line. They show broad wing emission components on the red side, and sharp cut-off on the blue side of the Ly-alpha line. Such asymmetry is often seen in superwind galaxies known to date, and also consistent with a theoretical prediction of superwind activity. There are eight objects (8/18 ~ 40%) that have large EWs exceeding 200 A, and no clear signature of superwind activities. Such large EWs cannot be explained in terms of photo-ionization by a moderately old (>10^7 yrs) stellar population, even with a top-heavy IMF or an extremely low metallicity. These eight objects clearly show a positive correlation between the Ly-alpha luminosity and the velocity width. This suggests that these eight objects are good candidates for forming-galaxies in a gas-cooling phase.

We present the results of spectroscopy of Lyman Break Galaxies (LBGs) at z~5 in the J0053+1234 field with the Faint Object Camera and Spectrograph on the Subaru telescope. Among 5 bright candidates with z' < 25.0 mag, 2 objects are confirmed to be at z~5 from their Ly alpha emission and the continuum depression shortward of Ly alpha. The EWs of Ly alpha emission of the 2 LBGs are not so strong to be detected as Ly alpha emitters, and one of them shows strong low-ionized interstellar (LIS) metal absorption lines. Two faint objects with z' \geq 25.0 mag are also confirmed to be at z~5, and their spectra show strong Ly alpha emission in contrast to the bright ones. These results suggest a deficiency of strong Ly alpha emission in bright LBGs at z~5, which has been discussed in our previous paper. Combined with our previous spectra of LBGs at z~5 obtained around the Hubble Deep Field-North (HDF-N), we made a composite spectrum of UV luminous (M_1400 \leq -21.5 mag) LBGs at z~5. The resultant spectrum shows a weak Ly alpha emission and strong LIS absorptions which suggests that the bright LBGs at z~5 have chemically evolved at least to ~0.1 solar metallicity. For a part of our sample in the HDF-N region, we obtained near-to-mid infrared data, which constraint stellar masses of these objects. With the stellar mass and the metallicity estimated from LIS absorptions, the metallicities of the LBGs at z~5 tend to be lower than those of the galaxies with the same stellar mass at z \lesssim 2, although the uncertainty is very large.

SDSS J083942.11+380526.3 is an Iron Low-ionization Broad Absorption Line (FeLoBAL) quasar at z = 2.3, and Aoki et al. (2006) recently found the presence of an H alpha absorption line in the broad H alpha emission line. Motivated by an idea that this quasar may be a huge molecular gas reservoir in the early phase of quasar evolution, we made CO(J=3-2) observations of it using the Nobeyama Millimeter Array. No significant CO emission was detected; although an emission-like feature (2.5 sigma) was seen close (~ 2") to the quasar, we regard it as a noise. The obtained 3 sigma upper limit on the CO luminosity is L'_{{\rm CO}(J=3-2)} = 4.5 \times 10^{10} K km/s pc^2, which corresponds to M({\rm H}_2) = 3.6 times 10^10 M_{\odot} if we adopt the CO-to-H_2 conversion factor of 0.8 M_{\odot} (K km/s pc^2)^{-1}. This upper limit is comparable to L'_{{\rm CO}(J=3-2)} (and thus the molecular gas mass) detected in quasars and BAL quasars at z=1-3, and no sign of the presence of the huge amount of molecular gas in this FeLoBAL quasar was obtained.

(abridged) We report the UV luminosity function (LF) of Lyman break galaxies at z~5 derived from a deep and wide survey using the Subaru/Suprime-Cam. Target fields consist of two blank regions of the sky (the HDF-N and J0053+1234), and the total effective surveyed area is 1290 sqarcmin. Applications of carefully determined colour selection criteria in V-I and I-z' yield a detection of 853 candidates with z'AB<26.5 mag. The derived UVLF at z~5 shows no significant change in the number density of bright (L>=L*_z=3) LBGs from that at z~3, while there is a significant decline in the LF's faint end with increasing lookback time. This result means that the evolution of the number densities is differential with UV luminosity: the number density of UV luminous objects remains almost constant from z~5 to 3 while the number density of fainter objects gradually increases with cosmic time. This trend becomes apparent thanks to the small uncertainties in number densities both in the bright and faint parts of LFs at different epochs that are made possible by the deep and wide surveys. We discuss the origins of this differential evolution and suggest that our observational findings are consistent with the biased galaxy evolution scenario: a galaxy population hosted by massive dark haloes starts active star formation preferentially at early cosmic time, while less massive galaxies increase their number density later. We also calculated the UV luminosity density by integrating the UVLF and at z~5 found it to be 38.8% of that at z~3 for the luminosity range L>0.1 L^*_z=3. By combining our results with those from the literature, we find that the cosmic UV luminosity density marks its peak at z=2-3 and then slowly declines toward higher redshift.

We present the first measurement of clustering properties of low mass galaxies with a stellar mass down to M_*~10^9 Msun at 1<z<4 in 24.4 arcmin^2 of the GOODS-North region with a depth of K_{AB}~25, based on the near infrared observations performed with MOIRCS at the Subaru Telescope. The correlation amplitude strongly depends on the K-band flux, color, and stellar mass of the galaxies. We find that K-band luminous galaxies have a larger correlation length than K-band faint galaxies. For color selected samples at 2<z<4, distant red galaxies with J-K>1.3 show a large bias of b~7.2+-1.3 on scales up to \theta~100" or 3.1 comoving Mpc, while blue galaxies with 0.5<J-K<1.3, in which most Lyman break galaxies are populated, have a weak clustering signal on large scales, but a possible strong small scale excess at \theta<10". For massive galaxies with M_*>~10^{10} Msun, we estimate a correlation length and bias to be r_0~4.5 h^{-1} Mpc and b=1.9-3.5, which are much larger than those of low mass (M_*~10^9-10^{10} Msun) galaxies. The comparison of our measurements with analytic CDM models constrains the properties of hosting dark halos, and indicates that the low mass galaxies would be progenitors of galaxies with a typical luminosity of L<~L_* in the local Universe. The blue galaxies in low mass samples are more strongly clustered in more massive halos with higher occupation numbers than low mass red galaxies. This fact suggests an environment effect due to the halo mass on star formation activity at high-z.

We discovered non-stellar Balmer absorption lines in two many-narrow-trough FeLoBALs (mntBALs) by the near-infrared spectroscopy with Subaru/CISCO. Presence of the non-stellar Balmer absorption lines is known to date only in the Seyfert galaxy NGC 4151, thus our discovery is the first cases for quasars. Since all known active galactic nuclei with Balmer absorption lines share characteristics, it is suggested that there is a population of BAL quasars which have unique structures at their nuclei or unique evolutionary phase.

We study the properties of very faint, sub-L* Lyman break galaxies at z~2-5 - thus far a largely neglected but numerically and energetically very important population. We find that the LBG luminosity function undergoes luminosity-dependent evolution: the number of luminous galaxies remains constant while the number of faint ones grows with time. The total UV luminosity density increases with cosmic time from at least z~5 until reaching a peak or a plateau around z~2 - behaviour that is governed by the sub-L* galaxies in the LF's "faint tail". Using broadband SED fitting we find a nearly-linear relationship between SFR and galaxy stellar mass at z~2. A typical L* LBG at z~2 shows a stellar mass of ~10^10M_sun, remarkably similar to the bimodality mass at low redshift. This similarity suggests that the mechanisms responsible for the galaxy bimodality at low-z may have also been at play at z~2.

In this contribution we briefly describe our recent results on the properties of Lyman break galaxies at z~5 obtained from deep and wide blank field surveys using Subaru telescope, and through the comparison with samples at lower redshift ranges we discuss the evolution of star-forming galaxies in the early universe.

We use very deep near-infrared imaging data taken with Multi-Object InfraRed Camera and Spectrograph (MOIRCS) on the Subaru Telescope to investigate the number counts of Distant Red Galaxies (DRGs). We have observed a 4x7 arcmin^2 field in the Great Observatories Origins Deep Survey North (GOODS-N), and our data reach J=24.6 and K=23.2 (5sigma, Vega magnitude). The surface density of DRGs selected by J-K>2.3 is 2.35+-0.31 arcmin^-2 at K<22 and 3.54+-0.38 arcmin^-2 at K<23, respectively. These values are consistent with those in the GOODS-South and FIRES. Our deep and wide data suggest that the number counts of DRGs turn over at K~22, and the surface density of the faint DRGs with K>22 is smaller than that expected from the number counts at the brighter magnitude. The result indicates that while there are many bright galaxies at 2<z<4 with the relatively old stellar population and/or heavy dust extinction, the number of the faint galaxies with the similar red color is relatively small. Different behaviors of the number counts of the DRGs and bluer galaxies with 2<z_phot<4 at K>22 suggest that the mass-dependent color distribution, where most of low-mass galaxies are blue while more massive galaxies tend to have redder colors, had already been established at that epoch.

We describe deep radio imaging at 1.4-GHz of the 1.3 square degree Subaru/XMM-Newton Deep Field (SXDF), made with the Very Large Array in B and C configurations. We present a radio map of the entire field, and a catalogue of 505 sources covering 0.8 square degrees to a peak flux density limit of 100 microJy. Robust optical identifications are provided for 90% of the sources, and suggested IDs are presented for all but 14 (of which 7 are optically blank, and 7 are close to bright contaminating objects). We show that the optical properties of the radio sources do not change with flux density, suggesting that AGNs continue to contribute significantly at faint flux densities. We test this assertion by cross-correlating our radio catalogue with the X-ray source catalogue and conclude that radio-quiet AGNs become a significant population at flux densities below 300 microJy, and may dominate the population responsible for the flattening of the radio source counts if a significant fraction of them are Compton-thick.

We discovered an H alpha absorption in a broad H alpha emission line of an unusual broad absorption line quasar, SDSS J083942.11+380526.3 at z=2.318, by near-infrared spectroscopy with the Cooled Infrared Spectrograph and Camera for OHS (CISCO) on the Subaru telescope. The Presence of non-stellar H alpha absorption is known only in the Seyfert galaxy NGC 4151 to date, thus our discovery is the first case for quasars. The H alpha absorption line is blueshifted by 520 km/s relative to the H alpha emission line, and its redshift almost coincides with those of UV low-ionization metal absorption lines. The width of the H alpha absorption (~ 340 km/s) is similar to those of the UV low-ionization absorption lines. These facts suggest that the H alpha and the low-ionization metal absorption lines are produced by the same low-ionization gas which has a substantial amount of neutral gas. The column density of the neutral hydrogen is estimated to be ~ 10^18 cm^-2 by assuming a gas temperature of 10,000 K from the analysis of the curve of growth. The continuum spectrum is reproduced by a reddened (E(B-V) ~ 0.15 mag for the SMC-like reddening law) composite quasar spectrum. Furthermore, the UV spectrum of SDSS J083942.11+380526.3 shows a remarkable similarity to that of NGC 4151 in its low state, suggesting the physical condition of the absorber in SDSS J083942.11+380526.3 is similar to that of NGC 4151 in the low state. As proposed for NGC 4151, SDSS J083942.11+380526.3 may be also seen through the close direction of the surface of the obscuring torus.

We report a deficiency of luminous Lyman break galaxies (LBGs) with a large rest-frame equivalent width (EW_rest) of Lyman-alpha emission at z~5-6. Combining our spectroscopic sample of LBGs at z~5 and those from the literature, we found that luminous LBGs at z~5-6 generally show weak Lyman-alpha emissions, while faint LBGs show a wide range of Lyman-alpha EW_rest and tend to have strong (EW_rest >20A) Lyman-alpha emissions; i.e., there is a deficiency of strong Lyman-alpha emission in luminous LBGs. There seems to be a threshold UV luminosity for the deficiency; it is M_1400 = -21.5 \~ -21.0 mag, which is close to or somewhat brighter than the M* of the UV luminosity function at z~5 and 6. Since the large EW_rest of Lyman-alpha emission can be seen among the faint LBGs, the fraction of Lyman-alpha emitters in LBGs may change rather abruptly with the UV luminosity. If the weakness of Lyman-alpha emission is due to dust absorption, the deficiency suggests that luminous LBGs at z=5-6 tend to be in dusty and more chemically evolved environments and started star formation earlier than faint ones, though other causes cannot be ruled out.

Spatially extended Ly-alpha sources which are faint and compact in coninuum are candidates for extremely young galaxies (age of ~<10e+7 yrs) at high redshifts. We carried out a systematic survey for extended Ly-alpha sources, using deep intermediate-band imaging data taken with the Suprime-Cam on the Subaru Telescope. Our survey covers a field of view of 33' * 25' and a redshift range of 3.24<z<4.95, down to a Ly-alpha flux of ~1e-17 ergs/cm^2/s. We identified 41 extended Ly-alpha sources. The redshift distribution of these sources shows that this kind of objects are common in the early universe over the surveyed redshift range. The identified objects have typical sizes of \~10-15 kpc and luminosities of ~10^42 ergs/s. Follow-up spectroscopy made for 7 of the 41 objects showed that our sample suffers from little contamination. All 7 objects have large equivalent widths of Ly-alpha emission line, all but one exceeding 240A in the rest frame. The large equivalent widths suggest that their extended Ly-alpha emissions are unlikely to be due to normal starbursts, but possibly originated from very young galaxies. All 41 objects in our sample have much smaller Ly-alpha luminosities than the two Ly-alpha Blobs (LABs) found by Steidel et al. (2000) in spite of our much larger survey volume. This suggests that large and luminous extended Ly-alpha objects like the two LABs are very rare and are clustered in overdense regions.

We study properties of the host galaxies of 15 hard X-ray selected type-2 active galactic nuclei (AGNs) at intermediate redshifts (0.05$<z<$0.6) detected in $ASCA$ surveys. The absorption corrected hard X-ray luminosities $L_{\rm 2-10 keV}$ range from 10$^{42}$ erg s$^{-1}$ to $10^{45}$ erg s$^{-1}$. We took the $R$-band image of these AGNs with the University of Hawaii 2.2 m telescope. Thanks to the intrinsic obscuration of nuclear light, we can decompose the galaxies with a spheroid component and a disk component. The resulting spheroid luminosities correlate with $L_{\rm 2-10 keV}$; higher (lower) X-ray luminosity AGNs tend to reside in luminous (less luminous) spheroids. It is also found that the hosts of luminous AGNs show a large spheroid-to-disk luminosity ratio ($\sim$1), while those of less luminous AGNs spread between 0 and 1. The correlation between $L_{\rm 2-10keV}$ and spheroid luminosity indicates that the relation between mass of a supermassive black hole (SMBH) and spheroid luminosity (BS-relation) at the intermediate redshifts. The BS-relation agrees with that in the local universe if the Eddington ratio of 0.24 is adopted, which is a mean value determined from our $ASCA$ type-1 AGN sample at similar redshifts through the broad-line width and continuum luminosity. The present study demonstrates the effectiveness of using type-2 AGNs at high redshifts to study their host properties.

We present the XMM-Newton and the optical-VLT spectra along with the optical and the near-infrared photometric data of one of the brightest X-ray (F(2-10 keV)~1e-13 erg/s cm^2) extremely red objects (R-K>=5) discovered so far. The source, XBSJ0216-0435, belongs to the XMM-Newton Bright Serendipitous Survey and it has extreme X-ray-to-optical (~220) and X-ray-to-near-infrared (~60) flux ratios. Thanks to its brightness, the X-ray statistics are good enough for an accurate spectral analysis by which the presence of an X-ray obscured (NH>1e22 cm^-2) QSO (L(2-10 keV)=4e45 erg/s) is determined. A statistically significant (~99%) excess around 2 keV in the observed-frame suggests the presence of an emission line. By assuming that this feature corresponds to the iron Kalpha line at 6.4 keV, a first estimate of the redshift of the source is derived (z_x~2). The presence of a high redshift QSO2 has been finally confirmed through dedicated VLT optical spectroscopic observations (z_o=1.985+/-0.002). This result yields to an optical validation of a new X-ray Line Emitting Object (XLEO) for which the redshift has been firstly derived from the X-ray data. XBSJ0216-0435 can be considered one of the few examples of X-ray obscured QSO2 at high redshift for which a detailed X-ray and optical spectral analysis has been possible. The spectral energy distribution from radio to X-rays is also presented. Finally from the near-infrared data the luminosity and the stellar mass of the host galaxy has been estimated finding a new example of the coexistence at high-z between massive galaxies and powerful QSOs.

We report angular correlation function (ACF) of Lyman Break Galaxies (LBGs) with unprecedented statistical quality on the basis of 16,920 LBGs at z=4 detected in the 1 deg^2 sky of the Subaru/XMM-Newton Deep Field. The ACF significantly departs from a power law, and shows an excess on small scale. Particularly, the ACF of LBGs with i'<27.5 have a clear break between the small and large-scale regimes at the angular separation of ~7'' whose projected length corresponds to the virial radius of dark halos with a mass of 10^11-12 Mo, indicating multiple LBGs residing in a single dark halo. Both on small (2''<theta<3'') and large (40''<theta<400'') scales, clustering amplitudes monotonically increase with luminosity for the magnitude range of i'=24.5-27.5, and the small-scale clustering shows a stronger luminosity dependence than the large-scale clustering. The small-scale bias reaches b~10-50, and the outskirts of small-scale excess extend to a larger angular separation for brighter LBGs. The ACF and number density of LBGs can be explained by the cold dark matter model.

We report results of deep optical spectroscopy with Subaru/FOCAS for Lyman Break Galaxy (LBG) candidates at z~5. So far, we made spectroscopic observations for 24 LBG candidates among ~200 bright (z'<25.0) LBG sample and confirmed 9 objects to be LBGs at z~5. Intriguingly, these bright LBGs show no or a weak Ly-alpha emission and relatively strong low ionization interstellar metal absorption lines. We also identified 2 faint (z'>25.0) objects to be at z~5 with their strong (EW_rest>20AA) Ly-alpha emission. Combining our results with other spectroscopic observations of galaxies at the similar redshift range, we found a clear luminosity dependence of EW_rest of Ly-alpha emission, i.e., the lack of strong Ly-alpha emission in bright LBGs. If the absence of Ly-alpha emission is due to dust absorption, these results suggest that bright LBGs at z~5 are in dusty and more chemically evolved environment than faint ones. This interpretation implies that bright LBGs started star formation earlier than faint ones, suggesting biased star formation.

We report the updated UV luminosity function (LF) of Lyman break galaxies (LBGs) at z~5. Combining with lower redshift data, we found that the evolution of UV LF is differential depending on the UV luminosity; i.e., the number density of luminous LBGs does not change significantly from z=5 to 3, while the number density of faint (L<L*) galaxies increases by a factor of five. We also found the signs of differential evolution within other observational results for LBGs at z~5. UV-luminous objects show strong clustering, and they have weak or no Lyman-alpha emissions while showing prominent metal absorption lines. There is also a clear deficiency of objects with blue UV-optical colors for those bright in rest-frame optical bands (from the Spitzer/IRAC data obtained by the GOODS project). By the SED fitting using rest-frame optical to near-infrared photometry data it is suggested that most luminous objects at z~5 have already accumulated stellar mass more than 10^10 Msun. We discuss that this differential evolution might be understood as a consequence of the biased galaxy evolution.

We obtained the number counts and the rest-frame B-band luminosity function of the color-selected old passively-evolving galaxies (OPEGs) at z=1 with very high statistical accuracy using a large and homogeneous sample of about 4000 such objects with z' <25 detected in the area of 1.03 deg^2 in the Subaru/XMM-Newton Deep Survey (SXDS) field. Our selection criteria are defined on the i'-z' and R-z' color-magnitude plane so that OPEGs at z=0.9-1.1 with formation redshift z_f=2-10 are properly sampled. The limiting magnitude corresponds to the luminosity of galaxies with M_*+3 at z=0. We made a pilot redshift observations for 99 OPEG candidates with 19 < z' < 22 and found that at least 78% (73/93) of the entire sample, or 95% (73/77) of these whose redshifts were obtained are indeed lie between z=0.87 and 1.12 and the most of their spectra show the continuum break and strong Ca H and K lines, indicating that these objects are indeed dominated by the old stellar populations. We then compare our results with the luminosity functions of the color- or the morphologically-selected early type galaxies at z=0 taking the evolutionary factor into account and found that the number density of old passive galaxies with sim M_* magnitude at z~1 averaged over the SXDS area is 40-60% of the equivalently red galaxies and 60-85% of the morphologically-selected E/S0 galaxies at z=0 depending on their luminosity evolution. It is revealed that more than half, but not all, of the present-day early-type galaxies had already been formed into quiescent passive galaxies at z=1.

The star-formation rates and the stellar masses of the host galaxies of AGNs at high-redshifts are keys to understanding the evolution of the relation between the mass of the spheroidal component of a galaxy and the mass of its central black hole. We investigate the host galaxies of 31 AGNs with spectroscopic redshifts between 2 and 4 found in the deep Chandra surveys of the GOODS fields with the ACS images. The sample can be divided into 17 ``extended'' AGNs and 14 ``compact'' AGNs based on the concentration parameter defined as the difference between the aperture magnitudes with 0.07" and 0.25" diameter. We derive upper limits of the UV luminosities of the host galaxies of the ``compact'' AGN sample, and upper and lower limits of the UV luminosities of the host galaxies of the ``extended'' AGN sample.These limits are consistent with the knee of the luminosity function of the Lyman Break Galaxies (LBGs) at z=3, suggesting moderate star-formation rates, less than 40 Msolar/yr, in the host galaxies of the AGNs at 2<z_sp<4 without correcting the dust extinction. By combining the limits of the UV luminosities with the observed K-band magnitudes for the ``extended'' AGNs, we derive the upper and lower limits of the stellar masses of their host galaxies. The derived upper limits on the stellar masses range from 10^10 Msolar to 10^12 Msolar. The upper limits imply that the Mbulge-MBH relation of the high-redshift AGNs is different from that of the galaxies in the nearby universe or the average Eddington ratio of the high-redshift AGNs is higher than that of low-redshift AGNs with lower-luminosity.

We report the serendipitous discovery of several unresolved X-ray sources lying in the prominent spiral arms of the galaxy KUG 0214-057 in XMM-Newton observations. The location of these X-ray sources strongly suggests that at least three, and possibly four, of these may be physically related to the galaxy. The luminosity of each of these sources at the distance of KUG 0214-057 is > 5 x 10^39 erg/s (0.3-10 keV), making each a strong candidate ultraluminous X-ray source (ULX). Using the ULXs as a metric implies that this relatively low-mass galaxy may be experiencing rather intense starburst activity. The serendipitous discovery of these ULXs suggests that such objects are not a negligible component of the overall extragalactic X-ray source population.

We report the discovery of primeval large-scale structures (LSSs) including two proto-clusters in a forming phase at z=5.7. We carried out extensive deep narrow-band imaging in the 1 deg^2 sky of the Subaru/XMM-Newton Deep Field, and obtained a cosmic map of 515 Ly alpha emitters (LAEs) in a volume with a transverse dimension of 180 Mpc x 180 Mpc and a depth of ~40 Mpc in comoving units. This cosmic map shows filamentary LSSs, including clusters and surrounding 10-40 Mpc scale voids, similar to the present-day LSSs. Our spectroscopic follow-up observations identify overdense regions in which two dense clumps of LAEs with a sphere of 1-Mpc diameter in physical units are included. These clumps show about 130 times higher star-formation rate density, mainly due to a large overdensity, ~80, of LAEs. These clumps would be clusters in a formation phase involving a burst of galaxy formation.

The Subaru Deep Field (SDF) project is a program of Subaru Observatory to carry out a deep galaxy survey over a blank field as large as 34'x27'. The program consists of very deep multi-band optical imaging, near infrared imaging for smaller portions of the field and follow-up optical spectroscopy. Major scientific goals of the project are to construct large samples of Lyman-break galaxies at z~4-5 and Lyman alpha emitters at z~5.7 and 6.6, and to make detailed studies these very high-redshift galaxy populations. In this paper, we describe the optical imaging observations and data reduction, presenting mosaicked images and object catalogs in seven bandpasses.The optical imaging was made through five broad-band filters, B, V, R, i', z', and two narrow-band filters, NB816 (lambda_c=8150A) and NB921 (lambda_c=9196A) with almost 10 hours long integrations for each band. The limiting magnitudes measured at 3-sigma on a 2" aperture are B=28.45, V=27.74, R=27.80, i'=27.43, z'=26.62, NB816=26.63, and NB921=26.54 in the AB system. The object catalog constructed for each of the seven bands contains more than 10^5 objects. The galaxy number counts corrected for detection incompleteness and star count contribution are found to be consistent with previous results in the literature. The mosaicked images and catalogs of all the bands have been made open to the public on Oct. 1, 2004 on the SDF project website at http://soaps.naoj.org/sdf/.

We report on the discovery of a star-forming galaxy at z=6.33 in the Subaru Deep Field. This object is selected as a candidate of an i'-dropout, high-redshift galaxy around z=6 because of its red i'-z' color in our deep optical imaging survey in the Subaru Deep Field. Our follow up optical spectroscopy reveals that this object is a strong Ly-alpha emitter with only very faint ultraviolet continuum. The rest-frame equivalent width of the detected Ly-alpha emission is as much as 130 A. Thus the light detected in our z' image is largely attributed to the Ly-alpha emission, i.e., ~40% of the z'-band flux is the strong Ly-alpha emission, giving a very red i'-z' color. This is consistent with the photometric property of this object because the narrow-band data obtained with the NB921 filter shows a significant depression, z'-NB921 = -0.54 mag. By using the photometric data, we show that some other objects among the 48 i'-dropout high-redshift galaxy candidates found in the Subaru Deep Field also show a significant NB921 depression. We briefly discuss the nature of these NB921-depressed objects.

We report initial results for spectroscopic observations of candidates of Lyman Break Galaxies (LBGs) at $z\sim5$ in a region centered on the Hubble Deep Field-North by using the Faint Object Camera and Spectrograph attached to the Subaru Telescope. Eight objects with $I_C\leq25.0$ mag, including one AGN, are confirmed to be at $4.5<z<5.2$. The rest-frame UV spectra of seven LBGs commonly show no or weak Lyalpha emission line (rest-frame equivalent width of 0-10\AA) and relatively strong low-ionization interstellar metal absorption lines of SiII $\lambda$1260, OI+SiII $\lambda$1303, and CII $\lambda$1334 (mean rest-frame equivalent widths of them are $-1.2 \sim -5.1 $\AA). These properties are significantly different from those of the mean rest-frame UV spectrum of LBGs at $z\sim3$, but are quite similar to those of subgroups of LBGs at $z\sim3$ with no or weak Lyalpha emission. The weakness of Lyalpha emission and strong low-ionization interstellar metal absorption lines may indicate that these LBGs at $z\sim5$ are chemically evolved to some degree and have a dusty environment. Since the fraction of such LBGs at $z\sim5$ in our sample is larger than that at $z\sim3$, we may witness some sign of evolution of LBGs from $z\sim5$ to $z\sim3$, though the present sample size is very small. It is also possible, however, that the brighter LBGs tend to show no or weak Lyalpha emission, because our spectroscopic sample is bright (brighter than $L^{\ast}$) among LBGs at $z\sim5$. More observations are required to establish spectroscopic nature of LBGs at $z\sim5$.

We present results of near-infrared photometry (J, H, K_S) for a sample of active galactic nuclei (AGNs) obtained from hard X-ray surveys with ASCA. The sample covers the AGNs at z=0.1-1 with L(2-10keV)=10^42-10^46 erg/s with very high completeness. The fraction of red (J-K_S>2 mag) AGNs in our sample is 2(+-1)%, which is comparable to that for optically- or UV-selected quasi-stellar objects (QSOs, i.e. luminous AGNs). The number of red AGNs found in our sample is also consistent with that expected from the surface density of red AGNs found in 2MASS by Cutri et al. (2001). We find that the anomalously-small dust-to-gas ratios in circumnuclear gas, which is seen in some AGNs with Seyfert-class luminosity, also occur in the QSOs (AGNs with luminosity of L(2-10keV) > 10^44.5 erg/s). For all the QSOs with an X-ray absorption of N_H > 10^22 /cm2 in our sample, the values of A_V/N_H are smaller than the Galactic value by a factor of 5 to 100. Since a fraction of this population among the QSOs in our sample is about 30%, such fraction of optical/UV-selected type 1 QSOs known to date may show type 2 nature in X-ray.

We investigate the cosmological evolution of the hard X-ray luminosity function (HXLF) of Active Galactic Nuclei (AGN) in the 2-10 keV luminosity range of 10^{41.5} - 10^{46.5} erg s^-1 as a function of redshift up to 3. From a combination of surveys conducted at photon energies above 2 keV with HEAO1, ASCA, and Chandra, we construct a highly complete (>96%) sample consisting of 247 AGNs over the wide flux range of 10^{-10} - 3.8*10^{-15} erg cm^-2 s^-1 (2-10 keV). For our purpose, we develop an extensive method of calculating the intrinsic (before-absorption) HXLF and the absorption (N_H) function. This utilizes the maximum likelihood method fully correcting for observational biases with consideration of the X-ray spectrum of each source. We find that (i) the fraction of X-ray absorbed AGNs decreases with the intrinsic luminosity and (ii) the evolution of the HXLF of all AGNs (including both type-I and type-II AGNs) is best described with a luminosity dependent density evolution (LDDE) where the cutoff redshift increases with the luminosity. Our results directly constrain the evolution of AGNs that produce a major part of the hard X-ray background, thus solving its origin quantitatively. A combination of the HXLF and the NH function enables us to construct a purely "observation based" population synthesis model. We present basic consequences of this model, and discuss the contribution of Compton-thick AGNs to the rest of the hard X-ray background.

Results of optical identification of the ASCA Lynx deep survey are presented. Six X-ray sources are detected in the 2-7 keV band using the SIS in a 20'x20' field of view with fluxes larger than ~4x10^{-14} erg s-1 cm-2 in the band. Follow-up optical spectroscopic observations were made, and five out of six sources are identified with AGNs/QSOs at redshifts of 0.5-1.3. We also identify two more additional X-ray sources detected in a soft X-ray band with AGNs/QSOs. It is found that three QSOs identified are located at z~1.3. Two rich clusters and several groups of galaxies are also placed at the same redshift in the surveyed field, and projected separations between the QSOs and the clusters are 3-8 Mpc at the redshift.

We present the results of optical spectroscopic identifications of a bright subsample of 2-10keV hard X-ray selected sources from the ASCA Medium Sensitivity Survey in the northern sky. The flux limit of the subsample is 3*10^-13 erg s^-1 cm^-2 in the 2-10keV band. All but one of the 87 hard X-ray selected sources are optically identified, with AGNs, 7 clusters of galaxies, and 1 galactic star. It is the largest complete sample of hard X-ray selected AGNs at the bright flux limit. Amounts of absorption to their nuclei are estimated to be hydrogen column densities (N_H) of up to ~3*10^23 cm^-2 from their X-ray spectra. Optical properties of X-ray absorbed AGNs with N_H > 1*10^22 cm^-2 indicate the effects of dust absorption: at redshifts, z<0.6, AGNs without broad H\beta emission lines have significantly larger N_H value than AGNs with broad H\beta emission lines. At z>0.6, the X-ray absorbed AGNs have a large hard X-ray to optical flux ratio (log f2-10keV/fR > +1). However, three X-ray absorbed z>0.6 AGNs show strong broad lines. In combination with hard X-ray selected AGN samples from the ASCA Large Sky Survey, the ASCA Deep Survey in the Lockman Hole and Chandra Deep Field North, the luminosity distributions of absorbed and less-absorbed AGNs are compared.

We have performed a detailed analysis of 3 optically normal galaxies extracted from the XMM Bright Serendipitous Source Sample. Thanks to the good statistics of the XMM-Newton data, we have unveiled the presence of an AGN in all of them. In particular, we detect both X-ray obscured (N_H>10^{22} cm^{-2}) and unobscured (N_H<10^{22} cm^{-2}) AGN with intrinsic 2--10 keV luminosities in the range between 10^{42} -- 10^{43} erg s^{-1}. We find that the X-ray and optical properties of the sources discussed here could be explained assuming a standard AGN hosted by galaxies with magnitudes M_R<M^*, taking properly into account the absorption associated with the AGN, the optical faintness of the nuclear emission with respect to the host galaxy, and the inadequate set--up and atmospheric conditions during the optical spectroscopic observations. Our new spectroscopic observations have revealed the expected AGN features also in the optical band. These results clearly show that optical spectroscopy sometimes can be inefficient in revealing the presence of an AGN, which instead is clearly found from an X-ray spectroscopic investigation. This remarks the importance of being careful in proposing the identification of X-ray sources (especially at faint fluxes) when only low quality optical spectra are in hand. This is particularly important for faint surveys (such as those with XMM-Newton and Chandra), in which optically dull but X-ray active objects are being found in sizeable numbers.

(abridged) We present results of a search for Lyman break galaxies (LBGs) at z ~ 5 in a 618 square-arcmin field including the HDF-N taken by Subaru Prime Focus Camera. Utilizing the published redshift data of the HDF-N and its flanking fields, the color selection criteria are chosen so that LBGs are picked out most efficiently and least contaminated by foreground objects. The numbers of LBG candidates detected are 310 in 23.0 < I_c < 25.5. The rest-frame UV luminosity function(LF) of LBGs at z ~ 5 is derived statistically. The fraction of contamination is estimated to be ~50% in the faintest magnitude range. The completeness of the survey is ~80% at the bright part of the sample, and ~20% in the faintest magnitude range (25.0 < I_c <= 25.5). The LF of LBG candidates at z ~ 5 does not show a significant difference from those at z ~ 3 and 4, though there might be a slight decrease in the fainter part. The UV luminosity density within the observational limit is 0.56 - 0.69 times smaller than that obtained for LBGs at z ~ 3, depending on the adopted cosmology and the integration range of the LF. The similarity of the LFs at redshifts 5 to 3 implies that most of LBGs at z ~ 5 should have faded out at z ~ 3 and LBGs at z \~ 5 are different galaxies from those seen at z ~ 3, if we take face values for ages of the LBGs at z ~ 3 obtained by the SED fitting in which a continuous star formation in an individual galaxy is assumed. However, if the star formation in LBGs is sporadic, the similarity of the LF at z ~ 3 and 5 would be explained. Such sporadic star formation has been suggested by hydrodynamical simulations and semi-analytic models with collisional starbursts, and the trend of the cosmic star formation history predicted by these studies resembles to that estimated from the UV luminosity density within the observational limit.

We present the first results of the Chandra and optical follow-up observations of hard X-ray sources detected in the ASCA Medium Sensitivity Survey (AMSS). Optical identifications are made for five objects. Three of them show either weak or absent optical narrow emission lines and are at low redshift <z>~0.06. One of them is a broad line object at z=0.910 and one is a z=0.460 object with only narrow lines. All the narrow line objects show strong evidence for absorption in their X-ray spectra. Their line ratios are consistent with a Seyfert II/LINER identification as are the line widths. The three low redshift objects have the colors of normal galaxies and apparently the light is dominated by stars. This could be due to the extinction of the underlying nuclear continuum by the same matter that absorbs X-rays and/or due to the dilution of the central source by starlight. These results suggest that X-ray sources that appear as ``normal'' galaxies in optical and near-IR bands significantly contribute to the hard X-ray background. This population of objects has a high space density and probably dominates the entire population of active galaxies.

To understand luminous AGNs in the z<1 universe, the ASCA AGN samples are the best at present. Combining the identified sample of AGNs from ASCA Large Sky Survey and Medium Sensitivity Survey, the sample of hard X-ray selected AGNs have been expanded up to 108 AGNs above the flux limit of 10^{-13} erg s^{-1} cm^{-2} in the 2-10 keV hard X-ray band. We discuss the fraction of absorbed AGNs in the hard X-ray selected AGN sample, and nature of absorbed luminous AGNs.

Deep hard X-ray and near infrared observations of a type-2 radio-quiet QSO candidate at z=0.9, AXJ08494+4454, are reported. The 0.5-10keV Chandra X-ray Observatory spectrum of AXJ08494+4454 is hard, and is explained well with a power-law continuum absorbed by a hydrogen column density of 2x10^23 cm^-2. The 2-10keV luminosity of the object reaches hard X-ray luminosity of QSOs. The large X-ray absorption and the large intrinsic luminosity support the original identification of AXJ08494+4454 as a type-2 radio-quiet QSO. Nevertheless, the deep Subaru/IRCS J-band spectroscopic observation suggests the presence of a strong broad H-alpha emission line from AXJ08494+4454. If real, the broad H-alpha emission line has a velocity width of 9400+-1000 km/s, which corresponds to a typical broad-Balmer line velocity width of a luminous QSO. The existence of the strong broad H-alpha line means that the object is not a type-2 QSO, but a luminous cousin of a Seyfert 1.9 galaxy in the source-frame optical spectrum. The Balmer decrement of broad lines, and the broad H-alpha emission to the hard X-ray luminosity ratio suggest that the nucleus is affected by dust extinction with A_V of 1-3 mag in the optical wavelength. The estimated amount of dust extinction is much smaller than that expected from the X-ray column density.

In this paper, results of optical identification of ASCA surveys are summarized. To understand luminous AGNs in the z<1 universe, the ASCA AGN sample is still better than samples of AGNs from deep Chandra or XMM-Newton surveys. Combining the identified sample of AGNs from ASCA Large Sky Survey and Medium Sensitivity Survey, the sample of hard X-ray selected AGNs have been expanded up to about 80 AGNs above the flux limit of 10^{-13} erg s^{-1} cm^{-2} in the 2-10keV band. One fifth of the ASCA AGNs are absorbed narrow-line AGNs. The luminosity distribution of the absorbed AGNs is limited below 10^{45} erg s^{-1}, though there are 14 non-absorbed broad-line AGNs are detected above the luminosity. The result suggests lack of absorbed luminous AGNs (deficiency of type 2 QSO).

The results of K-band-imaging observations of a candidate of an absorbed QSO at z=0.653, AX J131831+3341, are presented. The B-K color of the object is 4.85 mag, which is much redder than optically-selected QSOs. The K-band image shows nuclear and extended components, the same as in the optical V-, R-, and I-band images. The nuclear component (I-K = 4.29 mag) is much redder than the power-law models with energy indices of 0 to -1.0, which well reproduce the V-R and R-I optical colors of the nuclear component. A heavily absorbed (A_V = 3 mag) nucleus may emerge in the K-band, while optical light may originate from scattered nuclear light. The I-K color of the extended component is 2.2 mag, which is consistent with the post-starburst nature of the host galaxy, which is also suggested from the V-R and R-I colors of the extended component.

The results of deep-imaging observations of a candidate of an absorbed QSO at z=0.653, AX J131831+3341, are presented. AX J131831+3341 was found during the course of optical follow-up observations of the ASCA Large Sky Survey, and has an X-ray luminosity of 10^{45} erg s^{-1} (2-10keV), which corresponds to those of QSOs. Its optical spectrum shows no significant broad H_beta emission line, suggesting that the object is an absorbed QSO. Deep R and V band images reveal the presence of a point-like nucleus and an asymmetric extended component. The nuclear component has a blue color, and the optical magnitude is much fainter than that expected from the observed X-ray flux for typical type-1 AGNs. These photometric properties and the presence of broad MgII 2800A emission can be explained simultaneously if the observed nuclear light is dominated by scattered nuclear light, though there is a possibility that the nuclear component is a slightly absorbed nucleus if its intrinsic X-ray to optical flux ratio is the largest among X-ray selected AGNs. The size of the extended component, which is thought to be the host galaxy of the QSO, is larger than those of normal disk galaxies at z=0-0.75, and the absolute magnitude is similar to those of the brightest host galaxies of QSOs at redshifts smaller than 0.5. The V-R and R-I colors of the component are consistent with a 1 Gyr-old stellar population model without absorption.

We present results of optical identification of the X-ray sources detected in the ASCA Large Sky Survey. Optical spectroscopic observations were done for 34 X-ray sources which were detected with the SIS in the 2-7 keV band above 3.5 sigma. The sources are identified with 30 AGNs, 2 clusters of galaxies, and 1 galactic star. Only 1 source is still unidentified. The flux limit of the sample corresponds to 1 x 10^{-13} erg s^{-1} cm^{-2} in the 2-10 keV band. Based on the sample, the paper discusses optical and X-ray spectral properties of the AGNs, contribution of the sources to the Cosmic X-ray Background, and redshift and luminosity distributions of the AGNs. An interesting result is that the redshift distribution of the AGNs suggests a deficiency of high-redshift (0.5 < z < 2) and luminous (Lx(2-10 keV) > 10^{44} erg s^{-1}) absorbed narrow-line AGNs (so called type 2 QSOs).

We carried out the first wide-area unbiased survey with the ASCA satellite in the 0.7-10 keV band around a north Galactic-pole region covering a continuous area of 7 square degrees (Large Sky Survey; LSS). To make the best use of ASCA capability, we developed a new source-detection method where the complicated detector responses are fully taken into account. Applying this method to the entire LSS data independently in the total (0.7-7 keV), hard (2-10 keV), and soft (0.7-2 keV) band, we detected 107 sources altogether with sensitivity limits of 6 x 10E-14 (0.7-7 keV), 1 x 10E-13 (2-10 keV), and 2 x 10E-14 erg sE-1 cmE-2 (0.7-2 keV), respectively. A complete list of the detected sources is presented. Based on detailed studies by Monte Carlo simulations, we evaluated effects of the source confusion and accurately derived Log N - Log S relation in each survey band. The Log N - Log S relation in the hard band is located on the extrapolation from the GINGA and HEAO1 results with the Euclidean slope of -3/2, while that in the soft band is consistent with the results by ROSAT. At these flux limits, 30 (+/- 3) percent of the CXB in the 0.7-7 keV band and 23 (+/- 3) percent in the 2-10 keV band have been resolved into discrete sources. The average spectrum of faint sources detected in the total band shows a photon index of 1.63 +/- 0.07 in the 0.7-10 keV range, consistent with the comparison of source counts between the hard and the soft energy band. Those detected in the hard band show a photon index of 1.49 +/- 0.10 in the 2-10 keV range. These spectral properties suggest that contribution of sources with hard energy spectra become significant at a flux of 10E-13 erg sE-1 cmE-2 (2-10 keV). The most plausible candidates are type-II AGNs, as indicated by on-going optical identifications.

We present results of optical identification of the ASCA Large Sky Survey. X-ray sources which have hard X-ray spectra were identified with type-2 AGN at redshifts smaller than 0.5. It is supported that the absorbed X-ray spectrum of type-2 AGN makes the Cosmic X-ray Background harder in the hard X-ray band than type-1 AGN, which is main contributer in the soft X-ray band. Absence of type-2 AGN at redshift larger than 1 in the identified sample, which contrasts to the existence of 6 broad-line QSOs, suggests a deficiency of so-called ``type-2 QSO'' at high redshift.

We present results of ASCA deep exposure observations of the hardest X-ray source discovered in the ASCA Large Sky Survey (LSS) project, designated as AX J131501+3141. We extract its accurate X-ray spectrum, taking account of the contamination from a nearby soft source (AX J131502+3142), separated only by 1'. AX J131501+3141 exhibits a large absorption of NH = (6 +4 -2)x 10^22 H/cm^2 with a photon index \Gamma = 1.5 +0.7 -0.6. The 2--10 keV flux was about 5 x 10^-13 erg/s/cm^2 and was time variable by a factor of 30% in 0.5 year. From the highly absorbed X-ray spectrum and the time variability, as well as the results of the optical follow-up observations (Akiyama et al. 1998, astro-ph/9801173), we conclude that AX J131501+3141 is a type 2 Seyfert galaxy. Discovery of such a low flux and highly absorbed X-ray source could have a significant impact on the origin of the cosmic X-ray background.

We report the optical identification of the hardest X-ray source (AX J131501+3141) detected in an unbiased wide-area survey in the 0.5--10 keV band, the ASCA Large Sky Survey. The X-ray spectrum of the source is very hard and is well reproduced by a power-law component (Gamma = 1.5^+0.7_-0.6) with N_H = 6^+4_-2 *10^22 cm^-2 (Sakano et al. 1998). We have found a galaxy with R=15.62 mag near the center of the error circle for the X-ray source. The optical spectrum of the galaxy shows only narrow emission lines whose ratios correspond to those of a type 2 Seyfert galaxy at z = 0.072, implying an absorption-corrected X-ray luminosity of 2*10^43 erg sec^-1 (2--10 keV) and M_B = -20.93 mag. A radio point source is also associated with the center of the galaxy. We thus identify the X-ray source with this galaxy as an obscured AGN. The hidden nature of the nucleus of the galaxy in the optical band is consistent with the X-ray spectrum. These results support the idea that the obscured AGNs/QSOs contribute significantly to the cosmic X-ray background in the hard band at the faint flux level.