TY - JOUR

T1 - Grid-based calculation for perturbation theory of large-scale structure

AU - Taruya, Atsushi

AU - Nishimichi, Takahiro

AU - Jeong, Donghui

N1 - Funding Information:
This work was supported in part by MEXT/JSPS KAKENHI Grants No. JP15H05899, No. JP16H03977 (A. T.), and No. JP17K14273 (T. N.). T. N. also acknowledges financial support from Japan Science and Technology Agency (JST) CREST Grant No. JPMJCR1414. D. J. acknowledges support from NSF Grant (AST-1517363) and NASA 80NSSC18K1103. Numerical computation was partly carried out at the Yukawa Institute Computer Facility.

PY - 2018

Y1 - 2018

N2 - Perturbation theory calculation of large-scale structures has been used to interpret the observed nonlinear statistics of large-scale structures at the quasilinear regime. In particular, the so-called standard perturbation theory (SPT) provides a basis for the analytical computation of the higher-order quantities of large-scale structures. Here, we present a novel grid-based algorithm for the SPT calculation, hence named GridSPT, to generate the higher-order density and velocity fields from a given linear power spectrum. Taking advantage of the fast Fourier transform, the GridSPT quickly generates the nonlinear density fields at each order, from which we calculate the statistical quantities such as nonlinear power spectrum and bispectrum. Comparing the density fields (to fifth order) from GridSPT with those from the full N-body simulations in the configuration space, we find that GridSPT accurately reproduces the N-body result on large scales. The agreement worsens with smaller smoothing radius, particularly for the underdense regions where we find that the second-order Lagrangian perturbation theory algorithm performs well.

AB - Perturbation theory calculation of large-scale structures has been used to interpret the observed nonlinear statistics of large-scale structures at the quasilinear regime. In particular, the so-called standard perturbation theory (SPT) provides a basis for the analytical computation of the higher-order quantities of large-scale structures. Here, we present a novel grid-based algorithm for the SPT calculation, hence named GridSPT, to generate the higher-order density and velocity fields from a given linear power spectrum. Taking advantage of the fast Fourier transform, the GridSPT quickly generates the nonlinear density fields at each order, from which we calculate the statistical quantities such as nonlinear power spectrum and bispectrum. Comparing the density fields (to fifth order) from GridSPT with those from the full N-body simulations in the configuration space, we find that GridSPT accurately reproduces the N-body result on large scales. The agreement worsens with smaller smoothing radius, particularly for the underdense regions where we find that the second-order Lagrangian perturbation theory algorithm performs well.

UR - http://www.scopus.com/inward/record.url?scp=85057847899&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85057847899&partnerID=8YFLogxK

U2 - 10.1103/PhysRevD.98.103532

DO - 10.1103/PhysRevD.98.103532

M3 - Article

AN - SCOPUS:85057847899

VL - 98

JO - Physical Review D

JF - Physical Review D

SN - 2470-0010

IS - 10

M1 - 103532

ER -