• 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 present the structure of the Milky Way stellar halo beyond Galactocentric distances of $r = 50$ kpc traced by blue horizontal-branch (BHB) stars, which are extracted from the survey data in the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). We select BHB candidates based on $(g,r,i,z)$ photometry, where the $z$-band is on the Paschen series and the colors that involve the $z$-band are sensitive to surface gravity. About 450 BHB candidates are identified between $r = 50$ kpc and 300 kpc, most of which are beyond the reach of previous large surveys including the Sloan Digital Sky Survey. We find that the global structure of the stellar halo in this range has substructures, which are especially remarkable in the GAMA15H and XMM-LSS fields in the HSC-SSP. We find that the stellar halo can be fitted to a single power-law density profile with an index of $\alpha \simeq 3.3$ ($3.5$) with (without) these fields and its global axial ratio is $q \simeq 2.2$ ($1.3$). Thus, the stellar halo may be significantly disturbed and be made in a prolate form by halo substructures, perhaps associated with the Sagittarius stream in its extension beyond $r \sim 100$ kpc. For a broken power-law model allowing different power-law indices inside/outside a break radius, we obtain a steep power-law slope of $\alpha \sim 5$ outside a break radius of $\simeq 100$ kpc ($200$ kpc) for the case with (without) GAMA15H and XMM-LSS. This radius of $200$ kpc might be as close as a halo boundary if there is any, although larger BHB sample is required from further HSC-SSP survey to increase its statistical significance.
  • We report the result of optical identifications of FIRST radio sources with the Hyper Suprime-Cam Subaru Strategic Program survey (HSC-SSP). The positional cross-match within 1" between the FIRST and HSC-SSP catalogs (i ~< 26) produced more than 3600 optical counterparts in the 156 deg^2 of the HSC-SSP field. The matched counterparts account for more than 50 % of the FIRST sources in the search field, which substantially exceeds previously reported fractions of SDSS counterparts (i ~< 22) of ~ 30 %. Among the matched sample, 9 % are optically unresolved sources such as radio-loud quasars. The optically faint (i > 21) radio galaxies (RGs) show a flatter slope of a fitting linear function of the 1.4 GHz source counts than the bright RGs, while optically faint radio quasars show a steeper slope than the bright radio quasars. The optically faint RGs shows a flat slope in the i-band number counts down to 24 mag, implying either less-massive or distant radio-AGNs beyond 24 mag. The photometric redshift and the comparison of colors with the galaxy models show that most of the matched RGs are distributed at redshifts from 0 to 1.5. The optically faint sample includes the high radio-loudness sources that are not seen in the optically bright sample. Such sources are located at a redshift of larger than 1. This study provides a large number of radio-AGNs lying at the optically faint end and high redshift regime that are not probed by the previous searches.
  • March 19, 2018 cs.CV, cs.LG
    The network-based machine learning algorithm is very powerful tools. However, it requires huge training dataset. Researchers often meet privacy issues when they collect image dataset especially for surveillance applications. A learnable image encryption scheme is introduced. The key idea of this scheme is to encrypt images, so that human cannot understand images but the network can be train with encrypted images. This scheme allows us to train the network without the privacy issues. In this paper, a simple learnable image encryption algorithm is proposed. Then, the proposed algorithm is validated with cifar dataset.
  • We present a multi-wavelength catalog in the Subaru-XMM Deep Field (SXDF) as part of the Spitzer Large Area Survey with Hyper-Suprime-Cam (SPLASH). We include the newly acquired optical data from the Hyper-Suprime Cam Subaru Strategic Program, accompanied by IRAC coverage from the SPLASH survey. All available optical and near-infrared data is homogenized and resampled on a common astrometric reference frame. Source detection is done using a multi-wavelength detection image including the $u$-band to recover the bluest objects. We measure multi-wavelength photometry and compute photometric redshifts as well as stellar masses for $\sim$1.17 million objects over $\sim$4.2 deg$^2$ with $\sim$800,000 objects in the 2.4 deg$^2$ HSC-UltraDeep coverage. Using the available spectroscopic redshifts from various surveys over the range of $0<z<6$, we verify the performance of the photometric redshifts and we find a normalized median absolute deviation of 0.023 and outlier fraction of 3.2%. The SPLASH-SXDF catalog is a valuable, publicly available resource that is perfectly suited for studying galaxies in the early universe and tracing their evolution through cosmic time. The catalog is available for download from https://z.umn.edu/SXDF .
  • 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 report our first discoveries of high-redshift supernovae from the Subaru HIgh-Z sUpernova CAmpaign (SHIZUCA), a transient survey using Subaru/Hyper Suprime-Cam. We report the discovery of three supernovae at spectroscopically-confirmed redshifts of 2.399 (HSC16adga), 1.965 (HSC17auzg), and 1.851 (HSC17dbpf), and two supernova candidates with host-galaxy photometric redshifts of 3.2 (HSC16apuo) and 4.2 (HSC17dsid), respectively. In this paper, we present their photometric properties and the spectroscopic properties of the confirmed high-redshift supernovae are presented in the accompanying paper Curtin et al. (2018). The supernovae with the confirmed redshifts of z ~ 2 have rest ultraviolet peak magnitudes of around -21 mag, which make them superluminous supernovae. The discovery of three supernovae at z ~ 2 roughly corresponds to an event rate of ~ 900 Gpc-3 yr-1, which is already consistent with the total superluminous supernova rate estimated by extrapolating the local rate based on the cosmic star-formation history. Adding unconfirmed superluminous supernova candidates would increase the event rate. Our superluminous supernova candidates at the redshifts of around 3 and 4 indicate minimum superluminous supernova rates of ~ 400 Gpc-3 yr-1 (z ~ 3) and ~ 500 Gpc-3 yr-1 (z ~ 4). Because we have only performed a pilot search for high-redshift supernovae so far and have not completed selecting all the high-redshift supernova candidates, these rates are lower limits. Our initial results demonstrate the amazing capability of Hyper Suprime-Cam to discover high-redshift supernovae.
  • We present Keck spectroscopic confirmation of three superluminous supernovae (SLSNe) at z = 1.851, 1.965 and 2.399 detected as part of the Subaru HIgh-Z sUpernova CAmpaign (SHIZUCA). The host galaxies have multi-band photometric redshifts consistent with the spectroscopic values. The supernovae were detected during their rise, allowing the spectra to be taken near maximum light. The restframe far-ultraviolet (FUV; ~1000--2500A) spectra are made up in flux of approximately equal parts supernova and host galaxy. Weather conditions during observations were not ideal, and while the signal-to-noise ratios of the spectra are sufficient for redshift confirmation, the type of each event remains ambiguous. We compare our spectra to the few low redshift SLSN FUV spectra available to date and offer an interpretation as to the type of each supernova. We prefer SLSN-II classifications for all three events. The success of the first SHIZUCA Keck spectroscopic follow-up program is encouraging. It demonstrates that campaigns such as SHIZUCA are capable of identifying high redshift SLSNe with sufficient accuracy, speed and depth for rapid, well-cadenced and informative follow-up.
  • We report fourteen and twenty-eight protocluster candidates at z=5.7 and 6.6 over 14 and 19 deg^2 areas, respectively, selected from 2,230 (259) Lya emitters (LAEs) photometrically (spectroscopically) identified with Subaru/Hyper Suprime-Cam (HSC) deep images (Keck, Subaru, and Magellan spectra and the literature data). Six out of the 42 protocluster candidates include 1-12 spectroscopically confirmed LAEs at redshifts up to z=6.574. By the comparisons with the cosmological Lya radiative transfer (RT) model reproducing LAEs with the reionization effects, we find that more than a half of these protocluster candidates are progenitors of the present-day clusters with a mass of > 10^14 M_sun. We then investigate the correlation between LAE overdensity delta and Lya rest-frame equivalent width EW_Lya^rest, because the cosmological Lya RT model suggests that a slope of EW_Lya^rest-delta relation is steepened towards the epoch of cosmic reionization (EoR), due to the existence of the ionized bubbles around galaxy overdensities easing the escape of Lya emission from the partly neutral intergalactic medium (IGM). The available HSC data suggest that the slope of the EW_Lya^rest-delta correlation does not evolve from the post-reionization epoch z=5.7 to the EoR z=6.6 beyond the moderately large statistical errors. There is a possibility that we would detect the evolution of the EW_Lya^rest - delta relation from z=5.7 to 7.3 by the upcoming HSC observations providing large samples of LAEs at z=6.6-7.3.
  • We present optimized source galaxy selection schemes for measuring cluster weak lensing (WL) mass profiles unaffected by cluster member dilution from the Subaru Hyper Suprime-Cam Strategic Survey Program (HSC-SSP). The ongoing HSC-SSP survey will uncover thousands of galaxy clusters to $z\lesssim1.5$. In deriving cluster masses via WL, a critical source of systematics is contamination and dilution of the lensing signal by cluster {members, and by foreground galaxies whose photometric redshifts are biased}. Using the first-year CAMIRA catalog of $\sim$900 clusters with richness larger than 20 found in $\sim$140 deg$^2$ of HSC-SSP data, we devise and compare several source selection methods, including selection in color-color space (CC-cut), and selection of robust photometric redshifts by applying constraints on their cumulative probability distribution function (PDF; P-cut). We examine the dependence of the contamination on the chosen limits adopted for each method. Using the proper limits, these methods give mass profiles with minimal dilution in agreement with one another. We find that not adopting either the CC-cut or P-cut methods results in an underestimation of the total cluster mass ($13\pm4\%$) and the concentration of the profile ($24\pm11\%$). The level of cluster contamination can reach as high as $\sim10\%$ at $R\approx 0.24$ Mpc/$h$ for low-z clusters without cuts, while employing either the P-cut or CC-cut results in cluster contamination consistent with zero to within the 0.5% uncertainties. Our robust methods yield a $\sim60\sigma$ detection of the stacked CAMIRA surface mass density profile, with a mean mass of $M_\mathrm{200c} = (1.67\pm0.05({\rm {stat}}))\times 10^{14}\,M_\odot/h$.
  • We 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 performed a detailed study of the evolution of the star formation rate (SFR) and stellar mass of the brightest group galaxies (BGGs) and their relative contribution to the total baryon budget within $R_{200}$ ($f^{BGG}_{b,200}$). The sample comprises 407 BGGs selected from X-ray galaxy groups ($M_{200}=10^{12.8}-10^{14} \;M_{\odot}$) out to $z\sim1.3$ identified in the COSMOS, XMM-LSS, and AEGIS fields. We find that BGGs constitute two distinct populations of quiescent and star-forming galaxies and their mean SFR is $\sim2$ dex higher than the median SFR at $ z<1.3 $. Both the mean and the median SFRs decline with time by $>2$ dex. The mean (median) of stellar mass of BGGs has grown by $0.3$ dex since $z=1.3$ to the present day. We show that up to $\sim45\% $ of the stellar mass growth in a star-forming BGG can be due to its star-formation activity. With respect to $f^{BGG}_{b,200}$, we find it to increase with decreasing redshift by $\sim0.35$ dex while decreasing with halo mass in a redshift dependent manner. We show that the slope of the relation between $f^{BGG}_{b,200}$ and halo mass increases negatively with decreasing redshift. This trend is driven by an insufficient star-formation in BGGs, compared to the halo growth rate. We separately show the BGGs with the 20\% highest $f^{BGG}_{b,200}$ are generally non-star-forming galaxies and grow in mass by processes not related to star formation (e.g., dry mergers and tidal striping). We present the $ M_\star-M_h $ and $ M_\star/M_h-M_h $ relations and compare them with semi-analytic model predictions and a number of results from the literature. We quantify the intrinsic scatter in stellar mass of BGGs at fixed halo mass ($\sigma_{log M_{\star}}$) and find that $\sigma_{log M_{\star}}$ increases from 0.3 dex at $ z\sim0.2 $ to 0.5 dex at $ z\sim1.0 $ due to the bimodal distribution of stellar mass.
  • We present and characterize the catalog of galaxy shape measurements that will be used for cosmological weak lensing measurements in the Wide layer of the first year of the Hyper Suprime-Cam (HSC) survey. The catalog covers an area of 136.9 deg$^2$ split into six fields, with a mean $i$-band seeing of $0.58$ arcsec and $5\sigma$ point-source depth of $i\sim 26$. Given conservative galaxy selection criteria for first year science, the depth and excellent image quality results in unweighted and weighted source number densities of 24.6 and 21.8 arcmin$^{-2}$, respectively. We define the requirements for cosmological weak lensing science with this catalog, then focus on characterizing potential systematics in the catalog using a series of internal null tests for problems with point-spread function (PSF) modeling, shear estimation, and other aspects of the image processing. We find that the PSF models narrowly meet requirements for weak lensing science with this catalog, with fractional PSF model size residuals of approximately $0.003$ (requirement: 0.004) and the PSF model shape correlation function $\rho_1<3\times 10^{-7}$ (requirement: $4\times 10^{-7}$) at 0.5$^\circ$ scales. A variety of galaxy shape-related null tests are statistically consistent with zero, but star-galaxy shape correlations reveal additive systematics on $>1^\circ$ scales that are sufficiently large as to require mitigation in cosmic shear measurements. Finally, we discuss the dominant systematics and the planned algorithmic changes to reduce them in future data reductions.
  • We present wide and deep photometry of the northwest part of the halo of the Andromeda galaxy (M31) using Hyper Suprime-Cam on the Subaru Telescope. The survey covers 9.2 deg$^{2}$ field in the $g$, $i$, and $NB515$ bands and shows a clear red giant branch (RGB) of M31's halo stars and a pronounced red clump (RC) feature. The spatial distribution of RC stars shows a prominent stream feature, the North Western (NW) Stream, and a diffuse substructure in the south part of our survey field. We estimate the distances based on the RC method and obtain $(m-M)$ = 24.63$\pm 0.191$(random)$\pm0.057$(systematic) and 24.29$\pm 0.211$(random)$\pm0.057$(systematic) mag for the NW stream and diffuse substructure, respectively, implying that the NW Stream is located behind M31, whereas the diffuse substructure is located in front. We also estimate line-of-sight distances along the NW Stream and find that the south part of the stream is $\sim$20 kpc closer to us relative to the north part. The distance to the NW Stream inferred from the isochrone fitting to the color-magnitude diagram favors the RC-based distance, but the TRGB-based distance estimated for $NB515$-selected RGB stars does not agree with it. The surface number density distribution of RC stars across the NW Stream is found to be approximately Gaussian with a FWHM of $\sim$25 arcmin (5.7 kpc), with a slight skew to the south-west side. That along the NW Stream shows a complicated structure including variations in number density and a significant gap in the stream.
  • The 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 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.
  • 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 an unprecedentedly large catalog consisting of 2,230 > L* Lya emitters (LAEs) at z=5.7 and 6.6 on the 13.8 and 21.2 deg2 sky, respectively, that are identified by the SILVERRUSH program with the first narrowband imaging data of the Hyper Suprime-Cam (HSC) survey. We confirm that the LAE catalog is reliable on the basis of 96 LAEs whose spectroscopic redshifts are already determined by this program and the previous studies. This catalogue is also available on-line. Based on this catalogue, we derive the rest-frame Lya equivalent-width distributions of LAEs at z~5.7-6.6 that are reasonably explained by the exponential profiles with the scale lengths of ~120-170A, showing no significant evolution from z~5.7 to z~6.6. We find that 275 LAEs with a large equivalent width (LEW) of >240A are candidates of young-metal poor galaxies and AGNs. We also find that the fraction of LEW LAEs to all ones is 4% and 21% at z~5.7 and z~6.6, respectively. Our LAE catalog includes 11 Lya blobs (LABs) that are LAEs with spatially extended Lya emission whose profile is clearly distinguished from those of stellar objects at the >~ 3sigma level. The number density of the LABs at z=6-7 is ~ 10^-7-10^-6 Mpc^-3, being ~ 10-100 times lower than those claimed for LABs at z~ 2-3, suggestive of disappearing LABs at z>~6, albeit with the different selection methods and criteria for the low and high-z LABs.
  • We present Lya and UV-nebular emission line properties of bright Lya emitters (LAEs) at z=6-7 with a luminosity of log L_Lya/[erg s-1] = 43-44 identified in the 21-deg2 area of the SILVERRUSH early sample developed with the Subaru Hyper Suprime-Cam (HSC) survey data. Our optical spectroscopy newly confirm 21 bright LAEs with clear Lya emission, and contribute to make a spectroscopic sample of 96 LAEs at z=6-7 in SILVERRUSH. From the spectroscopic sample, we select 7 remarkable LAEs as bright as Himiko and CR7 objects, and perform deep Keck/MOSFIRE and Subaru/nuMOIRCS near-infrared spectroscopy reaching the 3sigma-flux limit of ~ 2x10^{-18} erg s-1 for the UV-nebular emission lines of He II1640, C IV1548,1550, and O III]1661,1666. Except for one tentative detection of C IV, we find no strong UV-nebular lines down to the flux limit, placing the upper limits of the rest-frame equivalent widths (EW_0) of ~2-4 A for He II, C IV, and O III] lines. Here we also investigate the VLT/X-SHOOTER spectrum of CR7 whose 6 sigma detection of He II is claimed by Sobral et al. Although two individuals and the ESO-archive service carefully re-analyze the X-SHOOTER data that are used in the study of Sobral et al., no He II signal of CR7 is detected, supportive of weak UV-nebular lines of the bright LAEs even for CR7. Spectral properties of these bright LAEs are thus clearly different from those of faint dropouts at z~7 that have strong UV-nebular lines shown in the various studies. Comparing these bright LAEs and the faint dropouts, we find anti-correlations between the UV-nebular line EW_0 and UV-continuum luminosity, which are similar to those found at z~2-3.
  • We 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 constrain the quasar contribution to cosmic reionization based on our deep optical survey of z~6 quasars down to z_R=24.15 using Subaru/Suprime-Cam in three UKIDSS-DXS fields covering 6.5 deg^2. In Kashikawa et al. (2015), we select 17 quasar candidates and report our initial discovery of two low-luminosity quasars (M_1450~ -23) from seven targets, one of which might be a Lyman alpha emitting galaxy. From an additional optical spectroscopy, none of the four candidates out of the remaining ten turn out to be genuine quasars. Moreover, the deeper optical photometry provided by the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) shows that, unlike the two already-known quasars, the i-z and z-y colors of the last six candidates are consistent with M- or L-type brown dwarfs. Therefore, the quasar luminosity function (QLF) in the previous paper is confirmed. Compiling QLF measurements from the literature over a wide magnitude range, including an extremely faint AGN candidate from Parsa et al. (2017}, to fit them with a double power-law, we find that the best-fit faint-end slope is alpha=-2.04^+0.33_-0.18 (-1.98^+0.48_-0.21) and characteristic magnitude is M_1450^*=-25.8^+1.1_-1.9 (-25.7^+1.0_-1.8) in the case of two (one) quasar detection. Our result suggests that, if the QLF is integrated down to M_1450=-18, quasars produce ~1-12% of the ionizing photons required to ionize the whole universe at z~6 with 2sigma confidence level, assuming that the escape fraction is f_esc=1 and the IGM clumpy factor is C=3. Even when the systematic uncertainties are taken into account, our result supports the scenario that quasars are the minor contributors of reionization.
  • 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 investigate the galaxy overdensity around proto-cluster scale quasar pairs at high (z>3) and low (z~1) redshift based on the unprecedentedly wide and deep optical survey of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). Using the first-year survey data covering effectively ~121 deg^2 with the 5sigma depth of i~26.4 and the SDSS DR12Q catalog, we find two luminous pairs at z~3.3 and 3.6 which reside in >5sigma overdense regions of g-dropout galaxies at i<25. The projected separations of the two pairs are R_perp=1.75 and 1.04 proper Mpc, and their velocity offsets are Delta V=692 and 1448 km s^{-1}, respectively. This result is in clear contrast to the average z~4 quasar environments as discussed in Uchiyama et al. (2017) and implies that the quasar activities of the pair members are triggered via major mergers in proto-clusters, unlike the vast majority of isolated quasars in general fields that may turn on via non-merger events such as bar and disk instabilities. At z~1, we find 37 pairs with R_perp<2 pMpc and Delta V<2300 km s^{-1} in the current HSC-Wide coverage, including four from Hennawi et al. (2006). The distribution of the peak overdensity significance within two arcminutes around the pairs has a long tail toward high density (>4sigma) regions. Thanks to the large sample size, we find a statistical evidence that this excess is unique to the pair environments when compared to single quasar and randomly selected galaxy environments at the same redshift range. Moreover, there are nine small-scale (R_perp<1 pMpc) pairs, two of which are found to reside in cluster fields. Our results demonstrate that <2 pMpc-scale quasar pairs at both redshift range tend to occur in massive haloes, although perhaps not the most massive ones, and that they are useful to search for rare density peaks.
  • We present initial results from the Subaru Strategic Program (SSP) with Hyper Suprime-Cam (HSC) on a comprehensive survey of emission-line galaxies at z<1.5 based on narrowband (NB) imaging. The first Public Data Release (PDR1) provides us with data from two NB filters, specifically NB816 and NB921 over 5.7 deg$^2$ and 16.2 deg$^2$ respectively. The $5 \sigma$ limiting magnitudes are 25.2 (UDeep layer, 1.4 deg$^2$) and 24.8 (Deep layer, 4.3 deg$^2$) mag in NB816, and 25.1 (UDeep, 2.9 deg$^2$) and 24.6--24.8 (Deep, 13.3 deg$^2$) mag in NB921. The wide-field imaging allows us to construct unprecedentedly large samples of 8,054 H$\alpha$ emitters at z ~ 0.25 and 0.40, 8,656 [OIII] emitters at z ~ 0.63 and 0.84, and 16,877 [OII] emitters at z ~ 1.19 and 1.47. We map the cosmic web on scales out to about 50 comoving Mpc that includes galaxy clusters, identified by red sequence galaxies, located at the intersection of filamentary structures of star-forming galaxies. The luminosity functions of emission-line galaxies are measured with precision and consistent with published studies. The wide field coverage of the data enables us to measure the luminosity functions up to brighter luminosities than previous studies. The comparison of the luminosity functions between the different HSC-SSP fields suggests that a survey volume of $>5\times10^5$ Mpc$^3$ is essential to overcome cosmic variance. Since the current data have not reached the full depth expected for the HSC-SSP, the color cut in i-NB816 or z-NB921 induces a bias towards star-forming galaxies with large equivalent widths, primarily seen in the stellar mass functions for the H$\alpha$ emitters at z ~ 0.25--0.40. Even so, the emission-line galaxies clearly cover a wide range of luminosity, stellar mass, and environment, thus demonstrating the usefulness of the NB data from the HSC-SSP to investigate star-forming galaxies at z<1.5.