We report on the detection of the [CII] 157.7 $\mu$m emission from the Lyman break galaxy (LBG) MACS0416_Y1 at z = 8.3113, by using the Atacama Large Millimeter/submillimeter Array (ALMA). The luminosity ratio of [OIII] 88 $\mu$m (from previous campaigns) to [CII] is 9.31 $\pm$ 2.6, indicative of hard interstellar radiation fields and/or a low covering fraction of photo-dissociation regions. The emission of [CII] is cospatial to the 850 $\mu$m dust emission (90 $\mu$m rest-frame, from previous campaigns), however the peak [CII] emission does not agree with the peak [OIII] emission, suggesting that the lines originate from different conditions in the interstellar medium. We fail to detect continuum emission at 1.5 mm (160 $\mu$m rest-frame) down to 18 $\mu$Jy (3$\sigma$). This nondetection places a strong limit on the dust spectrum, considering the 137 $\pm$ 26 $\mu$Jy continuum emission at 850 $\mu$m. This suggests an unusually warm dust component (T $>$ 80 K, 90% confidence limit), and/or a steep dust-emissivity index ($\beta_{\rm dust}$ $>$ 2), compared to galaxy-wide dust emission found at lower redshifts (typically T $\sim$ 30 - 50 K, $\beta_{\rm dust}$ $\sim$ 1 - 2). If such temperatures are common, this would reduce the required dust mass and relax the dust production problem at the highest redshifts. We therefore warn against the use of only single-wavelength information to derive physical properties, recommend a more thorough examination of dust temperatures in the early Universe, and stress the need for instrumentation that probes the peak of warm dust in the Epoch of Reionization.

We present near-infrared observations of 42 gravitationally lensed galaxies obtained in the framework of the KMOS Lensed Emission Lines and VElocity Review (KLEVER) Survey, a program aimed at investigating the spatially resolved properties of the ionised gas in 1.2<z<2.5 galaxies by means of a full coverage of the YJ, H and K near-infrared bands. Detailed metallicity maps and gradients are derived for a sub-sample of 28 galaxies from reconstructed source plane emission line maps, exploiting the variety of different emission line diagnostics provided by the broad wavelength coverage of the survey. About 85% of these galaxies are characterised by metallicity gradients shallower than 0.05 dex/kpc and 89% are consistent with flat slope within 3$\sigma$ (67% within 1$\sigma$), suggesting a mild evolution with cosmic time. In the context of cosmological simulations and chemical evolution models, the presence of efficient feedback mechanisms and/or extended star formation profiles on top of the classical "inside-out" scenario of mass assembly is generally required to reproduce the observed flatness of the metallicity gradients beyond z$\sim$1 . Three galaxies with significantly (> 3$\sigma$) "inverted" gradients are also found, showing an anti-correlation between metallicity and star formation rate density on local scales, possibly suggesting recent episodes of pristine gas accretion or strong radial flows in place. Nevertheless, the individual metallicity maps are characterised by a variety of different morphologies, with flat radial gradients sometimes hiding non-axisymmetric variations on kpc scales which are washed out by azimuthal averages, especially in interacting systems or in those undergoing local episodes of recent star formation.

We have carried out follow-up spectroscopy on three overdense regions of $g$- and $r$-dropout galaxies in the Canada-France-Hawaii Telescope Legacy Survey Deep Fields, finding two new protoclusters at $z=4.898$, 3.721 and a possible protocluster at $z=3.834$. The $z=3.721$ protocluster overlaps with a previously identified protocluster at $z=3.675$. The redshift separation between these two protoclusters is $\Delta z=0.05$, which is slightly larger than the size of typical protoclusters. Therefore, if they are not the progenitors of a $>10^{15}\,\mathrm{M_\odot}$ halo, they would grow into closely-located independent halos like a supercluster. The other protocluster at $z=4.898$ is also surrounded by smaller galaxy groups. These systems including protoclusters and neighboring groups are regarded as the early phase of superclusters. We quantify the spatial distribution of member galaxies of the protoclusters at $z=3.675$ and 3.721 by fitting triaxial ellipsoids, finding a tentative difference: one has a pancake-like shape while the other is filamentary. This could indicate that these two protoclusters are in different stages of formation. We investigate the relation between redshift and the velocity dispersion of protoclusters, including other protoclusters from the literature, in order to compare their dynamical states. Although there is no significant systematic trend in the velocity dispersions of protoclusters with redshift, the distribution is skewed to higher velocity dispersion over the redshift range of $z=2\mathrm{-}6$. This could be interpreted as two phases of cluster formation, one dominated by the steady accretion of galaxies, and the other by the merging between group-size halos, perhaps depending on the surrounding large-scale environments.

We present new ALMA observations and physical properties of a Lyman Break Galaxy at z=7.15. Our target, B14-65666, has a bright ultra-violet (UV) absolute magnitude, $M_{\rm UV}\approx-22.4$, and has been spectroscopically identified in Ly$\alpha$ with a small rest-frame equivalent width of $\approx4$ \AA. Previous HST image has shown that the target is comprised of two spatially separated clumps in the rest-frame UV. With ALMA, we have newly detected spatially resolved [OIII] 88 $\mu$m, [CII] 158 $\mu$m, and their underlying dust continuum emission. In the whole system of B14-65666, the [OIII] and [CII] lines have consistent redshifts of $7.1520\pm0.0003$, and the [OIII] luminosity, $(34.4\pm4.1)\times10^{8}L_{\rm \odot}$, is about three times higher than the [CII] luminosity, $(11.0\pm1.4)\times10^{8}L_{\rm \odot}$. With our two continuum flux densities, the dust temperature is constrained to be $T_{\rm d}\approx50-60$ K under the assumption of the dust emissivity index of $\beta_{\rm d}=2.0-1.5$, leading to a large total infrared luminosity of $L_{\rm TIR}\approx1\times10^{12}L_{\rm \odot}$. Owing to our high spatial resolution data, we show that the [OIII] and [CII] emission can be spatially decomposed into two clumps associated with the two rest-frame UV clumps whose spectra are kinematically separated by $\approx200$ km s$^{-1}$. We also find these two clumps have comparable UV, infrared, [OIII], and [CII] luminosities. Based on these results, we argue that B14-65666 is a starburst galaxy induced by a major-merger. The merger interpretation is also supported by the large specific star-formation rate (defined as the star-formation rate per unit stellar mass), sSFR$=260^{+119}_{-57}$ Gyr$^{-1}$, inferred from our SED fitting. Probably, a strong UV radiation field caused by intense star formation contributes to its high dust temperature and the [OIII]-to-[CII] luminosity ratio.

We present the Atacama Large Millimeter/submillimeter Array (ALMA) detection of the [O III] 88 $\mu$m line and rest-frame 90 $\mu$m dust continuum emission in a Y-dropout Lyman break galaxy (LBG), MACS0416_Y1, lying behind the Frontier Field cluster MACS J0416.1-2403. This [O III] detection confirms the LBG with a spectroscopic redshift of $z = 8.3118 \pm 0.0003$, making this object one of the furthest galaxies ever identified spectroscopically. The observed 850 $\mu$m flux density of $137 \pm 26$ $\mu$Jy corresponds to a de-lensed total infrared (IR) luminosity of $L_{\rm IR} = (1.7 \pm 0.3) \times 10^{11} L_{\odot}$ if assuming a dust temperature of $T_{\rm dust} = 50$ K and an emissivity index of $\beta = 1.5$, yielding a large dust mass of $4 \times 10^6 M_{\odot}$. The ultraviolet-to-far IR spectral energy distribution modeling where the [O III] emissivity model is incorporated suggests the presence of a young ($\tau_{\rm age} \approx 4$ Myr), star-forming (SFR $\approx 60 M_{\odot}$ yr$^{-1}$), moderately metal-polluted ($Z \approx 0.2 Z_{\odot}$) stellar component with a mass of $M_{\rm star} = 3 \times 10^8 M_{\odot}$. An analytic dust mass evolution model with a single episode of star-formation does not reproduce the metallicity and dust mass in $\tau_{\rm age} \approx 4$ Myr, suggesting a pre-existing evolved stellar component with $M_{\rm star} \sim 3 \times 10^9 M_{\odot}$ and $\tau_{\rm age} \sim 0.3$ Gyr as the origin of the dust mass.

Quasars (QSOs) hosting supermassive black holes are believed to reside in massive halos harboring galaxy overdensities. However, many observations revealed average or low galaxy densities around $z\gtrsim6$ QSOs. This could be partly because they measured galaxy densities in only tens of arcmin$^2$ around QSOs and might have overlooked potential larger scale galaxy overdensities. Some previous studies also observed only Lyman break galaxies (LBGs, massive older galaxies) and missed low mass young galaxies like Ly$\alpha$ emitters (LAEs) around QSOs. Here we present observations of LAE and LBG candidates in $\sim700$ arcmin$^2$ around a $z=6.61$ luminous QSO using Subaru Telescope Suprime-Cam with narrow/broadbands. We compare their sky distributions, number densities and angular correlation functions with those of LAEs/LBGs detected in the same manner and comparable data quality in our control blank field. In the QSO field, LAEs and LBGs are clustering in 4-20 comoving Mpc angular scales, but LAEs show mostly underdensity over the field while LBGs are forming $30\times60$ comoving Mpc$^2$ large scale structure containing 3-$7\sigma$ high density clumps. The highest density clump includes a bright (23.78 mag in the narrowband) extended ($\gtrsim 16$ kpc) Ly$\alpha$ blob candidate, indicative of a dense environment. The QSO could be part of the structure but is not located exactly at any of the high density peaks. Near the QSO, LAEs show underdensity while LBGs average to $4\sigma$ excess densities compared to the control field. If these environments reflect halo mass, the QSO may not be in the most massive halo, but still in a moderately massive one.

We detect 20 $z=7.0$ Ly$\alpha$ emitter (LAE) candidates to $L({\rm Ly}\alpha) \geq 2 \times 10^{42}$ erg s$^{-1}$ or $0.3$ $L^*_{z=7}$ and in $6.1\times 10^5$ Mpc$^3$ volume in the Subaru Deep Field and the Subaru/XMM-Newton Deep Survey field by 82 and 37 hours of Subaru Suprime-Cam narrowband NB973 and reddest optical $y$-band imaging. We compare their Ly$\alpha$ and UV luminosity functions (LFs) and densities and Ly$\alpha$ equivalent widths (EWs) to those of $z=5.7$, 6.6 and 7.3 LAEs from previous Suprime-Cam surveys. The Ly$\alpha$ LF (density) rapidly declines by a factor of $\times$1.5 (1.9) in $L({\rm Ly}\alpha)$ at $z=5.7-6.6$ (160 Myr), $\times$1.5 (1.6) at $z=6.6-7.0$ (60 Myr) at the faint end and $\times$2.0 (3.8) at $z=7.0-7.3$ (40 Myr). Also, in addition to the systematic decrease in EW at $z=5.7-6.6$ previously found, 2/3 of the $z=7.0$ LAEs detected in the UV continuum exhibit lower EWs than the $z=6.6$ ones. Moreover, while the UV LF and density do not evolve at $z=5.7-6.6$, they modestly decline at $z=6.6-7.0$, implying galaxy evolution contributing to the decline of the Ly$\alpha$ LF. Comparison of the $z=7.0$ Ly$\alpha$ LF to the one predicted by an LAE evolution model further reveals that galaxy evolution alone cannot explain all the decline of Ly$\alpha$ LF. If we attribute the discrepancy to Ly$\alpha$ attenuation by neutral hydrogen, the intergalactic medium transmission of Ly$\alpha$ photons at $z=7.0$ would be $T_{{\rm Ly}\alpha}^{\rm IGM} \leq 0.6-0.7$. It is lower (higher) than $T_{{\rm Ly}\alpha}^{\rm IGM}$ at $z=6.6$ (7.3) derived by previous studies, suggesting rapid increase in neutral fraction at $z > 6$.

In this paper we use ASPECS, the ALMA Spectroscopic Survey in the {\em Hubble} Ultra Deep Field (UDF) in band 3 and band 6, to place blind constraints on the CO luminosity function and the evolution of the cosmic molecular gas density as a function of redshift up to $z\sim 4.5$. This study is based on galaxies that have been solely selected through their CO emission and not through any other property. In all of the redshift bins the ASPECS measurements reach the predicted `knee' of the CO luminosity function (around $5\times10^{9}$ K km/s pc$^2$). We find clear evidence of an evolution in the CO luminosity function with respect to $z\sim 0$, with more CO luminous galaxies present at $z\sim 2$. The observed galaxies at $z\sim 2$ also appear more gas-rich than predicted by recent semi-analytical models. The comoving cosmic molecular gas density within galaxies as a function of redshift shows a factor 3-10 drop from $z \sim 2$ to $z \sim 0$ (with significant error bars), and possibly a decline at $z>3$. This trend is similar to the observed evolution of the cosmic star formation rate density. The latter therefore appears to be at least partly driven by the increased availability of molecular gas reservoirs at the peak of cosmic star formation ($z\sim2$).

We study the molecular gas properties of high-$z$ galaxies observed in the ALMA Spectroscopic Survey (ASPECS) that targets a $\sim1$ arcmin$^2$ region in the Hubble Ultra Deep Field (UDF), a blind survey of CO emission (tracing molecular gas) in the 3mm and 1mm bands. Of a total of 1302 galaxies in the field, 56 have spectroscopic redshifts and correspondingly well-defined physical properties. Among these, 11 have infrared luminosities $L_{\rm{}IR}>10^{11}$ L$_\odot$, i.e. a detection in CO emission was expected. Out these, 7 are detected at various significance in CO, and 4 are undetected in CO emission. In the CO-detected sources, we find CO excitation conditions that are lower than typically found in starburst/SMG/QSO environments. We use the CO luminosities (including limits for non-detections) to derive molecular gas masses. We discuss our findings in context of previous molecular gas observations at high redshift (star-formation law, gas depletion times, gas fractions): The CO-detected galaxies in the UDF tend to reside on the low-$L_{\rm{}IR}$ envelope of the scatter in the $L_{\rm{}IR}-L'_{\rm{}CO}$ relation, but exceptions exist. For the CO-detected sources, we find an average depletion time of $\sim$ 1 Gyr, with significant scatter. The average molecular-to-stellar mass ratio ($M_{\rm{}H2}$/$M_*$) is consistent with earlier measurements of main sequence galaxies at these redshifts, and again shows large variations among sources. In some cases, we also measure dust continuum emission. On average, the dust-based estimates of the molecular gas are a factor $\sim$2-5$\times$ smaller than those based on CO. Accounting for detections as well as non-detections, we find large diversity in the molecular gas properties of the high-redshift galaxies covered by ASPECS.

We present an analysis of a deep (1$\sigma$=13 $\mu$Jy) cosmological 1.2-mm continuum map based on ASPECS, the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field. In the 1 arcmin$^2$ covered by ASPECS we detect nine sources at $>3.5\sigma$ significance at 1.2-mm. Our ALMA--selected sample has a median redshift of $z=1.6\pm0.4$, with only one galaxy detected at z$>$2 within the survey area. This value is significantly lower than that found in millimeter samples selected at a higher flux density cut-off and similar frequencies. Most galaxies have specific star formation rates similar to that of main sequence galaxies at the same epoch, and we find median values of stellar mass and star formation rates of $4.0\times10^{10}\ M_\odot$ and $\sim40~M_\odot$ yr$^{-1}$, respectively. Using the dust emission as a tracer for the ISM mass, we derive depletion times that are typically longer than 300 Myr, and we find molecular gas fractions ranging from $\sim$0.1 to 1.0. As noted by previous studies, these values are lower than using CO--based ISM estimates by a factor $\sim$2. The 1\,mm number counts (corrected for fidelity and completeness) are in agreement with previous studies that were typically restricted to brighter sources. With our individual detections only, we recover $55\pm4\%$ of the extragalactic background light (EBL) at 1.2 mm measured by the Planck satellite, and we recover $80\pm7\%$ of this EBL if we include the bright end of the number counts and additional detections from stacking. The stacked contribution is dominated by galaxies at $z\sim1-2$, with stellar masses of (1-3)$\times$10$^{10}$ M$_\odot$. For the first time, we are able to characterize the population of galaxies that dominate the EBL at 1.2 mm.

We present Karl G Jansky Very Large Array (VLA) observations of CO(2-1) line emission and rest-frame 250GHz continuum emission of the Hyper-Luminous IR Galaxies (HyLIRGs) BRI1202-0725 (z=4.69) and BRI1335-0417 (z=4.41), with an angular resolution as high as 0.15". Our low order CO observations delineate the cool molecular gas, the fuel for star formation in the systems, in unprecedented detail. For BRI1202-0725, line emission is seen from both extreme starburst galaxies: the quasar host and the optically obscured submm galaxy (SMG), in addition to one of the Lyman-alpha emitting galaxies in the group. Line emission from the SMG shows an east-west extension of about 0.6". For Lyalpha-2, the CO emission is detected at the same velocity as [CII] and [NII], indicating a total gas mass ~4.0*10^10 solar masses. The CO emission from BRI1335-0417 peaks at the nominal quasar position, with a prominent northern extension (~1", a possible tidal feature). The gas depletion timescales are ~10^7 years for the three HyLIRGs, consistent with extreme starbursts, while that of Lyalpha-2 may be consistent with main sequence galaxies. We interpret these sources as major star formation episodes in the formation of massive galaxies and supermassive black holes (SMBHs) via gas rich mergers in the early Universe.

We present a search for [CII] line and dust continuum emission from optical dropout galaxies at $z>6$ using ASPECS, our ALMA Spectroscopic Survey in the Hubble Ultra-Deep Field (UDF). Our observations, which cover the frequency range $212-272$ GHz, encompass approximately the range $6<z<8$ for [CII] line emission and reach a limiting luminosity of L$_{\rm [CII]}\sim$(1.6-2.5)$\times$10$^{8}$ L$_{\odot}$. We identify fourteen [CII] line emitting candidates in this redshift range with significances $>$4.5 $\sigma$, two of which correspond to blind detections with no optical counterparts. At this significance level, our statistical analysis shows that about 60\% of our candidates are expected to be spurious. For one of our blindly selected [CII] line candidates, we tentatively detect the CO(6-5) line in our parallel 3-mm line scan. None of the line candidates are individually detected in the 1.2 mm continuum. A stack of all [CII] candidates results in a tentative detection with $S_{1.2mm}=14\pm5\mu$Jy. This implies a dust-obscured star formation rate (SFR) of $(3\pm1)$ M$_\odot$ yr$^{-1}$. We find that the two highest--SFR objects have candidate [CII] lines with luminosities that are consistent with the low-redshift $L_{\rm [CII]}$ vs. SFR relation. The other candidates have significantly higher [CII] luminosities than expected from their UV--based SFR. At the current sensitivity it is unclear whether the majority of these sources are intrinsically bright [CII] emitters, or spurious sources. If only one of our line candidates was real (a scenario greatly favored by our statistical analysis), we find a source density for [CII] emitters at $6<z<8$ that is significantly higher than predicted by current models and some extrapolations from galaxies in the local universe.

We make use of deep 1.2mm-continuum observations (12.7microJy/beam RMS) of a 1 arcmin^2 region in the Hubble Ultra Deep Field to probe dust-enshrouded star formation from 330 Lyman-break galaxies spanning the redshift range z=2-10 (to ~2-3 Msol/yr at 1sigma over the entire range). Given the depth and area of ASPECS, we would expect to tentatively detect 35 galaxies extrapolating the Meurer z~0 IRX-beta relation to z>~2 (assuming T_d~35 K). However, only 6 tentative detections are found at z>~2 in ASPECS, with just three at >3sigma. Subdividing z=2-10 galaxies according to stellar mass, UV luminosity, and UV-continuum slope and stacking the results, we only find a significant detection in the most massive (>10^9.75 Msol) subsample, with an infrared excess (IRX=L_{IR}/L_{UV}) consistent with previous z~2 results. However, the infrared excess we measure from our large selection of sub-L* (<10^9.75 Msol) galaxies is 0.11(-0.42)(+0.32) and 0.14(-0.14)(+0.15) at z=2-3 and z=4-10, respectively, lying below even an SMC IRX-beta relation (95% confidence). These results demonstrate the relevance of stellar mass for predicting the IR luminosity of z>~2 galaxies. We furthermore find that the evolution of the IRX-stellar mass relationship depends on the evolution of the dust temperature. If the dust temperature increases monotonically with redshift (as (1+z)^0.32) such that T_d~44-50 K at z>=4, current results are suggestive of little evolution in this relationship to z~6. We use these results to revisit recent estimates of the z>~3 SFR density. One less obvious implication is in interpreting the high Halpha EWs seen in z~5 galaxies: our results imply that star-forming galaxies produce Lyman-continuum photons at twice the efficiency (per unit UV luminosity) as implied in conventional models. Star-forming galaxies can then reionize the Universe, even if the escape fraction is <10%.

We present the rationale for and the observational description of ASPECS: The ALMA SPECtroscopic Survey in the Hubble Ultra-Deep Field (UDF), the cosmological deep field that has the deepest multi-wavelength data available. Our overarching goal is to obtain an unbiased census of molecular gas and dust continuum emission in high-redshift (z$>$0.5) galaxies. The $\sim$1$'$ region covered within the UDF was chosen to overlap with the deepest available imaging from HST. Our ALMA observations consist of full frequency scans in band 3 (84-115 GHz) and band 6 (212-272 GHz) at approximately uniform line sensitivity ($L'_{\rm CO}\sim$2$\times$10$^{9}$ K km/s pc$^2$), and continuum noise levels of 3.8 $\mu$Jy beam$^{-1}$ and 12.7 $\mu$Jy beam$^{-1}$, respectively. The molecular surveys cover the different rotational transitions of the CO molecule, leading to essentially full redshift coverage. The [CII] emission line is also covered at redshifts $6.0<z<8.0$. We present a customized algorithm to identify line candidates in the molecular line scans, and quantify our ability to recover artificial sources from our data. Based on whether multiple CO lines are detected, and whether optical spectroscopic redshifts as well as optical counterparts exist, we constrain the most likely line identification. We report 10 (11) CO line candidates in the 3mm (1mm) band, and our statistical analysis shows that $<$4 of these (in each band) are likely spurious. Less than 1/3 of the total CO flux in the low-J CO line candidates are from sources that are not associated with an optical/NIR counterpart. We also present continuum maps of both the band 3 and band 6 observations. The data presented here form the basis of a number of dedicated studies that are presented in subsequent papers.

The physical properties and elemental abundances of the interstellar medium in galaxies during cosmic reionization are important for understanding the role of galaxies in this process. We report the Atacama Large Millimeter/submillimeter Array detection of an oxygen emission line at a wavelength of 88 micrometers from a galaxy at an epoch about 700 million years after the Big Bang. The oxygen abundance of this galaxy is estimated at about one-tenth that of the Sun. The non-detection of far-infrared continuum emission indicates a deficiency of interstellar dust in the galaxy. A carbon emission line at a wavelength of 158 micrometers is also not detected, implying an unusually small amount of neutral gas. These properties might allow ionizing photons to escape into the intergalactic medium.

We present the discovery of three protoclusters at $z\sim3\mathrm{-}4$ with spectroscopic confirmation in the Canada-France-Hawaii Telescope (CFHT) Legacy Survey Deep Fields. In these fields, we investigate the large-scale projected sky distribution of $z\sim3\mathrm{-}6$ Lyman break galaxies and identify 21 protocluster candidates from regions that are overdense at more than $4\sigma$ overdensity significance. Based on cosmological simulations, it is expected that more than $76\%$ of these candidates will evolve into a galaxy cluster of at least a halo mass of $10^{14}\,\mathrm{M_\odot}$ at $z=0$. We perform follow-up spectroscopy for eight of the candidates using Subaru/FOCAS, KeckII/DEIMOS, and Gemini-N/GMOS. In total we target 462 dropout candidates and obtain 138 spectroscopic redshifts. We confirm three real protoclusters at $z=3\mathrm{-}4$ with more than five members spectroscopically identified, and find one to be an incidental overdense region by mere chance alignment. The other four candidate regions at $z\sim5\mathrm{-}6$ require more spectroscopic follow-up in order to be conclusive. A $z=3.67$ protocluster, which has eleven spectroscopically confirmed members, shows a remarkable core-like structure composed of a central small region ($<0.5\,\mathrm{physical\>Mpc}$) and an outskirts region ($\sim1.0\,\mathrm{physical\>Mpc}$). The Ly$\alpha$ equivalent widths of members of the protocluster are significantly smaller than those of field galaxies at the same redshift while there is no difference in the UV luminosity distributions. These results imply that some environmental effects start operating as early as at $z\sim4$ along with the growth of the protocluster structure.

We have analysed 18 ALMA continuum maps in Bands 6 and 7, with rms down to 7.8$\mu$Jy, to derive differential number counts down to 60$\mu$Jy and 100$\mu$Jy at $\lambda=$1.3 mm and $\lambda=$1.1 mm, respectively. The area covered by the combined fields is $\rm 9.5\times10^{-4}deg^2$ at 1.1mm and $\rm 6.6\times10^{-4}deg^{2}$ at 1.3mm. We improved the source extraction method by requiring that the dimension of the detected sources be consistent with the beam size. This method enabled us to remove spurious detections that have plagued the purity of the catalogues in previous studies. We detected 50 faint sources with S/N$>$3.5 down to 60$\mu$Jy, hence improving the statistics by a factor of four relative to previous studies. The inferred differential number counts are $\rm dN/d(Log_{10}S)=1\times10^5~deg^2$ at a 1.1 mm flux $S_{\lambda = 1.1~mm} = 130~\mu$Jy, and $\rm dN/d(Log_{10}S)=1.1\times10^5~deg^2$ at a 1.3 mm flux $\rm S_{\lambda = 1.3~mm} = 60~\mu$Jy. At the faintest flux limits, i.e. 30$\mu$Jy and 40$\mu$Jy, we obtain upper limits on the differential number counts of $\rm dN/d(Log_{10}S) < 7\times10^5~deg^2$ and $\rm dN/d(Log_{10}S)<3\times10^5~deg^2$, respectively. Our results provide a new lower limit to CIB intensity of 17.2${\rm Jy\ deg^{-2}}$ at 1.1mm and of 12.9${\rm Jy\ deg^{-2}}$ at 1.3mm. Moreover, the flattening of the integrated number counts at faint fluxes strongly suggests that we are probably close to the CIB intensity. Our data imply that galaxies with SFR$<40~M_{\odot}/yr$ certainly contribute less than 50% to the CIB while more than 50% of the CIB must be produced by galaxies with $\rm SFR>40~M_{\odot}/yr$. The differential number counts are in nice agreement with recent semi-analytical models of galaxy formation even as low as our faint fluxes. Consequently, this supports the galaxy evolutionary scenarios and assumptions made in these models.

We present new observations of the highest-redshift quasar known, ULAS J1120+0641, redshift $z=7.084$, obtained in the optical, at near-, mid-, and far-infrared wavelengths, and in the sub-mm. We combine these results with published X-ray and radio observations to create the multiwavelength spectral energy distribution (SED), with the goals of measuring the bolometric luminosity $L_{\rm bol}$, and quantifying the respective contributions from the AGN and star formation. We find three components are needed to fit the data over the wavelength range $0.12-1000\,\mu$m: the unobscured quasar accretion disk and broad-line region, a dusty clumpy AGN torus, and a cool 47K modified black body to characterise star formation. Despite the low signal-to-noise ratio of the new long-wavelength data, the normalisation of any dusty torus model is constrained within $\pm40\%$. We measure a bolometric luminosity $L_{\rm bol}=2.6\pm0.6\times10^{47}\,$erg$\,$s$^{-1}=6.7 \pm 1.6\times10^{13}L_{\odot}$, to which the three components contribute $31\%,32\%,3\%$, respectively, with the remainder provided by the extreme UV $<0.12\,\mu$m. We tabulate the best-fit model SED. We use local scaling relations to estimate a star formation rate (SFR) in the range $60-270\,{\rm M}_\odot$/yr from the [C$\,{\scriptsize \rm II}$] line luminosity and the $158\,\mu$m continuum luminosity. An analysis of the equivalent widths of the [C$\,{\scriptsize \rm II}$] line in a sample of $z>5.7$ quasars suggests that these indicators are promising tools for estimating the SFR in high-redshift quasars in general. At the time observed the black hole was growing in mass more than 100 times faster than the stellar bulge, relative to the mass ratio measured in the local universe, i.e. compared to ${M_{\rm BH}}/{M_{\rm bulge}} \simeq 1.4\times10^{-3}$, for ULAS J1120+0641 we measure ${\dot{M}_{\rm BH}}/{\dot{M}_{\rm bulge}} \simeq 0.2$.

We present ALMA observations of the [CII] line and far-infrared (FIR) continuum of a normally star-forming galaxy in the reionization epoch, the z=6.96 Ly-alpha emitter (LAE) IOK-1. Probing to sensitivities of sigma_line = 240 micro-Jy/beam (40 km/s channel) and sigma_cont = 21 micro-Jy/beam, we found the galaxy undetected in both [CII] and continuum. Comparison of UV - FIR spectral energy distribution (SED) of IOK-1, including our ALMA limit, with those of several types of local galaxies (including the effects of the cosmic microwave background, CMB, on the FIR continuum) suggests that IOK-1 is similar to local dwarf/irregular galaxies in SED shape rather than highly dusty/obscured galaxies. Moreover, our 3 sigma FIR continuum limit, corrected for CMB effects, implies intrinsic dust mass M_dust < 6.4 x 10^7 M_sun, FIR luminosity L_FIR < 3.7 x 10^{10} L_sun (42.5 - 122.5 micron), total IR luminosity L_IR < 5.7 x 10^{10} L_sun (8 - 1000 micron) and dust-obscured star formation rate (SFR) < 10 M_sun/yr, if we assume that IOK-1 has a dust temperature and emissivity index typical of local dwarf galaxies. This SFR is 2.4 times lower than one estimated from the UV continuum, suggesting that < 29% of the star formation is obscured by dust. Meanwhile, our 3 sigma [CII] flux limit translates into [CII] luminosity, L_[CII] < 3.4 x 10^7 L_sun. Locations of IOK-1 and previously observed LAEs on the L_[CII] vs. SFR and L_[CII]/L_FIR vs. L_FIR diagrams imply that LAEs in the reionization epoch have significantly lower gas and dust enrichment than AGN-powered systems and starbursts at similar/lower redshifts, as well as local star-forming galaxies.

We performed a systematic spectroscopic observation of a protocluster at $z=6.01$ in the Subaru Deep Field. We took spectroscopy for all 53 $i'$-dropout galaxies down to $z'=27.09\,\mathrm{mag}$ in/around the protocluster region. From these observations, we confirmed that 28 galaxies are at $z\sim6$, of which ten are clustered in a narrow redshift range of $\Delta z<0.06$. To trace the evolution of this primordial structure, we applied the same $i'$-dropout selection and the same overdensity measurements used in the observations to a semi-analytic model built upon the Millennium Simulation. We obtain a relation between the significance of overdensities observed at $z\sim6$ and the predicted dark matter halo mass at $z=0$. This protocluster with $6\sigma$ overdensity is expected to grow into a galaxy cluster with a mass of $\sim5\times10^{14}\,\mathrm{M_\odot}$ at $z=0$. Ten galaxies within $10\,\mathrm{comoving\>Mpc}$ of the overdense region can, with more than an 80% probability, merge into a single dark matter halo by $z=0$. No significant differences appeared in UV and Ly$\alpha$ luminosities between the protocluster and field galaxies, suggesting that this protocluster is still in the early phase of cluster formation before the onset of any obvious environmental effects. However, further observations are required to study other properties, such as stellar mass, dust, and age. We do find that galaxies tend to be in close pairs in this protocluster. These pair-like subgroups will coalesce into a single halo and grow into a more massive structure. We may witness an onset of cluster formation at $z\sim6$ toward a cluster as seen in local universe.

We have identified a very interesting Ly-alpha emitter, whose Ly-alpha emission line has an extremely large observed equivalent width of EW_0=436^{+422}_{-149}A, which corresponds to an extraordinarily large intrinsic rest-frame equivalent width of EW_0^{int}=872^{+844}_{-298}A after the average intergalactic absorption correction. The object was spectroscopically confirmed to be a real Ly-alpha emitter by its apparent asymmetric Ly-alpha line profile detected at z=6.538. The continuum emission of the object was definitely detected in our deep z'-band image; thus, its EW_0 was reliably determined. Follow-up deep near-infrared spectroscopy revealed emission lines of neither He II lambda1640 as an apparent signature of Population III, nor C IV lambda1549 as a proof of active nucleus. No detection of short-lived He II lambda1640 line is not necessarily inconsistent with the interpretation that the underlying stellar population of the object is dominated by Population III. We found that the observed extremely large EW_0 of the Ly-alpha emission and the upper limit on the EW_0 of the He II lambda1640 emission can be explained by population synthesis models favoring a very young age less than 2-4Myr and massive metal-poor (Z<10^{-5}) or even metal-free stars. The observed large EW_0 of Ly-alpha is hardly explained by Population I/II synthesis models with Z>10^{-3}. However, we cannot conclusively rule out the possibility that this object is composed of a normal stellar population with a clumpy dust distribution, which could enhance the Ly-alpha EW_0, though its significance is still unclear.

We present the first detailed study of the stellar populations of star-forming galaxies at z~1.5, which are selected by their [O II] emission line, detected in narrow-band surveys. We identified ~1,300 [O II] emitters at z=1.47 and z=1.62 in the Subaru Deep Field with rest-frame EWs above 13\AA. Optical and near-infrared spectroscopic observations for ~10% of our samples show that our separation of [O II] from [O III] emission-line galaxies in two-color space is 99% successful. We analyze the multi-wavelength properties of a subset of ~1,200 galaxies with the best photometry. They have average rest-frame EW of 45\AA, stellar mass of 3 x 10^9 M_sun, and stellar age of 100 Myr. In addition, our SED fitting and broad-band colors indicate that [O II] emitters span the full range of galaxy populations at z~1.5. We also find that 80% of [O II] emitters are also photometrically classified as "BX/BM" (UV) galaxies and/or the star-forming "BzK" (near-IR) galaxies. Our [O II] emission line survey produces a far more complete, and somewhat deeper sample of z~1.5 galaxies than either the BX/BM or sBzK selection alone. We constructed average SEDs and find that higher [O II] EW galaxies have somewhat bluer continua. SED model-fitting shows that they have on average half the stellar mass of galaxies with lower [O II] EW. The observed [O II] luminosity is well-correlated with the far-UV continuum with a logarithmic slope slightly 0f 0.89\pm0.22. The scatter of the [O II] luminosity against the far-UV continuum suggests that [O II] can be used as a SFR indicator with a reliability of 0.23 dex.

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 searched for z=7.3 Lya emitters (LAEs) behind two lensing clusters, Abell 2390 and CL 0024, with the Subaru Telescope Suprime-Cam and a narrowband NB1006 (FWHM ~ 21 nm centered at 1005 nm). We investigated if there exist objects consistent with the color of z=7.3 LAEs behind the clusters but could not detect any LAEs to the unlensed line limit F(Lya) ~ 6.9 x 10^{-18} erg/s/cm^2. Using several z=7 Lya luminosity functions (LFs) from the literature, we estimated and compared the expected detection numbers of z ~ 7 LAEs in lensing and blank field surveys in the case of using an 8m class ground based telescope. Given the steep bright-end slope of the LFs, when the detector field-of view (FOV) is comparable to the angular extent of a massive lensing cluster, imaging cluster(s) is more efficient in detecting z ~ 7 LAEs than imaging a blank field. However, the gain is expected to be modest, a factor of two at most and likely much less depending on the adopted LFs. The main advantage of lensing-cluster survey, therefore, remains to be the gain in depth and not necessarily in detection efficiency. For much larger detectors, the lensing effect becomes negligible and the efficiency of LAE detection is proportional to the instrumental FOV. We also inspected NB1006 images of three z ~ 7 z-dropouts previously detected in Abell 2390 and found that none of them are detected in NB1006. Two of them are consistent with predictions from the previous studies that they would be at lower redshifts. The other one has a photometric redshift of z ~ 7.3, and if it is at z=7.3, its unlensed Lya line flux would be very faint: F(Lya) < 4.4 x 10^{-18} erg/s/cm^2 (1 sigma upper limit) or rest frame equivalent width of W(Lya) < 26A. Its Lya emission might be attenuated by neutral hydrogen, as recent studies show that the fraction of Lyman break galaxies displaying strong Lya emission is lower at z ~ 7 than at z <~ 6.

We have performed deep imaging surveys for LyA emitters (LAEs) at redshift ~7.3 in two blank fields, the Subaru Deep Field (SDF) and the Subaru/XMM-Newton Deep survey Field (SXDF), using the Subaru/Suprime-Cam equipped with new red-sensitive CCDs and a new narrow-band filter, NB1006 (lambda_c=10052 Ang, FWHM=214 Ang). We identified four objects as LAE candidates that exhibit luminosity excess in NB1006. By carrying out deep follow-up spectroscopy for three of them using Subaru/FOCAS and Keck/DEIMOS, a definitively asymmetric emission line is detected for one of them, SXDF-NB1006-2. Assuming this line is LyA, this object is a LAE at z=7.215 which has luminosity of 1.2^{+1.5}_{-0.6} x 10^43 [erg s-1] and a weighted skewness S_w=4.90+-0.86. Another object, SDF-NB1006-2, shows variable photometry and is thus probably a quasar (QSO) or an active galactic nucleus (AGN). It shows an asymmetric emission line at 10076 Ang, which may be due to either LyA at z=7.288 or [OII] at z=1.703. The third object, SDF-NB1006-1, is likely a galaxy with temporal luminosity enhancement associated with a supernova explosion, as the brightness of this object varies between the observed epochs. Its spectrum does not show any emission lines. The inferred decrease in the number density of LAEs toward higher redshift is n_LyA(z=7.3)/n_LyA(z=5.7) = 0.05^+0.11_-0.05 from z=5.7 to 7.3 down to L(LyA)=1.0 x 10^43 [erg s-1]. The present result is consistent with the interpretation in previous studies that the neutral hydrogen fraction is rapidly increasing from z=5.7 to 7.3.

Recent observations of z >~ 7 Ly{\alpha} emitters (LAEs) have derived a variety of Ly{\alpha} luminosity functions (LFs) with contradictory results, evolution or non-evolution from z <~ 6, the epoch after reionization. This could be because most of z >~ 7 LFs comprise photometric candidates and might include some contaminations. We conducted the Subaru Telescope Faint Object Camera And Spectrograph narrowband NB980 ({\lambda}c ~ 9800A, FWHM ~ 100A) imaging and spectroscopy survey of z=7-7.1 LAEs to compare its "contamination-free" result with z >~ 7 photometric Ly{\alpha} LFs previously derived. We imaged the Subaru Deep Field and the sky around a cluster MS 1520.1+3002 and found one LAE candidate, but spectroscopy did not reveal Ly{\alpha} though deep enough to detect it. We calculated the expected number of LAEs in our survey, using five z=7 and three z=7.7 Ly{\alpha} LFs from recent surveys. Seven of them are consistent with null detection (0.1^{+1.8}_{-0.1}-1.1^{+2.2}_{-1.0} LAEs) within errors including Poisson statistics and cosmic variance, but average values (0.7-1.1 LAEs) predicted from one z=7 and two z=7.7 LFs among the seven indicate nearly a single detection. The remaining one z=7 LF predicts 3.0^{+3.2}_{-2.0} LAEs. As to z=7, the discrepancy likely comes from different LAE selection criteria. For z=7.7, there are two possibilities; (1) If z=7.7 LAEs are somehow brighter in Ly{\alpha} luminosity than lower redshift LAEs, z=7.7 LF is observed to be similar to or higher than lower redshift LFs even if attenuated by neutral hydrogen. (2) All/most of the z=7.7 candidates are not LAEs. This supports the decline of LF from z ~ 6 to 7.7 and reionization at z ~ 6-7.7.

We present a panoramic narrow-band imaging survey of [OII] emitters in and around the ClG J0218.3-0510 cluster at z=1.62 with Suprime-Cam on Subaru telescope. 352 [OII] emitters were identified on the basis of narrow-band excesses and photometric redshifts. We discovered a huge filamentary structure with some clumps traced by [OII] emitters and found that the ClG J0218.3-0510 cluster is embedded in an even larger super-structure than the one reported previously. 31 [OII] emitters were spectroscopically confirmed with the detection of H-alpha and/or [OIII] emission lines by FMOS observations. In the high density regions such as cluster core and clumps, star-forming [OII] emitters show a high overdensity by a factor of more than 10 compared to the field region. Although the star formation activity is very high even in the cluster core, some massive quiescent galaxies also exits at the same time. Furthermore, the properties of the individual [OII] emitters, such as star formation rates, stellar masses and specific star formation rates, do not show a significant dependence on the local density, either. Such lack of environmental dependence is consistent with our earlier result by Hayashi et al. (2011) on a z=1.5 cluster and its surrounding region. The fact that the star-forming activity of galaxies in the cluster core is as high as that in the field at z~1.6 may suggest that the star-forming galaxies are probably just in a transition phase from a starburst mode to a quiescent mode, and are thus showing comparable level of star formation rates to those in lower density environments. We may be witnessing the start of the reversal of the local SFR--density relation due to the "biased" galaxy formation and evolution in high density regions at high this redshift, beyond which massive galaxies would be forming vigorously in a more biased way in proto-cluster cores.

We report the discovery of a protocluster at z~6 containing at least eight cluster member galaxies with spectroscopic confirmations in the wide-field image of the Subaru Deep Field (SDF). The overdensity of the protocluster is significant at the 6 sigma level, based on the surface number density of i'-dropout galaxies. The overdense region covers ~36 sq. arcmin, and includes 30 i'-dropout galaxies. Follow-up spectroscopy revealed that 15 of these are real z~6 galaxies (5.7 < z < 6.3). Eight of the 15 are clustering in a narrow redshift range centered at z=6.01, corresponding to a seven-fold increase in number density over the average in redshift space. We found no significant difference in the observed properties, such as Ly-alpha luminosities and UV continuum magnitudes, between the eight protocluster members and the seven non-members. The velocity dispersion of the eight protocluster members is 647 km/s, which is about three times higher than that predicted by the standard cold dark matter model. This discrepancy could be attributed to the distinguishing three-dimensional distribution of the eight protocluster members. We discuss two possible explanations for this discrepancy: either the protocluster is already mature, with old galaxies at the center, or it is still immature and composed of three subgroups merging to become a larger cluster. In either case, this concentration of z=6.01 galaxies in the SDF may be one of the first sites of formation of a galaxy cluster in the universe.

Using deep narrow-band and broad-band imaging, we identify 401 z~0.40 and 249 z~0.49 H-alpha line-emitting galaxies in the Subaru Deep Field. Compared to other H-alpha surveys at similar redshifts, our samples are unique since they probe lower H-alpha luminosities, are augmented with multi-wavelength (rest-frame 1000AA--1.5 microns) coverage, and a large fraction (20%) of our samples has already been spectroscopically confirmed. Our spectra allow us to measure the Balmer decrement for nearly 60 galaxies with H-beta detected above 5-sigma. The Balmer decrements indicate an average extinction of A(H-alpha)=0.7^{+1.4}_{-0.7} mag. We find that the Balmer decrement systematically increases with higher H-alpha luminosities and with larger stellar masses, in agreement with previous studies with sparser samples. We find that the SFRs estimated from modeling the spectral energy distribution (SED) is reliable---we derived an "intrinsic" H-alpha luminosity which is then reddened assuming the color excess from SED modeling. The SED-predicted H-alpha luminosity agrees with H-alpha narrow-band measurements over 3 dex (rms of 0.25 dex). We then use the SED SFRs to test different statistically-based dust corrections for H-alpha and find that adopting one magnitude of extinction is inappropriate: galaxies with lower luminosities are less reddened. We find that the luminosity-dependent dust correction of Hopkins et al. yields consistent results over 3 dex (rms of 0.3 dex). Our comparisons are only possible by assuming that stellar reddening is roughly half of nebular reddening. The strong correspondence argue that with SED modeling, we can derive reliable intrinsic SFRs even in the absence of H-alpha measurements at z~0.5.

We present Keck spectroscopic observations of z>6 Lyman-break galaxy (LBG) candidates in the Subaru Deep Field (SDF). The candidates were selected as i'-dropout objects down to z'=27 AB magnitudes from an ultra-deep SDF z'-band image. With the Keck spectroscopy we identified 19 LBGs with prominent Ly_alpha emission lines at 6< z < 6.4. The median value of the Ly_alpha rest-frame equivalent widths (EWs) is ~50 A, with four EWs >100 A. This well-defined spectroscopic sample spans a UV-continuum luminosity range of -21.8< M_{UV}<-19.5 (0.6~5 L*_{UV}) and a Ly_alpha luminosity range of (0.3~3) x 10^{43} erg s^{-1} (0.3~3 L*_ {Ly_alpha}). We derive the UV and Ly_alpha luminosity functions (LFs) from our sample at <z>~6.2 after we correct for sample incompleteness. We find that our measurement of the UV LF is consistent with the results of previous studies based on photometric LBG samples at 5<z<7. Our Ly_alpha LF is also generally in agreement with the results of Ly_alpha-emitter surveys at z~5.7 and 6.6. This study shows that deep spectroscopic observations of LBGs can provide unique constraints on both the UV and Ly_alpha LFs at z>6.

We carried out extended spectroscopic confirmations of Ly-alpha emitters (LAEs) at z=6.5 and 5.7 in the Subaru Deep Field. Now, the total number of spectroscopically confirmed LAEs is 45 and 54 at z=6.5 and 5.7, respectively, and at least 81% (70%) of our photometric candidates at z=6.5 (5.7) have been spectroscopically identified as real LAEs. We made careful measurements of the Ly-alpha luminosity, both photometrically and spectroscopically, to accurately determine the Ly-alpha and rest-UV luminosity functions (LFs). The substantially improved evaluation of the Ly-alpha LF at z=6.5 shows an apparent deficit from z=5.7 at least at the bright end, and a possible decline even at the faint end, though small uncertainties remain. The rest-UV LFs at z=6.5 and 5.7 are in good agreement, at least at the bright end, in clear contrast to the differences seen in the Ly-alpha LF. These results imply an increase in the neutral fraction of the intergalactic medium from z=5.7 to 6.5. The rest-frame equivalent width (EW_0) distribution at z=6.5 seems to be systematically smaller than z=5.7, and it shows an extended tail toward larger EW_0. The bright end of the rest-UV LF can be reproduced from the observed Ly-alpha LF and a reasonable EW_0-UV luminosity relation. Integrating this rest-UV LF provides the first measurement of the contribution of LAEs to the photon budget required for reionization. The derived UV LF suggests that the fractional contribution of LAEs to the photon budget among Lyman break galaxies significantly increases towards faint magnitudes. Low-luminosity LAEs could dominate the ionizing photon budget, though this inference depends strongly on the uncertain faint-end slope of the Ly-alpha LF.

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

We conducted a deep narrowband NB973 (FWHM = 200 A centered at 9755 A) survey of z=7 Lyman alpha emitters (LAEs) in the Subaru/XMM-Newton Deep Survey Field, using the fully depleted CCDs newly installed on the Subaru Telescope Suprime-Cam, which is twice more sensitive to z=7 Lyman alpha at ~ 1 micron than the previous CCDs. Reaching the depth 0.5 magnitude deeper than our previous survey in the Subaru Deep Field that led to the discovery of a z=6.96 LAE, we detected three probable z=7 LAE candidates. Even if all the candidates are real, the Lyman alpha luminosity function (LF) at z=7 shows a significant deficit from the LF at z=5.7 determined by previous surveys. The LAE number and Lyman alpha luminosity densities at z=7 is ~ 7.7-54% and ~5.5-39% of those at z=5.7 to the Lyman alpha line luminosity limit of L(Ly-alpha) >~ 9.2 x 10^{42} erg s^{-1}. This could be due to evolution of the LAE population at these epochs as a recent galaxy evolution model predicts that the LAE modestly evolves from z=5.7 to 7. However, even after correcting for this effect of galaxy evolution on the decrease in LAE number density, the z=7 Lyman alpha LF still shows a deficit from z=5.7 LF. This might reflect the attenuation of Lyman alpha emission by neutral hydrogen remaining at the epoch of reionization and suggests that reionization of the universe might not be complete yet at z=7. If we attribute the density deficit to reionization, the intergalactic medium (IGM) transmission for Lyman alpha photons at z=7 would be 0.4 <= T_{Ly-alpha}^{IGM} <= 1, supporting the possible higher neutral fraction at the earlier epochs at z > 6 suggested by the previous surveys of z=5.7-7 LAEs, z ~ 6 quasars and z > 6 gamma-ray bursts.

We obtained mid-infrared 3.6 and 4.5 micron imaging of a z=6.96 Lyman alpha emitter (LAE) IOK-1 discovered in the Subaru Deep Field, using Spitzer Space Telescope Infrared Array Camera observations. After removal of a nearby bright source, we find that IOK-1 is not significantly detected in any of these infrared bands to m_3.6 ~ 24.00 and m_4.5 ~ 23.54 at 3 sigma. Fitting population synthesis models to the spectral energy distribution consisting of the upper limit fluxes of the optical to infrared non-detection images and fluxes in detection images, we constrain the stellar mass M* of IOK-1. This LAE could have either a mass as low as M* <~ 2-9 x 10^8 Msun for the young age (<~ 10 Myr) and the low dust reddening (A_V ~ 0) or a mass as large as M* <~ 1-4 x 10^{10} Msun for either the old age (> 100 Myr) or the high dust reddening (A_V ~ 1.5). This would be within the range of masses of z ~ 3-6.6 LAEs studied to date, ~ 10^6-10^{10} Msun. Hence, IOK-1 is not a particularly unique galaxy with extremely high mass or low mass but is similar to one of the LAEs seen at the later epochs.

We investigate the sky distribution of z ~ 6 Lyman break galaxies selected as i'-dropouts having i' - z' > 1.45 down to z' < 26.5 in the Subaru Deep Field (SDF). We discover 37 i'-dropouts clustered in a projected comoving 21.6 x 21.6 Mpc^2 region at z = 6, showing a local density excess. Carrying out follow-up spectroscopy, we identify four of them as Lyman-alpha emitters at z = 5.92, 6.01, 6.03 and 6.03 (spread over a distance of 46.6 Mpc). The number density of the cluster itself in SDF is ~ 2.2 x 10^{-7} Mpc^{-3}, smaller than those of protoclusters (i.e., forming galaxy clusters) at z ~ 2-5.7. Also, the structure shows ~4-21 times larger galaxy number density than those of z ~ 6 galaxies in a general field. It has a mass of M ~ 1.5^{+1.8}_{-0.5} x 10^{15}M_sun, comparable to those of z ~ 0-5 protoclusters. Since the contamination of our sample by interlopers is estimated to be quite low, 5.9%, most of the other unconfirmed i'-dropouts in the overdense region can be also z ~ 6 galaxies. Hence, it could be a candidate forming cluster at z ~ 6, representing a progenitor of galaxy clusters seen in the recent-day Universe.

We made a narrowband NB973 (bandwidth of 200A at 9755A) imaging of the Subaru Deep Field (SDF) and found two z=7 Lyman alpha emitter (LAE) candidates down to NB973=24.9. Carrying out deep follow-up spectroscopy, we identified one of them as a real z=6.96 LAE. This has shown that galaxy formation was in progress just 750 Myr after the Big Bang. Meanwhile, the Lyman alpha line luminosity function of LAE is known to decline from z=5.7 to 6.6 in the SDF. L* at z=6.6 is 40-60% of that at z=5.7. We also confirm that the number density of z=7 LAE is only 17% of the density at z=6.6 comparing the latest SDF LAE samples. This series of significant decreases in LAE density with increasing redshift can be the result of galaxy evolution during these epochs. However, using the UV continuum luminosity functions of LAEs, those of Lyman break galaxies, and a LAE evolution model based on the hierarchical clustering, we find that galaxy evolution alone cannot explain all the decrease in density. This extra density deficit can be interpreted as the attenuation of the Lyman alpha photons from LAEs due to a rapid evolution of neutral hydrogen fraction during the ongoing cosmic reionization at z~6.6-7.

We carried out a target survey for Lyman break galaxies (LBGs) and Lyman alpha emitters (LAEs) around QSO SDSS J0211-0009 at z=4.87. The deep and wide broadband and narrowband imaging simultaneously revealed the perspective structure of these two high-z populations. The LBGs without Ly-alpha emission form a filamentary structure including the QSO, while the LAEs are distributed around the QSO but avoid it within a distance of ~4.5Mpc. On the other hand, we serendipitously discovered a protocluster with a significant concentration of LBGs and LAEs where no strongly UV ionizing source such as a QSO or radio galaxy is known to exist. In this cluster field, two populations are spatially cross-correlated with each other. The relative spatial distribution of LAEs to LBGs is in stark contrast between the QSO and the cluster fields. We also found a weak trend showing that the number counts based on Ly-alpha and UV continuum fluxes of LAEs in the QSO field are slightly lower than in the cluster field, whereas the number counts of LBGs are almost consistent with each other. The LAEs avoid the nearby region around the QSO where the local UV background radiation could be ~100 times stronger than the average for the epoch. The clustering segregation between LBGs and LAEs seen in the QSO field could be due to either enhanced early galaxy formation in an overdense environment having caused all the LAEs to evolve into LBGs, or local photoionization due to the strong UV radiation from the QSO effectively causing a deficit in low-mass galaxies like LAEs.

When galaxy formation started in the history of the Universe remains unclear. Studies of the cosmic microwave background indicate that the Universe, after initial cooling (following the Big Bang), was reheated and reionized by hot stars in newborn galaxies at a redshift in the range 6 < z < 14 (ref. 1). Though several candidate galaxies at redshift z > 7 have been identified photometrically (refs 2,3), galaxies with spectroscopically confirmed redshifts have been confined to z < 6.6 (refs. 4-8). Here we report a spectroscopic redshift of z = 6.96 (corresponding to just 750 Myr after the Big Bang) for a galaxy whose spectrum clearly shows Lyman-alpha emission at 9,682 A, indicating active star formation at a rate of about 10 M_sun/yr, where M_sun us the mass of the Sun. This demonstrates that galaxy formation was under way when the Universe was only about 6 per cent of its present age. The number density of galaxies at z = 7 seems to be only 18-36 per cent of the density at z = 6.6.

We investigate the clustering properties of 2,600 Lyman Break Galaxies (LBGs) at z=3.5-5.2 in two large blank fields, the Subaru Deep Field and the Subaru/XMM Deep Field (600arcmin^2 each). The angular correlation functions of these LBGs show a clear clustering at both z~4 and 5. The correlation lengths are r_0= 4.1^{+0.2}_{-0.2} and 5.9^{+1.3}_{-1.7} h_{100}^{-1} Mpc (r_0= 5.1^{+1.0}_{-1.1} and 5.9^{+1.3}_{-1.7} h_{100}^{-1} Mpc) for all the detected LBGs (for L>L* LBGs) at z~4 and 5, respectively. These correlation lengths correspond to galaxy-dark matter biases of b_g= 2.9^{+0.1}_{-0.1} and 4.6^{+0.9}_{-1.2} (b_g=3.5^{+0.6}_{-0.7} and 4.6^{+0.9}_{-1.2}), for all the detected LBGs (for L>L^* LBGs) at z~4 and 5, respectively. These results, combined with estimates for z~3 LBGs in the literature, show that the correlation length of L>L^* LBGs is almost constant, ~5 h_{100}^{-1} Mpc, over z~3-5, while the bias monotonically increases with redshift at z>3. We also find that for LBGs at z~4 the clustering amplitude increases with UV-continuum luminosity and with the amount of dust extinction. We estimate the mass of dark halos hosting various kinds of high-z galaxies including LBGs with the analytic model given by Sheth & Tormen (1999). We find that the typical mass of dark halos hosting L>L^* LBGs is about 1x10^{12} h_{70}^{-1}Msol over z~3-5, which is comparable to that of the Milky Way Galaxy. A single dark halo with ~10^{12} h_{70}^{-1} Msol is found to host 0.1-0.3 LBG on average but host about four K-band selected galaxies.

(abridged) We investigate photometric properties of Lyman Break Galaxies (LBGs) at z=3.5-5.2 based on large samples of 2,600 LBGs detected in deep (i'~27) and wide-field (1,200 arcmin^2) images taken in the Subaru Deep Field (SDF) and the Subaru/XMM Deep Field (SXDF). The selection criteria for the LBG samples are examined with 85 spectroscopically identified objects and by Monte Carlo simulations. We find in the luminosity functions of LBGs (i) that the number density of bright galaxies (M_{1700}<-22; corresponding to SFR_{corr}>100 Msolar yr^{-1}) decreases significantly from z=4 to 5 and (ii) that the faint-end slope of the luminosity function may become steeper towards higher redshifts. We estimate dust extinction of z=4 LBGs with M<M^* from UV slopes, and obtain E(B-V)=0.15+/-0.03 as the mean value. The dust extinction remains constant with apparent luminosity, but increases with intrinsic luminosity. We find no evolution in dust extinction between LBGs at z=3 and 4. We investigate the evolution of UV-luminosity density at 1700A, rho, and find that rho does not significantly change from z=3 to z=5, i.e., rho(z=4)/rho(z=3)=1.0+/-0.2 and rho(z=5)/rho(z=3)=0.8+/-0.4, thus the cosmic star-formation rate (SFR) density remains constant. We find that the stellar mass density estimated from the cosmic SFR is consistent with those derived directly from the stellar mass function at z=0-1, but exceeds those at z~3 by a factor of 3. We find that the ratio of the UV-luminosity density of Ly-a emitters (LAEs) to that of LBGs is ~60% at z=5, and thus about a half of the star formation at z=5 probably occurs in LAEs. We obtain a constraint on the escape fraction of UV-ionizing photons produced by LBGs, f_{esc}>~0.13.

We present spectropolarimetry of the Type Ic supernova SN 2002ap and give a preliminary analysis: the data were taken at two epochs, close to and one month later than the visual maximum (2002 February 8). In addition we present June 9 spectropolarimetry without analysis. The data show the development of linear polarization. Distinct polarization profiles were seen only in the O I \lambda 7773 multiplet/Ca II IR triplet absorption trough at maximum light and in the Ca II IR triplet absorption trough a month later, with the latter showing a peak polarization as high as ~2 %. The intrinsic polarization shows three clear position angles: 80 degs for the February continuum, 120 degs for the February line feature, and 150 degs for the March data. We conclude that there are multiple asymmetric components in the ejecta. We suggest that the supernova has a bulk asymmetry with an axial ratio projected on the sky that is different from 1 by of order 10 %. Furthermore, we suggest very speculatively that a high velocity ejecta component moving faster than ~0.115c (e.g., a jet) contributes to polarization in the February epoch.