Subaru Hyper Suprime-Cam Survey Research Articles

Active Galactic Nucleus Properties of ∼1 Million Member Galaxies of Galaxy Groups and Clusters at z < 1.4 Based on the Subaru Hyper Suprime-Cam Survey

Herein, we present the statistical properties of active galactic nuclei (AGNs) for approximately 1 million member galaxies of galaxy groups and clusters with 0.1 < cluster redshift (z cl) < 1.4 selected using the Subaru Hyper Suprime-Cam, the so-called CAMIRA clusters. In this research, we focused on the AGN power fraction (f AGN), which is defined as the proportion of the contribution of AGNs to the total infrared (IR) luminosity, L IR (AGN)/L IR, and examined how f AGN depends on (i) z cl and (ii) the distance from the cluster center. We compiled multiwavelength data using the ultraviolet–mid-IR range. Moreover, we performed spectral energy distribution fits to determine f AGN using the CIGALE code with the SKIRTOR AGN model. We found that (i) the value of f AGN in the CAMIRA clusters is positively correlated with z cl, with the correlation slope being steeper than that for field galaxies, and (ii) f AGN exhibits a high value at the cluster outskirts. These results indicate that the emergence of the AGN population depends on the redshift and environment and that galaxy groups and clusters at high redshifts are important in AGN evolution. Additionally, we demonstrated that cluster–cluster mergers may enhance AGN activity at the outskirts of particularly massive galaxy clusters. Our findings are consistent with a related study on the CAMIRA clusters that was based on the AGN number fraction.

Cosmology from weak lensing peaks and minima with Subaru Hyper Suprime-Cam Survey first-year data

ABSTRACT We present cosmological constraints derived from peak counts, minimum counts, and the angular power spectrum of the Subaru Hyper Suprime-Cam first-year (HSC Y1) weak lensing shear catalogue. Weak lensing peak and minimum counts contain non-Gaussian information and hence are complementary to the conventional two-point statistics in constraining cosmology. In this work, we forward-model the three summary statistics and their dependence on cosmology, using a suite of N-body simulations tailored to the HSC Y1 data. We investigate systematic and astrophysical effects including intrinsic alignments, baryon feedback, multiplicative bias, and photometric redshift uncertainties. We mitigate the impact of these systematics by applying cuts on angular scales, smoothing scales, signal-to-noise ratio bins, and tomographic redshift bins. By combining peaks, minima, and the power spectrum, assuming a flat-ΛCDM model, we obtain $S_{8} \equiv \sigma _8\sqrt{\Omega _m/0.3}= 0.810^{+0.022}_{-0.026}$, a 35 per cent tighter constraint than that obtained from the angular power spectrum alone. Our results are in agreement with other studies using HSC weak lensing shear data, as well as with Planck 2018 cosmology and recent CMB lensing constraints from the Atacama Cosmology Telescope and the South Pole Telescope.

Galaxy–dark matter connection of photometric galaxies from the HSC-SSP Survey: galaxy–galaxy lensing and the halo model

ABSTRACT We infer the connection between stellar mass of galaxies from Subaru Hyper Suprime-Cam (HSC) survey, and their dark matter halo masses in two bins of redshifts between [0.3, 0.8]. We use measurements of the weak lensing of background galaxies from Year-1 shape catalogue from the HSC survey. We bin lens galaxies in stellar mass with varying thresholds ranging from $8.6\le \log [M_*/(h^{-2}{M_\odot })]\le \, 11.2$ and use stringent cuts in the selection of source galaxies to measure the lensing signal. We model these measurements of the lensing signals together with the abundance of galaxies in the halo occupation distribution framework. We obtain constraints on the halo occupation parameters of central galaxies Mmin and σlog M along with parameters that describe the occupation of the satellite galaxies. The measurements of abundance and lensing individually constrain different degeneracy directions in Mmin− σlog M plane, thus breaking the degeneracy in these parameters. We demonstrate that the lensing measurements are best able to constrain average central halo masses, 〈Mcen〉. We compare our measurements to those obtained using the abundance and clustering of these galaxies. We find that the galaxy–dark matter connection does not vary significantly between redshift bins we explore in this study. Uncertainties in the photometric redshift of the lens galaxies need to be studied to understand the true underlying stellar mass–halo mass relation and evolution of galaxies.

Subhalo abundance matching using progenitor mass at varying redshift: two modes of stellar mass growth imprinted into the Subaru HSC galaxy clustering

ABSTRACT We propose a novel subhalo abundance matching (SHAM) model that uses the virial mass of the main progenitor of each (sub)halo Mprog as a proxy of the galaxy stellar mass M* at the time of observation. This Mprog model predicts the two-point correlation functions depending on the choice of the epoch zprog at which Mprog is quoted. With zprog as a fitting parameter, we apply the Mprog model to the angular correlation functions measured with varying stellar mass thresholds from M*, lim/(h−2 M⊙) = 1011–108.6 using a sample of galaxies at z ≃ 0.4 from the Subaru Hyper Suprime-Cam survey. The Mprog model can reproduce the observations very well over 10 h−1kpc–10 h−1Mpc. We find that, for the samples of 109.2 ≤ M*, lim/(h−2 M⊙) ≤ 1010.2, the correlation functions predicted by the widely used Vpeak model lack amplitudes at ≲1 h−1 Mpc, suggesting that Mprog is a better proxy of the galaxy stellar mass than conventional Vpeak. The zprog parameter is highest (zprog ≃ 3) for intermediate mass galaxies at M* ≃ 109.9 h−2 M⊙, and becomes smaller down to zprog ≃ 1 for both lower and higher mass galaxies. We interpret these trends as reflecting the downsizing in the in situ star formation in lower mass galaxies and the larger contribution of the ex situ stellar mass growth in higher mass galaxies.

The nature and origins of the low surface brightness outskirts of massive, central galaxies in Subaru HSC

ABSTRACT We explore the stellar mass density and colour profiles of 118 low redshift, massive, central galaxies, selected to have assembled 90 per cent of their stellar mass 6 Gyr ago, finding evidence of the minor merger activity expected to be the driver behind the size growth of quiescent galaxies. We use imaging data in the g, r, i, z, y bands from the Subaru Hyper Suprime-Cam survey and perform SED fitting to construct spatially well-resolved radial profiles in colour and stellar mass surface density. Our visual morphological classification reveals that ∼42 per cent of our sample displays tidal features, similar to previous studies, ∼43 per cent of the remaining sample displays a diffuse stellar halo, and only ∼14 per cent displays no features, down to a limiting μr-band ∼ 28 mag arcsec−2. We find good agreement between the stacked colour profiles of our sample to those derived from previous studies and an expected smooth, declining stellar mass surface density profile in the central regions (&amp;lt; 3 Re). However, we also see a flattening of the profile (Σ* ∼ 107.5 M⊙ kpc−2) in the outskirts (up to 10 Re), which is revealed by our method of specifically targeting tidal/accretion features. We find similar levels of tidal features and behaviour in the stellar mass surface density profiles in a younger comparison sample, however, a lack of diffuse haloes. We also apply stacking techniques similar to those in previous studies, finding such procedures wash out tidal features and thereby produce smooth declining profiles. The stellar material in the outskirts contributes on average ∼1010 M⊙ or a few per cent of the total stellar mass and has similar colours to SDSS satellites of similar stellar mass.

The eROSITA Final Equatorial-Depth Survey (eFEDS)

Context. We present the results of a systematic X-ray analysis of optically rich galaxy clusters detected by the Subaru Hyper Suprime-Cam (HSC) survey in the eROSITA Final Equatorial-Depth Survey (eFEDS) field. Aims. Through a joint analysis of the SRG (Spectrum Roentgen Gamma)/eROSITA and Subaru/HSC surveys, we aim to investigate the dynamical status of the optically selected clusters and to derive the cluster scaling relations. Methods. The sample consists of 43 optically selected galaxy clusters with a richness &gt;40 in the redshift range of 0.16–0.89. We systematically analyzed the X-ray images and emission spectra using the eROSITA data. We identified the brightest cluster galaxy (BCG) using the optical and far-infrared databases. We evaluated the cluster’s dynamical status by measuring three quantities: offset between the X-ray peak and BCG position, the gas concentration parameter, and the number of galaxy-density peaks. We investigated the luminosity–temperature and mass–luminosity relations based on eROSITA X-ray spectra and HSC weak-lensing data analyses. Results. Based on these three measurements, we estimated the fraction of relaxed clusters to be 2(&lt; 39)%, which is smaller than that of the X-ray-selected cluster samples. After correcting for a selection bias due to the richness cut, we obtained a shallow L−T slope of 2.1 ± 0.5, which is consistent with the predictions of the self-similar model and the baseline model incorporating a mass–concentration relation. The L−M slope of 1.5 ± 0.3 is in agreement with the above-cited theoretical models as well as the data on the shear-selected clusters in the eFEDs field. Conclusions. Our analysis of high-richness optical clusters yields a small fraction of relaxed clusters and a shallow slope for the luminosity–temperature relation. This suggests that the average X-ray properties of the optical clusters are likely to be different from those observed in the X-ray samples. Thus, the joint eROSITA and HSC observations are a powerful tool in extending the analysis to a larger sample and understanding the selection effect with a view to establish cluster scaling relations.

A tale of a tail: a tidally disrupting ultra-diffuse galaxy in the M81 group

ABSTRACT We present the discovery of a giant tidal tail of stars associated with F8D1, the closest known example of an ultra-diffuse galaxy (UDG). F8D1 sits in a region of the sky heavily contaminated by Galactic cirrus and has been poorly studied since its discovery two decades ago. The tidal feature was revealed in a deep map of resolved red giant branch stars constructed using data from our Subaru Hyper Suprime-Cam survey of the M81 Group. It has an average surface brightness of μg ∼ 32 mag arcsec−2 and can be traced for over a degree on the sky (60 kpc at the distance of F8D1) with our current imagery. We revisit the main body properties of F8D1 using deep multiband imagery acquired with MegaCam on CFHT and measure effective radii of 1.7–1.9 kpc, central surface brightnesses of 24.7–25.7 mag, and a stellar mass of ∼7 × 107M⊙. Assuming a symmetric feature on the other side of the galaxy, we calculate that 30–36 per cent of F8D1’s present-day luminosity is contained in the tail. We argue that the most likely origin of F8D1’s disruption is a recent close passage to M81, which would have stripped its gas and quenched its star formation. As the only UDG that has so far been studied to such faint surface brightness depths, the unveiling of F8D1’s tidal disruption is important. It leaves open the possibility that many other UDGs could be the result of similar processes, with the most telling signatures of this lurking below current detection limits.

Subaru HSC weak lensing of SDSS redMaPPer cluster satellite galaxies: empirical upper limit on orphan fractions

ABSTRACT Weak gravitational lensing directly probes the matter distribution surrounding satellite galaxies in galaxy clusters. We measure the weak lensing signal induced on the shapes of background galaxies around the Sloan Digital Sky Survey (SDSS) redMaPPer cluster satellite galaxies, which have their central galaxies assigned with a probability Pcen &amp;gt; 0.95 in the redshift range, 0.1 ≤ z ≤ 0.33. We use the galaxy shapes from the Subaru Hyper Suprime-Cam (HSC) survey for this purpose. We bin satellite galaxies by their distance from the cluster centre and compare it to the signal around a control sample of galaxies that do not reside in clusters but have similar colours and magnitudes. We explore the effect of environmental processes on the dark matter mass around satellites. We see hints of a difference in the mass of the subhalo of the satellite compared to the halo masses of galaxies in our control sample, especially in the innermost cluster-centric radial bin (0.1 &amp;lt; r &amp;lt; 0.3 [${h^{-1}\,\rm Mpc}$]). For the first time, we put an upper limit on the prevalence of orphan galaxies that have entirely lost their dark matter haloes with cluster-centric distances with the help of our measurements. However, these upper limits could be relaxed if there is substantial contamination in the satellite galaxy sample.

The statistical properties of 28 IR-bright dust-obscured galaxies and SED modelling using CIGALE

Abstract The aim of this study is to characterize the physical and statistical properties of a sample of infrared-bright dust-obscured galaxies (DOGs) by fitting their spectral energy distributions (SEDs). We examined 28 DOGs at redshifts 0.47 ≤ z ≤ 1.63 discovered by combining images of the Subaru Hyper Suprime-Cam (HSC) survey, the VISTA Kilo-degree Infrared Galaxy (VIKING) survey, and the Wide-field Infrared Survey Explorer (WISE) all-sky survey, and detected at Herschel Spectral and Photometric Imaging Receiver (SPIRE) bands. We have detected a significant active galactic nucleus (AGN) contribution to the mid-infrared luminosity for $71\%$ of DOGs. Our DOGs contain several types of AGNs; the majority of AGN emission originates from Type 2 AGNs. Our DOG sample shows very high values of stellar mass [log (M*/M${}_\odot $) = 11.49 ± 1.61] compared with other samples of DOGs selected at infrared wavelengths. Our study is promising to identify a new type of DOGs called “overweight” DOGs (ODOGs). ODOGs may indicate the end of the DOG phase, and then they become visible quasars. Principal component (PC) analysis is applied to reduce the number of dimensions of our sample, removing the dependency on the observed variables. There are two significant PCs describing $72.7\%$ of the total variance. The first PC strongly correlates with redshift, dust luminosity, dust mass, and stellar mass, while far-ultraviolet (FUV) attenuation strongly correlates with the second PC, which is orthogonal to the first one. The partial correlation between the resulted physical parameters is tested, supporting the reliability of the correlations.

The Subaru HSC weak lensing mass-observable scaling relations of spectroscopic galaxy groups from the GAMA survey

ABSTRACT We utilize the galaxy shape catalogue from the first-year data release of the Subaru Hyper Suprime-Cam (HSC) survey to study the dark matter content of galaxy groups in the Universe using weak lensing. We use galaxy groups from the Galaxy Mass and Assembly galaxy survey in approximately 100 sq. degrees of the sky that overlap with the HSC survey as lenses. We restrict our analysis to the 1587 groups with at least five members. We divide these groups into six bins each of group luminosity and group member velocity dispersion and measure the lensing signal with a signal-to-noise ratio of 55 and 51 for these two different selections, respectively. We use a Bayesian halo model framework to infer the halo mass distribution of our groups binned in the two different observable properties and constrain the power-law scaling relation and the scatter between mean halo masses and the two-group observable properties. We obtain a 5 per cent constraint on the amplitude of the scaling relation between halo mass and group luminosity with 〈M〉 = (0.81 ± 0.04) × 1014 h−1 M⊙ for Lgrp = 1011.5 h−2 L⊙, and a power-law index of α = 1.01 ± 0.07. We constrain the amplitude of the scaling relation between halo mass and velocity dispersion to be 〈M〉 = (0.93 ± 0.05) × 1014 h−1 M⊙ for $\sigma = 500\, {\rm km\, s}^{-1}$ and a power-law index to be α = 1.52 ± 0.10. However, these scaling relations are sensitive to the exact cuts applied to the number of group members. Comparisons with similar scaling relations from the literature show that our results are consistent and have significantly reduced errors.

Can QCD axion stars explain Subaru HSC microlensing?

A non-negligible fraction of the QCD axion dark matter may form gravitationally bound Bose Einstein condensates, which are commonly known as axion stars or axion clumps. Such astrophysical objects have been recently proposed as the cause for the single candidate event reported by Subaru Hyper Suprime-Cam (HSC) microlensing search in the Andromeda galaxy. Depending on the breaking scale of the Peccei-Quinn symmetry and the details of the dark matter scenario, QCD axion clumps may form via gravitational condensation during radiation domination, in the dense core of axion miniclusters, or within axion minihalos around primordial black holes. We analyze all these scenarios and conclude that the microlensing candidate detected by the Subaru HSC survey is likely not caused by QCD axion stars.

Subaru Hyper Suprime-Cam Survey of Cygnus OB2 Complex – I. Introduction, photometry, and source catalogue

ABSTRACT Low-mass star formation inside massive clusters is crucial to understand the effect of cluster environment on processes like circumstellar disc evolution, planet, and brown dwarf formation. The young massive association of Cygnus OB2, with a strong feedback from massive stars, is an ideal target to study the effect of extreme environmental conditions on its extensive low-mass population. We aim to perform deep multiwavelength studies to understand the role of stellar feedback on the IMF, brown dwarf fraction and circumstellar disc properties in the region. We introduce here, the deepest and widest optical photometry of 1.5○ diameter region centred at Cygnus OB2 in r2, i2, z, and Y-filters, using Subaru Hyper Suprime-Cam (HSC). This work presents the data reduction, source catalogue generation, data quality checks, and preliminary results about the pre-main sequence sources. We obtain 713 529 sources in total, with detection down to ∼28, 27, 25.5, and 24.5 mag in r2, i2, z, and Y-band, respectively, which is ∼3 – 5 mag deeper than the existing Pan-STARRS and GTC/OSIRIS photometry. We confirm the presence of a distinct pre-main sequence branch by statistical field subtraction of the central 18 arcmin region. We find the median age of the region as ∼5 ± 2 Myr with an average disc fraction of ∼9 per cent. At this age, combined with A $_V\, \sim$ 6 – 8 mag, we detect sources down to a mass range of ∼0.01–0.17 M⊙. The deep HSC catalogue will serve as the groundwork for further studies on this prominent active young cluster.

Ionizing radiation from AGNs at z &amp;gt; 3.3 with the Subaru Hyper Suprime-Cam Survey and the CFHT Large Area U-band Deep Survey (CLAUDS)

ABSTRACT We use deep and wide imaging data from the CFHT Large Area U-band Deep Survey (CLAUDS) and the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) to constrain the ionizing radiation (Lyman Continuum, LyC) escape fraction from active galactic nuclei (AGNs) at z ∼ 3–4. For 94 AGNs with spectroscopic redshifts at 3.3 &amp;lt; z &amp;lt; 4.0, we use their U/i-band flux ratios to estimate LyC transmission of individual AGNs. The distribution of their LyC transmission shows values lower than the range of LyC transmission values for IGM of the same redshift range, which suggests that LyC escape fraction of AGNs at z &amp;gt; 3.3 is considerably lower than unity in most cases. We do not find any trend in LyC transmission values depending on their UV luminosities. Based on the photometry of stacked images we find the average flux ratio of LyC and non-ionizing UV photons escaping from the objects (fLyC/fUV)out = 0.182 ± 0.043 for AGNs at 3.3 &amp;lt; z &amp;lt; 3.6, which corresponds to LyC escape fraction fesc = 0.303 ± 0.072 if we assume a fiducial intrinsic SED of AGN. Based on the estimated LyC escape fraction and the UV luminosity function (UVLF) of AGNs, we argue that UV-selected AGNs’ contribution to the LyC emissivity at the epoch is minor, although the size of their contribution largely depends on the shape of the UVLF.

Subaru High-z Exploration of Low-luminosity Quasars (SHELLQs). XIV. A Candidate Type II Quasar at z = 6.1292

Abstract We present deep Keck/MOSFIRE near-infrared spectroscopy of a strong Lyα emitting source at z = 6.1292, HSC J142331.71−001809.1, which was discovered by the SHELLQS program from imaging data of the Subaru Hyper Suprime-Cam (HSC) survey. This source is one of five objects that show narrow (FWHM &lt; 500 km s−1) and prominent (L Lyα &gt; 1044 erg s−1) Lyα emission lines at absolute 1450 Å continuum magnitudes of M 1450 ∼ −22 mag. Its rest-frame Lyα equivalent width (EW) is 370 ± 30 Å. In the 2 hr Keck/MOSFIRE spectrum in Y band, the high-ionization C iv λλ1548,1550 doublet emission line was clearly detected with FWHM = 120 − 20 + 20 km s−1 and a total rest-frame EW of 37 − 5 + 6 Å. We also report the detection of weak continuum emission, and the tentative detection of O iii] λλ1661,1666 in the 4 hr J-band spectrum. Judging from the UV magnitude, line widths, luminosities, and EWs of Lyα and C iv, we suggest that this source is a reionization-era analog of classical type-II AGNs, although there is a possibility that it represents a new population of AGN/galaxy composite objects in the early universe. We compare the properties of J1423−0018 to intermediate-redshift type-II AGNs and C iv emitters seen in z = 6–7 galaxy samples. Further observations of other metal emission lines in the rest-frame UV or optical, along with X-ray follow-up observations of the z = 6–7 narrow-line quasars, are needed for more robust diagnostics and to determine their nature.

The physical properties of IR-bright Dust-Obscured Galaxies

AbstractThis work aims to represent the physical properties of a sample of infrared-bright dust-obscured galaxies (DOGs) studied by Suleiman et al. (2022) by fitting the spectral energy distributions (SEDs). Twenty-eight DOGs were examined at redshifts 0.47 ≤ z ≤ 1.63 discovered by combining images of the Subaru Hyper Suprime-Cam (HSC) survey, VISTA Kilo-degree Infrared Galaxy (VIKING) survey, and the Wide-field Infrared Survey Explorer (WISE) all-sky survey, and detected at Herschel Spectral and Photometric Imaging Receiver (SPIRE) bands. The results show a correlation between the star formation rate (SFR) and the dust luminosity of Suleiman et al. (2022) DOG sample, the SFR ranges of the sample according to different redshifts, and a comparison between Suleiman et al. (2022) sample and other samples of DOGs.

Three-dimensional Reconstruction of Weak-lensing Mass Maps with a Sparsity Prior. I. Cluster Detection

Abstract We propose a novel method to reconstruct high-resolution three-dimensional mass maps using data from photometric weak-lensing surveys. We apply an adaptive LASSO algorithm to perform a sparsity-based reconstruction on the assumption that the underlying cosmic density field is represented by a sum of Navarro–Frenk–White halos. We generate realistic mock galaxy shear catalogs by considering the shear distortions from isolated halos for the configurations matched to the Subaru Hyper Suprime-Cam Survey with its photometric redshift estimates. We show that the adaptive method significantly reduces line-of-sight smearing that is caused by the correlation between the lensing kernels at different redshifts. Lensing clusters with lower mass limits of 1014.0 h−1 M ⊙, 1014.7 h−1 M ⊙, 1015.0 h−1 M ⊙ can be detected with 1.5σ confidence at the low (z &lt; 0.3), median (0.3 ≤ z &lt; 0.6), and high (0.6 ≤ z &lt; 0.85) redshifts, respectively, with an average false detection rate of 0.022 deg−2. The estimated redshifts of the detected clusters are systematically lower than the true values by Δz ∼ 0.03 for halos at z ≤ 0.4, but the relative redshift bias is below 0.5% for clusters at 0.4 &lt; z ≤ 0.85. The standard deviation of the redshift estimation is 0.092. Our method enables direct three-dimensional cluster detection with accurate redshift estimates.

The Evolution of the Ultraluminous Lyα Luminosity Function over z = 5.7–6.6

Abstract Ultraluminous Lyα emitting galaxies (ULLAEs) with erg s−1 near the epoch of reionization (z &gt; 5) make up the bright end of the LAE luminosity function (LF) and may provide insight into the process of reionization, including the formation of ionized bubbles around these extreme systems. We present a spectroscopic LF for ULLAEs at z = 5.7. We used data from the HEROES ∼45 deg2 Subaru Hyper Suprime-Cam survey, which is centered on the North Ecliptic Pole and has both broadband (grizY) and narrowband (NB816 and NB921) imaging, to select candidate ULLAEs based on a NB816 excess and a strong broadband Lyman break. We spectroscopically observed 17 ULLAE candidates with DEIMOS on Keck II. We confirmed 12 as LAEs at z = 5.7, 9 of which are ULLAEs. The remaining sources are an active galactic nuclei at z = 5.7, an [O III]λ5007 emitter at z = 0.63, a red star, and two spectroscopic nondetections. Using the nine confirmed ULLAEs, we construct a ULLAE LF at z = 5.7. After applying a comprehensive incompleteness correction, we compare our new z = 5.7 LF with our recent z = 6.6 LF and with other LFs from the literature to look for evolution at the ultraluminous end. We find the overall ratio of the z = 5.7 to z = 6.6 ULLAE comoving number densities to be , which corresponds to an LF offset of dex.

Noise reduction for weak lensing mass mapping: an application of generative adversarial networks to Subaru Hyper Suprime-Cam first-year data

ABSTRACT We propose a deep-learning approach based on generative adversarial networks (GANs) to reduce noise in weak lensing mass maps under realistic conditions. We apply image-to-image translation using conditional GANs to the mass map obtained from the first-year data of Subaru Hyper Suprime-Cam (HSC) Survey. We train the conditional GANs by using 25 000 mock HSC catalogues that directly incorporate a variety of observational effects. We study the non-Gaussian information in denoised maps using one-point probability distribution functions (PDFs) and also perform matching analysis for positive peaks and massive clusters. An ensemble learning technique with our GANs is successfully applied to reproduce the PDFs of the lensing convergence. About $60{{\ \rm per\ cent}}$ of the peaks in the denoised maps with height greater than 5σ have counterparts of massive clusters within a separation of 6 arcmin. We show that PDFs in the denoised maps are not compromised by details of multiplicative biases and photometric redshift distributions, nor by shape measurement errors, and that the PDFs show stronger cosmological dependence compared to the noisy counterpart. We apply our denoising method to a part of the first-year HSC data to show that the observed mass distribution is statistically consistent with the prediction from the standard ΛCDM model.

The Saga of M81: Global View of a Massive Stellar Halo in Formation

Abstract Recent work has shown that Milky Way–mass galaxies display an incredible range of stellar halo properties, yet the origin of this diversity is unclear. The nearby galaxy M81—currently interacting with M82 and NGC 3077—sheds unique light on this problem. We present a Subaru Hyper Suprime-Cam survey of the resolved stellar populations around M81, revealing M81's stellar halo in never-before-seen detail. We resolve the halo to unprecedented V-band equivalent surface brightnesses of 33 mag arcsec −2 and produce the first-ever global stellar mass density map for a Milky Way–mass stellar halo outside of the Local Group. Using the minor axis, we confirm M81's halo as one of the lowest mass and metal poorest known (M ⋆ ≃ 1.16 × 109 M ⊙, [Fe/H] ≃ −1.2)—indicating a relatively quiet prior accretion history. Yet, our global halo census finds that tidally unbound material from M82 and NGC 3077 provides a substantial infusion of metal-rich material (M ⋆ ≃ 5.4 × 108 M ⊙, [Fe/H] ≃−0.9). We further show that, following the accretion of its massive satellite M82 (and the LMC-like NGC 3077), M81 will host one of the most massive and metal-rich stellar halos in the nearby universe. Thus, the saga of M81: following a passive history, M81's merger with M82 will completely transform its halo from a low-mass, anemic halo rivaling the Milky Way, to a metal-rich behemoth rivaled only by systems such as M31. This dramatic transformation indicates that the observed diversity in stellar halo properties is primarily driven by diversity in the largest mergers these galaxies have experienced.

The Faint End of the Quasar Luminosity Function at z ∼ 5 from the Subaru Hyper Suprime-Cam Survey

Abstract We present the quasar luminosity function at z ∼ 5 derived from the optical wide-field survey data obtained as a part of the Subaru strategic program (SSP) with the Hyper Suprime-Cam (HSC). From a ∼81.8 deg2 area in the Wide layer of the HSC-SSP survey, we selected 224 candidates of low-luminosity quasars at z ∼ 5 by adopting the Lyman-break method down to i = 24.1 mag. Based on our candidates and spectroscopically confirmed quasars from the Sloan Digital Sky Survey (SDSS), we derived the quasar luminosity function at z ∼ 5, covering a wide luminosity range of −28.76 &lt; M 1450 &lt; −22.32 mag. We found that the quasar luminosity function is fitted by a double power-law model with a break magnitude of mag. The inferred number density of low-luminosity quasars is lower, and the derived faint-end slope, , is flatter than those of previous studies at z ∼ 5. A compilation of the quasar luminosity function at 4 ≤ z ≤ 6 from the HSC-SSP suggests that there is little redshift evolution in the break magnitude and in the faint-end slope within this redshift range, although previous studies suggest that the faint-end slope becomes steeper at higher redshifts. The number density of low-luminosity quasars decreases more rapidly from z ∼ 5 to z ∼ 6 than from z ∼ 4 to z ∼ 5.