We have initiated a new survey for local extremely metal-poor galaxies (EMPGs) with Subaru/Hyper Suprime-Cam (HSC) large-area (~500 deg^2) optical images reaching a 5 sigma limit of ~26 magnitude, about 100 times deeper than the Sloan Digital Sky Survey (SDSS). To select Z/Z_sun<0.1 EMPGs from ~40 million sources detected in the Subaru images, we first develop a machine-learning (ML) classifier based on a deep neural network algorithm with a training data set consisting of optical photometry of galaxy, star, and QSO models. We test our ML classifier with SDSS objects having spectroscopic metallicity measurements, and confirm that our ML classifier accomplishes 86%-completeness and 46%-purity EMPG classifications with photometric data. Applying our ML classifier to the photometric data of the Subaru sources as well as faint SDSS objects with no spectroscopic data, we obtain 27 and 86 EMPG candidates from the Subaru and SDSS photometric data, respectively. We conduct optical follow-up spectroscopy for 10 out of our EMPG candidates with Magellan/LDSS-3+MagE, Keck/DEIMOS, and Subaru/FOCAS, and find that the 10 EMPG candidates are star-forming galaxies at z=0.007-0.03 with large H_beta equivalent widths of 104-265 A, stellar masses of log(M*/M_sun)=5.0-7.1, and high specific star-formation rates of ~300 Gyr^{-1}, which are similar to those of early galaxies at z>6 reported recently. We spectroscopically confirm that 3 out of 10 candidates are truly EMPGs with Z/Z_sun<0.1, one of which is HSC J1631+4426, the most metal-poor galaxy with Z/Z_sun=0.016 reported ever.

We present element abundance ratios and ionizing radiation of local young low-mass (~$10^{6}$ M_sun) extremely metal poor galaxies (EMPGs) with a 2% solar oxygen abundance (O/H)_sun and a high specific star-formation rate (sSFR~300 Gyr$^{-1}$), and other (extremely) metal poor galaxies, which are compiled from Extremely Metal-Poor Representatives Explored by the Subaru Survey (EMPRESS) and the literature. Weak emission lines such as [FeIII]4658 and HeII4686 are detected in very deep optical spectra of the EMPGs taken with 8m-class telescopes including Keck and Subaru (Kojima et al. 2019, Izotov et al. 2018), enabling us to derive element abundance ratios with photoionization models. We find that neon- and argon-to-oxygen ratios are comparable to those of known local dwarf galaxies, and that the nitrogen-to-oxygen abundance ratios (N/O) are lower than 20% (N/O)_sun consistent with the low oxygen abundance. However, the iron-to-oxygen abundance ratios (Fe/O) of the EMPGs are generally high; the EMPGs with the 2%-solar oxygen abundance show high Fe/O ratios of ~90-140% (Fe/O)_sun, which are unlikely explained by suggested scenarios of Type Ia supernova iron productions, iron's dust depletion, and metal-poor gas inflow onto previously metal-riched galaxies with solar abundances. Moreover, these EMPGs have very high HeII4686/H$\beta$ ratios of ~1/40, which are not reproduced by existing models of high-mass X-ray binaries whose progenitor stellar masses are less than 120 M_sun. Comparing stellar-nucleosynthesis and photoionization models with a comprehensive sample of EMPGs identified by this and previous EMPG studies, we propose that both the high Fe/O ratios and the high HeII4686/H$\beta$ ratios are explained by the past existence of super massive ($>$300 M_sun) stars, which may evolve into intermediate-mass black holes ($\gtrsim$100 M_sun).

We present the first stellar velocity dispersion measurement of a massive quenching galaxy at z=4.01. The galaxy is first identified as a massive z>~4 galaxy with suppressed star formation from photometric redshifts based on deep multi-band data in the UKIDSS Ultra Deep Survey field. A follow-up spectroscopic observation with MOSFIRE on Keck revealed strong multiple absorption features, which are identified as Balmer absorption lines, giving a secure redshift of z=4.01. Thanks to the high S/N of the spectrum, we are able to estimate the stellar velocity dispersion, sigma=268+/-59 km/s. This velocity dispersion is consistent with that of massive galaxies today, implying no significant evolution in stellar velocity dispersion over the last 12 Gyr. Based on an upper limit on its physical size from deep optical images (r_eff<1.3 kpc), we find that its dynamical mass is consistent with the stellar mass inferred from photometry. Furthermore, the galaxy is located on the mass fundamental plane extrapolated from lower redshift galaxies. Combining all these results, we find that the velocity dispersion does not significantly evolve with redshift, although the size and mass of massive quenched galaxies do. This suggests that the mass in the core of massive galaxies does not evolve significantly, while most of the mass growth occurs in the outskirts of the galaxies, which also increases the size. This picture is consistent with a two-phase formation scenario in which mass and size growth is due to accretion in the outskirts of galaxies via mergers.

We present the morphology and stellar population of 27 extremely metal poor galaxies (EMPGs) at $z\sim0$ with metallicities of $0.01-0.1$ Z$_{\odot}$. We conduct multi-component surface brightness (SB) profile fitting for the deep Subaru/HSC $i$-band images of the EMPGs with the Galfit software, carefully removing the SB contributions of potential associated galaxies (PAGs). We find that the EMPGs with a median stellar mass of $\log(M_{*}/{\rm M}_{\odot})=5.77$ have a median S\'ersic index of $n=1.08$ and a median effective radius of $r_{\rm e}=176$ pc, suggesting that typical EMPGs are very compact disk galaxies. We compare the EMPGs with $z\sim 6$ galaxies and local galaxies on the size-mass ($r_{\rm e}$--$M_*$) diagram, and identify that the majority of the EMPGs have a $r_{\rm e}$--$M_*$ relation similar to $z\sim0$ star-forming galaxies rather than $z\sim6$ galaxies. Not every EMPG is a local analog of high-$z$ young galaxies in the $r_{\rm e}$--$M_*$ relation. We also study the PAGs of our galaxies, and find that 23 out of the 27 EMPGs show detectable PAGs within the projected distance of 10 kpc. The PAGs have median values of $n=0.93$, $r_{\rm e}=1.41$ kpc, and $\log(M_{*}/{\rm M}_{\odot})=7.47$ that are similar to those of local dwarf irregulars and ultra-diffuse galaxies. A spectrum of one EMPG-PAG system, so far available, indicates that the PAG is dynamically related to the EMPG with a median velocity difference of $\Delta V=96.1$ km s$^{-1}$. This moderately-large $\Delta V$ cannot be explained by the dynamics of the PAG, but likely by the infall on the PAG. EMPGs may form in infalling gas, and become part of PAG disks undergoing further metal enrichment.

Star-forming galaxies with strong nebular and collisional emission lines are privileged target galaxies in forthcoming cosmological large galaxy redshift surveys. We use the COSMOS2015 photometric catalog to model galaxy spectral energy distributions and emission-line fluxes. We adopt an empirical but physically-motivated model that uses information from the best-fitting spectral energy distribution of stellar continuum to each galaxy. The emission-line flux model is calibrated and validated against direct flux measurements in subsets of galaxies that have 3D-HST or zCOSMOS-Bright spectra. We take a particular care in modelling dust attenuation such that our model can explain both H$\alpha$ and [OII] observed fluxes at different redshifts. We find that a simple solution to this is to introduce a redshift evolution in the dust attenuation fraction parameter, $f=E_{\rm star}(B-V)/E_{\rm gas}(B-V)$, as $f(z)=0.44+0.2z$. From this catalog, we derive the H$\alpha$ and [OII] luminosity functions up to redshifts of about 2.5 after carefully accounting for emission line flux and redshift errors. This allows us to make predictions for H$\alpha$ and [OII] galaxy number counts in next-generation cosmological redshift surveys. Our modeled emission lines and spectra in the COSMOS2015 catalog shall be useful to study the target selection for planned next-generation galaxy redshift surveys and we make them publicly available as `EL-COSMOS' on the ASPIC database.

We report two secure ($z=3.775, 4.012$) and one tentative ($z\approx3.767$) spectroscopic confirmations of massive and quiescent galaxies through $K$-band observations with Keck/MOSFIRE and VLT/X-Shooter. The stellar continuum emission, the absence of strong nebular emission lines and the lack of significant far-infrared detections confirm the passive nature of these objects, disfavoring the alternative solution of low-redshift dusty star-forming interlopers. We derive stellar masses of $\mathrm{log}(M_{\star}/M_\odot)\sim11$ and ongoing star formation rates placing these galaxies $\gtrsim 1-2$ dex below the main sequence at their redshifts. The adopted parametrization of the star formation history suggests that these sources experienced a strong ($\langle \rm SFR \rangle \sim 1200-3500\,M_\odot\,\mathrm{yr}^{-1}$) and short ($\sim 50$ Myr) burst of star formation, peaking $\sim 150-500$ Myr before the time of observation, all properties reminiscent of the characteristics of sub-millimeter galaxies (SMGs) at $z>4$. We investigate this connection by comparing the comoving number densities and the properties of these two populations. We find a fair agreement only with the deepest sub-mm surveys detecting not only the most extreme starbursts, but also more normal galaxies. We support these findings by further exploring the Illustris-TNG cosmological simulation, retrieving populations of both fully quenched massive galaxies at $z\sim3-4$ and SMGs at $z\sim4-5$, with number densities and properties in agreement with the observations at $z\sim3$, but in increasing tension at higher redshift. Nevertheless, as suggested by the observations, not all the progenitors of quiescent galaxies at these redshifts shine as bright SMGs in their past and, similarly, not all bright SMGs quench by $z\sim3$, both fractions depending on the threshold assumed to define the SMGs themselves.

We examine the impact of fiber assignment on clustering measurements from fiber-fed spectroscopic galaxy surveys. We identify new effects which were absent in previous, relatively shallow galaxy surveys such as Baryon Oscillation Spectroscopic Survey . Specifically, we consider deep surveys covering a wide redshift range from z=0.6 to z=2.4, as in the Subaru Prime Focus Spectrograph survey. Such surveys will have more target galaxies than we can place fibers on. This leads to two effects. First, it eliminates fluctuations with wavelengths longer than the size of the field of view, as the number of observed galaxies per field is nearly fixed to the number of available fibers. We find that we can recover the long-wavelength fluctuation by weighting galaxies in each field by the number of target galaxies. Second, it makes the preferential selection of galaxies in under-dense regions. We mitigate this effect by weighting galaxies using the so-called individual inverse probability. Correcting these two effects, we recover the underlying correlation function at better than 1 percent accuracy on scales greater than 10 Mpc/h.

We present the rest-frame optical sizes of massive quiescent galaxies (QGs) at $z\sim4$ measured at $K'$-band with the Infrared Camera and Spectrograph (IRCS) and AO188 on the Subaru telescope. Based on a deep multi-wavelength catalog in the Subaru XMM-Newton Deep Survey Field (SXDS), covering a wide wavelength range from the $u$-band to the IRAC $8.0\mu m$ over 0.7 deg$^2$, we evaluate photometric redshift to identify massive ($M_{\star}\sim10^{11}\ M_\odot$) galaxies with suppressed star formation. These galaxies show a prominent 4000$\rm \AA$ break feature at $z\sim4$, suggestive of an evolved stellar population. We then conduct follow-up $K'$-band imaging with adaptive optics for the five brightest galaxies ($K_{AB,total}=22.5\sim23.4$). Compared to lower redshift ones, QGs at $z\sim4$ have smaller physical sizes of effective radii $r_{eff}=0.2$ to $1.8$ kpc. The mean size measured by stacking the four brightest objects is $r_{eff}=0.7\rm\ kpc$. This is the first measurement of the rest-frame optical sizes of QGs at $z\sim4$. We evaluate the robustness of our size measurements using simulations and find that our size estimates are reasonably accurate with an expected systematic bias of $\sim0.2$ kpc. If we account for the stellar mass evolution, massive QGs at $z\sim4$ are likely to evolve into the most massive galaxies today. We find their size evolution with cosmic time in a form of $\log(r_e/{\rm kpc})= -0.44+1.77 \log(t/\rm Gyr)$. Their size growth is proportional to the square of stellar mass, indicating the size-stellar mass growth driven by minor dry mergers.

Intrinsic alignments (IA), the coherent alignment of intrinsic galaxy orientations, can be a source of a systematic error of weak lensing surveys. The redshift evolution of IA also contains information about the physics of galaxy formation and evolution. This paper presents the first measurement of IA at high redshift, $z\sim 1.4$, using the spectroscopic catalog of blue star-forming galaxies of the FastSound redshift survey, with the galaxy shape information from the Canada-Hawaii-France telescope lensing survey. The IA signal is consistent with zero with power-law amplitudes fitted to the projected correlation functions for density-shape and shape-shape correlation components, $A_{\delta+}=-0.0071\pm 0.1340$ and $A_{++}=-0.0505\pm 0.0848$, respectively. These results are consistent with those obtained from blue galaxies at lower redshifts (e.g., $A_{\delta+}=0.0035_{-0.0389}^{+0.0387}$ and $A_{++}=0.0045_{-0.0168}^{+0.0166}$ at $z=0.51$ from the WiggleZ survey). The upper limit of the constrained IA amplitude corresponds to a few percent contamination to the weak-lensing shear power spectrum, resulting in systematic uncertainties on the cosmological parameter estimations by $-0.052<\Delta \sigma_8<0.039$ and $-0.039<\Delta \Omega_m<0.030$.

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

We report the study of far-IR sizes of submillimeter galaxies (SMGs) in relation to their dust-obscured star formation rate (SFR) and active galactic nuclei (AGN) presence, determined using mid-IR photometry. We determined the millimeter-wave ($\lambda_{\rm obs}=1100 \mu$m) sizes of 69 ALMA-identified SMGs, selected with $\geq10$$\sigma$ confidence on ALMA images ($F_{\rm 1100 \mu m}=1.7$--7.4 mJy). We found that all the SMGs are located above an avoidance region in the millimeter size-flux plane, as expected by the Eddington limit for star formation. In order to understand what drives the different millimeter-wave sizes in SMGs, we investigated the relation between millimeter-wave size and AGN fraction for 25 of our SMGs at $z=1$--3. We found that the SMGs for which the mid-IR emission is dominated by star formation or AGN have extended millimeter-sizes, with respective median $R_{\rm c,e} = 1.6^{+0.34}_{-0.21}$ and 1.5$^{+0.93}_{-0.24}$ kpc. Instead, the SMGs for which the mid-IR emission corresponds to star-forming/AGN composites have more compact millimeter-wave sizes, with median $R_{\rm c,e}=1.0^{+0.20}_{-0.20}$ kpc. The relation between millimeter-wave size and AGN fraction suggests that this size may be related to the evolutionary stage of the SMG. The very compact sizes for composite star-forming/AGN systems could be explained by supermassive black holes growing rapidly during the SMG coalescing, star-formation phase.

We present initial results from the Subaru Strategic Program (SSP) with Hyper Suprime-Cam (HSC) on a comprehensive survey of emission-line galaxies at z<1.5 based on narrowband (NB) imaging. The first Public Data Release (PDR1) provides us with data from two NB filters, specifically NB816 and NB921 over 5.7 deg$^2$ and 16.2 deg$^2$ respectively. The $5 \sigma$ limiting magnitudes are 25.2 (UDeep layer, 1.4 deg$^2$) and 24.8 (Deep layer, 4.3 deg$^2$) mag in NB816, and 25.1 (UDeep, 2.9 deg$^2$) and 24.6--24.8 (Deep, 13.3 deg$^2$) mag in NB921. The wide-field imaging allows us to construct unprecedentedly large samples of 8,054 H$\alpha$ emitters at z ~ 0.25 and 0.40, 8,656 [OIII] emitters at z ~ 0.63 and 0.84, and 16,877 [OII] emitters at z ~ 1.19 and 1.47. We map the cosmic web on scales out to about 50 comoving Mpc that includes galaxy clusters, identified by red sequence galaxies, located at the intersection of filamentary structures of star-forming galaxies. The luminosity functions of emission-line galaxies are measured with precision and consistent with published studies. The wide field coverage of the data enables us to measure the luminosity functions up to brighter luminosities than previous studies. The comparison of the luminosity functions between the different HSC-SSP fields suggests that a survey volume of $>5\times10^5$ Mpc$^3$ is essential to overcome cosmic variance. Since the current data have not reached the full depth expected for the HSC-SSP, the color cut in i-NB816 or z-NB921 induces a bias towards star-forming galaxies with large equivalent widths, primarily seen in the stellar mass functions for the H$\alpha$ emitters at z ~ 0.25--0.40. Even so, the emission-line galaxies clearly cover a wide range of luminosity, stellar mass, and environment, thus demonstrating the usefulness of the NB data from the HSC-SSP to investigate star-forming galaxies at z<1.5.

We present the environmental dependence of colour, stellar mass, and star formation (SF) activity in H-alpha-selected galaxies along the large-scale structure at z=0.4 hosting twin clusters in DEEP2-3 field, discovered by Subaru Strategic Programme of Hyper Suprime-Cam (HSC SSP). By combining photo-z selected galaxies and H-alpha emitters selected with broad-band and narrow-band (NB) data of the recent data release of HSC SSP (DR1), we confirm that galaxies in higher-density environments or galaxies in the cluster central regions show redder colours. We find that there still remains a possible colour-density and colour-radius correlation even if we restrict the sample to H-alpha-selected galaxies, likely due to the presence of massive H-alpha emitters in denser regions. We also find a hint of increased star formation rates (SFR) amongst H-alpha emitters towards the highest-density environment, again primarily driven by the excess of red/massive H-alpha emitters in high-density environment, while their specific SFR does not significantly change with environment. This work demonstrates the power of the HSC SSP NB data to study SF galaxies across environment in the distant universe.

We discuss the host galaxy metallicity distribution of all long gamma-ray bursts (GRBs) whose redshifts are known to be $< 0.4$, including newly obtained spectroscopic datasets of the host galaxies of GRB 060614, 090417B, and 130427A. We compare the metallicity distribution of the low-redshift sample to the model predictions, and constrain the relation between metallicity and GRB occurrence. We take account of spatial variation of metallicities among star forming regions within a galaxy. We found that the models, in which only low-metallicity stars produce GRBs with a sharp cutoff of GRB production efficiency around 12+log(O/H) $\sim$ 8.3, can well reproduce the observed distribution, while the models with no metallicity dependence are not consistent with the observations. We also discuss possible sampling biases we may suffer by collecting long GRBs whose redshifts are known, presenting the photometric observations of the host galaxy of GRB 111225A at $z = 0.297$ whose redshift has been undetermined until $\sim$ 2.3 years after the burst.

We present the detailed characterization of two extremely red submillimeter galaxies (SMGs), ASXDF1100.053.1 and 231.1, with the Atacama Large Millimeter/submillimeter Array (ALMA) and the Jansky Very Large Array (VLA). These SMGs were selected originally using AzTEC at 1100 micron, and are observed by Herschel to be faint at 100--500 micron. Their (sub)millimeter colors are as red as -- or redder -- than known z>~5 SMGs; indeed, ASXDF1100.053.1 is redder than HFLS 3, which lies at z=6.3. They are also faint and red in the near-/mid-infrared: ~1 microJy at IRAC 4.5 micron and <0.2 microJy in the Ks filter. These SMGs are also faint in the radio waveband, where F_6GHz=4.5 microJy for ASXDF1100.053.1 and F_1.4GHz=28 microJy for ASXDF1100.231.1, suggestive of z=6.5^{+1.4}_{-1.1} and z=4.1^{+0.6}_{-0.7} for ASXDF1100.053.1 and 231.1, respectively. ASXDF1100.231.1 has a flux excess in the 3.6-micron filter, probably due to H$\alpha$ emission at z=4--5. Derived properties of ASXDF1100.053.1 for z=5.5--7.5 and 231.1 for z=3.5--5.5 are as follows: their infrared luminosities are [6.5-7.4]x10^{12} and [4.2-4.5]x10^{12} L_sun; their stellar masses are [0.9-2]x10^{11} and [0.4-3]x10^{10} M_sun; their circularized half-light radii in the ALMA maps are ~1 and <~0.2 kpc (~2--3 kpc for 90% of the total flux). Lastly, their surface infrared luminosity densities, Sigma_IR, are ~1x10^{12} and >~1.5x10^{13} L_sun kpc^{-2}, similar to values seen for local (U)LIRGs. These data suggest that ASXDF1100.053.1 and 231.1 are compact SMGs at z>~4 and can plausibly evolve into z>~3 compact quiescent galaxies.

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.

The Prime Focus Spectrograph (PFS) is a wide-field, multi-object spectrograph accommodating 2394 fibers to observe the sky at the prime focus of the Subaru telescope. The software system to operate a spectroscopic survey is structured by the four packages: Instrument control software, exposure targeting software, data reduction pipeline, and survey planning and tracking software. In addition, we operate a database system where various information such as properties of target objects, instrument configurations, and observation conditions is stored and is organized via a standardized data model for future references to update survey plans and to scientific researches. In this article, we present an overview of the software system and describe the workflows that need to be performed in the PFS operation, with some highlights on the database that organizes various information from sub-processes in the survey operation, and on the process of fiber configuration from the software perspectives.

PFS (Prime Focus Spectrograph), a next generation facility instrument on the 8.2-meter Subaru Telescope, is a very wide-field, massively multiplexed, optical and near-infrared spectrograph. Exploiting the Subaru prime focus, 2394 reconfigurable fibers will be distributed over the 1.3 deg field of view. The spectrograph has been designed with 3 arms of blue, red, and near-infrared cameras to simultaneously observe spectra from 380nm to 1260nm in one exposure at a resolution of ~1.6-2.7A. An international collaboration is developing this instrument under the initiative of Kavli IPMU. The project is now going into the construction phase aiming at undertaking system integration in 2017-2018 and subsequently carrying out engineering operations in 2018-2019. This article gives an overview of the instrument, current project status and future paths forward.

We present a multi-wavelength analysis of five submillimeter sources (S_1.1mm = 0.54-2.02 mJy) that were detected during our 1.1-mm-deep continuum survey in the SXDF-UDS-CANDELS field (2 arcmin^2, 1sigma = 0.055 mJy beam^-1) using the Atacama Large Millimeter/submillimeter Array (ALMA). The two brightest sources correspond to a known single-dish (AzTEC) selected bright submillimeter galaxy (SMG), whereas the remaining three are faint SMGs newly uncovered by ALMA. If we exclude the two brightest sources, the contribution of the ALMA-detected faint SMGs to the infrared extragalactic background light is estimated to be ~ 4.1^{+5.4}_{-3.0} Jy deg^{-2}, which corresponds to ~ 16^{+22}_{-12}% of the infrared extragalactic background light. This suggests that their contribution to the infrared extragalactic background light is as large as that of bright SMGs. We identified multi-wavelength counterparts of the five ALMA sources. One of the sources (SXDF-ALMA3) is extremely faint in the optical to near-infrared region despite its infrared luminosity (L_IR ~ 1e12 L_sun or SFR ~ 100 M_sun yr^{-1}). By fitting the spectral energy distributions (SEDs) at the optical-to-near-infrared wavelengths of the remaining four ALMA sources, we obtained the photometric redshifts (z_photo) and stellar masses (M_*): z_photo ~ 1.3-2.5, M_* ~ (3.5-9.5)e10 M_sun. We also derived their star formation rates (SFRs) and specific SFRs (sSFRs) as ~ 30-200 M_sun yr^{-1} and ~ 0.8-2 Gyr^{-1}, respectively. These values imply that they are main-sequence star-forming galaxies.

We present results of 12CO(J=2-1) observations toward four massive star-forming galaxies at z~1.4 with the Nobeyama 45~m radio telescope. The galaxies are detected with Spitzer/MIPS in 24 um, Herschel/SPIRE in 250 um, and 350 um and they mostly reside in the main sequence. Their gas-phase metallicities derived with N2 method by using the Ha and [NII]6584 emission lines are near the solar value. CO lines are detected toward three galaxies. The molecular gas masses obtained are (9.6-35) x 10^{10} Msun by adopting the Galactic CO-to-H2 conversion factor and the CO(2-1)/CO(1-0) flux ratio of 3. The dust masses derived with the modified blackbody model (assuming the dust temperature of 35 K and the emissivity index of 1.5) are (2.4-5.4) x 10^{8} Msun. The resulting gas-to-dust ratios (not accounting for HI mass) at z~1.4 are 220-1450, which are several times larger than those in local star-forming galaxies. A dependence of the gas-to-dust ratio on the far-infrared luminosity density is not clearly seen.

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.

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.

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.

We report the source size distribution, as measured by ALMA millimetric continuum imaging, of a sample of 13 AzTEC-selected submillimeter galaxies (SMGs) at z_photo ~ 3-6. Their infrared luminosities and star-formation rates (SFR) are L_IR ~ 2-6 x 10^12 L_sun and ~ 200-600 M_sun yr-1, respectively. The size of z ~ 3-6 SMGs ranges from 0".10 to 0".38 with a median of 0".20+0".03-0".05 (FWHM), corresponding to a median circularized effective radius (Rc,e) of 0.67+0.13-0.14 kpc, comparable to the typical size of the stellar component measured in compact quiescent galaxies at z ~ 2 (cQGs) --- R ~ 1 kpc. The median surface SFR density of our z ~ 3-6 SMGs is 100+42-26 M_sun yr-1 kpc-2, comparable to that seen in local merger-driven (U)LIRGsrather than in extended disk galaxies at low and high redshifts. The discovery of compact starbursts in z >~ 3 SMGs strongly supports a massive galaxy formation scenario wherein z ~ 3-6 SMGs evolve into the compact stellar components of z ~ 2 cQGs. These cQGs are then thought to evolve into the most massive ellipticals in the local Universe, mostly via dry mergers. Our results thus suggest that z >~ 3 SMGs are the likely progenitors of massive local ellipticals, via cQGs, meaning that we can now trace the evolutionary path of the most massive galaxies over a period encompassing ~ 90% of the age of the Universe.

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.

When and how galaxy morphology such as disk and bulge seen in the present-day universe emerged is still not clear. In the universe at $z\gtrsim 2$, galaxies with various morphology are seen, and star-forming galaxies at $z\sim2$ show an intrinsic shape of bar-like structure. Then, when did round disk structure form? Here we take a simple and straightforward approach to see the epoch when a round disk galaxy population emerged by constraining the intrinsic shape statistically based on apparent axial ratio distribution of galaxies. We derived the distributions of the apparent axial ratios in the rest-frame optical light ($\sim 5000$ \AA) of star-forming main sequence galaxies at $2.5>z>1.4$, $1.4>z>0.85$, and $0.85>z>0.5$, and found that the apparent axial ratios of them show peaky distributions at $z\gtrsim0.85$, while a rather flat distribution at the lower redshift. By using a tri-axial model ($A>B>C$) for the intrinsic shape, we found the best-fit models give the peaks of the $B/A$ distribution of $0.81\pm0.04$, $0.84\pm0.04$, and $0.92\pm0.05$ at $2.5>z>1.4$, $1.4>z>0.85$, and $0.85>z>0.5$, respectively. The last value is close to the local value of 0.95. Thickness ($C/A$) is $\sim0.25$ at all the redshifts and is close to the local value (0.21). The results indicate the shape of the star-forming galaxies in the main sequence changes gradually, and the round disk is established at around $z\sim0.9$. Establishment of the round disk may be due to a cease of violent interaction of galaxies or a growth of a bulge and/or a super-massive black hole resides at the center of a galaxy which dissolves the bar structure.

We present near-infrared spectroscopy of the host galaxy of dark GRB 080325 using Subaru/MOIRCS. The obtained spectrum provides a clear detection of H$\alpha$ emission and marginal [NII]$\lambda$6584. The host is a massive (M$_{*}\sim10^{11}$M$_{\odot}$), dusty ($A_{V}\sim 1.2$) star-forming galaxy at z=1.78. The star formation rate calculated from the H$\alpha$ luminosity (35.6-47.0 M$_{\odot}$ yr$^{-1}$) is typical among GRB host galaxies (and star-forming galaxies generally) at z $>$1; however, the specific star formation rate is lower than normal star-forming galaxies at redshift $\sim$ 1.6, in contrast to the high specific star formation rates measured for many of other GRB hosts. The metallicity of the host is estimated to be 12+log(O/H)$_{\rm KK04}$$=$8.88. We emphasize that this is one of the most massive distant host galaxies for which metallcity is measured with emission-line diagnostics. The metallicity is fairly high among GRB hosts. However, this is still lower than the metallicity of normal star-forming galaxies of the same mass at z$\sim$1.6. The metallicity offset from normal star-forming galaxies is close to a typical value of other GRB hosts and indicates that GRB host galaxies are uniformly biased toward low metalicity over a wide range of redshift and stellar mass. The low-metallicity nature of the GRB 080325 host is likely not attributable to the fundamental metallicity relation of star-forming galaxies beacuse it is a metal-poor outlier from the relation and has a low sSFR. Thus we conclude that metallicity is important to the mechanism that produced this GRB.

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 carried out 12CO(J=2-1) observations toward three star-forming galaxies on the main sequence at z~1.4 with the Nobeyama 45m radio telescope. These galaxies are detected with Spitzer/MIPS in 24 um, Herschel/SPIRE in 250 um and 350 um, and their gas metallicity, derived from optical emission line ratios based on near infrared spectroscopic observations, is close to the solar metallicity. Although weak signal-like features of CO were seen, we could not detect significant CO emission. The dust mass and the upper limits on the molecular gas mass are (3.4-6.7) x 10^{8} Msun and (9.7-14) x 10^{10} Msun, respectively. The upper limits on the gas-to-dust ratios at z~1.4 are 150-410 which are comparable to the gas-to-dust ratios in local galaxies with similar gas metallicity. A line stacking analysis enables us to detect a significant CO emission and to derive an average molecular gas mass of 1.3 x 10^{11} Msun and gas-to-dust ratio of 250. This gas-to-dust ratio is also near that in local galaxies with solar metallicity. These results suggest that the gas-to-dust ratio in star-forming galaxies with solar metallicity does not evolve significantly up to z~1.4. By comparing to a theoretical calculation, a rapid increase of the dust mass in an earlier epoch of galaxy evolution is suggested.

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

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

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.

(Abridge) We study statistical intrinsic shape of star-forming BzK galaxies (sBzK galaxies) at z~2 in both rest-frame UV and rest-frame optical wavelengths. The sBzK galaxies are selected down to K(AB)=24.0 mag in the GOODS-South and SXDS fields, where high-resolution images from Hubble Space Telescope are publicly available. 57% (583) of all 1028 galaxies in GOODS-S show a single component in the ACS/F850LP image. As WFC3/F160W images cover only some part of GOODS-S and SXDS, 724/1028 and 2500/29835 sBzK galaxies in the GOODS-S and SXDS have the WFC3 coverage. 86% (626) and 82% (2044) of the sBzK galaxies in WFC3/F160W images appear as a single component in the GOODS-S and SXDS, respectively. Larger fraction of single-component objects in F850LP images represents multiple star-forming regions in galaxies, while they are not so obvious in the F160W image which appears smoother. Most of the single-component sBzK galaxies show S\'ersic indices of n=0.5-2.5, in agreement with those of local disk galaxies. Their effective radii are 1.0-3.0 kpc and 1.5-4.0 kpc in F850LP and F160W images, respectively, regardless of the observed fields. Stellar surface mass density of the sBzK galaxies is also comparable to that of the local disk galaxies. However, the intrinsic shape of sBzK galaxies is not a round disk as seen in the local disk galaxies. By comparing apparent axial ratio (b/a) distributions of the sBzK galaxies with those by assuming tri-axial model with axes A>B>C, we found their intrinsic face-on B/A ratios peak at B/A=0.70 and B/A=0.77-0.79 in the rest-frame UV and optical, respectively and are statistically more bar-like than that of the local disk galaxies. The intrinsic edge-on C/A ratios in both rest-frame UV and optical wavelengths peak at 0.26, which is slightly larger than that of the local disk galaxies.

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 initial results from the CO survey toward high redshift galaxies using the Nobeyama 45m telescope. Using the new wide bandwidth spectrometer equipped with a two-beam SIS receiver, we have robust new detections of three high redshift (z=1.6-3.4) submillimeter galaxies (SXDF 1100.001, SDP9, and SDP17), one tentative detection (SDSS J160705+533558), and one non-detection (COSMOS-AzTEC1). The galaxies observed during the commissioning phase are sources with known spectroscopic redshifts from previous optical or from wide-band submm spectroscopy. The derived molecular gas mass and line widths from Gaussian fits are ~10^11 Msun and 430-530 km/s, which are consistent with previous CO observations of distant submm galaxies and quasars. The spectrometer that allows a maximum of 32 GHz instantaneous bandwidth will provide new science capabilities at the Nobeyama 45m telescope, allowing us to determine redshifts of bright submm selected galaxies without any prior redshift information.

We present results of a search for a supernova (SN) component associated with GRB 100418A at the redshift of 0.624. The field of GRB 100418A was observed with FOCAS on Subaru 8.2m telescope under a photometric condition (seeing 0.3"-0.4") on 2010 May 14 (UT). The date corresponds to 25.6 days after the burst trigger (15.8 days in the restframe). We did imaging observations in V, Rc, and Ic bands, and two hours of spectrophotometric observations. We got the resolved host galaxy image which elongated 1.6" (= 11 kpc) from north to south. No point source was detected on the host galaxy. The time variation of Rc-band magnitude shows that the afterglow of GRB 100418A has faded to Rc \sim > 24 without SN like rebrightening, when we compare our measurement to the reports in GCN circulars. We could not identify any SN feature such as broad emission-lines or bumps in our spectrum. Assuming the SN is fainter than the 3{\sigma} noise spectrum of our observation, we estimate the upper limit on the SN absolute magnitude MIc,obs > -17.2 in observer frame Ic-band. This magnitude is comparable to the faintest type Ic SNe. We also estimate host galaxy properties from the spectrum. The host galaxy of GRB 100418A is relatively massive (log M_{star}/M_{sun} = 9.54) compared to typical long GRB host galaxies, and has 12+log(O/H) = 8.75.

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 structure of star-forming galaxies at z~2 in GOODS-N field selected as sBzK galaxies down to K_{AB} <24.0 mag. Among 1029 sBzK galaxies, 551 galaxies (54%) show a single component in ACS/F850LP image obtained with the Hubble Space Telescope; the rest show multiple components. We fit the single-component sBzK galaxies with the single S\'ersic profile using the ACS/F850LP image and find that a majority of them (64%) show S\'ersic index of n=0.5-2.5, indicating that they have a disk-like structure. The resulting effective radii typically range from 1.0 to 3.0 kpc in the rest-frame UV wavelength. After correcting the effective radii to those in the rest-frame optical wavelength, we find that the single-component sBzK galaxies locate in the region where the local and z~1 disk galaxies distribute in the stellar mass-size diagram, suggesting comparable surface stellar mass density between the sBzK and z~0-1 disk galaxies. All these properties suggest that the single-component sBzK galaxies are progenitors of the present-day disk galaxies. However, by studying their intrinsic shape through comparison between the observed distribution of apparent axial ratios and the distribution for triaxial models with axes (A>B>C), we find that the mean B/A ratio is 0.61^{+0.05}_{-0.08} and disk thickness C/A is 0.28^{+0.03}_{-0.04}. This indicates that the single-component sBzK galaxies at z~2 have a bar-like or oval shape rather than a round disk shape. The shape seems to resemble to a bar/oval structure that form through bar instability; if it is the case, the intrinsic shape may give us a clue to understand dynamical evolution of baryonic matter in a dark matter halo.

(abridged) We present a study of stellar population of LAEs at z=4.86 in GOODS-N and its flanking field. With the publicly available IRAC data in GOODS-N and further IRAC observations in the flanking fields, we select five LAEs which are not contaminated by neighboring objects in IRAC images and construct their observed SEDs with I_c, z', IRAC 3.6micron, and 4.5micron band photometry. The SEDs cover the rest-frame UV to optical wavelengths. We derive stellar masses, ages, color excesses, and star formation rates of five LAEs using SED fitting method. Assuming the constant star formation history, we find that the stellar masses range from 10^8 to $10^{10} Msun with the median value of 2.5x10^9 Msun. The derived ages range from very young ages (7.4 Myr) to 437 Myr with a median age of 25 Myr. The color excess E(B-V) are between 0.1-0.4 mag. Star formation rates are 55-209 Msun/yr. A comparison of the stellar populations is made between three LAEs and 88 LBGs selected at the same redshift, in the same observed field, and down to the same limit of the rest-frame UV luminosity. These three LAEs are the brightest and reddest samples among the whole LAE samples at z=4.86. The LAEs distribute at the relatively faint part of UV-luminosity distribution of LBGs. Deriving the stellar properties of the LBGs by fitting their SEDs with the same model ensures that model difference does not affect the comparison. It is found that the stellar properties of the LAEs lie on distributions of those of LBGs. On average, the LAEs show less dust extinction, and lower star formation rates than LBGs, while the stellar mass of LAEs nearly lies in the middle part of the mass distribution of LBGs. However, the stellar properties of LAEs and LBGs are similar at the fixed UV or optical luminosity. We also examine the relations between the output properties from the SED fitting and the rest-frame Lya equivalent width.

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 optical and near infrared observations of GRB 080325 classified as a "Dark GRB". Near-infrared observations with Subaru/MOIRCS provided a clear detection of afterglow in Ks band, although no optical counterpart was reported. The flux ratio of rest-wavelength optical to X-ray bands of the afterglow indicates that the dust extinction along the line of sight to the afterglow is Av = 2.7 - 10 mag. This large extinction is probably the major reason for optical faintness of GRB 080325. The J - Ks color of the host galaxy, (J - Ks = 1.3 in AB magnitude), is significantly redder than those for typical GRB hosts previously identified. In addition to J and Ks bands, optical images in B, Rc, i', and z' bands with Subaru/Suprime-Cam were obtained at about one year after the burst, and a photometric redshift of the host is estimated to be z_{photo} = 1.9. The host luminosity is comparable to L^{*} at z \sim 2 in contrast to the sub-L^{*} property of typical GRB hosts at lower redshifts. The best-fit stellar population synthesis model for the host shows that a large dust extinction (Av = 0.8 mag) attributes to the red nature of the host and that the host galaxy is massive (M_{*} = 7.0 \times 10^{10} Msun) which is one of the most massive GRB hosts previously identified. By assuming that the mass-metallicity relation for star-forming galaxies at z \sim 2 is applicable for the GRB host, this large stellar mass suggests the high metallicity environment around GRB 080325, consistent with inferred large extinction.

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

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