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

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

We perform a stacking analysis of {\it Planck}, {\it AKARI}, Infrared Astronomical Satellite ($IRAS$), Wide-field Infrared Survey Eplorer ($WISE$), and {\it Herschel} images of the largest number of (candidate) protoclusters at $z\sim3.8$ selected from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). Stacking the images of the $179$ candidate protoclusters, the combined infrared (IR) emission of the protocluster galaxies in the observed $12-850~\mu$m wavelength range is successfully detected with $>5\sigma$ significance (at $Planck$). This is the first time that the average IR spectral energy distribution (SED) of a protocluster has been constrained at $z\sim4$. The observed IR SEDs of the protoclusters exhibit significant excess emission in the mid-IR compared to that expected from typical star-forming galaxies (SFGs). They are reproduced well using SED models of intense starburst galaxies with warm/hot dust heated by young stars, or by a population of active galactic nuclei (AGN)/SFG composites. For the pure star-forming model, a total IR (from 8 to 1000 $\mu$m) luminosity of $19.3_{-4.2}^{+0.6}\times10^{13}~L_{\odot}$ and a star formation rate (SFR) of $16.3_{-7.8}^{+1.0}\times10^3~M_{\odot}$ yr$^{-1}$ are found whereas for the AGN/SFG composite model, $5.1_{-2.5}^{+2.5}\times10^{13}~L_{\odot}$ and $2.1^{+6.3}_{-1.7}\times10^3~M_{\odot}$ yr$^{-1}$ are found. Uncertainty remaining in the total SFRs; however, the IR luminosities of the most massive protoclusters are likely to continue increasing up to $z\sim4$. Meanwhile, no significant IR flux excess is observed around optically selected QSOs at similar redshifts, which confirms previous results. Our results suggest that the $z\sim4$ protoclusters trace dense, intensely star-forming environments that may also host obscured AGNs missed by the selection in the optical.

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 present deep near-infrared spectroscopy of six quasars at 6.1<z<6.7 with VLT/X-Shooter and Gemini-N/GNIRS. Our objects, originally discovered through a wide-field optical survey with the Hyper Suprime-Cam (HSC) Subaru Strategic Program (HSC-SSP), have the lowest luminosities (-25.5< M1450<-23.1 mag) of the z>5.8 quasars with measured black hole masses. From single-epoch mass measurements based on MgII2798, we find a wide range in black hole masses, from M_BH=10^7.6 to 10^9.3 Msun. The Eddington ratios L_bol/L_Edd range from 0.16 to 1.1, but the majority of the HSC quasars are powered by M_BH=10^9 Msun supermassive black holes (SMBHs) accreting at sub-Eddington rates. The Eddington ratio distribution of the HSC quasars is inclined to lower accretion rates than those of Willott et al. (2010a), who measured the black hole masses for similarly faint z=6 quasars. This suggests that the global Eddington ratio distribution is wider than has previously been thought. The presence of M_BH=10^9 Msun SMBHs at z=6 cannot be explained with constant sub-Eddington accretion from stellar remnant seed black holes. Therefore, we may be witnessing the first buildup of the most massive black holes in the first billion years of the universe, the accretion activity of which is transforming from active growth to a quiescent phase. Measurements of a larger complete sample of z>6 low-luminosity quasars, as well as deeper observations with future facilities will enable us to better understand the early SMBH growth in the reionization epoch.

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

We have carried out deep and wide field imaging observations with narrow bands targeting 11 quasar fields to systematically study the possible photoevaporation effect of quasar radiation on surrounding low mass galaxies at $z\sim2-3$. We focused on Lyman alpha emitters (LAEs) at the same redshifts as the quasars that lie within the quasar proximity zones, where the ultra-violet radiation from the quasars is higher than the average background at that epoch. We found that LAEs with high rest-frame equivalent width of Ly$\alpha$ emission ($EW_0$) of $\gtrsim 150$\AA$~$ with low stellar mass ($\lesssim 10^8 M_{\odot}$), are predominantly scarce in the quasar proximity zones, suggesting that quasar photoevaporation effects may be taking place. The halo mass of LAEs with $EW_0>150$\AA$~$ is estimated to be $3.6^{+12.7}_{-2.3}\times10^9 M_{\odot}$ either from the Spectral Energy Distribution (SED) fitting or the main sequence. Based on a hydrodynamical simulation, the predicted delay in star formation under a local UV background intensity with $J (\nu_L)\gtrsim10^{-21}$ erg s$^{-1}$ cm$^{-2}$ Hz$^{-1}$ sr$^{-1}$ for galaxies having less than this halo mass is about $>20$ Myr, which is longer than the expected age of LAEs with $EW_0>150$\AA. On the other hand, the photoevaporation seems to be less effective around very luminous quasars, which is consistent with the idea that these quasars are still in an early stage of quasar activity.

We report the results of a neutrino search in Super-Kamiokande for coincident signals with the first detected gravitational wave produced by a binary neutron star merger, GW170817, which was followed by a short gamma-ray burst, GRB170817A, and a kilonova/macronova. We searched for coincident neutrino events in the range from 3.5 MeV to $\sim$100 PeV, in a time window $\pm$500 seconds around the gravitational wave detection time, as well as during a 14-day period after the detection. No significant neutrino signal was observed for either time window. We calculated 90% confidence level upper limits on the neutrino fluence for GW170817. From the upward-going-muon events in the energy region above 1.6 GeV, the neutrino fluence limit is $16.0^{+0.7}_{-0.6}$ ($21.3^{+1.1}_{-0.8}$) cm$^{-2}$ for muon neutrinos (muon antineutrinos), with an error range of $\pm5^{\circ}$ around the zenith angle of NGC4993, and the energy spectrum is under the assumption of an index of $-2$. The fluence limit for neutrino energies less than 100 MeV, for which the emission mechanism would be different than for higher-energy neutrinos, is also calculated. It is $6.6 \times 10^7$ cm$^{-2}$ for anti-electron neutrinos under the assumption of a Fermi-Dirac spectrum with average energy of 20 MeV.

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 the results of a search in the Super-Kamiokande (SK) detector for excesses of neutrinos with energies above a few GeV that are in the direction of the track events reported in IceCube. Data from all SK phases (SK-I through SK-IV) were used, spanning a period from April 1996 to April 2016 and corresponding to an exposure of 225 kilotonne-years . We considered the 14 IceCube track events from a data set with 1347 livetime days taken from 2010 to 2014. We use Poisson counting to determine if there is an excess of neutrinos detected in SK in a 10 degree search cone (5 degrees for the highest energy data set) around the reconstructed direction of the IceCube event. No significant excess was found in any of the search directions we examined. We also looked for coincidences with a recently reported IceCube multiplet event. No events were detected within a $\pm$ 500 s time window around the first detected event, and no significant excess was seen from that direction over the lifetime of SK.

Vacuum-ultraviolet (VUV) photoabsorption spectra were recorded of the $A\,{}^2\Sigma^+(v'=0)\leftarrow{}X\,{}^2\Pi(v''=0)$, $D\,{}^2\Sigma^-(v'=0)\leftarrow{}X\,{}^2\Pi(v''=0)$ and $D\,{}^2\Sigma^-(v'=1)\leftarrow{}X\,{}^2\Pi(v''=0)$ bands of the OH and OD radicals generated in a plasma-discharge source using synchrotron radiation as a background continuum coupled with the VUV Fourier-transform spectrometer on the DESIRS beamline of synchrotron SOLEIL. High-resolution spectra permitted the quantification of transition frequencies, relative $f$-values, and natural line broadening. The $f$-values were absolutely calibrated with respect to a previous measurement of $A\,{}^2\Sigma^+(v'=0)\leftarrow{}X\,{}^2\Pi(v''=0)$ ([wang1979]). Lifetime broadening of the excited $D\,{}^2\Sigma^-(v=0)$ and $D\,{}^2\Sigma^-(v=1)$ levels is measured for the first time and compared with previous experimental limits, and implies a lifetime 5 times shorter than a theoretical prediction ([van_der_loo2005]). A local perturbation of the $D\,{}^2\Sigma^-(v=0)$ level in OH was found.