Extensive air showers induced from high-energy cosmic rays provide a window into understanding the most energetic phenomena in the universe. We present a new method for observing these showers using the silicon imaging detector Subaru Hyper Suprime-Cam (HSC). This method has the advantage of being able to measure individual secondary particles. When paired with a surface detector array, silicon imaging detectors like Subaru HSC will be useful for studying the properties of extensive air showers in detail. The following report outlines the first results of observing extensive air showers with Subaru HSC. The potential for reconstructing the incident direction of primary cosmic rays is demonstrated and possible interdisciplinary applications are discussed.

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 ∼ 4 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 ∼ 2 to z ∼ 4. We carefully remove possible contaminants from foreground galaxies and, for each structure, 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 ∼28% 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 ∼ 4. These results suggest that specific environmental effects or assembly biases have already emerged in some protoclusters as early as z ∼ 4, and we suggest that proto-BCGs have different star formation histories than other galaxies in the same epoch.

We report the rest-frame ultraviolet luminosity function of g-dropout galaxies in 177 protocluster candidates (PC UVLF) at z ∼ 4 selected in the Hyper Suprime-Cam Subaru Strategic Program. Comparing it with the UVLF of field galaxies at the same redshift, we find that the PC UVLF shows a significant excess toward the bright end. This excess cannot be explained by the contribution of only active galactic nuclei, and we also find that this excess is more significant in higher density 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 ∼ 4. We also estimate the star formation rate density (SFRD) in protocluster regions as ≃6%-20% of the cosmic SFRD, based on the measured PC UVLF after correction for the selection incompleteness in our protocluster sample. This high value suggests that protoclusters make a nonnegligible contribution to the cosmic SFRD at z ∼ 4, as previously suggested by simulations. Our results suggest that protoclusters are essential components for galaxy evolution at z ∼ 4.

We present Lyα and UV-nebular emission line properties of bright Lyα emitters (LAEs) at z = 6-7 with a luminosity of log LLyα/[erg s-1] = 43-44 identified in the 21 deg2 area of the SILVERRUSH early sample developed with the Subaru Hyper Suprime-Cam survey data. Our optical spectroscopy newly confirms 21 bright LAEs with clear Lyα emission, and contributes to making a spectroscopic sample of 96 LAEs at z = 6-7 in SILVERRUSH. From the spectroscopic sample, we select seven remarkable LAEs as bright as Himiko and CR7 objects, and perform deep Keck/MOSFIRE and Subaru/nuMOIRCS near-infrared spectroscopy reaching the 3 σ flux limit of ∼2 × 10-18 erg s-1 for the UV-nebular emission lines of He II λ1640, CIV λλ1548,1550, and OIII]λλ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 (EW0) of ∼2-4 A for CIV, He II, and OIII] lines. We also investigate the VLT/X-SHOOTER spectrum of CR7 whose 6 σ detection of He II is claimed by Sobral et al. Although two individuals and the ESO archive service carefully reanalyzed 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. The 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 EW0 and the UV-continuum luminosity, which are similar to those found at z ∼ 2-3.

We present high-fidelity cosmology results from a blinded joint analysis of galaxy-galaxy weak lensing (Δς) and projected galaxy clustering (wp) measured from the Hyper Suprime-Cam Year-1 (HSC-Y1) data and spectroscopic Sloan Digital Sky Survey (SDSS) galaxy catalogs in the redshift range 0.15

We explored the clustering properties of Lyman break galaxies at z = 4 and 5 with an angular two-point correlation function on the basis of the very deep and wide Subaru Deep Field data. We confirmed the previous result that the clustering strength of LBGs depends on the UV luminosity in the sense that brighter LBGs are more strongly clustered. In addition, we found an apparent dependence of the correlation function slope on UV luminosity for LBGs at both z = 4 and 5. More luminous LBGs have a steeper correlation function. The bias parameter was found to be a scale-dependent function for bright LBGs, whereas it appears to be almost scale-independent for faint LBGs. Luminous LBGs have a higher bias at smaller angular scales, which decreases as the scale increases. To compare these observational results, we constructed numerical mock LBG catalogs based on a semianalytic model of hierarchical clustering combined with high-resolution N-body simulation, carefully mimicking the observational selection effects. The luminosity functions and the overall correlation functions for LBGs at z = 4 and 5 predicted by this mock catalog were found to be almost consistent with the observation. The observed dependence of the clustering on UV luminosity was not reproduced by the model, unless subsamples of distinct halo mass were considered. That is, LBGs belonging to more massive dark halos had steeper and larger amplitude correlation functions. With this model, we found that LBG multiplicity in massive dark halos amplifies the clustering strength at small scales, which steepens the correlation function. The hierarchical clustering model could therefore be reconciled with the observed luminosity dependence of the correlation function if there is a tight correlation between UV luminosity and halo mass. Our finding that the slope of the correlation function depends on luminosity could be an indication that massive dark halos hosted multiple bright LBGs. © 2006. The American Astronomical Society. All rights reserved.

We perform a blinded cosmology analysis with cosmic shear two-point correlation functions measured from more than 25 million galaxies in the Hyper Suprime-Cam three-year shear catalog in four tomographic redshift bins ranging from 0.3 to 1.5. After conservative masking and galaxy selection, the survey covers 416 deg2 of the northern sky with an effective galaxy number density of 15 arcmin-2 over the four redshift bins. The 2PCFs adopted for cosmology analysis are measured in the angular range; 7.1<θ/arcmin<56.6 for ζ+ and 31.2<θ/arcmin<248 for ζ-, with a total signal-to-noise ratio of 26.6. We apply a conservative, wide, flat prior on the photometric redshift errors on the last two tomographic bins, and the relative magnitudes of the cosmic shear amplitude across four redshift bins allow us to calibrate the photometric redshift errors. With this flat prior on redshift errors, we find ωm=0.256-0.044+0.056 and S8σ8ωm/0.3=0.769-0.034+0.031 (both 68% C.I.) for a flat Λ cold dark matter cosmology. We find, after unblinding, that our constraint on S8 is consistent with the Fourier space cosmic shear and the 3×2 pt analyses on the same HSC dataset. We carefully study the potential systematics from astrophysical and systematic model uncertainties in our fiducial analysis using synthetic data, and report no biases (including projection bias in the posterior space) greater than 0.5σ in the estimation of S8. Our analysis hints that the mean redshifts of the two highest tomographic bins are higher than initially estimated. In addition, a number of consistency tests are conducted to assess the robustness of our analysis. Comparing our result with Planck-2018 cosmic microwave background observations, we find a ∼2σ tension for the ΛCDM model.

This paper presents the third data release of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP), a wide-field multi-band imaging survey with the Subaru 8.2 m telescope. HSC-SSP has three survey layers (Wide, Deep, and UltraDeep) with different area coverages and depths, designed to address a wide array of astrophysical questions. This third release from HSC-SSP includes data from 278 nights of observing time and covers about 670 deg2 in all five broad-band filters (grizy) at the full depth (∼26 mag at 5σ depending on filter) in the Wide layer. If we include partially observed areas, the release covers 1470 deg2. The Deep and UltraDeep layers have ∼ 80% of the originally planned integration times, and are considered done, as we have slightly changed the observing strategy in order to compensate for various time losses. There are a number of updates in the image processing pipeline. Of particular importance is the change in the sky subtraction algorithm; we subtract the sky on small scales before the detection and measurement stages, which has significantly reduced the number of false detections. Thanks to this and other updates, the overall quality of the processed data has improved since the previous release. However, there are limitations in the data (for example, the pipeline is not optimized for crowded fields), and we encourage the user to check the quality assurance plots as well as a list of known issues before exploiting the data.