We study the Effective Field Theory of Large Scale Structure for cosmic density and momentum fields. We show that the finite part of the two-loop calculation and its counterterms introduce an apparent scale dependence for the leading order parameter $c_\text{s}^2$ of the EFT starting at k=0.1 h/Mpc. These terms limit the range over which one can trust the one-loop EFT calculation at the 1 % level to k<0.1 h/Mpc at redshift z=0. We construct a well motivated one parameter ansatz to fix the relative size of the one- and two-loop counterterms using their high-k sensitivity. Although this one parameter model is a very restrictive choice for the counterterms, it explains the apparent scale dependence of $c_\text{s}^2$ seen in simulations. It is also able to capture the scale dependence of the density power spectrum up to k$\approx$ 0.3 h/Mpc at the 1 % level at redshift $z=0$. Considering a simple scheme for the resummation of large scale motions, we find that the two loop calculation reduces the need for this IR-resummation at k<0.2 h/Mpc. Finally, we extend our calculation to momentum statistics and show that the same one parameter model can also describe density-momentum and momentum-momentum statistics.

We study the bispectrum in the Effective Field Theory of Large Scale Structure, consistently accounting for the effects of short-scale dynamics. We begin by proving that, as long as the theory is perturbative, it can be formulated to arbitrary order using only operators that are local in time. We then derive all the new operators required to cancel the UV-divergences and obtain a physically meaningful prediction for the one-loop bispectrum. In addition to new, subleading stochastic noises and the viscosity term needed for the one-loop power spectrum, we find three new effective operators. The three new parameters can be constrained by comparing with N-body simulations. The best fit is precisely what is suggested by the structure of UV-divergences, hence justifying a formula for the EFTofLSS bispectrum whose only fitting parameter is already fixed by the power spectrum. This result predicts the bispectrum of N-body simulations up to $k \approx 0.22\, h\, \text{Mpc}^{-1}$ at $z=0$, an improvement by nearly a factor of two as compared to one-loop standard perturbation theory.

We study the effective Lagrangian, at leading order in derivatives, that describes the propagation of density and metric fluctuations in a fluid composed by an arbitrary number of interacting components. Our results can be applied to any situation in cosmology where different species have non-gravitational interactions. In time-dependent backgrounds, such as FLRW, the quadratic action can only be diagonalized at fixed time slices and flavour mixing is unavoidable as time flows. In spite of this, the fluid can be interpreted at the level of the linear equations of motion as an ensemble of individual interacting species. We show that interactions lead to anisotropic stresses that depend on the mixing terms in the quadratic action for perturbations. In addition to the standard entrainment among different components, we find a new operator in the effective action that behaves as a cosmological constant when it dominates the dynamics.

We compute the renormalized two-point functions of density, divergence and vorticity of the velocity in the Effective Field Theory of Large Scale Structures. Because of momentum and mass conservation, the corrections from short scales to the large-scale power spectra of density, divergence and vorticity must start at order $k^{4}$. For the vorticity this constitutes one of the two leading terms. Exact (approximated) self-similarity of an Einstein-de Sitter ($\Lambda$CDM) background fixes the time dependence so that the vorticity power spectrum at leading order is determined by the symmetries of the problem and the power spectrum around the non-linear scale. We show that to cancel all divergences in the velocity correlators one needs new counterterms. These fix the definition of velocity and do not represent new properties of the system. For an Einstein-de Sitter universe, we show that all three renormalized cross- and auto-correlation functions have the same structure but different numerical coefficients, which we compute. We elucidate the differences between using momentum and velocity.