Showing 1–50 of 194 results for author: Takada, M

Weak-lensing Shear-selected Galaxy Clusters from the Hyper Suprime-Cam Subaru Strategic Program: II. Cosmological Constraints from the Cluster Abundance

Authors: I-Non Chiu, Kai-Feng Chen, Masamune Oguri, Markus M. Rau, Hironao Miyatake, Satoshi Miyazaki, Surhud More, Takashi Hamana, Tomomi Sunayama, Sunao Sugiyama, Masahiro Takada

Abstract: We present cosmological constraints using the abundance of weak-lensing shear-selected galaxy clusters in the Hyper Suprime-Cam (HSC) Subaru Strategic Program. The clusters are selected on the mass maps constructed using the three-year (Y3) weak-lensing data with an area of $\approx500~$deg$^2$, resulting in a sample size of $129$ clusters with high signal-to-noise ratios $ν$ of $ν\geq4.7$. Owing… ▽ More We present cosmological constraints using the abundance of weak-lensing shear-selected galaxy clusters in the Hyper Suprime-Cam (HSC) Subaru Strategic Program. The clusters are selected on the mass maps constructed using the three-year (Y3) weak-lensing data with an area of $\approx500~$deg$^2$, resulting in a sample size of $129$ clusters with high signal-to-noise ratios $ν$ of $ν\geq4.7$. Owing to the deep, wide-field, and uniform imaging of the HSC survey, this is by far the largest sample of shear-selected clusters, in which the selection solely depends on gravity and is free from any assumptions about the dynamical state. Informed by the optical counterparts, the shear-selected clusters span a redshift range of $z\lesssim0.7$ with a median of $z\approx0.3$. The lensing sources are securely selected at $z\gtrsim0.7$ with a median of $z\approx1.3$, leading to nearly zero cluster member contamination. We carefully account for (1) the bias in the photometric redshift of sources, (2) the bias and scatter in the weak-lensing mass using a simulation-based calibration, and (3) the measurement uncertainty that is directly estimated on the mass maps using an injection-based method developed in a companion paper (Chen et al. submitted). In a blind analysis, the fully marginalized posteriors of the cosmological parameters are obtained as $Ω_{\mathrm{m}} = 0.50^{+0.28}_{-0.24}$, $σ_8 = 0.685^{+0.161}_{-0.088}$, $\hat{S}_{8}\equivσ_8\left(Ω_{\mathrm{m}}/0.3\right)^{0.25} = 0.835^{+0.041}_{-0.044}$, and $σ_8\sqrt{Ω_{\mathrm{m}}/0.3} = 0.993^{+0.084}_{-0.126}$ in a flat $Λ$CDM model. We compare our cosmological constraints with other studies, including those based on cluster abundances, galaxy-galaxy lensing and clustering, and Cosmic Microwave Background observed by $Planck$, and find good agreement at levels of $\lesssim2σ$. [abridged] △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: 26 pages, 14 Figures, 1 Table. To be submitted to the Open Journal of Astrophysics. Main results are Figures 11 and 12. The cluster catalog and the cosmological chain are accessible at https://github.com/inonchiu/hsc_shear_selected_clusters

Weak-Lensing Shear-Selected Galaxy Clusters from the Hyper Suprime-Cam Subaru Strategic Program: I. Cluster Catalog, Selection Function and Mass--Observable Relation

Authors: Kai-Feng Chen, I-Non Chiu, Masamune Oguri, Yen-Ting Lin, Hironao Miyatake, Satoshi Miyazaki, Surhud More, Takashi Hamana, Markus M. Rau, Tomomi Sunayama, Sunao Sugiyama, Masahiro Takada

Abstract: We present the first step towards deriving cosmological constraints through the abundances of galaxy clusters selected in a $510\,\mathrm{deg}^2$ weak-lensing aperture mass map, constructed with the Year-Three shear catalog from the Hyper Suprime-Cam Subaru Strategic Program. We adopt a conservative source galaxy selection to construct a sample of $129$ weak-lensing peaks with a signal-to-noise ra… ▽ More We present the first step towards deriving cosmological constraints through the abundances of galaxy clusters selected in a $510\,\mathrm{deg}^2$ weak-lensing aperture mass map, constructed with the Year-Three shear catalog from the Hyper Suprime-Cam Subaru Strategic Program. We adopt a conservative source galaxy selection to construct a sample of $129$ weak-lensing peaks with a signal-to-noise ratio above $4.7$. We use semi-analytical injection simulations to derive the selection function and the mass--observable relation of our sample. These results take into account complicated uncertainties associated with weak-lensing measurements, such as the non-uniform survey depth and the complex survey geometry, projection effects from uncorrelated large-scale structures, and the intrinsic alignment of source galaxies. We also propose a novel modeling framework to make parts of the mass--observable relation insensitive to assumed cosmological parameters. Such a framework not only offers a great computational advantage to cosmological studies, but can also benefit future astrophysical studies using shear-selected clusters. Our results are an important step towards utilizing these cluster samples that are constructed nearly independent of any baryonic assumptions in upcoming deep-and-wide lensing surveys from the Vera Rubin Observatory, Euclid, and the Nancy Grace Roman Space Telescope. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: 17 pages, 9 figures, submitted to OJA. The cluster catalog of the cosmological sample (signal-to-noise ratio above $4.7$) is accessible at https://github.com/inonchiu/hsc_shear_selected_clusters

Exploring the baryonic effect signature in the Hyper Suprime-Cam Year 3 cosmic shear two-point correlations on small scales: the $S_8$ tension remains present

Authors: Ryo Terasawa, Xiangchong Li, Masahiro Takada, Takahiro Nishimichi, Satoshi Tanaka, Sunao Sugiyama, Toshiki Kurita, Tianqing Zhang, Masato Shirasaki, Ryuichi Takahashi, Hironao Miyatake, Surhud More, Atsushi J. Nishizawa

Abstract: The baryonic feedback effect is considered as a possible solution to the so-called $S_8$ tension indicated in cosmic shear cosmology. The baryonic effect is more significant on smaller scales, and affects the cosmic shear two-point correlation functions (2PCFs) with different scale- and redshift-dependencies from those of the cosmological parameters. In this paper, we use the Hyper Suprime-Cam Yea… ▽ More The baryonic feedback effect is considered as a possible solution to the so-called $S_8$ tension indicated in cosmic shear cosmology. The baryonic effect is more significant on smaller scales, and affects the cosmic shear two-point correlation functions (2PCFs) with different scale- and redshift-dependencies from those of the cosmological parameters. In this paper, we use the Hyper Suprime-Cam Year 3 (HSC-Y3) data to measure the cosmic shear 2PCFs ($ξ_{\pm}$) down to 0.28 arcminutes, taking full advantage of the high number density of source galaxies in the deep HSC data, to explore a possible signature of the baryonic effect. While the published HSC analysis used the cosmic shear 2PCFs on angular scales, which are sensitive to the matter power spectrum at $k\lesssim 1~h{\rm Mpc}^{-1}$, the smaller scale HSC cosmic shear signal allows us to probe the signature of matter power spectrum up to $k\simeq 20~h{\rm Mpc}^{-1}$. Using the accurate emulator of the nonlinear matter power spectrum, DarkEmulator2, we show that the dark matter-only model can provide an acceptable fit to the HSC-Y3 2PCFs down to the smallest scales. In other words, we do not find any clear signature of the baryonic effects or do not find a systematic shift in the $S_8$ value with the inclusion of the smaller-scale information as would be expected if the baryonic effect is significant. Alternatively, we use a flexible 6-parameter model of the baryonic effects, which can lead to both enhancement and suppression in the matter power spectrum compared to the dark matter-only model, to perform the parameter inference of the HSC-Y3 2PCFs. We find that the small-scale HSC data allow only a fractional suppression of up to 5 percent in the matter power spectrum at $k\sim 1~h{\rm Mpc}^{-1}$, which is not sufficient to reconcile the $S_8$ tension. △ Less

Submitted 29 March, 2024; originally announced March 2024.

Comments: 30 pages, 16 figures

Exploring faint white dwarfs and the luminosity function with Subaru HSC and SDSS in Stripe 82

Authors: Tian Qiu, Masahiro Takada, Naoki Yasuda, Akira Tokiwa, Kazumi Kashiyama, Yoshihisa Suzuki, Kenta Hotokezaka

Abstract: We present 4,987 white dwarf (WD) candidates selected from matched stars between the multi-band imaging datasets of Subaru Hyper Suprime-Cam (HSC) Survey and SDSS in the Stripe82 region covering about 165 deg$^2$. We first select WD candidates from the "reduced proper motion" diagram that is obtained by combining the apparent magnitude in the range $i=19$ -- 24 and the proper motion measured by co… ▽ More We present 4,987 white dwarf (WD) candidates selected from matched stars between the multi-band imaging datasets of Subaru Hyper Suprime-Cam (HSC) Survey and SDSS in the Stripe82 region covering about 165 deg$^2$. We first select WD candidates from the "reduced proper motion" diagram that is obtained by combining the apparent magnitude in the range $i=19$ -- 24 and the proper motion measured by comparing the astrometric positions of each object between the two datasets over about 14 yr time baseline. We refine the WD candidates by fitting blackbody and template WD atmosphere models to HSC photometries for each candidate, enabling the estimation of photometric distance and tangential velocity ($v_{\rm t}$) with respect to the Sun. The deep HSC data allows us to identify low-temperature ($<4000$ K) and faint WD candidates down to absolute magnitude, $M_{\rm bol}\simeq 17$. We evaluate the selection function of our WD candidates using the mock catalogue for spatial and kinematic distributions of WDs in the (thin and thick) disc and halo regions based on the standard Milky Way model. We construct the samples of disc and halo WD candidates by selecting WDs with the cuts of tangential velocity, $40<v_{\rm t}/[{\rm km}~{\rm s}^{-1}]<80$ and $200<v_{\rm t}/[{\rm km}~{\rm s}^{-1}]<500$, respectively. The total number densities of the disc and halo WDs are $(9.45 \pm 0.94) \times 10^{-3}$ pc$^{-3}$ and $(4.20 \pm 1.74) \times 10^{-4}$ pc$^{-3}$, respectively. Our LFs extend down to fainter absolute magnitudes compared with previous work. The faint WDs could represent the oldest generation of building blocks in the tens of billions of years of our Milky Way's assembly history. △ Less

Submitted 31 January, 2024; originally announced January 2024.

Comments: 18 pages, 18 figures, submitted to MNRAS

Prospects for constraining quasar ages with fiber spectrographs: quasar-induced Ly$α$ emission from the intergalactic medium

Authors: Ryuichiro Hada, Masahiro Takada, Akio K. Inoue

Abstract: We present a theoretical framework for linking quasar properties, such as quasar age, to the surrounding Ly$α$ emission intensity. In particular, we focus on a method for mapping the large-scale structure of Ly$α$ emission intensity with galaxy spectra from wide-field spectroscopic surveys, e.g., the Subaru Prime Focus Spectrograph (PFS) or the Dark Energy Spectroscopic Instrument (DESI), and cons… ▽ More We present a theoretical framework for linking quasar properties, such as quasar age, to the surrounding Ly$α$ emission intensity. In particular, we focus on a method for mapping the large-scale structure of Ly$α$ emission intensity with galaxy spectra from wide-field spectroscopic surveys, e.g., the Subaru Prime Focus Spectrograph (PFS) or the Dark Energy Spectroscopic Instrument (DESI), and consider the quasar-induced Ly$α$ emission from the intergalactic medium (IGM). To do this, we construct a theoretical model based on two physical processes: resonant scattering of quasar Ly$α$ photons and fluorescence due to quasar ionizing photons, finding that the fluorescence contribution due to optically thick gas clouds is dominant. Taking into account the light cone effect and assuming a typical quasar spectrum, we calculate the fluorescence contribution to the spectrum stacked within each bin of the separation angle from the quasar as a function of quasar age. Furthermore, we compute the quasar-Ly$α$ emission cross-correlation and its SNR for the planned PFS survey. The predicted signal can only account for $\sim10\%$ of the measurements indicated from the BOSS and eBOSS surveys in the outer region of $> 10\ \rm{cMpc}\ \rm{h}^{-1}$, and the predicted SNR is not sufficient to detect the quasar-induced contribution. However, we found that our model, combined with the contribution of star-forming galaxies, is not in conflict with these measurements. We discuss other possible contributions to the Ly$α$ emission excess around quasars, the efficiency of using spectroscopic fibers, and the redshift dependence of our model. △ Less

Submitted 9 June, 2024; v1 submitted 4 December, 2023; originally announced December 2023.

Comments: 20 pages, 12 figures, accepted for publication in MNRAS; moderate changes to match published version

The Milky Way Tomography with Subaru Hyper Suprime-Cam. I. Halo substructures

Authors: Yoshihisa Suzuki, Masashi Chiba, Yutaka Komiyama, Kohei Hayashi, Masayuki Tanaka, Tetsuya Fukushima, Scott G. Carlsten, Akira Tokiwa, Tian Qiu, Masahiro Takada

Abstract: We analyze the photometric data in the Wide layer of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) over $\sim 1,200$ deg$^{2}$ to uncover new halo substructures beyond the distance, $D_{\odot}\sim$ 30 kpc, from the Sun. For this purpose, we develop an isochrone filter for an old, metal-poor stellar system to extract the faint main-sequence stars at a range of distances. With this method… ▽ More We analyze the photometric data in the Wide layer of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) over $\sim 1,200$ deg$^{2}$ to uncover new halo substructures beyond the distance, $D_{\odot}\sim$ 30 kpc, from the Sun. For this purpose, we develop an isochrone filter for an old, metal-poor stellar system to extract the faint main-sequence stars at a range of distances. With this method, we detect, not only the previously discovered substructures such as the Orphan Stream, but also the new overdensity toward Boötes at about $D_{\odot}\sim$ 60 kpc and the new stream-like feature toward Pisces at around $D_{\odot}\sim$ 60 kpc. It has been suggested that a small-scale overdensity exists in this direction of Pisces (the so-called Pisces Overdensity), but our results show that the overdensity is widely spread with a tidally elongated feature. Combining our results with the ongoing Hyper Suprime-Cam narrow-band survey and the near-future spectroscopic survey with Prime Focus Spectrograph (PFS) will allow us to place strong constraints on the origin of these halo substructures. △ Less

Submitted 10 January, 2024; v1 submitted 29 November, 2023; originally announced November 2023.

Comments: 16 pages, 11 figures, 1 table, accepted for publication in PASJ

Flat to nonflat: Calculating nonlinear power spectra of biased tracers for nonflat $Λ$CDM model

Authors: Ryo Terasawa, Ryuichi Takahashi, Takahiro Nishimichi, Masahiro Takada

Abstract: The growth of large-scale structure, together with the geometrical information of cosmic expansion history and cosmological distances, can be used to obtain constraints on the spatial curvature of the universe that probes the early universe physics, whereas modeling the nonlinear growth in a nonflat universe is still challenging due to computational expense of simulations in a high-dimensional cos… ▽ More The growth of large-scale structure, together with the geometrical information of cosmic expansion history and cosmological distances, can be used to obtain constraints on the spatial curvature of the universe that probes the early universe physics, whereas modeling the nonlinear growth in a nonflat universe is still challenging due to computational expense of simulations in a high-dimensional cosmological parameter space. In this paper, we develop an approximate method to compute the halo-matter and halo-auto power spectra for nonflat $Λ$CDM model, from quantities representing the nonlinear evolution of the corresponding flat $Λ$CDM model, based on the separate universe (SU) method. By utilizing the fact that the growth response to long-wavelength fluctuations (equivalently the curvature), $T_{δ_{\rm b}}(k)$, is approximated by the response to the Hubble parameter, $T_h(k)$, our method allows one to estimate the nonlinear power spectra in a nonflat universe efficiently from the power spectra of the flat universe. We use $N$-body simulations to show that the estimator can provide the halo-matter (halo-auto) power spectrum at $\sim 1\%$ ($\sim 2\%$ ) accuracy up to $k \simeq 3 (1) \, h {\rm Mpc}^{-1}$ even for a model with large curvature $Ω_K = \pm 0.1$. Using the estimator we can extend the prediction of the existing emulators such as Dark Emulator to nonflat models without degrading their accuracy. Since the response to long-wavelength fluctuations is also a key quantity for estimating the super sample covariance (SSC), we discuss that the approximate identity $T_{δ_{\rm b}}(k) \approx T_h(k)$ can be used to calculate the SSC terms analytically. △ Less

Submitted 20 October, 2023; originally announced October 2023.

Comments: 18 pages, 7 figures. arXiv admin note: text overlap with arXiv:2205.10339

The Intrinsic Alignment of Galaxy Clusters and Impact of Projection Effects

Authors: Jingjing Shi, Tomomi Sunayama, Toshiki Kurita, Masahiro Takada, Sunao Sugiyama, Rachel Mandelbaum, Hironao Miyatake, Surhud More, Takahiro Nishimichi, Harry Johnston

Abstract: Galaxy clusters, being the most massive objects in the Universe, exhibit the strongest alignment with the large-scale structure. However, mis-identification of members due to projection effects from the large scale structure can occur. We studied the impact of projection effects on the measurement of the intrinsic alignment of galaxy clusters, using galaxy cluster mock catalogs. Our findings showe… ▽ More Galaxy clusters, being the most massive objects in the Universe, exhibit the strongest alignment with the large-scale structure. However, mis-identification of members due to projection effects from the large scale structure can occur. We studied the impact of projection effects on the measurement of the intrinsic alignment of galaxy clusters, using galaxy cluster mock catalogs. Our findings showed that projection effects result in a decrease of the large scale intrinsic alignment signal of the cluster and produce a bump at $r_p\sim 1h^{-1}/Mpc$, most likely due to interlopers and missed member galaxies. This decrease in signal explains the observed similar alignment strength between bright central galaxies and clusters in the SDSS redMaPPer cluster catalog. The projection effect and cluster intrinsic alignment signal are coupled, with clusters having lower fractions of missing members or having higher fraction of interlopers exhibiting higher alignment signals in their projected shapes. We aim to use these findings to determine the impact of projection effects on galaxy cluster cosmology in future studies. △ Less

Submitted 10 January, 2024; v1 submitted 16 June, 2023; originally announced June 2023.

Comments: 15 pages, 14 figures, Accepted by MNRAS

Quasar Luminosity Function at z = 7

Authors: Yoshiki Matsuoka, Masafusa Onoue, Kazushi Iwasawa, Michael A. Strauss, Nobunari Kashikawa, Takuma Izumi, Tohru Nagao, Masatoshi Imanishi, Masayuki Akiyama, John D. Silverman, Naoko Asami, James Bosch, Hisanori Furusawa, Tomotsugu Goto, James E. Gunn, Yuichi Harikane, Hiroyuki Ikeda, Kohei Inayoshi, Rikako Ishimoto, Toshihiro Kawaguchi, Satoshi Kikuta, Kotaro Kohno, Yutaka Komiyama, Chien-Hsiu Lee, Robert H. Lupton , et al. (19 additional authors not shown)

Abstract: We present the quasar luminosity function (LF) at $z = 7$, measured with 35 spectroscopically confirmed quasars at $6.55 < z < 7.15$. The sample of 22 quasars from the Subaru High-$z$ Exploration of Low-Luminosity Quasars (SHELLQs) project, combined with 13 brighter quasars in the literature, covers an unprecedentedly wide range of rest-frame ultraviolet magnitudes over $-28 < M_{1450} < -23$. We… ▽ More We present the quasar luminosity function (LF) at $z = 7$, measured with 35 spectroscopically confirmed quasars at $6.55 < z < 7.15$. The sample of 22 quasars from the Subaru High-$z$ Exploration of Low-Luminosity Quasars (SHELLQs) project, combined with 13 brighter quasars in the literature, covers an unprecedentedly wide range of rest-frame ultraviolet magnitudes over $-28 < M_{1450} < -23$. We found that the binned LF flattens significantly toward the faint end populated by the SHELLQs quasars. A maximum likelihood fit to a double power-law model has a break magnitude $M^*_{1450} = -25.60^{+0.40}_{-0.30}$, a characteristic density $Φ^* = 1.35^{+0.47}_{-0.30}$ Gpc$^{-3}$ mag$^{-1}$, and a bright-end slope $β= -3.34^{+0.49}_{-0.57}$, when the faint-end slope is fixed to $α= -1.2$ as observed at $z \le 6$. The overall LF shape remains remarkably similar from $z = 4$ to $7$, while the amplitude decreases substantially toward higher redshifts, with a clear indication of an accelerating decline at $z \ge 6$. The estimated ionizing photon density, $10^{48.2 \pm 0.1}$ s$^{-1}$ Mpc$^{-3}$, is less than 1 % of the critical rate to keep the intergalactic medium ionized at $z = 7$, and thus indicates that quasars are not a major contributor to cosmic reionization. △ Less

Submitted 18 May, 2023; originally announced May 2023.

Comments: The Astrophysical Journal Letters, in press

Hyper Suprime-Cam Year 3 Results: Cosmology from Galaxy Clustering and Weak Lensing with HSC and SDSS using the Minimal Bias Model

Authors: Sunao Sugiyama, Hironao Miyatake, Surhud More, Xiangchong Li, Masato Shirasaki, Masahiro Takada, Yosuke Kobayashi, Ryuichi Takahashi, Takahiro Nishimichi, Atsushi J. Nishizawa, Markus M. Rau, Tianqing Zhang, Roohi Dalal, Rachel Mandelbaum, Michael A. Strauss, Takashi Hamana, Masamune Oguri, Ken Osato, Arun Kannawadi, Robert Armstrong, Yutaka Komiyama, Robert H. Lupton, Nate B. Lust, Satoshi Miyazaki, Hitoshi Murayama , et al. (5 additional authors not shown)

Abstract: We present cosmological parameter constraints from a blind joint analysis of three two-point correlation functions measured from the Year 3 Hyper Suprime-Cam (HSC-Y3) imaging data, covering 416 deg$^2$, and the SDSS DR11 spectroscopic galaxies spanning the redshift range $[0.15, 0.70]$. We subdivide the SDSS galaxies into three volume-limited samples separated in redshift, each of which acts as a… ▽ More We present cosmological parameter constraints from a blind joint analysis of three two-point correlation functions measured from the Year 3 Hyper Suprime-Cam (HSC-Y3) imaging data, covering 416 deg$^2$, and the SDSS DR11 spectroscopic galaxies spanning the redshift range $[0.15, 0.70]$. We subdivide the SDSS galaxies into three volume-limited samples separated in redshift, each of which acts as a large-scale structure tracer characterized by the measurement of the projected correlation function, $w_{\rm p}(R)$. We also use the measurements of the galaxy-galaxy weak lensing signal $ΔΣ(R)$ for each of these SDSS samples which act as lenses for a secure sample of source galaxies selected from the HSC-Y3 shape catalog based on their photometric redshifts. We combine these measurements with the cosmic shear correlation functions, $ξ_{\pm}(\vartheta)$, measured for our HSC source sample. We model these observables with the minimal bias model of the galaxy clustering observables in the context of a flat $Λ$CDM cosmology. We use conservative scale cuts, $R>12$ and $8~h^{-1}$Mpc, for $ΔΣ$ and $w_{\rm p}$, respectively, where the minimal bias model is valid, in addition to conservative prior on the residual bias in the mean redshift of the HSC photometric source galaxies. Our baseline analysis yields $S_8=0.775^{+0.043}_{-0.038}$ (68% C.I.) for the $Λ$CDM model, after marginalizing over uncertainties in other parameters. Our value of $S_8$ is consistent with that from the Planck 2018 data, but the credible interval of our result is still relatively large. Our results are statistically consistent with those of a companion paper, which extends this analysis to smaller scales with an emulator-based halo model. △ Less

Submitted 27 December, 2023; v1 submitted 2 April, 2023; originally announced April 2023.

Comments: 29 pages, 17 figures, 4 tables, version accepted for publication in PRD together with other HSC Y3 weak lensing cosmology papers - see https://hsc-release.mtk.nao.ac.jp/doc/index.php/wly3/

Journal ref: Phys. Rev. D 108, 123521 (2023)

Hyper Suprime-Cam Year 3 Results: Cosmology from Galaxy Clustering and Weak Lensing with HSC and SDSS using the Emulator Based Halo Model

Authors: Hironao Miyatake, Sunao Sugiyama, Masahiro Takada, Takahiro Nishimichi, Xiangchong Li, Masato Shirasaki, Surhud More, Yosuke Kobayashi, Atsushi J. Nishizawa, Markus M. Rau, Tianqing Zhang, Ryuichi Takahashi, Roohi Dalal, Rachel Mandelbaum, Michael A. Strauss, Takashi Hamana, Masamune Oguri, Ken Osato, Wentao Luo, Arun Kannawadi, Bau-Ching Hsieh, Robert Armstrong, Yutaka Komiyama, Robert H. Lupton, Nate B. Lust , et al. (9 additional authors not shown)

Abstract: We present cosmology results from a blinded joint analysis of cosmic shear, $ξ_{\pm}(\vartheta)$, galaxy-galaxy weak lensing, $Δ\!Σ(R)$, and projected galaxy clustering, $w_{\rm p}(R)$, measured from the Hyper Suprime-Cam three-year (HSC-Y3) shape catalog and the Sloan Digital Sky Survey (SDSS) DR11 spectroscopic galaxy catalog - a 3$\times$2pt cosmology analysis. We define luminosity-cut samples… ▽ More We present cosmology results from a blinded joint analysis of cosmic shear, $ξ_{\pm}(\vartheta)$, galaxy-galaxy weak lensing, $Δ\!Σ(R)$, and projected galaxy clustering, $w_{\rm p}(R)$, measured from the Hyper Suprime-Cam three-year (HSC-Y3) shape catalog and the Sloan Digital Sky Survey (SDSS) DR11 spectroscopic galaxy catalog - a 3$\times$2pt cosmology analysis. We define luminosity-cut samples of SDSS galaxies to serve as the tracers of $w_{\rm p}$ and as the lens samples for $Δ\!Σ$ in three spectroscopic redshift bins spanning the range $0.15<z<0.7$. For the $ξ_{\pm}$ and $Δ\!Σ$ measurements, we use a single source sample over 416 deg$^2$, selected from HSC-Y3 based on having photometric redshifts (photo-$z$) greater than 0.75. For cosmological parameter inference, we use Dark Emulator combined with a halo occupation distribution prescription to model $w_{\rm p}$ and $Δ\!Σ$ down to quasi-nonlinear scales. In our baseline analysis we employ an uninformative flat prior of the residual photo-$z$ error to model a residual bias in the mean redshift of HSC source galaxies. We obtain a robust constraint on the cosmological parameters for the flat $Λ$CDM model: $S_8=σ_8(Ω_{\rm m}/0.3)^{0.5}=0.763^{+0.040}_{-0.036}$ (68% C.I.), or the best-constrained parameter given by $S'_8=σ_8(Ω_{\rm m}/0.3)^{0.22}=0.721\pm 0.028$, determined with about 4% fractional precision. Our HSC-Y3 data exhibits about 2.5$σ$ tension with the Planck inferred $S_8$ value for the $Λ$CDM model, and hints at a non-zero residual photo-$z$ bias implying that the true mean redshift of the HSC galaxies at $z\gtrsim 0.75$ is higher than that implied by the original photo-$z$ estimates. △ Less

Submitted 6 April, 2023; v1 submitted 2 April, 2023; originally announced April 2023.

Comments: 39 pages, 25 figures, 5 tables, for coordinated submission to PRD with other HSC Y3 weak lensing cosmology papers - see https://hsc-release.mtk.nao.ac.jp/doc/index.php/wly3/. arXiv admin note: text overlap with arXiv:2111.02419

Hyper Suprime-Cam Year 3 Results: Measurements of Clustering of SDSS-BOSS Galaxies, Galaxy-Galaxy Lensing and Cosmic Shear

Authors: Surhud More, Sunao Sugiyama, Hironao Miyatake, Markus Michael Rau, Masato Shirasaki, Xiangchong Li, Atsushi J. Nishizawa, Ken Osato, Tianqing Zhang, Masahiro Takada, Takashi Hamana, Ryuichi Takahashi, Roohi Dalal, Rachel Mandelbaum, Michael A. Strauss, Yosuke Kobayashi, Takahiro Nishimichi, Masamune Oguri, Wentao Luo, Arun Kannawadi, Bau-Ching Hsieh, Robert Armstrong, James Bosch, Yutaka Komiyama, Robert H. Lupton , et al. (9 additional authors not shown)

Abstract: We use the Sloan Digital Sky Survey (SDSS) BOSS galaxies and their overlap with approximately 416 sq. degree of deep $grizy$-band imaging from the Subaru Hyper Suprime-Cam Survey (HSC). We measure three two-point correlations that form the basis of the cosmological inference presented in our companion papers, Miyatake et al. and Sugiyama et al. We use three approximately volume limited subsamples… ▽ More We use the Sloan Digital Sky Survey (SDSS) BOSS galaxies and their overlap with approximately 416 sq. degree of deep $grizy$-band imaging from the Subaru Hyper Suprime-Cam Survey (HSC). We measure three two-point correlations that form the basis of the cosmological inference presented in our companion papers, Miyatake et al. and Sugiyama et al. We use three approximately volume limited subsamples of spectroscopic galaxies by their $i$-band magnitude from the SDSS-BOSS: LOWZ (0.1<z<0.35), CMASS1 (0.43<z<0.55) and CMASS2 (0.55<z<0.7), respectively. We present high signal-to-noise ratio measurements of the projected correlation functions of these galaxies, which is expected to be proportional to the matter correlation function times the bias of galaxies on large scales. In order to break the degeneracy between the amplitude of the matter correlation and the bias of these galaxies, we use the distortions of the shapes of galaxies in HSC due to weak gravitational lensing, to measure the galaxy-galaxy lensing signal, which probes the galaxy-matter cross-correlation of the SDSS-BOSS galaxies. We also measure the cosmic shear correlation functions from HSC galaxies which is related to the projected matter correlation function. We demonstrate the robustness of our measurements with a variety of systematic tests. Our use of a single sample of HSC source galaxies is crucial to calibrate any residual systematic biases in the inferred redshifts of our galaxies. We also describe the construction of a suite of mocks: i) spectroscopic galaxy catalogs which obey the clustering and abundance of each of the three SDSS-BOSS subsamples, and ii) galaxy shape catalogs which obey the footprint of the HSC survey and have been appropriately sheared by the large-scale structure expected in a $Λ$-CDM model. We use these mock catalogs to compute the covariance of each of our observables. △ Less

Submitted 16 November, 2023; v1 submitted 2 April, 2023; originally announced April 2023.

Comments: 31 pages, 24 figures, version accepted for publication in PRD together with other HSC Y3 weak lensing cosmology papers - see https://hsc-release.mtk.nao.ac.jp/doc/index.php/wly3/

Hyper Suprime-Cam Year 3 Results: Cosmology from Cosmic Shear Two-point Correlation Functions

Authors: Xiangchong Li, Tianqing Zhang, Sunao Sugiyama, Roohi Dalal, Ryo Terasawa, Markus M. Rau, Rachel Mandelbaum, Masahiro Takada, Surhud More, Michael A. Strauss, Hironao Miyatake, Masato Shirasaki, Takashi Hamana, Masamune Oguri, Wentao Luo, Atsushi J. Nishizawa, Ryuichi Takahashi, Andrina Nicola, Ken Osato, Arun Kannawadi, Tomomi Sunayama, Robert Armstrong, James Bosch, Yutaka Komiyama, Robert H. Lupton , et al. (10 additional authors not shown)

Abstract: We perform a blinded cosmology analysis with cosmic shear two-point correlation functions (2PCFs) 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 deg$^2$ of the northern sky with an effective galaxy number density of 15 arcm… ▽ More We perform a blinded cosmology analysis with cosmic shear two-point correlation functions (2PCFs) 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 deg$^2$ 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 < θ/{\rm arcmin} < 56.6$ for $ξ_+$ and $31.2 <θ/{\rm 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 $Ω_{\rm m}=0.256_{-0.044}^{+0.056}$ and $S_8\equiv σ_8 \sqrt{Ω_{\rm m}/0.3}=0.769_{-0.034}^{+0.031}$ (both 68\% CI) for a flat $Λ$ cold dark matter cosmology. We find, after unblinding, that our constraint on $S_8$ is consistent with the Fourier space cosmic shear and the 3$\times$2pt 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 $S_8$. 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. △ Less

Submitted 30 November, 2023; v1 submitted 2 April, 2023; originally announced April 2023.

Comments: 38 pages, 32 figures, 4 tables (PRD in press.)

Hyper Suprime-Cam Year 3 Results: Cosmology from Cosmic Shear Power Spectra

Authors: Roohi Dalal, Xiangchong Li, Andrina Nicola, Joe Zuntz, Michael A. Strauss, Sunao Sugiyama, Tianqing Zhang, Markus M. Rau, Rachel Mandelbaum, Masahiro Takada, Surhud More, Hironao Miyatake, Arun Kannawadi, Masato Shirasaki, Takanori Taniguchi, Ryuichi Takahashi, Ken Osato, Takashi Hamana, Masamune Oguri, Atsushi J. Nishizawa, Andrés A. Plazas Malagón, Tomomi Sunayama, David Alonso, Anže Slosar, Robert Armstrong , et al. (13 additional authors not shown)

Abstract: We measure weak lensing cosmic shear power spectra from the three-year galaxy shear catalog of the Hyper Suprime-Cam (HSC) Subaru Strategic Program imaging survey. The shear catalog covers $416 \ \mathrm{deg}^2$ of the northern sky, with a mean $i$-band seeing of 0.59 arcsec and an effective galaxy number density of 15 $\mathrm{arcmin}^{-2}$ within our adopted redshift range. With an $i$-band magn… ▽ More We measure weak lensing cosmic shear power spectra from the three-year galaxy shear catalog of the Hyper Suprime-Cam (HSC) Subaru Strategic Program imaging survey. The shear catalog covers $416 \ \mathrm{deg}^2$ of the northern sky, with a mean $i$-band seeing of 0.59 arcsec and an effective galaxy number density of 15 $\mathrm{arcmin}^{-2}$ within our adopted redshift range. With an $i$-band magnitude limit of 24.5 mag, and four tomographic redshift bins spanning $0.3 \leq z_{\mathrm{ph}} \leq 1.5$ based on photometric redshifts, we obtain a high-significance measurement of the cosmic shear power spectra, with a signal-to-noise ratio of approximately 26.4 in the multipole range $300<\ell<1800$. The accuracy of our power spectrum measurement is tested against realistic mock shear catalogs, and we use these catalogs to get a reliable measurement of the covariance of the power spectrum measurements. We use a robust blinding procedure to avoid confirmation bias, and model various uncertainties and sources of bias in our analysis, including point spread function systematics, redshift distribution uncertainties, the intrinsic alignment of galaxies and the modeling of the matter power spectrum. For a flat $Λ$CDM model, we find $S_8 \equiv σ_8 (Ω_m/0.3)^{0.5} =0.776^{+0.032}_{-0.033}$, which is in excellent agreement with the constraints from the other HSC Year 3 cosmology analyses, as well as those from a number of other cosmic shear experiments. This result implies a $\sim$$2σ$-level tension with the Planck 2018 cosmology. We study the effect that various systematic errors and modeling choices could have on this value, and find that they can shift the best-fit value of $S_8$ by no more than $\sim$$0.5σ$, indicating that our result is robust to such systematics. △ Less

Submitted 4 April, 2023; v1 submitted 2 April, 2023; originally announced April 2023.

Comments: 35 pages, 18 figures, 6 tables, for coordinated submission to PRD with other HSC Y3 weak lensing cosmology papers - see https://hsc-release.mtk.nao.ac.jp/doc/index.php/wly3/

Journal ref: Physical Review D, Volume 108, Issue 12, December 2023, article id.123519

Study of structural parameters and systemic proper motion of Sextans dwarf spheroidal galaxy with Subaru Hyper Suprime-Cam data

Authors: Akira Tokiwa, Masahiro Takada, Tian Qiu, Naoki Yasuda, Yutaka Komiyama, Masashi Chiba, Kohei Hayashi

Abstract: We use the Subaru Hyper Suprime-Cam (HSC) data to study structural parameters and systemic proper motion of the Sextans dwarf spheroidal galaxy at the heliocentric distance of 86 kpc, which is one of the most important targets for studies of dark matter nature and galaxy formation physics. Thanks to the superb image quality and wide area coverage of the Sextans field, the HSC data enables a secure… ▽ More We use the Subaru Hyper Suprime-Cam (HSC) data to study structural parameters and systemic proper motion of the Sextans dwarf spheroidal galaxy at the heliocentric distance of 86 kpc, which is one of the most important targets for studies of dark matter nature and galaxy formation physics. Thanks to the superb image quality and wide area coverage of the Sextans field, the HSC data enables a secure selection of member star candidates based on the colour-magnitude cut, yielding about 10,000 member candidates at magnitudes down to $i\sim 24$. We use a likelihood analysis of the two-dimensional distribution of stars to estimate the structural parameters of Sextans taking into account the contamination of foreground halo stars in the Milky Way, and find that the member star distribution is well-fitted by an elliptical King profile with ellipticity $ε\simeq 0.25$ and the core and tidal radii of $R_c=(368.4\pm 8.5)$ pc and $R_t=(2.54\pm 0.046)$ kpc, respectively. Then using the two HSC datasets of 2.66 years time baseline on average, we find the systemic proper motions of Sextans to be $(μ_α, μ_δ)=(-0.448\pm{0.075},0.083\pm{0.078})$ $\mathrm{mas}\ \mathrm{yr}^{-1}$, which is consistent with some of the previous works using the $Gaia$ data of relatively bright member stars in Sextans. Thus, our results give a demonstration that a ground-based, large-aperture telescope data which covers a wide solid angle of the sky and has a long time baseline, such as the upcoming LSST data, can be used to study systemic proper motions of dwarf galaxies. △ Less

Submitted 13 July, 2023; v1 submitted 12 February, 2023; originally announced February 2023.

Comments: 14 pages, 15 figures

Constraints on anisotropic primordial non-Gaussianity from intrinsic alignments of SDSS-III BOSS galaxies

Authors: Toshiki Kurita, Masahiro Takada

Abstract: We measure the three-dimensional cross-power spectrum of galaxy density and intrinsic alignment (IA) fields for the first time from the spectroscopic and imaging data of SDSS-III BOSS galaxies, for each of the four samples in the redshift range $0.2 < z < 0.75$. In the measurement we use the power spectrum estimator, developed in our previous work, to take into account the line-of-sight dependent… ▽ More We measure the three-dimensional cross-power spectrum of galaxy density and intrinsic alignment (IA) fields for the first time from the spectroscopic and imaging data of SDSS-III BOSS galaxies, for each of the four samples in the redshift range $0.2 < z < 0.75$. In the measurement we use the power spectrum estimator, developed in our previous work, to take into account the line-of-sight dependent projection of galaxy shapes onto the sky coordinate and the $E/B$-mode decomposition of the spin-2 shape field. Our method achieves a significant detection of the $E$-mode power spectrum with the total signal-to-noise ratio comparable with that of the quadrupole moment of the galaxy density power spectrum, while the measured $B$-mode power spectra are consistent with a null signal to within the statistical errors for all the galaxy samples. We also show that, compared to the previous results based on the two-dimensional projected correlation function, our method improves the precision of the linear shape bias parameter estimation by up to a factor of two thanks to the three-dimensional information. By performing a joint analysis of the galaxy density and IA power spectra in the linear regime, we constrain the isotropic and anisotropic local primordial non-Gaussianities (PNGs) parameters, $f_\mathrm{NL}^{s=0}$ and $f_\mathrm{NL}^{s=2}$, simultaneously, where the two types of PNGs induce characteristic scale-dependent biases at very large scales in the density and IA power spectra, respectively. We do not find any significant detection for both PNGs: the constraints $f^{s=0}_\mathrm{NL}=57^{+30}_{-29}$ and $f^{s=2}_\mathrm{NL} = -67_{-269}^{+285}$ ($68\%$ C.L.), respectively. Our method paves the way for using the IA power spectrum as a cosmological probe for current and future galaxy surveys. △ Less

Submitted 6 February, 2023; originally announced February 2023.

Comments: 46 pages, 15 figures

Constraining cosmological vorticity modes with CMB secondary anisotropies

Authors: William R. Coulton, Kazuyuki Akitsu, Masahiro Takada

Abstract: Observational searches for large-scale vorticity modes in the late time Universe are underexplored. Within the standard $Λ$CDM model, this is well motivated given the observed properties of the cosmic microwave background (CMB). However, this means that searches for cosmic vorticity modes can serve as a powerful consistency test of our cosmological model. We show that through combining CMB measure… ▽ More Observational searches for large-scale vorticity modes in the late time Universe are underexplored. Within the standard $Λ$CDM model, this is well motivated given the observed properties of the cosmic microwave background (CMB). However, this means that searches for cosmic vorticity modes can serve as a powerful consistency test of our cosmological model. We show that through combining CMB measurements of the kinetic Sunyaev-Zel'dovich and the moving lens effects with galaxy survey data we can constrain vorticity fields independently from the large scale cosmic velocity field. This approach can provide stringent constraints on the largest scale modes and can be achieved by a simple change in the standard estimators. Alternatively if one assumes there are no cosmic vorticity modes, this estimator can be used to test for systematic biases in existing analyses of kinetic Sunyaev-Zel'dovich effect in a manner analogous to curl-lensing. △ Less

Submitted 26 January, 2023; originally announced January 2023.

Comments: 6 pages, 1 figure

A General Framework for Removing Point Spread Function Additive Systematics in Cosmological Weak Lensing Analysis

Authors: Tianqing Zhang, Xiangchong Li, Roohi Dalal, Rachel Mandelbaum, Michael A. Strauss, Arun Kannawadi, Hironao Miyatake, Andrina Nicola, Andrés A. Plazas Malagón, Masato Shirasaki, Sunao Sugiyama, Masahiro Takada, Surhud More

Abstract: Cosmological weak lensing measurements rely on a precise measurement of the shear two-point correlation function (2PCF) along with a deep understanding of systematics that affect it. In this work, we demonstrate a general framework for detecting and modeling the impact of PSF systematics on the cosmic shear 2PCF, and mitigating its impact on cosmological analysis. Our framework can describe leakag… ▽ More Cosmological weak lensing measurements rely on a precise measurement of the shear two-point correlation function (2PCF) along with a deep understanding of systematics that affect it. In this work, we demonstrate a general framework for detecting and modeling the impact of PSF systematics on the cosmic shear 2PCF, and mitigating its impact on cosmological analysis. Our framework can describe leakage and modeling error from all spin-2 quantities contributed by the PSF second and higher moments, rather than just the second moments, using the cross-correlations between galaxy shapes and PSF moments. We interpret null tests using the HSC Year 3 (Y3) catalogs with this formalism, and find that leakage from the spin-2 combination of PSF fourth moments is the leading contributor to additive shear systematics, with total contamination that is an order of magnitude higher than that contributed by PSF second moments alone. We conducted a mock cosmic shear analysis for HSC Y3, and find that, if uncorrected, PSF systematics can bias the cosmological parameters $Ω_m$ and $S_8$ by $\sim$0.3$σ$. The traditional second moment-based model can only correct for a 0.1$σ$ bias, leaving the contamination largely uncorrected. We conclude it is necessary to model both PSF second and fourth moment contamination for HSC Y3 cosmic shear analysis. We also reanalyze the HSC Y1 cosmic shear analysis with our updated systematics model, and identify a 0.07$σ$ bias on $Ω_m$ when using the more restricted second moment model from the original analysis. We demonstrate how to self-consistently use the method in both real space and Fourier space, assess shear systematics in tomographic bins, and test for PSF model overfitting. △ Less

Submitted 15 November, 2023; v1 submitted 6 December, 2022; originally announced December 2022.

Comments: 33 pages, 28 figures, match the refereed version in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, 525, 2441-2471, 2023

Weak Lensing Tomographic Redshift Distribution Inference for the Hyper Suprime-Cam Subaru Strategic Program three-year shape catalogue

Authors: Markus Michael Rau, Roohi Dalal, Tianqing Zhang, Xiangchong Li, Atsushi J. Nishizawa, Surhud More, Rachel Mandelbaum, Hironao Miyatake, Michael A. Strauss, Masahiro Takada

Abstract: We present posterior sample redshift distributions for the Hyper Suprime-Cam Subaru Strategic Program Weak Lensing three-year (HSC Y3) analysis. Using the galaxies' photometry and spatial cross-correlations, we conduct a combined Bayesian Hierarchical Inference of the sample redshift distributions. The spatial cross-correlations are derived using a subsample of Luminous Red Galaxies (LRGs) with ac… ▽ More We present posterior sample redshift distributions for the Hyper Suprime-Cam Subaru Strategic Program Weak Lensing three-year (HSC Y3) analysis. Using the galaxies' photometry and spatial cross-correlations, we conduct a combined Bayesian Hierarchical Inference of the sample redshift distributions. The spatial cross-correlations are derived using a subsample of Luminous Red Galaxies (LRGs) with accurate redshift information available up to a photometric redshift of $z < 1.2$. We derive the photometry-based constraints using a combination of two empirical techniques calibrated on spectroscopic- and multiband photometric data that covers a spatial subset of the shear catalog. The limited spatial coverage induces a cosmic variance error budget that we include in the inference. Our cross-correlation analysis models the photometric redshift error of the LRGs to correct for systematic biases and statistical uncertainties. We demonstrate consistency between the sample redshift distributions derived using the spatial cross-correlations, the photometry, and the posterior of the combined analysis. Based on this assessment, we recommend conservative priors for sample redshift distributions of tomographic bins used in the three-year cosmological Weak Lensing analyses. △ Less

Submitted 20 November, 2023; v1 submitted 29 November, 2022; originally announced November 2022.

Comments: 23 pages, 11 figures, 1 table, updated to match the accepted version in the MNRAS

Galaxy clustering from the bottom up: A Streaming Model emulator I

Authors: Carolina Cuesta-Lazaro, Takahiro Nishimichi, Yosuke Kobayashi, Cheng-Zong Ruan, Alexander Eggemeier, Hironao Miyatake, Masahiro Takada, Naoki Yoshida, Pauline Zarrouk, Carlton M. Baugh, Sownak Bose, Baojiu Li

Abstract: In this series of papers, we present a simulation-based model for the non-linear clustering of galaxies based on separate modelling of clustering in real space and velocity statistics. In the first paper, we present an emulator for the real-space correlation function of galaxies, whereas the emulator of the real-to-redshift space mapping based on velocity statistics is presented in the second pape… ▽ More In this series of papers, we present a simulation-based model for the non-linear clustering of galaxies based on separate modelling of clustering in real space and velocity statistics. In the first paper, we present an emulator for the real-space correlation function of galaxies, whereas the emulator of the real-to-redshift space mapping based on velocity statistics is presented in the second paper. Here, we show that a neural network emulator for real-space galaxy clustering trained on data extracted from the Dark Quest suite of N-body simulations achieves sub-per cent accuracies on scales $1 < r < 30 $ $h^{-1} \,\mathrm{Mpc}$, and better than $3\%$ on scales $r < 1$ $h^{-1}\mathrm{Mpc}$ in predicting the clustering of dark-matter haloes with number density $10^{-3.5}$ $(h^{-1}\mathrm{Mpc})^{-3}$, close to that of SDSS LOWZ-like galaxies. The halo emulator can be combined with a galaxy-halo connection model to predict the galaxy correlation function through the halo model. We demonstrate that we accurately recover the cosmological and galaxy-halo connection parameters when galaxy clustering depends only on the mass of the galaxies' host halos. Furthermore, the constraining power in $σ_8$ increases by about a factor of $2$ when including scales smaller than $5$ $h^{-1} \,\mathrm{Mpc}$. However, when mass is not the only property responsible for galaxy clustering, as observed in hydrodynamical or semi-analytic models of galaxy formation, our emulator gives biased constraints on $σ_8$. This bias disappears when small scales ($r < 10$ $h^{-1}\mathrm{Mpc}$) are excluded from the analysis. This shows that a vanilla halo model could introduce biases into the analysis of future datasets. △ Less

Submitted 10 August, 2022; originally announced August 2022.

Separate universe approach to evaluate nonlinear matter power spectrum for non-flat $Λ$CDM model

Authors: Ryo Terasawa, Ryuichi Takahashi, Takahiro Nishimichi, Masahiro Takada

Abstract: The spatial curvature ($Ω_K$) of the Universe is one of the most fundamental quantities that could give a link to the early universe physics. In this paper we develop an approximate method to compute the nonlinear matter power spectrum, $P(k)$, for "non-flat" $Λ$CDM models using the separate universe (SU) ansatz which states that the effect of the curvature on structure formation is equivalent to… ▽ More The spatial curvature ($Ω_K$) of the Universe is one of the most fundamental quantities that could give a link to the early universe physics. In this paper we develop an approximate method to compute the nonlinear matter power spectrum, $P(k)$, for "non-flat" $Λ$CDM models using the separate universe (SU) ansatz which states that the effect of the curvature on structure formation is equivalent to that of long-wavelength density fluctuation ($δ_{\rm b}$) in a local volume in the "flat" $Λ$CDM model, via the specific mapping between the background cosmological parameters and redshifts in the non-flat and flat models. By utilizing the fact that the normalized response of $P(k)$ to $δ_{\rm b}$ (equivalently $Ω_K$), which describes how the non-zero $Ω_K$ alters $P(k)$ as a function of $k$, is well approximated by the response to the Hubble parameter $h$ within the flat model, our method allows one to generalize the prediction of $P(k)$ for flat cosmologies via fitting formulae or emulators to that for non-flat cosmologies. We use $N$-body simulations for the non-flat $Λ$CDM models with $|Ω_K|\leq 0.1$ to show that our method can predict $P(k)$ for non-flat models up to $k \simeq 6\,h{\rm Mpc}^{-1}$ in the redshift range $z\simeq [0,1.5]$, to the fractional accuracy within $\sim 1$% that roughly corresponds to requirements for weak lensing cosmology with upcoming surveys. We find that the emulators, those built for flat cosmologies such as EuclidEmulator, can predict the non-flat $P(k)$ with least degradation. △ Less

Submitted 3 June, 2022; v1 submitted 20 May, 2022; originally announced May 2022.

Comments: 18 pages, 11 figures

Report number: IPMU22-0027,YITP-22-49

Analysis method for 3D power spectrum of projected tensor field with fast estimator and window convolution modelling: an application to intrinsic alignments

Authors: Toshiki Kurita, Masahiro Takada

Abstract: Rank-2 tensor fields of large-scale structure, e.g. a tensor field inferred from shapes of galaxies, open up a window to directly access 2-scalar, 2-vector and 2-tensor modes, where the scalar fields can be measured independently from the standard density field that is traced by distribution of galaxies. Here we develop an estimator of the multipole moments of power spectra for the three-dimension… ▽ More Rank-2 tensor fields of large-scale structure, e.g. a tensor field inferred from shapes of galaxies, open up a window to directly access 2-scalar, 2-vector and 2-tensor modes, where the scalar fields can be measured independently from the standard density field that is traced by distribution of galaxies. Here we develop an estimator of the multipole moments of power spectra for the three-dimensional tensor field, taking into account the projection onto plane perpendicular to the line-of-sight direction. To do this, we find that a convenient representation of the power spectrum multipoles can be obtained by the use of the associated Legendre polynomials in the form which allows for the fast Fourier transform estimations under the local plane-parallel (LPP) approximation. The formulation also allows us to obtain the Hankel transforms to connect the two-point statistics in Fourier and configuration space, which are needed to derive theoretical templates of the power spectrum including convolution of a survey window. To validate our estimators, we use the simulation data of the projected tidal field assuming a survey window that mimics the BOSS-like survey footprint. We show that the LPP estimators fairly well recover the multipole moments that are inferred from the global plane-parallel approximation. We find that the survey window causes a more significant change in the multipole moments of projected tensor power spectrum at $k\lesssim 0.1\,h{\rm Mpc}^{-1}$ from the input power spectrum, than in the density power spectrum. Nevertheless, our method to compute the theory template including the survey window effects successfully reproduces the window-convolved multipole moments measured from the simulations. The analysis method presented here paves the way for a cosmological analysis using three-dimensional tensor-type tracers of large-scale structure for current and future surveys. △ Less

Submitted 10 June, 2022; v1 submitted 23 February, 2022; originally announced February 2022.

Comments: 29 pages, 6 figures; accepted for publication in PRD

Report number: IPMU 22-0001

Constraining Dark matter annihilation with Dark Energy Survey Y3 LSBG sample

Authors: Daiki Hashimoto, Atsushi J. Nishizawa, Masahiro Takada

Abstract: To reveal natures of the dark matter (DM) particles, a gamma-ray signal produced in annihilation processes of DM into the standard model particles has been one of the major probes. The cross-correlation between highly DM dominated structures, such as local dwarf galaxies, and observed photons in the direction of the structures, has been explored and provided stringent constraints on the annihilati… ▽ More To reveal natures of the dark matter (DM) particles, a gamma-ray signal produced in annihilation processes of DM into the standard model particles has been one of the major probes. The cross-correlation between highly DM dominated structures, such as local dwarf galaxies, and observed photons in the direction of the structures, has been explored and provided stringent constraints on the annihilation rate. In our previous work, we have shown that it is sufficient to know the distance distribution of the galaxy sample and individual distance measurement is not required for constraining the annihilation rate. In this work, we apply the method to low-surface brightness galaxies (LSBGs) with unknown individual redshifts provided from the Dark Energy Survey (DES) Year 3 data. With all DES-LSBGs of about 24,000 objects, we find that the upper limits of the cross section for bb channel with 95% C.L. is 3*10^-25 cm^3/s at DM mass of 100 GeV. To be more conservative, we remove about 7000 LSBGs within 1 degree from resolved gamma-ray point sources and then the constraint becomes 30 % weaker than the one with all samples in all DM mass ranges. △ Less

Submitted 2 February, 2022; originally announced February 2022.

Comments: 8 pages, 3 figures, submitted to PRD

Cluster cosmology with anisotropic boosts: Validation of a novel forward modeling analysis and application on SDSS redMaPPer clusters

Authors: Youngsoo Park, Tomomi Sunayama, Masahiro Takada, Yosuke Kobayashi, Hironao Miyatake, Surhud More, Takahiro Nishimichi, Sunao Sugiyama

Abstract: We present a novel analysis for cluster cosmology that fully forward models the abundances, weak lensing, and the clustering of galaxy clusters. Our analysis notably includes an empirical model for the anisotropic boosts impacting the lensing and clustering signals of optical clusters. These boosts arise from a preferential selection of clusters surrounded by anisotropic large scale structure, a c… ▽ More We present a novel analysis for cluster cosmology that fully forward models the abundances, weak lensing, and the clustering of galaxy clusters. Our analysis notably includes an empirical model for the anisotropic boosts impacting the lensing and clustering signals of optical clusters. These boosts arise from a preferential selection of clusters surrounded by anisotropic large scale structure, a consequence of the limited discrimination between line-of-sight interlopers and true cluster members offered by photometric surveys. We validate our analysis via a blind cosmology challenge on mocks, and find that we can obtain tight and unbiased cosmological constraints without informative priors or external calibrations on any of our model parameters. We then apply our analysis on the SDSS redMaPPer clusters, and find results favoring low $Ω_\mathrm{m}$ and high $σ_8$, combining to yield the lensing strength constraint $S_8 = 0.718_{-0.021}^{+0.024}$. We investigate potential drivers behind these results through a series of post-unblinding tests, noting that our results are consistent with existing cluster cosmology constraints but clearly inconsistent with other CMB/LSS based cosmology results. From these tests, we find hints that a suppression in the cluster lensing signal may be driving our results. △ Less

Submitted 1 April, 2023; v1 submitted 16 December, 2021; originally announced December 2021.

Comments: v1: 19 pages, 14 figures. Comments welcome. v2: Added clarifications in text and reflected a minor bug fix for Section 3; version submitted to MNRAS. v3: version accepted and published at MNRAS

Report number: YITP-21-160

Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs). XVI. 69 New Quasars at 5.8 < z < 7.0

Authors: Yoshiki Matsuoka, Kazushi Iwasawa, Masafusa Onoue, Takuma Izumi, Nobunari Kashikawa, Michael A. Strauss, Masatoshi Imanishi, Tohru Nagao, Masayuki Akiyama, John D. Silverman, Naoko Asami, James Bosch, Hisanori Furusawa, Tomotsugu Goto, James E. Gunn, Yuichi Harikane, Hiroyuki Ikeda, Rikako Ishimoto, Toshihiro Kawaguchi, Nanako Kato, Satoshi Kikuta, Kotaro Kohno, Yutaka Komiyama, Chien-Hsiu Lee, Robert H. Lupton , et al. (19 additional authors not shown)

Abstract: We present the spectroscopic discovery of 69 quasars at 5.8 < z < 7.0, drawn from the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP) imaging survey data. This is the 16th publication from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, and completes identification of all but the faintest candidates (i.e., i-band dropouts with zAB < 24 and y-band detections, and z… ▽ More We present the spectroscopic discovery of 69 quasars at 5.8 < z < 7.0, drawn from the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP) imaging survey data. This is the 16th publication from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, and completes identification of all but the faintest candidates (i.e., i-band dropouts with zAB < 24 and y-band detections, and z-band dropouts with yAB < 24) with Bayesian quasar probability Pq > 0.1 in the HSC-SSP third public data release (PDR3). The sample reported here also includes three quasars with Pq < 0.1 at z ~ 6.6, which we selected in an effort to completely cover the reddest point sources with simple color cuts. The number of high-z quasars discovered in SHELLQs has now grown to 162, including 23 type-II quasar candidates. This paper also presents identification of seven galaxies at 5.6 < z < 6.7, an [O III] emitter at z = 0.954, and 31 Galactic cool stars and brown dwarfs. High-z quasars and galaxies comprise 75 % and 16 % respectively of all the spectroscopic SHELLQs objects that pass our latest selection algorithm with the PDR3 photometry. That is, a total of 91 % of the objects lie at z > 5.6. This demonstrates that the algorithm has very high efficiency, even though we are probing an unprecedentedly low-luminosity population down to M1450 ~ -21 mag. △ Less

Submitted 24 November, 2021; originally announced November 2021.

Comments: Accepted for publication in The Astrophysical Journal Supplement Series

HSC Year 1 cosmology results with the minimal bias method: HSC$\times$BOSS galaxy-galaxy weak lensing and BOSS galaxy clustering

Authors: Sunao Sugiyama, Masahiro Takada, Hironao Miyatake, Takahiro Nishimichi, Masato Shirasaki, Yosuke Kobayashi, Surhud More, Ryuichi Takahashi, Ken Osato, Masamune Oguri, Jean Coupon, Chiaki Hikage, Bau-Ching Hsieh, Yotaka Komiyama, Alexie Leauthaud, Xiangchong Li, Wentao Luo, Robert H. Lupton, Hitoshi Murayama, Atsushi J. Nishizawa, Youngsoo Park, Paul A. Price, Melanie Simet, Joshua S. Speagle, Michael A. Strauss , et al. (1 additional authors not shown)

Abstract: We present cosmological parameter constraints from a blinded joint analysis of galaxy-galaxy weak lensing, $Δ\!Σ(R)$, and projected correlation function, $w_\mathrm{p}(R)$, measured from the first-year HSC (HSC-Y1) data and SDSS spectroscopic galaxies over $0.15<z<0.7$. We use luminosity-limited samples as lens samples for $Δ\!Σ$ and as large-scale structure tracers for $w_\mathrm{p}$ in three red… ▽ More We present cosmological parameter constraints from a blinded joint analysis of galaxy-galaxy weak lensing, $Δ\!Σ(R)$, and projected correlation function, $w_\mathrm{p}(R)$, measured from the first-year HSC (HSC-Y1) data and SDSS spectroscopic galaxies over $0.15<z<0.7$. We use luminosity-limited samples as lens samples for $Δ\!Σ$ and as large-scale structure tracers for $w_\mathrm{p}$ in three redshift bins, and use the HSC-Y1 galaxy catalog to define a secure sample of source galaxies at $z_\mathrm{ph}>0.75$ for the $Δ\!Σ$ measurements, selected based on their photometric redshifts. For theoretical template, we use the "minimal bias" model for the cosmological clustering observables for the flat $Λ$CDM cosmological model. We compare the model predictions with the measurements in each redshift bin on large scales, $R>12$ and $8~h^{-1}\mathrm{Mpc}$ for $Δ\!Σ(R)$ and $w_\mathrm{p}(R)$, respectively, where the perturbation theory-inspired model is valid. When we employ weak priors on cosmological parameters, without CMB information, we find $S_8=0.936^{+0.092}_{-0.086}$, $σ_8=0.85^{+0.16}_{-0.11}$, and $Ω_\mathrm{m}=0.283^{+0.12}_{-0.035}$ for the flat $Λ$CDM model. Although the central value of $S_8$ appears to be larger than those inferred from other cosmological experiments, we find that the difference is consistent with expected differences due to sample variance, and our results are consistent with the other results to within the statistical uncertainties. (abriged) △ Less

Submitted 21 November, 2021; originally announced November 2021.

Comments: 24 pages, 19 figures, 4 tables, to be submitted to Phys. Rev. D

Report number: IPMU21-0078, YITP-21-126

Cosmological inference from the emulator based halo model II: Joint analysis of galaxy-galaxy weak lensing and galaxy clustering from HSC-Y1 and SDSS

Authors: Hironao Miyatake, Sunao Sugiyama, Masahiro Takada, Takahiro Nishimichi, Masato Shirasaki, Yosuke Kobayashi, Rachel Mandelbaum, Surhud More, Masamune Oguri, Ken Osato, Youngsoo Park, Ryuichi Takahashi, Jean Coupon, Chiaki Hikage, Bau-Ching Hsieh, Alexie Leauthaud, Xiangchong Li, Wentao Luo, Robert H. Lupton, Satoshi Miyazaki, Hitoshi Murayama, Atsushi J. Nishizawa, Paul A. Price, Melanie Simet, Joshua S. Speagle , et al. (3 additional authors not shown)

Abstract: We present high-fidelity cosmology results from a blinded joint analysis of galaxy-galaxy weak lensing ($Δ\!Σ$) and projected galaxy clustering ($w_{\rm p}$) 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<z<0.7$. We define luminosity-limited samples of SDSS galaxies to serve as the tracers of… ▽ More We present high-fidelity cosmology results from a blinded joint analysis of galaxy-galaxy weak lensing ($Δ\!Σ$) and projected galaxy clustering ($w_{\rm p}$) 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<z<0.7$. We define luminosity-limited samples of SDSS galaxies to serve as the tracers of $w_{\rm p}$ in three spectroscopic redshift bins, and as the lens samples for $Δ\!Σ$. For the $Δ\!Σ$ measurements, we select a single sample of 4 million source galaxies over 140 deg$^2$ from HSC-Y1 with photometric redshifts (photo-$z$) greater than 0.75, enabling a better handle of photo-$z$ errors by comparing the $Δ\!Σ$ amplitudes for the three lens redshift bins. For cosmological parameter inference, we use an input galaxy-halo connection model built on the {\tt Dark Emulator} package with a halo occupation distribution that includes nuisance parameters to marginalize over modeling uncertainties. We model the $Δ\!Σ$ and $w_{\rm p}$ measurements on scales from $R\simeq 3$ and $2\,h^{-1}{\rm Mpc}$, respectively, up to $30\,h^{-1}{\rm Mpc}$ assuming a flat $Λ$CDM cosmology. With various tests using mock catalogs described in Miyatake et al. (2021), we show that any bias in the clustering amplitude $S_8\equiv σ_8(Ω_{\rm m}/0.3)^{0.5}$ due to uncertainties in the galaxy-halo connection is less than $\sim50$\% of the statistical uncertainty of $S_8$, {\it unless} the assembly bias effect is unexpectedly large. Our best-fit models have $S_8=0.795^{+0.049}_{-0.042}$ (mode and 68\% credible interval) for the flat $Λ$CDM model; we find tighter constraints on the quantity $S_8(α=0.17)\equivσ_8(Ω_{\rm m}/0.3)^{0.17} =0.745^{+0.039}_{-0.031}$. (abriged) △ Less

Submitted 29 November, 2021; v1 submitted 3 November, 2021; originally announced November 2021.

Comments: 44 pages, 30 figures, 6 tables, submitted to Phys. Rev. D, YouTube video is available at https://www.youtube.com/watch?v=jRhQvuee97Q

Report number: IPMU21-0069, YITP-21-125

Precision analysis of the redshift-space galaxy bispectrum

Authors: Mikhail M. Ivanov, Oliver H. E. Philcox, Takahiro Nishimichi, Marko Simonović, Masahiro Takada, Matias Zaldarriaga

Abstract: We study the information content of the angle-averaged (monopole) redshift space galaxy bispectrum. The main novelty of our approach is the use of a systematic tree-level perturbation theory model that includes galaxy bias, IR resummation, and also accounts for nonlinear redshift space distortions, binning, and projection effects. We analyze data from the PT challenge simulations, whose cumulative… ▽ More We study the information content of the angle-averaged (monopole) redshift space galaxy bispectrum. The main novelty of our approach is the use of a systematic tree-level perturbation theory model that includes galaxy bias, IR resummation, and also accounts for nonlinear redshift space distortions, binning, and projection effects. We analyze data from the PT challenge simulations, whose cumulative volume of 566 $h^{-3}$Gpc$^3$ allows for a precise comparison to theoretical predictions. Fitting the power spectrum and bispectrum of our simulated data, and varying all necessary cosmological and nuisance parameters in a consistent Markov chain Monte Carlo analysis, we find that our tree-level bispectrum model is valid up to $k_{\max}=0.08~h{\rm Mpc}^{-1}$ (at $z=0.61$). We also find that inclusion of the bispectrum monopole improves constraints on cosmological parameters by $(5-15)\%$ relative to the power spectrum. The improvement is more significant for the quadratic bias parameters of our simulated galaxies, which we also show to deviate from biases of the host dark matter halos at the $\sim 3σ$ level. Finally, we adjust the covariance and scale cuts to match the volume of the BOSS survey, and estimate that within the minimal $Λ$CDM model the bispectrum data can tighten the constraint on the mass fluctuation amplitude $σ_8$ by roughly $10\%$. △ Less

Submitted 19 October, 2021; originally announced October 2021.

Comments: 52 pages, 11 figures

Report number: YITP-21-120, CERN-TH-2021-155

Full-shape cosmology analysis of SDSS-III BOSS galaxy power spectrum using emulator-based halo model: a $5\%$ determination of $σ_8$

Authors: Yosuke Kobayashi, Takahiro Nishimichi, Masahiro Takada, Hironao Miyatake

Abstract: We present the results obtained from the full-shape cosmology analysis of the redshift-space power spectra for 4 galaxy samples of the SDSS-III BOSS DR12 galaxy catalog over $0.2 < z < 0.75$. For the theoretical template, we use an emulator that was built from an ensemble set of $N$-body simulations, which enables fast and accurate computation of the redshift-space power spectrum of halos. Combini… ▽ More We present the results obtained from the full-shape cosmology analysis of the redshift-space power spectra for 4 galaxy samples of the SDSS-III BOSS DR12 galaxy catalog over $0.2 < z < 0.75$. For the theoretical template, we use an emulator that was built from an ensemble set of $N$-body simulations, which enables fast and accurate computation of the redshift-space power spectrum of halos. Combining with the halo occupation distribution to model the halo-galaxy connection, we can compute the redshift-space power spectrum of BOSS-like galaxies in the flat $Λ$CDM cosmology. In our cosmology inference, we use the power spectrum monopole, quadrupole and hexadecapole and include 7 nuisance parameters to model uncertainties in the halo-galaxy connection for each galaxy sample, but do not use any information on the abundance of galaxies. We demonstrate a validation of our analysis pipeline using the mock catalogs of BOSS-like galaxies, generated using different recipes of the halo-galaxy connection and including the assembly bias effect. Assuming weak priors on cosmological parameters, except for $Ω_{\rm b}h^2$ and $n_{\rm s}$, we show that our model well reproduces the BOSS power spectra. Including the power spectrum information up to $k_{\rm max}=0.25\,h{\rm Mpc}^{-1}$, we find $Ω_\mathrm{m}=0.301^{+0.012}_{-0.011}$, $H_0=68.2 \pm 1.4~\mathrm{km\,s}^{-1}\mathrm{Mpc}^{-1}$, and $σ_8=0.786^{+0.036}_{-0.037}$, for the mode and 68\% credible interval, after marginalization over nuisance parameters. We find little improvement in the cosmological parameters beyond a maximum wavelength $k_{\rm max}\simeq 0.2\,h\,{\rm Mpc}^{-1}$ due to the shot noise domination and marginalization of the halo-galaxy connection parameters. Our results are consistent with the Planck CMB results within $1σ$ statistical uncertainties. △ Less

Submitted 10 March, 2022; v1 submitted 13 October, 2021; originally announced October 2021.

Comments: 27 pages, 18 figures, accepted by Phys. Rev. D. The main results are replaced with the analysis of the updated measurement of BOSS DR12 power spectrum provided by Beutler & McDonald (2021, arXiv:2106.06324)

Report number: IPMU21-0058, YITP-21-112

Cosmological constraints without fingers of God

Authors: Mikhail M. Ivanov, Oliver H. E. Philcox, Marko Simonović, Matias Zaldarriaga, Takahiro Nishimichi, Masahiro Takada

Abstract: Non-linear redshift-space distortions ("fingers of God") are challenging to model analytically, a fact that limits the applicability of perturbation theory in redshift space as compared to real space. We show how this problem can be mitigated using a new observable, $Q_0$, which can be easily estimated from the redshift space clustering data and is approximately equal to the real space power spect… ▽ More Non-linear redshift-space distortions ("fingers of God") are challenging to model analytically, a fact that limits the applicability of perturbation theory in redshift space as compared to real space. We show how this problem can be mitigated using a new observable, $Q_0$, which can be easily estimated from the redshift space clustering data and is approximately equal to the real space power spectrum. The new statistic does not suffer from fingers of God and can be accurately described with perturbation theory down to $k_{\rm max}\simeq 0.4~h~\text{Mpc}^{-1}$. It can be straightforwardly included in the likelihood at negligible additional computational cost, and yields noticeable improvements on cosmological parameters compared to standard power spectrum multipole analyses. Using both simulations and observational data from the Baryon Oscillation Spectroscopic Survey, we show that improvements vary from $10\%$ to $100\%$ depending on the cosmological parameter considered, the galaxy sample and the survey volume. △ Less

Submitted 23 February, 2022; v1 submitted 30 September, 2021; originally announced October 2021.

Comments: 29 pages, 6 figures, minor edits, new BOSS likelihood used, available at github.com/oliverphilcox/full_shape_likelihoods

Report number: YITP-21-106, CERN-TH-2021-143

Journal ref: Phys. Rev. D 105, 043531 (2022)

Dark matter constraints with stacked gamma rays scales with the number of galaxies

Authors: Daiki Hashimoto, Atsushi J. Nishizawa, Masahiro Takada, Oscar Macias

Abstract: Low-surface-brightness galaxies (LSBGs) are interesting targets for searches of dark matter emission due to their low baryonic content. However, predicting their expected dark matter emissivities is difficult because of observational challenges in their distance measurements. Here we present a stacking method that makes use of catalogs of LSBGs and maps of unresolved gamma-ray emission measured by… ▽ More Low-surface-brightness galaxies (LSBGs) are interesting targets for searches of dark matter emission due to their low baryonic content. However, predicting their expected dark matter emissivities is difficult because of observational challenges in their distance measurements. Here we present a stacking method that makes use of catalogs of LSBGs and maps of unresolved gamma-ray emission measured by the Fermi Gamma-Ray Space Telescope. We show that, for relatively large number of LSBGs, individual distance measurements to the LSBGs are not necessary, instead the overall distance distribution of the population is sufficient in order to impose dark matter constraints. Further, we demonstrate that the effect of the covariance between two galaxies located closely -- at an angular distance comparable to the size of the Fermi point spread function -- is negligibly small. As a case in point, we apply our pipeline to a sample of 800 faint LSBGs discovered by Hyper Suprime-Cam and find that, the 95 per cent confidence level upper limits on the dark matter annihilation cross-section scales with inverse of the number LSBGs. In light of this linear dependence with the number of objects, we argue this methodology could provide extremely powerful limits if it is applied to the more than 105 LSBGs readily available with the Legacy Survey of Space and Time. △ Less

Submitted 29 April, 2022; v1 submitted 18 September, 2021; originally announced September 2021.

Comments: 12 pages, 7 figures, submitted to PRD

Third Data Release of the Hyper Suprime-Cam Subaru Strategic Program

Authors: Hiroaki Aihara, Yusra AlSayyad, Makoto Ando, Robert Armstrong, James Bosch, Eiichi Egami, Hisanori Furusawa, Junko Furusawa, Sumiko Harasawa, Yuichi Harikane, Bau-Ching Hsieh, Hiroyuki Ikeda, Kei Ito, Ikuru Iwata, Tadayuki Kodama, Michitaro Koike, Mitsuru Kokubo, Yutaka Komiyama, Xiangchong Li, Yongming Liang, Yen-Ting Lin, Robert H. Lupton, Nate B Lust, Lauren A. MacArthur, Ken Mawatari , et al. (42 additional authors not shown)

Abstract: The paper presents the third data release of Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP), a wide-field multi-band imaging survey with the Subaru 8.2m 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 ob… ▽ More The paper presents the third data release of Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP), a wide-field multi-band imaging survey with the Subaru 8.2m 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 square degrees in all five broad-band filters at the full depth ($\sim26$~mag at $5σ$) in the Wide layer. If we include partially observed area, the release covers 1,470 square degrees. The Deep and UltraDeep layers have $\sim80\%$ 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 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. The data release website is https://hsc-release.mtk.nao.ac.jp/. △ Less

Submitted 30 August, 2021; originally announced August 2021.

Comments: 25 pages, 19 figures, submitted to PASJ. Data available at https://hsc-release.mtk.nao.ac.jp/

Brightest Cluster Galaxies Are Statistically Special From $z=0.3$ to $z=1$

Authors: Roohi Dalal, Michael A. Strauss, Tomomi Sunayama, Masamune Oguri, Yen-Ting Lin, Song Huang, Youngsoo Park, Masahiro Takada

Abstract: We study Brightest Cluster Galaxies (BCGs) in $\sim5000$ galaxy clusters from the Hyper Suprime-Cam (HSC) Subaru Strategic Program. The sample is selected over an area of 830 $\textrm{deg}^2$ and is uniformly distributed in redshift over the range $z=0.3-1.0$. The clusters have stellar masses in the range $10^{11.8} - 10^{12.9} M_{\odot}$. We compare the stellar mass of the BCGs in each cluster to… ▽ More We study Brightest Cluster Galaxies (BCGs) in $\sim5000$ galaxy clusters from the Hyper Suprime-Cam (HSC) Subaru Strategic Program. The sample is selected over an area of 830 $\textrm{deg}^2$ and is uniformly distributed in redshift over the range $z=0.3-1.0$. The clusters have stellar masses in the range $10^{11.8} - 10^{12.9} M_{\odot}$. We compare the stellar mass of the BCGs in each cluster to what we would expect if their masses were drawn from the mass distribution of the other member galaxies of the clusters. The BCGs are found to be "special", in the sense that they are not consistent with being a statistical extreme of the mass distribution of other cluster galaxies. This result is robust over the full range of cluster stellar masses and redshifts in the sample, indicating that BCGs are special up to a redshift of $z=1.0$. However, BCGs with a large separation from the center of the cluster are found to be consistent with being statistical extremes of the cluster member mass distribution. We discuss the implications of these findings for BCG formation scenarios. △ Less

Submitted 25 August, 2021; originally announced August 2021.

Comments: 14 pages, 12 figures. This is a pre-copyedited, author-produced PDF of an article accepted for publication in MNRAS following peer review. The version of record is available online at: https://doi.org/10.1093/mnras/stab2363

Possible evidence of QCD axion stars in HSC and OGLE microlensing events

Authors: Sunao Sugiyama, Masahiro Takada, Alexander Kusenko

Abstract: Dark matter in the form of axions is expected to form miniclusters, and their dense regions can harbor compact axion stars. Such axion stars could be discovered by microlensing events. In particular, some candidate events reported by Subaru HSC and OGLE can be explained simultaneously if the axion stars with masses of the order of the Earth mass make up about $\sim27^{+7}_{013}$ percent of dark ma… ▽ More Dark matter in the form of axions is expected to form miniclusters, and their dense regions can harbor compact axion stars. Such axion stars could be discovered by microlensing events. In particular, some candidate events reported by Subaru HSC and OGLE can be explained simultaneously if the axion stars with masses of the order of the Earth mass make up about $\sim27^{+7}_{013}$ percent of dark matter. For QCD axions, this corresponds to the axion mass in the range $10^{-9}-10^{-6}$ eV, which is consistent with the experimental constraints, as well as the cosmological anthropic window of parameters. △ Less

Submitted 7 April, 2023; v1 submitted 6 August, 2021; originally announced August 2021.

Comments: 8 pages, 2 figures; latest version for the publication in Physics Letters B

Report number: IPMU21-0050

Mock catalogues of emission line galaxies based on the local mass density in dark-matter only simulations

Authors: Ken Osato, Takahiro Nishimichi, Masahiro Takada

Abstract: The high-precision measurement of spatial clustering of emission line galaxies (ELGs) is a primary objective for upcoming cosmological spectroscopic surveys. The source of strong emission of ELGs is nebular emission from surrounding ionized gas irradiated by massive short-lived stars in star-forming galaxies. As a result, ELGs are more likely to reside in newly-formed halos and this leads to a non… ▽ More The high-precision measurement of spatial clustering of emission line galaxies (ELGs) is a primary objective for upcoming cosmological spectroscopic surveys. The source of strong emission of ELGs is nebular emission from surrounding ionized gas irradiated by massive short-lived stars in star-forming galaxies. As a result, ELGs are more likely to reside in newly-formed halos and this leads to a nonlinear relation between ELG number density and matter density fields. In order to estimate the covariance matrix of cosmological observables, it is essential to produce many independent realisations to simulate ELG distributions for large survey volumes. To this end, we present a novel and fast scheme to populate ELGs in dark-matter only $N$-body simulations based on local density field. This method enables fast production of mock ELG catalogues suitable for verifying analysis methods and quantifying observational systematics in upcoming spectroscopic surveys and can populate ELGs in moderately high-density regions even though the halo structure cannot be resolved due to low resolution. The power spectrum of simulated ELGs is consistent with results of hydrodynamical simulations up to fairly small scales ($\lesssim 1 h \, \mathrm{Mpc}^{-1}$), and the simulated ELGs are more likely to be found in filamentary structures, which is consistent with results of semi-analytic and hydrodynamical simulations. Furthermore, we address the redshift-space power spectrum of simulated ELGs. The measured multipole moments of simulated ELGs clearly exhibit a weaker Finger-of-God effect than those of matter due to infalling motions towards halo centre, rather than random virial motions inside halos. △ Less

Submitted 13 January, 2022; v1 submitted 28 July, 2021; originally announced July 2021.

Comments: 10 pages, 9 figures, MNRAS in press

Report number: YITP-21-73

Journal ref: MNRAS, 511, 1131-1140 (2022)

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

Authors: Divya Rana, Surhud More, Hironao Miyatake, Takahiro Nishimichi, Masahiro Takada, Aaron S. G. Robotham, Andrew M. Hopkins, Benne W. Holwerda

Abstract: We utilize the galaxy shape catalogue from the first-year data release of the Subaru Hyper Suprime-cam Survey (HSC) 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$… ▽ More We utilize the galaxy shape catalogue from the first-year data release of the Subaru Hyper Suprime-cam Survey (HSC) 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 percent constraint on the amplitude of the scaling relation between halo mass and group luminosity with $\avg{M} = (0.81\pm 0.04)\times10^{14}\hinvMsun$ for $L_{\rm grp}=10^{11.5}\hinvsqLsun$, and a power-law index of $α=1.01\pm 0.07$. We constrain the amplitude of the scaling relation between halo mass and velocity dispersion to be $\avg{M}=(0.93\pm 0.05)\times10^{14}\hinvMsun$ for $σ=500 \kms$ and a power-law index to be $α=1.52\pm0.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. △ Less

Submitted 9 January, 2022; v1 submitted 12 July, 2021; originally announced July 2021.

Comments: 17 pages, 11 figures, Accepted for publication in MNRAS

Report number: YITP-21-74

The three-year shear catalog of the Subaru Hyper Suprime-Cam SSP Survey

Authors: Xiangchong Li, Hironao Miyatake, Wentao Luo, Surhud More, Masamune Oguri, Takashi Hamana, Rachel Mandelbaum, Masato Shirasaki, Masahiro Takada, Robert Armstrong, Arun Kannawadi, Satoshi Takita, Satoshi Miyazaki, Atsushi J. Nishizawa, Andrés A. Plazas Malagón, Michael A. Strauss, Masayuki Tanaka, Naoki Yoshida

Abstract: We present the galaxy shear catalog that will be used for the three-year cosmological weak gravitational lensing analyses using data from the Wide layer of the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP) Survey. The galaxy shapes are measured from the $i$-band imaging data acquired from 2014 to 2019 and calibrated with image simulations that resemble the observing conditions of the surv… ▽ More We present the galaxy shear catalog that will be used for the three-year cosmological weak gravitational lensing analyses using data from the Wide layer of the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP) Survey. The galaxy shapes are measured from the $i$-band imaging data acquired from 2014 to 2019 and calibrated with image simulations that resemble the observing conditions of the survey based on training galaxy images from the Hubble Space Telescope in the COSMOS region. The catalog covers an area of 433.48 deg$^2$ of the northern sky, split into six fields. The mean $i$-band seeing is 0.59 arcsec. With conservative galaxy selection criteria (e.g., $i$-band magnitude brighter than 24.5), the observed raw galaxy number density is 22.9 arcmin$^{-2}$, and the effective galaxy number density is 19.9 arcmin$^{-2}$. The calibration removes the galaxy property-dependent shear estimation bias to a level: $|δm|<9\times 10^{-3}$. The bias residual $δm$ shows no dependence on redshift in the range $0<z\leq 3$. We define the requirements for cosmological weak lensing science for this shear catalog, and quantify potential systematics in the catalog using a series of internal null tests for systematics related to point-spread function modelling and shear estimation. A variety of the null tests are statistically consistent with zero or within requirements, but (i) there is evidence for PSF model shape residual correlations; and (ii) star-galaxy shape correlations reveal additive systematics. Both effects become significant on $>1$ degree scales and will require mitigation during the inference of cosmological parameters using cosmic shear measurements. △ Less

Submitted 1 February, 2022; v1 submitted 30 June, 2021; originally announced July 2021.

Comments: 39 pages, 31 figures, 5 tables; PASJ (in press)

Hunting isolated neutron stars with proper motions from wide-area optical surveys

Authors: Daisuke Toyouchi, Kenta Hotokezaka, Masahiro Takada

Abstract: High-velocity neutron stars (HVNSs) that were kicked out from their birth location can be potentially identified with their large proper motions, and possibly with large parallax, when they come across the solar neighborhood. In this paper, we study the feasibility of hunting isolated HVNSs in wide-area optical surveys by modeling the evolution of NS luminosity taking into account spin-down and th… ▽ More High-velocity neutron stars (HVNSs) that were kicked out from their birth location can be potentially identified with their large proper motions, and possibly with large parallax, when they come across the solar neighborhood. In this paper, we study the feasibility of hunting isolated HVNSs in wide-area optical surveys by modeling the evolution of NS luminosity taking into account spin-down and thermal radiation. Assuming the upcoming 10-year VRO LSST observation, our model calculations predict that about 10 HVNSs mainly consisting of pulsars with ages of $10^4$--$10^5$ yr and thermally emitting NSs with $10^5$--$10^6$ yr are detectable. We find that a few NSs with effective temperature $< 5 \times 10^5$ K, which are likely missed in the current and future X-ray surveys, are also detectable. In addition to the standard neutron star cooling models, we consider a dark matter heating model. If such a strong heating exists we find that the detectable HVNSs would be significantly cooler, i.e., $\lesssim 5\times 10^5$ K. Thus, the future optical observation will give an unique NS sample, which can provide essential constraints on the NS cooling and heating mechanisms. Moreover, we suggest that providing HVNS samples with optical surveys is helpful for understanding the intrinsic kick-velocity distribution of NSs. △ Less

Submitted 9 June, 2021; originally announced June 2021.

Comments: 8 pages, 6 figures, submitted to MNRAS

An Optimal Estimator of Intrinsic Alignments for Star-forming Galaxies in IllustrisTNG Simulation

Authors: Jingjing Shi, Ken Osato, Toshiki Kurita, Masahiro Takada

Abstract: Emission line galaxies (ELGs), more generally star-forming galaxies, are valuable tracers of large-scale structure and therefore main targets of upcoming wide-area spectroscopic galaxy surveys. We propose a fixed-aperture shape estimator of each ELG for extracting the intrinsic alignment (IA) signal, and assess the performance of the method using image simulations of ELGs generated from the Illust… ▽ More Emission line galaxies (ELGs), more generally star-forming galaxies, are valuable tracers of large-scale structure and therefore main targets of upcoming wide-area spectroscopic galaxy surveys. We propose a fixed-aperture shape estimator of each ELG for extracting the intrinsic alignment (IA) signal, and assess the performance of the method using image simulations of ELGs generated from the IllustrisTNG simulation including observational effects such as the sky background noise. We show that our method enables a significant detection of the IA power spectrum with the linear-scale coefficient $A_{\rm IA}\simeq (13$--$15)\pm 3.0$ up to $z=2$, even from the small simulation volume $\sim0.009\,(h^{-1}{\rm Gpc})^3$, in contrast to the null detection with the standard method. Thus the ELG IA signal, measured with our method, opens up opportunities to exploit cosmology and galaxy physics in high-redshift universe. △ Less

Submitted 6 August, 2021; v1 submitted 25 April, 2021; originally announced April 2021.

Comments: accepted for publication in ApJ

Report number: YITP-21-42

The stellar mass in and around isolated central galaxies: connections to the total mass distribution through galaxy-galaxy lensing in the Hyper Suprime-Cam survey

Authors: Wenting Wang, Xiangchong Li, Jingjing Shi, Jiaxin Han, Naoki Yasuda, Yipeng Jing, Surhud More, Masahiro Takada, Hironao Miyatake, Atsushi J. Nishizawa

Abstract: Using photometric galaxies from the HSC survey, we measure the stellar mass density profiles for satellite galaxies as a function of the projected distance, $r_p$, to isolated central galaxies (ICGs) selected from SDSS/DR7 spectroscopic galaxies at $z\sim0.1$. By stacking HSC images, we also measure the projected stellar mass density profiles for ICGs and their stellar halos. The total mass distri… ▽ More Using photometric galaxies from the HSC survey, we measure the stellar mass density profiles for satellite galaxies as a function of the projected distance, $r_p$, to isolated central galaxies (ICGs) selected from SDSS/DR7 spectroscopic galaxies at $z\sim0.1$. By stacking HSC images, we also measure the projected stellar mass density profiles for ICGs and their stellar halos. The total mass distributions are further measured from HSC weak lensing signals. ICGs dominate within $\sim$0.15 times the halo virial radius ($0.15R_{200}$). The stellar mass versus total mass fractions drop with the increase in $r_p$ up to $\sim0.15R_{200}$, beyond which they are less than 1\% while stay almost constant, indicating the radial distribution of satellites trace dark matter. The total stellar mass in satellites is proportional to the virial mass of the host halo, $M_{200}$, for ICGs more massive than $10^{10.5}M_\odot$, i.e., $M_{\ast,\mathrm{sat}} \propto M_{200}$, whereas the relation between the stellar mass of ICGs $+$ stellar halos and $M_{200}$ is close to $M_{\ast,\mathrm{ICG+diffuse}}\propto M_{200}^{1/2}$. Below $10^{10.5}M_\odot$, the change in $M_{200}$ is much slower with the decrease in $M_{\ast,\mathrm{ICG+diffuse}}$. At fixed stellar mass, red ICGs are hosted by more massive dark matter halos and have more satellites. At $M_{200}\sim10^{12.7}M_\odot$, both $M_{\ast,\mathrm{sat}}$ and the fraction of stellar mass in satellites versus total stellar mass, $f_\mathrm{sat}$, tend to be slightly higher around blue ICGs, perhaps implying the late formation of blue galaxies. $f_\mathrm{sat}$ increases with the increase in both $M_{\ast,\mathrm{ICG+diffuse}}$ and $M_{200}$, and scales more linearly with $M_{200}$. We provide best-fitting formulas for these scaling relations and for red and blue ICGs separately. △ Less

Submitted 8 August, 2021; v1 submitted 12 April, 2021; originally announced April 2021.

Comments: accepted by ApJ - comments welcome - data available upon request

Finding gravitational-wave black holes with parallax microlensing

Authors: Satoshi Toki, Masahiro Takada

Abstract: The LIGO-Virgo gravitational-wave (GW) observation unveiled the new population of black holes (BHs) that appears to have an extended mass spectrum up to around $70M_\odot$, much heavier than the previously-believed mass range ($\sim 8M_\odot$). In this paper, we study the capability of a microlensing observation of stars in the Milky Way (MW) bulge region to identify BHs of GW mass scales, taking… ▽ More The LIGO-Virgo gravitational-wave (GW) observation unveiled the new population of black holes (BHs) that appears to have an extended mass spectrum up to around $70M_\odot$, much heavier than the previously-believed mass range ($\sim 8M_\odot$). In this paper, we study the capability of a microlensing observation of stars in the Milky Way (MW) bulge region to identify BHs of GW mass scales, taking into account the microlensing parallax characterized by the parameter $π_{\rm E}\propto M^{-1/2}$ ($M$ is the mass of a lens), which is a dimension-less quantity defined by the ratio of the astronomical unit to the projected Einstein radius. First, assuming that BHs follow the same spatial and velocity distributions of stars as predicted by the standard MW model, we show that microlensing events with long light curve timescales, $t_{\rm E}\gtrsim 100~{\rm days}$, and small parallax effects, $π_{\rm E}\sim 10^{-2}$, are dominated by BH lenses compared to stellar-mass lenses. Second, using a Markov chain Monte Carlo analysis of the simulated light curve, we show that BH lens candidates are securely identified on individual basis, if the parallax effect is detected or well constrained to the precision of a percent level in $π_{\rm E}$. We also discuss that a microlensing event of an intermediate-mass BH of $\sim 1000M_\odot$, if it occurs, can be identified in a distinguishable way from stellar-mass BHs. △ Less

Submitted 25 March, 2021; v1 submitted 24 March, 2021; originally announced March 2021.

Comments: 18 pages, 14 figures, 1 table

Cosmological inference from emulator based halo model I: Validation tests with HSC and SDSS mock catalogs

Authors: Hironao Miyatake, Yosuke Kobayashi, Masahiro Takada, Takahiro Nishimichi, Masato Shirasaki, Sunao Sugiyama, Ryuichi Takahashi, Ken Osato, Surhud More, Youngsoo Park

Abstract: We present validation tests of emulator-based halo model method for cosmological parameter inference, assuming hypothetical measurements of the projected correlation function of galaxies, $w_{\rm p}(R)$, and the galaxy-galaxy weak lensing, $Δ\!Σ(R)$, from the spectroscopic SDSS galaxies and the Hyper Suprime-Cam Year1 (HSC-Y1) galaxies. To do this, we use \textsc{Dark Emulator} developed in Nishim… ▽ More We present validation tests of emulator-based halo model method for cosmological parameter inference, assuming hypothetical measurements of the projected correlation function of galaxies, $w_{\rm p}(R)$, and the galaxy-galaxy weak lensing, $Δ\!Σ(R)$, from the spectroscopic SDSS galaxies and the Hyper Suprime-Cam Year1 (HSC-Y1) galaxies. To do this, we use \textsc{Dark Emulator} developed in Nishimichi et al. based on an ensemble of $N$-body simulations, which is an emulation package enabling a fast, accurate computation of halo clustering quantities for flat-geometry $w$CDM cosmologies. Adopting the halo occupation distribution, the emulator allows us to obtain model predictions of $Δ\!Σ$ and $w_{\rm p}$ for the SDSS-like galaxies at a few CPU seconds for an input set of parameters. We present performance and validation of the method by carrying out Markov Chain Monte Carlo analyses of the mock signals measured from a variety of mock catalogs that mimic the SDSS and HSC-Y1 galaxies. We show that the halo model method can recover the underlying true cosmological parameters to within the 68\% credible interval, except for the mocks including the assembly bias effect (although we consider the unrealistically-large amplitude of assembly bias effect). Even for the assembly bias mock, we demonstrate that the cosmological parameters can be recovered {\it if} the analysis is restricted to scales $R\gtrsim 10~h^{-1}{\rm Mpc}$. We also show that, by using a single population of source galaxies to infer the relative strengths of $Δ\!Σ$ for multiple lens samples at different redshifts, the joint probes method allows for self-calibration of photometric redshift errors and multiplicative shear bias. Thus we conclude that the emulator-based halo model method can be safely applied to the HSC-Y1 dataset, achieving a precision of $σ(S_8)\simeq 0.04$. △ Less

Submitted 31 December, 2020; originally announced January 2021.

Comments: 39 pages, 22 figures and 5 tables

Report number: IPMU21-0002, YITP-21-03

Super-sample covariance of the thermal Sunyaev-Zel'dovich effect

Authors: Ken Osato, Masahiro Takada

Abstract: The thermal Sunyaev-Zel'dovich (tSZ) effect is a powerful probe of cosmology. The statistical errors in the tSZ power spectrum measurements are dominated by the presence of massive clusters in a survey volume that are easy to identify on individual cluster basis. First, we study the impact of super sample covariance (SSC) on the tSZ power spectrum measurements, and find that the sample variance is… ▽ More The thermal Sunyaev-Zel'dovich (tSZ) effect is a powerful probe of cosmology. The statistical errors in the tSZ power spectrum measurements are dominated by the presence of massive clusters in a survey volume that are easy to identify on individual cluster basis. First, we study the impact of super sample covariance (SSC) on the tSZ power spectrum measurements, and find that the sample variance is dominated by the connected non-Gaussian (cNG) covariance arising mainly from Poisson number fluctuations of massive clusters in the survey volume. Second, we find that removing such individually-detected, massive clusters from the analysis significantly reduces the cNG contribution, thereby leading the SSC to be a leading source of the sample variance. We then show, based on Fisher analysis, that the power spectrum measured from the remaining diffuse tSZ effects can be used to obtain tight constraints on cosmological parameters as well as the hydrostatic mass bias parameter. Our method offers complementary use of individual tSZ cluster counts and the power spectrum measurements of diffuse tSZ signals for cosmology and intracluster gas physics. △ Less

Submitted 2 February, 2021; v1 submitted 8 October, 2020; originally announced October 2020.

Comments: 15 pages, 9 figures, PRD in press

Journal ref: Phys. Rev. D 103, 063501 (2021)

Testing Stochastic Gravitational Wave Signals from Primordial Black Holes with Optical Telescopes

Authors: Sunao Sugiyama, Volodymyr Takhistov, Edoardo Vitagliano, Alexander Kusenko, Misao Sasaki, Masahiro Takada

Abstract: Primordial black holes (PBHs) can constitute the predominant fraction of dark matter (DM) if PBHs reside in the currently unconstrained "sublunar" mass range. PBHs originating from scalar perturbations generated during inflation can naturally appear with a broad spectrum in a class of models. The resulting stochastic gravitational wave (GW) background generated from such PBH production can account… ▽ More Primordial black holes (PBHs) can constitute the predominant fraction of dark matter (DM) if PBHs reside in the currently unconstrained "sublunar" mass range. PBHs originating from scalar perturbations generated during inflation can naturally appear with a broad spectrum in a class of models. The resulting stochastic gravitational wave (GW) background generated from such PBH production can account for the recently reported North American Nanohertz Observatory for Gravitational Waves (NANOGrav) pulsar timing array data signal, and will be testable in future GW observations by interferometer-type experiments such as Laser Interferometer Space Antenna (LISA). We show that the broad mass function of such PBH DM is already being probed by Subaru Hyper Suprime-Cam (HSC) microlensing data and is consistent with a detected candidate event. Upcoming observations of HSC will be able to provide an independent definitive test of the stochastic GW signals originating from such PBH DM production scenarios. △ Less

Submitted 27 January, 2021; v1 submitted 5 October, 2020; originally announced October 2020.

Comments: 6 pages, 1 figure; matches published version

Report number: IPMU20-0105, YITP-20-125

Journal ref: Phys.Lett.B 814 (2021) 136097

A comparative study of satellite galaxies in Milky Way-like galaxies from HSC, DECaLS and SDSS

Authors: Wenting Wang, Masahiro Takada, Xiangchong Li, Scott G. Carlsten, Ting-Wen Lan, Jingjing Shi, Hironao Miyatake, Surhud More, Rachael L. Beaton, Robert Lupton, Yen-Ting Lin, Tian Qiu, Wentao Luo

Abstract: We conduct a comprehensive and statistical study of the luminosity functions (LFs) for satellite galaxies, by counting photometric galaxies from HSC, DECaLS and SDSS around isolated central galaxies (ICGs) and paired galaxies from the SDSS/DR7 spectroscopic sample. Results of different surveys show very good agreement. The satellite LFs can be measured down to $M_V\sim-10$, and for central primary… ▽ More We conduct a comprehensive and statistical study of the luminosity functions (LFs) for satellite galaxies, by counting photometric galaxies from HSC, DECaLS and SDSS around isolated central galaxies (ICGs) and paired galaxies from the SDSS/DR7 spectroscopic sample. Results of different surveys show very good agreement. The satellite LFs can be measured down to $M_V\sim-10$, and for central primary galaxies as small as $8.5<\log_{10}M_\ast/M_\odot<9.2$ and $9.2<\log_{10}M_\ast/M_\odot<9.9$, implying there are on average 3--8 satellites with $M_V<-10$ around LMC-mass ICGs. The bright end cutoff of satellite LFs and the satellite abundance are both sensitive to the magnitude gap between the primary and its companions, indicating galaxy systems with larger magnitude gaps are on average hosted by less massive dark matter haloes. By selecting primaries with stellar mass similar to our MW, we discovered that i) the averaged satellite LFs of ICGs with different magnitude gaps to their companions and of galaxy pairs with different colour or colour combinations all show steeper slopes than the MW satellite LF; ii) there are on average more satellites with $-15<M_V<-10$ than those in our MW; iii) there are on average 1.5 to 2.5 satellites with $M_V<-16$ around ICGs, consistent with our MW; iv) even after accounting for the large scatter predicted by numerical simulations, the MW satellite LF is uncommon at $M_V>-12$. Hence the MW and its satellite system are statistically atypical of our sample of MW-mass systems. In consequence, our MW is not a good representative of other MW-mass galaxies. Strong cosmological implications based on only MW satellites await additional discoveries of fainter satellites in extra-galactic systems. Interestingly, the MW satellite LF is typical among other MW-mass systems within 40~Mpc in the local Universe, perhaps implying the Local Volume is an under-dense region. △ Less

Submitted 6 November, 2020; v1 submitted 15 September, 2020; originally announced September 2020.

Comments: accepted by MNRAS, with more detailed discussions about the scatter

Power Spectrum of Intrinsic Alignments of Galaxies in IllustrisTNG

Authors: Jingjing Shi, Toshiki Kurita, Masahiro Takada, Ken Osato, Yosuke Kobayashi, Takahiro Nishimichi

Abstract: We present the 3-{\it dimensional} intrinsic alignment power spectra between the projected 2d galaxy shape/spin and the 3d tidal field across $0.1<k/h{\rm Mpc}^{-1}<60$ using cosmological hydrodynamical simulation, Illustris-TNG300, at redshifts ranging from $0.3$ to $2$. The shape-tidal field alignment increases with galaxy mass and the linear alignment coefficient $A_{\rm IA}$, defined with resp… ▽ More We present the 3-{\it dimensional} intrinsic alignment power spectra between the projected 2d galaxy shape/spin and the 3d tidal field across $0.1<k/h{\rm Mpc}^{-1}<60$ using cosmological hydrodynamical simulation, Illustris-TNG300, at redshifts ranging from $0.3$ to $2$. The shape-tidal field alignment increases with galaxy mass and the linear alignment coefficient $A_{\rm IA}$, defined with respect to the primordial tidal field, is found to have weak redshift dependence. We also show a promising detection of the shape/spin-tidal field alignments for stellar mass limited samples and a weak or almost null signal for star-forming galaxies for the TNG300 volume, $\sim 0.01~(h^{-1}{\rm Gpc})^3$. We further study the morphology and environmental dependence of the intrinsic alignment power spectra. The shape of massive disk- and spheroid-galaxies tend to align with the tidal field. The spin of low mass disks (and spheroids at low redshifts) tend to be parallel with the tidal field, while the spin of massive spheroids and disks tend to be perpendicular to tidal field. The shape and spin of massive centrals align with the tidal field at both small and large scales. Satellites show a radial alignment within the one-halo term region, and low mass satellites have an intriguing alignment signal in the two-halo term region. We also forecast a feasibility to measure the intrinsic alignment power spectrum for spectroscopic and imaging surveys such as Subaru HSC/PFS and DESI. Our results thus suggest that galaxy intrinsic alignment can be used as a promising tool for constraining the galaxy formation models. △ Less

Submitted 28 January, 2021; v1 submitted 1 September, 2020; originally announced September 2020.

Comments: accepted for publication in JCAP, major change made after first version

Report number: YITP-20-113

Validating a minimal galaxy bias method for cosmological parameter inference using HSC-SDSS mock catalogs

Authors: Sunao Sugiyama, Masahiro Takada, Yosuke Kobayashi, Hironao Miyatake, Masato Shirasaki, Takahiro Nishimichi, Youngsoo Park

Abstract: We assess the performance of a perturbation theory inspired method for inferring cosmological parameters from the joint measurements of galaxy-galaxy weak lensing ($ΔΣ$) and the projected galaxy clustering ($w_{\rm p}$). To do this, we use a wide variety of mock galaxy catalogs constructed based on a large set of $N$-body simulations that mimic the Subaru HSC-Y1 and SDSS galaxies, and apply the me… ▽ More We assess the performance of a perturbation theory inspired method for inferring cosmological parameters from the joint measurements of galaxy-galaxy weak lensing ($ΔΣ$) and the projected galaxy clustering ($w_{\rm p}$). To do this, we use a wide variety of mock galaxy catalogs constructed based on a large set of $N$-body simulations that mimic the Subaru HSC-Y1 and SDSS galaxies, and apply the method to the mock signals to address whether to recover the underlying true cosmological parameters in the mocks. We find that, as long as the appropriate scale cuts, $12$ and $8~h^{-1}{\rm Mpc}$ for $ΔΣ$ and $w_{\rm p}$ respectively, are adopted, a "minimal-bias" model using the linear bias parameter $b_1$ alone and the nonlinear matter power spectrum can recover the true cosmological parameters (here focused on $Ω_{\rm m}$ and $σ_8$) to within the 68% credible interval, for all the mocks we study including one in which an assembly bias effect is implemented. This is as expected if physical processes inherent in galaxy formation/evolution are confined to local, small scales below the scale cut, and thus implies that real-space observables have an advantage in filtering out the impact of small-scale nonlinear effects in parameter estimation, compared to their Fourier-space counterparts. In addition, we find that a theoretical template including the higher-order bias contributions such as nonlinear bias parameter $(b_2)$ does not improve the cosmological constraints, but rather leads to a larger parameter bias compared to the baseline $b_1$-method. △ Less

Submitted 16 August, 2020; originally announced August 2020.

Comments: Submitted to PRD, A revised version in response to the referee's comments

Report number: IPMU20-0063, YITP-20-80

Journal ref: Phys. Rev. D 102, 083520 (2020)

Impacts of pre-initial conditions on anisotropic separate universe simulations: a boosted tidal response in the epoch of reionization

Authors: Shogo Masaki, Takahiro Nishimichi, Masahiro Takada

Abstract: To generate initial conditions for cosmological $N$-body simulations, one needs to prepare a uniform distribution of simulation particles, so-called the pre-initial condition (pre-IC). The standard method to construct the pre-IC is to place the particles on the lattice grids evenly spaced in the three-dimensional spatial coordinates. However, even after the initial displacement of each particle ac… ▽ More To generate initial conditions for cosmological $N$-body simulations, one needs to prepare a uniform distribution of simulation particles, so-called the pre-initial condition (pre-IC). The standard method to construct the pre-IC is to place the particles on the lattice grids evenly spaced in the three-dimensional spatial coordinates. However, even after the initial displacement of each particle according to cosmological perturbations, the particle distribution remains to display an artificial anisotropy. Such an artifact causes systematic effects in simulations at later time until the evolved particle distribution sufficiently erases the initial anisotropy. In this paper, we study the impacts of the pre-IC on the anisotropic separate universe simulation, where the effect of large-scale tidal field on structure formation is taken into account using the anisotropic expansion in a local background (simulation volume). To quantify the impacts, we compare the simulations employing the standard grid pre-IC and the glass one, where the latter is supposed to suppress the initial anisotropy. We show that the artificial features in the grid pre-IC simulations are seen until $z\sim 9$, while the glass pre-IC simulations appear to be stable and accurate over the range of scales we study. From these results we find that a coupling of the large-scale tidal field with matter clustering is enhanced compared to the leading-order prediction of perturbation theory in the quasi non-linear regime in the redshift range $5\lesssim z\lesssim 15$, indicating the importance of tidal field on structure formation at such high redshifts, e.g. during the epoch of reionization. △ Less

Submitted 21 October, 2020; v1 submitted 16 July, 2020; originally announced July 2020.

Comments: 11 pages, 8 figures, 1 table, accepted for publication in MNRAS, added a subsection and Appendix

Report number: YITP-20-92, IPMU20-0080

Imprint of anisotropic primordial non-Gaussianity on halo intrinsic alignments in simulations

Authors: Kazuyuki Akitsu, Toshiki Kurita, Takahiro Nishimichi, Masahiro Takada, Satoshi Tanaka

Abstract: Using $N$-body simulations of cosmological large-scale structure formation, for the first time, we show that the anisotropic primordial non-Gaussianity (PNG) causes a scale-dependent modification, given by $1/k^2$ at small $k$ limit, in the three-dimensional power spectra of halo shapes (intrinsic alignments), whilst the conventional power spectrum of halo number density field remains unaffected.… ▽ More Using $N$-body simulations of cosmological large-scale structure formation, for the first time, we show that the anisotropic primordial non-Gaussianity (PNG) causes a scale-dependent modification, given by $1/k^2$ at small $k$ limit, in the three-dimensional power spectra of halo shapes (intrinsic alignments), whilst the conventional power spectrum of halo number density field remains unaffected. We discuss that wide-area imaging and spectrocopic surveys observing the same region of the sky allow us to constrain the quadrupole PNG coefficient $f_{\rm NL}^{s=2}$ at a precision comparable with or better than that of the cosmic microwave background. △ Less

Submitted 22 March, 2021; v1 submitted 7 July, 2020; originally announced July 2020.

Comments: 9 pages, 4 figures; an accepted version for publication in PRD

Report number: IPMU20-0075, YITP-20-85

Journal ref: Phys. Rev. D 103, 083508 (2021)

Hunting gravitational wave black holes with microlensing

Authors: Natasha S. Abrams, Masahiro Takada

Abstract: Gravitational microlensing is a powerful tool to search for a population of invisible black holes (BHs) in the Milky Way (MW), including isolated BHs and binary BHs at wide orbits that are complementary to gravitational wave observations. By monitoring highly populated regions of source stars like the MW bulge region, one can pursue microlensing events due to these BHs. We find that if BHs have a… ▽ More Gravitational microlensing is a powerful tool to search for a population of invisible black holes (BHs) in the Milky Way (MW), including isolated BHs and binary BHs at wide orbits that are complementary to gravitational wave observations. By monitoring highly populated regions of source stars like the MW bulge region, one can pursue microlensing events due to these BHs. We find that if BHs have a Salpeter-like mass function extended beyond $30M_\odot$ and a similar velocity and spatial structure to stars in the Galactic bulge and disk regions, the BH population is a dominant source of the microlensing events at long timescales of the microlensing light curve $\gtrsim 100~$days. This is due to a boosted sensitivity of the microlensing event rate to lens mass, given as $M^2$, for such long-timescale events. A monitoring observation of $2 \times 10^{10}$ stars in the bulge region over 10 years with the Rubin Observatory Legacy Survey of Space and Time (LSST) would enable one to find about $6\times 10^5$ BH microlensing events. We evaluate the efficiency of potential LSST cadences for characterizing the light curves of BH microlensing and find that nearly all events of long timescales can be detected. △ Less

Submitted 5 November, 2020; v1 submitted 9 June, 2020; originally announced June 2020.

Comments: 16 pages, 8 figures: accepted for publication in ApJ