Perturbation theory remixed. II. Improved modeling of nonlinear bispectrum

We present the application of the nth order Eulerian perturbation theory (nEPT) for modeling the matter bispectrum in real space as an advancement over the standard perturbation theory (SPT). The nEPT method, detailed in Wang et al. [Phys. Rev. D 107, 103534 (2023)PRVDAQ2470-001010.1103/PhysRevD.107.103534], sums up the density perturbations up to the nth order before computing summary statistics such as bispectrum. Taking advantage of grid-based calculation of SPT (GridSPT), we make a realization-based comparison of the analytical nonlinear bispectrum predictions from nEPT and SPT against a suite of N-body simulations. Using a spherical-bispectrum visualization scheme, we show that nEPT bispectrum matches better than SPT bispectrum over a wide range of scales in general wCDM cosmologies. Like the power spectrum case, we find that nEPT bispectrum modeling accuracy is controlled by σ8(z)σ8D(z), where D(z) is the linear growth factor at a redshift z. Notably, the 6EPT doubles the bispectrum model's validity range compared to the one-loop SPT for σ8(z)<0.5, corresponding to redshifts z≥1 for the best-fitting Planck-2018 cosmology. For n≥5, however, nEPT bispectrum depends sensitively on the cut-off scale or the grid resolution. The percent-level modeling accuracy achieved for the spherical bispectrum (where we average over all triangular configurations) becomes much degraded when fixing configurations. Thus, we show that the validity range of the field-level cosmological inferences must be different from that derived from averaged summary statistics such as n-point correlation functions.