We present the results of the deep Subaru/FOCAS and Keck/MOSFIRE spectroscopy for four spatially extended [O III]$\lambda\lambda$4959,5007 sources, dubbed [O III] blobs, at $z=0.6-0.8$ that are originally pinpointed by large-area Subaru imaging surveys. The line diagnostics of the rest-frame optical lines suggests that only one [O III] blob, OIIIB-3, presents an AGN signature, indicating that hot gas of the rest of the [O III] blobs is heated by star formation. One of such star-forming [O III] blobs, OIIIB-4, at $z=0.838$ has an [O III] equivalent width of $845\pm27$ \r{A} and an [O III] to [O II]$\lambda\lambda$3726,3729 ratio of [O III]/[O II]= $6.5\pm2.7$ that are as high as those of typical green peas (Cardamone et al. 2009). The spatially resolved spectrum of OIIIB-4 shows [O III]/[O II]= $5-10$ over $14$ kpc in the entire large [O III] extended regions of OIIIB-4, unlike the known green peas whose strong [O III] emission region is compact. Moreover, OIIIB-4 presents no high ionization emission lines, unlike green beans that have extended [O III] emission with a type-2 AGN. OIIIB-4 is thus a giant green pea, which is a low stellar mass ($7\times10^7$ $M_\odot$) galaxy with a very high specific star formation rate (sSFR = $2\times10^2\, {\rm Gyr}^{-1}$), a high ionization parameter ($q_{ion} \sim 3\times10^8\, {\rm cm\,s^{-1}}$), and a low metallicity similar to those of green peas. Neither an AGN-light echo nor a fast radiative shock likely takes place due to the line diagnostics for spatially-resolved components of OIIIB-4 and no detections of He II$\lambda$4686 or [Ne V]$\lambda\lambda$3346,3426 lines that are fast-radiative shock signatures. There is a possibility that the spatially-extended [O III] emission of OIIIB-4 is originated from outflowing gas produced by the intense star formation in a density-bounded ionization state.

We present 45 ground-based photometric observations of the K2-22 system collected between December 2016 and May 2017, which we use to investigate the evolution of the transit of the disintegrating planet K2-22b. Last observed in early 2015, in these new observations we recover the transit at multiple epochs and measure a typical depth of <1.5%. We find that the distribution of our measured transit depths is comparable to the range of depths measured in observations from 2014 and 2015. These new observations also support ongoing variability in the K2-22b transit shape and time, although the overall shallowness of the transit makes a detailed analysis of these transit parameters difficult. We find no strong evidence of wavelength-dependent transit depths for epochs where we have simultaneous coverage at multiple wavelengths, although our stacked Las Cumbres Observatory data collected over days-to-months timescales are suggestive of a deeper transit at blue wavelengths. We encourage continued high-precision photometric and spectroscopic monitoring of this system in order to further constrain the evolution timescale and to aid comparative studies with the other few known disintegrating planets.

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

We present an unprecedentedly large catalog consisting of 2,230 > L* Lya emitters (LAEs) at z=5.7 and 6.6 on the 13.8 and 21.2 deg2 sky, respectively, that are identified by the SILVERRUSH program with the first narrowband imaging data of the Hyper Suprime-Cam (HSC) survey. We confirm that the LAE catalog is reliable on the basis of 96 LAEs whose spectroscopic redshifts are already determined by this program and the previous studies. This catalogue is also available on-line. Based on this catalogue, we derive the rest-frame Lya equivalent-width distributions of LAEs at z~5.7-6.6 that are reasonably explained by the exponential profiles with the scale lengths of ~120-170A, showing no significant evolution from z~5.7 to z~6.6. We find that 275 LAEs with a large equivalent width (LEW) of >240A are candidates of young-metal poor galaxies and AGNs. We also find that the fraction of LEW LAEs to all ones is 4% and 21% at z~5.7 and z~6.6, respectively. Our LAE catalog includes 11 Lya blobs (LABs) that are LAEs with spatially extended Lya emission whose profile is clearly distinguished from those of stellar objects at the >~ 3sigma level. The number density of the LABs at z=6-7 is ~ 10^-7-10^-6 Mpc^-3, being ~ 10-100 times lower than those claimed for LABs at z~ 2-3, suggestive of disappearing LABs at z>~6, albeit with the different selection methods and criteria for the low and high-z LABs.

We present Lya and UV-nebular emission line properties of bright Lya emitters (LAEs) at z=6-7 with a luminosity of log L_Lya/[erg s-1] = 43-44 identified in the 21-deg2 area of the SILVERRUSH early sample developed with the Subaru Hyper Suprime-Cam (HSC) survey data. Our optical spectroscopy newly confirm 21 bright LAEs with clear Lya emission, and contribute to make a spectroscopic sample of 96 LAEs at z=6-7 in SILVERRUSH. From the spectroscopic sample, we select 7 remarkable LAEs as bright as Himiko and CR7 objects, and perform deep Keck/MOSFIRE and Subaru/nuMOIRCS near-infrared spectroscopy reaching the 3sigma-flux limit of ~ 2x10^{-18} erg s-1 for the UV-nebular emission lines of He II1640, C IV1548,1550, and O III]1661,1666. Except for one tentative detection of C IV, we find no strong UV-nebular lines down to the flux limit, placing the upper limits of the rest-frame equivalent widths (EW_0) of ~2-4 A for He II, C IV, and O III] lines. Here we also investigate the VLT/X-SHOOTER spectrum of CR7 whose 6 sigma detection of He II is claimed by Sobral et al. Although two individuals and the ESO-archive service carefully re-analyze the X-SHOOTER data that are used in the study of Sobral et al., no He II signal of CR7 is detected, supportive of weak UV-nebular lines of the bright LAEs even for CR7. Spectral properties of these bright LAEs are thus clearly different from those of faint dropouts at z~7 that have strong UV-nebular lines shown in the various studies. Comparing these bright LAEs and the faint dropouts, we find anti-correlations between the UV-nebular line EW_0 and UV-continuum luminosity, which are similar to those found at z~2-3.

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.

The Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) is a three-layered imaging survey aimed at addressing some of the most outstanding questions in astronomy today, including the nature of dark matter and dark energy. The survey has been awarded 300 nights of observing time at the Subaru Telescope and it started in March 2014. This paper presents the first public data release of HSC-SSP. This release includes data taken in the first 1.7 years of observations (61.5 nights) and each of the Wide, Deep, and UltraDeep layers covers about 108, 26, and 4 square degrees down to depths of i~26.4, ~26.5, and ~27.0 mag, respectively (5sigma for point sources). All the layers are observed in five broad bands (grizy), and the Deep and UltraDeep layers are observed in narrow bands as well. We achieve an impressive image quality of 0.6 arcsec in the i-band in the Wide layer. We show that we achieve 1-2 per cent PSF photometry (rms) both internally and externally (against Pan-STARRS1), and ~10 mas and 40 mas internal and external astrometric accuracy, respectively. Both the calibrated images and catalogs are made available to the community through dedicated user interfaces and database servers. In addition to the pipeline products, we also provide value-added products such as photometric redshifts and a collection of public spectroscopic redshifts. Detailed descriptions of all the data can be found online. The data release website is https://hsc-release.mtk.nao.ac.jp/.

We conduct a systematic search for galaxies at $z=0.1-1.5$ with [OII]$\lambda3727$, [OIII]$\lambda5007$, or H$\alpha\lambda6563$ emission lines extended over at least 30 kpc by using deep narrowband and broadband imaging in Subaru-XMM Deep Survey (SXDS) field. These extended emission-line galaxies are dubbed [OII], [OIII], or H$\alpha$ blobs. Based on a new selection method that securely selects extended emission-line galaxies, we find 77 blobs at $z=0.40-1.46$ with the isophotal area of emission lines down to $1.2\times10^{-18}$ erg s$^{-1}$ cm$^{-2}$ kpc$^{-2}$. Four of them are spectroscopically confirmed to be [OIII] blobs at $z=0.83$. We identify AGN activities in 8 blobs with X-ray and radio data and find that the fraction of AGN contribution increases with increasing isophotal area of the extended emission. With the Kolmogorov-Smirnov (KS) and Anderson-Darling tests, we confirm that the stellar-mass distributions of H$\alpha$ and [OII] blobs are not drawn from those of the emitters at the $>90$% confidence level in that H$\alpha$ and [OII] blobs are located at the massive end of the distributions, but cannot reject null hypothesis of being the same distributions in terms of the specific star formation rates. It is suggested that galactic-scale outflows tend to be more prominent in more massive star-forming galaxies. Exploiting our sample homogeneously selected over the large area, we derive the number densities of blobs at each epoch. The number densities of blobs decrease drastically with redshifts at the rate that is larger than that of the decrease of cosmic star formation densities.

We present statistical properties of diffuse Lyman-alpha halos (LAHs) around high-$z$ star-forming galaxies with large Subaru samples of Lyman-alpha emitters (LAEs) at $z=2.2$. We make subsamples defined by the physical quantities of LAEs' central Lyman-alpha luminosities, UV magnitudes, Lyman-alpha equivalent widths, and UV slopes, and investigate LAHs' radial surface brightness (SB) profiles and scale lengths $r_n$ as a function of these physical quantities. We find that there exist prominent LAHs around LAEs with faint Lyman-alpha luminosities, bright UV luminosities, and small Lyman-alpha equivalent widths in cumulative radial Lyman-alpha SB profiles. We confirm this trend with the anti-correlation between $r_n$ and Lyman-alpha luminosities (equivalent widths) based on the Spearman's rank correlation coefficient that is $\rho=-0.9$ ($-0.7$) corresponding to the $96\%$ ($93\%$) confidence level, although the correlation between $r_n$ and UV magnitudes is not clearly found in the rank correlation coefficient. Our results suggest that LAEs with properties similar to typical Lyman-break galaxies (with faint Lyman-alpha luminosities and small equivalent widths) possess more prominent LAHs. We investigate scenarios for the major physical origins of LAHs with our results, and find that the cold stream scenario is not preferred, due to the relatively small equivalent widths up to $77$\AA\ in LAHs that include LAEs' central components. There remain two possible scenarios of Lyman-alpha scattering in circum-galactic medium and satellite galaxies that cannot be tested with our observational data.

We present new results on the evolution of the cosmic star formation rate as a function of stellar mass in the SXDS-UDS field. We make use of narrow-band selected emission line galaxies in four redshift slices between z = 1.46 and z = 0.63, and compute stellar masses by fitting a series of templates to recreate each galaxy's star formation history. We determine mass-binned luminosity functions in each redshift slice, and derive the star formation rate density (rhoSFR) as a function of mass using the [OIII] or [OII] emission lines. We calculate dust extinction and metallicity as a function of stellar mass, and investigate the effect of these corrections on the shape of the overall rhoSFR(M). We find that both these corrections are crucial for determining the shape of the rhoSFR(M), and its evolution with redshift. The fully corrected rhoSFR(M) is a relatively flat distribution, with the normalisation moving towards lower values of rhoSFR with increasing cosmic time/decreasing redshift, and requiring star formation to be truncated across all masses studied here. The peak of rhoSFR(M) is found in the 10^10.5<Msun<10^11.0 mass bin at z = 1.46. In the lower redshift slices the location of the peak is less certain, however low mass galaxies in the range 10^7.0<Msun<10^8.0 play an important part in the overall rhoSFR(M) out to at least z ~ 1.2.

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 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 present the ultra-deep Subaru narrowband imaging survey for Lya emitters (LAEs) at $z=7.3$ in SXDS and COSMOS fields with a total integration time of 106 hours. Exploiting our new sharp bandwidth filter, NB101, installed on Suprime-Cam, we have reached $L(Lya)=2.4\times10^{42} \ erg \ s^{-1}$ ($5\sigma$) for $z=7.3$ LAEs, about 4 times deeper than previous Subaru $z \gtrsim 7$ studies, which allows us to reliably investigate evolution of Lya luminosity function (LF), for the first time, down to the luminosity limit same as those of Subaru $z=3.1-6.6$ LAE samples. Surprisingly, we only find three and four LAEs in SXDS and COSMOS fields, respectively, while one expects a total of $\sim 65$ LAEs by our survey in the case of no Lya LF evolution from $z=6.6$ to $7.3$.We identify a decrease of Lya LF from $z=6.6$ to $7.3$ at the $>90\%$ confidence level from our $z=7.3$ Lya LF.Moreover, the evolution of Lya LF is clearly accelerated at $z>6.6$ beyond the measurement uncertainties including cosmic variance. Because no such accelerated evolution of UV-continuum LF or cosmic star-formation rate (SFR) is found at $z\sim 7$, but suggested only at $z>8$ (Oesch et al. 2013, Bouwens et al. 2014), this accelerated Lya LF evolution is explained by physical mechanisms different from a pure SFR decrease but related to Lya production and escape in the process of cosmic reionization. Because a simple accelerating increase of IGM neutral hydrogen absorbing Lya would not reconcile with Thomson scattering optical depth measurements from WMAP and Planck, our findings may support new physical pictures suggested by recent theoretical studies, such as the existence of HI clumpy clouds within cosmic ionized bubbles selectively absorbing Lya and the large ionizing photon escape fraction of galaxies making weak Lya emission.

We report our Keck/MOSFIRE and Magellan/LDSS3 spectroscopy for an [OII] Blob, OIIB10, that is a high-$z$ galaxy with spatially extended [OII]$\lambda\lambda3726,3729$ emission over 30 kpc recently identified by a Subaru large-area narrowband survey. The systemic redshift of OIIB10 is $z=1.18$ securely determined with [OIII]$\lambda\lambda4959,5007$ and H$\beta$ emission lines. We identify FeII$\lambda$2587 and MgII$\lambda\lambda$2796,2804 absorption lines blueshifted from the systemic redshift by $80\pm50$ and $260\pm40$ km s$^{-1}$, respectively, which indicate gas outflow from OIIB10 with the velocity of $\sim 80-260$ km s$^{-1}$. This outflow velocity is comparable with the escape velocity, $250\pm140$ km s$^{-1}$, estimated under the assumption of a singular isothermal halo potential profile. Some fraction of the outflowing gas could escape from the halo of OIIB10, suppressing OIIB10's star-formation activity. We estimate a mass loading factor, $\eta$, that is a ratio of mass outflow rate to star-formation rate, and obtain $\eta>0.8\pm 0.1$ which is relatively high compared with low-$z$ starbursts including U/LIRGs and AGNs. The major energy source of the outflow is unclear with the available data. Although no signature of AGN is found in the X-ray data, OIIB10 falls in the AGN/star-forming composite region in the line diagnostic diagrams. It is possible that the outflow is powered by star formation and a type-2 AGN with narrow FWHM emission line widths of $70-130$ km s$^{-1}$. This is the first detailed spectroscopic study of oxygen-line blobs, which includes the analyses of the escape velocity, the mass loading factor, and the presence of an AGN, and a significant step to understanding the nature of oxygen-line blobs and the relation with gas outflow and star-formation quenching at high redshift.

We present diffuse Lyman-alpha halos (LAHs) identified in the composite Subaru narrowband images of 100-3600 Lyman-alpha emitters (LAEs) at z=2.2, 3.1, 3.7, 5.7, and 6.6. First, we carefully examine potential artifacts mimicking LAHs that include a large-scale point-spread function (PSF) made by instrumental and atmospheric effects. Based on our critical test with composite images of non-LAE samples whose narrowband-magnitude and source-size distributions are the same as our LAE samples, we confirm that no artifacts can produce a diffuse extended feature similar to our LAHs. After this test, we measure the scale lengths of exponential profile for the LAHs estimated from our z=2.2-6.6 LAE samples of L(Lyman-alpha) > 2 x 10^42 erg s^-1. We obtain the scale lengths of ~ 5-10 kpc at z=2.2-5.7, and find no evolution of scale lengths in this redshift range beyond our measurement uncertainties. Combining this result and the previously-known UV-continuum size evolution, we infer that the ratio of LAH to UV-continuum sizes is nearly constant at z=2.2-5.7. The scale length of our z=6.6 LAH is larger than 5-10 kpc just beyond the error bar, which is a hint that the scale lengths of LAHs would increase from z=5.7 to 6.6. If this increase is confirmed by future large surveys with significant improvements of statistical and systematical errors, this scale length change at z > 6 would be a signature of increasing fraction of neutral hydrogen scattering Lyman-alpha photons, due to cosmic reionization.

We present the results of structure analyses for a large sample of 426 Lya emitters (LAEs) at z~2.2 that are observed with HST/ACS and WFC3-IR by deep extra-galactic legacy surveys. We confirm that the merger fraction and the average ellipticity of LAE's stellar component are 10-30 % and 0.4-0.6, respectively, that are comparable with previous study results. We successfully identify that some LAEs have a spatial offset between Lya and stellar-continuum emission peaks, d_Lya, by ~2.5-4 kpc beyond our statistical errors. To uncover the physical origin of strong Lya emission found in LAEs, we investigate Lya equivalent width (EW) dependences of these three structural parameters, merger fraction, d_Lya, and ellipticity of stellar distribution in the range of EW(Lya)=20-250A. Contrary to expectations, we find that merger fraction does not significantly increase with Lya EW. We reveal an anti-correlation between d_Lya and EW(Lya) by Kolmogorov-Smirnov (KS) test. There is a trend that the LAEs with a large Lya EW have a small ellipticity. This is consistent with the recent theoretical claims that Lya photons can more easily escape from face-on disks having a small ellipticity, due to less inter-stellar gas along the line of sight, although our KS test indicates that this trend is not statistically significant. Our results of Lya-EW dependence generally support the idea that an HI column density is a key quantity determining Lya emissivity.

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.

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

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.

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.