Showing 1–50 of 140 results for author: Nishimichi, T

A Parameter-Masked Mock Data Challenge for Beyond-Two-Point Galaxy Clustering Statistics

Authors: Beyond-2pt Collaboration, :, Elisabeth Krause, Yosuke Kobayashi, Andrés N. Salcedo, Mikhail M. Ivanov, Tom Abel, Kazuyuki Akitsu, Raul E. Angulo, Giovanni Cabass, Sofia Contarini, Carolina Cuesta-Lazaro, ChangHoon Hahn, Nico Hamaus, Donghui Jeong, Chirag Modi, Nhat-Minh Nguyen, Takahiro Nishimichi, Enrique Paillas, Marcos Pellejero Ibañez, Oliver H. E. Philcox, Alice Pisani, Fabian Schmidt, Satoshi Tanaka, Giovanni Verza , et al. (2 additional authors not shown)

Abstract: The last few years have seen the emergence of a wide array of novel techniques for analyzing high-precision data from upcoming galaxy surveys, which aim to extend the statistical analysis of galaxy clustering data beyond the linear regime and the canonical two-point (2pt) statistics. We test and benchmark some of these new techniques in a community data challenge "Beyond-2pt", initiated during the… ▽ More The last few years have seen the emergence of a wide array of novel techniques for analyzing high-precision data from upcoming galaxy surveys, which aim to extend the statistical analysis of galaxy clustering data beyond the linear regime and the canonical two-point (2pt) statistics. We test and benchmark some of these new techniques in a community data challenge "Beyond-2pt", initiated during the Aspen 2022 Summer Program "Large-Scale Structure Cosmology beyond 2-Point Statistics," whose first round of results we present here. The challenge dataset consists of high-precision mock galaxy catalogs for clustering in real space, redshift space, and on a light cone. Participants in the challenge have developed end-to-end pipelines to analyze mock catalogs and extract unknown ("masked") cosmological parameters of the underlying $Λ$CDM models with their methods. The methods represented are density-split clustering, nearest neighbor statistics, BACCO power spectrum emulator, void statistics, LEFTfield field-level inference using effective field theory (EFT), and joint power spectrum and bispectrum analyses using both EFT and simulation-based inference. In this work, we review the results of the challenge, focusing on problems solved, lessons learned, and future research needed to perfect the emerging beyond-2pt approaches. The unbiased parameter recovery demonstrated in this challenge by multiple statistics and the associated modeling and inference frameworks supports the credibility of cosmology constraints from these methods. The challenge data set is publicly available and we welcome future submissions from methods that are not yet represented. △ Less

Submitted 3 May, 2024; originally announced May 2024.

Comments: New submissions welcome! Challenge data available at https://github.com/ANSalcedo/Beyond2ptMock

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

Emulating power spectra for pre- and post-reconstructed galaxy samples

Authors: Yuting Wang, Ruiyang Zhao, Zhongxu Zhai, Kazuya Koyama, Will J. Percival, Hong Guo, Yin Li, Gong-Bo Zhao, Takahiro Nishimichi, Héctor Gil-Marín, Yonghao Feng, Hanyu Zhang, Yi Wu

Abstract: The small-scale linear information in galaxy samples typically lost during non-linear growth can be restored to a certain level by the density field reconstruction, which has been demonstrated for improving the precision of the baryon acoustic oscillations (BAO) measurements. As proposed in the literature, a joint analysis of the power spectrum before and after the reconstruction enables an effici… ▽ More The small-scale linear information in galaxy samples typically lost during non-linear growth can be restored to a certain level by the density field reconstruction, which has been demonstrated for improving the precision of the baryon acoustic oscillations (BAO) measurements. As proposed in the literature, a joint analysis of the power spectrum before and after the reconstruction enables an efficient extraction of information carried by high-order statistics. However, the statistics of the post-reconstruction density field are difficult to model. In this work, we circumvent this issue by developing an accurate emulator for the pre-reconstructed, post-reconstructed, and cross power spectra ($P_{\rm pre}$, $P_{\rm post}$, $P_{\rm cross}$) up to $k=0.5~h~{\rm Mpc^{-1}}$ based on the \textsc{Dark Quest} N-body simulations. The accuracy of the emulator is at percent level, namely, the error of the emulated monopole and quadrupole of the power spectra is less than $1\%$ and $10\%$ of the ground truth, respectively. A fit to an example power spectra using the emulator shows that the constraints on cosmological parameters get largely improved using $P_{\rm pre}$+$P_{\rm post}$+$P_{\rm cross}$ with $k_{\rm max}=0.25~h~{\rm Mpc^{-1}}$, compared to that derived from $P_{\rm pre}$ alone, namely, the constraints on ($Ω_m$, $H_0$, $σ_8$) are tightened by $\sim41 \%-55\%$, and the uncertainties of the derived BAO and RSD parameters ($α_{\perp}$, $α_{||}$, $fσ_8$) shrink by $\sim 28\%-54\%$, respectively. This highlights the complementarity among $P_{\rm pre}$, $P_{\rm post}$ and $P_{\rm cross}$, which demonstrates the efficiency and practicability of a joint $P_{\rm pre}$, $P_{\rm post}$ and $P_{\rm cross}$ analysis for cosmological implications. △ Less

Submitted 25 February, 2024; v1 submitted 9 November, 2023; originally announced November 2023.

Comments: 19 pages, 11 figures, 2 tables; accepted for publication in ApJ

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

Nonlinear redshift space distortion in halo ellipticity correlations: Analytical model and N-body simulations

Authors: Teppei Okumura, Atsushi Taruya, Toshiki Kurita, Takahiro Nishimichi

Abstract: We present an analytic model of nonlinear correlators of galaxy/halo ellipticities in redshift space. The three-dimensional ellipticity field is not affected by the redshift-space distortion (RSD) at linear order, but by the nonlinear one, known as the Finger-of-God effect, caused by the coordinate transformation from real to redshift space. Adopting a simple Gaussian damping function to describe… ▽ More We present an analytic model of nonlinear correlators of galaxy/halo ellipticities in redshift space. The three-dimensional ellipticity field is not affected by the redshift-space distortion (RSD) at linear order, but by the nonlinear one, known as the Finger-of-God effect, caused by the coordinate transformation from real to redshift space. Adopting a simple Gaussian damping function to describe the nonlinear RSD effect and the nonlinear alignment model for the relation between the observed ellipticity and underlying tidal fields, we derive analytic formulas for the multipole moments of the power spectra of the ellipticity field in redshift space expanded in not only the associated Legendre basis, a natural basis for the projected galaxy shape field, but also the standard Legendre basis, conventionally used in literature. The multipoles of the corresponding correlation functions of the galaxy shape field are shown to be expressed by a simple Hankel transform, as is the case for those of the conventional galaxy density correlations. We measure these multipoles of the power spectra and correlation functions of the halo ellipticity field using large-volume N-body simulations. We then show that the measured alignment signals can be better predicted by our nonlinear model than the existing linear alignment model. The formulas derived here have already been used to place cosmological constraints using from the redshift-space correlation functions of the galaxy shape field measured from the Sloan Digital Sky Survey (Okumura and Taruya, 2023). △ Less

Submitted 31 March, 2024; v1 submitted 11 October, 2023; originally announced October 2023.

Comments: 26 pages, 9 figures, 3 tables, accepted to PRD

Report number: YITP-23-128

Journal ref: Phys. Rev. D 109, 103501 (2024)

Multi-stream radial structure of cold dark matter haloes from particle trajectories: deep inside splashback radius

Authors: Yohsuke Enomoto, Takahiro Nishimichi, Atsushi Taruya

Abstract: By tracking trajectories of dark matter (DM) particles accreting onto haloes in cosmological $N$-body simulations, we investigate the radial phase-space distribution of cold dark matter (CDM) haloes, paying attention to their inner regions deep inside the halo boundary called the splashback radius, where the particles undergo multi-stream flows. Improving the analysis by Sugiura et al., we classif… ▽ More By tracking trajectories of dark matter (DM) particles accreting onto haloes in cosmological $N$-body simulations, we investigate the radial phase-space distribution of cold dark matter (CDM) haloes, paying attention to their inner regions deep inside the halo boundary called the splashback radius, where the particles undergo multi-stream flows. Improving the analysis by Sugiura et al., we classify DM particles by the number of apocenter passages, $p$, and count it up to $p=40$ for each halo over a wide mass range. Quantifying the radial density profile for particles having the same value of $p$, we find that it generally exhibits a double-power law feature, whose indices of inner and outer slopes are well-described by $-1$ and $-8$, respectively. Its characteristic scale and density are given as a simple fitting function of $p$, with a weak halo mass dependence. Interestingly, summing up these double-power law profiles beyond $p=40$ reproduces well the total density profile of simulated haloes. The double-power law nature is persistent and generic not only in mass-selected haloes but also in haloes selected in different criteria. Our results are compared with self-similar solutions that describe the stationary and spherical accretion of DM. We find that even when introducing a non-zero angular momentum, none of them explain the radial multi-stream structure. The analysis with particle trajectories tracing back to higher redshifts suggests that the double-power law nature has been established during an early accretion phase and remains stable. △ Less

Submitted 8 December, 2023; v1 submitted 24 September, 2023; originally announced September 2023.

Comments: 24 pages, 20 figures. Version consistent with that accepted by MNRAS

Report number: YITP-23-117

Mock Observatory: two thousand lightcone mock catalogues of luminous red galaxies from the Hyper Suprime-Cam Survey for the cosmological large-scale analysis

Authors: Shogo Ishikawa, Teppei Okumura, Takahiro Nishimichi

Abstract: Estimating a reliable covariance matrix for correlation functions of galaxies is a crucial task to obtain accurate cosmological constraints from galaxy surveys. We generate 2,000 independent lightcone mock luminous red galaxy (LRGs) catalogues at $0.3 \leq z \leq 1.25$, designed to cover CAMIRA LRGs observed by the Subaru Hyper Suprime-Cam Subaru Strategic Programme (HSC SSP). We first produce ful… ▽ More Estimating a reliable covariance matrix for correlation functions of galaxies is a crucial task to obtain accurate cosmological constraints from galaxy surveys. We generate 2,000 independent lightcone mock luminous red galaxy (LRGs) catalogues at $0.3 \leq z \leq 1.25$, designed to cover CAMIRA LRGs observed by the Subaru Hyper Suprime-Cam Subaru Strategic Programme (HSC SSP). We first produce full-sky lightcone halo catalogues using a COmoving Lagrangian Acceleration (COLA) technique, and then trim them to match the footprints of the HSC SSP S20A Wide layers. The mock LRGs are subsequently populated onto the trimmed halo catalogues according to the halo occupation distribution model constrained by the observed CAMIRA LRGs. The stellar mass ($M_{\star}$) is assigned to each LRG by the subhalo abundance-matching technique using the observed stellar-mass functions of CAMIRA LRGs. We evaluate photometric redshifts (photo-$z$) of mock LRGs by incorporating the photo-$z$ scatter, which is derived from the observed $M_{\star}$--photo-$z$-scatter relations of the CAMIRA LRGs. We validate the constructed full-sky halo and lightcone LRG mock catalogues by comparing their angular clustering statistics (i.e., power spectra and correlation functions) with those measured from the halo catalogues of full $N$-body simulations and the CAMIRA LRG catalogues from the HSC SSP, respectively. We detect clear signatures of baryon acoustic oscillations (BAOs) from our mock LRGs, whose angular scales are well consistent with theoretical predictions. These results demonstrate that our mock LRGs can be used to evaluate covariance matrices at large scales and provide predictions for the BAO detectability and cosmological constraints. △ Less

Submitted 1 March, 2024; v1 submitted 7 August, 2023; originally announced August 2023.

Comments: 14 pages, 11 figures, accepted for publication in MNRAS

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

Robustness of Baryon Acoustic Oscillations Measurements with Photometric Redshift Uncertainties

Authors: Keitaro Ishikawa, Tomomi Sunayama, Atsushi J. Nishizawa, Hironao Miyatake, Takahiro Nishimichi

Abstract: We investigate the robustness of baryon acoustic oscillations (BAO) measurements with a photometric galaxy sample using mock galaxy catalogues with various sizes of photometric redshift (photo-$z$) uncertainties. We first conduct the robustness of BAO measurements, assuming we have a perfect knowledge of photo-$z$ uncertainties. We find that the BAO shift parameter $α$ can be constrained in an unb… ▽ More We investigate the robustness of baryon acoustic oscillations (BAO) measurements with a photometric galaxy sample using mock galaxy catalogues with various sizes of photometric redshift (photo-$z$) uncertainties. We first conduct the robustness of BAO measurements, assuming we have a perfect knowledge of photo-$z$ uncertainties. We find that the BAO shift parameter $α$ can be constrained in an unbiased manner even for 3% photometric redshift uncertainties up to $z\sim 1$. For instance, $α=1.006 \pm 0.078$ with 95% confidence level is obtained from 3% photo-$z$ uncertainty data at $z=1.03$ using the sample of $M_* \ge 10^{10.25} M_{\odot}/h^2$. We also find that a sparse galaxy sample, e.g. $<2\times10^{-4}$ [$h$ Mpc$^{-1}]^3$ causes additional noise in the covariance matrix calculation and can bias the constraint on $α$. Following this, we look into the scenario where incorrect photometric redshift uncertainties are assumed in the fitting model. We find that underestimating the photo-$z$ uncertainty leads to a degradation in the constraining power on $α$. However, the constrained value of $α$ is not biased. We also quantify the constraining power on $Ω_{\rm m0}$ assuming the LSST-like covariance and find that the 95% confidence level is $σ(Ω_{\rm m0})\sim0.03$-$0.05$ corresponding to the photo-$z$ uncertainties of 1% to 3% respectively. Finally, we examine whether the skewness in the photometric redshift can bias the constraint on $α$ and confirm that the constraint on $α$ is unbiased, even assuming a Gaussian photo-$z$ uncertainty in our model. △ Less

Submitted 13 July, 2023; v1 submitted 2 June, 2023; originally announced June 2023.

Comments: 11 pages, 11 figures, 3 tables, submitted to MNRAS

Perturbation theory challenge for cosmological parameters estimation II.: Matter power spectrum in redshift space

Authors: Ken Osato, Takahiro Nishimichi, Atsushi Taruya, Francis Bernardeau

Abstract: Constraining cosmological parameters from large-scale structure observations requires precise and accurate tools to compute its properties. While perturbation theory (PT) approaches can serve this purpose, exploration of large parameter space is challenging due to the potentially large computational cost of such calculations. In this study, we show that a response function approach applied to the… ▽ More Constraining cosmological parameters from large-scale structure observations requires precise and accurate tools to compute its properties. While perturbation theory (PT) approaches can serve this purpose, exploration of large parameter space is challenging due to the potentially large computational cost of such calculations. In this study, we show that a response function approach applied to the regularized PT (RegPT) model at 2-loop order, plus correction terms induced by redshift space distortion effects, can reduce the runtime by a factor of 50 compared to direct integration. We illustrate the performance of this approach by performing the parameter inference of five fundamental cosmological parameters from the redshift space power spectrum measured from $N$-body simulations as mock measurements, and inferred cosmological parameters are directly compared with parameters used to generate initial conditions of the simulations. From this \textit{PT challenge} analysis, the constraining power of cosmological parameters and parameter biases are quantified with the survey volume and galaxy number density expected for the \textit{Euclid} mission at the redshift $z = 1$ as a function of the maximum wave-number of data points $k_\mathrm{max}$. We find that RegPT with correction terms reproduces the input cosmological parameters without bias up to maximum wave-number $k_\mathrm{max}= 0.18 \, h\,\mathrm{Mpc}^{-1}$. Moreover, RegPT+, which introduces one free parameter to RegPT to handle the damping feature on small scales, delivers the best performance among the examined models and achieves tighter constraints without significant parameter bias for higher maximum wave-number $k_\mathrm{max} = 0.21 \, h\,\mathrm{Mpc}^{-1}$. △ Less

Submitted 9 March, 2024; v1 submitted 2 May, 2023; originally announced May 2023.

Comments: 32 pages, 21 figures, PRD in press, codes are available at https://github.com/0satoken/Eclairs

Report number: YITP-23-32

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

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

Machine Learning for Observational Cosmology

Authors: Kana Moriwaki, Takahiro Nishimichi, Naoki Yoshida

Abstract: An array of large observational programs using ground-based and space-borne telescopes is planned in the next decade. The forthcoming wide-field sky surveys are expected to deliver a sheer volume of data exceeding an exabyte. Processing the large amount of multiplex astronomical data is technically challenging, and fully automated technologies based on machine learning and artificial intelligence… ▽ More An array of large observational programs using ground-based and space-borne telescopes is planned in the next decade. The forthcoming wide-field sky surveys are expected to deliver a sheer volume of data exceeding an exabyte. Processing the large amount of multiplex astronomical data is technically challenging, and fully automated technologies based on machine learning and artificial intelligence are urgently needed. Maximizing scientific returns from the big data requires community-wide efforts. We summarize recent progress in machine learning applications in observational cosmology. We also address crucial issues in high-performance computing that are needed for the data processing and statistical analysis. △ Less

Submitted 25 May, 2023; v1 submitted 28 March, 2023; originally announced March 2023.

Comments: 55 pages, 8 figures, accepted for publication in Reports on Progress in Physics

Report number: YITP-23-40

Cosmology with the Galaxy Bispectrum Multipoles: Optimal Estimation and Application to BOSS Data

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

Abstract: We present a framework for self-consistent cosmological analyses of the full-shape anisotropic bispectrum, including the quadrupole $(\ell=2)$ and hexadecapole $(\ell=4)$ moments. This features a novel window-free algorithm for extracting the latter quantities from data, derived using a maximum-likelihood prescription. Furthermore, we introduce a theoretical model for the bispectrum multipoles (wh… ▽ More We present a framework for self-consistent cosmological analyses of the full-shape anisotropic bispectrum, including the quadrupole $(\ell=2)$ and hexadecapole $(\ell=4)$ moments. This features a novel window-free algorithm for extracting the latter quantities from data, derived using a maximum-likelihood prescription. Furthermore, we introduce a theoretical model for the bispectrum multipoles (which does not introduce new free parameters), and test both aspects of the pipeline on several high-fidelity mocks, including the PT Challenge suite of gigantic cumulative volume. This establishes that the systematic error is significantly below the statistical threshold, both for the measurement and modeling. As a realistic example, we extract the large-scale bispectrum multipoles from BOSS DR12 and analyze them in combination with the power spectrum data. Assuming a minimal $Λ$CDM model, with a BBN prior on the baryon density and a \textit{Planck} prior on $n_s$, we can extract the remaining cosmological parameters directly from the clustering data. The inclusion of the unwindowed higher-order $(\ell>0)$ large-scale bispectrum multipoles is found to moderately improve one-dimensional cosmological parameter posteriors (at the $5\%-10\%$ level), though these multipoles are detected only in three out of four BOSS data segments at $\approx 5σ$. Combining information from the power spectrum and bispectrum multipoles, the real space power spectrum, and the post-reconstructed BAO data, we find $H_0 = 68.2\pm 0.8~\mathrm{km}\,\mathrm{s}^{-1}\mathrm{Mpc}^{-1}$, $Ω_m =0.33\pm 0.01$ and $σ_8 = 0.736\pm 0.033$ (the tightest yet found in perturbative full-shape analyses). Our estimate of the growth parameter $S_8=0.77\pm 0.04$ agrees with both weak lensing and CMB results. △ Less

Submitted 8 February, 2023; originally announced February 2023.

Comments: 40 pages, 8 figures, 6 tables; estimators and data are publicly available at https://github.com/oliverphilcox/Spectra-Without-Windows

Report number: YITP-23-13, CERN-TH-2023-022

Universal multi-stream radial structures of cold dark matter halos

Authors: Yohsuke Enomoto, Takahiro Nishimichi, Atsushi Taruya

Abstract: Virialized halos of cold dark matter generically exhibit multi-stream structures of accreted dark matter within an outermost radial caustic known as the splashback radius. By tracking the particle trajectories that accrete onto the halos in cosmological $N$-body simulations, we count their number of apocenter passages ($p$), and use them to characterize the multi-stream structure of dark matter pa… ▽ More Virialized halos of cold dark matter generically exhibit multi-stream structures of accreted dark matter within an outermost radial caustic known as the splashback radius. By tracking the particle trajectories that accrete onto the halos in cosmological $N$-body simulations, we count their number of apocenter passages ($p$), and use them to characterize the multi-stream structure of dark matter particles. We find that the radial density profile for each stream, classified by the number of apocenter passages, exhibits universal features, and can be described by a double power-law function comprising inner shallow and outer steep slopes of indices of $-1$ and $-8$, respectively. Surprisingly, these properties hold over a wide range of halo masses. The double-power law feature is persistent when dividing the sample by concentration or accretion rate. The dependence of the characteristic scale and amplitude of the profile on $p$ cannot be replicated by known self-similar solutions, requiring consideration of complexities such as the distribution of angular momentum or mergers. △ Less

Submitted 21 June, 2023; v1 submitted 2 February, 2023; originally announced February 2023.

Comments: 7 pages, 5 figures, updated version matches to publication in ApJL

Report number: YITP-22-163

Journal ref: ApJL 950 L13 (2023)

Cosmological gravity probes: connecting recent theoretical developments to forthcoming observations

Authors: Shun Arai, Katsuki Aoki, Yuji Chinone, Rampei Kimura, Tsutomu Kobayashi, Hironao Miyatake, Daisuke Yamauchi, Shuichiro Yokoyama, Kazuyuki Akitsu, Takashi Hiramatsu, Shin'ichi Hirano, Ryotaro Kase, Taishi Katsuragawa, Yosuke Kobayashi, Toshiya Namikawa, Takahiro Nishimichi, Teppei Okumura, Maresuke Shiraishi, Masato Shirasaki, Tomomi Sunayama, Kazufumi Takahashi, Atsushi Taruya, Junsei Tokuda

Abstract: Since the discovery of the accelerated expansion of the present Universe, significant theoretical developments have been made in the area of modified gravity. In the meantime, cosmological observations have been providing more high-quality data, allowing us to explore gravity on cosmological scales. To bridge the recent theoretical developments and observations, we present an overview of a variety… ▽ More Since the discovery of the accelerated expansion of the present Universe, significant theoretical developments have been made in the area of modified gravity. In the meantime, cosmological observations have been providing more high-quality data, allowing us to explore gravity on cosmological scales. To bridge the recent theoretical developments and observations, we present an overview of a variety of modified theories of gravity and the cosmological observables in the cosmic microwave background and large-scale structure, supplemented with a summary of predictions for cosmological observables derived from cosmological perturbations and sophisticated numerical studies. We specifically consider scalar-tensor theories in the Horndeski and DHOST family, massive gravity/bigravity, vector-tensor theories, metric-affine gravity, and cuscuton/minimally-modified gravity, and discuss the current status of those theories with emphasis on their physical motivations, validity, appealing features, the level of maturity, and calculability. We conclude that the Horndeski theory is one of the most well-developed theories of modified gravity, although several remaining issues are left for future observations. The paper aims to help to develop strategies for testing gravity with ongoing and forthcoming cosmological observations. △ Less

Submitted 18 December, 2022; originally announced December 2022.

Comments: 115 pages, 8 figures

KiDS-1000: Combined halo-model cosmology constraints from galaxy abundance, galaxy clustering and galaxy-galaxy lensing

Authors: Andrej Dvornik, Catherine Heymans, Marika Asgari, Constance Mahony, Benjamin Joachimi, Maciej Bilicki, Elisa Chisari, Hendrik Hildebrandt, Henk Hoekstra, Harry Johnston, Konrad Kuijken, Alexander Mead, Hironao Miyatake, Takahiro Nishimichi, Robert Reischke, Sandra Unruh, Angus H. Wright

Abstract: We present constraints on the flat $Λ$CDM cosmological model through a joint analysis of galaxy abundance, galaxy clustering and galaxy-galaxy lensing observables with the Kilo-Degree Survey. Our theoretical model combines a flexible conditional stellar mass function, to describe the galaxy-halo connection, with a cosmological N-body simulation-calibrated halo model to describe the non-linear matt… ▽ More We present constraints on the flat $Λ$CDM cosmological model through a joint analysis of galaxy abundance, galaxy clustering and galaxy-galaxy lensing observables with the Kilo-Degree Survey. Our theoretical model combines a flexible conditional stellar mass function, to describe the galaxy-halo connection, with a cosmological N-body simulation-calibrated halo model to describe the non-linear matter field. Our magnitude-limited bright galaxy sample combines 9-band optical-to-near-infrared photometry with an extensive and complete spectroscopic training sample to provide accurate redshift and stellar mass estimates. Our faint galaxy sample provides a background of accurately calibrated lensing measurements. We constrain the structure growth parameter $S_8=σ_8\sqrt{Ω_{\mathrm{m}}/0.3}=0.773^{+0.028}_{-0.030}$, and the matter density parameter $Ω_{\mathrm{m}}=0.290^{+0.021}_{-0.017}$. The galaxy-halo connection model adopted in the work is shown to be in agreement with previous studies. Our constraints on cosmological parameters are comparable to, and consistent with, joint $3\times2{\mathrm{pt}}$ clustering-lensing analyses that additionally include a cosmic shear observable. This analysis therefore brings attention to the significant constraining power in the often-excluded non-linear scales for galaxy clustering and galaxy-galaxy lensing observables. By adopting a theoretical model that accounts for non-linear halo bias, halo exclusion, scale-dependent galaxy bias and the impact of baryon feedback, this work demonstrates the potential and a way forward to include non-linear scales in cosmological analyses. Varying the width of the satellite galaxy distribution with an additional parameter yields a strong preference for sub-Poissonian variance, improving the goodness of fit by 0.18 in reduced $χ^{2}$ value compared to a fixed Poisson distribution. △ Less

Submitted 13 June, 2023; v1 submitted 6 October, 2022; originally announced October 2022.

Comments: 25 pages, 16 figures, accepted for publication in A&A

Journal ref: A&A 675, A189 (2023)

Perturbation Theory Remixed: Improved Nonlinearity Modeling beyond Standard Perturbation Theory

Authors: Zhenyuan Wang, Donghui Jeong, Atsushi Taruya, Takahiro Nishimichi, Ken Osato

Abstract: We present a novel $n$EPT ($n$th-order Eulerian Perturbation Theory) scheme to model the nonlinear density field by the summation up to $n$th-order density fields in perturbation theory. The obtained analytical power spectrum shows excellent agreement with the results from all 20 Dark-Quest suites of $N$-body simulations spreading over a broad range of cosmologies. The agreement is much better tha… ▽ More We present a novel $n$EPT ($n$th-order Eulerian Perturbation Theory) scheme to model the nonlinear density field by the summation up to $n$th-order density fields in perturbation theory. The obtained analytical power spectrum shows excellent agreement with the results from all 20 Dark-Quest suites of $N$-body simulations spreading over a broad range of cosmologies. The agreement is much better than the conventional two-loop Standard Perturbation Theory and would reach out to $k_{\rm max}\simeq 0.4~h/{\rm Mpc}$ at $z=3$ for the best-fitting Planck cosmology, without any free parameters. The method can accelerate the forward modeling of the non-linear cosmological density field, an indispensable probe of cosmic mysteries such as inflation, dark energy, and dark matter. △ Less

Submitted 31 August, 2022; originally announced September 2022.

Comments: 9 pages, 4 figures; to be submitted to PRL

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.

Cosmology with the Redshift-Space Galaxy Bispectrum Monopole at One-Loop Order

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

Abstract: We study the cosmological information content of the redshift-space galaxy bispectrum monopole at one-loop order in perturbation theory. We incorporate all effects necessary for comparison to data: fourth-order galaxy bias, infrared resummation (accounting for the non-linear evolution of baryon acoustic oscillations), ultraviolet counterterms, non-linear redshift-space distortions, stochastic cont… ▽ More We study the cosmological information content of the redshift-space galaxy bispectrum monopole at one-loop order in perturbation theory. We incorporate all effects necessary for comparison to data: fourth-order galaxy bias, infrared resummation (accounting for the non-linear evolution of baryon acoustic oscillations), ultraviolet counterterms, non-linear redshift-space distortions, stochastic contributions, projection, and binning effects. The model is implemented using FFTLog, and validated with the PT Challenge suite of $N$-body simulations, whose large volume allows for high-precision tests. Focusing on the mass fluctuation amplitude, $σ_8$, and galaxy bias parameters, we find that including one-loop corrections allow us to significantly extend the range of scales over which the bispectrum can be modeled, and greatly tightens constraints on bias parameters. However, this does not lead to noticeable improvements in the $σ_8$ errorbar due to the necessary marginalization over a large number of nuisance parameters with conservative priors. Analyzing a BOSS-volume likelihood, we find that the addition of the one-loop bispectrum may lead to improvements on primordial non-Gaussianity constraints by $\lesssim 30\%$ and on $σ_8$ by $\approx 10\%$, though we caution that this requires pushing the analysis to short scales where the galaxy bias parameters may not be correctly recovered; this may lead to biases in the recovered parameter values. We conclude that restrictive priors from simulations or higher-order statistics such as the bispectrum multipoles will be needed in order to realize the full information content of the galaxy bispectrum. △ Less

Submitted 13 June, 2022; v1 submitted 6 June, 2022; originally announced June 2022.

Comments: 51 pages, 6 figures, one-loop redshift-space galaxy bispectrum code is publicly available at https://github.com/oliverphilcox/OneLoopBispectrum

Report number: CERN-TH-2022-092, YITP-22-60

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

Dark matter halos of luminous AGNs from galaxy-galaxy lensing with the HSC Subaru Strategic Program

Authors: Wentao Luo, John D. Silverman, Surhud More, Andy Goulding, Hironao Miyatake, Takahiro Nishimichi, Chiaki Hikage, Lalitwadee Kawinwanichakij, Junyao Li, Xiangchong Li, Elinor Medezinski, Masamune Oguri, Taira Oogi, Cristobal Sifon

Abstract: We assess the dark matter halo masses of luminous AGNs over the redshift range 0.2 to 1.2 using galaxy-galaxy lensing based on imaging data from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). We measure the weak lensing signal of a sample of 8882 AGNs constructed using HSC and WISE photometry. The lensing detection around AGNs has a signal-to-noise ratio of 15. As expected, we find that… ▽ More We assess the dark matter halo masses of luminous AGNs over the redshift range 0.2 to 1.2 using galaxy-galaxy lensing based on imaging data from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). We measure the weak lensing signal of a sample of 8882 AGNs constructed using HSC and WISE photometry. The lensing detection around AGNs has a signal-to-noise ratio of 15. As expected, we find that the lensing mass profile is consistent with that of massive galaxies ($M_{*}\sim 10.8~M_\odot$). Surprisingly, the lensing signal remains unchanged when the AGN sample is split into low and high stellar mass hosts. Specifically, we find that the excess surface density (ESD) of AGNs, residing in galaxies with high stellar masses, significantly differs from that of the control sample. We further fit a halo occupation distribution model to the data to infer the posterior distribution of parameters including the average halo mass. We find that the characteristic halo mass of the full AGN population lies near the knee ($\rm log(M_h/h^{-1}M_{\odot})=12.0$) of the stellar-to-halo mass relation (SHMR). Illustrative of the results given above, the halo masses of AGNs residing in host galaxies with high stellar masses (i.e., above the knee of the SHMR) falls below the calibrated SHMR while the halo mass of the low stellar mass sample is more consistent with the established SHMR. These results indicate that massive halos with higher clustering bias tend to suppress AGN activity, probably due to the lack of available gas. △ Less

Submitted 11 April, 2022; v1 submitted 7 April, 2022; originally announced April 2022.

Comments: 17 pages,10 figures, accepted for publication in the Astrophysical Journal

Consistent lensing and clustering in a low-$S_8$ Universe with BOSS, DES Year 3, HSC Year 1 and KiDS-1000

Authors: A. Amon, N. C. Robertson, H. Miyatake, C. Heymans, M. White, J. DeRose, S. Yuan, R. H. Wechsler, T. N. Varga, S. Bocquet, A. Dvornik, S. More, A. J. Ross, H. Hoekstra, A. Alarcon, M. Asgari, J. Blazek, A. Campos, R. Chen, A. Choi, M. Crocce, H. T. Diehl, C. Doux, K. Eckert, J. Elvin-Poole , et al. (83 additional authors not shown)

Abstract: We evaluate the consistency between lensing and clustering probes of large-scale structure based on measurements of projected galaxy clustering from BOSS combined with overlapping galaxy-galaxy lensing from three surveys: DES Y3, HSC Y1, and KiDS-1000. An intra-lensing-survey study finds good agreement between these lensing data. We model the observations using the Dark Emulator and fit the data a… ▽ More We evaluate the consistency between lensing and clustering probes of large-scale structure based on measurements of projected galaxy clustering from BOSS combined with overlapping galaxy-galaxy lensing from three surveys: DES Y3, HSC Y1, and KiDS-1000. An intra-lensing-survey study finds good agreement between these lensing data. We model the observations using the Dark Emulator and fit the data at two fixed cosmologies: Planck, with $S_8=0.83$, and a Lensing cosmology with $S_8=0.76$. For a joint analysis limited to scales with $R>5.25h^{-1}$Mpc, we find that both cosmologies provide an acceptable fit to the data. Full utilisation of the small-scale clustering and lensing measurements is hindered by uncertainty in the impact of baryon feedback and assembly bias, which we account for with a reasoned theoretical error budget. We incorporate a systematic scaling parameter for each redshift bin, $A$, that decouples the lensing and clustering to capture any inconsistency. When a wide range of scales ($0.15<R<60h^{-1}$Mpc) are incorporated, we find different results for the consistency of clustering and lensing between the two cosmologies. Limiting the analysis to the bins for which the impact of the selection of the lens sample is expected to be minimal, for the low-$S_8$ Lensing cosmology, the measurements are consistent with $A$=1; $A=0.91\pm0.04$ using DES+KiDS and $A=0.97\pm0.06$ using HSC. For the Planck cosmology case, we find a discrepancy: $A=0.79\pm0.03$ using DES+KiDS and $A=0.84\pm0.05$ using HSC. We demonstrate that a kSZ-based estimate for baryonic effects alleviates some of the discrepancy in the Planck cosmology. This analysis demonstrates the statistical power of these small-scale measurements, but also indicates that caution is still warranted given current uncertainties in modelling baryonic effects, assembly bias, and selection effects in the foreground sample. △ Less

Submitted 13 October, 2022; v1 submitted 15 February, 2022; originally announced February 2022.

Comments: 28 pages, 11 figures

The halo model with beyond-linear halo bias: unbiasing cosmological constraints from galaxy-galaxy lensing and clustering

Authors: Constance Mahony, Andrej Dvornik, Alexander Mead, Catherine Heymans, Marika Asgari, Hendrik Hildebrandt, Hironao Miyatake, Takahiro Nishimichi, Robert Reischke

Abstract: We determine the error introduced in a joint halo model analysis of galaxy-galaxy lensing and galaxy clustering observables when adopting the standard approximation of linear halo bias. Considering the Kilo-Degree Survey, we forecast that ignoring the non-linear halo bias would result in up to 5$σ$ offsets in the recovered cosmological parameters describing structure growth, $S_8$, and the matter… ▽ More We determine the error introduced in a joint halo model analysis of galaxy-galaxy lensing and galaxy clustering observables when adopting the standard approximation of linear halo bias. Considering the Kilo-Degree Survey, we forecast that ignoring the non-linear halo bias would result in up to 5$σ$ offsets in the recovered cosmological parameters describing structure growth, $S_8$, and the matter density parameter, $Ω_{\mathrm{m}}$. We include the scales $10^{-1.3}<r_{\rm{p}} \ / h^{-1}\, \mathrm{Mpc}<10$ in the data vector, and the direction of these offsets are shown to depend on the freedom afforded to the halo model through other nuisance parameters. We conclude that a beyond-linear halo bias correction must therefore be included in future cosmological halo model analyses of large-scale structure observables on non-linear scales. △ Less

Submitted 30 June, 2022; v1 submitted 3 February, 2022; originally announced February 2022.

Comments: v2: Version accepted in MNRAS. Minor clarifications and added robustness test, conclusions unchanged

Report number: YITP-22-07

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

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

Taming redshift-space distortion effects in the EFTofLSS and its application to data

Authors: Guido D'Amico, Leonardo Senatore, Pierre Zhang, Takahiro Nishimichi

Abstract: Former analyses of the BOSS data using the Effective Field Theory of Large-Scale Structure (EFTofLSS) have measured that the largest counterterms are the redshift-space distortion ones. This allows us to adjust the power-counting rules of the theory, and to explicitly identify that the leading next-order terms have a specific dependence on the cosine of the angle between the line-of-sight and the… ▽ More Former analyses of the BOSS data using the Effective Field Theory of Large-Scale Structure (EFTofLSS) have measured that the largest counterterms are the redshift-space distortion ones. This allows us to adjust the power-counting rules of the theory, and to explicitly identify that the leading next-order terms have a specific dependence on the cosine of the angle between the line-of-sight and the wavenumber of the observable, $μ$. Such a specific $μ$-dependence allows us to construct a linear combination of the data multipoles, $\slashed{P}$, where these contributions are effectively projected out, so that EFTofLSS predictions for $\slashed{P}$ have a much smaller theoretical error and so a much higher $k$-reach. The remaining data are organized in wedges in $μ$ space, have a $μ$-dependent $k$-reach because they are not equally affected by the leading next-order contributions, and therefore can have a higher $k$-reach than the multipoles. Furthermore, by explicitly including the highest next-order terms, we define a `one-loop+' procedure, where the wedges have even higher $k$-reach. We study the effectiveness of these two procedures on several sets of simulations and on the BOSS data. The resulting analysis has identical computational cost as the multipole-based one, but leads to an improvement on the determination of some of the cosmological parameters that ranges from $10\%$ to $100\%$, depending on the survey properties. △ Less

Submitted 30 September, 2021; originally announced October 2021.

Comments: 26 pages, 6 figures

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)

Grid-based calculations of redshift-space matter fluctuations from perturbation theory: UV sensitivity and convergence at the field level

Authors: Atsushi Taruya, Takahiro Nishimichi, Donghui Jeong

Abstract: Perturbation theory (PT) has been used to interpret the observed nonlinear large-scale structure statistics at the quasi-linear regime. To facilitate the PT-based analysis, we have presented the GridSPT algorithm, a grid-based method to compute the nonlinear density and velocity fields in standard perturbation theory (SPT) from a given linear power spectrum. Here, we further put forward the approa… ▽ More Perturbation theory (PT) has been used to interpret the observed nonlinear large-scale structure statistics at the quasi-linear regime. To facilitate the PT-based analysis, we have presented the GridSPT algorithm, a grid-based method to compute the nonlinear density and velocity fields in standard perturbation theory (SPT) from a given linear power spectrum. Here, we further put forward the approach by taking the redshift-space distortions into account. With the new implementation, we have, for the first time, generated the redshift-space density field to the fifth order and computed the next-to-next-to-leading order (2 loop) power spectrum and the next-to-leading order (1 loop) bispectrum of matter clustering in redshift space. By comparing the result with corresponding analytical SPT calculation and $N$-body simulations, we find that the SPT calculation (A) suffers much more from the UV sensitivity due to the higher-derivative operators and (B) deviates from the $N$-body results from the Fourier wavenumber smaller than real space $k_{\rm max}$. Finally, we have shown that while Padé approximation removes spurious features in morphology, it does not improve the modeling of power spectrum and bispectrum. △ Less

Submitted 14 September, 2021; originally announced September 2021.

Comments: 23 pages, 18 figures

Report number: YITP-21-95

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

Implementing spectra response function approaches for fast calculation of power spectra and bispectra

Authors: Ken Osato, Takahiro Nishimichi, Atsushi Taruya, Francis Bernardeau

Abstract: Perturbation theory of large-scale structures of the Universe at next-to-leading order and next-to-next-to-leading order provides us with predictions of cosmological statistics at sub-percent level in the mildly non-linear regime. Its use to infer cosmological parameters from spectroscopic surveys, however, is hampered by the computational cost of making predictions for a large number of parameter… ▽ More Perturbation theory of large-scale structures of the Universe at next-to-leading order and next-to-next-to-leading order provides us with predictions of cosmological statistics at sub-percent level in the mildly non-linear regime. Its use to infer cosmological parameters from spectroscopic surveys, however, is hampered by the computational cost of making predictions for a large number of parameters. In order to reduce the running time of the codes, we present a fast scheme in the context of the regularized perturbation theory approach and applied it to power spectra at 2-loop level and bispectra at 1-loop level, including the impact of binning. This method utilizes a Taylor expansion of the power spectrum as a functional of the linear power spectrum around fiducial points at which costly direct evaluation of perturbative diagrams is performed and tabulated. The computation of the predicted spectra for arbitrary cosmological parameters then requires only one-dimensional integrals that can be done within a few minutes. It makes this method suitable for Markov chain Monte-Carlo analyses for cosmological parameter inference. △ Less

Submitted 29 October, 2021; v1 submitted 9 July, 2021; originally announced July 2021.

Comments: 32 pages, 21 figures, PRD in press

Report number: YITP-21-71

Journal ref: Phys. Rev. D 104, 103501 (2021)

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

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

Covariance of the matter power spectrum including the survey window function effect: N-body simulations vs. fifth-order perturbation theory on grid

Authors: Atsushi Taruya, Takahiro Nishimichi, Donghui Jeong

Abstract: We present a Next-to-next-to-leading (fifth or NNLO) order calculation for the covariance matrix of the matter power spectrum, taking into account the effect of survey window functions. Using the grid-based calculation scheme for the standard perturbation theory, GridSPT, we quickly generate multiple realizations of the nonlinear density fields to fifth order in perturbation theory, then estimate… ▽ More We present a Next-to-next-to-leading (fifth or NNLO) order calculation for the covariance matrix of the matter power spectrum, taking into account the effect of survey window functions. Using the grid-based calculation scheme for the standard perturbation theory, GridSPT, we quickly generate multiple realizations of the nonlinear density fields to fifth order in perturbation theory, then estimate the power spectrum and the covariance matrix from the sample. To the end, we have obtained the non-Gaussian covariance originated from the one-loop trispectrum without explicitly computing the trispectrum. By comparing the GridSPT calculations with the N-body results, we show that NNLO GridSPT result reproduces the N-body results on quasi-linear scales where SPT accurately models nonlinear matter power spectrum. Incorporating the survey window function effect to GridSPT is rather straightforward, and the resulting NNLO covariance matrix also matches well with the N-body results. △ Less

Submitted 6 January, 2021; v1 submitted 10 July, 2020; originally announced July 2020.

Comments: 20 pages, 16 figures, updated to match the version published in PRD

Report number: YITP-20-88

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

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)

Shapes and alignments of dark matter haloes and their brightest cluster galaxies in 39 strong lensing clusters

Authors: Taizo Okabe, Masamune Oguri, Sébastien Peirani, Yasushi Suto, Yohan Dubois, Christophe Pichon, Tetsu Kitayama, Shin Sasaki, Takahiro Nishimichi

Abstract: We study shapes and alignments of 45 dark matter (DM) haloes and their brightest cluster galaxies (BCGs) using a sample of 39 massive clusters from Hubble Frontier Field (HFF), Cluster Lensing And Supernova survey with Hubble (CLASH), and Reionization Lensing Cluster Survey (RELICS). We measure shapes of the DM haloes by strong gravitational lensing, whereas BCG shapes are derived from their light… ▽ More We study shapes and alignments of 45 dark matter (DM) haloes and their brightest cluster galaxies (BCGs) using a sample of 39 massive clusters from Hubble Frontier Field (HFF), Cluster Lensing And Supernova survey with Hubble (CLASH), and Reionization Lensing Cluster Survey (RELICS). We measure shapes of the DM haloes by strong gravitational lensing, whereas BCG shapes are derived from their light profiles in Hubble Space Telescope images. Our measurements from a large sample of massive clusters presented here provide new constraints on dark matter and cluster astrophysics. We find that DM haloes are on average highly elongated with the mean ellipticity of $0.482\pm 0.028$, and position angles of major axes of DM haloes and their BCGs tend to be aligned well with the mean value of alignment angles of $22.2\pm 3.9$ deg. We find that DM haloes in our sample are on average more elongated than their BCGs with the mean difference of their ellipticities of $0.11\pm 0.03$. In contrast, the Horizon-AGN cosmological hydrodynamical simulation predicts on average similar ellipticities between DM haloes and their central galaxies. While such a difference between the observations and the simulation may well be explained by the difference of their halo mass scales, other possibilities include the bias inherent to strong lensing measurements, limited knowledge of baryon physics, or a limitation of cold dark matter. △ Less

Submitted 23 May, 2020; originally announced May 2020.

Comments: 18 pages, 16 figures, 6 Tables, accepted by MNRAS

Report number: YITP-20-76

Model independent measurement of the growth rate from the consistency relations of the LSS

Authors: Marco Marinucci, Takahiro Nishimichi, Massimo Pietroni

Abstract: The Consistency Relations for the Large Scale Structure provide a link between the amplitude of Baryonic Acoustic Oscillations in the squeezed bispectrum (BS) and in the power spectrum (PS). This relation depends on the large scale bias of the considered tracer, $b_α$, and on the growth rate of structures, $f$. Remarkably, originating from basic symmetry principles, this relation is exact and inde… ▽ More The Consistency Relations for the Large Scale Structure provide a link between the amplitude of Baryonic Acoustic Oscillations in the squeezed bispectrum (BS) and in the power spectrum (PS). This relation depends on the large scale bias of the considered tracer, $b_α$, and on the growth rate of structures, $f$. Remarkably, originating from basic symmetry principles, this relation is exact and independent on the underlying cosmological model. By analysing data from large volume simulations, both for dark matter and for haloes, we illustrate how BS and PS measurements can be used to extract $b_α$ and $f$ without the need of any theoretical approximation scheme for the computation of the BS and the PS. We show that, combining measurements of the squeezed BS with the quadrupole to monopole ratios for the PS at large scales can successfully break the $b_α-f$ degeneracy. We forecast that this method, applied to a Euclid-like survey, will be able to measure bias, and then the growth rate, at better than $10\%$ level, with no extra assumption. △ Less

Submitted 19 May, 2020; originally announced May 2020.

Comments: 22 pages, 7 figures

Report number: YITP-20-49

Accurate emulator for the redshift-space power spectrum of dark matter halos and its application to galaxy power spectrum

Authors: Yosuke Kobayashi, Takahiro Nishimichi, Masahiro Takada, Ryuichi Takahashi, Ken Osato

Abstract: An accurate theoretical template of the redshift-space galaxy power spectrum, if applicable out to nonlinear scales, enables us to extract more stringent and robust constraints on cosmological parameters from the measured galaxy clustering. In this work, we develop a simulation-based template, so-called emulator, for the redshift-space power spectrum of dark matter halos. Using the redshift-space… ▽ More An accurate theoretical template of the redshift-space galaxy power spectrum, if applicable out to nonlinear scales, enables us to extract more stringent and robust constraints on cosmological parameters from the measured galaxy clustering. In this work, we develop a simulation-based template, so-called emulator, for the redshift-space power spectrum of dark matter halos. Using the redshift-space halo power spectra measured from the Dark Quest $N$-body simulation suite that covers 101 flat-geometry $w$CDM cosmologies around the Planck $Λ$CDM model, we feed these data into a feed-forward neural network to build the fast and accurate emulation of the power spectrum from the linear to nonlinear scales up to $k \simeq 0.6 \,h \,{\rm Mpc}^{-1}$. Our emulator achieves about 1% and 5% fractional accuracies in predicting the monopole and quadrupole moments of the power spectrum, respectively, for halos of $\sim 10^{13}h^{-1}M_\odot$ that correspond to host halos of the SDSS LOWZ- and CMASS-like galaxies, where the achieved accuracies are sufficient compared to the statistical errors of SDSS volume. The validation and performance of the emulator are given by the comparison of the emulator predictions with the power spectra directly measured from the simulations for validation sets that are not used in the training. We demonstrate that the emulator outputs can be used to make model predictions for the redshift-space power spectrum of galaxies by employing user-fed models for the halo-galaxy connection, such as the halo occupation distribution. The emulator allows us to easily incorporate the Finger-of-God effect due to the virial motions of galaxies and the Alcock-Paczyński distortions. Our code can compute the redshift-space galaxy power spectrum in a CPU subseconds and is ready to perform the emulator-based cosmological analysis for the exiting and upcoming galaxy redshift surveys. △ Less

Submitted 18 October, 2020; v1 submitted 12 May, 2020; originally announced May 2020.

Comments: 26 pages, 19 figures, published by Phys. Rev. D

Report number: IPMU20-0046, YITP-20-41

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

Lagrangian-space Gaussian ansatz for the matter redshift-space power spectrum and correlation function

Authors: Patrick Valageas, Takahiro Nishimichi

Abstract: We study the predictions for the matter redshift-space power spectrum and correlation function of a Lagrangian-space Gaussian ansatz introduced in a previous work. This model is a natural extension of the Zeldovich approximation, where the displacement and velocity power spectra are determined by the equations of motion, instead of being set equal to the linear power spectrum. It does not contain… ▽ More We study the predictions for the matter redshift-space power spectrum and correlation function of a Lagrangian-space Gaussian ansatz introduced in a previous work. This model is a natural extension of the Zeldovich approximation, where the displacement and velocity power spectra are determined by the equations of motion, instead of being set equal to the linear power spectrum. It does not contain any free parameter. As for the real-space statistics, we find that this Lagrangian-space approach is much more efficient for the correlation functions than for the power spectra. The damping of the BAO oscillations is well recovered but there is a large smooth drift from the simulations in the power spectra. The multipoles of the correlation functions are well recovered on BAO scales, with an accuracy of $2\%$ for $ξ^s_0$ down to $10 h^{-1}$ Mpc, and of $3\%$ for $ξ^s_2$ down to $26 h^{-1}$ Mpc, at $z \geq 0.35$. △ Less

Submitted 29 April, 2020; originally announced April 2020.

Comments: 15 pages

Report number: YITP-20-50

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

Power spectrum of halo intrinsic alignments in simulations

Authors: Toshiki Kurita, Masahiro Takada, Takahiro Nishimichi, Ryuichi Takahashi, Ken Osato, Yosuke Kobayashi

Abstract: We use a suite of $N$-body simulations to study intrinsic alignments (IA) of halo shapes with the surrounding large-scale structure in the $Λ$CDM model. For this purpose, we develop a novel method to measure multipole moments of the three-dimensional power spectrum of the $E$-mode field of halo shapes with the matter/halo distribution, $P_{δE}^{(\ell)}(k)$ (or $P^{(\ell)}_{{\rm h}E}$), and those o… ▽ More We use a suite of $N$-body simulations to study intrinsic alignments (IA) of halo shapes with the surrounding large-scale structure in the $Λ$CDM model. For this purpose, we develop a novel method to measure multipole moments of the three-dimensional power spectrum of the $E$-mode field of halo shapes with the matter/halo distribution, $P_{δE}^{(\ell)}(k)$ (or $P^{(\ell)}_{{\rm h}E}$), and those of the auto-power spectrum of the $E$ mode, $P^{(\ell)}_{EE}(k)$, based on the $E$/$B$-mode decomposition. The IA power spectra have non-vanishing amplitudes over the linear to nonlinear scales, and the large-scale amplitudes at $k\lesssim 0.1~h~{\rm Mpc}^{-1}$ are related to the matter power spectrum via a constant coefficient ($A_{\rm IA}$), similar to the linear bias parameter of galaxy or halo density field. We find that the cross- and auto-power spectra $P_{δE}$ and $P_{EE}$ at nonlinear scales, $k\gtrsim 0.1~h~{\rm Mpc}^{-1}$, show different $k$-dependences relative to the matter power spectrum, suggesting a violation of the nonlinear alignment model commonly used to model contaminations of cosmic shear signals. The IA power spectra exhibit baryon acoustic oscillations, and vary with halo samples of different masses, redshifts and cosmological parameters ($Ω_{\rm m}, S_8$). The cumulative signal-to-noise ratio for the IA power spectra is about 60% of that for the halo density power spectrum, where the super-sample covariance is found to give a significant contribution to the total covariance. Thus our results demonstrate that the IA power spectra of galaxy shapes, measured from imaging and spectroscopic surveys for an overlapping area of the sky, can be used to probe the underlying matter power spectrum, the primordial curvature perturbations, and cosmological parameters, in addition to the standard galaxy density power spectrum. △ Less

Submitted 19 November, 2020; v1 submitted 27 April, 2020; originally announced April 2020.

Comments: 20 pages, 21 figures; updated to reflect accepted version

Report number: IPMU20-0044, YITP-20-48

Observing Baryonic Acoustic Oscillations in tomographic cosmic shear surveys

Authors: Francis Bernardeau, Takahiro Nishimichi, Atsushi Taruya

Abstract: We show that it is possible to build effective matter density power spectra in tomographic cosmic shear observations that exhibit the Baryonic Acoustic Oscillations (BAO) features once a nulling transformation has been applied to the data. The precision with which the amplitude and position of these features can be reconstructed is quantified in terms of sky coverage, intrinsic shape noise, median… ▽ More We show that it is possible to build effective matter density power spectra in tomographic cosmic shear observations that exhibit the Baryonic Acoustic Oscillations (BAO) features once a nulling transformation has been applied to the data. The precision with which the amplitude and position of these features can be reconstructed is quantified in terms of sky coverage, intrinsic shape noise, median source redshift and number density of sources. BAO detection in Euclid or LSST like wide surveys will be possible with a modest signal-to-noise ratio. It would improve dramatically for slightly deeper surveys. △ Less

Submitted 7 April, 2020; originally announced April 2020.

Comments: 5 pages, 4 figures

Report number: YITP-20-15