We present the physical properties of AKARI sources without optical counterparts in optical images from the Hyper Suprime-Cam (HSC) on the Subaru telescope. Using the AKARI infrared (IR) source catalog and HSC optical catalog, we select 583 objects that do not have HSC counterparts in the AKARI North Ecliptic Pole (NEP) wide survey field ($\sim 5$ deg$^{2}$). Because the HSC limiting magnitude is deep ($g_{\rm AB}$ $\sim 28.6$), these are good candidates for extremely red star-forming galaxies (SFGs) and/or active galactic nuclei (AGNs), possibly at high redshifts. We compile multi-wavelength data out to 500 $\mu$m and use it for Spectral Energy Distribution (SED) fitting with CIGALE to investigate the physical properties of AKARI galaxies without optical counterparts. We also compare their physical quantities with AKARI mid-IR selected galaxies with HSC counterparts. The estimated redshifts of AKARI objects without HSC counterparts range up to $z\sim 4$, significantly higher than that of AKARI objects with HSC counterparts. We find that: (i) 3.6 $-$ 4.5 $\mu$m color, (ii) AGN luminosity, (iii) stellar mass, (iv) star formation rate, and (v) $V$-band dust attenuation in the interstellar medium of AKARI objects without HSC counterparts are systematically larger than those of AKARI objects with counterparts. These results suggest that our sample includes luminous, heavily dust-obscured SFGs/AGNs at $z\sim 1-4$ that are missed by previous optical surveys, providing very interesting targets for the coming James Webb Space Telescope era.

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 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 use Hyper Suprime-Cam on the Subaru Telescope to investigate the structural and photometric properties of early-type dwarf galaxies and young stellar systems at the center of the M81 Group. We have mapped resolved stars to $\sim2$ magnitudes below the tip of the red giant branch over almost 6.5 square degrees, corresponding to a projected area of $160\times160 \rm{kpc}$ at the distance of M81. The resulting stellar catalogue enables a homogeneous analysis of the member galaxies with unprecedented sensitivity to low surface brightness emission. The radial profiles of the dwarf galaxies are well-described by Sersic and King profiles, and show no obvious signatures of tidal disruption. The measured radii for most of these systems are larger than the existing literature values and we find the total luminosity of IKN ($\rm{M_{V,0}}=-14.29$) to be almost 3 magnitudes brighter than previously-thought. We identify new dwarf satellite candidates, d1006+69 and d1009+68, which we estimate to lie at a distance of $4.3\pm0.2$ Mpc and $3.5\pm0.5$ Mpc. With $\rm{M_{V,0}}=-8.91\pm0.40$ and $\rm{[M/H]}=-1.83\pm0.28$, d1006+69 is one of the faintest and most metal-poor dwarf satellites currently-known in the M81 Group. The luminosity functions of young stellar systems in the outlying tidal HI debris imply continuous star formation in the recent past and the existence of populations as young as 30 Myr old. We find no evidence for old RGB stars coincident with the young MS/cHeB stars which define these objects, supporting the idea that they are genuinely new stellar systems resulting from triggered star formation in gaseous tidal debris.

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.

Much of the star formation is obscured by dust. For the complete understanding of the cosmic star formation history (CSFH), infrared (IR) census is indispensable. AKARI carried out deep mid-infrared observations using its continuous 9-band filters in the North Ecliptic Pole (NEP) field (5.4 deg$^2$). This took significant amount of satellite's lifetime, $\sim$10\% of the entire pointed observations. By combining archival Spitzer (5 bands) and WISE (4 bands) mid-IR photometry, we have, in total, 18 band mid-IR photometry, which is the most comprehensive photometric coverage in mid-IR for thousands of galaxies. However previously, we only had shallow optical imaging ($\sim$25.9ABmag) in a small area of 1.0 deg$^2$. As a result, there remained thousands of AKARI's infrared sources undetected in optical. Using the new Hyper Suprime-Cam on Subaru telescope, we obtained deep enough optical images of the entire AKARI NEP field in 5 broad bands ($g\sim$27.5mag). These provided photometric redshift, and thereby IR luminosity for the previously undetected faint AKARI IR sources. Combined with the accurate mid-IR luminosity measurement, we constructed mid-IR LFs, and thereby performed a census of dust-obscured CSFH in the entire AKARI NEP field. We have measured restframe 8$\mu$m, 12$\mu$m luminosity functions (LFs), and estimated total infrared LFs at 0.35$<$z$<$2.2. Our results are consistent with our previous work, but with much reduced statistical errors thanks to the large area coverage of the new data. We have possibly witnessed the turnover of CSFH at $z\sim$2.

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

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

We explore the relations between size, stellar mass and average stellar population age (indicated by D$_n4000$ indices) for a sample of $\sim11000$ intermediate-redshift galaxies from the SHELS spectroscopic survey (Geller et al. 2014) augmented by high-resolution Subaru Telescope Hyper Suprime-Cam imaging. In the redshift interval $0.1<z<0.6$, star forming galaxies are on average larger than their quiescent counterparts. The mass-complete sample of $\sim3500$ $M_*>10^{10}\, M_\odot$ quiescent galaxies shows that the average size of a $10^{11}\, M_\odot$ quiescent galaxy increases by $\lesssim25\%$ from $z\sim0.6$ to $z\sim0.1$. This growth rate is a function of stellar mass: the most massive ($M_*>10^{11}\, M_\odot$) galaxies grow significantly more slowly in size than an order of magnitude less massive quiescent systems that grow by 70\% in the $0.1\lesssim z\lesssim0.3$ redshift interval. For $M_*<10^{11}\, M_\odot$ galaxies age and size are anti-correlated at fixed mass; more massive quiescent systems show no significant trend in size with average stellar population age. The evolution in absolute and fractional abundances of quiescent systems at intermediate redshift are also a function of galaxy stellar mass. The suite of evolutionary trends suggests that galaxies more massive than $\sim10^{11}\, M_\odot$ have mostly assembled their mass by $z\sim0.6$. Quiescent galaxies with lower stellar masses show more complex evolution that is characterized by a combination of individual quiescent galaxy size growth (through mergers) and an increase in the size of newly quenched galaxies joining the population at later times (progenitor bias). The $M_*\sim10^{10}\, M_\odot$ population grows predominantly as a result of progenitor bias. For $M_*\sim5\times10^{10}\, M_\odot$ quiescent galaxies, mergers and progenitor bias make more comparable contributions to the size growth.[abridged]

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

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

We present 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 use the dense HectoMAP redshift survey to explore the properties of 104 redMaPPer cluster candidates. The redMaPPer systems in HectoMAP cover the full range of richness and redshift (0.08 $< z <$ 0.60). Fifteen systems included in the Subaru/Hyper Suprime-Cam public data release are bona fide clusters. The median number of spectroscopic members per cluster is $\sim20$. We include redshifts of 3547 member candidates listed in the redMaPPer catalog whether they are cluster members or not. We evaluate the redMaPPer membership probability spectroscopically. The scaled richness ({\lambda}rich/S) provided by redMaPPer correlates tightly with the spectroscopic richness regardless of the cluster redshift and appears to be a better mass proxy than the original richness, {\lambda}rich. The purity (number of real systems) in redMaPPer exceeds 90% even at the lowest richness; however, there is some incompleteness. Five massive galaxy clusters (M $\gtrsim 2 \times 10^{13}$ M$_{\odot}$) associated with X-ray emission in the HectoMAP region are missing from the catalog.

We apply a friends-of-friends algorithm to the HectoMAP redshift survey and cross-identify associated X-ray emission in the ROSAT All-Sky Survey data (RASS). The resulting flux limited catalog of X-ray cluster survey is complete to a limiting flux of $\sim3 \times10^{-13}$ erg s$^{-1}$ cm$^{-2}$ and includes 15 clusters (7 newly discovered) with redshift $z \leq 0.4$. HectoMAP is a dense survey ($\sim1200$ galaxies deg$^{-2}$) that provides $\sim50$ members (median) in each X-ray cluster. We provide redshifts for the 1036 cluster members. Subaru/Hyper Suprime-Cam imaging covers three of the X-ray systems and confirms that they are impressive clusters. The HectoMAP X-ray clusters have an $L_{X} - {\sigma}_{cl}$ scaling relation similar to that of known massive X-ray clusters. The HectoMAP X-ray cluster sample predicts $\sim 12000 \pm3000$ detectable X-ray clusters in the RASS to the limiting flux, comparable with previous estimates.

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.

Understanding infrared (IR) luminosity is fundamental to understanding the cosmic star formation history and AGN evolution. Japanese infrared satellite, AKARI, provided unique data sets to probe this both at low and high redshift; the AKARI all sky survey in 6 bands (9-160 $\mu$m), and the AKARI NEP survey in 9 bands (2-24$\mu$m). The AKARI performed all sky survey in 6 IR bands (9, 18, 65, 90, 140, and 160 $\mu$m) with 3-10 times better sensitivity than IRAS, covering the crucial far-IR wavelengths across the peak of the dust emission. Combined with a better spatial resolution, we measure the total infrared luminosity ($L_{TIR}$) of individual galaxies, and thus, the total infrared luminosity density of the local Universe much more precisely than previous work. In the AKARI NEP wide field, AKARI has obtained deep images in the mid-infrared (IR), covering 5.4 deg$^2$. However, our previous work was limited to the central area of 0.25 deg$^2$ due to the lack of deep optical coverage. To rectify the situation, we used the newly advent Subaru telescope's Hyper Suprime-Cam to obtain deep optical images over the entire 5.4 deg$^2$ of the AKARI NEP wide field. With this deep and wide optical data, we, for the first time, can use the entire AKARI NEP wide data to construct restframe 8$\mu$m, 12$\mu$m, and total infrared (TIR) luminosity functions (LFs) at 0.15$<z<$2.2. A continuous 9-band filter coverage in the mid-IR wavelength (2.4, 3.2, 4.1, 7, 9, 11, 15, 18, and 24$\mu$m) by the AKARI satellite allowed us to estimate restframe 8$\mu$m and 12$\mu$m luminosities without using a large extrapolation based on a SED fit, which was the largest uncertainty in previous work. By combining these two results, we reveal dust-hidden cosmic star formation history and AGN evolution from z=0 to z=2.2, all probed by the AKARI satellite.

The AKARI North Ecliptic Pole (NEP) survey consists of two survey projects: NEP-Deep (0.5 sq.deg) and NEP-Wide (5.4 sq.deg), providing with tens of thousands of galaxies. A continuous filter coverage in the mid-infrared wavelengths (7, 9, 11, 15, 18 and 24 $\mu$m) is unique to diagnose the contributions from dusty star-formation activity and AGNs. Here we present current status focused on the newly obtained optical images and near-future prospects with a new X-ray telescope. Hyper Suprime-Cam on Subaru telescope is a gigantic optical camera with huge Field of View (FoV). Thanks to the wide FoV, we successfully obtained deep optical images at g, r, i, z and Y-bands covering most of the NEP-Wide field. Using the deep optical images, we identified over 5000 optical counterparts of the mid-IR sources, presumably deeply obscured galaxies in NEP-Wide field. We also investigated properties of these infrared sources with SED-fitting. eROSITA, to be launched early 2018, is a new all-sky X-ray survey telescope, and expected to conduct ultra deep 2-10 keV imaging toward NEP. We expect unprecedentedly numerous Compton-thick AGN candidates when combined with the multi-wavelength data in NEP region.

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.

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.

We present the results of the detailed analysis of an optical imaging survey conducted using the Subaru / Hyper Suprime-Cam (HSC), which aims to identify an optical counterpart to the gravitational wave event GW151226. In half a night, the $i$- and $z$-band imaging survey by HSC covers 63.5deg$^2$ of the error region, which contains about 7\% of the LIGO localization probability, and the same field is observed in three different epochs. The detectable magnitude of the candidates in a differenced image is evaluated as $i \sim 23.2$ mag for the requirement of at least two 5$\sigma$ detections, and 1744 candidates are discovered. Assuming a kilonova as an optical counterpart, we compared the optical properties of the candidates with model predictions. A red and rapidly declining light curve condition enables the discrimination of a kilonova from other transients, and a small number of candidates satisfy this condition. The presence of stellar-like counterparts in the reference frame suggests that the surviving candidates are likely to be flare stars. The fact that most of those candidates are in galactic plane, $|b|<5^{\circ}$, supports this interpretation. We also checked whether the candidates are associated with the nearby GLADE galaxies, which reduces the number of contaminants even with a looser color cut. When a better probability map (with localization accuracy of $\sim50{\rm deg}^2$) is available, kilonova searches of up to approximately $200$ Mpc will become feasible by conducting immediate follow-up observations with an interval of 3--6 days.

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

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 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 report the results of optical--infrared follow-up observations of the gravitational wave (GW) event GW151226 detected by the Advanced LIGO in the framework of J-GEM (Japanese collaboration for Gravitational wave ElectroMagnetic follow-up). We performed wide-field optical imaging surveys with Kiso Wide Field Camera (KWFC), Hyper Suprime-Cam (HSC), and MOA-cam3. The KWFC survey started at 2.26 days after the GW event and covered 778 deg$^2$ centered at the high Galactic region of the skymap of GW151226. We started the HSC follow-up observations from 12 days after the event and covered an area of 63.5 deg$^2$ of the highest probability region of the northern sky with the limiting magnitudes of 24.6 and 23.8 for i band and z band, respectively. MOA-cam3 covered 145 deg$^2$ of the skymap with MOA-red filter 2.5 months after the GW alert. Total area covered by the wide-field surveys was 986.5 deg$^2$. The integrated detection probability of all the observed area was $\sim$29%. We also performed galaxy-targeted observations with six optical and near-infrared telescopes from 1.61 days after the event. Total of 238 nearby (<100 Mpc) galaxies were observed with the typical I band limiting magnitude of $\sim$19.5. We detected 13 supernova candidates with the KWFC survey, and 60 extragalactic transients with the HSC survey. Two third of the HSC transients were likely supernovae and the remaining one third were possible active galactic nuclei. With our observational campaign, we found no transients that are likely to be associated with GW151226.

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

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

We present our optical follow-up observations to search for an electromagnetic counterpart of the first gravitational wave source GW150914 in the framework of the Japanese collaboration for Gravitational wave ElectroMagnetic follow-up (J-GEM), which is an observing group utilizing optical and radio telescopes in Japan, as well as those in New Zealand, China, South Africa, Chile, and Hawaii. We carried out a wide-field imaging survey with Kiso Wide Field Camera (KWFC) on the 1.05-m Kiso Schmidt telescope in Japan and a galaxy-targeted survey with Tripole5 on the B&C 61-cm telescope in New Zealand. Approximately 24 deg2 regions in total were surveyed in i-band with KWFC and 18 nearby galaxies were observed with Tripole5 in g-, r-, and i-bands 4-12 days after the gravitational wave detection. Median 5-sigma depths are i~18.9 mag for the KWFC data and g~18.9 mag, r~18.7 mag, and i~18.3 mag for the Tripole5 data. Probability for a counterpart to be in the observed area is 1.2% in the initial skymap and 0.1% in the final skymap. We do not find any transient source associated to an external galaxy with spatial offset from its center, which is consistent with the local supernova rate. We summarize future prospects and ongoing efforts to pin down electromagnetic counterparts of binary black hole mergers as well as neutron star mergers.

We present the first results of a wide-field mapping survey of the M81 group conducted with Hyper Suprime-Cam on the Subaru Telescope. Our deep photometry reaches $\sim2$ magnitudes below the tip of the red giant branch (RGB) and reveals the spatial distribution of both old and young stars over an area of $\sim 100\times115$ kpc at the distance of M81. The young stars ($\sim30-160$ Myr old) closely follow the neutral hydrogen distribution and can be found in a stellar stream between M81 and NGC\,3077 and in numerous outlying stellar associations, including the known concentrations of Arp's Loop, Holmberg\,IX, an arc in the halo of M82, BK3N, and the Garland. Many of these groupings do not have counterparts in the RGB maps, suggesting they may be genuinely young systems. Our survey also reveals for the first time the very extended ($\geq 2\times \rm{R_{25}}$) halos of RGB stars around M81, M82 and NGC\,3077, as well as faint tidal streams that link these systems. The halos of M82 and NGC\,3077 exhibit highly disturbed morphologies, presumably a consequence of the recent gravitational encounter and their ongoing disruption. While the halos of M81, NGC\,3077 and the inner halo of M82 have the similar $(g-i)_{0}$ colors, the outer halo of M82 is significantly bluer indicating it is more metal-poor. Remarkably, our deep panoramic view of the M81 group demonstrates that the complexity long-known to be present in HI is equally matched in the low surface brightness stellar component.

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.

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.

The Leonids show meteor storms in a period of 33 years, and known as one of the most active meteor showers. It has recently shown a meteor stream consisting of several narrow dust trails made by meteoroids ejected from a parent comet. Hence, an analysis of the temporal behavior of the meteor flux is important to study the structure of the trails. However, statistical inference for the count data is not an easy task, because of its Poisson characteristics. We carried out a wide-field video observation of the Leonid meteor storm in 2001. We formulated a state-of-the-art statistical analysis, which is called a self-organizing state space model, to infer the true behavior of the dust density of the trails properly from the meteor count data. {}From this analysis, we found that the trails have a fairly smooth spatial structure, with small and dense clumps that cause a temporal burst of meteor flux. We also proved that the time behavior (trend) of the fluxes of bright meteors and that of faint meteors are significantly different. In addition we comment on some other application of the self-organizing state-space model in fields related to astronomy and astrophysics.

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.

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

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