We report the rest-frame ultraviolet luminosity function of $g$-dropout galaxies in 177 protocluster candidates (PC UVLF) at $z\sim4$ selected in the Hyper Suprime-Cam Subaru Strategic Program. Comparing with the UVLF of field galaxies at the same redshift, we find that the PC UVLF shows a significant excess towards the bright-end. This excess can not be explained by the contribution of only active galactic nuclei, and we also find that this is more significant in higher dense regions. Assuming that all protocluster members are located on the star formation main sequence, the PC UVLF can be converted into a stellar mass function. Consequently, our protocluster members are inferred to have a 2.8 times more massive characteristic stellar mass than that of the field Lyman break galaxies at the same redshift. This study, for the first time, clearly shows that the enhancement in star formation or stellar mass in overdense regions can generally be seen as early as at $z\sim4$. We also estimate the star formation rate density (SFRD) in protocluster regions as $\simeq 6-20\%$ of the cosmic SFRD, based on the measured PC UVLF after correcting for the selection incompleteness in our protocluster sample. This high value suggests that protoclusters make a non-negligible contribution to the cosmic SFRD at $z\sim4$, as previously suggested by simulations. Our results suggest that protoclusters are essential components for the galaxy evolution at $z\sim4$.

We present the galaxy luminosity functions (LFs) of four Hickson Compact Groups using image data from the Subaru Hyper Suprime-Cam. A distinct dip appeared in the faint-ends of all the LFs at $M_g\sim-12$. A similar dip was observed in the LFs of the galaxy clusters Coma and Centaurus. However, LFs in the Virgo, Hydra, and the field had flatter slopes and no dips. As the relative velocities among galaxies are lower in compact groups than in clusters, the effect of galaxy-galaxy interactions would be more significant in compact groups. The $M_g\sim-12$ dip of compact groups may imply that frequent galaxy-galaxy interactions would affect the evolution of galaxies, and the dip in LF could become a boundary between different galaxy populations.

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

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

We use MMT spectroscopy and deep Subaru Hyper Suprime-Cam (HSC) imaging to compare the spectroscopic central stellar velocity dispersion of quiescent galaxies with the effective dispersion of the dark matter halo derived from the stacked lensing signal. The spectroscopic survey (the Smithsonian Hectospec Lensing Survey) provides a sample of 4585 quiescent galaxy lenses with measured line-of-sight central stellar velocity dispersion ($\sigma_{\rm SHELS}$) that is more than 85% complete for $R < 20.6$, $D_{n}4000> 1.5$ and $M_* > 10^{9.5}{\rm M}_{\odot}$. The median redshift of the sample of lenses is 0.32. We measure the stacked lensing signal from the HSC deep imaging. The central stellar velocity dispersion is directly proportional to the velocity dispersion derived from the lensing $\sigma_{\rm Lens}$, $\sigma_{\rm Lens} = (0.98\pm 0.14)\sigma_{\rm SHELS}+(19.89\pm31.89)$. The independent spectroscopic and weak lensing velocity dispersions probe different scales, $\sim3$kpc and $\gtrsim$ 100 kpc, respectively, and strongly indicate that the observable central stellar velocity dispersion for quiescent galaxies is a good proxy for the velocity dispersion of the dark matter halo. We thus demonstrate the power of combining high-quality imaging and spectroscopy to shed light on the connection between galaxies and their dark matter halos.

Green valley galaxies represent the population that is likely to transition from the star-forming to the quiescent phases. To investigate the role of the environment in quenching star formation, we use the wide-field data from the Hyper Suprime-Cam Strategic Subaru Proposal survey to quantify the frequency of green valley galaxies in different environments and their redshift evolution. We find that the green valley fraction, in general, is less than 20% in any redshift and environment. The green valley fraction, when normalized to the total population, is higher in the field than that in groups or clusters and decreases with a decreasing redshift and increasing mass. The lower fraction of transitional galaxies in denser environments could be a consequence of the lack of star-forming galaxies, which could be the progenitors of green valley galaxies. To assess the effect of the environment on star formation quenching, we define the effective green valley fraction as the ratio of the number of green valley galaxies to that of nonquiescent galaxies only. The effective green valley fraction for field galaxies is lower than that for group or cluster galaxies, which reveals a strong positive mass dependence and mild redshift evolution. Moreover, the specific star formation rate (sSFR) is reduced by 0.1-0.3 dex in groups or clusters. Our results thus imply that an ongoing slow quenching process has been acting in the dense environment since z~1.

We present proper motion measurements for more than $0.55$ million main-sequence stars, by comparing astrometric positions of matched stars between the multi-band imaging datasets from the Hyper Suprime-Cam (HSC) Survey and the SDSS Stripe 82. In doing this we use $3$ million galaxies to recalibrate the astrometry and set up a common reference frame between the two catalogues. The exquisite depth and the nearly $12$ years of time baseline between HSC and SDSS enable high-precision measurements of statistical proper motions for stars down to $i\simeq 24$. A validation of our method is demonstrated by the agreement with the $Gaia$ proper motions, to the precision better than $0.1$ mas yr$^{-1}$. To retain the precision, we make a correction of the subtle effects due to the differential chromatic refraction in the SDSS images based on the comparison with the $Gaia$ proper motions against colour of stars, which is validated using the SDSS spectroscopic quasars. Combining with the photometric distance estimates for individual stars based on the precise HSC photometry, we show a significant detection of the net proper motions for stars in each bin of distance out to $100$ kpc. The two-component tangential velocities after subtracting the apparent motions due to our own motion display rich phase-space structures including a clear signature of the Sagittarius stream in the halo region of distance range $[10,\ 35]$ kpc. We also measure the tangential velocity dispersion in the distance range $5-20$ kpc and find that the data are consistent with a constant isotropic dispersion of $80\pm 10\ {\rm km/s}$. More distant stars appear to have random motions with respect to the Galactic centre on average.

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

We present measurements of cosmic shear two-point correlation functions (TPCFs) from Hyper Suprime-Cam Subaru Strategic Program (HSC SSP) first-year data, and derived cosmological constraints based on a blind analysis. The HSC first-year shape catalog is divided into four tomographic redshift bins ranging from $z=0.3$ to 1.5 with equal widths of $\Delta z =0.3$. The unweighted galaxy number densities in each tomographic bin are 5.9, 5.9, 4.3, and 2.4 arcmin$^{-2}$ from lower to higher redshifts, respectively. We adopt the standard TPCF estimators, $\xi_\pm$, for our cosmological analysis, given that we find no evidence of the significant B-mode shear. The TPCFs are detected at high significance for all ten combinations of auto- and cross-tomographic bins over a wide angular range, yielding a total signal-to-noise ratio of 19 in the angular ranges adopted in the cosmological analysis, $7'<\theta<56'$ for $\xi_+$ and $28'<\theta<178'$ for $\xi_-$. We perform the standard Bayesian likelihood analysis for cosmological inference from the measured cosmic shear TPCFs, including contributions from intrinsic alignment of galaxies as well as systematic effects from PSF model errors, shear calibration uncertainty, and source redshift distribution errors. We adopt a covariance matrix derived from realistic mock catalogs constructed from full-sky gravitational lensing simulations that fully account for survey geometry and measurement noise. For a flat $\Lambda$ cold dark matter model, we find $S_8 \equiv \sigma_8\sqrt{\Omega_m/0.3}=0.804_{-0.029}^{+0.032}$, and $\Omega_m=0.346_{-0.100}^{+0.052}$. We carefully check the robustness of the cosmological results against astrophysical modeling uncertainties and systematic uncertainties in measurements, and find that none of them has a significant impact on the cosmological constraints.

Like massive galaxies, dwarf galaxies are expected to undergo major mergers with other dwarfs. However, the end state of these mergers and the role that merging plays in regulating dwarf star formation is uncertain. Using imaging from the Hyper Suprime-Cam Subaru Strategic program, we construct a sample of dwarf-dwarf mergers and examine the star formation and host properties of the merging systems. These galaxies are selected via an automated detection algorithm from a sample of 6875 spectroscopically selected isolated dwarf galaxies at $z<0.12$ and $\log(M_\star/M_\odot)<9.6$ from the Galaxy and Mass Assembly (GAMA) and Sloan Digital Sky Survey (SDSS) spectroscopic campaigns. We find a total tidal feature detection fraction of 3.29% (6.1% when considering only galaxies at $z<0.05$). The tidal feature detection fraction rises strongly as a function of star formation activity; 15%-20% of galaxies with extremely high H$\alpha$ equivalent width (H$\alpha$ EW > 250 Angstrom) show signs of tidal debris. Galaxies that host tidal debris are also systematically bluer than the average galaxy at fixed stellar mass. These findings extend the observed dwarf-dwarf merger sequence with a significant sample of dwarf galaxies, indicating that star formation triggered in mergers between dwarf galaxies continues after coalescence.

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

We report the discovery of 28 quasars and 7 luminous galaxies at 5.7 $\le$ z $\le$ 7.0. This is the tenth in a series of papers from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The total number of spectroscopically identified objects in SHELLQs has now grown to 93 high-z quasars, 31 high-z luminous galaxies, 16 [O III] emitters at z ~ 0.8, and 65 Galactic cool dwarfs (low-mass stars and brown dwarfs). These objects were found over 900 deg2, surveyed by HSC between 2014 March and 2018 January. The full quasar sample includes 18 objects with very strong and narrow Ly alpha emission, whose stacked spectrum is clearly different from that of other quasars or galaxies. While the stacked spectrum shows N V 1240 emission and resembles that of lower-z narrow-line quasars, the small Ly alpha width may suggest a significant contribution from the host galaxies. Thus these objects may be composites of quasars and star-forming galaxies.

We present the large-scale structure over more than 50 comoving Mpc scale at z $\sim$ 0.9 where the CL1604 supercluster, which is one of the largest structures ever known at high redshifts, is embedded. The wide-field deep imaging survey by the Subaru Strategic Program with Hyper Suprime-Cam reveals that the already-known CL1604 supercluster is a mere part of larger-scale structure extending to both the north and the south. We confirm that there are galaxy clusters at three slightly different redshifts in the northern and southern sides of the supercluster by determining the redshifts of 55 red-sequence galaxies and 82 star-forming galaxies in total by the follow-up spectroscopy with Subaru/FOCAS and Gemini-N/GMOS. This suggests that the structure ever known as the CL1604 supercluster is the tip of the iceberg. We investigate stellar population of the red-sequence galaxies using 4000 A break and Balmer H$\delta$ absorption line. Almost all of the red-sequence galaxies brighter than 21.5 mag in $z$-band show an old stellar population with $\gtrsim2$ Gyr. The comparison of composite spectra of the red-sequence galaxies in the individual clusters show that the galaxies at a similar redshift have similar stellar population age, even if they are located $\sim$50 Mpc apart from each other. However, there could be a large variation in the star formation history. Therefore, it is likely that galaxies associated with the large-scale structure at 50 Mpc scale formed at almost the same time, have assembled into the denser regions, and then have evolved with different star formation history along the hierarchical growth of the cosmic web.

We report on the discovery of a new Milky Way (MW) satellite in Bo\"otes based on data from the on-going Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP). This satellite, named Bo\"otes IV, is the third ultra-faint dwarf that we have discovered in the HSC-SSP. We have identified a statistically significant (32.3$\sigma$) overdensity of stars having characteristics of a metal-poor, old stellar population. The distance to this stellar system is $D_{\odot}=209^{+20}_{-18}$ kpc with a $V$-band absolute magnitude of $M_V=-4.53^{+0.23}_{-0.21}$ mag. Bo\"otes IV has a half-light radius of $r_h=462^{+98}_{-84}$ pc and an ellipticity of $0.64^{+0.05}_{-0.05}$, which clearly suggests that this is a dwarf satellite galaxy. We also found another overdensity that appears to be a faint globular cluster with $M_V=-0.20^{+0.59}_{-0.83}$ mag and $r_h=5.9^{+1.5}_{-1.3}$ pc located at $D_{\odot}=46^{+4}_{-4}$ kpc. Adopting the recent prediction for the total population of satellites in a MW-sized halo by Newton et al. (2018), which combined the characteristics of the observed satellites by SDSS and DES with the subhalos obtained in $\Lambda$CDM models, we estimate that there should be about two MW satellites at $M_V\le0$ in the $\sim676$ deg$^2$ covered by HSC-SSP, whereas that area includes six satellites. Thus, the observed number of satellites is larger than the theoretical prediction. On the face of it, we have a problem of too many satellites, instead of the well-known missing satellites problem whereby the $\Lambda$CDM theory overpredicts the number of satellites in a MW-sized halo. This may imply that the models need more refinements for the assignment of subhalos to satellites such as considering those found by the current deeper survey. [abridged]

We report on the global structure of the Milky Way (MW) stellar halo up to its outer boundary based on the analysis of blue-horizontal branch stars (BHBs). These halo tracers are extracted from the $(g,r,i,z)$ band multi-photometry in the internal data release of the on-going Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) surveyed over $\sim550$~deg$^2$ area. In order to select most likely BHBs by removing blue straggler stars (BSs) and other contamination in a statistically significant manner, we have developed and applied an extensive Bayesian method, instead of the simple color cuts adopted in our previous work, where each of the template BHBs and non-BHBs obtained from the available catalogs is represented as a mixture of multiple Gaussian distributions in the color-color diagrams. We found from the candidate BHBs in the range of 18.5<g<23.5 mag that the radial density distribution over a Galactocentric radius of r=36-360 kpc can be approximated as a single power-law profile with an index of $\alpha=3.74^{+0.21}_{-0.22}$ or a broken power-law profile with an index of $\alpha_{\rm in}=2.92^{+0.33}_{-0.33}$ at $r$ below a broken radius of $r_{\rm b}=160^{+18}_{-19}$~kpc and a very steep slope of $\alpha_{\rm out}=15.0^{+3.7}_{-4.5}$ at $r>r_{\rm b}$. The latter profile with a prolate shape having an axial ratio of $q=1.72^{+0.44}_{-0.28}$ is most likely and this halo may hold a rather sharp boundary at r=160kpc. The slopes of the halo density profiles are compared with those from the suite of hydrodynamical simulations for the formation of stellar halos. This comparison suggests that the MW stellar halo may consist of the two overlapping components: the in situ. inner halo as probed by RR Lyrae stars showing a relatively steep radial density profile and the ex situ. outer halo with a shallow profile probed by BHBs here, which is made by accretion of small stellar systems.

We present results of the cross-correlation Lya intensity mapping with Subaru/Hyper Suprime-Cam (HSC) ultra-deep narrowband images and Lya emitters (LAEs) at $z=5.7$ and $6.6$ in a total area of $4$ deg$^2$. Although overwhelming amount of data quality controls have been performed for the narrowband images and the LAE samples, we further conduct extensive analysis evaluating systematics of large-scale point spread function wings, sky subtractions, and unknown errors on the basis of physically uncorrelated signals and sources found in real HSC images and object catalogs, respectively. Removing the systematics, we carefully calculate cross-correlations between Lya intensity of the narrowband images and the LAEs. We identify very diffuse Lya emission with the $3\sigma$ ($2\sigma$) significance at $>150$ comoving kpc (ckpc) far from the LAEs at $z = 5.7$ ($6.6$), beyond a virial radius of star-forming galaxies with $M_\mathrm{h} \sim 10^{11} M_\odot$. The diffuse Lya emission possibly extends up to $1,000$ ckpc with the surface brightness of $10^{-20}-10^{-19}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$. We confirm that the small-scale ($<150$ ckpc) Lya radial profiles of LAEs in our Lya intensity maps are consistent with those obtained by recent MUSE observations (Leclercq et al. 2017). Comparisons with numerical simulations suggest that the large-scale ($\sim 150-1,000$ ckpc) Lya emission are not explained by unresolved faint sources of neighboring galaxies including satellites, but by a combination of Lya photons emitted from the central LAE and other unknown sources, such as a cold-gas stream and galactic outflow. We find no evolution in the Lya radial profiles of our LAEs from $z=5.7$ to $6.6$, where theoretical models predict a flattening of the profile slope made by cosmic reionization, albeit with our moderately large observational errors.

We study the faint stellar halo of isolated central galaxies, by stacking galaxy images in the HSC survey and accounting for the residual sky background sampled with random points. The surface brightness profiles in HSC $r$-band are measured for a wide range of galaxy stellar masses ($9.2<\log_{10}M_\ast/M_\odot<11.4$) and out to 120 kpc. Failing to account for the stellar halo below the noise level of individual images will lead to underestimates of the total luminosity by $\leq 15\%$. Splitting galaxies according to the concentration parameter of their light distributions, we find that the surface brightness profiles of low concentration galaxies drop faster between 20 and 100 kpc than those of high concentration galaxies. Albeit the large galaxy-to-galaxy scatter, we find a strong self-similarity of the stellar halo profiles. They show unified forms once the projected distance is scaled by the halo virial radius. The colour of galaxies is redder in the centre and bluer outside, with high concentration galaxies having redder and more flattened colour profiles. There are indications of a colour minimum, beyond which the colour of the outer stellar halo turns red again. This colour minimum, however, is very sensitive to the completeness in masking satellite galaxies. We also examine the effect of the extended PSF in the measurement of the stellar halo, which is particularly important for low mass or low concentration galaxies. The PSF-corrected surface brightness profile can be measured down to $\sim$31 $\mathrm{mag}/\mathrm{arcsec}^2$ at 3-$\sigma$ significance. PSF also slightly flattens the measured colour profiles.

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

We present an overview of a deep transient survey of the COSMOS field with the Subaru Hyper Suprime-Cam (HSC). The survey was performed for the 1.77 deg$^2$ ultra-deep layer and 5.78 deg$^2$ deep layer in the Subaru Strategic Program over 6- and 4-month periods from 2016 to 2017, respectively. The ultra-deep layer shows a median depth per epoch of 26.4, 26.3, 26.0, 25.6, and 24.6 mag in $g$, $r$, $i$, $z$, and $y$ bands, respectively; the deep layer is $\sim0.6$ mag shallower. In total, 1,824 supernova candidates were identified. Based on light curve fitting and derived light curve shape parameter, we classified 433 objects as Type Ia supernovae (SNe); among these candidates, 129 objects have spectroscopic or COSMOS2015 photometric redshifts and 58 objects are located at $z > 1$. Our unique dataset doubles the number of Type Ia SNe at $z > 1$ and enables various time-domain analyses of Type II SNe, high redshift superluminous SNe, variable stars, and active galactic nuclei.

We report on the results of deep and wide-field (1.1 deg$^2$) narrow-band observations with Subaru/Hyper Suprime-Cam (HSC) of a field around a hyperluminous QSO (HLQSO), HS1549+1919, residing in a protocluster at $z=2.84$, to map the large-scale structure of Ly$\alpha$ emitters (LAEs). One HSC pointing enables us to detect 3490 LAEs and 76 extended Ly$\alpha$ blobs (LABs), probing diverse environments from voids to protoclusters. The HLQSO is found to be near the center of the protocluster, which corresponds to the intersection of $\sim$100 cMpc-scale structures of LAEs. LABs are basically distributed along the large-scale structure, with larger ones particularly clustered around the HLQSO, confirming a previously noted tendency of LABs to prefer denser environments. Moreover, the shapes of LABs near the HLQSO appear to be aligned with the large-scale structure. Finally, a deep Ly$\alpha$ image reveals a diffuse Ly$\alpha$ nebula along a filamentary structure with no luminous UV/sub-mm counterpart. We suggest that the diffuse nebula is due to a cold filament with high clumping factor illuminated by the QSO, with a required high clumpiness provided by unresolved residing halos of mass $\leq 10^{9-10}M_\odot$.

We report on the optical properties of infrared (IR)-bright dust-obscured galaxies (DOGs) that are defined as (i-[22])_AB > 7.0. Because supermassive black holes (SMBHs) in IR-bright DOGs are expected to be rapidly growing in the major merger scenario, they provide useful clues for understanding the co-evolution of SMBHs and their host galaxies. However, the optical properties of IR-bright DOGs remain unclear because the optical emission of a DOG is very faint. By combining ~ 105 deg^2 images of the optical, near-IR, and mid-IR data obtained from the Subaru Hyper Suprime-Cam (HSC) survey, the VISTA VIKING survey, and the WISE all-sky survey, respectively, 571 IR-bright DOGs were selected. We found that IR-bright DOGs show a redder (g-z)_AB color than other populations of dusty galaxies, such as ultra-luminous IR galaxies (ULIRGs) at a similar redshift, with a significantly large dispersion. Among the selected DOGs, star-formation (SF) dominated DOGs show a relatively red color, while active galactic nucleus (AGN) dominated DOGs show a rather blue color in optical. This result is consistent with the idea that the relative AGN contribution in the optical emission becomes more significant at a later stage in the major merger scenario. We discovered eight IR-bright DOGs showing a significant blue excess in blue HSC bands (BluDOGs). This blue excess can be interpreted as a leaked AGN emission that is either a directly leaking or a scattered AGN emission, as proposed for some blue-excess Hot DOGs in earlier studies.

We present the results of a survey of the brightest UV-selected galaxies in protoclusters. These proto-brightest cluster galaxy (proto-BCG) candidates are drawn from 179 overdense regions of $g$-dropout galaxies at $z\sim4$ from the Hyper Suprime-Cam Subaru Strategic Program identified previously as good protocluster candidates. This study is the first to extend the systematic study of the progenitors of BCGs from $z\sim2$ to $z\sim4$. We carefully remove possible contaminants from foreground galaxies and, for each structure, we select the brightest galaxy that is at least 1 mag brighter than the fifth brightest galaxy. We select 63 proto-BCG candidates and compare their properties with those of galaxies in the field and those of other galaxies in overdense structures. The proto-BCG candidates and their surrounding galaxies have different rest-UV color $(i - z)$ distributions to field galaxies and other galaxies in protoclusters that do not host proto-BCGs. In addition, galaxies surrounding proto-BCGs are brighter than those in protoclusters without proto-BCGs. The image stacking analysis reveals that the average effective radius of proto-BCGs is $\sim28\%$ larger than that of field galaxies. The $i-z$ color differences suggest that proto-BCGs and their surrounding galaxies are dustier than other galaxies at $z\sim4$. These results suggest that specific environmental effects or assembly biasses have already emerged in some protoclusters as early as $z \sim 4$, and we suggest that proto-BCGs have different star formation histories than other galaxies in the same epoch.

We report our first discoveries of high-redshift supernovae from the Subaru HIgh-Z sUpernova CAmpaign (SHIZUCA), the transient survey using Subaru/Hyper Suprime-Cam. We report the discovery of three supernovae at the spectroscopically-confirmed redshifts of 2.399 (HSC16adga), 1.965 (HSC17auzg), and 1.851 (HSC17dbpf), and two supernova candidates with the host-galaxy photometric redshifts of 3.2 (HSC16apuo) and 4.2 (HSC17dsid), respectively. In this paper, we present their photometric properties and the spectroscopic properties of the confirmed high-redshift supernovae are presented in the accompanying paper (Curtin et al. 2019). The supernovae with the confirmed redshifts of z ~ 2 have the rest ultraviolet peak magnitudes close to -21 mag and they are likely superluminous supernovae. The discovery of three supernovae at z ~ 2 roughly corresponds to the approximate event rate of ~ 900 +/- 520 Gpc-3 yr-1 with Poisson error, which is consistent with the total superluminous supernova rate estimated by extrapolating the local rate based on the cosmic star-formation history. Adding unconfirmed superluminous supernova candidates would increase the event rate. Our superluminous supernova candidates at the redshifts of around 3 and 4 indicate the approximate superluminous supernova rates of ~ 400 +/- 400 Gpc-3 yr-1 (z ~ 3) and ~ 500 +/- 500 Gpc-3 yr-1 (z ~ 4) with Poisson errors. Our initial results demonstrate the outstanding capability of Hyper Suprime-Cam to discover high-redshift supernovae.

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

We measure cosmic weak lensing shear power spectra with the Subaru Hyper Suprime-Cam (HSC) survey first-year shear catalog covering 137deg$^2$ of the sky. Thanks to the high effective galaxy number density of $\sim$17 arcmin$^{-2}$ even after conservative cuts such as magnitude cut of $i<24.5$ and photometric redshift cut of $0.3\leq z \leq 1.5$, we obtain a high significance measurement of the cosmic shear power spectra in 4 tomographic redshift bins, achieving a total signal-to-noise ratio of 16 in the multipole range $300 \leq \ell \leq 1900$. We carefully account for various uncertainties in our analysis including the intrinsic alignment of galaxies, scatters and biases in photometric redshifts, residual uncertainties in the shear measurement, and modeling of the matter power spectrum. The accuracy of our power spectrum measurement method as well as our analytic model of the covariance matrix are tested against realistic mock shear catalogs. For a flat $\Lambda$ cold dark matter ($\Lambda$CDM) model, we find $S_8\equiv \sigma_8(\Omega_{\rm m}/0.3)^\alpha=0.800^{+0.029}_{-0.028}$ for $\alpha=0.45$ ($S_8=0.780^{+0.030}_{-0.033}$ for $\alpha=0.5$) from our HSC tomographic cosmic shear analysis alone. In comparison with Planck cosmic microwave background constraints, our results prefer slightly lower values of $S_8$, although metrics such as the Bayesian evidence ratio test do not show significant evidence for discordance between these results. We study the effect of possible additional systematic errors that are unaccounted in our fiducial cosmic shear analysis, and find that they can shift the best-fit values of $S_8$ by up to $\sim 0.6\sigma$ in both directions. The full HSC survey data will contain several times more area, and will lead to significantly improved cosmological constraints.

We present new measurements of the quasar luminosity function (LF) at $z \sim 6$, over an unprecedentedly wide range of the rest-frame ultraviolet luminosity $M_{1450}$ from $-30$ to $-22$ mag. This is the fifth in a series of publications from the Subaru High-$z$ Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The LF was calculated with a complete sample of 110 quasars at $5.7 \le z \le 6.5$, which includes 48 SHELLQs quasars discovered over 650 deg$^2$, and 63 brighter quasars discovered by the Sloan Digital Sky Survey and the Canada-France-Hawaii Quasar Survey (including one overlapping object). This is the largest sample of $z \sim 6$ quasars with a well-defined selection function constructed to date, and has allowed us to detect significant flattening of the LF at its faint end. A double power-law function fit to the sample yields a faint-end slope $\alpha = -1.23^{+0.44}_{-0.34}$, a bright-end slope $\beta = -2.73^{+0.23}_{-0.31}$, a break magnitude $M_{1450}^* = -24.90^{+0.75}_{-0.90}$, and a characteristic space density $\Phi^* = 10.9^{+10.0}_{-6.8}$ Gpc$^{-3}$ mag$^{-1}$. Integrating this best-fit model over the range $-18 < M_{1450} < -30$ mag, quasars emit ionizing photons at the rate of $\dot{n}_{\rm ion} = 10^{48.8 \pm 0.1}$ s$^{-1}$ Mpc$^{-3}$ at $z = 6.0$. This is less than 10 % of the critical rate necessary to keep the intergalactic medium ionized, which indicates that quasars are not a major contributor to cosmic reionization.

We present optical long-slit spectroscopy and far-ultraviolet to near-infrared spectral energy distribution fitting of two diffuse dwarf galaxies, LSBG-285 and LSBG-750, which were recently discovered by the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). We measure redshifts using H$\alpha$ line emission, and find that these galaxies are at comoving distances of ${\approx}25$ and ${\approx}41$ Mpc, respectively, after correcting for the local velocity field. They have effective radii of $r_\mathrm{eff}=1.2$ and 1.8 kpc and stellar masses of $M_\star\approx2$-$3\times10^{7}~M_\odot$. There are no massive galaxies ($M_\star>10^{10} M_\odot$) within a comoving separation of at least 1.5 Mpc from LSBG-285 and 2 Mpc from LSBG-750. These sources are similar in size and surface brightness to ultra-diffuse galaxies, except they are isolated, star-forming objects that were optically selected in an environmentally blind survey. Both galaxies likely have low stellar metallicities $[Z_\star/Z_\odot] < -1.0$ and are consistent with the stellar mass-metallicity relation for dwarf galaxies. We set an upper limit on LSBG-750's rotational velocity of ${\sim}50$ km s$^{-1}$, which is comparable to dwarf galaxies of similar stellar mass with estimated halo masses $<10^{11}~M_\odot$. We find tentative evidence that the gas-phase metallicities in both of these diffuse systems are high for their stellar mass, though a statistically complete, optically-selected galaxy sample at very low surface brightness will be necessary to place these results into context with the higher-surface-brightness galaxy population.

We investigate the local and line-of-sight overdensities of strong gravitational lens galaxies using wide-area multiband imaging from the Hyper Suprime-Cam Subaru Strategic Program. We present 41 new definite or probable lens candidates discovered in Data Release 2 of the survey. Using a combined sample of 87 galaxy-scale lenses out to a lens redshift of $z_{\mathrm{L}} \sim 0.8$, we compare galaxy number counts in lines of sight toward known and newly-discovered lenses in the survey to those of a control sample consisting of random lines of sight. We also compare the local overdensity of lens galaxies to a sample of "twin" galaxies that have a similar redshift and velocity dispersion to test whether lenses lie in different environments from similar non-lens galaxies. We find that lens fields contain higher number counts of galaxies compared to the control fields, but this effect arises from the local environment of the lens. Once galaxies in the lens plane are removed, the lens lines of sight are consistent with the control sample. The local environments of the lenses are overdense compared to the control sample, and are slightly overdense compared to those of the twin sample, although the significance is marginal. There is no significant evidence of the evolution of the local overdensity of lens environments with redshift.

We present the Ly$\alpha$ luminosity function (LF) derived from 34 Ly$\alpha$ emitters (LAEs) at $z=7.0$ on the sky of $3.1$ deg$^2$, the largest sample compared to those in the literature obtained at a redshift $z\gtrsim7$. The LAE sample is made by deep large-area Subaru narrowband observations conducted by the Cosmic HydrOgen Reionization Unveiled with Subaru (CHORUS) project. The $z=7.0$ Ly$\alpha$ LF of our project is consistent with those of the previous DECam and Subaru studies at the bright and faint ends, respectively, while our $z=7.0$ Ly$\alpha$ LF has uncertainties significantly smaller than those of the previous study results. Exploiting the small errors of our measurements, we investigate the shape of the faint to bright-end Ly$\alpha$ LF. We find that the $z=7.0$ Ly$\alpha$ LF shape can be explained by the steep slope of $\alpha \simeq -2.5$ suggested at $z=6.6$, and that there is no clear signature of a bright-end excess at $z\simeq 7$ claimed by the previous work, which was thought to be made by the ionized bubbles around bright LAEs whose Ly$\alpha$ photons could easily escape from the partly neutral IGM at $z \simeq 7$. We estimate the Ly$\alpha$ luminosity densities (LDs) with Ly$\alpha$ LFs at $z\simeq 6-8$ given by our and the previous studies, and compare the evolution of the UV-continuum LD estimated with dropouts. The Ly$\alpha$ LD monotonically decreases from $z\sim 6$ to $8$, and evolves stronger than the UV-continuum LD, indicative of the Ly$\alpha$ damping wing absorption of the IGM towards the heart of the reionization epoch.

The Initial Mass Function (IMF) for massive galaxies can be constrained by combining stellar dynamics with strong gravitational lensing. However, this method is limited by degeneracies between the density profile of dark matter and the stellar mass-to-light ratio. In this work we reduce this degeneracy by combining weak lensing together with strong lensing and stellar kinematics. Our analysis is based on two galaxy samples: 45 strong lenses from the SLACS survey and 1,700 massive quiescent galaxies from the SDSS main spectroscopic sample with weak lensing measurements from the Hyper Suprime-Cam survey. We use a Bayesian hierarchical approach to jointly model all three observables. We fit the data with models of varying complexity and show that a model with a radial gradient in the stellar mass-to-light ratio is required to simultaneously describe both galaxy samples. This result is driven by a subset of strong lenses with very steep total density profile, that cannot be fitted by models with no gradient. Our measurements are unable to determine whether $M_*/L$ gradients are due to variations in stellar population parameters at fixed IMF, or to gradients in the IMF itself. The inclusion of $M_*/L$ gradients decreases dramatically the inferred IMF normalisation, compared to previous lensing-based studies, with the exact value depending on the assumed dark matter profile. The main effect of strong lensing selection is to shift the stellar mass distribution towards the high mass end, while the halo mass and stellar IMF distribution at fixed stellar mass are not significantly affected.

We present results from a pilot observation of nearby (~20 Mpc) galaxies with mass similar to that of the Milky Way (MW) to address the missing satellite problem. This is the first paper from an on-going project to address the problem with a statistical sample of galaxies outside of the Local Group (LG) without employing an assumption that the LG is a typical halo in the Universe. Thanks to the close distances of our targets, dwarf galaxies around them can be identified as extended, diffuse galaxies. By applying a surface brightness cut together with a careful visual screening to remove artifacts and background contamination, we construct a sample of dwarf galaxies. The luminosity function (LF) of one of the targets is broadly consistent with that of the MW, but the other has a more abundant dwarf population. Numerical simulations by Okamoto (2013) seem to overpredict the number of dwarfs on average, while more recent predictions from Copernicus Complexio are in a better agreement. In both observations and simulations, there is a large diversity in the LFs, demonstrating the importance of addressing the missing satellite problem with a statistically representative sample. We also characterize the projected spatial distributions of the satellites and do not observe strong evidence for alignments around the central galaxies. Based on this successful pilot observation, we are carrying out further observations to increase the sample of nearby galaxies, which we plan to report in our future paper.

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

We present 1,201 galaxies at $0.05<z<0.45$ that host tidal features, detected from the first $\sim\! 200$ deg$^2$ of imaging from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). All galaxies in the present sample have spectroscopic observations from the Sloan Digital Sky Survey (SDSS) spectroscopic campaigns, generating a sample of 21208 galaxies. Of these galaxies, we identify 214 shell systems and 987 stream systems. For 575 of these systems, we are additionally able to measure the $(g-i)$ colors of the tidal features. We find evidence for star formation in a subset of the streams, with the exception of streams around massive ellipticals, and find that stream host galaxies span the full range of stellar masses in our sample. Galaxies which host shells are predominantly red and massive: we find that observable shells form more frequently around ellipticals than around disc galaxies of the same stellar mass. Although the majority of the shells in our sample are consistent with being formed by minor mergers, $15\% \pm 4.4\%$ of shell host galaxies have $(g-i)$ colors as red as their host galaxy, consistent with being formed by major mergers. These "red shells" are additionally preferentially aligned with the major axis of the host galaxy, as previously predicted from simulations. We suggest that although the bulk of the observable shell population originates from fairly minor mergers, which preferentially form shells that are not aligned with the major axis of the galaxy, major mergers produce a significant number of observable shells.

We present a catalog of extended low-surface-brightness galaxies (LSBGs) identified in the Wide layer of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). Using the first ${\sim}$200 deg$^2$ of the survey, we have uncovered 781 LSBGs, spanning red ($g-i\geq0.64$) and blue ($g-i<0.64$) colors and a wide range of morphologies. Since we focus on extended galaxies ($r_\mathrm{eff}=2.5$-$14^{\prime\prime}$), our sample is likely dominated by low-redshift objects. We define LSBGs to have mean surface brightnesses $\bar{\mu}_\mathrm{eff}(g)>24.3$ mag arcsec$^{-2}$, which allows nucleated galaxies into our sample. As a result, the central surface brightness distribution spans a wide range of $\mu_0(g)=18$-$27.4$ mag arcsec$^{-2}$, with 50% and 95% of galaxies fainter than 24.3 and 22 mag arcsec$^{-2}$, respectively. Furthermore, the surface brightness distribution is a strong function of color, with the red distribution being much broader and generally fainter than that of the blue LSBGs, and this trend shows a clear correlation with galaxy morphology. Red LSBGs typically have smooth light profiles that are well-characterized by single-component S\'{e}rsic functions. In contrast, blue LSBGs tend to have irregular morphologies and show evidence for ongoing star formation. We crossmatch our sample with existing optical, HI, and ultraviolet catalogs to gain insight into the physical nature of the LSBGs. We find that our sample is diverse, ranging from dwarf spheroidals and ultra-diffuse galaxies in nearby groups to gas-rich irregulars to giant LSB spirals, demonstrating the potential of the HSC-SSP to provide a truly unprecedented view of the LSBG population.

We present the structure of the Milky Way stellar halo beyond Galactocentric distances of $r = 50$ kpc traced by blue horizontal-branch (BHB) stars, which are extracted from the survey data in the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). We select BHB candidates based on $(g,r,i,z)$ photometry, where the $z$-band is on the Paschen series and the colors that involve the $z$-band are sensitive to surface gravity. About 450 BHB candidates are identified between $r = 50$ kpc and 300 kpc, most of which are beyond the reach of previous large surveys including the Sloan Digital Sky Survey. We find that the global structure of the stellar halo in this range has substructures, which are especially remarkable in the GAMA15H and XMM-LSS fields in the HSC-SSP. We find that the stellar halo can be fitted to a single power-law density profile with an index of $\alpha \simeq 3.3$ ($3.5$) with (without) these fields and its global axial ratio is $q \simeq 2.2$ ($1.3$). Thus, the stellar halo may be significantly disturbed and be made in a prolate form by halo substructures, perhaps associated with the Sagittarius stream in its extension beyond $r \sim 100$ kpc. For a broken power-law model allowing different power-law indices inside/outside a break radius, we obtain a steep power-law slope of $\alpha \sim 5$ outside a break radius of $\simeq 100$ kpc ($200$ kpc) for the case with (without) GAMA15H and XMM-LSS. This radius of $200$ kpc might be as close as a halo boundary if there is any, although larger BHB sample is required from further HSC-SSP survey to increase its statistical significance.

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 use ~100 square deg of deep (>28.5 mag arcsec$^{-2}$ in i-band), high-quality (median 0.6 arcsec seeing) imaging data from the Hyper Suprime-Cam (HSC) survey to reveal the halo mass dependence of the surface mass density profiles and outer stellar envelopes of massive galaxies. The i-band images from the HSC survey reach ~4 magnitudes deeper than Sloan Digital Sky Survey and enable us to directly trace stellar mass distributions to 100 kpc without requiring stacking. We conclusively show that, at fixed stellar mass, the stellar profiles of massive galaxies depend on the masses of their dark matter haloes. On average, massive central galaxies with $\log M_{\star, 100\ \mathrm{kpc}}>11.6$ in more massive haloes at 0.3 < z < 0.5 have shallower inner stellar mass density profiles (within ~10-20 kpc) and more prominent outer envelopes. These differences translate into a halo mass dependence of the mass-size relation. Central galaxies in haloes with $\log M_{\rm{Halo}}>14.0$ are ~20% larger in $R_{\mathrm{50}}$ at fixed stellar mass. Such dependence is also reflected in the relationship between the stellar mass within 10 and 100 kpc. Comparing to the mass--size relation, the $\log M_{\star, 100\ \rm{kpc}}$-$\log M_{\star, 10\ \rm{kpc}}$ relation avoids the ambiguity in the definition of size, and can be straightforwardly compared with simulations. Our results demonstrate that, with deep images from HSC, we can quantify the connection between halo mass and the outer stellar halo, which may provide new constraints on the formation and assembly of massive central galaxies.

We present the result of searching for clusters of galaxies based on weak gravitational lensing analysis of the $\sim 160$~deg$^2$ area surveyed by Hyper Suprime-Cam (HSC) as a Subaru Strategic Program. HSC is a new prime focus optical imager with a 1.5 diameter field of view on the 8.2-meter Subaru telescope. The superb median seeing on the HSC $i$-band images of $0.56$ arcsec allows the reconstruction of high angular resolution mass maps via weak lensing, which is crucial for the weak lensing cluster search. We identify 65 mass map peaks with signal-to-noise (SN) ratio larger than 4.7, and carefully examine their properties by cross-matching the clusters with optical and X-ray cluster catalogs. We find that all the 39 peaks with SN$>5.1$ have counterparts in the optical cluster catalogs, and only 2 out of the 65 peaks are probably false positives. The upper limits of X-ray luminosities from ROSAT All Sky Survey (RASS) imply the existence of an X-ray under-luminous cluster population. We show that the X-rays from the shear selected clusters can be statistically detected by stacking the RASS images. The inferred average X-ray luminosity is about half that of the X-ray selected clusters of the same mass. The radial profile of the dark matter distribution derived from the stacking analysis is well modeled by the Navarro-Frenk-White profile with a small concentration parameter value of $c_{500}\sim 2.5$, which suggests that the selection bias on the orientation or the internal structure for our shear selected cluster sample is not strong.

We present optimized source galaxy selection schemes for measuring cluster weak lensing (WL) mass profiles unaffected by cluster member dilution from the Subaru Hyper Suprime-Cam Strategic Survey Program (HSC-SSP). The ongoing HSC-SSP survey will uncover thousands of galaxy clusters to $z\lesssim1.5$. In deriving cluster masses via WL, a critical source of systematics is contamination and dilution of the lensing signal by cluster {members, and by foreground galaxies whose photometric redshifts are biased}. Using the first-year CAMIRA catalog of $\sim$900 clusters with richness larger than 20 found in $\sim$140 deg$^2$ of HSC-SSP data, we devise and compare several source selection methods, including selection in color-color space (CC-cut), and selection of robust photometric redshifts by applying constraints on their cumulative probability distribution function (PDF; P-cut). We examine the dependence of the contamination on the chosen limits adopted for each method. Using the proper limits, these methods give mass profiles with minimal dilution in agreement with one another. We find that not adopting either the CC-cut or P-cut methods results in an underestimation of the total cluster mass ($13\pm4\%$) and the concentration of the profile ($24\pm11\%$). The level of cluster contamination can reach as high as $\sim10\%$ at $R\approx 0.24$ Mpc/$h$ for low-z clusters without cuts, while employing either the P-cut or CC-cut results in cluster contamination consistent with zero to within the 0.5% uncertainties. Our robust methods yield a $\sim60\sigma$ detection of the stacked CAMIRA surface mass density profile, with a mean mass of $M_\mathrm{200c} = (1.67\pm0.05({\rm {stat}}))\times 10^{14}\,M_\odot/h$.

We report fourteen and twenty-eight protocluster candidates at z=5.7 and 6.6 over 14 and 19 deg^2 areas, respectively, selected from 2,230 (259) Lya emitters (LAEs) photometrically (spectroscopically) identified with Subaru/Hyper Suprime-Cam (HSC) deep images (Keck, Subaru, and Magellan spectra and the literature data). Six out of the 42 protocluster candidates include 1-12 spectroscopically confirmed LAEs at redshifts up to z=6.574. By the comparisons with the cosmological Lya radiative transfer (RT) model reproducing LAEs with the reionization effects, we find that more than a half of these protocluster candidates are progenitors of the present-day clusters with a mass of > 10^14 M_sun. We then investigate the correlation between LAE overdensity delta and Lya rest-frame equivalent width EW_Lya^rest, because the cosmological Lya RT model suggests that a slope of EW_Lya^rest-delta relation is steepened towards the epoch of cosmic reionization (EoR), due to the existence of the ionized bubbles around galaxy overdensities easing the escape of Lya emission from the partly neutral intergalactic medium (IGM). The available HSC data suggest that the slope of the EW_Lya^rest-delta correlation does not evolve from the post-reionization epoch z=5.7 to the EoR z=6.6 beyond the moderately large statistical errors. There is a possibility that we would detect the evolution of the EW_Lya^rest - delta relation from z=5.7 to 7.3 by the upcoming HSC observations providing large samples of LAEs at z=6.6-7.3.

Using $\sim$140 deg$^2$ Subaru Hyper Suprime-Cam (HSC) survey data, we stack the weak lensing (WL) signal around five Planck clusters found within the footprint. This yields a 15$\sigma$ detection of the mean Planck cluster mass density profile. The five Planck clusters span a relatively wide mass range, $M_{\rm WL,500c} = (2-30)\times10^{14}\,M_\odot/h$ with a mean mass of $M_{\rm WL,500c} = (4.15\pm0.61)\times10^{14}\,M_\odot/h$. The ratio of the stacked Planck Sunyaev-Zel'dovich (SZ) mass to the stacked WL mass is $ \langle M_{\rm SZ}\rangle/\langle M_{\rm WL}\rangle = 1-b = 0.80\pm0.14$. This mass bias is consistent with previous WL mass calibrations of Planck clusters within the errors. We discuss the implications of our findings for the calibration of SZ cluster counts and the much discussed tension between Planck SZ cluster counts and Planck $\Lambda$CDM cosmology.

We present and characterize the catalog of galaxy shape measurements that will be used for cosmological weak lensing measurements in the Wide layer of the first year of the Hyper Suprime-Cam (HSC) survey. The catalog covers an area of 136.9 deg$^2$ split into six fields, with a mean $i$-band seeing of $0.58$ arcsec and $5\sigma$ point-source depth of $i\sim 26$. Given conservative galaxy selection criteria for first year science, the depth and excellent image quality results in unweighted and weighted source number densities of 24.6 and 21.8 arcmin$^{-2}$, respectively. We define the requirements for cosmological weak lensing science with this catalog, then focus on characterizing potential systematics in the catalog using a series of internal null tests for problems with point-spread function (PSF) modeling, shear estimation, and other aspects of the image processing. We find that the PSF models narrowly meet requirements for weak lensing science with this catalog, with fractional PSF model size residuals of approximately $0.003$ (requirement: 0.004) and the PSF model shape correlation function $\rho_1<3\times 10^{-7}$ (requirement: $4\times 10^{-7}$) at 0.5$^\circ$ scales. A variety of galaxy shape-related null tests are statistically consistent with zero, but star-galaxy shape correlations reveal additive systematics on $>1^\circ$ scales that are sufficiently large as to require mitigation in cosmic shear measurements. Finally, we discuss the dominant systematics and the planned algorithmic changes to reduce them in future data reductions.

We present wide and deep photometry of the northwest part of the halo of the Andromeda galaxy (M31) using Hyper Suprime-Cam on the Subaru Telescope. The survey covers 9.2 deg$^{2}$ field in the $g$, $i$, and $NB515$ bands and shows a clear red giant branch (RGB) of M31's halo stars and a pronounced red clump (RC) feature. The spatial distribution of RC stars shows a prominent stream feature, the North Western (NW) Stream, and a diffuse substructure in the south part of our survey field. We estimate the distances based on the RC method and obtain $(m-M)$ = 24.63$\pm 0.191$(random)$\pm0.057$(systematic) and 24.29$\pm 0.211$(random)$\pm0.057$(systematic) mag for the NW stream and diffuse substructure, respectively, implying that the NW Stream is located behind M31, whereas the diffuse substructure is located in front. We also estimate line-of-sight distances along the NW Stream and find that the south part of the stream is $\sim$20 kpc closer to us relative to the north part. The distance to the NW Stream inferred from the isochrone fitting to the color-magnitude diagram favors the RC-based distance, but the TRGB-based distance estimated for $NB515$-selected RGB stars does not agree with it. The surface number density distribution of RC stars across the NW Stream is found to be approximately Gaussian with a FWHM of $\sim$25 arcmin (5.7 kpc), with a slight skew to the south-west side. That along the NW Stream shows a complicated structure including variations in number density and a significant gap in the stream.

The Hyper Suprime-Cam Subaru Strategic Program (HSC SSP) is an excellent survey for the search for strong lenses, thanks to its area, image quality and depth. We use three different methods to look for lenses among 43,000 luminous red galaxies from the Baryon Oscillation Spectroscopic Survey (BOSS) sample with photometry from the S16A internal data release of the HSC SSP. The first method is a newly developed algorithm, named YATTALENS, which looks for arc-like features around massive galaxies and then estimates the likelihood of an object being a lens by performing a lens model fit. The second method, CHITAH, is a modeling-based algorithm originally developed to look for lensed quasars. The third method makes use of spectroscopic data to look for emission lines from objects at a different redshift from that of the main galaxy. We find 15 definite lenses, 36 highly probable lenses and 282 possible lenses. Among the three methods, YATTALENS, which was developed specifically for this problem, performs best in terms of both completeness and purity. Nevertheless five highly probable lenses were missed by YATTALENS but found by the other two methods, indicating that the three methods are highly complementary. Based on these numbers we expect to find $\sim$300 definite or probable lenses by the end of the HSC SSP.

We present the Ly$\alpha$ luminosity functions (LFs) at $z=$5.7 and 6.6 derived from a new large sample of 1,266 Ly$\alpha$ emitters (LAEs) identified in total areas of 14 and 21 deg$^2$, respectively, based on the early narrowband data of the Subaru/Hyper Suprime-Cam (HSC) survey. Together with careful Monte-Carlo simulations that account for the incompleteness of the LAE selection and the flux estimate systematics in the narrowband imaging, we have determined the Ly$\alpha$ LFs with the unprecedentedly small statistical and systematic uncertainties in a wide Ly$\alpha$ luminosity range of $10^{42.8-43.8}$ erg s$^{-1}$. We obtain the best-fit Schechter parameters of $L^{*}_{{\rm Lya}}=1.6^{+2.2}_{-0.6} (1.7^{+0.3}_{-0.7}) \times10^{43}$ erg s$^{-1}$, $\phi^{*}_{{\rm Lya}}=0.85^{+1.87}_{-0.77}\ (0.47^{+1.44}_{-0.44})\times10^{-4}$ Mpc$^{-3}$, and $\alpha=-2.6^{+0.6}_{-0.4}\ (-2.5^{+0.5}_{-0.5})$ at $z=5.7$ ($6.6$). We confirm that our best-estimate Ly$\alpha$ LFs are consistent with the majority of the previous studies, but find that our Ly$\alpha$ LFs do not agree with the high number densities of LAEs recently claimed by Matthee/Santos et al.'s studies that may overcorrect the incompleteness and the flux systematics. Our Ly$\alpha$ LFs at $z=5.7$ and $6.6$ show an indication that the faint-end slope is very steep ($\alpha \simeq -2.5$), although it is also possible that the bright-end LF results are enhanced by systematic effects such as the contribution from AGNs, blended merging galaxies, and/or large ionized bubbles around bright LAEs. Comparing our Ly$\alpha$ LF measurements with four independent reionization models, we estimate the neutral hydrogen fraction of the IGM to be $x_{\rm HI}=0.3\pm0.2$ at $z=$6.6 that is consistent with the small Thomson scattering optical depth obtained by Planck 2016.

We present a joint X-ray, optical and weak-lensing analysis for X-ray luminous galaxy clusters selected from the MCXC (Meta-Catalog of X-Ray Detected Clusters of Galaxies) cluster catalog in the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) survey field with S16A data, As a pilot study of our planned series papers, we measure hydrostatic equilibrium (H.E.) masses using XMM-Newton data for four clusters in the current coverage area out of a sample of 22 MCXC clusters. We additionally analyze a non-MCXC cluster associated with one MCXC cluster. We show that H.E. masses for the MCXC clusters are correlated with cluster richness from the CAMIRA catalog (Oguri et al. 2017), while that for the non-MCXC cluster deviates from the scaling relation. The mass normalization of the relationship between the cluster richness and H.E. mass is compatible with one inferred by matching CAMIRA cluster abundance with a theoretical halo mass function. The mean gas mass fraction based on H.E. masses for the MCXC clusters is $\langle f_{\rm gas} \rangle = 0.125\pm0.012$ at spherical overdensity $\Delta=500$, which is $\sim80-90$ percent of the cosmic mean baryon fraction, $\Omega_b/\Omega_m$, measured by cosmic microwave background experiments. We find that the mean baryon fraction estimated from X-ray and HSC-SSP optical data is comparable to $\Omega_b/\Omega_m$. A weak-lensing shear catalog of background galaxies, combined with photometric redshifts, is currently available only for three clusters in our sample. Hydrostatic equilibrium masses roughly agree with weak-lensing masses, albeit with large uncertainty. This study demonstrates that further multiwavelength study for a large sample of clusters using X-ray, HSC-SSP optical and weak lensing data will enable us to understand cluster physics and utilize cluster-based cosmology.

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

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

Although Type Ia supernova cosmology has now reached a mature state, it is important to develop as many independent methods as possible to understand the true nature of dark energy. Recent studies have shown that Type II supernovae (SNe II) offer such a path and could be used as alternative distance indicators. However, the majority of these studies were unable to extend the Hubble diagram above redshift $z=0.3$ because of observational limitations. Here, we show that we are now ready to move beyond low redshifts and attempt high-redshift ($z \gtrsim 0.3$) SN~II cosmology as a result of new-generation deep surveys such as the Subaru/Hyper Suprime-Cam (HSC) survey. Applying the "standard candle method" to SN$\sim$2016jhj ($z=0.3398 \pm 0.0002$; discovered by HSC) together with a low-redshift sample, we are able to construct the highest-redshift SN II Hubble diagram to date with an observed dispersion of 0.27 mag (i.e., 12-13% in distance). This work demonstrates the bright future of SN~II cosmology in the coming era of large, wide-field surveys like that of the Large Synoptic Survey Telescope.

We present a statistical study of the redshift evolution of the cluster galaxy population over a wide redshift range from 0.1 to 1.1, using $\sim 1900$ optically-selected CAMIRA clusters from $\sim 232$~deg$^2$ of the Hyper Suprime-Cam (HSC) Wide S16A data. Our stacking technique with a statistical background subtraction reveals color-magnitude diagrams of red-sequence and blue cluster galaxies down to faint magnitudes of $m_z\sim 24$. We find that the linear relation of red-sequence galaxies in the color-magnitude diagram extends down to the faintest magnitudes we explore with a small intrinsic scatter $\sigma_{\rm int}(g-r)<0.1$. The scatter does not evolve significantly with redshift. The stacked color-magnitude diagrams are used to define red and blue galaxies in clusters for studying their radial number density profiles without resorting to photometric redshifts of individual galaxies. We find that red galaxies are significantly more concentrated toward cluster centers and blue galaxies dominate the outskirt of clusters. We explore the fraction of red galaxies in clusters as a function of redshift, and find that the red fraction decreases with increasing distances from cluster centers. The red fraction exhibits a moderate decrease with increasing redshift. The radial number density profiles of cluster member galaxies are also used to infer the location of the steepest slope in the three dimensional galaxy density profiles. For a fixed threshold in richness, we find little redshift evolution in this location.

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

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

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

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

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

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

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

We present the SILVERRUSH program strategy and clustering properties investigated with $\sim 2,000$ Ly$\alpha$ emitters at $z=5.7$ and $6.6$ found in the early data of the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey exploiting the carefully designed narrowband filters. We derive angular correlation functions with the unprecedentedly large samples of LAEs at $z=6-7$ over the large total area of $14-21$ deg$^2$ corresponding to $0.3-0.5$ comoving Gpc$^2$. We obtain the average large-scale bias values of $b_{\rm avg}=4.1\pm 0.2$ ($4.5\pm 0.6$) at $z=5.7$ ($z=6.6$) for $\gtrsim L^*$ LAEs, indicating the weak evolution of LAE clustering from $z=5.7$ to $6.6$. We compare the LAE clustering results with two independent theoretical models that suggest an increase of an LAE clustering signal by the patchy ionized bubbles at the epoch of reionization (EoR), and estimate the neutral hydrogen fraction to be $x_{\rm HI}=0.15^{+0.15}_{-0.15}$ at $z=6.6$. Based on the halo occupation distribution models, we find that the $\gtrsim L^*$ LAEs are hosted by the dark-matter halos with the average mass of $\log (\left < M_{\rm h} \right >/M_\odot) =11.1^{+0.2}_{-0.4}$ ($10.8^{+0.3}_{-0.5}$) at $z=5.7$ ($6.6$) with a Ly$\alpha$ duty cycle of 1 % or less, where the results of $z=6.6$ LAEs may be slightly biased, due to the increase of the clustering signal at the EoR. Our clustering analysis reveals the low-mass nature of $\gtrsim L^*$ LAEs at $z=6-7$, and that these LAEs probably evolve into massive super-$L^*$ galaxies in the present-day universe.

Massive galaxies display extended light profiles that can reach several hundreds of kilo parsecs. These stellar halos provide a fossil record of galaxy assembly histories. Using data that is both wide (~100 square degree) and deep (i>28.5 mag/arcsec^2 in i-band), we present a systematic study of the stellar halos of a sample of more than 3000 galaxies at 0.3 < z < 0.5 with $\log M_{\star}/M_{\odot} > 11.4$. Our study is based on high-quality (0.6 arcsec seeing) imaging data from the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP), which enables us to individually estimate surface mass density profiles to 100 kpc without stacking. As in previous work, we find that more massive galaxies exhibit more extended outer profiles. When this extended light is not properly accounted for as a result of shallow imaging or inadequate profile modeling, the derived stellar mass function can be significantly underestimated at the highest masses. Across our sample, the ellipticity of outer light profiles increases substantially as we probe larger radii. We show for the first time that these ellipticity gradients steepen dramatically as a function of galaxy mass, but we detect no mass-dependence in outer color gradients. Our results support the two-phase formation scenario for massive galaxies in which outer envelopes are built up at late times from a series of merging events. We provide surface mass surface mass density profiles in a convenient tabulated format to facilitate comparisons with predictions from numerical simulations of galaxy formation.

We present wide-field (167 deg$^2$) weak lensing mass maps from the Hyper Supreme-Cam Subaru Strategic Program (HSC-SSP). We compare these weak lensing based dark matter maps with maps of the distribution of the stellar mass associated with luminous red galaxies. We find a strong correlation between these two maps with a correlation coefficient of $\rho=0.54\pm0.03$ (for a smoothing size of $8'$). This correlation is detected even with a smaller smoothing scale of $2'$ ($\rho=0.34\pm 0.01$). This detection is made uniquely possible because of the high source density of the HSC-SSP weak lensing survey ($\bar{n}\sim 25$ arcmin$^{-2}$). We also present a variety of tests to demonstrate that our maps are not significantly affected by systematic effects. By using the photometric redshift information associated with source galaxies, we reconstruct a three-dimensional mass map. This three-dimensional mass map is also found to correlate with the three-dimensional galaxy mass map. Cross-correlation tests presented in this paper demonstrate that the HSC-SSP weak lensing mass maps are ready for further science analyses.

We present an optically-selected cluster catalog from the Hyper Suprime-Cam (HSC) Subaru Strategic Program. The HSC images are sufficiently deep to detect cluster member galaxies down to $M_*\sim 10^{10.2}M_\odot$ even at $z\sim 1$, allowing a reliable cluster detection at such high redshifts. We apply the CAMIRA algorithm to the HSC Wide S16A dataset covering $\sim 232$ deg$^2$ to construct a catalog of 1921 clusters at redshift $0.1<z<1.1$ and richness $\hat{N}_{\rm mem}>15$ that roughly corresponds to $M_{\rm 200m}\gtrsim 10^{14}h^{-1}M_\odot$. We confirm good cluster photometric redshift performance, with the bias and scatter in $\Delta z/(1+z)$ being better than 0.005 and 0.01 over most of the redshift range, respectively. We compare our cluster catalog with large X-ray cluster catalogs from XXL and XMM-LSS surveys and find good correlation between richness and X-ray properties. We also study the miscentering effect from the distribution of offsets between optical and X-ray cluster centers. We confirm the high ($>0.9$) completeness and purity for high mass clusters by analyzing mock galaxy catalogs.

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

We present the procedure to build and validate the bright-star masks for the Hyper-Suprime-Cam Strategic Subaru Proposal (HSC-SSP) survey. To identify and mask the saturated stars in the full HSC-SSP footprint, we rely on the Gaia and Tycho-2 star catalogues. We first assemble a pure star catalogue down to $G_{\rm Gaia} < 18$ after removing $\sim1.5\%$ of sources that appear extended in the Sloan Digital Sky Survey (SDSS). We perform visual inspection on the early data from the S16A internal release of HSC-SSP, finding that our star catalogue is $99.2\%$ pure down to $G_{\rm Gaia} < 18$. Second, we build the mask regions in an automated way using stacked detected source measurements around bright stars binned per $G_{\rm Gaia}$ magnitude. Finally, we validate those masks from visual inspection and comparison with the literature of galaxy number counts and angular two-point correlation functions. This version (Arcturus) supersedes the previous version (Sirius) used in the S16A internal and DR1 public releases. We publicly release the full masks and tools to flag objects in the entire footprint of the planned HSC-SSP observations at this address: ftp://obsftp.unige.ch/pub/coupon/brightStarMasks/HSC-SSP/.

We present the deepest optical images of the COSMOS field based on a joint dataset taken with Hyper Suprime-Cam (HSC) by the HSC Subaru Strategic Program (SSP) team and the University of Hawaii (UH). The COSMOS field is one of the key extragalactic fields with a wealth of deep, multi-wavelength data. However, the current optical data are not sufficiently deep to match with, e.g., the UltraVista data in the near-infrared. The SSP team and UH have joined forces to produce very deep optical images of the COSMOS field by combining data from both teams. The coadd images reach depths of g=27.8, r=27.7, i=27.6, z=26.8, and y=26.2 mag at 5 sigma for point sources based on flux uncertainties quoted by the pipeline and they cover essentially the entire COSMOS 2 square degree field. The seeing is between 0.6 and 0.9 arcsec on the coadds. We perform several quality checks and confirm that the data are of science quality; ~2% photometry and 30 mas astrometry. This accuracy is identical to the Public Data Release 1 from HSC-SSP. We make the joint dataset including fully calibrated catalogs of detected objects available to the community at https://hsc-release.mtk.nao.ac.jp/.

We present the results from a search for new Milky Way (MW) satellites from the first two years of data from the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP) $\sim 300$~deg$^2$ and report the discovery of a highly compelling ultra-faint dwarf galaxy candidate in Cetus. This is the second ultra-faint dwarf we have discovered after Virgo~I reported in our previous paper. This satellite, Cetus~III, has been identified as a statistically significant (10.7$\sigma$) spatial overdensity of star-like objects, which are selected from a relevant isochrone filter designed for a metal-poor and old stellar population. This stellar system is located at a heliocentric distance of 251$^{+24}_{-11}$~kpc with a most likely absolute magnitude of $M_V = -2.4 \pm 0.6$~mag estimated from a Monte Carlo analysis. Cetus~III is extended with a half-light radius of $r_h = 90^{+42}_{-17}$~pc, suggesting that this is a faint dwarf satellite in the MW located beyond the detection limit of the Sloan Digital Sky Survey. Further spectroscopic studies are needed to assess the nature of this stellar system. We also revisit and update the parameters for Virgo~I finding $M_V = -0.33^{+0.75}_{-0.87}$~mag and $r_h = 47^{+19}_{-13}$~pc. Using simulations of $\Lambda$-dominated cold dark matter models, we predict that we should find one or two new MW satellites from $\sim 300$~deg$^2$ HSC-SSP data, in rough agreement with the discovery rate so far. The further survey and completion of HSC-SSP over $\sim 1,400$~deg$^2$ will provide robust insights into the missing satellites problem.

The Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) is currently the deepest wide- field survey in progress. The 8.2 m aperture of Subaru telescope is very powerful in detect- ing faint/small moving objects, including near-Earth objects, asteroids, centaurs and Tran- Neptunian objects (TNOs). However, the cadence and dithering pattern of the HSC-SSP are not designed for detecting moving objects, making it difficult to do so systematically. In this paper, we introduce a new pipeline for detecting moving objects (specifically TNOs) in a non-dedicated survey. The HSC-SSP catalogs are re-arranged into the HEALPix architecture. Then, the stationary detections and false positive are removed with a machine learning al- gorithm to produce a list of moving object candidates. An orbit linking algorithm and visual inspections are executed to generate the final list of detected TNOs. The preliminary results of a search for TNOs using this new pipeline on data from the first HSC-SSP data release (Mar 2014 to Nov 2015) are also presented.

We present a list of galaxy-scale lens candidates including a highly probable interacting galaxy-scale lens in the Hyper Suprime-Cam (HSC) imaging survey. We combine HSC imaging with the blended-spectra catalog from the Galaxy And Mass Assembly (GAMA) survey to identify lens candidates, and use lens mass modeling to confirm the candidates. We find 46 matches between the HSC S14A_0b imaging data release and the GAMA catalog. Ten of them are probable lens systems according to their morphology and redshifts. There is one system with an interacting galaxy pair, HSC J084928+000949, that has a valid mass model. We predict the total mass enclosed by the Einstein radius of $\sim0.72$" ($\sim1.65$kpc) for this new expected lens system to be $\sim10^{10.59}M_{\odot}$. Using the photometry in the {\it grizy} bands of the HSC survey and stellar population synthesis modeling with a Salpeter stellar initial mass function, we estimate the stellar mass within the Einstein radius to be $\sim10^{10.46}\,M_{\odot}$. We thus find a dark matter mass fraction within the Einstein radius of $\sim25\%$. Further spectroscopy or high-resolution imaging would allow confirmation of the nature of these lens candidates. The particular system with the interacting galaxy pair, if confirmed, would provide an opportunity to study the interplay between dark matter and stars as galaxies build up through hierarchical mergers.

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

We report the discovery of a diffuse stellar cloud with an angular extent $\gtrsim30^{\prime\prime}$, which we term "Sumo Puff", in data from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). While we do not have a redshift for this object, it is in close angular proximity to a post-merger galaxy at redshift $z=0.0431$ and is projected within a few virial radii (assuming similar redshifts) of two other ${\sim}L_\star$ galaxies, which we use to bracket a potential redshift range of $0.0055 < z < 0.0431$. The object's light distribution is flat, as characterized by a low Sersic index ($n\sim0.3$). It has a low central $g$-band surface brightness of ${\sim}26.4$ mag arcsec$^{-2}$, large effective radius of ${\sim}13^{\prime\prime}$ (${\sim}11$ kpc at $z=0.0431$ and ${\sim}1.5$ kpc at $z=0.0055$), and an elongated morphology ($b/a\sim0.4$). Its red color ($g-i\sim1$) is consistent with a passively evolving stellar population and similar to the nearby post-merger galaxy, and we may see tidal material connecting Sumo Puff with this galaxy. We offer two possible interpretations for the nature of this object: (1) it is an extreme, galaxy-size tidal feature associated with a recent merger event, or (2) it is a foreground dwarf galaxy with properties consistent with a quenched, disturbed ultra-diffuse galaxy. We present a qualitative comparison with simulations that demonstrates the feasibility of forming a structure similar to this object in a merger event. Follow-up spectroscopy and/or deeper imaging to confirm the presence of the bridge of tidal material will be necessary to reveal the true nature of this object.

We utilize the HSC CAMIRA cluster catalog and the photo-$z$ galaxy catalog constructed in the HSC wide field (S16A), covering $\sim$ 174 deg$^{2}$, to study the star formation activity of galaxies in different environments over 0.2 $<$ $z$ $<$ 1.1. We probe galaxies down to $i \sim$ 26, corresponding to a stellar mass limit of log$_{10}$(M$_*$/M$_{\odot}$) $\sim$ 8.2 and $\sim$ 8.6 for star-forming and quiescent populations, respectively, at $z$ $\sim$ 0.2. The existence of the red sequence for low stellar mass galaxies in clusters suggests that the environmental quenching persists to halt the star formation in the low-mass regime. In addition, star-forming galaxies in groups or clusters are systematically biased toward lower values of specific star formation rate by 0.1 -- 0.3 dex with respect to those in the field and the offsets shows no strong redshift evolution over our redshift range, implying a universal slow quenching mechanism acting in the dense environments since $z$ $\sim$ 1.1. Moreover, the environmental quenching dominates the mass quenching in low mass galaxies, and the quenching dominance reverses in high mass ones. The transition mass is greater in clusters than in groups, indicating that the environmental quenching is more effective for massive galaxies in clusters compared to groups.

We present a visible multi-band photometry of trans-Neptunian objects (TNOs) observed by the Subaru Telescope in the framework of Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) from March in 2014 to September in 2016. We measured the five broad-band (g, r, i, z, and Y) colors over the wavelength range from 0.4 um to 1.0 um for 30 known TNOs using the HSC-SSP survey data covering ~500 deg2 of sky within +/-30 deg of ecliptic latitude. This dataset allows us to characterize the dynamical classes based on visible reflectance spectra as well as to examine the relationship between colors and the other parameters such as orbital elements. Our results show that the hot classical and scattered populations share similar color distributions, while the cold classical population has a reflective decrease toward shorter wavelength below the i band. Based on the obtained color properties, we found that the TNO sample examined in the present work can be separated into two groups by inclination (I), the low-I population consisting of cold classical objects and high-I population consisting of hot classical and scattered objects. The whole sample exhibits an anti-correlation between colors and inclination, but no significant correlation between colors and semi-major axis, perihelion distance, eccentricity, or absolute magnitude. The color-inclination correlation does not seem to be continuous over the entire inclination range. Rather, it is seen only in the high-I population. We found that the low- and high-I populations are distinguishable in the g-i vs. eccentricity plot, but four high-I objects show g-i colors similar to those of the low-I population. If we exclude these four objects, the high-I objects show a positive correlation between g-i and eccentricity and a negative correlation between g-i and inclination with high significance levels.

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

We present the first results of the Subaru/Hyper Suprime-Cam (HSC) survey of the interacting galaxy system, NGC4631 and NGC4656. From the maps of resolved stellar populations, we identify 11 dwarf galaxies (including already-known dwarfs) in the outer region of NGC4631 and the two tidal stellar streams around NGC4631, named Stream SE and Stream NW, respectively. This paper describes the fundamental properties of these tidal streams. Based on the tip of red giant branch method and the Bayesian statistics, we find that StreamSE (7.10 Mpc in Expected a posteriori, EAP, with the 90% credible intervals of [6.22, 7.29] Mpc) and StreamNW (7.91 Mpc in EAP with the 90% credible intervals of [6.44, 7.97] Mpc) are located in front of and behind NGC4631, respectively. We also calculate the metallicity distribution of stellar streams by comparing the member stars with theoretical isochrones on the color-magnitude diagram. We find that both streams have the same stellar population based on the Bayesian model selection method, suggesting that they originated from a tidal interaction between NGC4631 and a single dwarf satellite. The expected progenitor has a positively skewed metallicity distribution function with [M/H]_EAP=-0.92 with the 90% credible intervals of [-1.46, -0.51]. The stellar mass of the progenitor is estimated as 3.7 x 10e+8 Msun with the 90% credible intervals of [5.8 x 10e+6, 8.6 x 10e+9] Msun based on the mass-metallicity relation for Local group dwarf galaxies. This is in good agreement with an initial stellar mass of the progenitor presumed in the previous N-body simulation.

We surveyed a central 0.6 deg^2 region of Abell 1367 cluster for extended ionized gas clouds (EIGs) using the Subaru prime-focus camera (Suprime-Cam) with a narrow-band filter that covers Halpha. We discovered six new EIGs in addition to five known EIGs. We also found that the Halpha tail from the blue infalling group (BIG) is extended to about 330 kpc in projected distance, which is about twice longer than previously reported. Candidates of star-forming blobs in the tail are detected. The properties of the EIG parent galaxies in Abell 1367 basically resemble those in the Coma cluster. A noticeable difference is that the number of detached EIGs is significantly fewer in Abell 1367, while the fraction of blue member galaxies is higher. The results suggest a difference in the evolutionary stage of the clusters; Abell 1367 is at an earlier stage than the Coma cluster.

We report the serendipitous discovery of a quadruply lensed source at $z_{\rm s}=3.76$, HSC~J115252+004733, from the Hyper Suprime-Cam (HSC) Survey. The source is lensed by an early-type galaxy at $z_{\rm l}=0.466$ and a satellite galaxy. Here, we investigate the properties of the source by studying its size and luminosity from the imaging and the luminosity and velocity width of the Ly-$\alpha$ line from the spectrum. Our analyses suggest that the source is most probably a low-luminosity active galactic nucleus (LLAGN) but the possibility of it being a compact bright galaxy (e.g., a Lyman-$\alpha$ emitter or Lyman Break Galaxy) cannot be excluded. The brighter pair of lensed images appears point-like except in the HSC $i$-band (with a seeing $\sim0.5"$). The extended emission in the $i$-band image could be due to the host galaxy underneath the AGN, or alternatively, due to a highly compact lensed galaxy (without AGN) which appears point-like in all bands except in $i$-band. We also find that the flux ratio of the brighter pair of images is different in the Ks-band compared to optical wavelengths. Phenomena such as differential extinction and intrinsic variability cannot explain this chromatic variation. While microlensing from stars in the foreground galaxy is less likely to be the cause, it cannot be ruled out completely. If the galaxy hosts an AGN, then this represents the highest redshift quadruply imaged AGN known to date, enabling study of a distant LLAGN. Discovery of this unusually compact and faint source demonstrates the potential of the HSC survey.

We present the deep and wide $V$ and $I_c$ photometry of the Sextans dwarf spheroidal galaxy (dSph) taken by Suprime-Cam imager on the Subaru Telescope, which extends out to the tidal radius. The colour-magnitude diagram (CMD) reaches two magnitudes below the main sequence (MS) turn-off, showing a steep red giant branch, blue and red horizontal branch (HB), sub-giant branch (SGB), MS, and blue stragglers (BS). We construct the radial profile of each evolutionary phase and demonstrate that blue HB stars are more spatially extended, while red HB stars are more centrally concentrated than the other components. The colour distribution of SGB stars also varies with the galactocentric distance; the inner SGB stars shift bluer than those in the outskirt. The radial differences in the CMD morphology indicate the existence of the age gradient. The relatively younger stars ($\sim10$ Gyr) are more centrally concentrated than the older ones ($\sim13$ Gyr). The spatial contour maps of stars in different age bins also show that the younger population has higher concentration and higher ellipticity than the older one. We also detect the centrally concentrated bright BS stars, the number of which is consistent with the idea that a part of these stars belongs to the remnant of a disrupted star cluster discovered in the previous spectroscopic studies.

We present measurements of the clustering properties of a sample of infrared (IR) bright dust-obscured galaxies (DOGs). Combining 125 deg$^2$ of wide and deep optical images obtained with the Hyper Suprime-Cam on the Subaru Telescope and all-sky mid-IR (MIR) images taken with Wide-Field Infrared Survey Explorer, we have discovered 4,367 IR-bright DOGs with $(i - [22])_{\rm AB}$ $>$ 7.0 and flux density at 22 $\mu$m $>$ 1.0 mJy. We calculate the angular autocorrelation function (ACF) for a uniform subsample of 1411 DOGs with 3.0 mJy $<$ flux (22 $mu$m) $<$ 5.0 mJy and $i_{\rm AB}$ $<$ 24.0. The ACF of our DOG subsample is well-fit with a single power-law, $\omega (\theta)$ = (0.010 $\pm$ 0.003) $\theta^{-0.9}$, where $\theta$ in degrees. The correlation amplitude of IR-bright DOGs is larger than that of IR-faint DOGs, which reflects a flux-dependence of the DOG clustering, as suggested by Brodwin et al. (2008). We assume that the redshift distribution for our DOG sample is Gaussian, and consider 2 cases: (1) the redshift distribution is the same as IR-faint DOGs with flux at 22 $\mu$m $<$ 1.0 mJy, mean and sigma $z$ = 1.99 $\pm$ 0.45, and (2) $z$ = 1.19 $\pm$ 0.30, as inferred from their photometric redshifts. The inferred correlation length of IR-bright DOGs is $r_0$ = 12.0 $\pm$ 2.0 and 10.3 $\pm$ 1.7 $h^{-1}$ Mpc, respectively. IR-bright DOGs reside in massive dark matter halos with a mass of $\log [\langle M_{\mathrm{h}} \rangle / (h^{-1} M_{\odot})]$ = 13.57$_{-0.55}^{+0.50}$ and 13.65$_{-0.52}^{+0.45}$ in the two cases, respectively.

We describe a weak lensing view of the downsizing of star forming galaxies based on cross correlating a weak lensing ($\kappa$) map with a predicted map constructed from a redshift survey. Moderately deep and high resolution images with Subaru/Hyper Suprime-Cam covering the 4 deg^2 DLS F2 field provide a $\kappa$ map with 1 arcmin resolution. A dense complete redshift survey of the F2 field including 12,705 galaxies with $R\leq20.6$ is the basis for construction of the predicted map. The zero-lag cross-correlation between the \kappa and predicted maps is significant at the $30\sigma$ level. The width of the cross-correlation peak is comparable with the angular scale of rich cluster at $z\sim0.3$, the median depth of the redshift survey. Slices of the predicted map in $\delta{z} = 0.05$ redshift bins enable exploration of the impact of structure as a function of redshift. The zero-lag normalised cross-correlation has significant local maxima at redshifts coinciding with known massive X-ray clusters. Even in slices where there are no known massive clusters, there is significant signal in the cross-correlation originating from lower mass groups that trace the large-scale of the universe. Spectroscopic $D_n4000$ measurements enable division of the sample into star-forming and quiescent populations. The significance of the cross-correlation with structure containing star-forming galaxies increases with redshift from $5\sigma$ at $z = 0.3$ to $7 \sigma$ at $z = 0.5$. The weak lensing results are consistent with the downsizing view of galaxy evolution established on the basis of many other independent studies.

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

We report the discovery of 15 quasars and bright galaxies at 5.7 < z < 6.9. This is the initial result from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the exquisite multiband imaging data produced by the Subaru Hyper Suprime-Cam (HSC) Strategic Program survey. The candidate selection is performed by combining several photometric approaches including a Bayesian probabilistic algorithm to reject stars and dwarfs. The spectroscopic identification was carried out with the Gran Telescopio Canarias and the Subaru Telescope for the first 80 deg2 of the survey footprint. The success rate of our photometric selection is quite high, approaching 100 % at the brighter magnitudes (zAB < 23.5 mag). Our selection also recovered all the known high-z quasars on the HSC images. Among the 15 discovered objects, six are likely quasars, while the other six with interstellar absorption lines and in some cases narrow emission lines are likely bright Lyman-break galaxies. The remaining three objects have weak continua and very strong and narrow Ly alpha lines, which may be excited by ultraviolet light from both young stars and quasars. These results indicate that we are starting to see the steep rise of the luminosity function of z > 6 galaxies, compared with that of quasars, at magnitudes fainter than M1450 ~ -22 mag or zAB ~24 mag. Follow-up studies of the discovered objects as well as further survey observations are ongoing.

We report the serendipitous discovery of HSC J142449-005322, a double source plane lens system in the Hyper Suprime-Cam Subaru Strategic Program. We dub the system Eye of Horus. The lens galaxy is a very massive early-type galaxy with stellar mass of ~7x10^11 Msun located at z_L=0.795. The system exhibits two arcs/rings with clearly different colors, including several knots. We have performed spectroscopic follow-up observations of the system with FIRE on Magellan. The outer ring is confirmed at z_S2=1.988 with multiple emission lines, while the inner arc and counterimage is confirmed at z_S1=1.302. This makes it the first double source plane system with spectroscopic redshifts of both sources. Interestingly, redshifts of two of the knots embedded in the outer ring are found to be offset by delta_z=0.002 from the other knots, suggesting that the outer ring consists of at least two distinct components in the source plane. We perform lens modeling with two independent codes and successfully reproduce the main features of the system. However, two of the lensed sources separated by ~0.7 arcsec cannot be reproduced by a smooth potential, and the addition of substructure to the lens potential is required to reproduce them. Higher-resolution imaging of the system will help decipher the origin of this lensing feature and potentially detect the substructure.

We present clustering analysis results from 10,381 Lyman break galaxies (LBGs) at z~ 4-7, identified in the Hubble legacy deep imaging and new complimentary large-area Subaru/Hyper Suprime-Cam data. We measure the angular correlation functions (ACFs) of these LBGs at z~4, 5, 6, and 7, and fit these measurements using halo occupation distribution (HOD) models that provide an estimate of halo masses, M_h~(1-20)x10^11 Msun. Our M_h estimates agree with those obtained by previous clustering studies in a UV-magnitude vs. M_h plane, and allow us to calculate stellar-to-halo mass ratios (SHMRs) of LBGs. By comparison with the z~0 SHMR, we identify evolution of the SHMR from z~0 to z~4, and z~4 to z~7 at the >98% confidence levels. The SHMR decreases by a factor of ~2 from z~0 to 4, and increases by a factor of ~4 from z~4 to 7. We compare our SHMRs with results of a hydrodynamic simulation and a semi-analytic model, and find that these theoretical studies do not predict the SHMR increase from z~4 to 7. We obtain the baryon conversion efficiency (BCE) of LBGs at z~4, and find that the BCE increases with increasing dark matter halo mass. Finally, we compare our clustering+HOD estimates with results from abundance matching techniques, and conclude that the M_h estimates of the clustering+HOD analyses agree with those of the simple abundance matching within a factor of 3, and that the agreement improves when using more sophisticated abundance matching techniques that include subhalos, incompleteness, and/or evolution in the star formation and stellar mass functions.

We revealed the detailed structure of a vastly extended H$\alpha$-emitting nebula ("H$\alpha$ nebula") surrounding the starburst/merging galaxy NGC 6240 by deep narrow-band imaging observations with the Subaru Suprime-Cam. The extent of the nebula is $\sim$90 kpc in diameter and the total H$\alpha$ luminosity amounts to $\approx 1.6 \times 10^{42}$ erg s$^{-1}$. The volume filling factor and the mass of the warm ionized gas are $\sim$10$^{-4}$--10$^{-5}$ and $\sim$$5 \times 10^8$ $M_\odot$, respectively. The nebula has a complicated structure, which includes numerous filaments, loops, bubbles, and knots. We found that there is a tight spatial correlation between the H$\alpha$ nebula and the extended soft X-ray-emitting gas, both in large and small scales. The overall morphology of the nebula is dominated by filamentary structures radially extending from the center of the galaxy. A large-scale bi-polar bubble extends along the minor axis of the main stellar disk. The morphology strongly suggests that the nebula was formed by intense outflows - superwinds - driven by starbursts. We also found three bright knots embedded in a looped filament of ionized gas that show head-tail morphologies in both emission-line and continuum, suggesting close interactions between the outflows and star forming regions. Based on the morphology and surface brightness distribution of the H$\alpha$ nebula, we propose the scenario that three major episodes of starburst/superwind activities which were initiated $\sim$10$^2$ Myr ago formed the extended ionized gas nebula of NGC 6240.

We present both the observations and the data reduction procedures of the Subaru COSMOS 20 project that is an optical imaging survey of the HST COSMOS field, carried out by using Suprime-Cam on the Subaru Telescope with the following 20 optical filters: 6 broad-band (B, g', V, r', i', and z'), 2 narrow-band (NB711 and NB816), and 12 intermediate-band filters (IA427, IA464, IA484, IA505, IA527, IA574, IA624, IA679, IA709, IA738, IA767, and IA827). A part of this project is described in Taniguchi et al. (2007) and Capak et al. (2007) for the six broad-band and one narrow-band (NB816) filter data. In this paper, we present details of the observations and data reduction for remaining 13 filters (the 12 IA filters and NB711). In particular, we describe the accuracy of both photometry and astrometry in all the filter bands. We also present optical properties of the Suprime-Cam IA filter system in Appendix.

We report the discovery of 854 ultra diffuse galaxies (UDGs) in the Coma cluster using deep R band images, with partial B, i, and Halpha band coverage, obtained with the Subaru telescope. Many of them (332) are Milky Way-sized with very large effective radii of r_e>1.5kpc. This study was motivated by the recent discovery of 47 UDGs by van-Dokkum et al. (2015); our discovery suggests >1,000 UDGs after accounting for the smaller Subaru field. The new UDGs show a distribution concentrated around the cluster center, strongly suggesting that the great majority are (likely longtime) cluster members. They are a passively evolving population, lying along the red sequence in the CM diagram with no Halpha signature. Star formation was, therefore, quenched in the past. They have exponential light profiles, effective radii re ~ 800 pc- 5 kpc, effective surface brightnesses mu_e(R)=25-28 mag arcsec-2, and stellar masses ~1x10^7 - 5x10^8Msun. There is also a population of nucleated UDGs. Some MW-sized UDGs appear closer to the cluster center than previously reported; their survival in the strong tidal field, despite their large sizes, possibly indicates a large dark matter fraction protecting the diffuse stellar component. The indicated baryon fraction ~<1% is less than the cosmic average, and thus the gas must have been removed from the possibly massive dark halo. The UDG population appears to be elevated in the Coma cluster compared to the field, indicating that the gas removal mechanism is related primarily to the cluster environment.

We present the photometric properties of a sample of infrared (IR) bright dust obscured galaxies (DOGs). Combining wide and deep optical images obtained with the Hyper Suprime-Cam (HSC) on the Subaru Telescope and all-sky mid-IR (MIR) images taken with Wide-Field Infrared Survey Explorer (WISE), we discovered 48 DOGs with $i - K_\mathrm{s} > 1.2$ and $i - [22] > 7.0$, where $i$, $K_\mathrm{s}$, and [22] represent AB magnitude in the $i$-band, $K_\mathrm{s}$-band, and 22 $\mu$m, respectively, in the GAMA 14hr field ($\sim$ 9 deg$^2$). Among these objects, 31 ($\sim$ 65 %) show power-law spectral energy distributions (SEDs) in the near-IR (NIR) and MIR regime, while the remainder show a NIR bump in their SEDs. Assuming that the redshift distribution for our DOGs sample is Gaussian, with mean and sigma $z$ = 1.99 $\pm$ 0.45, we calculated their total IR luminosity using an empirical relation between 22 $\mu$m luminosity and total IR luminosity. The average value of the total IR luminosity is (3.5 $\pm$ 1.1) $\times$ $10^{13}$ L$_{\odot}$, which classifies them as hyper-luminous infrared galaxies (HyLIRGs). We also derived the total IR luminosity function (LF) and IR luminosity density (LD) for a flux-limited subsample of 18 DOGs with 22 $\mu$m flux greater than 3.0 mJy and with $i$-band magnitude brighter than 24 AB magnitude. The derived space density for this subsample is log $\phi$ = -6.59 $\pm$ 0.11 [Mpc$^{-3}$]. The IR LF for DOGs including data obtained from the literature is well fitted by a double-power law. The derived lower limit for the IR LD for our sample is $\rho_{\mathrm{IR}}$ $\sim$ 3.8 $\times$ 10$^7$ [L$_{\odot}$ Mpc$^{-3}$] and its contributions to the total IR LD, IR LD of all ultra-luminous infrared galaxies (ULIRGs), and that of all DOGs are $>$ 3 %, $>$ 9 %, and $>$ 15 %, respectively.

We present properties of moderately massive clusters of galaxies detected by the newly developed Hyper Suprime-Cam on the Subaru telescope using weak gravitational lensing. Eight peaks exceeding a S/N ratio of 4.5 are identified on the convergence S/N map of a 2.3 square degree field observed during the early commissioning phase of the camera. Multi-color photometric data is used to generate optically selected clusters using the CAMIRA algorithm. The optical cluster positions were correlated with the peak positions from the convergence map. All eight significant peaks have optical counterparts. The velocity dispersion of clusters are evaluated by adopting the Singular Isothemal Sphere (SIS) fit to the tangential shear profiles, yielding virial mass estimates, M500c, of the clusters which range from 2.7x10^13 to 4.4x10^14 solar mass. The number of peaks is considerably larger than the average number expected from LambdaCDM cosmology but this is not extremely unlikely if one takes the large sample variance in the small field into account. We could, however, safely argue that the peak count strongly favours the recent Planck result suggesting high sigma8$value of 0.83. The ratio of stellar mass to the dark matter halo mass shows a clear decline as the halo mass increases. If the gas mass fraction, fg, in halos is universal, as has been suggested in the literature, the observed baryon mass in stars and gas shows a possible deficit compared with the total baryon density estimated from the baryon oscillation peaks in anisotropy of the cosmic microwave background.

We report the discovery of a new dwarf galaxy (NGC6503-d1) during the Subaru extended ultraviolet (XUV) disk survey. It is a likely companion of the spiral galaxy NGC6503. The resolved images, in B, V, R, i, and Halpha, show an irregular appearance due to bright stars with underlying, smooth and unresolved stellar emission. It is classified as the transition type (dIrr/dSph). Its structural properties are similar to those of the dwarfs in the Local Group, with a V absolute magnitude ~ -10.5, half-light radius ~400 pc, and central surface brightness ~25.2. Despite the low stellar surface brightness environment, one HII region was detected, though its Halpha luminosity is low, indicating an absence of any appreciable O-stars at the current epoch. The presence of multiple stellar populations is indicated by the color-magnitude diagram of ~300 bright resolved stars and the total colors of the dwarf, with the majority of its total stellar mass ~4x10^6 Msun in an old stellar population.

We present photometric recalibration of the Subaru Deep Field (SDF) and Subaru/XMM-Newton Deep Survey (SXDS). Recently, Yamanoi et al. (2012) suggested the existence of a discrepancy between the SDF and SXDS catalogs. We have used the Sloan Digital Sky Survey (SDSS) Data Release 8 (DR8) catalog and compared stars in common between SDF/SXDS and SDSS. We confirmed that there exists a 0.12 mag offset in B-band between the SDF and SXDS catalogs. Moreover, we found that significant zero point offsets in i-band (~ 0.10 mag) and z-band (~ 0.14 mag) need to be introduced to the SDF/SXDS catalogs to make it consistent with the SDSS catalog. We report the measured zero point offsets of five filter bands of SDF/SXDS catalogs. We studied the potential cause of these offsets, but the origins are yet to be understood.

The Andromeda galaxy (M31) shows many tidal features in its halo, including the Giant Southern Stream (GSS) and a sharp ledge in surface density on its western side (the W Shelf). Using DEIMOS on the Keck telescope, we obtain radial velocities of M31's giant stars along its NW minor axis, in a radial range covering the W Shelf and extending beyond its edge. In the space of velocity versus radius, the sample shows the wedge pattern expected from a radial shell, which is detected clearly here for the first time. This confirms predictions from an earlier model of formation of the GSS, which proposed that the W Shelf is a shell from the third orbital wrap of the same tidal debris stream that produces the GSS, with the main body of the progenitor lying in the second wrap. We calculate the distortions in the shelf wedge pattern expected from its outward expansion and angular momentum, and show that these effects are echoed in the data. In addition, a hot, relatively smooth spheroid population is clearly present. We construct a bulge-disk-halo N-body model that agrees with surface brightness and kinematic constraints, and combine it with a simulation of the GSS. From the contrasting kinematic signatures of the hot spheroid and shelf components, we decompose the observed stellar metallicity distribution into contributions from each component using a non-parametric mixture model. The shelf component's metallicity distribution matches previous observations of the GSS superbly, further strengthening the evidence they are connected and bolstering the case for a massive progenitor of this stream.

We derived the luminosity function (LF) of dwarf galaxies in the Coma Cluster down to M_R=-10 at three fields located at the center, intermediate, and outskirt. The LF (-19<M_R<-10) shows no significant differences among the three fields. It shows a clear dip at M_R\sim-13, and is composed of two distinct components of different slopes; the bright component with -19<M_R<-13 has a flatter slope than the faint component with -13<M_R<-10 which has a steep slope. The bright component (-19<M_R<-13) consists of mostly red extended galaxies including few blue galaxies whose colors are typical of late-type galaxies. On the other hand, the faint component (-13<M_R<-10) consists of largely PSF-like compact galaxies. We found that both these compact galaxies and some extended galaxies are present in the center while only compact galaxies are seen in the outskirt. In the faint component, the fraction of blue galaxies is larger in the outskirt than in the center. We suggest that the dwarf galaxies in the Coma Cluster, which make up the two components in the LF, are heterogeneous with some different origins.

From H\alpha narrow band observations, we identified three H\alpha emitting regions in the direction of Magellanic Stream IV (MS IV). They consist of three parallel filaments with 2 arcmin width and 6 -- 30 arcmin length at 12 arcmin intervals. The mean surface brightness of them is $\sim 2 \times 10^{-18}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$. Because of their low surface brightness, the regions were not detected in previous H\alpha surveys. In HI map, the position of the filaments overlap MS, suggesting that they are parts of MS, but there also exists a local HI structure. If the filaments associate with MS, the sizes are 30 pc $\times$ 100 -- 500 pc. The filaments lie at the leading edge of a downstream cloud, which supports a shock heating and its propagation (shock cascade) model for the ionizing source. If they are local objects, on the other hand, Fossil Str\"omgren Trails of more than two stars is a possible interpretation, and the sizes would be 0.1 pc $\times$ 0.3 -- 1.5 pc at 180 pc distance. The positional information of the H\alpha filaments presented in this letter enables us future spectroscopic observations to clarify their nature.

We present the results of deep imaging and spectroscopic observations of very extended ionized gas (EIG) around four member galaxies of the Coma cluster of galaxies: RB199, IC4040, GMP2923 and GMP3071. The EIGs were serendipitously found in an H-alpha narrow band imaging survey of the central region of the Coma cluster. The relative radial velocities of the EIGs with respect to the systemic velocities of the parent galaxies from which they emanate increase almost monotonically with the distance from the nucleus of the respective galaxies, reaching -400 - -800 km/s at around 40 - 80 kpc from the galaxies. The one-sided morphologies and the velocity fields of the EIGs are consistent with the predictions of numerical simulations of ram pressure stripping. We found a very low-velocity filament (v_rel = -1300 km/s) at the southeastern edge of the disk of IC4040. Some bright compact knots in the EIGs of RB199 and IC4040 exhibit blue continuum and strong H-alpha emission. The equivalent widths of the H-alpha emission exceed 200 A, and are greater than 1000 A for some knots. The emission line intensity ratios of the knots are basically consistent with those of sub-solar abundance HII regions. These facts indicate that intensive star formation occurs in the knots. Some filaments, including the low velocity filament of the IC4040 EIG, exhibit shock-like emission line spectra, suggesting that shock heating plays an important role in ionization and excitation of the EIGs.

As part of our survey of galactic stellar halos, we investigate the structure and stellar populations of the northern outer part of the stellar halo in NGC55, a member galaxy of the Sculptor Group, using deep and wide-field V- and I-band images taken with Subaru/Suprime-Cam. Based on the analysis of the color-magnitude diagrams (CMDs) for red-giant-branch (RGB) stars, we derive a tip of RGB (TRGB)-based distance modulus to the galaxy of (m-M)_0 = 26.58 +/- 0.11 (d = 2.1 +/- 0.1 Mpc). From the stellar density maps, we detect the asymmetrically disturbed, thick disk structure and two metal-poor overdense substructures in the north region of NGC55, which may correspond to merger remnants associated with hierarchical formation of NGC55's halo. In addition, we identify a diffuse metal-poor halo extended out to at least z ~ 16 kpc from the galactic plane. The surface-brightness profiles toward the z-direction perpendicular to the galactic plane suggest that the stellar density distribution in the northern outer part of NGC55 is described by a locally isothermal disk at z <~ 6 kpc and a likely diffuse metal-poor halo with V-band surface brightness of mu_V >~ 32 mag arcsec^{-2}, where old RGB stars dominate. We derive the metallicity distributions (MDs) of these structures on the basis of the photometric comparison of RGB stars with the theoretical stellar evolutionary models. The MDs of the thick disk structures show the peak and mean metallicity of [Fe/H]peak ~ -1.4 and [Fe/H]mean ~ -1.7, respectively, while the outer substructures show more metal-poor features than the thick disk structure. Combined with the current results with our previous study for M31's halo, we discuss the possible difference in the formation process of stellar halos among different Hubble types.

The Coma cluster was the target of a HST-ACS Treasury program designed for deep imaging in the F475W and F814W passbands. Although our survey was interrupted by the ACS instrument failure in 2007, the partially completed survey still covers ~50% of the core high-density region in Coma. Observations were performed for 25 fields that extend over a wide range of cluster-centric radii (~1.75 Mpc) with a total coverage area of 274 arcmin^2. The majority of the fields are located near the core region of Coma (19/25 pointings) with six additional fields in the south-west region of the cluster. In this paper we present reprocessed images and SExtractor source catalogs for our survey fields, including a detailed description of the methodology used for object detection and photometry, the subtraction of bright galaxies to measure faint underlying objects, and the use of simulations to assess the photometric accuracy and completeness of our catalogs. We also use simulations to perform aperture corrections for the SExtractor Kron magnitudes based only on the measured source flux and half-light radius. We have performed photometry for ~73,000 unique objects; one-half of our detections are brighter than the 10-sigma point-source detection limit at F814W=25.8 mag (AB). The slight majority of objects (60%) are unresolved or only marginally resolved by ACS. We estimate that Coma members are 5-10% of all source detections, which consist of a large population of unresolved objects (primarily GCs but also UCDs) and a wide variety of extended galaxies from a cD galaxy to dwarf LSB galaxies. The red sequence of Coma member galaxies has a constant slope and dispersion across 9 magnitudes (-21<M_F814W<-13). The initial data release for the HST-ACS Coma Treasury program was made available to the public in 2008 August. The images and catalogs described in this study relate to our second data release.

We present images of extended H-alpha clouds associated with 14 member galaxies in the Coma cluster obtained from deep narrow band imaging observations with Suprime-Cam at the Subaru Telescope. The parent galaxies of the extended H-alpha clouds are distributed farther than 0.2 Mpc from the peak of X-ray emission of the cluster. Most of the galaxies have colors bluer than g-r approx 0.5 and they account for 57% of the blue (g-r<0.5) bright (r<17.8 mag) galaxies in the central region of the Coma cluster. They reside near the red- and blue-shifted edges of the Coma cluster's radial velocity distribution. These findings suggest that the most of the parent galaxies were recently captured by the Coma cluster potential and are now infalling toward the cluster center with their disk gas being stripped off and producing the observed H-alpha clouds.

QSOs have been thought to be important for tracing highly biased regions in the early universe, from which the present-day massive galaxies and galaxy clusters formed. While overdensities of star-forming galaxies have been found around QSOs at 2<z<5, the case for excess galaxy clustering around QSOs at z>6 is less clear. Previous studies with HST have reported the detection of small excesses of faint dropout galaxies in some QSO fields, but these surveys probed a relatively small region surrounding the QSOs. To overcome this problem, we have observed the most distant QSO at z=6.4 using the large field of view of the Suprime-Cam (34' x 27'). Newly-installed CCDs allowed us to select Lyman break galaxies (LBG) at z~6.4 more efficiently. We found seven LBGs in the QSO field, whereas only one exists in a comparison field. The significance of this apparent excess is difficult to quantify without spectroscopic confirmation and additional control fields. The Poisson probability to find seven objects when one expects four is ~10%, while the probability to find seven objects in one field and only one in the other is less than 0.4%, suggesting that the QSO field is significantly overdense relative to the control field. We find some evidence that the LBGs are distributed in a ring-like shape centered on the QSO with a radius of ~3 Mpc. There are no candidate LBGs within 2 Mpc from the QSO, i.e., galaxies are clustered around the QSO but appear to avoid the very center. These results suggest that the QSO is embedded in an overdense region when defined on a sufficiently large scale. This suggests that the QSO was indeed born in a massive halo. The central deficit of galaxies may indicate that (1) the strong UV radiation from the QSO suppressed galaxy formation in its vicinity, or (2) that star-formation closest to the QSO occurs mostly in an obscured mode that is missed by our UV selection.

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

We search for stars with proper motions in a set of deep Subaru images, covering about 0.48 square degrees to a depth of $i' \simeq 26$, taken over a span of five and a half years. We follow the methods described in \citet{Richmond2009} to reduce and analyze this dataset. We present a sample of 69 stars with motions of high significance, and discuss briefly the populations from which they are likely drawn. Based on photometry and motions alone, we expect that 14 of the candidates may be white dwarfs. Our candidate with the largest proper motion is surprisingly faint and likely to prove interesting: its colors and motions suggest that it might be an M dwarf moving at over 500 km/sec or an L dwarf in the halo.

We present a photometric survey of the stellar halo of the Andromeda galaxy, using Suprime-Cam on the Subaru Telescope. A detailed analysis of VI color-magnitude diagrams of the resolved stellar population is used to measure properties such as line-of-sight distance, surface brightness, metallicity, and age, and these are used to isolate and characterize different components of the M31 halo: (1) several substructures, and (2) the smooth halo. First, we study M31's halo substructure along the north-west/south-east minor axis out to R ~ 100 kpc and the south-west major axis region at R ~ 60 kpc. We confirm two substructures in the south-east halo reported by Ibata et al. (2007) and discover two overdense substructures in the north-west halo. We investigate the properties of these four substructures as well as other structures including the western shelf and find that differences in stellar populations among these systems, thereby suggesting each has a different origin. Our statistical analysis implies that the M31 halo as a whole may contain at least 16 substructures, each with a different origin. Second, we investigate the properties of an underlying, smooth and extended halo component out to R > 100 kpc. We find that the surface density of this smooth halo can be fitted to a Hernquist model of scale radius ~ 17 kpc or a power-law profile with ~ R^{-2.17 +/- 0.15}. In contrast to the relative smoothness of the halo density profile, its metallicity distribution appears to be spatially non-uniform with non-monotonic variations with radius, suggesting that the halo population has not had sufficient time to dynamically homogenize the accreted populations. Further implications for the formation of the M31 halo are discussed.

We report an optical detection of an extended structure around a QSO at z=6.43 (CFHQSJ2329-0301, the highest redshift QSO currently known) in deep z' and z_r-band images of the Subaru/Suprime-Cam. After a careful PSF (QSO) subtraction, a structure in the z'-band extends more than 4'' on the sky (R_e=11 kpc), and thus, is well-resolved (16sigma detection). The PSF-subtracted z_r-band structure is in a similar shape to that in the z'-band, but less significant with a 3 sigma detection. In the z'-band, a radial profile of the QSO+host shows a clear excess over that of the averaged PSF in 0.8-3'' radius. Since the z'-band includes a Lya emission at z=6.43, the z' flux is perhaps a mixture of the host (continuum light) and its Lya emission, whereas the z_r-band flux is from the host. Through a SED modeling, we estimate 40% of the PSF-subtracted z'-band light is from the host (continuum) and 60% is from Lya emission. The absolute magnitude of the host is M_{1450}=-23.9 (c.f. M_{1450}=-26.4 for the QSO). A lower limit of the SFR(Lya) is 1.6 Msun yr^{-1} with stellar mass ranging 6.2 x 10^8 to 1.1 x 10^10 Msun when 100 Myrs of age is assumed. The detection shows that a luminous QSO is already harbored by a large, star-forming galaxy in the early Universe only after ~840 Myr after the big bang. The host may be a forming giant galaxy, co-evolving with a super massive black hole.

We introduce the Dichroic-Mirror Camera (DMC), an instrument developed at the University of Tokyo which is capable of performing simultaneous imaging in fifteen bands. The main feature of the DMC is the dichroic mirrors, which split incoming light into red and blue components. Combination of dichroic mirrors split light from the telescope focus into fifteen intermediate-width bands across 390 -- 950 nm. The fifteen bands of DMC provide measurements of the object's spectral energy distribution (SED) at fifteen wavelength points. During May -- June 2007 and March 2008, observing run of the DMC was carried out at Higashi-Hiroshima Observatory, Japan. We observed several objects i.e. planets, asteroids, standard stars & star clusters, planetary nebulae, galaxies, and supernovae. We describe several early scientific results from the DMC.

We search for stars with proper motions in a set of twenty deep Subaru images, covering about 0.28 square degrees to a depth of i ~ 25, taken over a span of six years. In this paper, we describe in detail our reduction and techniques to identify moving objects. We present a first sample of 99 stars with motions of high significance, and discuss briefly the populations from which they are likely drawn. Based on photometry and motions alone, we expect that 9 of the candidates may be white dwarfs. We also find a group of stars which may be extremely metal-poor subdwarfs in the halo.

Using the HST ACS, we have obtained deep optical images reaching well below the oldest main sequence turnoff in fields on the southeast minor-axis of the Andromeda Galaxy, 35 kpc from the nucleus. These data probe the star formation history in the extended halo of Andromeda -- that region beyond 30 kpc that appears both chemically and morphologically distinct from the metal-rich, highly-disturbed inner spheroid. The present data, together with our previous data for fields at 11 and 21 kpc, do not show a simple trend toward older ages and lower metallicities, as one might expect for populations further removed from the obvious disturbances of the inner spheroid. Specifically, the mean ages and [Fe/H] values at 11 kpc, 21 kpc, and 35 kpc are 9.7 Gyr and -0.65, 11.0 Gyr and -0.87, and 10.5 Gyr and -0.98, respectively. In the best-fit model of the 35 kpc population, one third of the stars are younger than 10 Gyr, while only ~10% of the stars are truly ancient and metal-poor. The extended halo thus exhibits clear evidence of its hierarchical assembly, and the contribution from any classical halo formed via early monolithic collapse must be small.

We found an unusual complex of narrow blue filaments, bright blue knots, and H-alpha emitting filaments and clouds, which morphologically resembled a complex of ``fireballs,'' extending up to 80 kpc south from an E+A galaxy RB199 in the Coma cluster. The galaxy has a highly disturbed morphology indicative of a galaxy--galaxy merger remnant. The narrow blue filaments extend in straight shapes toward the south from the galaxy, and several bright blue knots are located at the southern ends of the filaments. The Rc band absolute magnitudes, half light radii and estimated masses of the bright knots are -12 - -13 mag, 200 - 300 pc and 10^6-7 Msolar, respectively. Long, narrow H-alpha emitting filaments are connected at the south edge of the knots. The average color of the fireballs is B - Rc = 0.5, which is bluer than RB199 (B - R = 0.99), suggesting that most of the stars in the fireballs were formed within several times 10^8 yr. The narrow blue filaments exhibit almost no H-alpha emission. Strong H-alpha and UV emission appear in the bright knots. These characteristics indicate that star formation recently ceased in the blue filaments and now continues in the bright knots. The gas stripped by some mechanism from the disk of RB199 may be traveling in the intergalactic space, forming stars left along its trajectory. The most plausible fireball formation mechanism is ram pressure stripping by high-speed collision between the galaxy and the hot intra-cluster medium. The fireballs may be a snapshot of diffuse intra-cluster population formation, or halo star population formation in a cluster galaxy.

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

We describe the HST ACS Coma cluster Treasury survey, a deep two-passband imaging survey of one of the nearest rich clusters of galaxies, the Coma cluster (Abell 1656). The survey was designed to cover an area of 740 square arcmin in regions of different density of both galaxies and intergalactic medium within the cluster. The ACS failure of January 27th 2007 leaves the survey 28% complete, with 21 ACS pointings (230 square arcmin) complete, and partial data for a further 4 pointings (44 square arcmin). Predicted survey depth for 10 sigma detections for optimal photometry of point sources is g' = 27.6 in the F475W filter, and IC=26.8 mag in F814 (AB magnitudes). Initial simulations with artificially injected point sources show 90% recovered at magnitude limits of g' = 27.55 and IC = 26.65. For extended sources, the predicted 10 sigma limits for a 1 square arcsecond region are g' = 25.8 mag/sq. arcsec and IC = 25.0 mag/sq. arcsec. We highlight several motivating science goals of the survey, including study of the faint end of the cluster galaxy luminosity function, structural parameters of dwarf galaxies, stellar populations and their effect on colors and color gradients, evolution of morphological components in a dense environment, the nature of ultra compact dwarf galaxies, and globular cluster populations of cluster galaxies of a range of luminosities and types. This survey will also provide a local rich cluster benchmark for various well known global scaling relations and explore new relations pertaining to the nuclear properties of galaxies.

We report our first results of photometric and spectroscopic observations for Tycho's supernova remnant (SNR Tycho) to search for the companion star of a type Ia supernova (SN Ia). From photometric observations using Suprime-Cam on the Subaru Telescope, we have picked up stars brighter than 22 mag (in $V$-band) for spectroscopy, which are located within a circular region with the radius of 30" around the center of SNR Tycho. If the ejecta of young supernova remnants, such as SNR Tycho, have a sufficient amount of Fe I, we should be able to detect absorption lines at 3720 \AA and 3860 \AA associated with transitions from the ground state of Fe I in the spectrum of the companion star. To identify the companion star of a SN Ia using these characteristic absorption lines of Fe I, we made optical low-resolution spectroscopy of their targets using FOCAS on the Subaru Telescope. In our spectroscopic observations, we obtained spectra of 17 stars in the SNR Tycho region and compare them with template stellar spectra. We detect significant absorption lines from two stars at 3720 \AA. Since widths of their absorption lines are broad, it is likely that the detected absorptions are due to Fe I in the expanding ejecta of SNR Tycho. However, none of stars exhibits a clear red wing in the observed profiles of the absorption, though a star in the background of the SNR should show it. Hence, we suggest another interpretation that the detected absorption lines might be caused by the peculiarity of stars. A star named Tycho(E) has the absorption line at 3720 \AA and its projected position is close to the center of SNR Tycho. Based on our observations, Tycho(E) is a new candidate as the companion star of Tycho's supernova.

We carried out a wide-field V, I imaging survey of the Local Group dwarf spheroidal galaxy Leo II using the Subaru Prime Focus Camera on the 8.2-m Subaru Telescope. The survey covered an area of 26.67 x 26.67 arcmin^2, far beyond the tidal radius of Leo II (8.63 arcmin), down to the limiting magnitude of V ~26, which is roughly 1 mag deeper than the turn-off point of the main sequence stars of Leo II. Radial number density profiles of bright and faint red giant branch (RGB) stars were found to change their slopes at around the tidal radius, and extend beyond the tidal radius with shallower slopes. A smoothed surface brightness map of Leo II suggests the existence of a small substructure of globular cluster luminosity beyond the tidal radius. We investigated the properties of the stellar population by means of the color-magnitude diagram. The horizontal branch (HB) morphology index shows a radial gradient in which red HB stars are more concentrated than blue HB stars, which is common to many Local Group dwarf spheroidal galaxies. The color distribution of RGB stars around the mean RGB sequence shows a larger dispersion at the center than in the outskirts, indicating a mixture of stellar populations at the center and a more homogeneous population in the outskirts. Based on the age estimation using subgiant branch (SGB) stars, we found that although the major star formation took place ~8 Gyr ago, a considerable stellar population younger than 8 Gyr is found at the center; such a younger population is insignificant in the outskirts.

We report the serendipitous findings of 13 faint meteors and 44 artificial space objects by Subaru SuprimeCam imaging observations during 11-16 August 2004. The meteors, at about 100km altitude, and artificial satellites/debris in orbit, at 500km altitude or higher, were clearly discriminated by their apparent defocused image sizes. CCD photometry of the 13 meteors, including 1 Perseid, 1 Aquarid, and 11 sporadic meteors, was performed. We defined a peak video-rate magnitude by comparing the integrated photon counts from the brightest portion of the track traversed within 33ms to those from a 0-mag star during the same time duration. This definition gives magnitudes in the range 4.0< V_{vr} <6.4 and 4.1< I_{vr}<5.9 for these 13 meteors. The corresponding magnitude for virtual naked-eye observers could be somewhat fainter especially for the V-band observation, in which the [OI] 5577 line lasting about 1 sec as an afterglow could contribute to the integrated flux of the present 5-10 min CCD exposures. Although the spatial resolution is insufficient to resolve the source size of anything smaller than about 1 m, we developed a new estimate of the collisionally excited column diameter of these meteors. A diameter as small as a few mm was derived from their collisionally excited photon rates, meteor speed, and the volume density of the oxygen atoms at the 100km altitude. The actual column diameter of the radiating zone, however, could be as large as few 100m because the excited atoms travel that distance before they emit forbidden lines in 0.7 sec of its average lifetime. Among the 44 artificial space objects, we confirmed that 17 were cataloged satellites/space debris.

We have used Suprime-Cam on the Subaru Telescope to conduct a V- and I-band imaging survey of fields sampling the spheroid of the Andromeda galaxy along its south-east minor axis. Our photometric data are deep enough to resolve stars down to the red clump. Based on a large and reliable sample of red giant stars available from this deep wide-field imager, we have derived metallicity distributions vs. radius and a surface brightness profile over projected distances of R=23-66 kpc from the galaxy's center. The metallicity distributions across this region shows a clear high mean metallicity and a broad distribution ([Fe/H] ~ -0.6 +/- 0.5), and indicates no metallicity gradient within our observed range. The surface brightness profile at R>40 kpc is found to be flatter than previously thought. It is conceivable that this part of the halo samples as yet unidentified, metal-rich substructure.

We present imaging data and photometry for the COSMOS survey in 15 photometric bands between 0.3um and 2.4um. These include data taken on the Subaru 8.3m telescope, the KPNO and CTIO 4m telescopes, and the CFHT 3.6m telescope. Special techniques are used to ensure that the relative photometric calibration is better than 1% across the field of view. The absolute photometric accuracy from standard star measurements is found to be 6%. The absolute calibration is corrected using galaxy spectra, providing colors accurate to 2% or better. Stellar and galaxy colors and counts agree well with the expected values. Finally, as the first step in the scientific analysis of these data we construct panchromatic number counts which confirm that both the geometry of the universe and the galaxy population are evolving.

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.

We present results from a narrow-band optical survey of a contiguous area of 1.95 deg^2, covered by the Cosmic Evolution Survey (COSMOS). Both optical narrow-band (lambda_c = 8150 AA and Delta_lambda = 120 AA) and broad-band (B, V, g', r', i', and z') imaging observations were performed with the Subaru prime-focus camera, Suprime-Cam on the Subaru Telescope. We provide the largest contiguous narrow-band survey, targetting Ly alpha emitters (LAEs) at z~5.7. We find a total of 119 LAE candidates at z~5.7. Over the wide-area covered by this survey, we find no strong evidence for large scale clustering of LAEs. We estimate a star formation rate (SFR) density of ~7*10^-4 M_sun yr^-1 Mpc^-3 for LAEs at z~5.7, and compare it with previous measurements.

We study the galaxy population in the central region and a region about 0.6 Mpc away from the center of the Hydra I cluster in B- and Rc-bands down to M~-10 using the Subaru Suprime-Cam photometry. We find that the luminosity function of the entire population has a slightly steeper slope (alpha~-1.6) in the range of -20<M<-10 than those reported for other clusters in slightly brighter ranges. The slope appears to be steeper in poorer clusters. The number of faint galaxies (M>-14) increases in the cluster center as well as bright galaxies. The Hydra I cluster is dominated by red dwarfs and the luminosity function shows a significant upturn at M~-16 as is seen in several other nearby clusters, but not in the field. This upturn and the variation in the faint-end slope of the luminosity function may be caused by the cluster environment or the evolution history of individual clusters.

In the deep narrow band image of the Coma Cluster taken with Suprime-Cam of the Subaru telescope, we found an extremely long and narrow (~ 60 kpc x 2 kpc) Halpha emitting region associated with a poststarburst galaxy (D100). Follow up spectroscopy shows that the region has the same redshift as D100. The surface brightness of the region is typically 25 mag(AB) arcsec^-2 in Halpha, which corresponds to 0.5--4x10^-17 erg s^-1 cm^-2 arcsec^-2. We set two possible explanations for the origin of the region; gas stripped off from a merged dwarf, or gas stripped off from D100 by ram pressure. Either scenario has difficulty to fully explain all the observed characteristics of the region.

We report on the discovery of a new potential galaxy threshing system in the COSMOS 2 square degree field using the prime-focus camera, Suprime-Cam, on the 8.2 m Subaru Telescope. This system consists of a giant elliptical galaxy with $M_V \approx -21.6$ and a tidally disrupted satellite galaxy with $M_V \approx -17.7$ at a photometric redshift of $z \approx 0.08$. This redshift is consistent with the spectroscopic redshift of 0.079 for the giant elliptical galaxy obtained from the Sloan Digital Sky Survey (SDSS) archive. The luminosity masses of the two galaxies are $3.7 \times 10^{12} \cal{M}_{\odot}$ and $3.1 \times 10^{9} \cal{M}_{\odot}$, respectively. The distance between the two galaxies is greater than 100 kpc. The two tidal tails emanating from the satellite galaxy extend over 150 kpc. This system would be the second well-defined galaxy threshing system found so far.

We present deep optical imaging observations of 2 square degree area, covered by the Cosmic Evolution Survey (COSMOS), made by the prime-focus Camera (Supreme-Cam) on the 8.2m Subaru Telescope. Observations were done in six broad-band [B (4459.7 AA), g' (4723.1 AA), V (5483.8 AA), r' (6213.0 AA), i' (7640.8 AA), z' (8855.0 AA)], and one narrow-band (NB816) filters. A total of 10^6 galaxies were detected to i'~26.5 mag. These data, combined with observations at u* and K-band are used to construct the photometric catalogs for the COSMOS and to measure their photometric redshifts, multi-band spectral energy distributions, stellar masses and identification of high redshift candidates. This catalog provides multi-waveband data for scientific analysis of the COSMOS survey.

Gamma-ray bursts (GRBs) and their afterglows have been proposed as an excellent probe to study the evolution of cosmic star formation, the reionization of the intergalactic medium, and the metal enrichment history of the universe, since the prompt gamma-ray emission of GRBs should be detectable out to distances z>10. Hitherto, the highest measured redshift for a GRB has been z=4.50. Here we report the optical spectrum of the afterglow of GRB 050904 obtained 3.4 days after the burst. The spectrum shows a clear continuum at the long wavelength end of the spectrum with a sharp cutoff at around 9000 A due to Ly alpha absorption at a redshift of 6.3 with a damping wing. Little flux is present in the waveband shortward of the Ly alpha break. A system of absorption lines of heavy elements at redshift z=6.295 +- 0.002 were also detected, yielding a precise measurement of the largest known redshift of a GRB. Analysis of the Si II fine structure lines suggest a dense metal-enriched environment around the GRB progenitor, providing unique information on the properties of the gas in a galaxy when the universe was younger than one billion years.

The Cosmic Evolution Survey (COSMOS) is a Hubble Space Telescope (HST) treasury project.The COSMOS aims to perform a 2 square degree imaging survey of an equatorial field in $I$(F814W) band, using the Advanced Camera for Surveys (ACS). Such a wide field survey, combined with ground-based photometric and spectroscopic data, is essential to understand the interplay between large scale structure, evolution and formation of galaxies and dark matter. In 2004, we have obtained high-quality, broad band images of the COSMOS field ($B, V, r^\prime, i^\prime,$ and $ z^\prime$) using Suprime-Cam on the Subaru Telescope, and we have started our new optical multi-band program, COSMOS-21 in 2005. Here, we present a brief summary of the current status of the COSMOS project together with contributions from the Subaru Telescope. Our future Subaru program, COSMOS-21, is also discussed briefly.

We report the discovery of a probable large scale structure composed of many galaxy clumps around the known twin clusters at z=1.26 and z=1.27 in the Lynx region. Our analysis is based on deep, panoramic, and multi-colour imaging (26.4'x24.1') in VRi'z'-bands with the Suprime-Cam on the 8.2m Subaru telescope. This unique, deep and wide-field imaging data-set allows us for the first time to map out the galaxy distribution in the highest redshift supercluster known. We apply a photometric redshift technique to extract plausible cluster members at z~1.27 down to i'=26.15 (5sigma) corresponding to \~M*+2.5 at this redshift. From the 2-D distribution of these photometrically selected galaxies, we newly identify seven candidates of galaxy groups or clusters where the surface density of red galaxies is significantly high (>5sigma), in addition to the two known clusters. These candidates show clear red colour-magnitude sequences consistent with a passive evolution model, which suggests the existence of additional high density regions around the Lynx superclusters.

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

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

We present the results of a survey for Ly-alpha emitters at z ~ 5.8 using a new intermediate-band filter centered at lambda_c = 8275 AA with Delta-lambda_FWHM = 340 AA (i.e., the spectroscopic resolution is R ~ 23) with a combination with a traditional narrow-band centered at lambda_c = 8150 AA with Delta-lambda_FWHM = 120 AA (R ~ 68). Our observations were made with use of the Subaru Prime Focus Camera, Suprime-Cam, on the 8.2 m Subaru telescope in a sky area surrounding the high redshift quasar, SDSSp J104433.04-012522.2 at z=5.74, covering an effective sky area with ~ 720 arcmin^2. In this survey, we have found four Ly-alpha-emitter candidates from the intermediate-band image (z \~ 5.8 with Delta z ~ 0.3). Combined with our previous results based on the NB816 imaging, we discuss the star formation activity in galaxies between z ~ 5.7 and z ~ 5.9.

A main topic at this meeting is how galaxies are affected when they enter for the first time the cluster environment from the outskirts. Most of the times we are forced to infer the environmental effects indirectly, relying on systematic variations of galaxy properties with environment, but there aren't many examples of direct observations able to unveil ongoing transformations taking place, and the corresponding mechanism producing it. We present a case in which it is possible to identify the cluster environment, and in particular the intracluster medium and the recent infall history of galaxies onto the cluster, as the cause for a recent, abrupt change in the evolutionary history of galaxies.

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 have conducted a deep multi-color imaging survey of 0.2 degrees^2 centered on the Hubble Deep Field North (HDF-N). We shall refer to this region as the Hawaii-HDF-N. Deep data were collected in U, B, V, R, I, and z' bands over the central 0.2 degrees^2 and in HK' over a smaller region covering the Chandra Deep Field North (CDF-N). The data were reduced to have accurate relative photometry and astrometry across the entire field to facilitate photometric redshifts and spectroscopic followup. We have compiled a catalog of 48,858 objects in the central 0.2 degrees^2 detected at 5 sigma significance in a 3" aperture in either R or z' band. Number counts and color-magnitude diagrams are presented and shown to be consistent with previous observations. Using color selection we have measured the density of objects at 3<z<7. Our multi-color data indicates that samples selected at z>5.5 using the Lyman break technique suffer from more contamination by low redshift objects than suggested by previous studies.

We report results of a deep wide-field narrowband survey for redshift z~5.7 Ly alpha emitters carried out with SuprimeCam on Subaru 8.3-m telescope. Deep narrowband imaging of the SSA22 field through a 120 A bandpass filter centered at 8150 A was combined with deep multicolor RIz' SuprimeCam broadband imaging, and BVRZ imaging taken with CFHT's CFH12K camera to select high-redshift galaxy candidates. Spectroscopic observations were made using the new wide-field multi-object DEIMOS spectrograph on Keck for 22 of the 26 candidate objects. Eighteen objects were identified as z~5.7 Lyman alpha emitters, and a further nineteenth candidate was identified based on an LRIS spectrum. At the 3.3 A resolution of the DEIMOS spectra the asymmetric profile for Ly alpha emission with its steep blue fall-off can be clearly seen. We use this to describe the distribution of equivalent widths and the continuum color break properties for z~5.7 Ly alpha galaxies compared with foreground objects. The large majority (>75%) of Ly alpha lines have rest frame equivalent widths less than 240 A and can be understood in terms of young star forming galaxies with a Salpeter initial mass function for the stars. With narrowband selection criteria of I-N > 0.7 and N<25.05 (AB mags) we find a surface density of Ly alpha emitters of 0.03 per square arcminute per (deltaz=0.1) to a limiting flux just under 2e-17 erg/cm2/s. The luminosity function of the Ly alpha emitters is similar to that at lower redshifts to the lowest measurable luminosity of 1e43 ergs/s as is the universal star formation rate based on their continuum properties. We note that the objects are highly structured in both spatial and spectral properties on the angular scale of the fields (~60 Mpc), and that multiple fields will have to be averaged to accurately measure their ensemble properties.

A satisfactory understanding of the origin of the dependence of galaxy properties on their environment has remained, so far, out of reach. In the light of numerous observational results and substantial theoretical progress obtained for clusters of galaxies in the last years, a primary goal is to understand how the star formation activity depends on cluster substructure, i.e. on the merging/accretion history of a cluster. In this contribution we present a case in which it is possible to identify the cluster environment, and in particular the intracluster medium and the recent infall history of galaxies onto the cluster, as the cause for an abrupt change in the star formation histories of a subset of galaxies in the Coma cluster.

A surprising result of our recent spectroscopic survey of galaxies in the Coma cluster has been the discovery of a possible bimodal distribution in the metallicities of faint galaxies at $M_B>-17$. We identified a group of dwarfs with luminosity-weighted metallicities around solar and a group with [M/H] around -1.5. A metallicity bimodality among galaxies of similar luminosities is unexpected and suggests that faint cluster galaxies could be an heterogeneous population that formed through more than one evolutionary path, possibly as a consequence of the cluster environment.

We present an optical multicolor-imaging study of the galaxy cluster Abell 521 at $z = 0.25$, using Suprime-Cam on the Subaru Telescope, covering an area of $32 \times 20$ arcmin$^2$ ($9.4 \times 5.8 h_{50}^{-2}$ Mpc$^2$ at $z = 0.25$). Our imaging data taken with both a narrow-band filter, $NB816$ ($\lambda_0 = 8150$\AA and $\Delta \lambda = 120$\AA), and broad-band filters, $B,V,R_{\rm C}, i^\prime$, and $z^\prime$ allow us to find 165 H$\alpha$ emitters. We obtain the H$\alpha$ luminosity function (LF) for the cluster galaxies within 2 Mpc; the Schechter parameters are $\alpha = -0.75 \pm 0.23$, $\phi^\star = 10^{-0.25 \pm 0.20}$ Mpc$^{-3}$, and $L^\star = 10^{42.03 \pm 0.17}$ erg s$^{-1}$. Although the faint end slope, $\alpha$, is consistent with that of the local cluster H$\alpha$ LFs, the characteristic luminosity, $L^\star$, is about 6 times (or $\approx 2$ mag) brighter. This strong evolution implies that Abell 521 contains more active star-forming galaxies than the local clusters, being consistent with the observed Butcher-Oemler effect. However, the bright $L^\star$ of Abell 521 may be, at least in part, due to the dynamical condition of this cluster.

We report here the results of deep optical spectroscopy of the very extended emission-line region (VEELR) found serendipitously around the Seyfert 2 galaxy NGC 4388 in the Virgo cluster using the Subaru Telescope. The H-alpha recession velocities of most of the filaments of the region observed are highly blue-shifted with respect to the systemic velocity of the galaxy. The velocity field is complicated, and from the kinematic and morphological points of view, there seem to be several streams of filaments: low velocity (v = -100 km/s) filaments, high velocity (v = -300 km/s) filaments, and a very high velocity (v > -500 km/s) cloud complex. The emission-line ratios of the VEELR filaments are well explained by power-law photoionization models with solar abundances, suggesting that the Seyfert nucleus of NGC 4388 is the dominant ionization source of the VEELR and that the VEELR gas has moderate metallicity. In addition to photoionization, shock heating probably contributes to the ionization of the gas. In particular, the filaments outside the ionization cone of the Seyfert nucleus are mainly excited by shocks. We conclude that the VEELR was formerly the disk gas of NGC 4388, which has been stripped by ram pressure due to the interaction between the hot intra-cluster medium (ICM) and the galaxy. The velocity field and the morphology of the VEELR closely resemble snapshots from some numerical simulations of this process. In the case of NGC 4388, the ram pressure-stripped gas, which is normally seen as extended HI filaments, happens to be exposed and ionized by the radiation from the AGN, and so can be seen as optical emission-line gas.

The spectroscopic properties of galaxies in the Coma cluster are compared with those of galaxies in rich clusters at $z \sim 0.5$, to investigate the evolution of the star formation history in clusters. Luminous galaxies with $M_V \leq -20$ and post-starburst/post-starforming (k+a) spectra which constitute a significant fraction of galaxies in distant cluster samples are absent in Coma, where spectacular cases of k+a spectra are found instead at $M_V>-18.5$ and represent a significant proportion of the cluster dwarf galaxy population. A simple inspection of their positions on the sky indicates that this type of galaxy does not show a preferential location within the cluster, but the bluest and strongest-lined group of k+a's lies in projection towards the central 1.4 Mpc of Coma and have radial velocities significantly higher than the cluster mean. We find a striking correlation between the positions of these young and strong post-starburst galaxies and substructure in the hot intracluster medium (ICM) identified from {\it XMM-Newton} data, with these galaxies lying close to the edges of two infalling substructures. This result strongly suggests that the interaction with the dense ICM could be responsible for the quenching of the star formation (thus creating the k+a spectrum), and possibly, for any previous starburst. The evolution with redshift of the luminosity distribution of k+a galaxies can be explained by a ``downsizing effect'', with the maximum luminosity/mass of actively star-forming galaxies infalling onto clusters decreasing at lower redshift. We discuss the possible physical origin of this downsizing effect and the implications of our results for current scenarios of environmental effects on the star formation in galaxies.

We present the results of a survey for Ly$\alpha$ emitters at $z\approx 5.7$ based on optical narrow-band ($\lambda_{\rm c} = 8150$ \AA ~ and $\Delta\lambda = 120$ \AA), and broad-band ($B$, $R_{\rm C}$, $I_{\rm C}$, and $z^\prime$) observations of the field surrounding the high redshift quasar, SDSSp J104433.04$-$012522.2, on the 8.2 m Subaru Telescope with the Subaru Prime Focus Camera, Suprime-Cam. This survey covers a sky area of $\approx 720$ arcmin$^2$ and a co-moving volume of $\simeq 2 \times 10^5 ~ h_{0.7}^{-3}$ Mpc$^3$. We have found 20 candidates of Ly$\alpha$ emitters at $z \approx$ 5.7 with $\Delta z \approx 0.1$. Two of them have been confirmed star-forming galaxies at $z=5.655$ and $z=5.687$ from our follow-up optical spectroscopy. We discuss star-formation properties of the 20 objects from a statistical point of view. Our survey leads to a new estimate of the star formation rate density at $z \approx 5.7$, $\sim 1.2 \times 10^{-3} h_{0.7} M_\odot$ yr$^{-1}$ Mpc$^{-3}$.

We present measurements of the cosmic shear correlation in the shapes of galaxies in the Suprime-Cam 2.1 deg^2 R_c-band imaging data. As an estimator of the shear correlation originated from the gravitational lensing, we adopt the aperture mass variance. We detect a non-zero E mode variance on scales between 2 and 40arcmin. We also detect a small but non-zero B mode variance on scales larger than 5arcmin. We compare the measured E mode variance to the model predictions in CDM cosmologies using maximum likelihood analysis. A four-dimensional space is explored, which examines sigma_8, Omega_m, Gamma and zs (a mean redshift of galaxies). We include three possible sources of error: statistical noise, the cosmic variance estimated using numerical experiments, and a residual systematic effect estimated from the B mode variance. We derive joint constraints on two parameters by marginalizing over the two remaining parameters. We obtain an upper limit of Gamma<0.5 for zs>0.9 (68% confidence). For a prior Gamma\in[0.1,0.4] and zs\in[0.6,1.4], we find sigma_8=(0.50_{-0.16}^{+0.35})Omega_m^{-0.37} for flat cosmologies and sigma_8=(0.51_{-0.16}^{+0.29})Omega_m^{-0.34}$ for open cosmologies (95% confidence). If we take the currently popular LCDM model, we obtain a one-dimensional confidence interval on sigma_8 for the 95.4% level, 0.62<\sigma_8<1.32 for zs\in[0.6,1.4]. Information on the redshift distribution of galaxies is key to obtaining a correct cosmological constraint. An independent constraint on Gamma from other observations is useful to tighten the constraint.

Spectroscopic and spectropolarimetric observations of SN 2003dh/GRB 030329 obtained in 2003 May using the Subaru 8.2 m telescope are presented. The properties of the SN are investigated through a comparison with spectra of the Type Ic hypernovae SNe 1997ef and 1998bw. (Hypernovae being a tentatively defined class of SNe with very broad absorption features: these features suggest a large velocity of the ejected material and possibly a large explosion kinetic energy.) Comparison with spectra of other hypernovae shows that the spectrum of SN 2003dh obtained on 2003 May 8 and 9, i.e., 34-35 rest-frame days after the GRB (for z=0.1685), are similar to those of SN 1997ef obtained ~34-42 days after the fiducial time of explosion of that SN. The match with SN 1998bw spectra is not as good (at rest 7300-8000 A, but again spectra obtained ~33-43 days after GRB 980425 are preferred. This indicates that the SN may have intermediate properties between SNe 1997ef and 1998bw. Based on the analogy with the other hypernovae, the time of explosion of SN 2003dh is then constrained to be between -8 and +2 days of the GRB. The Si and O P-Cygni lines of SN 2003dh seem comparable to those of SN 1997ef, which suggests that the ejected mass in SN 2003dh may match that in SN 1997ef. Polarization was marginally detected at optical wavelengths. This is consistent with measurements of the late afterglow, implying that it mostly originated in the interstellar medium of the host galaxy.

Several high-z (z > 5.7) quasars have been found in the course of Sloan Digital Sky Survey. The presence of such very high-z quasars is expected to give constraints on early structure formation. On one hand, it is suggested that these most luminous objects at high redshift are biased toward the highly magnified objects by gravitational lensing. To clarify the effect of gravitational lensing on the high-z quasars, we began the imaging survey of intervening lensing galaxies. Indeed our previous optical image showed that SDSSp J104433.04+012502.2 at z=5.74 is gravitationally magnified by a factor 2. In this paper, we report our new optical imaging of other two high-z quasars, SDSSp J103027.10+052455.0 at z=6.28 and SDSSp J130608.26+035626.3 at z=5.99. Since we find neither intervening galaxy nor counter image with i^{\prime} < 25.4-25.8 around each quasar, we conclude that they are not strongly magnified regardless that a lens galaxy is dusty.

We have carried out a deep imaging survey for H$\alpha$ emitting galaxies at $z\approx$0.24 using a narrowband filter tuned with the redshifted line. The total sky area covered is 706 arcmin$^2$ within a redshift range from 0.234 to 0.252 ($\delta z$=0.018). This corresponds to a volume of 3.9$\times10^3$ Mpc$^3$ when $\Omega_{\rm matter}=0.3$, $\Omega_\Lambda=0.7$, and $H_{\mathrm{0}}$=70 km s$^{-1}$ Mpc$^{-1}$ are adopted. We obtain a sample of 348 H$\alpha$ emitting galaxies whose observed emission-line equivalent widths are greater than 12 \AA. We find an extinction-corrected H$\alpha$ luminosity density of $10^{39.65^{+0.08}_{-0.12}}$ ergs s$^{-1}$ Mpc$^{-3}$. Using the Kennicutt relation between the H$\alpha$ luminosity and star formation rate, the star formation rate density in the covered volume is estimated as $0.036^{+0.006}_{-0.012}$ $M_{\sun}$ yr$^{-1}$ Mpc$^{-3}$. This value is higher by a factor of 3 than the local SFR density.

We report the discovery of a large-scale structure of Lyman alpha emitters (LAEs) at z=4.86 based on wide-field imaging with the prime-focus camera (Suprime-Cam) on the Subaru telescope. We observed a 25' x 45' area of the Subaru Deep Field in a narrow band (NB, lambdaC=7126 A and FWHM=73 A) together with R and i'. We isolate from these data 43 LAE candidates down to NB=25.5 mag using color criteria. Follow-up spectroscopy of five candidates suggests the contamination by low-z objects to be ~ 20%. We find that the LAE candidates are clustered in an elongated region on the sky of 20 Mpc in width and 50 Mpc in length at z=4.86, which is comparable in size to present-day large-scale structures (we adopt H0=70 km/s/Mpc, Omega0=0.3, lambda0=0.7). This elongated region includes a circular region of 12 Mpc radius of higher surface overdensity (delta=2), which may be the progenitor of a cluster of galaxies. Assuming this circular region to be a sphere with a spatial overdensity of 2, we compare our observation with predictions by Cold Dark Matter models. We find that an Omega0=0.3 flat model with sigma8=0.9 predicts the number of such spheres consistent with the observed number (one sphere in our survey volume) if the bias parameter of LAEs is b \simeq 6. This value suggests that the typical mass of dark haloes hosting LAEs at z \simeq 5 is of the order of 10^12 Msolar. Such a large mass poses an interesting question about the nature of LAEs.

We report on the discovery of a very narrow-line star forming object beyond redshift of 5. Using the prime-focus camera, Suprime-Cam, on the 8.2 m Subaru telescope together with a narrow-passband filter centered at $\lambda_{\rm c}$ = 8150 \AA with passband of $\Delta\lambda$ = 120 \AA, we have obtained a very deep image of the field surrounding the quasar SDSSp J104433.04$-$012502.2 at a redshift of 5.74. Comparing this image with optical broad-band images, we have found an object with a very strong emission line. Our follow-up optical spectroscopy has revealed that this source is at a redshift of $z=5.655\pm0.002$, forming stars at a rate $\sim 13 ~ h_{0.7}^{-2} ~ M_\odot$ yr$^{-1}$. Remarkably, the velocity dispersion of Ly$\alpha$-emitting gas is only 22 km s$^{-1}$. Since a blue half of the Ly$\alpha$ emission could be absorbed by neutral hydrogen gas, perhaps in the system, a modest estimate of the velocity dispersion may be $\gtrsim$ 44 km s$^{-1}$. Together with a linear size of 7.7 $h_{0.7}^{-1}$ kpc, we estimate a lower limit of the dynamical mass of this object to be $\sim 2 \times 10^9 M_\odot$. It is thus suggested that LAE J1044$-$0123 is a star-forming dwarf galaxy (i.e., a subgalactic object or a building block) beyond redshift 5 although we cannot exclude a possibility that most Ly$\alpha$ emission is absorbed by the red damping wing of neutral intergalactic matter.

A large spectroscopic survey is constructed of galaxies in the Coma cluster. The survey covers a wide area (1 deg$^2$) to deep magnitudes (R ~ 19.5), covering both the core (high density) and outskirts (intermediate to low density) of the cluster. The spectroscopic sample consists of a total of 1191 galaxies, of which, 760 galaxies are confirmed members of the Coma cluster. A statistical technique is developed to correct the spectroscopic sample for incompletness. The corrected sample is then used to construct R-band luminosity function (LF) spanning a range of 7 magnitudes (-23 < M_R < -16) both at the core and outskirts of the cluster. Dependence of the LF on local environment in Coma is explored. The LFs are found to be the same, within the errors, between the inner and outer regions and close to those from recent measurements for field galaxies. The total B-band LF for the Coma cluster, fitted to a Schecter form is also measured and shows a dip at M(B) = -18 mag., in agreement with previous studies. The implications of this feature are discussed. The LF is studied in B-R color intervals and shows a steep faint-end slope for red (B-R > 1.35) galaxies, both at the core and outskirts of the cluster. This population of low luminosity red galaxies has a higher surface density than the blue (B-R < 1.35) star-forming population, dominating the faint-end of the Coma cluster LF. It is found that relative number of high surface brightness galaxies is larger at the cluster core, implying destruction of low surface brightness galaxies in dense core environment.

During the course of our optical deep survey program on L$\alpha$ emitters at $z \approx 5.7$ in the sky area surrounding the quasar SDSSp J104433.04$-$012502.2 at $z=5.74$, we found that a faint galaxy with $m_B$(AB) $\approx 25$ is located at \timeform{1".9} southwest of the quasar. Its broad-band color properties from $B$ to $z^\prime$ suggest that the galaxy is located at a redshift of $z \sim 1.5$ -- 2.5. This is consistent with no strong emission line in our optical spectroscopy. Since the counter image of the quasar cannot be seen in our deep optical images, the magnification factor seems not to be very high. Our modest estimate is that this quasar is gravitationally magnified by a factor of 2.

We have built an 80-mega pixels (10240 X 8192) mosaic CCD camera, called Suprime-Cam, for the wide-field prime focus of the 8.2 m Subaru telescope. Suprime-Cam covers a field of view 34' X 27', a unique facility among the the 8-10 m class telescopes, with a resolution of 0.''202 per pixel. The focal plane consists of ten high-resistivity 2k X 4k CCDs developed by MIT Lincoln Laboratory, which are cooled by a large stirling-cycle cooler. The CCD readout electronics was designed to be scalable, which allows the multiple read-out of tens of CCDs. It takes 50 seconds to readout entire arrays. We designed a filter-exchange mechanism of the jukebox type that can hold up to ten large filters (205 X 170 X 15 mm^3). The wide-field corrector is basically a three-lens Wynne-type, but has a new type of atmospheric dispersion corrector. The corrector provides a flat focal plane and an un-vignetted field of view of 30' in diameter. The achieved co-planarity of the focal array mosaic is smaller than 30 micron peak-to-peak, which realizes mostly the seeing limited image over the entire field. The median seeing in the Ic-band, measured over one year and a half, is 0.''61. The PSF anisotropy in Suprime-Cam images, estimated by stellar ellipticities, is about 2% under this median seeing condition. At the time of commissioning, Suprime-Cam had the largest survey speed, which is defined as the field of view multiplied by the primary mirror area of the telescope, among those cameras built for sub-arcsecond imaging.

We have identified intracluster planetary nebulae in a Virgo cluster core field by imaging with the Subaru Suprime-Cam through two narrow band filters centered at the redshifted wavelengths of the [OIII] $\lambda=5007$ \AA and the H$\alpha$ $\lambda = 6563$ \AA lines; broad-band images in V and R were acquired to check for emission in the adjacent continuum. Emission line objects in Virgo are then selected from the two-color diagram ([OIII] -- Halpha) vs. ([OIII] -- (V+R)), which is calibrated using PNe in M84 (Jacoby et al. 1990). Using both [OIII] and Halpha allows us to distinguish bona-fide planetary nebulae from high redshift emission-line galaxies at the bright end of the [OIII] luminosity function. Spectroscopic observations of a subsample of these objects were made at the TNG and at the VLT, in a region around M84 and in an intracluster field respectively. The observations confirm the efficiency of the combined [OIII]+H$\alpha$ imaging to identify true PNe. We also obtained the first spectrum of an intracluster PN which shows the [OIII] doublet with S/N > 10 and its Halpha emission. From the results based on the spectroscopic follow-up, we derive a lower limit to the fraction of the Virgo cluster light contributed by the intracluster stars at the surveyed position in the cluster core: it amounts to 10%.

We have identified 38 candidates of intracluster planetary nebulae (ICPNe) in a $34'\times27'$ field in the core of the Virgo cluster based on the Suprime-Cam imaging through two narrow-band filters centered at the redshifted wavelengths of the [OIII]$\lambda=5007\AA$ and the Halpha \lambda=6563\AA lines. Broad-band images in V and R bands are used to check for any emissions in the adjacent continuum. We describe the method briefly and present the list of intracluster planetary nebulae candidates, together with their finding charts. The ICPN candidates show a highly inhomogeneous distribution, which may suggest an association with the M86-M84 subcluster. Fraction of diffuse intracluster light with respect to total light in galaxies is estimated to be about 10%, leading to an estimate of about 20% for the baryon fraction. Spectroscopic follow up and a wider survey are critical to reveal the nature of intracluster stellar population.

We study the evolution of elliptical galaxies at z >~ 1 on the basis of a sample of 247 Extremely Red Objects (EROs) with R-Ks >= 3.35 (AB) and Ks <= 22.1 (AB) constructed from BVRi'z'JHKs multicolor data of a 114 arcmin2 area in the Subaru/XMM Deep Survey Field. By fitting template spectra of old ellipticals (OEs) and young, dusty starbursts (DSs) to the multicolor data, we classify EROs into these two classes and estimate their redshifts. We find that 58% of the EROs in our sample belong to the OE class and that 24% of the OEs are fit by a spectrum having a disk component with the B-band bulge-to-total luminosity ratio of <= 0.9. We compare the observed surface density of OEs with predictions from passive evolution models. We then derive rest-frame B-band luminosity functions (LFs) of OEs in our sample at z=1-1.5 and 1.5-2.5. We find that the LF at z=1-1.5 roughly agrees with the LF of local ellipticals if a dimming of 1.3 mag from z=1.25 to the present epoch is assumed. On the other hand, the amplitude of the LF at z=1.5-2.5 is found to be lower than that of the local LF by a factor of ~ 3 over the whole range of magnitude observed. Taking account of a strong decrease in the number density of morphologically classified early-type galaxies at z >~ 1.5 found by several authors, we conclude that the majority of ellipticals seen at present have not established either a red color or a smooth 1/4-law profile before z ~ 1.5$. (abridged)

We report the first result of weak gravitational lensing survey on a 2.1 sq deg Rc-band image taken with a wide field camera (Suprime-Cam) on the prime focus of 8.2 m Subaru Telescope. The weak lensing mass reconstruction is applied to the data to search for dark matter halos of cluster scale; M >= 10^14 solar mass. The reconstructed convergence field is divided by 1-sigma noise to obtain the signal-to-noise ratio map (S/N-map) of the detection. Local maxima and minima are searched on the S/N-map and the probability distribution function (PDF) of the peaks are created to compare with model predictions. We found excess over noise PDF created from the randomized realization on both positive and negative sides. Negative peaks imply the presence of voids in the dark matter distribution and this is the first report of the detection. Positive peaks, on the other hand, represent the dark matter halos and the number count of the halos on the 2.1 sq deg image is 4.9 +- 2.3 for S/N > 5 where the Gaussian smoothing radius of the convergence map is 1'. The result is consistent with the prediction assuming the Press-Schechter mass function and the NFW halo profile under the cluster normalized CDM cosmology. Although the present statistics is not enough due to the limited field of view, this work demonstrates that dark matter halo count via weak lensing is a promising way to test the paradigm of structure formation and cosmological model.

We present the results of a survey for emission-line objects based on optical intermediate-band ($\lambda_{\rm c}$ = 5736 \AA ~ and $\Delta\lambda$ = 280 \AA) and broad-band ($B$, $V$, $R$, and $i^\prime$) observations of the Subaru/XMM-Newton Deep Field on the 8.2 m Subaru telescope with the Subaru Prime Focus Camera, Suprime-Cam. All the data were obtained during the guaranteed time observations of the Suprime-Cam instrument. The intermediate-band image covered a sky area with 10\minpoint62 $\times 12\minpoint40 \approx 132$ arcmin$^2$ in the Subaru/XMM-Newton Deep Field (Ouchi et al.). Using this image, we have found 23 emission-line sources whose observed emission-line equivalent widths are greater than 250 \AA. Their optical multicolor properties indicate that six emission-line sources are Ly$\alpha$ emitters at $z \approx$ 3.7 ($\Delta z \approx 0.22$). They are either intense starburst galaxies or active galactic nuclei like quasars at $z \approx$ 3.7. Two more emission-line sources may also be Ly$\alpha$ emitters at $z \approx$ 3.7 although their multicolor properties are marginal. Among the remaining 15 emission-line objects, eight objects appear strong emission-line galaxies at lower redshift; e.g., [O {\sc ii}] $\lambda$3727 emitters at $z \approx 0.54$, H$\beta$ at $z \approx 0.18$, or [O {\sc iii}]$\lambda$5007 emitters at $z \approx 0.15$. The remaining seven objects are unclassified because they are too faint to be detected in broad-band images. We discuss observational properties of these strong emission-line sources. In particular, our data allow us to estimate the star formation density at $z \approx 3.7$ for the first time.

We report on early results from deep and wide-field narrow-band imaging on a 543 arcmin^2 area of the Subaru Deep Field. We find 87 Lya emitters (LAEs) at z=4.86(+/-0.03) which are photometrically selected by a combination of two broad bands (R and i') and one narrow band (NB711; 7126A FWHM=73A). We derive the luminosity functions (LFs) of the LAEs at Lya luminosity and at UV-continuum (1700A) luminosity. The LFs show little evolution between z=3.4 and z=4.86 either in Lya or UV-continuum emission. The amplitude of the LAE LF tends to decline at the bright magnitudes more rapidly than that of the LBG LF at similar redshifts. We calculate the angular correlation function of our LAEs up to ~15 arcmin separations. The angular correlation w(theta) is found to increase with decreasing angular separations, showing a clear signal of clustering. It is also found that the distribution of LAEs shows a large density gradient with a scale of >~15 arcmin, which would indicate the existence of a large-scale structure of LAEs on >~20h^(-1) Mpc scales.We fit the observed correlation function by Aw theta^(-0.8) to find Aw=29arcsec^(0.8). The estimated correlation length is r0=3.5(+/-0.3)h^(-1)Mpc in comoving units (Omega_m=0.3 and Omega_Lambda=0.7), which is slightly larger than the value for z~4 LBGs with i'<26. We calculate the angular correlation function for two sub-samples of the 87 LAEs divided by Lya luminosity, UV-continuum luminosity, and Lya equivalent width (EW). The Lya bright sub-sample shows a larger correlation amplitude than the Lya faint sub-sample, while no significant difference is found for the sub-samples divided by UV-continuum luminosity or EW.

During the course of our deep optical imaging survey for Ly alpha emitters at z approximately 5.7 in the field around the z=5.74 quasar SDSSp J104433.04-012502.2, we have found a candidate strong emission-line source. Follow-up optical spectroscopy shows that the emission line profile of this object is asymmetric, showing excess red-wing emission. These properties are consistent with an identification of Ly alpha emission at a redshift of z=5.687 +/- 0.002. The observed broad line width, Delta V_{FWHM} ~= 340 km s^{-1} and excess red-wing emission also suggest that this object hosts a galactic superwind.

We found a very large, 35 kpc, emission-line region around the Seyfert type 2 galaxy NGC 4388, using deep narrow-band imaging with the prime focus camera (Suprime-Cam) of the Subaru telescope. This region consists of many faint gas clouds or filaments, and extends northeastwards from the galaxy. The typical H-alpha luminosity L(H-alpha) of the clouds is 10^37 erg s^-1, and the total L(H-alpha) of the region within 10 kpc from the nucleus is 2 x 10^38 erg s^-1, which corresponds to an ionized gas mass of 10^5 M_solar. The map of the emission-line intensity ratio I([OIII])/I(H-alpha) indicates that the inner (r < 12 kpc) region of the VEELR may be excited by nuclear ionizing radiation. The excitation mechanism of the outer (r > 12 kpc) region is unclear, but it is likely that the nuclear radiation is also a dominant source of its ionization. We discuss the origin of the ionized gas. Two plausible origins of the gas in the VEELR are (i) the tidal debris resulting from a past interaction with a gas-rich dwarf galaxy, i.e., a minor merger, or (ii) the interstellar medium of NGC 4388, stripped by the ram pressure of the hot intracluster medium of the Virgo cluster.

We investigate the radial dependence of the spectroscopic properties, in particular the Mg2, <Fe> and H beta spectroscopic indices, in a sample of galaxies spanning a wide range of absolute luminosity in the Coma cluster. After allowing for the magnitude dependence of these indices, we find a significant gradient in Mg2, in the sense that galaxies in the core of the cluster have stronger Mg2. We find only weak gradients in <Fe> and H beta. Using the model grids presented in an earlier paper in this series, we attribute the Mg2 gradient to changes in metal abundance. One possible mechanism to create this abundance gradient is pressure confinement by the intracluster medium of material from Supernova driven winds early in the history of the galaxies.

We discuss the properties of galaxies around the radio galaxy 3C 324 at z=1.2 based on BVRIK' multi-band imaging data. We have applied a photometric-redshift technique to objects in the 3C 324 field, and identified 35 objects as plausible cluster members. We have found that red and luminous members are concentrated in a small region enclosed by a circle of 40" radius (0.33 Mpc at z=1.2 for \Omega_0=0.3, \lambda_0=0.7, H_0=70km/s/Mpc cosmology) from the 3C 324 galaxy. The 3C 324 cluster is probably much more compact in size compared with the local clusters. We constructed a K'-band luminosity function of the cluster members and fit a Schechter function, and found the characteristic magnitude to be K*_{AB}=20.2+-0.6. This value is consistent with the extrapolation of the pure passive evolution seen for z<1 clusters. We have identified eight bright galaxies which form a red color-magnitude sequence. The slope of the sequence is consistent with the passive evolution model down to K'_{AB}<22; we also found that there is no clear age variation in these bright red galaxies. However, seven out of these eight galaxies exhibit a significant excess in the rest UV light with respect to the passive evolution model. This may suggest that the massive early-type galaxies in this high-redshift cluster are still forming stars to some extent. We have confirmed a truncation of the color-magnitude sequence at K'_{AB}\sim22; faint passively-evolving galaxies may not yet be present in this cluster at z\sim1.2. The overall color distribution of the cluster members, selected by the photometric redshift technique, is found to be very broad. We derived the fraction of blue galaxies in this cluster following a

(abr.) We present the analysis of the spectroscopic catalog of galaxies in the Coma cluster from Mobasher et al. This catalog comprises ~300 spectra of cluster members with magnitudes M_B=-20.5 to -14 in two areas of ~1 X 1.5 Mpc towards the center and the SW region of the cluster. In this paper the analysis is restricted to the 257 galaxies with no emission lines in their spectra. The strength of the age-sensitive indices is found to correlate with galaxy magnitude over the whole magnitude range. Similarly, the metallicity-sensitive indices anticorrelate with magnitude. By comparing the observed indices with model grids based on the Padova isochrones, we derive luminosity-weighted ages and metallicities. We present the distributions of ages and metallicities for galaxies in various magnitude bins. The mean metallicity decreases with galaxy magnitude and, at a given luminosity, appears to be generally lower for galaxies in the SW region of Coma as compared to the center of the cluster. A broad range of ages, from younger than 3 Gyr to older than 9 Gyr, is found in galaxies of any magnitude. However, systematic trends of age with luminosity are present among galaxies in the central field. In the central Mpc of Coma, a large fraction of galaxies at any luminosity (50-60% of the giants, >30 % of the dwarfs) show no evidence in their central regions of star formation occurred at redshift z<2, while the proportion of galaxies with significant star formation occurring at intermediate (0.35<z<2) and low (z<0.35) redshifts is found to depend on galaxy luminosity. Faint galaxies with young luminosity-weighted ages appear to have a bimodal metallicity distribution that, if confirmed, would point to a composite formation scenario. Finally, we present an interpretation of the index-magnitude relations observed.

We study the luminosity function and the correlation function of about 1200 z~4 Lyman break galaxies (LBGs) with i'<26 that are photometrically selected from deep BRi' imaging data of a 618 arcmin^2 area in the Subaru/XMM-Newton Deep Field taken with Subaru Prime Focus Camera. The contamination and completeness of our LBG sample are evaluated, on the basis of the Hubble Deep Field-North (HDF-N) objects, to be 17% and 45%, respectively. We derive the UV (rest 1700A) luminosity functions (LFs) and find a large population of UV-luminous galaxies at z~4. The LFs of the red and blue subsamples imply that the bright LBGs are redder in the UV continuum than the average color of the LBGs. Then we calculate the correlation function over theta = 2''-1000'' and find that it is fitted fairly well by a power law, omega(theta)=A_omega theta^(-0.8), with A_omega=0.71 +/- 0.26. We estimate the correlation length r_0 (in comoving units) of the two-point spatial correlation function xi(r) = (r/r_0)^(-1.8) to be r_0=2.7 +0.5/-0.6 h^(-1) Mpc (Omega_m=0.3 and Omega_Lambda=0.7). The correlation function shows an excess of omega (theta) on small scales (theta < 5''), departing from the power-law fit at > 3 sigma significance level. Interpreting this as being due to galaxy mergers, we evaluate the fraction of galaxies undergoing mergers to be 3.0 +/- 0.9%, which is significantly smaller than those of galaxies at intermediate redshifts.

The ages of stellar populations in 52 elliptical and S0 galaxies in the Coma cluster are investigated using a new spectroscopic sample. More than 40% of the S0s are found to have experienced recent star formation in their central regions during the last ~5 Gyrs, while such activity is absent in the ellipticals. Galaxies in this sample have absolute magnitudes in the range -20.5 < M_B < -17.5, and the fraction of S0 galaxies with recent star formation is higher at fainter luminosities. The observed luminosity range of S0 galaxies with signs of recent star formation activity is consistent with them being the descendants of typical star-forming spirals at intermediate redshift whose star formation has been halted as a consequence of the dense environment.

We study the clustering properties of about 1200 z~4 Lyman Break Galaxy (LBG) candidates with i'<26 which are selected by color from deep BRi' imaging data of a 618 arcmin^2 area in the Subaru/XMM-Newton Deep Field taken with Subaru Prime Focus Camera. The contamination and completeness of our LBG sample are evaluated, on the basis of the Hubble Deep Field North (HDFN) objects, to be 17% and 45%, respectively. We derive the angular correlation function over theta = 2''-1000'', and find that it is fitted fairly well by a power law, omega(theta)=A_omega theta^{-0.8}, with A_omega = 0.71 +/- 0.26. We then calculate the correlation length r0 (in comoving units) of the two-point spatial correlation function xi(r) = (r/r0)^{-1.8} from A_omega using the redshift distribution of LBGs derived from the HDFN, and find r0=2.7 (+0.5/-0.6) h^{-1} Mpc in a Lambda-dominated universe (Omega_m=0.3 and Omega_Lambda=0.7). This is twice larger than the correlation length of the dark matter at z~4 predicted from an analytic model by Peacock & Dodds but about twice smaller than that of bright galaxies predicted by a semi-analytic model of Baugh et al. We find an excess of omega(theta) on small scales (theta < 5'') departing from the power law fit over 3 sigma significance levels. Interpreting this as due to galaxy mergers, we estimate the fraction of galaxies undergoing mergers in our LBG sample to be 3.0 +/- 0.9%, which is significantly smaller than those of galaxies at intermediate redshifts.

This is the second paper in a series studying the photometric and spectroscopic properties of galaxies of different luminosities in the Coma cluster. The sample selection, spectroscopic observations and completeness functions are presented here. To study the spectral properties of galaxies as a function of their local environment, two fields were selected for spectroscopic observations to cover both the core (Coma1) and outskirts (ie. south-west of the core and centered on NGC4839)- (Coma3) of the cluster. Medium resolution spectroscopy (6-9 \AA) was carried out for a total of 490 galaxies in both fields (302 in Coma1 and 188 in Coma3), using the WYFFOS multi-fiber spectrograph on the William Herschel Telescope. The galaxies cover a range of $12 < R < 20$, corresponding to -23 < M_R < -15 (H0=65 km/sec/Mpc). The redshifts are measured with an accuracy of 100 km/sec. The spectral line strengths and equivalent widths are also measured for the same galaxies and analysed in Poggianti et al (2001- paper III). A total of 189 (Coma1) and 90 (Coma3) galaxies are identified as members of the Coma cluster. An analysis of the colors show that only two members of the Coma cluster in our sample have B-R > 2. The completeness functions for the spectroscopic sample is presented.

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

We investigate pure luminosity evolution models for early-type (elliptical and S0) galaxies (i.e., no number density change or morphology transition), and examine whether these models are consistent with observed number counts in the B, I and K bands and redshift distributions of two samples of faint galaxies selected in the I and K bands. The models are characterized by the star formation time scale $\tau_{SF}$ and the time $t_{gw}$ when galactic wind blows in addition to several other conventional parameters. We find the single-burst model ($\tau_{SF}$=0.1 Gyr and $t_{gw}$=0.353 Gyr), which is known to reproduce the photometric properties of early-type galaxies in clusters, is inconsistent with redshift distributions of early-type galaxies in the field environment due to overpredictions of galaxies at $z\gsim1.4$ even with strong extinction which is at work until $t_{gw}$. In order for dust extinction to be more effective, we change $\tau_{SF}$ and $t_{gw}$ as free parameters, and find that models with $\tau_{SF}\gsim0.5$ Gyr and $t_{gw}>1.0$ Gyr can be made consistent with both the observed redshift distributions and number counts, if we introduce strong extinction $(E(B-V)\geq1$ as a peak value). These results suggest that early-type galaxies in the field environment do not have the same evolutionary history as described by the single-burst model.

We carry out CCD photometry of galaxies in the 5.25 square region centered on Coma cluster down to $M_R=-16.0$, beyond the limit of conventional morphological classification. We use the angular two-point correlation function as well as radial profiles in order to characterize the luminosity segregation. We find strong luminosity segregation for our total sample over the magnitude range of $-20 \leq M_R \leq -16$, which is not entirely accounted for in terms of the morphology-density relation that is known to exist only for bright galaxies. We use a single consistent parameter, the degree of luminosity concentration, to parameterize the morphology of galaxies over the wide magnitude range, where both giant and dwarf galaxies are included. Galaxies with high central concentration (HCC) show strong luminosity segregation, i.e. their clustering strength depends strongly on luminosity while those with low central concentration (LCC) show almost no luminosity segregation. Radial density profile shows that brighter HCC-type galaxies tend to more strongly concentrate near the cluster center while LCC-type galaxies do not show such a dependence on luminosity. We show that these results are tenable against the contamination by field galaxies and uncertainties in our method of classification.