We present the analysis of a very high-magnification ($A\sim 900$) microlensing event KMT-2019-BLG-1953. A single-lens single-source (1L1S) model appears to approximately delineate the observed light curve, but the residuals from the model exhibit small but obvious deviations in the peak region. A binary lens (2L1S) model with a mass ratio $q\sim 2\times 10^{-3}$ improves the fits by $\Delta\chi^2=181.8$, indicating that the lens possesses a planetary companion. From additional modeling by introducing an extra planetary lens component (3L1S model) and an extra source companion (2L2S model), it is found that the residuals from the 2L1S model further diminish, but claiming these interpretations is difficult due to the weak signals with $\Delta\chi^2=16.0$ and $13.5$ for the 3L1S and 2L2L models, respectively. From a Bayesian analysis, we estimate that the host of the planets has a mass of $M_{\rm host}=0.31^{+0.37}_{-0.17}~M_\odot$ and that the planetary system is located at a distance of $D_{\rm L}=7.04^{+1.10}_{-1.33}~{\rm kpc}$ toward the Galactic center. The mass of the securely detected planet is $M_{\rm p}=0.64^{+0.76}_{-0.35}~M_{\rm J}$. The signal of the potential second planet could have been confirmed if the peak of the light curve had been more densely observed by followup observations, and thus the event illustrates the need for intensive followup observations for very high-magnification events even in the current generation of high-cadence surveys.

We report the discovery of a planet in the microlensing event OGLE-2018-BLG-1269, with planet-host mass ratio $q \sim 6\times10^{-4}$, i.e., $0.6$ times smaller than the Jupiter/Sun mass ratio. Combined with the $Gaia$ parallax and proper motion, a strong one-dimensional constraint on the microlens parallax vector allows us to significantly reduce the uncertainties of lens physical parameters. A Bayesian analysis that ignores any information about light from the host yields that the planet is a cold giant $(M_{2} = 0.69_{-0.22}^{+0.44}\,M_{\rm J})$ orbiting a Sun-like star $(M_{1} = 1.13_{-0.35}^{+0.72}\,M_{\odot})$ at a distance of $D_{\rm L} = 2.56_{-0.62}^{+0.92}\,{\rm kpc}$. The projected planet-host separation is $a_{\perp} = 4.61_{-1.17}^{+1.70}\,{\rm au}$. Using {\it Gaia} astrometry, we show that the blended light lies $\lesssim 12\,$mas from the host and therefore must be either the host star or a stellar companion to the host. An isochrone analysis favors the former possibility at $>99.6\%$. The host is therefore a subgiant. For host metallicities in the range of $0.0 \leq {\rm [Fe/H]} \leq +0.3$, the host and planet masses are then in the range of $1.16 \leq M_{1}/M_{\odot} \leq 1.38$ and $0.74 \leq M_{2}/M_{\rm J} \leq 0.89$, respectively. Low host metallicities are excluded. The brightness and proximity of the lens make the event a strong candidate for spectroscopic followup both to test the microlensing solution and to further characterize the system.

We report the discovery and analysis of the planetary microlensing event OGLE-2017-BLG-0406, which was observed both from the ground and by the ${\it Spitzer}$ satellite in a solar orbit. At high magnification, the anomaly in the light curve was densely observed by ground-based-survey and follow-up groups, and it was found to be explained by a planetary lens with a planet/host mass ratio of $q=7.0 \times 10^{-4}$ from the light-curve modeling. The ground-only and ${\it Spitzer}$-"only" data each provide very strong one-dimensional (1-D) constraints on the 2-D microlens parallax vector $\bf{\pi_{\rm E}}$. When combined, these yield a precise measurement of $\bf{\pi_{\rm E}}$, and so of the masses of the host $M_{\rm host}=0.56\pm0.07\,M_\odot$ and planet $M_{\rm planet} = 0.41 \pm 0.05\,M_{\rm Jup}$. The system lies at a distance $D_{\rm L}=5.2 \pm 0.5 \ {\rm kpc}$ from the Sun toward the Galactic bulge, and the host is more likely to be a disk population star according to the kinematics of the lens. The projected separation of the planet from the host is $a_{\perp} = 3.5 \pm 0.3 \ {\rm au}$, i.e., just over twice the snow line. The Galactic-disk kinematics are established in part from a precise measurement of the source proper motion based on OGLE-IV data. By contrast, the ${\it Gaia}$ proper-motion measurement of the source suffers from a catastrophic $10\,\sigma$ error.

We present the analysis of microlensing event OGLE-2006-BLG-284, which has a lens system that consists of two stars and a gas giant planet with a mass ratio of $q_p = (1.26\pm 0.19) \times 10^{-3}$ to the primary. The mass ratio of the two stars is $q_s = 0.289\pm 0.011$, and their projected separation is $s_s = 2.1\pm 0.7\,$AU, while the projected separation of the planet from the primary is $s_p = 2.2\pm 0.8\,$AU. For this lens system to have stable orbits, the three-dimensional separation of either the primary and secondary stars or the planet and primary star must be much larger than that these projected separations. Since we do not know which is the case, the system could include either a circumbinary or a circumstellar planet. Because there is no measurement of the microlensing parallax effect or lens system brightness, we can only make a rough Bayesian estimate of the lens system masses and brightness. We find host star and planet masses of $M_{L1} = 0.35^{+0.30}_{-0.20}\,M_\odot$, $M_{L2} = 0.10^{+0.09}_{-0.06}\,M_\odot$, and $m_p = 144^{+126}_{-82}\,M_\oplus$, and the $K$-band magnitude of the combined brightness of the host stars is $K_L = 19.7^{+0.7}_{-1.0}$. The separation between the lens and source system will be $\sim 90\,$mas in mid-2020, so it should be possible to detect the host system with follow-up adaptive optics or Hubble Space Telescope observations.

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.

We present a weak-lensing analysis of X-ray galaxy groups and clusters selected from the XMM-XXL survey using the first-year data from the Hyper Suprime-Cam (HSC) Subaru Strategic Program. Our joint weak-lensing and X-ray analysis focuses on 136 spectroscopically confirmed X-ray-selected systems at 0.031 < z < 1.033 detected in the 25sqdeg XXL-N region. We characterize the mass distributions of individual clusters and establish the concentration-mass (c-M) relation for the XXL sample, by accounting for selection bias and statistical effects, and marginalizing over the remaining mass calibration uncertainty. We find the mass-trend parameter of the c-M relation to be \beta = -0.07 \pm 0.28 and the normalization to be c200 = 4.8 \pm 1.0 (stat) \pm 0.8 (syst) at M200=10^{14}Msun/h and z = 0.3. We find no statistical evidence for redshift evolution. Our weak-lensing results are in excellent agreement with dark-matter-only c-M relations calibrated for recent LCDM cosmologies. The level of intrinsic scatter in c200 is constrained as \sigma(\ln[c200]) < 24% (99.7% CL), which is smaller than predicted for the full population of LCDM halos. This is likely caused in part by the X-ray selection bias in terms of the relaxation state. We determine the temperature-mass (Tx-M500) relation for a subset of 105 XXL clusters that have both measured HSC lensing masses and X-ray temperatures. The resulting Tx-M500 relation is consistent with the self-similar prediction. Our Tx-M500 relation agrees with the XXL DR1 results at group scales, but has a slightly steeper mass trend, implying a smaller mass scale in the cluster regime. The overall offset in the Tx-M500 relation is at the $1.5\sigma$ level, corresponding to a mean mass offset of (34\pm 20)%. We also provide bias-corrected, weak-lensing-calibrated M200 and M500 mass estimates of individual XXL clusters based on their measured X-ray temperatures.

We first present new Subaru narrow-band observations of the Ly$\alpha$ halo around the quasi-stellar object (QSO) CFHQ J232908$-$030158 at $z=6.42$, which appears the most luminous and extended halo at $z>5$ ($L_{Ly\alpha}=9.8\times10^{43}$ erg s$^{-1}$ within $37$ pkpc diameter). Then, combining these measurements with available data in the literature, we find two different evolutions of QSOs' Ly$\alpha$ halos. First is a possible short-term evolution with QSO age seen in four $z>6$ QSOs. We find the anti-correlation between the Ly$\alpha$ halo scales with QSOs' IR luminosity, with J2329-0301's halo being the brightest and largest. It indicates that ionizing photons escape more easily out to circum-galactic regions when host galaxies are less dusty. We also find a positive correlation between IR luminosity and black hole mass ($M_\text{BH}$). Given $M_\text{BH}$ as an indicator of QSO age, we propose a hypothesis that a large Ly$\alpha$ halo mainly exists around QSOs in the young phase of their activity due to a small amount of dust. The second is an evolution with cosmic time seen over $z\sim2-5$. We find the increase of surface brightness toward lower-redshift with a similar growth rate to that of dark matter halos (DHs) which evolve to $M_\text{DH}=10^{12}-10^{13}$ M$_\odot$ at $z=2$. The extent of Ly$\alpha$ halos is also found to increase at a rate scaling with the virial radius of growing DHs, $r_\text{vir} \propto M_\text{DH}^{1/3}$($1+z$)$^{-1}$. These increases are consistent with a scenario that the CGM around QSOs evolves in mass and size keeping pace with hosting DHs.

We estimate the bias on weak lensing mass measurements of shear-selected galaxy cluster samples. The mass bias is expected to be significant because constructions of cluster samples from peaks in weak lensing mass maps and measurements of cluster masses from their tangential shear profiles share the same noise. We quantify this mass bias from large sets of mock cluster samples with analytical density profiles and realistic large-scale structure noise from ray-tracing simulations. We find that, even for peaks with signal-to-noise ratio larger than $4.0$ in weak lensing mass maps constructed in a deep survey with a high source galaxy number density of $30~\mathrm{arcmin}^{-2}$, derived weak lensing masses for these shear-selected clusters are still biased high by $\sim55\%$ on average. Such a large bias mainly originates from up-scattered low mass objects, which is an inevitable consequence of selecting clusters with a noisy observable directly linked to the mass measurement. We also investigate the dependence of the mass bias on different physical and observational parameters, finding that the mass bias strongly correlates with cluster redshifts, true halo masses, and selection signal-to-noise thresholds, but having moderate dependence on observed weak lensing masses and survey depths. This bias, albeit considerable, can still be modeled accurately in statistical studies of shear-selected clusters as the intrinsic scatter around the mean bias is found to be reasonable in size. We demonstrate that such a bias can explain the deviation in X-ray properties previously found on a shear-selected cluster sample. Our result will be useful for turning large samples of shear-selected clusters available in future surveys into potential probes of cosmology and cluster astrophysics.

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.

We present the analysis of the binary-lens microlensing event OGLE-2013-BLG-0911. The best-fit solutions indicate the binary mass ratio of q~0.03 which differs from that reported in Shvartzvald+2016. The event suffers from the well-known close/wide degeneracy, resulting in two groups of solutions for the projected separation normalized by the Einstein radius of s~0.15 or s~7. The finite source and the parallax observations allow us to measure the lens physical parameters. The lens system is an M-dwarf orbited by a massive Jupiter companion at very close (M_{host}=0.30^{+0.08}_{-0.06} M_{Sun}, M_{comp}=10.1^{+2.9}_{-2.2} M_{Jup}, a_{exp}=0.40^{+0.05}_{-0.04} au) or wide (M_{host}=0.28^{+0.10}_{-0.08} M_{Sun}, M_{comp}=9.9^{+3.8}_{-3.5}M_{Jup}, a_{exp}=18.0^{+3.2}_{-3.2} au) separation. Although the mass ratio is slightly above the planet-brown dwarf (BD) mass-ratio boundary of q=0.03 which is generally used, the median physical mass of the companion is slightly below the planet-BD mass boundary of 13M_{Jup}. It is likely that the formation mechanisms for BDs and planets are different and the objects near the boundaries could have been formed by either mechanism. It is important to probe the distribution of such companions with masses of ~13M_{Jup} in order to statistically constrain the formation theories for both BDs and massive planets. In particular, the microlensing method is able to probe the distribution around low-mass M-dwarfs and even BDs which is challenging for other exoplanet detection methods.

With $q=(4.09\pm 0.27) \times 10^{-5}$, KMT-2019-BLG-0842Lb is the second cold planet with planet-host mass ratio $q$ below that of Uranus relative to the Sun ($q=4.37\times 10^{-5}$), the first having been discovered only in 2018. Hence, such cold, sub-Uranus-ratio planets are probably not rare. We discuss possible improvements in the search for these planets. The Bayesian estimates for the host mass, planet mass, system distance, and planet-host projected separation are $M_{\rm host}=0.76\pm 0.40 M_\odot$, $M_{\rm planet}=10.3\pm 5.5 M_\oplus$, $D_L = 3.3\pm 1.3$ kpc, and $a_\perp = 3.3\pm 1.4$ AU, respectively. The blended light is plausibly due to the host star, an identification that could be tested by high resolution images taken during the next (2020) bulge season.

Scaling relations trace the formation and evolution of galaxy clusters. We exploited multi-wavelength surveys -- the XXL survey at \emph{XMM-Newton} in the X-ray band, and the Hyper Suprime-Cam Subaru Strategic Program for optical weak lensing -- to study an X-ray selected, complete sample of clusters and groups. The scalings of gas mass, temperature, and soft-band X-ray luminosity with the weak lensing mass show imprints of radiative cooling and AGN feedback in groups. From the multi-variate analysis, we found some evidence for steeper than self-similar slopes for gas mass ($\beta_{m_\text{g}|m}=1.73 \pm0.80$) and luminosity ($\beta_{l|m}=1.91\pm0.94$) and a nearly self-similar slope for the temperature ($\beta_{t|m}=0.78\pm0.43$). Intrinsic scatters of X-ray properties appear to be positively correlated at a fixed mass (median correlation factor $\rho_{X_1X_2|m}\sim0.34$) due to dynamical state and merger history of the halos. Positive correlations with the weak lensing mass (median correlation factor $\rho_{m_\text{wl}X|m}\sim0.35$) can be connected to triaxiality and orientation. Comparison of weak lensing and hydrostatic masses suggests a small role played by non-thermal pressure support ($9\pm17\%$).

We present results from simultaneous modeling of high angular resolution GBT/MUSTANG-2 90 GHz Sunyaev-Zel'dovich effect (SZE) measurements and XMM-XXL X-ray images of three rich galaxy clusters selected from the HSC-SSP Survey. The combination of high angular resolution SZE and X-ray imaging enables a spatially resolved multi-component analysis, which is crucial to understand complex distributions of cluster gas properties. The targeted clusters have similar optical richnesses and redshifts, but exhibit different dynamical states in their member galaxy distributions; a single-peaked cluster, a double-peaked cluster, and a cluster belonging to a supercluster. A large-scale residual pattern in both regular Compton-parameter $y$ and X-ray surface brightness distributions is found in the single-peaked cluster, indicating a sloshing mode. The double-peaked cluster shows an X-ray remnant cool core between two SZE peaks associated with galaxy concentrations. The temperatures of the two peaks reach $\sim20-30$ keV in contrast to the cool core component of $\sim2$ keV, indicating a violent merger. The highest SZE signal for the supercluster is elongated along a direction perpendicular to the major axis of the X-ray core, suggesting a minor merger before core passage. The $S_X$ and $y$ distributions are thus perturbed at some level, regardless of the optical properties. We find that an integrated Compton $y$ parameter and a temperature for the major merger are significantly boosted from those expected by the weak-lensing mass and those for the other two clusters show no significant deviations, which is consistent with predictions of numerical simulations.

We report the discovery of a gas giant planet orbiting a low-mass host star in the microlensing event MOA-bin-29 that occurred in 2006. We find five degenerate solutions with the planet/host-star mass ratio of $q \sim 10^{-2}$. The Einstein radius crossing time of all models are relatively short ($\sim 4-7$ days), which indicates that the mass of host star is likely low. The measured lens-source proper motion is $5-9$ ${\rm mas}\ {\rm yr}^{-1}$ depending on the models. Since only finite source effects are detected, we conduct a Bayesian analysis in order to obtain the posterior probability distribution of the lens physical properties. As a result, we find the lens system is likely to be a gas giant orbiting a brown dwarf or a very late M-dwarf in the Galactic bulge. The probability distributions of the physical parameters for the five degenerate models are consistent within the range of error. By combining these probability distributions, we conclude that the lens system is a gas giant with a mass of $M_{\rm p} = 0.63^{+1.13}_{-0.39}\ M_{\rm Jup}$ orbiting a brown dwarf with a mass of $M_{\rm h} = 0.06^{+0.11}_{-0.04}\ M_\odot$ at a projected star-planet separation of $r_\perp = 0.53^{+0.89}_{-0.18}\ {\rm au}$. The lens distance is $D_{\rm L} = 6.89^{+1.19}_{-1.19}\ {\rm kpc}$, i.e., likely within the Galactic bulge.

Short-timescale microlensing events are likely to be produced by substellar brown dwarfs (BDs), but it is difficult to securely identify BD lenses based on only event timescales $t_{\rm E}$ because short-timescale events can also be produced by stellar lenses with high relative lens-source proper motions. In this paper, we report three strong candidate BD-lens events found from the search for lensing events not only with short timescales ($t_{\rm E} \lesssim 6~{\rm days}$) but also with very small angular Einstein radii ($\theta_{\rm E}\lesssim 0.05~{\rm mas}$) among the events that have been found in the 2016--2019 observing seasons. These events include MOA-2017-BLG-147, MOA-2017-BLG-241, and MOA-2019-BLG-256, in which the first two events are produced by single lenses and the last event is produced by a binary lens. From the Bayesian analysis conducted with the combined $t_{\rm E}$ and $\theta_{\rm E}$ constraint, it is estimated that the lens masses of the individual events are $0.051^{+0.100}_{-0.027}~M_\odot$, $0.044^{+0.090}_{-0.023}~M_\odot$, and $0.046^{+0.067}_{-0.023}~M_\odot/0.038^{+0.056}_{-0.019}~M_\odot$ and the probability of the lens mass smaller than the lower limit of stars is $\sim 80\%$ for all events. We point out that routine lens mass measurements of short time-scale lensing events require survey-mode space-based observations.

We present the first measurement of cross-correlation between the lensing potential, reconstructed from cosmic microwave background (CMB) {\it polarization} data, and the cosmic shear field from galaxy shapes. This measurement is made using data from the POLARBEAR CMB experiment and the Subaru Hyper Suprime-Cam (HSC) survey. By analyzing an 11~deg$^2$ overlapping region, we reject the null hypothesis at 3.5$\sigma$\ and constrain the amplitude of the {\bf cross-spectrum} to $\widehat{A}_{\rm lens}=1.70\pm 0.48$, where $\widehat{A}_{\rm lens}$ is the amplitude normalized with respect to the Planck~2018{} prediction, based on the flat $\Lambda$ cold dark matter cosmology. The first measurement of this {\bf cross-spectrum} without relying on CMB temperature measurements is possible due to the deep POLARBEAR map with a noise level of ${\sim}$6\,$\mu$K-arcmin, as well as the deep HSC data with a high galaxy number density of $n_g=23\,{\rm arcmin^{-2}}$. We present a detailed study of the systematics budget to show that residual systematics in our results are negligibly small, which demonstrates the future potential of this cross-correlation technique.

We present the first results of a pilot X-ray study of 37 rich galaxy clusters at $0.1<z<1.1$ in the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) field. Diffuse X-ray emissions from these clusters were serendipitously detected in the XMM-Newton fields of view. We systematically analyze X-ray images of 37 clusters and emission spectra of a subsample of 17 clusters with high photon statistics by using the XMM-Newton archive data. The frequency distribution of the offset between the X-ray centroid or peak and the position of the brightest cluster galaxy was derived for the optical cluster sample. The fraction of relaxed clusters estimated from the X-ray peak offsets in 17 clusters is $29\pm11(\pm13)$\%, which is smaller than that of the X-ray cluster samples such as HIFLUGCS. Since the optical cluster search is immune to the physical state of X-ray-emitting gas, it is likely to cover a larger range of the cluster morphology. We also derived the luminosity-temperature relation and found that the slope is marginally shallower than those of X-ray-selected samples and consistent with the self-similar model prediction of 2. Accordingly, our results show that the X-ray properties of the optical clusters are marginally different from those observed in the X-ray samples. The implication of the results and future prospects are briefly discussed.

We report the discovery of an exoplanet in microlensing event OGLE-2015-BLG-1649. The planet/host-star mass ratio is $q =7.2 \times 10^{-3}$ and the projected separation normalized by the Einstein radius is $s = 0.9$. The upper limit of the lens flux is obtained from adaptive optics observations by IRCS/Subaru, which excludes the probability of a G-dwarf or more massive host star and helps to put a tighter constraint on the lens mass as well as commenting on the formation scenarios of giant planets orbiting low-mass stars. We conduct a Bayesian analysis including constraints on the lens flux to derive the probability distribution of the physical parameters of the lens system. We thereby find that the masses of the host star and planet are $M_{L} = 0.34 \pm 0.19 M_{\odot}$ and $M_{p} = 2.5^{+1.5}_{-1.4} M_{Jup}$, respectively. The distance to the system is $D_{L} = 4.23^{+1.51}_{-1.64}$kpc. The projected star-planet separation is $a_{\perp} = 2.07^{+0.65}_{-0.77}$AU. The lens-source relative proper motion of the event is quite high, at $\sim 7.1 \, {\rm mas/yr}$. Therefore, we may be able to determine the lens physical parameters uniquely or place much stronger constraints on them by measuring the color-dependent centroid shift and/or the image elongation with additional high resolution imaging already a few years from now.

We report two microlensing planet candidates discovered by the KMTNet survey in $2017$. However, both events have the 2L1S/1L2S degeneracy, which is an obstacle to claiming the discovery of the planets with certainty unless the degeneracy can be resolved. For KMT-2017-BLG-0962, the degeneracy cannot be resolved. If the 2L1S solution is correct, KMT-2017-BLG-0962 might be produced by a super Jupiter-mass planet orbiting a mid-M dwarf host star. For KMT-2017-BLG-1119, the light curve modeling favors the 2L1S solution but higher-resolution observations of the baseline object tend to support the 1L2S interpretation rather than the planetary interpretation. This degeneracy might be resolved by a future measurement of the lens-source relative proper motion. This study shows the problem of resolving 2L1S/1L2S degeneracy exists over a much wider range of conditions than those considered by the theoretical study of Gaudi (1998).

The CFHT Large Area U-band Deep Survey (CLAUDS) uses data taken with the MegaCam mosaic imager on CFHT to produce images of 18.60 deg2 with median seeing of FWHM=0.92 arcsec and to a median depth of U = 27.1 AB (5 sigma in 2 arcsec apertures), with selected areas that total 1.36 deg2 reaching a median depth of U=27.7 AB. These are the deepest U-band images assembled to date over this large an area. These data are located in four fields also imaged to comparably faint levels in grizy and several narrowband filters as part of the Hyper Suprime-Cam (HSC) Subaru Strategic Program (HSC-SSP). These CFHT and Subaru datasets will remain unmatched in their combination of area and depth until the advent of the Large Synoptic Survey Telescope (LSST). This paper provides an overview of the scientific motivation for CLAUDS and gives details of the observing strategy, observations, data reduction, and data merging with the HSC-SSP. Three early applications of these deep data are used to illustrate the potential of the dataset: deep U-band galaxy number counts, z~3 Lyman break galaxy (LBG) selection, and photometric redshifts improved by adding CLAUDS U to the Subaru HSC grizy photometry.

We report the discovery of a planet in a binary that was discovered from the analysis of the microlensing event OGLE-2018-BLG-1700. We identify the triple nature of the lens from the fact that the complex anomaly pattern can be decomposed into two parts produced by two binary-lens events, in which one binary pair has a very low mass ratio of $\sim 0.01$ between the lens components and the other pair has a mass ratio of $\sim 0.3$. We find two sets of degenerate solutions, in which one solution has a projected separation between the primary and its stellar companion less than the angular Einstein radius $\thetae$ (close solution), while the other solution has a separation greater than $\thetae$ (wide solution). From the Bayesian analysis with the constraints of the event time scale and angular Einstein radius together with the location of the source lying in the far disk behind the bulge, we find that the planet is a super-Jupiter with a mass of $4.4^{+3.0}_{-2.0}~M_{\rm J}$ and the stellar binary components are early and late M-type dwarfs with masses $0.42^{+0.29}_{-0.19}~M_\odot$ and $0.12^{+0.08}_{-0.05}~M_\odot$, respectively, and the planetary system is located at a distance of $D_{\rm L}=7.6^{+1.2}_{-0.9}~{\rm kpc}$. The planet is a circumstellar planet according to the wide solution, while it is a circumbinary planet according to the close solution. The projected primary-planet separation is $2.8^{+3.2}_{-2.5}~{\rm au}$ commonly for the close and wide solutions, but the primary-secondary binary separation of the close solution, $0.75^{+0.87}_{-0.66}~{\rm au}$, is widely different from the separation, $10.5^{+12.1}_{-9.2}~{\rm au}$, of the wide solution.

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 present the analysis of the event OGLE-2017-BLG-1186 from the 2017 Spitzer microlensing campaign. This is a remarkable microlensing event because its source is photometrically bright and variable, which makes it possible to perform an asteroseismic analysis using ground-based data. We find that the source star is an oscillating red giant with average timescale of $\sim 9$ d. The asteroseismic analysis also provides us source properties including the source angular size ($\sim 27~\mu{\rm as}$) and distance ($\sim 11.5$ kpc), which are essential for inferring the properties of the lens. When fitting the light curve, we test the feasibility of Gaussian Processes (GPs) in handling the correlated noise caused by the variable source. We find that the parameters from the GP model are generally more loosely constrained than those from the traditional $\chi^2$ minimization method. We note that this event is the first microlensing system for which asteroseismology and GPs have been used to account for the variable source. With both finite-source effect and microlens parallax measured, we find that the lens is likely a $\sim 0.045~M_{\odot}$ brown dwarf at distance $\sim 9.0$ kpc, or a $\sim 0.073~M_{\odot}$ ultracool dwarf at distance $\sim 9.8$ kpc. Combining the estimated lens properties with a Bayesian analysis using a Galactic model, we find a $\sim 35$ per cent probability for the lens to be a bulge object and $\sim 65$ per cent to be a background disc object.

We report the discovery of the microlensing planet OGLE-2018-BLG-0740Lb. The planet is detected with a very strong signal of $\Delta\chi^2\sim 4630$, but the interpretation of the signal suffers from two types of degeneracies. One type is caused by the previously known close/wide degeneracy, and the other is caused by an ambiguity between two solutions, in which one solution requires to incorporate finite-source effects, while the other solution is consistent with a point-source interpretation. Although difficult to be firmly resolved based on only the photometric data, the degeneracy is resolved in strong favor of the point-source solution with the additional external information obtained from astrometric and spectroscopic observations. The small astrometric offset between the source and baseline object supports that the blend is the lens and this interpretation is further secured by the consistency of the spectroscopic distance estimate of the blend with the lensing parameters of the point-source solution. The estimated mass of the host is $1.0\pm 0.1~M_\odot$ and the mass of the planet is $4.5\pm 0.6~M_{\rm J}$ (close solution) or $4.8\pm 0.6~M_{\rm J}$ (wide solution) and the lens is located at a distance of $3.2\pm 0.5$~kpc. The bright nature of the lens, with $I\sim 17.1$ ($V\sim 18.2$), combined with its dominance of the observed flux suggest that radial-velocity (RV) follow-up observations of the lens can be done using high-resolution spectrometers mounted on large telescopes, e.g., VLT/ESPRESSO, and this can potentially not only measure the period and eccentricity of the planet but also probe for close-in planets. We estimate that the expected RV amplitude would be $\sim 60\sin i ~{\rm m~s}^{-1}$.

We report a multiplanetary system found from the analysis of microlensing event OGLE-2018-BLG-1011, for which the light curve exhibits a double-bump anomaly around the peak. We find that the anomaly cannot be fully explained by the binary-lens or binary-source interpretations and its description requires the introduction of an additional lens component. The 3L1S (3 lens components and a single source) modeling yields three sets of solutions, in which one set of solutions indicates that the lens is a planetary system in a binary, while the other two sets imply that the lens is a multiplanetary system. By investigating the fits of the individual models to the detailed light curve structure, we find that the multiple-planet solution with planet-to-host mass ratios $\sim 9.5\times 10^{-3}$ and $\sim 15\times 10^{-3}$ are favored over the other solutions. From the Bayesian analysis, we find that the lens is composed of two planets with masses $1.8^{+3..4}_{-1.1}~M_{\rm J}$ and $2.8^{+5.1}_{-1.7}~M_{\rm J}$ around a host with a mass $0.18^{+0.33}_{-0.10}~M_\odot$ and located at a distance $7.1^{+1.1}_{-1.5}~{\rm kpc}$. The estimated distance indicates that the lens is the farthest system among the known multiplanetary systems. The projected planet-host separations are $a_{\perp,2}=1.8^{+2.1}_{-1.5}~{\rm au}$ ($0.8^{+0.9}_{-0.6}~{\rm au}$) and $a_{\perp,3}=0.8^{+0.9}_{-0.6}~{\rm au}$, where the values of $a_{\perp,2}$ in and out the parenthesis are the separations corresponding to the two degenerate solutions, indicating that both planets are located beyond the snow line of the host, as with the other four multiplanetary systems previously found by microlensing.

We present the analysis of the microlensing event OGLE-2015-BLG-1670, detected in a high-extinction field, very close to the Galactic plane. Due to the dust extinction along the line of sight, this event was too faint to be detected before it reached the peak of magnification. The microlensing light-curve models indicate a high-magnification event with a maximum of $A_\mathrm{max}\gtrsim200$, very sensitive to planetary deviations. An anomaly in the light curve has been densely observed by the microlensing surveys MOA, KMTNet, and OGLE. From the light-curve modeling, we find a planetary anomaly characterized by a planet-to-host mass ratio, $q=\left(1.00^{+0.18}_{-0.16}\right)\times 10^{-4}$, at the peak recently identified in the mass-ratio function of microlensing planets. Thus, this event is interesting to include in future statistical studies about planet demography. We have explored the possible degeneracies and find two competing planetary models resulting from the $s\leftrightarrow1/s$ degeneracy. However, because the projected separation is very close to $s=1$, the physical implications for the planet for the two solutions are quite similar, except for the value of $s$. By combining the light-curve parameters with a Galactic model, we have estimated the planet mass $M_2=17.9^{+9.6}_{-8.8}\,\mathrm{M}_\oplus$ and the lens distance $D_\mathrm{L}=6.7^{+1.0}_{-1.3}\,\mathrm{kpc}$, corresponding to a Neptune-mass planet close to the Galactic bulge. Such events with a low absolute latitude ($|b|\approx 1.1\,\mathrm{deg}$) are subject to both high extinction and more uncertain source distances, two factors that may affect the mass measurements in the provisional Wide Field Infrared Survey Telescope fields. More events are needed to investigate the potential trade-off between the higher lensing rate and the difficulty in measuring masses in these low-latitude fields.

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 present the analysis of stellar binary microlensing event OGLE-2015-BLG-0060 based on observations obtained from 13 different telescopes. Intensive coverage of the anomalous parts of the light curve was achieved by automated follow-up observations from the robotic telescopes of the Las Cumbres Observatory. We show that, for the first time, all main features of an anomalous microlensing event are well covered by follow-up data, allowing us to estimate the physical parameters of the lens. The strong detection of second-order effects in the event light curve necessitates the inclusion of longer-baseline survey data in order to constrain the parallax vector. We find that the event was most likely caused by a stellar binary-lens with masses $M_{\star1} = 0.87 \pm 0.12 M_{\odot}$ and $M_{\star2} = 0.77 \pm 0.11 M_{\odot}$. The distance to the lensing system is 6.41 $\pm 0.14$ kpc and the projected separation between the two components is 13.85 $\pm 0.16$ AU. Alternative interpretations are also considered.

We report the discovery of a $Spitzer$ microlensing planet OGLE-2018-BLG-0596Lb, with preferred planet-host mass ratio $q \sim 2\times10^{-4}$. The planetary signal, which is characterized by a short $(\sim 1~{\rm day})$ "bump" on the rising side of the lensing light curve, was densely covered by ground-based surveys. We find that the signal can be explained by a bright source that fully envelops the planetary caustic, i.e., a "Hollywood" geometry. Combined with the source proper motion measured from $Gaia$, the $Spitzer$ satellite parallax measurement makes it possible to precisely constrain the lens physical parameters. The preferred solution, in which the planet perturbs the minor image due to lensing by the host, yields a Uranus-mass planet with a mass of $M_{\rm p} = 13.9\pm1.6~M_{\oplus}$ orbiting a mid M-dwarf with a mass of $M_{\rm h} = 0.23\pm0.03~M_{\odot}$. There is also a second possible solution that is substantially disfavored but cannot be ruled out, for which the planet perturbs the major image. The latter solution yields $M_{\rm p} = 1.2\pm0.2~M_{\oplus}$ and $M_{\rm h} = 0.15\pm0.02~M_{\odot}$. By combining the microlensing and $Gaia$ data together with a Galactic model, we find in either case that the lens lies on the near side of the Galactic bulge at a distance $D_{\rm L} \sim 6\pm1~{\rm kpc}$. Future adaptive optics observations may decisively resolve the major image/minor image degeneracy.

We report the mass and distance measurements of two single-lens events from the 2017 Spitzer microlensing campaign. The ground-based observations yield the detection of finite-source effects, and the microlens parallaxes are derived from the joint analysis of ground-based observations and Spitzer observations. We find that the lens of OGLE-2017-BLG-1254 is a $0.60 \pm 0.03 M_{\odot}$ star with $D_{\rm LS} = 0.53 \pm 0.11~\text{kpc}$, where $D_{\rm LS}$ is the distance between the lens and the source. The second event, OGLE-2017-BLG-1161, is subject to the known satellite parallax degeneracy, and thus is either a $0.51^{+0.12}_{-0.10} M_{\odot}$ star with $D_{\rm LS} = 0.40 \pm 0.12~\text{kpc}$ or a $0.38^{+0.13}_{-0.12} M_{\odot}$ star with $D_{\rm LS} = 0.53 \pm 0.19~\text{kpc}$. Both of the lenses are therefore isolated stars in the Galactic bulge. By comparing the mass and distance distributions of the eight published Spitzer finite-source events with the expectations from a Galactic model, we find that the Spitzer sample is in agreement with the probability of finite-source effects occurrence in single lens events.

In this work, we present the analysis of the binary microlensing event OGLE-2018-BLG-0022 that is detected toward the Galactic bulge field. The dense and continuous coverage with the high-quality photometry data from ground-based observations combined with the space-based {\it Spitzer} observations of this long time-scale event enables us to uniquely determine the masses $M_1=0.40 \pm 0.05~M_\odot$ and $M_2=0.13\pm 0.01~M_\odot$ of the individual lens components. Because the lens-source relative parallax and the vector lens-source relative proper motion are unambiguously determined, we can likewise unambiguously predict the astrometric offset between the light centroid of the magnified images (as observed by the {\it Gaia} satellite) and the true position of the source. This prediction can be tested when the individual-epoch {\it Gaia} astrometric measurements are released.

We analyze the binary microlensing event MOA-2016-BLG-231, which was observed from the ground and from Spitzer. The lens is composed of very low-mass brown dwarfs (BDs) with $M_1 = 21^{+12}_{-5} \ M_J$ and $M_2 = 9^{+5}_{-2}\ M_J$, and it is located in the Galactic disk $D_{\rm L} = 2.85^{+0.88}_{-0.50}\ {\rm kpc}$. This is the fifth binary brown dwarf discovered by microlensing, and the BD binary is moving counter to the orbital motion of disk stars. Constraints on the lens physical properties come from late time, non-caustic-crossing features of the Spitzer light curve. Thus, MOA-2016-BLG-231 shows how Spitzer plays a crucial role in resolving the nature of BDs in binary BD events with short timescale ($\lesssim 10$ days).

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.

Planet formation theories predict the existence of free-floating planets that have been ejected from their parent systems. Although they emit little or no light, they can be detected during gravitational microlensing events. Microlensing events caused by rogue planets are characterized by very short timescales $t_{\rm E}$ (typically below two days) and small angular Einstein radii $\theta_{\rm E}$ (up to several uas). Here we present the discovery and characterization of two ultra-short microlensing events identified in data from the Optical Gravitational Lensing Experiment (OGLE) survey, which may have been caused by free-floating or wide-orbit planets. OGLE-2012-BLG-1323 is one of the shortest events discovered thus far ($t_{\rm E}$=0.155 +/- 0.005 d, $\theta_{\rm E}$=2.37 +/- 0.10 uas) and was caused by an Earth-mass object in the Galactic disk or a Neptune-mass planet in the Galactic bulge. OGLE-2017-BLG-0560 ($t_{\rm E}$=0.905 +/- 0.005 d, $\theta_{\rm E}$=38.7 +/- 1.6 uas) was caused by a Jupiter-mass planet in the Galactic disk or a brown dwarf in the bulge. We rule out stellar companions up to a distance of 6.0 and 3.9 au, respectively. We suggest that the lensing objects, whether located on very wide orbits or free-floating, may originate from the same physical mechanism. Although the sample of ultrashort microlensing events is small, these detections are consistent with low-mass wide-orbit or unbound planets being more common than stars in the Milky Way.

We analyze the gravitational binary-lensing event OGLE-2016-BLG-0156, for which the lensing light curve displays pronounced deviations induced by microlens-parallax effects. The light curve exhibits 3 distinctive widely-separated peaks and we find that the multiple-peak feature provides a very tight constraint on the microlens-parallax effect, enabling us to precisely measure the microlens parallax $\pi_{\rm E}$. All the peaks are densely and continuously covered from high-cadence survey observations using globally located telescopes and the analysis of the peaks leads to the precise measurement of the angular Einstein radius $\theta_{\rm E}$. From the combination of the measured $\pi_{\rm E}$ and $\theta_{\rm E}$, we determine the physical parameters of the lens. It is found that the lens is a binary composed of two M dwarfs with masses $M_1=0.18\pm 0.01\ M_\odot$ and $M_2=0.16\pm 0.01\ M_\odot$ located at a distance $D_{\rm L}= 1.35\pm 0.09\ {\rm kpc}$. According to the estimated lens mass and distance, the flux from the lens comprises an important fraction, $\sim 25\%$, of the blended flux. The bright nature of the lens combined with the high relative lens-source motion, $\mu=6.94\pm 0.50\ {\rm mas}\ {\rm yr}^{-1}$, suggests that the lens can be directly observed from future high-resolution follow-up observations.

We compare the planet-to-star mass-ratio distribution measured by gravitational microlensing to core accretion theory predictions from population synthesis models. The core accretion theory's runaway gas accretion process predicts a dearth of intermediate-mass giant planets that is not seen in the microlensing results. In particular, the models predict $\sim10\,\times$ fewer planets at mass ratios of $10^{-4} \leq q \leq 4 \times 10^{-4}$ than inferred from microlensing observations. This tension implies that gas giant formation may involve processes that have hitherto been overlooked by existing core accretion models or that the planet-forming environment varies considerably as a function of host-star mass. Variation from the usual assumptions for the protoplanetary disk viscosity and thickness could reduce this discrepancy, but such changes might conflict with microlensing results at larger or smaller mass ratios, or with other observations. The resolution of this discrepancy may have important implications for planetary habitability because it has been suggested that the runaway gas accretion process may have triggered the delivery of water to our inner solar system. So, an understanding of giant planet formation may help us to determine the occurrence rate of habitable planets.

We present an analysis of the microlensing event OGLE-2015-BLG-0232. This event is challenging to characterize for two reasons. First, the light curve is not well sampled during the caustic crossing due to the proximity of the full Moon impacting the photometry quality. Moreover, the source brightness is difficult to estimate because this event is blended with a nearby K dwarf star. We found that the light curve deviations are likely due to a close brown dwarf companion (i.e., s = 0.55 and q = 0.06), but the exact nature of the lens is still unknown. We finally discuss the potential of follow-up observations to estimate the lens mass and distance in the future.

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 the analysis of the caustic-crossing binary microlensing event OGLE-2017-BLG-0039. Thanks to the very long duration of the event, with an event time scale $t_{\rm E}\sim 130$ days, the microlens parallax is precisely measured despite its small value of $\pie\sim 0.06$. The analysis of the well-resolved caustic crossings during both the source star's entrance and exit of the caustic yields the angular Einstein radius $\thetae\sim 0.6$~mas. The measured $\pie$ and $\thetae$ indicate that the lens is a binary composed of two stars with masses $\sim 1.0~M_\odot$ and $\sim 0.15~M_\odot$, and it is located at a distance of $\sim 6$ kpc. From the color and brightness of the lens estimated from the determined lens mass and distance, it is expected that $\sim 2/3$ of the $I$-band blended flux comes from the lens. Therefore, the event is a rare case of a bright lens event for which high-resolution follow-up observations can confirm the nature of the lens.

We present the analysis of the planetary microlensing event MOA-2016-BLG-319. The event light curve is characterized by a brief ($\sim 3$ days) anomaly near the peak produced by minor-image perturbations. From modeling, we find two distinct solutions that describe the observed light curve almost equally well. From the investigation of the lens-system configurations, we find that the confusion in the lensing solution is caused by the degeneracy between the two solutions resulting from the source passages on different sides of the planetary caustic. These degeneracies can be severe for major-image perturbations but it is known that they are considerably less severe for minor-image perturbations. From the comparison of the lens-system configuration with those of two previously discovered planetary events, for which similar degeneracies were reported, we find that the degeneracies are caused by the special source trajectories that passed the star-planet axes at approximately right angles. By conducting a Bayesian analysis, it is estimated that the lens is a planetary system in which a giant planet with a mass $M_{\rm p}=0.62^{+1.16}_{-0.33}~M_{\rm J}$ ($0.65^{+1.21}_{-0.35}~M_{\rm J}$) is orbiting a low-mass M-dwarf host with a mass $M_{\rm h}=0.15^{+0.28}_{-0.08}~M_\odot$. Here the planet masses in and out of the parentheses represent the masses for the individual degenerate solutions. The projected host-planet separations are $a_\perp\sim 0.95$ au and $\sim 1.05$ au for the two solutions. The identified degeneracy indicates the need to check similar degeneracies in future analyses of planetary lensing events with minor-image perturbations.

We present results from a set of simulations designed to constrain the weak lensing shear calibration for the Hyper Suprime-Cam (HSC) survey. These simulations include HSC observing conditions and galaxy images from the Hubble Space Telescope (HST), with fully realistic galaxy morphologies and the impact of nearby galaxies included. We find that the inclusion of nearby galaxies in the images is critical to reproducing the observed distributions of galaxy sizes and magnitudes, due to the non-negligible fraction of unrecognized blends in ground-based data, even with the excellent typical seeing of the HSC survey (0.58" in the $i$-band). Using these simulations, we detect and remove the impact of selection biases due to the correlation of weights and the quantities used to define the sample (S/N and apparent size) with the lensing shear. We quantify and remove galaxy property-dependent multiplicative and additive shear biases that are intrinsic to our shear estimation method, including a $\sim 10$ per cent-level multiplicative bias due to the impact of nearby galaxies and unrecognized blends. Finally, we check the sensitivity of our shear calibration estimates to other cuts made on the simulated samples, and find that the changes in shear calibration are well within the requirements for HSC weak lensing analysis. Overall, the simulations suggest that the weak lensing multiplicative biases in the first-year HSC shear catalog are controlled at the 1 per cent level.

We report the discovery and the analysis of the short timescale binary-lens microlensing event, MOA-2015-BLG-337. The lens system could be a planetary system with a very low mass host, around the brown dwarf/planetary mass boundary, or a brown dwarf binary. We found two competing models that explain the observed light curves with companion/host mass ratios of q~0.01 and ~0.17, respectively. From the measurement of finite source effects in the best-fit planetary model, we find a relatively small angular Einstein radius of theta_E ~ 0.03 mas which favors a low mass lens. We conduct a Bayesian analysis to obtain the probability distribution of the lens properties. The results for the planetary models strongly depend on the minimum mass, M_min, in the assumed mass function. In summary, there are two solutions of the lens system: (1) a brown dwarf/planetary mass boundary object orbited by a super-Neptune (the planetary model with M_min=0.001 M_sun) and (2) a brown dwarf binary (the binary model). If the planetary models is correct, this system can be one of a new class of planetary system, having a low host mass and also a planetary mass ratio (q <0.03) between the host and its companion. The discovery of the event is important for the study of planetary formation in very low mass objects. In addition, it is important to consider all viable solutions in these kinds of ambiguous events in order for the future comprehensive statistical analyses of planetary/binary microlensing events.

We present the analysis of planetary microlensing event MOA-2011-BLG-291, which has a mass ratio of $q=(3.8\pm0.7)\times10^{-4}$ and a source star that is redder (or brighter) than the bulge main sequence. This event is located at a low Galactic latitude in the survey area that is currently planned for NASA's WFIRST exoplanet microlensing survey. This unusual color for a microlensed source star implies that we cannot assume that the source star is in the Galactic bulge. The favored interpretation is that the source star is a lower main sequence star at a distance of $D_S=4.9\pm1.3\,$kpc in the Galactic disk. However, the source could also be a turn-off star on the far side of the bulge or a sub-giant in the far side of the Galactic disk if it experiences significantly more reddening than the bulge red clump stars. However, these possibilities have only a small effect on our mass estimates for the host star and planet. We find host star and planet masses of $M_{\rm host} =0.15^{+0.27}_{-0.10}M_\odot$ and $m_p=18^{+34}_{-12}M_\oplus$ from a Bayesian analysis with a standard Galactic model under the assumption that the planet hosting probability does not depend on the host mass or distance. However, if we attempt to measure the host and planet masses with host star brightness measurements from high angular resolution follow-up imaging, the implied masses will be sensitive to the host star distance. The WFIRST exoplanet microlensing survey is expected to use this method to determine the masses for many of the planetary systems that it discovers, so this issue has important design implications for the WFIRST exoplanet microlensing survey.

OGLE-2014-BLG-0962 (OB140962) is a stellar binary microlensing event that was well-covered by observations from the Spitzer satellite as well as ground-based surveys. Modelling yields a unique physical solution: a mid-M+M-dwarf binary with $M_{\rm prim} = 0.20 \pm 0.01 M_\odot$ and $M_{\rm sec} = 0.16 \pm 0.01 M_\odot$, with projected separation of $2.0 \pm 0.3$ AU. The lens is only $D_{\rm LS} = 0.41 \pm 0.06$ kpc in front of the source, making OB140962 a bulge lens and the most distant Spitzer binary lens to date. In contrast, because the Einstein radius ($\theta_{\rm E} = 0.143 \pm 0.007$ mas) is unusually small, a standard Bayesian analysis, conducted in the absence of parallax information, would predict a brown dwarf binary. We test the accuracy of Bayesian analysis over a set of Spitzer lenses, finding overall good agreement throughout the sample. We also illustrate the methodology for probing the Galactic distribution of planets by comparing the cumulative distance distribution of the Spitzer 2-body lenses to that of the Spitzer single lenses.

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 show that dense OGLE and KMTNet $I$-band survey data require four bodies (sources plus lenses) to explain the microlensing light curve of OGLE-2015-BLG-1459. However, these can equally well consist of three lenses and one source (3L1S), two lenses and two sources (2L2S) or one lens and three sources (1L3S). In the 3L1S and 2L2S interpretations, the host is a brown dwarf and the dominant companion is a Neptune-class planet, with the third body (in the 3L1S case) being a Mars-class object that could have been a moon of the planet. In the 1L3S solution, the light curve anomalies are explained by a tight (five stellar radii) low-luminosity binary source that is offset from the principal source of the event by $\sim 0.17\,\au$. These degeneracies are resolved in favor of the 1L3S solution by color effects derived from comparison to MOA data, which are taken in a slightly different ($R/I$) passband. To enable current and future ($WFIRST$) surveys to routinely characterize exomoons and distinguish among such exotic systems requires an observing strategy that includes both a cadence faster than 9 min$^{-1}$ and observations in a second band on a similar timescale.

We present weak-lensing measurements using the first-year data from the Hyper Suprime-Cam Strategic Survey Program on the Subaru telescope for eight galaxy clusters selected through their thermal Sunyaev-Zel'dovich (SZ) signal measured at 148 GHz with the Atacama Cosmology Telescope Polarimeter experiment. The overlap between the two surveys in this work is 33.8 square degrees, before masking bright stars. The signal-to-noise ratio of individual cluster lensing measurements ranges from 2.2 to 8.7, with a total of 11.1 for the stacked cluster weak-lensing signal. We fit for an average weak-lensing mass distribution using three different profiles, a Navarro-Frenk-White profile, a dark-matter-only emulated profile, and a full cosmological hydrodynamic emulated profile. We interpret the differences among the masses inferred by these models as a systematic error of 10\%, which is currently smaller than the statistical error. We obtain the ratio of the SZ-estimated mass to the lensing-estimated mass (the so-called hydrostatic mass bias $1-b$) of $0.74^{+0.13}_{-0.12}$, which is comparable to previous SZ-selected clusters from the Atacama Cosmology Telescope and from the {\sl Planck} Satellite. We conclude with a discussion of the implications for cosmological parameters inferred from cluster abundances compared to cosmic microwave background primary anisotropy measurements.

We measure and model the wavelength dependence of the PSF in the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP) survey. We find that PSF chromaticity is present in that redder stars appear smaller than bluer stars in the $g, r,$ and $i$-bands at the 1-2 per cent level and in the $z$ and $y$-bands at the 0.1-0.2 per cent level. From the color dependence of the PSF, we fit a model between the monochromatic PSF trace radius, $R$, and wavelength of the form $R(\lambda)\propto \lambda^{b}$. We find values of $b$ between -0.2 and -0.5, depending on the epoch and filter. This is consistent with the expectations of a turbulent atmosphere with an outer scale length of $\sim 10-100$ m, indicating that the atmosphere is dominating the chromaticity. We find evidence in the best seeing data that the optical system and detector also contribute some wavelength dependence. Meyers and Burchat (2015) showed that $b$ must be measured to an accuracy of $\sim 0.02$ not to dominate the systematic error budget of the Large Synoptic Survey Telescope (LSST) weak lensing (WL) survey. Using simple image simulations, we find that $b$ can be inferred with this accuracy in the $r$ and $i$-bands for all positions in the LSST field of view, assuming a stellar density of 1 star arcmin$^{-2}$ and that the optical PSF can be accurately modeled. Therefore, it is possible to correct for most, if not all, of the bias that the wavelength-dependent PSF will introduce into an LSST-like WL survey.

We report the discovery of a planetary system in which a super-earth orbits a late M-dwarf host. The planetary system was found from the analysis of the microlensing event OGLE-2017-BLG-0482, wherein the planet signal appears as a short-term anomaly to the smooth lensing light curve produced by the host. Despite its weak signal and short duration, the planetary signal was firmly detected from the dense and continuous coverage by three microlensing surveys. We find a planet/host mass ratio of $q\sim 1.4\times 10^{-4}$. We measure the microlens parallax $\pi_{\rm E}$ from the long-term deviation in the observed lensing light curve, but the angular Einstein radius $\theta_{\rm E}$ cannot be measured because the source trajectory did not cross the planet-induced caustic. Using the measured event timescale and the microlens parallax, we find that the masses of the planet and the host are $M_{\rm p}=9.0_{-4.5}^{+9.0}\ M_\oplus$ and $M_{\rm host}=0.20_{-0.10}^{+0.20}\ M_\odot$, respectively, and the projected separation between them is $a_\perp=1.8_{-0.7}^{+0.6}$ au. The estimated distance to the lens is $D_{\rm L}=5.8_{-2.1}^{+1.8}$ kpc. The discovery of the planetary system demonstrates that microlensing provides an important method to detect low-mass planets orbiting low-mass stars.

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 report the discovery of a sub-Jupiter mass planet orbiting beyond the snow line of an M-dwarf most likely in the Galactic disk as part of the joint Spitzer and ground-based monitoring of microlensing planetary anomalies toward the Galactic bulge. The microlensing parameters are strongly constrained by the light curve modeling and in particular by the Spitzer-based measurement of the microlens parallax, $\pi_\mathrm{E}$. However, in contrast to many planetary microlensing events, there are no caustic crossings, so the angular Einstein radius, $\theta_\mathrm{E}$ has only an upper limit based on the light curve modeling alone. Additionally, the analysis leads us to identify 8 degenerate configurations: the four-fold microlensing parallax degeneracy being doubled by a degeneracy in the caustic structure present at the level of the ground-based solutions. To pinpoint the physical parameters, and at the same time to break the parallax degeneracy, we make use of a series of arguments: the $\chi^2$ hierarchy, the Rich argument, and a prior Galactic model. The preferred configuration is for a host at $D_L=3.73_{-0.67}^{+0.66}~\mathrm{kpc}$ with mass $M_\mathrm{L}=0.30_{-0.12}^{+0.15}~\mathrm{M_\odot}$, orbited by a Saturn-like planet with $M_\mathrm{planet}=0.43_{-0.17}^{+0.21}~\mathrm{M_\mathrm{Jup}}$ at projected separation $a_\perp = 1.70_{-0.39}^{+0.38}~\mathrm{au}$, about 2.1 times beyond the system snow line. Therefore, it adds to the growing population of sub-Jupiter planets orbiting near or beyond the snow line of M-dwarfs discovered by microlensing. Based on the rules of the real-time protocol for the selection of events to be followed up with Spitzer, this planet will not enter the sample for measuring the Galactic distribution of planets.

We present the analysis of the binary-microlensing event OGLE-2014-BLG-0289. The event light curve exhibits very unusual five peaks where four peaks were produced by caustic crossings and the other peak was produced by a cusp approach. It is found that the quintuple-peak features of the light curve provide tight constraints on the source trajectory, enabling us to precisely and accurately measure the microlensing parallax $\pi_{\rm E}$. Furthermore, the three resolved caustics allow us to measure the angular Einstein radius $\thetae$. From the combination of $\pi_{\rm E}$ and $\thetae$, the physical lens parameters are uniquely determined. It is found that the lens is a binary composed of two M dwarfs with masses $M_1 = 0.52 \pm 0.04\ M_\odot$ and $M_2=0.42 \pm 0.03\ M_\odot$ separated in projection by $a_\perp = 6.4 \pm 0.5$ au. The lens is located in the disk with a distance of $D_{\rm L} = 3.3 \pm 0.3$~kpc. It turns out that the reason for the absence of a lensing signal in the {\it Spitzer} data is that the time of observation corresponds to the flat region of the light curve.

We perform a $z$-band survey for an optical counterpart of a binary neutron star coalescence GW170817 with Subaru/Hyper Suprime-Cam. Our untargeted transient search covers $23.6$ deg$^2$ corresponding to the $56.6\%$ credible region of GW170817 and reaches the $50\%$ completeness magnitude of $20.6$ mag on average. As a result, we find 60 candidates of extragalactic transients, including J-GEM17btc (a.k.a. SSS17a/DLT17ck). While J-GEM17btc is associated with NGC 4993 that is firmly located inside the 3D skymap of GW170817, the other 59 candidates do not have distance information in the GLADE v2 catalog or NASA/IPAC Extragalactic Database (NED). Among 59 candidates, 58 are located at the center of extended objects in the Pan-STARRS1 catalog, while one candidate has an offset. We present location, $z$-band apparent magnitude, and time variability of the candidates and evaluate the probabilities that they are located inside of the 3D skymap of GW170817. The probability for J-GEM17btc is $64\%$ being much higher than those for the other 59 candidates ($9.3\times10^{-3}-2.1\times10^{-1}\%$). Furthermore, the possibility, that at least one of the other 59 candidates is located within the 3D skymap, is only $3.2\%$. Therefore, we conclude that J-GEM17btc is the most-likely and distinguished candidate as the optical counterpart of GW170817.

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.

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 report the discovery of OGLE-2016-BLG-1190Lb, which is likely to be the first Spitzer microlensing planet in the Galactic bulge/bar, an assignation that can be confirmed by two epochs of high-resolution imaging of the combined source-lens baseline object. The planet's mass M_p= 13.4+-0.9 M_J places it right at the deuterium burning limit, i.e., the conventional boundary between "planets" and "brown dwarfs". Its existence raises the question of whether such objects are really "planets" (formed within the disks of their hosts) or "failed stars" (low mass objects formed by gas fragmentation). This question may ultimately be addressed by comparing disk and bulge/bar planets, which is a goal of the Spitzer microlens program. The host is a G dwarf M_host = 0.89+-0.07 M_sun and the planet has a semi-major axis a~2.0 AU. We use Kepler K2 Campaign 9 microlensing data to break the lens-mass degeneracy that generically impacts parallax solutions from Earth-Spitzer observations alone, which is the first successful application of this approach. The microlensing data, derived primarily from near-continuous, ultra-dense survey observations from OGLE, MOA, and three KMTNet telescopes, contain more orbital information than for any previous microlensing planet, but not quite enough to accurately specify the full orbit. However, these data do permit the first rigorous test of microlensing orbital-motion measurements, which are typically derived from data taken over <1% of an orbital period.

The unprecedented depth and area surveyed by the Subaru Strategic Program with the Hyper Suprime-Cam (HSC-SSP) have enabled us to construct and publish the largest distant cluster sample out to $z\sim 1$ to date. In this exploratory study of cluster galaxy evolution from $z=1$ to $z=0.3$, we investigate the stellar mass assembly history of brightest cluster galaxies (BCGs), evolution of stellar mass and luminosity distributions, stellar mass surface density profile, as well as the population of radio galaxies. Our analysis is the first high redshift application of the top N richest cluster selection, which is shown to allow us to trace the cluster galaxy evolution faithfully. Over the 230 deg$^2$ area of the current HSC-SSP footprint, selecting the top 100 clusters in each of the 4 redshift bins allows us to observe the buildup of galaxy population in descendants of clusters whose $z\approx 1$ mass is about $2\times 10^{14}\,M_\odot$. Our stellar mass is derived from a machine-learning algorithm, which is found to be unbiased and accurate with respect to the COSMOS data. We find very mild stellar mass growth in BCGs (about 35% between $z=1$ and 0.3), and no evidence for evolution in both the total stellar mass-cluster mass correlation and the shape of the stellar mass surface density profile. We also present the first measurement of the radio luminosity distribution in clusters out to $z\sim 1$, and show hints of changes in the dominant accretion mode powering the cluster radio galaxies at $z\sim 0.8$.

We have developed a system combining a back-illuminated Complementary-Metal-Oxide-Semiconductor (CMOS) imaging sensor and Xilinx Zynq System-on-Chip (SoC) device for a soft X-ray (0.5-10 keV) imaging spectroscopy observation of the Sun to investigate the dynamics of the solar corona. Because typical timescales of energy release phenomena in the corona span a few minutes at most, we aim to obtain the corresponding energy spectra and derive the physical parameters, i.e., temperature and emission measure, every few tens of seconds or less for future solar X-ray observations. An X-ray photon-counting technique, with a frame rate of a few hundred frames per second or more, can achieve such results. We used the Zynq SoC device to achieve the requirements. Zynq contains an ARM processor core, which is also known as the Processing System (PS) part, and a Programmable Logic (PL) part in a single chip. We use the PL and PS to control the sensor and seamless recording of data to a storage system, respectively. We aim to use the system for the third flight of the Focusing Optics Solar X-ray Imager (FOXSI-3) sounding rocket experiment for the first photon-counting X-ray imaging and spectroscopy of the Sun.

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

Recent detection of gravitational waves from a neutron star (NS) merger event GW170817 and identification of an electromagnetic counterpart provide a unique opportunity to study the physical processes in NS mergers. To derive properties of ejected material from the NS merger, we perform radiative transfer simulations of kilonova, optical and near-infrared emissions powered by radioactive decays of r-process nuclei synthesized in the merger. We find that the observed near-infrared emission lasting for > 10 days is explained by 0.03 Msun of ejecta containing lanthanide elements. However, the blue optical component observed at the initial phases requires an ejecta component with a relatively high electron fraction (Ye). We show that both optical and near-infrared emissions are simultaneously reproduced by the ejecta with a medium Ye of ~ 0.25. We suggest that a dominant component powering the emission is post-merger ejecta, which exhibits that mass ejection after the first dynamical ejection is quite efficient. Our results indicate that NS mergers synthesize a wide range of r-process elements and strengthen the hypothesis that NS mergers are the origin of r-process elements in the Universe.

The first detected gravitational wave from a neutron star merger was GW170817. In this study, we present J-GEM follow-up observations of SSS17a, an electromagnetic counterpart of GW170817. SSS17a shows a 2.5-mag decline in the $z$-band from 1.7 days to 7.7 days after the merger. Such a rapid decline is not comparable with supernovae light curves at any epoch. The color of SSS17a also evolves rapidly and becomes redder for later epochs; the $z-H$ color changed by approximately 2.5 mag in the period of 0.7 days to 7.7 days. The rapid evolution of both the optical brightness and the color are consistent with the expected properties of a kilonova that is powered by the radioactive decay of newly synthesized $r$-process nuclei. Kilonova models with Lanthanide elements can reproduce the aforementioned observed properties well, which suggests that $r$-process nucleosynthesis beyond the second peak takes place in SSS17a. However, the absolute magnitude of SSS17a is brighter than the expected brightness of the kilonova models with the ejecta mass of 0.01 $\Msun$, which suggests a more intense mass ejection ($\sim 0.03 \Msun$) or possibly an additional energy source.

Type Ia supernovae (SNe Ia) arise from the thermonuclear explosion of carbon-oxygen white dwarfs. Though the uniformity of their light curves makes them powerful cosmological distance indicators, long-standing issues remain regarding their progenitors and explosion mechanisms. Recent detection of the early ultraviolet pulse of a peculiar subluminous SN Ia has been claimed as new evidence for the companion-ejecta interaction through the single-degenerate channel. Here, we report the discovery of a prominent but red optical flash at $\sim$ 0.5 days after the explosion of a SN Ia which shows hybrid features of different SN Ia sub-classes: a light curve typical of normal-brightness SNe Ia, but with strong titanium absorptions, commonly seen in the spectra of subluminous ones. We argue that the early flash of such a hybrid SN Ia is different from predictions of previously suggested scenarios such as the companion-ejecta interaction. Instead it can be naturally explained by a SN explosion triggered by a detonation of a thin helium shell either on a near-Chandrasekhar-mass white dwarf ($\gtrsim$ 1.3 M$_{\odot}$) with low-yield $^{56}$Ni or on a sub-Chandrasekhar-mass white dwarf ($\sim$ 1.0 M$_{\odot}$) merging with a less massive white dwarf. This finding provides compelling evidence that one branch of the previously proposed explosion models, the helium-ignition scenario, does exist in nature, and such a scenario may account for explosions of white dwarfs in a wider mass range in contrast to what was previously supposed.

We present the result of microlensing event MOA-2016-BLG-290, which received observations from the two-wheel Kepler (K2), Spitzer, as well as ground-based observatories. A joint analysis of data from K2 and the ground leads to two degenerate solutions of the lens mass and distance. This degeneracy is effectively broken once the (partial) Spitzer light curve is included. Altogether, the lens is found to be an extremely low-mass star located in the Galactic bulge. MOA-2016-BLG-290 is the first microlensing event for which we have signals from three well-separated ($\sim1$ AU) locations. It demonstrates the power of two-satellite microlensing experiment in reducing the ambiguity of lens properties, as pointed out independently by S. Refsdal and A. Gould several decades ago.

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

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

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

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

Photometric redshifts are a key component of many science objectives in the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). In this paper, we describe and compare the codes used to compute photometric redshifts for HSC-SSP, how we calibrate them, and the typical accuracy we achieve with the HSC five-band photometry (grizy). We introduce a new point estimator based on an improved loss function and demonstrate that it works better than other commonly used estimators. We find that our photo-z's are most accurate at 0.2<~zphot<~1.5, where we can straddle the 4000A break. We achieve sigma(d_zphot/(1+zphot))~0.05 and an outlier rate of about 15% for galaxies down to i=25 within this redshift range. If we limit to a brighter sample of i<24, we achieve sigma~0.04 and ~8% outliers. Our photo-z's should thus enable many science cases for HSC-SSP. We also characterize the accuracy of our redshift probability distribution function (PDF) and discover that some codes over/under-estimate the redshift uncertainties, which have implications for N(z) reconstruction. Our photo-z products for the entire area in the Public Data Release 1 are publicly available, and both our catalog products (such as point estimates) and full PDFs can be retrieved from the data release site, https://hsc-release.mtk.nao.ac.jp/.

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 analysis of the binary gravitational microlensing event MOA-2015-BLG-020. The event has a fairly long timescale (about 63 days) and thus the light curve deviates significantly from the lensing model that is based on the rectilinear lens-source relative motion. This enables us to measure the microlensing parallax through the annual parallax effect. The microlensing parallax parameters constrained by the ground-based data are confirmed by the Spitzer observations through the satellite parallax method. By additionally measuring the angular Einstein radius from the analysis of the resolved caustic crossing, the physical parameters of the lens are determined. It is found that the binary lens is composed of two dwarf stars with masses $M_1 = 0.606 \pm 0.028M_\odot$ and $M_2 = 0.125 \pm 0.006M_\odot$ in the Galactic disk. Assuming the source star is at the same distance as the bulge red clump stars, we find the lens is at a distance $D_L = 2.44 \pm 0.10 kpc$. In the end, we provide a summary and short discussion of all published microlensing events in which the annual parallax effect is confirmed by other independent observations.

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.

The XXL survey currently covers two 25 sq. deg. patches with XMM observations of ~10ks. We summarise the scientific results associated with the first release of the XXL data set, that occurred mid 2016. We review several arguments for increasing the survey depth to 40 ks during the next decade of XMM operations. X-ray (z<2) cluster, (z<4) AGN and cosmic background survey science will then benefit from an extraordinary data reservoir. This, combined with deep multi-$\lambda$ observations, will lead to solid standalone cosmological constraints and provide a wealth of information on the formation and evolution of AGN, clusters and the X-ray background. In particular, it will offer a unique opportunity to pinpoint the z>1 cluster density. It will eventually constitute a reference study and an ideal calibration field for the upcoming eROSITA and Euclid missions.

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.

Machine learning techniques are widely applied in many modern optical sky surveys, e.q. Pan-STARRS1, PTF/iPTF and Subaru/Hyper Suprime-Cam survey, to reduce human intervention for data verification. In this study, we have established a machine learning based real-bogus system to reject the false detections in the Subaru/Hyper-Suprime-Cam StrategicSurvey Program (HSC-SSP) source catalog. Therefore the HSC-SSP moving object detection pipeline can operate more effectively due to the reduction of false positives. To train the real-bogus system, we use the stationary sources as the real training set and the "flagged" data as the bogus set. The training set contains 47 features, most of which are photometric measurements and shape moments generated from the HSC image reduction pipeline (hscPipe). Our system can reach a true positive rate (tpr) ~96% with a false positive rate (fpr) ~ 1% or tpr ~99% at fpr ~5%. Therefore we conclude that the stationary sources are decent real training samples, and using photometry measurements and shape moments can reject the false positives effectively.

Understanding how QSO's UV radiation affects galaxy formation is vital to our understanding of reionization era. Using a custom made narrow-band filter, $NB906$, on Subaru/Suprime-Cam, we investigated the number density of Ly$\alpha$ emitters (LAE) around a QSO at z=6.4. To date, this is the highest redshift narrow-band observation, where LAEs around a luminous QSO are investigated. Due to the large field-of-view of Suprime-Cam, our survey area is $\sim$5400~cMpc$^2$, much larger than previously studies at z=5.7 ($\sim$200 cMpc$^2$). In this field, we previously found a factor of 7 overdensity of Lyman break galaxies (LBGs). Based on this, we expected to detect $\sim$100 LAEs down to $NB906$=25 ABmag. However, our 6.4 hour exposure found none. The obtained upper limit on the number density of LAEs is more than an order lower than the blank fields. Furthermore, this lower density of LAEs spans a large scale of 10 $p$Mpc across. A simple argument suggests a strong UV radiation from the QSO can suppress star-formation in halos with $M_{vir}<10^{10}M_{\odot}$ within a $p$Mpc from the QSO, but the deficit at the edge of the field (5 $p$Mpc) remains to be explained.

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

In this paper, we describe the optical imaging data processing pipeline developed for the Subaru Telescope's Hyper Suprime-Cam (HSC) instrument. The HSC Pipeline builds on the prototype pipeline being developed by the Large Synoptic Survey Telescope's Data Management system, adding customizations for HSC, large-scale processing capabilities, and novel algorithms that have since been reincorporated into the LSST codebase. While designed primarily to reduce HSC Subaru Strategic Program (SSP) data, it is also the recommended pipeline for reducing general-observer HSC data. The HSC pipeline includes high level processing steps that generate coadded images and science-ready catalogs as well as low-level detrending and image characterizations.

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.

The Subaru Strategic Program (SSP) is an ambitious multi-band survey using the Hyper Suprime-Cam (HSC) on the Subaru telescope. The Wide layer of the SSP is both wide and deep, reaching a detection limit of i~26.0 mag. At these depths, it is challenging to achieve accurate, unbiased, and consistent photometry across all five bands. The HSC data are reduced using a pipeline that builds on the prototype pipeline for the Large Synoptic Survey Telescope. We have developed a Python-based, flexible framework to inject synthetic galaxies into real HSC images called SynPipe. Here we explain the design and implementation of SynPipe and generate a sample of synthetic galaxies to examine the photometric performance of the HSC pipeline. For stars, we achieve 1% photometric precision at i~19.0 mag and 6% precision at i~25.0 in the i-band. For synthetic galaxies with single-Sersic profiles, forced CModel photometry achieves 13% photometric precision at i~20.0 mag and 18% precision at i~25.0 in the i-band. We show that both forced PSF and CModel photometry yield unbiased color estimates that are robust to seeing conditions. We identify several caveats that apply to the version of HSC pipeline used for the first public HSC data release (DR1) that need to be taking into consideration. First, the degree to which an object is blended with other objects impacts the overall photometric performance. This is especially true for point sources. Highly blended objects tend to have larger photometric uncertainties, systematically underestimated fluxes and slightly biased colors. Second, >20% of stars at 22.5< i < 25.0 mag can be misclassified as extended objects. Third, the current CModel algorithm tends to strongly underestimate the half-light radius and ellipticity of galaxy with i>21.5 mag.

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 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 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 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 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 perform a weak-lensing study of the nearby cool-core galaxy clusters, Hydra A ($z=0.0538$) and A478 ($z=0.0881$), of which brightest cluster galaxies (BCGs) host powerful activities of active galactic nuclei (AGNs). For each cluster, the observed tangential shear profile is well described either by a single Navarro--Frenk--White model or a two-component model including the BCG as an unresolved point mass. For A478, we determine the BCG and its host-halo masses from a joint fit to weak-lensing and stellar photometry measurements. We find that the choice of initial mass functions (IMFs) can introduce a factor of two uncertainty in the BCG mass, whereas the BCG host halo mass is well constrained by data. We perform a joint analysis of weak-lensing and stellar kinematics data available for the Hydra A cluster, which allows us to constrain the central mass profile without assuming specific IMFs. We find that the central mass profile ($r<300$kpc) determined from the joint analysis is in excellent agreement with those from independent measurements, including dynamical masses estimated from the cold gas disk component, X-ray hydrostatic total mass estimates, and the central stellar mass estimated based on the Salpeter IMF. The observed dark-matter fraction around the BCG for Hydra A is found to be smaller than those predicted by adiabatic contraction models, suggesting the importance of other physical processes, such as the the AGN feedback and/or dissipationless mergers.

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

The Kiso Supernova Survey (KISS) is a high-cadence optical wide-field supernova (SN) survey. The primary goal of the survey is to catch the very early light of a SN, during the shock breakout phase. Detection of SN shock breakouts combined with multi-band photometry obtained with other facilities would provide detailed physical information on the progenitor stars of SNe. The survey is performed using a 2.2x2.2 deg field-of-view instrument on the 1.05-m Kiso Schmidt telescope, the Kiso Wide Field Camera (KWFC). We take a three-minute exposure in g-band once every hour in our survey, reaching magnitude g~20-21. About 100 nights of telescope time per year have been spent on the survey since April 2012. The number of the shock breakout detections is estimated to be of order of 1 during our 3-year project. This paper summarizes the KISS project including the KWFC observing setup, the survey strategy, the data reduction system, and CBET-reported SNe discovered so far by KISS.

We present our discovery of dramatic variability in SDSS J1100+4421 by the high-cadence transient survey Kiso Supernova Survey (KISS). The source brightened in the optical by at least a factor of three within about half a day. Spectroscopic observations suggest that this object is likely a narrow-line Seyfert 1 galaxy (NLS1) at z=0.840, however with unusually strong narrow emission lines. The estimated black hole mass of ~ 10^7 Msun implies bolometric nuclear luminosity close to the Eddington limit. SDSS J1100+4421 is also extremely radio-loud, with a radio loudness parameter of R ~ 4 x 10^2 - 3 x 10^3, which implies the presence of relativistic jets. Rapid and large-amplitude optical variability of the target, reminiscent of that found in a few radio- and gamma-ray loud NLS1s, is therefore produced most likely in a blazar-like core. The 1.4 GHz radio image of the source shows an extended structure with a linear size of about 100 kpc. If SDSS J1100+4421 is a genuine NLS1, as suggested here, this radio structure would then be the largest ever discovered in this type of active galaxies.

We conduct a joint X-ray and weak-lensing study of four relaxed galaxy clusters (Hydra A, A478, A1689 and A1835) observed by both Suzaku and Subaru out to virial radii, with an aim to understand recently-discovered unexpected feature of the ICM in cluster outskirts. We show that the average hydrostatic-to-lensing total mass ratio for the four clusters decreases from \sim 70% to \sim 40% as the overdensity contrast decreases from 500 to the virial value.The average gas mass fraction from lensing total mass estimates increases with cluster radius and agrees with the cosmic mean baryon fraction within the virial radius, whereas the X-ray-based gas fraction considerably exceeds the cosmic values due to underestimation of the hydrostatic mass. We also develop a new advanced method for determining normalized cluster radial profiles for multiple X-ray observables by simultaneously taking into account both their radial dependence and multivariate scaling relations with weak-lensing masses. Although the four clusters span a range of halo mass, concentration, X-ray luminosity and redshift, we find that the gas entropy, pressure, temperature and density profiles are all remarkably self-similar when scaled with the lensing M_200 mass and r_200 radius.The entropy monotonically increases out to \sim 0.5r_200 following the accretion shock heating model K(r)\propto r^1.1, and flattens at \simgt 0.5r_200.The universality of the scaled entropy profiles indicates that the thermalization mechanism over the entire cluster region (>0.1r_200) is controlled by gravitation in a common to all clusters, although the heating efficiency in the outskirts needs to be modified from the standard law.The bivariate scaling functions of the gas density and temperature reveal that the flattening of the outskirts entropy profile is caused by the steepening of the temperature, rather than the flattening of the gas density.

We present the results of an X-ray follow-up campaign targeting 10 Weak Lensing (WL) selected galaxy clusters from a Subaru WL survey. Eight clusters were studied with dedicated Chandra pointings, whereas archival X-ray data were used for the remaining two clusters. The WL clusters appear to fit the same scaling relation between X-ray luminosity and temperature as X-ray selected clusters. However, when we consider the luminosity-mass relation, the WL selected clusters appear underluminous by a factor 3.9$\pm$0.9 (or, alternatively, more massive by 2.9$\pm$0.2), compared to X-ray selected clusters. Only by considering various observational effects that could potentially bias WL masses, can this difference be reconciled. We used X-ray imaging data to quantify the dynamical state of the clusters and found that one of the clusters appears dynamically relaxed, and two of the clusters host a cool core, consistent with Sunyaev-Zel'dovich effect selected clusters. Our results suggest that regular, cool core clusters may be over-represented in X-ray selected samples.

We examined the anisotropic point spread function (PSF) of Suprime-Cam data utilizing dense star field data. We decomposed the PSF ellipticities into three components, the optical aberration, atmospheric turbulence, and chip-misalignment in an empirical manner, and evaluated the amplitude of each component. We found that, for long-exposure data, the optical aberration has the largest contribution to the PSF ellipticities, which could be modeled well by a simple analytic function based on the lowest-order aberration theory. The statistical properties of PSF ellipticities resulting from the atmospheric turbulence were investigated by using the numerical simulations. The simulation results are in a reasonable agreement with the observed data. It is also found that the optical PSF can be well corrected by the standard correction method with a polynomial fitting function. However, for the atmospheric PSF, its correction is affected by the common limitation caused by sparse sampling of PSFs due to a limited number of stars. We also examined the effects of the residual PSF anisotropies on Suprime-Cam cosmic shear data. We found that the shape and amplitude of the B-mode shear variance are broadly consistent with those of the residual PSF ellipticities measured from the dense star field data. This indicates that most of the sources of residual systematic are understood, which is an important step for cosmic shear statistics to be a practical tool of the precision cosmology.

Weak lensing provides an important route toward collecting samples of clusters of galaxies selected by mass. Subtle systematic errors in image reduction can compromise the power of this technique. We use the B-mode signal to quantify this systematic error and to test methods for reducing this error. We show that two procedures are efficient in suppressing systematic error in the B-mode: (1) refinement of the mosaic CCD warping procedure to conform to absolute celestial coordinates and (2) truncation of the smoothing procedure on a scale of 10$^{\prime}$. Application of these procedures reduces the systematic error to 20% of its original amplitude. We provide an analytic expression for the distribution of the highest peaks in noise maps that can be used to estimate the fraction of false peaks in the weak lensing $\kappa$-S/N maps as a function of the detection threshold. Based on this analysis we select a threshold S/N = 4.56 for identifying an uncontaminated set of weak lensing peaks in two test fields covering a total area of $\sim 3$deg$^2$. Taken together these fields contain seven peaks above the threshold. Among these, six are probable systems of galaxies and one is a superposition. We confirm the reliability of these peaks with dense redshift surveys, x-ray and imaging observations. The systematic error reduction procedures we apply are general and can be applied to future large-area weak lensing surveys. Our high peak analysis suggests that with a S/N threshold of 4.5, there should be only 2.7 spurious weak lensing peaks even in an area of 1000 deg$^2$ where we expect $\sim$ 2000 peaks based on our Subaru fields.

We present a Subaru weak lensing measurement of ACT-CL J0022.2-0036, one of the most luminous, high-redshift (z=0.81) Sunyaev-Zel'dovich (SZ) clusters discovered in the 268 deg^2 equatorial region survey of the Atacama Cosmology Telescope. For the weak lensing analysis using i'-band images, we use a model-fitting (Gauss-Laguerre shapelet) method to measure shapes of galaxy images, where we fit galaxy images in different exposures simultaneously to obtain best-fit ellipticities taking into account the different PSFs in each exposure. We also take into account the astrometric distortion effect on galaxy images by performing the model fitting in the world coordinate system. To select background galaxies behind the cluster at z=0.81, we use photometric redshift (photo-z) estimates for every galaxy derived from the co-added images of multi-passband Br'i'z'Y, with PSF matching/homogenization. After a photo-z cut for background galaxy selection, we detect the tangential weak lensing distortion signal with a total signal-to-noise ratio of about 3.7. By fitting a Navarro-Frenk-White model to the measured shear profile, we find the cluster mass to be M_200\bar{\rho}_m = [7.5^+3.2_-2.8(stat.)^+1.3_-0.6(sys.)] x 10^14 M_\odot/h. The weak lensing-derived mass is consistent with previous mass estimates based on the SZ observation, with assumptions of hydrostatic equilibrium and virial theorem, as well as with scaling relations between SZ signal and mass derived from weak lensing, X-ray, and velocity dispersion, within the measurement errors. We also show that the existence of ACT-CL J0022.2-0036 at z=0.81 is consistent with the cluster abundance prediction of the \Lambda-dominated cold dark matter structure formation model. We thus demonstrate the capability of Subaru-type ground-based images for studying weak lensing of high-redshift clusters.

We present a method for attaining sub-arcsecond pointing stability during sub- orbital balloon flights, as designed for in the High Altitude Lensing Observatory (HALO) concept. The pointing method presented here has the potential to perform near-space quality optical astronomical imaging at 1-2% of the cost of space-based missions. We also discuss an architecture that can achieve sufficient thermomechanical stability to match the pointing stability. This concept is motivated by advances in the development and testing of Ultra Long Duration Balloon (ULDB) flights which promise to allow observation campaigns lasting more than three months. The design incorporates a multi-stage pointing architecture comprising: a gondola coarse azimuth control system, a multi-axis nested gimbal frame structure with arcsecond stability, a telescope de-rotator to eliminate field rotation, and a fine guidance stage consisting of both a telescope mounted angular rate sensor and guide CCDs in the focal plane to drive a fast-steering mirror. We discuss the results of pointing tests together with a preliminary thermo-mechanical analysis required for sub-arcsecond pointing at high altitude. Possible future applications in the areas of wide-field surveys and exoplanet searches are also discussed.

We use a dense redshift survey in the foreground of the Subaru GTO2deg^2 weak lensing field (centered at $\alpha_{2000}$ = 16$^h04^m44^s$;$\delta_{2000}$ =43^\circ11^{\prime}24^{\prime\prime}$) to assess the completeness and comment on the purity of massive halo identification in the weak lensing map. The redshift survey (published here) includes 4541 galaxies; 4405 are new redshifts measured with the Hectospec on the MMT. Among the weak lensing peaks with a signal-to-noise greater that 4.25, 2/3 correspond to individual massive systems; this result is essentially identical to the Geller et al. (2010) test of the Deep Lens Survey field F2. The Subaru map, based on images in substantially better seeing than the DLS, enables detection of less massive halos at fixed redshift as expected. We demonstrate that the procedure adopted by Miyazaki et al. (2007) for removing some contaminated peaks from the weak lensing map improves agreement between the lensing map and the redshift survey in the identification of candidate massive systems.

We present DLSCL J0916.2+2951 (z=0.53), a newly discovered major cluster merger in which the collisional cluster gas has become dissociated from the collisionless galaxies and dark matter. We identified the cluster using optical and weak lensing observations as part of the Deep Lens Survey. Our follow-up observations with Keck, Subaru, Hubble Space Telescope, and Chandra show that the cluster is a dissociative merger and constrain the dark matter self-interaction cross-section {\sigma}_{DM}m_{DM}^{-1}\leq7 cm^2g^{-1}. The system is observed at least 0.7\pm0.2 Gyr since first pass-through, thus providing a picture of cluster mergers 2-5 times further progressed than similar systems observed to date. This improved temporal leverage has implications for our understanding of merging clusters and their impact on galaxy evolution.

We have taken a deep, moderate-resolution Keck/Deimos spectra of QSO, CFHQS2329, at z=6.4. At the wavelength of Lya, the spectrum shows a spatially-extended component, which is significantly more extended than a stellar spectrum, and also a continuum part of the spectrum. The restframe line width of the extended component is 21+-7 A, and thus smaller than that of QSO (52+-4 A), where they should be identical if the light is incomplete subtraction of the QSO component. Therefore, these comparisons argue for the detection of a spatially extended Lya nebulae around this QSO. This is the first z>6 QSO that an extended Lya halo has been observed around. Careful subtraction of the central QSO spectrum reveals a lower limit to the Lya luminosity of (1.7+-0.1)x 10^43 erg s^-1. This emission may be from the theoretically predicted infalling gas in the process of forming a primordial galaxy that is ionized by a central QSO. On the other hand, if it is photoionized by the host galaxy, an estimated star-formation rate of >3.0 Msun yr^-1 is required. If we assume the gas is virialized, we obtain dynamical mass estimate of Mdyn=1.2x10^12 Msun. The derived MBH/Mhost is 2.1x10^-4, which is two orders smaller than those from more massive z~6 QSOs, and places this galaxy in accordance with the local M-sigma relation, in contrast to a previous claim on the evolution of M-sigma relation at z~6. We do not claim evolution or non-evolution of the M-sigma relation based on a single object, but our result highlights the importance of investigating fainter QSOs at z~6.

Understanding the cosmic re-ionization is one of the key goals of the modern observational cosmology. High redshift QSO spectra can be used as background light sources to measure absorption by intervening neutral hydrogen. We investigate neutral hydrogen absorption in a deep, moderate-resolution Keck/Deimos spectrum of QSO CFHQSJ2329-0301 at z=6.4. This QSO is one of the highest redshift QSOs presently known at z=6.4 but is 2.5 mag fainter than a previously well-studied QSO SDSSJ1148+5251 at z=6.4. Therefore, it has a smaller Stromgren sphere, and allows us to probe the highest redshift hydrogen absorption to date. The average transmitted flux at 5.915<z_abs<6.365 (200 comoving Mpc) is consistent with zero, in Ly_alpha, Ly_beta, and Ly_gamma absorption measurements. This corresponds to the lower limit of optical depth, tau_eff>4.9. These results are consistent with strong evolution of the optical depth at z>5.7.

Based on the ultraviolet to far-infrared photometry already compiled and presented in a companion paper (Barro et al. 2011a, Paper I), we present a detailed SED analysis of nearly 80,000 IRAC 3.6+4.5 micron selected galaxies in the Extended Groth Strip. We estimate photometric redshifts, stellar masses, and star formation rates separately for each galaxy in this large sample. The catalog includes 76,936 sources with [3.6] < 23.75 (85% completeness level of the IRAC survey) over 0.48 square degrees. The typical photometric redshift accuracy is Delta z/(1+z)=0.034, with a catastrophic outlier fraction of just 2%. We quantify the systematics introduced by the use of different stellar population synthesis libraries and IMFs in the calculation of stellar masses. We find systematic offsets ranging from 0.1 to 0.4 dex, with a typical scatter of 0.3 dex. We also provide UV- and IR-based SFRs for all sample galaxies, based on several sets of dust emission templates and SFR indicators. We evaluate the systematic differences and goodness of the different SFR estimations using the deep FIDEL 70 micron data available in the EGS. Typical random uncertainties of the IR-bases SFRs are a factor of two, with non-negligible systematic effects at z$\gtrsim$1.5 observed when only MIPS 24 micron data is available. All data products (SEDs, postage stamps from imaging data, and different estimations of the photometric redshifts, stellar masses, and SFRs of each galaxy) described in this and the companion paper are publicly available, and they can be accessed through our the web-interface utility Rainbow-navigator

We present an IRAC 3.6+4.5 microns selected catalog in the Extended Groth Strip (EGS) containing photometry from the ultraviolet to the far-infrared and stellar parameters derived from the analysis of the multi-wavelength data. In this paper, we describe the method used to build coherent spectral energy distributions (SEDs) for all the sources. In a companion paper, we analyze those SEDs to obtain robust estimations of stellar parameters such as photometric redshifts, stellar masses, and star formation rates. The catalog comprises 76,936 sources with [3.6]<23.75 mag (85% completeness level of the IRAC survey in the EGS) over 0.48 square degrees. For approximately 16% of this sample, we are able to deconvolve the IRAC data to obtain robust fluxes for the multiple counterparts found in ground-based optical images. Typically, the SEDs of the IRAC sources in our catalog count with more than 15 photometric data points, spanning from the UV to the FIR. Approximately 95% and 90% of all IRAC sources are detected in the deepest optical and near-infrared bands. Only 10% of the sources have optical spectroscopy and redshift estimations. Almost 20% and 2% of the sources are detected by MIPS at 24 and 70 microns, respectively. We also cross-correlate our catalog with public X-ray and radio catalogs. Finally, we present the Rainbow Navigator public web-interface utility designed to browse all the data products resulting from this work, including images, spectra, photometry, and stellar parameters.

Galaxy clusters are usually detected in blind optical surveys via suitable filtering methods. We present an optimal matched filter which maximizes their signal-to-noise ratio by taking advantage of the knowledge we have of their intrinsic physical properties and of the data noise properties. In this paper we restrict our application to galaxy magnitudes, positions and photometric redshifts if available, and we also apply the filter separately to weak lensing data. The method is suitable to be naturally extended to a multi-band approach which could include not only additional optical bands but also observables with different nature such as X-rays. For each detection, the filter provides its significance, an estimate for the richness and for the redshift even if photo-z are not given. The provided analytical error estimate is tested against numerical simulations. We finally apply our method to the COSMOS field and compare the results with previous cluster detections obtained with different methods. Our catalogue contains 27 galaxy clusters with minimal threshold at 3-sigma level including both optical and weak-lensing information.

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.

This paper analyzes a sample of 489 Spitzer/IRAC sources in the Extended Groth Strip whose spectral energy distributions fit a red power law from 3.6 to 8.0 \micron. The median for sources with known redshift is <z>=1.6. Though all or nearly all of the sample are likely to be active galactic nuclei, only 33% were detected in the EGS X-ray survey (AEGIS-X) using 200 ks Chandra observations. The detected sources are X-ray luminous with L_X > 10^43 erg/s and moderately to heavily obscured with N_H > 10^22 cm^-2. Stacking the X-ray-undetected sample members yields a statistically significant X-ray signal, suggesting that they are on average more distant or more obscured than sources with X-ray detections. The ratio of X-ray to mid-infrared fluxes suggests that a substantial fraction of the sources undetected in X-rays are obscured at the Compton-thick level, in contrast to the X-ray-detected sources, all of which appear to be Compton-thin. For the X-ray-detected power-law sources with redshifts, an X-ray luminosity L_X ~ 10^44 erg/s marks a transition between low-luminosity, blue sources dominated by the host galaxy to high-luminosity, red power-law sources dominated by nuclear activity. X-ray-to-optical ratios, infrared variability, and 24 micron properties of the sample are consistent with the identification of infrared power-law sources as active nuclei, but a rough estimate is that only 22% of AGNs are selected by the power law criteria. Comparison of the power-law selection technique and various IRAC color criteria for identifying AGNs confirms that high-redshift samples selected via simple IRAC colors may be heavily contaminated by starlight-dominated objects.

(Abridged) We present results of Suzaku observations of the intracluster medium (ICM) in Abell 1689, combined with complementary analysis of the SDSS data and weak and strong lensing analysis of Subaru/Suprime-Cam and HST/ACS observations. Faint X-ray emission from the ICM around the virial radius is detected at 4.0 sigma significance. We find anisotropic gas temperature and entropy distributions in cluster outskirts correlated with large-scale structure of galaxies. The high temperature and entropy region in the northeastern (NE) outskirts is connected to an overdense filamentary structure. The outskirt regions in contact with low density void environments have low gas temperatures and entropies, deviating from hydrostatic equilibrium. These results suggest that thermalization of the ICM occurs faster along the filamentary structures than the void regions. A joint X-ray and lensing analysis shows that the hydrostatic mass is $\sim60-90%$ of spherical lensing one but comparable to a triaxial halo mass within errors in $0.6r_{2500} \simlt r \simlt 0.8r_{500}$, and that it is significantly biased as low as $\simlt60%$ within $0.4r_{2500}$, irrespective of mass models. The thermal gas pressure within $r_{500}$ is, at most, $\sim50$--60% of the total pressure to balance fully the gravity of the spherical lensing mass, and $\sim30$--40% around the virial radius. Although these constitute lower limits when one considers the possible halo triaxiality, these small relative contributions of thermal pressure would require additional sources of pressure, such as bulk and/or turbulent motions.

We present a detailed analysis of the baryonic and dark matter distribution in the lensing cluster Abell 611 (z=0.288), with the goal of determining the dark matter profile over an unprecedented range of cluster-centric distance. By combining three complementary probes of the mass distribution, weak lensing from deep multi-color imaging, strong lensing constraints based on the identification of multiply-imaged sources, and resolved stellar velocity dispersion measures for the brightest cluster galaxy (BCG), we extend the methodology for separating the dark and baryonic mass components introduced by Sand et al. (2008). Our resulting dark matter profile samples the cluster from ~3 kpc to 3.25 Mpc, thereby providing an excellent basis for comparisons with recent numerical models. We demonstrate that only by combining our three observational techniques can degeneracies in constraining the form of the dark matter profile be broken on scales crucial for detailed comparisons with numerical simulations. Our analysis reveals that a simple Navarro, Frenk, and White (NFW) profile is an unacceptable fit to our data. We confirm earlier claims that the inner profile of the dark matter profile deviates significantly from the NFW form and find a inner logarithmic slope \beta flatter than 0.3 (68%; where rho_DM ~ r^{-\beta} at small radii). In order to reconcile our data with cluster formation in a LambdaCDM cosmology, we speculate that it may be necessary to revise our understanding of the nature of baryon--dark matter interactions in cluster cores. Comprehensive weak and strong lensing data, when coupled with kinematic information on the brightest cluster galaxy, can readily be applied to a larger sample of clusters to test the universality of these results.

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 present early phase observations in optical and near-infrared wavelengths for the extremely luminous Type Ia supernova (SN Ia) 2009dc. The decline rate of the light curve is $\Delta m_{15}(B)=0.65\pm 0.03$, which is one of the slowest among SNe Ia. The peak $V$-band absolute magnitude is $M_{V}=-19.90\pm 0.15$ mag even if the host extinction is $A_{V}=0$ mag. It reaches $M_{V}=-20.19\pm 0.19$ mag for the host extinction of $A_{V}=0.29$ mag as inferred from the observed Na {\sc i} D line absorption in the host. Our $JHK_{s}$-band photometry shows that the SN is one of the most luminous SNe Ia also in near-infrared wavelengths. These results indicate that SN 2009dc belongs to the most luminous class of SNe Ia, like SN 2003fg and SN 2006gz. We estimate the ejected $^{56}$Ni mass of $1.2\pm 0.3$ $\Msun$ for no host extinction case (or 1.6$\pm$ 0.4 M$_{\odot}$ for the host extinction of $A_{V}=0.29$ mag). The C {\sc ii} $\lambda$6580 absorption line keeps visible until a week after maximum, which diminished in SN 2006gz before its maximum brightness. The line velocity of Si {\sc ii} $\lambda$6355 is about 8000 km s$^{-1}$ around the maximum, being considerably slower than that of SN 2006gz, while comparable to that of SN 2003fg. The velocity of the C {\sc ii} line is almost comparable to that of the Si {\sc ii}. The presence of the carbon line suggests that thick unburned C+O layers remain after the explosion. SN 2009dc is a plausible candidate of the super-Chandrasekhar mass SNe Ia.

We present $u'g'R$ optical images taken with the MMT/Megacam and the Subaru/Suprime of the Extended Groth Strip survey. The total survey covers an area of about $\sim 1$ degree$^2$, including four sub-fields and is optimized for the study of galaxies at $z\sim3$. Our methods for photometric calibration in AB magnitudes, the limiting magnitude and the galaxy number count are described. A sample of 1642 photometrically selected candidate LBGs to an apparent $R_{AB}$ magnitude limit of 25.0 is present. The average sky surface density of our LBGs sample is $\sim$ 1.0 arcmin$^{-2}$, slightly higher than the previous finding.

We present the first results of a MOS campaign to follow up cluster candidates located via weak lensing. Our main goals are to search for spatial concentrations of galaxies that are plausible optical counterparts of the weak lensing signals, and to determine the cluster redshifts from those of member galaxies. Around each of 36 targeted cluster candidates, we obtain 15-32 galaxy redshifts. For 28 of these targets, we confirm a secure cluster identification. This includes three cases where two clusters at different redshifts are projected along the same line-of-sight. In 6 of the 8 unconfirmed targets, we find multiple small galaxy concentrations at different redshifts. In both the remaining two targets, a single small galaxy concentration is found. We evaluate the weak lensing mass of confirmed clusters. For a subsample of our most cleanly measured clusters, we investigate the statistical relation between their weak lensing mass and the velocity dispersion of their member galaxies, comparing our sample with optically and X-ray selected samples from the literature. Our lensing-selected clusters are consistent with sigma_v=sigma_sis, with a similar scatter to the optically and X-ray selected clusters. We thus find no evidence of selection bias compared to these other techniques. We also derive an empirical relation between the cluster mass and the galaxy velocity dispersion, which is in reasonable agreement with the prediction of N-body simulations in the LCDM cosmology.

We qualitatively examine properties of artificial deformation in shapes of objects (galaxies and stars) induced by the pixelization effects (also called as the aliasing effects) using toy mock simulation images. We pay a special attention to the the second pixelization because it might be a potential source of a systematic noise in a weak lensing analysis. In particular, it is found that resampling with rotation induces artificial ellipticities in object shapes having a periodic concentric-circle-shaped pattern. Our major findings are as follows. (1) Root-mean-square (RMS) of artificial ellipticities in object shapes induced by the first pixelization effect can be as large as RMS > 10^-2 if a characteristic size of objects (e.g., the FWHM) is smaller than twice of the pixel size. While for larger objects, it quickly becomes very small (RMS < 10^-5). (2) The amplitude of the shape deformation induced by the second pixelization effect depends on the object size. It also depends strongly on an interpolation scheme adopted to carry out resampling and on the grid size of the output pixels. The RMS of ellipticities in object shapes induced by the second pixelization effect can be suppressed to well below 10^{-2} if one adopts a proper interpolation scheme (implemented in popular image processing softwares). We also discuss an impact of the pixelization effects on a weak lensing analysis.

The application of multi-wavelength selection techniques is crucial for discovering a complete and unbiased set of Active Galactic Nuclei (AGNs). Here, we select a sample of 72 AGN candidates in the Extended Groth Strip (EGS) using deep radio and mid-infrared data from the Very Large Array (VLA) and the Spitzer Space Telescope, and analyze their properties across other wavelengths. Only 30% of these sources are detected in deep 200 ks Chandra X-ray Observatory pointings. The X-ray detected sources demonstrate moderate obscuration with column densities of N_H > 10^22 cm^-2. A stacked image of sources undetected by Chandra shows low levels of X-ray activity, suggesting they may be faint or obscured AGNs. Less than 40% of our sample are selected as AGNs with optical broad lines, mid-infrared power laws, or X-ray detections. Thus, if our candidates are indeed AGNs, the radio/mid-infrared selection criteria we use provide a powerful tool for identifying sources missed by other surveys.

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 assess the relative merits of weak lensing surveys, using overlapping imaging data from the ground-based Subaru telescope and the Hubble Space Telescope (HST). Our tests complement similar studies undertaken with simulated data. From observations of 230,000 matched objects in the 2 square degree COSMOS field, we identify the limit at which faint galaxy shapes can be reliably measured from the ground. Our ground-based shear catalog achieves sub-percent calibration bias compared to high resolution space-based data, for galaxies brighter than i'~24.5 and with half-light radii larger than 1.8". This selection corresponds to a surface density of ~15 galaxies per sq arcmin compared to ~71 per sq arcmin from space. On the other hand the survey speed of current ground-based facilities is much faster than that of HST, although this gain is mitigated by the increased depth of space-based imaging desirable for tomographic (3D) analyses. As an independent experiment, we also reconstruct the projected mass distribution in the COSMOS field using both data sets, and compare the derived cluster catalogs with those from X-ray observations. The ground-based catalog achieves a reasonable degree of completeness, with minimal contamination and no detected bias, for massive clusters at redshifts 0.2<z<0.5. The space-based data provide improved precision and a greater sensitivity to clusters of lower mass or at higher redshift.

We present the results of an ongoing weak lensing survey conducted with the Subaru telescope whose initial goal is to locate and study the distribution of shear-selected structures or halos. Using a Suprime-cam imaging survey spanning 21.82 square degree, we present a catalog of 100 candidate halos located from lensing convergence maps. Our sample is reliably drawn from that subset of our survey area, (totaling 16.72 square degree) uncontaminated by bright stars and edge effects and limited at a convergence signal to noise ratio of 3.69. To validate the sample detailed spectroscopic measures have been made for 26 candidates using the Subaru multi-object spectrograph, FOCAS. All are confirmed as clusters of galaxies but two arise as the superposition of multiple clusters viewed along the line of sight. Including data available in the literature and an ongoing Keck spectroscopic campaign, a total of 41 halos now have reliable redshifts. For one of our survey fields, the XMM LSS (Pierre et al. 2004) field, we compare our lensing-selected halo catalog with its X-ray equivalent. Of 15 halos detected in the XMM-LSS field, 10 match with published X-ray selected clusters and a further 2 are newly-detected and spectroscopically confirmed in this work. Although three halos have not yet been confirmed, the high success rate within the XMM-LSS field (12/15) confirms that weak lensing provides a reliable method for constructing cluster catalogs, irrespective of the nature of the constituent galaxies or the intracluster medium.

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

The Cosmic Evolution Survey (COSMOS) was initiated with an extensive allocation (590 orbits in Cycles 12-13) using the Hubble Space Telescope (HST) for high resolution imaging. Here we review the characteristics of the HST imaging with the Advanced Camera for Surveys (ACS) and parallel observations with NICMOS and WFPC2. A square field (1.8$\sq$\deg) has been imaged with single-orbit ACS I-F814W exposures with 50% completeness for sources 0.5\arcsec in diameter at I$_{AB} $ = 26.0 mag. The ACS imaging is a key part of the COSMOS survey, providing very high sensitivity and high resolution (0.09\arcsec FWHM, 0.05\arcsec pixels) imaging and detecting 1.2 million objects to a limiting magnitude of 26.5 (AB). These images yield resolved morphologies for several hundred thousand galaxies. The small HST PSF also provides greatly enhanced sensitivity for weak lensing investigations of the dark matter distribution.

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.

We study 87 Extremely Red Objects (EROs), selected both to have color redder than R-[3.6] = 4.0, and to have confirmed spectroscopic redshifts. Together, these two constraints result in this sample populating a fairly narrow redshift range at 0.76 < z < 1.42. The key new ingredient included here is deep Spitzer Space Telescope InfraRed Array Camera (IRAC) data. Based on [3.6]-[8.0] color, we demonstrate that it is possible to classify EROs into early-type, dusty starburst, or power-law (AGN) types. We present ultraviolet to mid-infrared spectral energy distributions (SEDs) and Advanced Camera for Surveys (ACS) images, both of which support our simple IRAC color classification.

A new 6-cm survey of almost 0.6 square degrees to a limit of 0.55-mJy/beam (10-sigma) finds 37 isolated radio sources and 7 radio source pairs (not necessarily physical companions). IRAC counterparts are identified for at least 92% of the radio sources within the area of deep IRAC coverage, which includes 31 isolated sources and 6 pairs. This contrasts with an identification rate of <74% to R<23.95 in visible light. Eight of the IRAC galaxies have power law spectral energy distributions, implying that the mid-infrared emission comes from a powerful AGN. The remaining 26 IRAC galaxies show stellar emission in the mid-infrared, probably in most of these galaxies because the stellar emission is bright enough to outshine an underlying AGN. The infrared colors suggest that the majority of these galaxies are bulge-dominated and have redshifts between approximately 0.5 and 1.0. Visible spectra from the DEEP2 redshift survey, available for 11 galaxies, are consistent with this suggestion. The IRAC galaxies fall into two distinct groups in a color-magnitude diagram, one group (the "stripe") includes all the AGN. The other group (the "blue clump") has blue 3.6 to 8 micron colors and a small range of 8 micron magnitudes. This separation should be useful in classifying galaxies found in other radio surveys.

We identify 17 possible 8.0 micron-selected counterparts to the submillimeter galaxies in the CUDSS 14-hour field, derived from deep imaging carried out with the IRAC and MIPS instruments aboard the Spitzer Space Telescope. Ten of the 17 counterparts are not the same as those previously identified at shorter wavelengths. We argue that 8.0 micron selection offers a better means for identifying counterparts to submillimeter galaxies than near-infrared or optical selection. Based on the panchromatic SEDs, most counterparts appear to be powered by ongoing star formation. Power-law fits to the SEDs suggest that five objects in the 8.0 micron-selected counterpart sample harbor dominant AGNs; a sixth object is identified as a possible AGN. The 3.6 to 8.0 micron colors of the infrared-selected counterparts are significantly redder than the general IRAC galaxy population in the CUDSS 14-hour field.

Mid-infrared observations of Active Galactic Nuclei (AGN) are important for understanding of the physical conditions around the central accretion engines. Chandra and XMM-Newton X-ray observations of a 300 arcmin^2 region in the Extended Groth Strip are used to select a sample of~150 AGN. The Spitzer instruments IRAC and MIPS detect 68-80% of these sources, which show a wide range of mid-infrared properties. About 40% of the sources have red power-law spectral energy distributions (f_nu ~ nu^alpha, alpha<0) in the 3.6-8 um IRAC bands. In these sources the central engine dominates the emission at both X-ray and IR wavelengths. Another 40% of the sources have blue mid-IR spectral energy distributions (alpha>0) with their infrared emission dominated by the host galaxy; the remaining 20% are not well-fit by a power law. Published IRAC color criteria for AGN select most of the red sources, but only some of the blue sources. As with all other known methods, selecting AGN with mid-IR colors will not produce a sample that is simultaneously complete and reliable. The IRAC SED type does not directly correspond to X-ray spectral type (hard/soft). The mid-IR properties of X-ray-detected Lyman-break, radio, submillimeter, and optically-faint sources vary widely and, for the most part, are not distinct from those of the general X-ray/infrared source population. X-ray sources emit 6-11% of the integrated mid--IR light, making them significant contributors to the cosmic infrared background.

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 Spitzer IRAC (3.6, 4.5, 5.8, 8.0 um) and MIPS (24 um) observations of the SCUBA submillimeter sources and uJy VLA radio sources in a 5'x5' area in the Lockman Hole East region. Out of the ~40 SCUBA/VLA sources in the field, Spitzer counterparts were detected for nearly all except for the few low-weight SCUBA detections. We show that the majority (80-90%) of the detected sources are cold (i.e., starburst-like) infrared-luminous galaxies (L_{IR} > 10^{11} Lsun) at redshift 0.5<z<3.5, whose star-formation rate density (SFRD) is comparable to that of the optically-selected star-forming galaxies.

We investigate the luminous X-ray sources in the Lockman Hole (LH) and the Extended Groth Strip (EGS) detected at 24microns using MIPS and also with IRAC on board Spitzer. We assemble optical/infrared spectral energy distributions (SEDs) for 45 X-ray/24micron sources in the EGS and LH. Only about 1/4 of the hard X-ray/24micron sources show pure type 1 AGN SEDs. More than half of the X-ray/24micron sources have stellar-emission-dominated or obscured SEDs, similar to those of local type 2 AGN and spiral/starburst galaxies. One-third of the sources detected in hard X-rays do not have a 24micron counterpart. Two such sources in the LH have SEDs resembling those of S0/elliptical galaxies. The broad variety of SEDs in the optical-to-Spitzer bands of X-ray selected AGN means that AGN selected according to the behavior in the optical/infrared will have to be supplemented by other kinds of data (e.g., X-ray) to produce unbiased samples of AGN.

We present Spitzer observations in five wavebands between 3.6 and 24um of an unbiased sample of 9 luminous, dusty galaxies selected at 1200um by the MAMBO camera on the IRAM 30-m telescope, a population akin to the well-known submm or `SCUBA' galaxies (hereafter SMGs). Owing to the coarse resolution of submm/mm instrumentation, SMGs have traditionally been difficult to identify at other wavelengths. We compare our multi-wavelength catalogs to show that the overlap between 24 and 1200um must be close to complete at these flux levels. We find that all (4/4) of the most secure >=4sigma SMGs have robust >=4sigma counterparts at 1.4GHz, while the fraction drops to 7/9 using all >=3sigma SMGs. We show that combining mid-IR and marginal (>=3sigma) radio detections provides plausible identifications in the remaining cases, enabling us to identify the complete sample. Accretion onto an obscured central engine is betrayed by the shape of the mid-IR continuum emission for several sources, confirming Spitzer's potential to weed out active galaxies. We demonstrate the power of a S(24um)/S(8um) vs S(8um)/S(4.5um) color-color plot as a diagnostic for this purpose. However, we conclude that the majority (~75%) of SMGs have rest-frame mid-/far-IR SEDs commensurate with obscured starbursts. Sensitive 24-um observations are clearly a useful route to identify and characterize reliable counterparts to high-redshift far-IR-bright galaxies, complementing what is possible via deep radio imaging.

We present preliminary results on 24micron detections of luminous infrared galaxies at z>1 with the Multiband Imaging Photometer for Spitzer (MIPS). Observations were performed in the Lockman Hole and the Extended Groth Strip (EGS), and were supplemented by data obtained with the Infrared Array Camera (IRAC) between 3 and 9microns. The positional accuracy of ~2arcsec for most MIPS/IRAC detections provides unambiguous identifications of their optical counterparts. Using spectroscopic redshifts from the Deep Extragalactic Evolutionary Probe survey, we identify 24micron sources at z>1 in the EGS, while the combination of the MIPS/IRAC observations with $BVRIJHK$ ancillary data in the Lockman Hole also shows very clear cases of galaxies with photometric redshifts at 1<z<2.5. The observed 24micron fluxes indicate infrared luminosities greater than 10^11 L_sol, while the data at shorter wavelengths reveal rather red and probably massive (M>=M*) galaxy counterparts. It is the first time that this population of luminous objects is detected up to z~2.5 in the infrared. Our work demonstrates the ability of the MIPS instrument to probe the dusty Universe at very high redshift, and illustrates how the forthcoming Spitzer deep surveys will offer a unique opportunity to illuminate a dark side of cosmic history not explored by previous infrared experiments.

We present sub-millimeter statistical detections of galaxies discovered in the 5'x5' Spitzer Early Release Observations (to 4-15 microJy 5 sigma at 3.6-8 microns, 170 microJy at 24 microns) through a stacking analysis of our reanalysed SCUBA 8mJy survey maps, and a Spitzer identification of a new sub-millimeter point source in the 8mJy survey region. For sources detected at 5.8 or 8 microns (154 and 111 sources respectively), we detect positive skews in the sub-millimeter flux distributions at 99.2-99.8% confidence using Kolmogorov-Smirnov tests, at both 850 microns and 450 microns. We also marginally detect the Spitzer 24 micron galaxies at 850 microns at 97% confidence, and place limits on the mean sub-millimeter fluxes of the 3.6 and 4.5 micron sources. Integrating the sub-millimeter fluxes of the Spitzer populations, we find the 5.8 micron galaxies contribute 0.12 +/- 0.05 nW/m^2/sr to the 850 micron background, and 2.4 +/- 0.7 nW/m^2/sr to the 450 micron background; similar contributions are made by the 8 micron-selected sample. We infer that the populations dominating the 5.8 and 8 micron extragalactic background light also contribute around a quarter of the 850 micron background and the majority of the 450 micron background.

We investigate Extremely Red Objects (EROs) using near- and mid-infrared observations in five passbands (3.6 to 24 micron) obtained from the Spitzer Space Telescope, and deep ground-based R and K imaging. The great sensitivity of the IRAC camera allows us to detect 64 EROs in only 12 minutes of IRAC exposure time, by means of an R-[3.6] color cut (analogous to the traditional red R-K cut). A pure infrared K-[3.6] red cut detects a somewhat different population and may be more effective at selecting z > 1.3 EROs. We find 17% of all galaxies detected by IRAC at 3.6 or 4.5 micron to be EROs. These percentages rise to about 40% at 5.8 micron, and about 60% at 8.0 micron. We utilize the spectral bump at 1.6 micron to divide the EROs into broad redshift slices using only near-infrared colors (2.2/3.6/4.5 micron). We conclude that two-thirds of all EROs lie at redshift z > 1.3. Detections at 24 micron imply that at least 11% of 0.6 < z < 1.3 EROs and at least 22% of z > 1.3 EROs are dusty star-forming 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 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.

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.

We present further observations of the extremely red object ERO J094258+4659.2, identified by \citet{iye00} as ERO R1 in their deep images of the cluster A851. We estimate its redshift independently by eight-band photometric redshift determination and cross-correlation of a new H-band spectrum with the optical spectra of local E/S0 galaxies, and conclude that it lies at $z \sim 1.5$. Although its colors are consistent both with an elliptical galaxy and an S0 galaxy at that redshift, its elongated shape and exponential luminosity profile suggest the presence of an evolved stellar disk component. We rule out the possibility that these properties are strongly influenced by gravitational lensing by the foreground cluster, and therefore conclude that this object is more likely to be an S0-like galaxy, rather than a lensed elliptical. The H-band spectrum does not show strong H$\alpha$ emission and the star formation rate therefore appears to be very modest. The presence of such a galaxy with an apparently relaxed disk of stars at this high redshift provides a new and strong constraint on theoretical models which aim to explain the formation and evolution of galaxies.

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

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

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.

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.

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.

Deep near-infrared images of a blank 2'x2' section of sky near the Galactic north pole taken by Subaru Telescope are presented. The total integration times of the J and K' bands are 12.1 hours and 9.7 hours, resulting in 5-sigma limiting magnitudes of 25.1 and 23.5 mag, respectively. The numbers of sources within these limiting magnitudes found with an automated detection procedure are 385 in the J band and 350 in K'. Based on photometric measurements of these sources, we present number count vs. magnitude relations, color vs. magnitude diagrams, size vs. color relationships, etc. The slope of the galaxy number count plotted against the AB magnitude scale is about 0.23 in the 22 to 26 AB magnitude range of both bands. The spatial number density of galaxies as well as the slopes in the faint-end region given by the Subaru Deep Field (SDF) survey is consistent with those given by HST-NICMOS surveys as expressed on the AB magnitude diagram. Several sources having very large J-K' color are found including a few K' objects without detection at J. In addition, a number of faint Galactic stars are also detected, most of which are assigned to M-subdwarfs, together with a few brown dwarf candidates.

A study of the stellar population of the M31 spiral arm around OB association A24 was carried out based on the photometric data obtained from deep V and JHK imaging. The luminosity function was obtained for -7 <~ Mbol <~ -3.5 by applying the extinction correction corresponding to Av=1 and the bolometric correction BC(K) as an empirical function of (J-K)o. In comparing the observed color-luminosity diagrams with semitheoretical isochrones modified for the dust-shell effects, we found the young population of t <~ 30 Myr with supergiants of Mbol <~ -5, the bulk of the intermediate-age population of t ~ 0.2 - 2.5 Gyr with bright asymptotic giant branch (AGB) stars of -5 <~ Mbol <~ -4, and old populations of t ~> 3 Gyr with AGB and red giant branch (RGB) stars of Mbol ~> -4. The average star formation rate was estimated to be ~1.8x10^4 M_o/Myr and ~0.7x10^4 M_o/Myr per deprojected disk area of 1 kpc^2 from the number density of B0 V stars around Mv=-4.0 (age ~10 Myr) and the number density of bright AGB stars around Mbol = -4.3 (age ~1 Gyr), respectively. A study of the local variation in the V and the J and H luminosity functions revealed a kind of anticorrelation between the population of the young component and that of the intermediate-age component when subdomains of ~100 pc scales were concerned. This finding suggests that the disk domain around the A24 area experienced a series of star formation episodes alternatively among different subdomains with a timescale of a few spiral passage periods. Brief discussions are given about the interstellar extinction and about the lifetimes of bright AGB stars and the highly red objects (HROs) in the same area.