In this paper we test the perturbative halo bias model at the field level. The advantage of this approach is that any analysis can be done without sample variance if the same initial conditions are used in simulations and perturbation theory calculations. We write the bias expansion in terms of modified bias operators in Eulerian space, designed such that the large bulk flows are automatically resummed and not treated perturbatively. Using these operators, the bias model accurately matches the Eulerian density of halos in N-body simulations. The mean-square model error is close to the Poisson shot noise for a wide range of halo masses and it is rather scale-independent, with scale-dependent corrections becoming relevant at the nonlinear scale. In contrast, for linear bias the mean-square model error can be higher than the Poisson prediction by factors of up to a few on large scales, and it becomes scale dependent already in the linear regime. We show that by weighting simulated halos by their mass, the mean-square error of the model can be further reduced by up to an order of magnitude, or by a factor of two when including $60\%$ mass scatter. We also test the Standard Eulerian bias model using the nonlinear matter field measured from simulations and show that it leads to a larger and more scale-dependent model error than the bias expansion based on perturbation theory. These results may be of particular relevance for cosmological inference methods that use a likelihood of the biased tracer at the field level, or for initial condition and BAO reconstruction that requires a precise estimate of the large-scale potential from the biased tracer density.

Angular statistics of cosmological observables are hard to compute. The main difficulty is due to the presence of highly-oscillatory Bessel functions which need to be integrated over. In this paper, we provide a simple and fast method to compute the angular power spectrum and bispectrum of any observable. The method is based on using an FFTlog algorithm to decompose the momentum-space statistics onto a basis of power-law functions. For each power law, the integrals over Bessel functions have a simple analytical solution. This allows us to efficiently evaluate these integrals, independently of the value of the multipole $\ell$. In particular, this method significantly speeds up the evaluation of the angular bispectrum compared to existing methods. To illustrate our algorithm, we compute the galaxy, lensing and CMB temperature angular power spectrum and bispectrum.

We present a systematic study of galaxy biasing in the presence of primordial non-Gaussianity. For a large class of non-Gaussian initial conditions, we define a general bias expansion and prove that it is closed under renormalization, thereby showing that the basis of operators in the expansion is complete. We then study the effects of primordial non-Gaussianity on the statistics of galaxies. We show that the equivalence principle enforces a relation between the scale-dependent bias in the galaxy power spectrum and that in the dipolar part of the bispectrum. This provides a powerful consistency check to confirm the primordial origin of any observed scale-dependent bias. Finally, we also discuss the imprints of anisotropic non-Gaussianity as motivated by recent studies of higher-spin fields during inflation.

We develop the effective theory of large-scale structure for non-Gaussian initial conditions. The effective stress tensor in the dark matter equations of motion contains new operators, which originate from the squeezed limit of the primordial bispectrum. Parameterizing the squeezed limit by a scaling and an angular dependence, captures large classes of primordial non-Gaussianity. Within this parameterization, we classify the possible contributions to the effective theory. We show explicitly how all terms consistent with the symmetries arise from coarse graining the dark matter equations of motion and its initial conditions. We also demonstrate that the system is closed under renormalization and that the basis of correction terms is therefore complete. The relevant corrections to the matter power spectrum and bispectrum are computed numerically and their relative importance is discussed.

The statistics of primordial curvature fluctuations are our window into the period of inflation, where these fluctuations were generated. To date, the cosmic microwave background has been the dominant source of information about these perturbations. Large scale structure is however from where drastic improvements should originate. In this paper, we explain the theoretical motivations for pursuing such measurements and the challenges that lie ahead. In particular, we discuss and identify theoretical targets regarding the measurement of primordial non-Gaussianity. We argue that when quantified in terms of the local (equilateral) template amplitude $f_{\rm NL}^{\rm loc}$ ($f_{\rm NL}^{\rm eq}$), natural target levels of sensitivity are $\Delta f_{\rm NL}^{\rm loc, eq.} \simeq 1$. We highlight that such levels are within reach of future surveys by measuring 2-, 3- and 4-point statistics of the galaxy spatial distribution. This paper summarizes a workshop held at CITA (University of Toronto) on October 23-24, 2014.

This paper provides a systematic study of renormalization in models of halo biasing. Building on work of McDonald, we show that Eulerian biasing is only consistent with renormalization if non-local terms and higher-derivative contributions are included in the biasing model. We explicitly determine the complete list of required bias parameters for Gaussian initial conditions, up to quartic order in the dark matter density contrast and at leading order in derivatives. At quadratic order, this means including the gravitational tidal tensor, while at cubic order the velocity potential appears as an independent degree of freedom. Our study naturally leads to an effective theory of biasing in which the halo density is written as a double expansion in fluctuations and spatial derivatives. We show that the bias expansion can be organized in terms of Galileon operators which aren't renormalized at leading order in derivatives. Finally, we discuss how the renormalized bias parameters impact the statistics of halos.

We show that the recent Planck limits on primordial non-Gaussianity impose strong constraints on light hidden sector fields coupled to the inflaton via operators suppressed by a high mass scale \Lambda. We study a simple effective field theory in which a hidden sector field is coupled to a shift-symmetric inflaton via arbitrary operators up to dimension five. Self-interactions in the hidden sector lead to non-Gaussianity in the curvature perturbations. To be consistent with the Planck limit on local non-Gaussianity, the coupling to any hidden sector with light fields and natural cubic couplings must be suppressed by a very high scale \Lambda > 10^5 H. Even if the hidden sector has Gaussian correlations, nonlinearities in the mixing with the inflaton still lead to non-Gaussian curvature perturbations. In this case, the non-Gaussianity is of the equilateral or orthogonal type, and the Planck data requires \Lambda > 10^2 H.

In this paper, we prove that the superhorizon conservation of the curvature perturbation zeta in single-field inflation holds as an operator statement. This implies that all zeta-correlators are time independent at all orders in the loop expansion. Our result follows directly from locality and diffeomorphism invariance of the underlying theory. We also explore the relationship between the conservation of zeta, the single-field consistency relation and the renormalization of composite operators.

Soft limits of inflationary correlation functions are both observationally relevant and theoretically robust. Various theorems can be proven about them that are insensitive to detailed model-building assumptions. In this paper, we re-derive several of these theorems in a universal way. Our method makes manifest why soft limits are such an interesting probe of the spectrum of additional light fields during inflation. We illustrate these abstract results with a detailed case study of the soft limits of quasi-single-field inflation.