Large scale structure surveys promise to be the next leading probe of cosmological information. It is therefore crucial to reliably predict their observables. The Effective Field Theory of Large Scale Structures (EFTofLSS) provides a manifestly convergent perturbation theory for the weakly non-linear regime of dark matter, where correlation functions are computed in an expansion of the wavenumber k of a mode over the wavenumber associated with the non-linear scale k NL . Since most of the information is contained at high wavenumbers, it is necessary to compute higher order corrections to correlation functions. After the one-loop correction to the matter power spectrum, we estimate that the next leading one is the two-loop contribution, which we compute here. At this order in k/k NL , there is only one counterterm in the EFTofLSS that must be included, though this term contributes both at tree-level and in several one-loop diagrams. We also discuss correlation functions involving the velocity and momentum fields. We find that the EFTofLSS prediction at two loops matches to percent accuracy the non-linear matter power spectrum at redshift zero up to k ∼ 0.6 h Mpc −1 , requiring just one unknown coefficient that needs to be fit to observations. Given that Standard Perturbation Theory stops converging at redshift zero at k ∼ 0.1 h Mpc −1 , our results demonstrate the possibility of accessing a factor of order 200 more dark matter quasi-linear modes than naively expected. If the remaining observational challenges to accessing these modes can be addressed with similar success, our results show that there is tremendous potential for large scale structure surveys to explore the primordial universe.

Large scale structure surveys are likely the next leading probe of cosmological information. It is therefore crucial to reliably predict their observables. The Effective Field Theory of Large Scale Structures (EFTofLSS) provides a manifestly convergent perturbation theory for the weakly non-linear regime, where dark matter correlation functions are computed in an expansion of the wavenumber k over the wavenumber associated to the non-linear scale kNL. To push the predictions to higher wavenumbers, it is necessary to compute the 2-loop matter power spectrum. For equal-time correlators, exactly as with standard perturturbation theory, there are IR divergences present in each diagram that cancel completely in the final result. We develop a method by which all 2-loop diagrams are computed as one integral, with an integrand that is manifestly free of any IR divergences. This allows us to compute the 2-loop power spectra in a reliable way that is much less numerically challenging than standard techniques. We apply our method to scaling universes where the linear power spectrum is a single power law of k, and where IR divergences can particularly easily interfere with accurate evaluation of loop corrections if not handled carefully. We show that our results are independent of IR cutoff and, after renormalization, of the UV cutoff, and comment how the method presented here naturally generalizes to higher loops.

Given the importance of future large scale structure surveys for delivering new cosmological information, it is crucial to reliably predict their observables. The Effective Field Theory of Large Scale Structures (EFTofLSS) provides a manifestly convergent perturbative scheme to compute the clustering of dark matter in the weakly nonlinear regime in an expansion in k/k NL , where k is the wavenumber of interest and k NL is the wavenumber associated to the nonlinear scale. It has been recently shown that the EFTofLSS matches to 1% level the dark matter power spectrum at redshift zero up to k 0.3h Mpc −1 and k 0.6h Mpc −1 at one and two loops respectively, using only one counterterm that is fit to data. Similar results have been obtained for the momentum power spectrum at one loop. This is a remarkable improvement with respect to former analytical techniques. Here we study the prediction for the equal-time dark matter bispectrum at one loop. We find that at this order it is sufficient to consider the same counterterm that was measured in the power spectrum. Without any remaining free parameter, and in a cosmology for which k NL is smaller than in the previously considered cases (σ 8 = 0.9), we find that the prediction from the EFTofLSS agrees very well with N -body simulations up to k 0.25h Mpc −1 , given the accuracy of the measurements, which is of order a few percent at the highest k's of interest. While the fit is very good on average up to k 0.25h Mpc −1 , the fit performs slightly worse on equilateral configurations, in agreement with expectations that for a given maximum k, equilateral triangles are the most nonlinear.