The disagreement between direct late-time measurements of the Hubble constant from the SH0ES collaboration, and early-universe measurements based on the $\Lambda$CDM model from the Planck collaboration might, at least in principle, be explained by new physics in the early universe. Recently, the application of the Effective Field Theory of Large-Scale Structure to the full shape of the power spectrum of the SDSS/BOSS data has revealed a new, rather powerful, way to measure the Hubble constant and the other cosmological parameters from Large-Scale Structure surveys. In light of this, we analyze two models for early universe physics, Early Dark Energy and Rock 'n' Roll, that were designed to significantly ameliorate the Hubble tension. Upon including the information from the full shape to the Planck, BAO, and Supernovae measurements, we find that the degeneracies in the cosmological parameters that were introduced by these models are well broken by the data, so that these two models do not significantly ameliorate the tension.

We study the effect of the Einstein de-Sitter (EdS) approximation on the one-loop power spectrum of galaxies in redshift space in the Effective Field Theory of Large-Scale Structure. The dark matter density perturbations and velocity divergence are treated with exact time dependence. Splitting the density perturbation into its different temporal evolutions naturally gives rise to an irreducible basis of biases. While, as in the EdS approximation, at each time this basis spans a seven-dimensional space, this space is a slightly different one, and the difference is captured by a single calculable time- and $\vec k$-dependent function. We then compute the redshift-space galaxy one-loop power spectrum with the EdS approximation ($P^{\text{EdS-approx}}$) and without ($P^{\text{Exact}}$). For the monopole we find $P_{\text{0}}^{\text{Exact}}/P_{\text{0}}^{\text{EdS-approx}}\sim 1.003$ and for the quadrupole $P_{\text{2}}^{\text{Exact}}/P_{\text{2}}^{\text{EdS-approx}}\sim 1.007$ at $z=0.57$, and sharply increasing at lower redshifts. Finally, we show that a substantial fraction of the effect remains even after allowing the bias coefficients to shift within a physically allowed range. This suggests that the EdS approximation can only fit the data to a level of precision that is roughly comparable to the precision of the next generation of cosmological surveys. Furthermore, we find that implementing the exact time dependence formalism is not demanding and is easily applicable to data. Both of these points motivate a direct study of this effect on the cosmological parameters.

We study, using Mean Curvature Flow methods, 3+1 dimensional cosmologies with a positive cosmological constant, matter satisfying the dominant and the strong energy conditions, and with spatial slices that can be foliated by 2-dimensional surfaces that are the closed orbits of a symmetry group. If these surfaces have non-positive Euler characteristic (or in the case of 2-spheres, if the initial 2-spheres are large enough) and also if the initial spatial slice is expanding everywhere, then we prove that asymptotically the spacetime becomes physically indistinguishable from de Sitter space on arbitrarily large regions of spacetime. This holds true notwithstanding the presence of initial arbitrarily-large density fluctuations.

An accurate theoretical template for the galaxy power spectrum is a key for the success of ongoing and future spectroscopic surveys. We examine to what extent the Effective Field Theory of Large Scale Structure is able to provide such a template and correctly estimate cosmological parameters. To that end, we initiate a blinded challenge to infer cosmological parameters from the redshift-space power spectrum of high-resolution mock catalogs mimicking the BOSS galaxy sample but covering a hundred times larger cumulative volume. This gigantic simulation volume allows us to separate systematic bias due to theoretical modeling from the statistical error due to sample variance. The challenge task was to measure three unknown input parameters used in the simulation: the Hubble constant, the matter density fraction, and the clustering amplitude. We present analyses done by two independent teams, who have fitted the mock simulation data generated by yet another independent group. This allows us to avoid any confirmation bias by analyzers and pin down possible tuning of the specific EFT implementations. Both independent teams have recovered the true values of the input parameters within sub-percent statistical errors corresponding to the total simulation volume.

We apply the Effective Field Theory of Large-Scale Structure to analyze the $w$CDM cosmological model. By using the full shape of the power spectrum and the BAO post-reconstruction measurements from BOSS, the Supernovae from Pantheon, and a prior from BBN, we set the competitive CMB-independent limit $w=-1.046_{-0.052}^{+0.055}$ at $68\%$ C.L.. After adding the Planck CMB data, we find $w=-1.023_{-0.030}^{+0.033}$ at $68\%$ C.L.. Our results are obtained using PyBird, a new, fast Python-based code which we make publicly available.

We study, using Mean Curvature Flow methods, 2+1 dimensional cosmologies with a positive cosmological constant and matter satisfying the dominant and the strong energy conditions. If the spatial slices are compact with non-positive Euler characteristic and are initially expanding everywhere, then we prove that the spatial slices reach infinite volume, asymptotically converge on average to de Sitter and they become, almost everywhere, physically indistinguishable from de Sitter. This holds true notwithstanding the presence of initial arbitrarily-large density fluctuations and the formation of black holes.

The precision of the cosmological data allows us to accurately approximate the predictions for cosmological observables by Taylor expanding up to a low order the dependence on the cosmological parameters around a reference cosmology. By applying this observation to the redshift-space one-loop galaxy power spectrum of the Effective Field Theory of Large-Scale Structure, we analyze the BOSS DR12 data by scanning over all the parameters of $\Lambda$CDM cosmology with massive neutrinos. We impose several sets of priors, the widest of which is just a Big Bang Nucleosynthesis prior on the current fractional energy density of baryons, $\Omega_b h^2$, and a bound on the sum of neutrino masses to be less than 0.9 eV. In this case we measure the primordial amplitude of the power spectrum, $A_s$, the abundance of matter, $\Omega_m$, the Hubble parameter, $H_0$, and the tilt of the primordial power spectrum, $n_s$, to about $19\%$, $5.7\%$, $2.2\%$ and $7.3\%$ respectively, obtaining $\ln ( 10^{10} A_s) =2.91\pm 0.19$, $\Omega_m=0.314\pm 0.018$, $H_0=68.7\pm 1.5$ km/(s Mpc) and $n_s=0.979\pm 0.071$ at $68\%$ confidence level. A public code is released with this preprint.

The Effective Field Theory of Large-Scale Structure (EFTofLSS) is a formalism that allows us to predict the clustering of Cosmological Large-Scale Structure in the mildly non-linear regime in an accurate and reliable way. After validating our technique against several sets of numerical simulations, we perform the analysis for the cosmological parameters of the DR12 BOSS data. We assume $\Lambda$CDM, a fixed value of the baryon/dark-matter ratio, $\Omega_b/\Omega_c$, and of the tilt of the primordial power spectrum, $n_s$, and no significant input from numerical simulations. By using the one-loop power spectrum multipoles, we measure the primordial amplitude of the power spectrum, $A_s$, the abundance of matter, $\Omega_m$, and the Hubble parameter, $H_0$, to about $13\%$, $3.2\%$ and $3.2\%$ respectively, obtaining $\ln(10^{10}As)=2.72\pm 0.13$, $\Omega_m=0.309\pm 0.010$, $H_0=68.5\pm 2.2$ km/(s Mpc) at 68\% confidence level. If we then add a CMB prior on the sound horizon, the error bar on $H_0$ is reduced to $1.6\%$. These results are a substantial qualitative and quantitative improvement with respect to former analyses, and suggest that the EFTofLSS is a powerful instrument to extract cosmological information from Large-Scale Structure.

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space.

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale.

The hot dense environment of the early universe is known to have produced large numbers of baryons, photons, and neutrinos. These extreme conditions may have also produced other long-lived species, including new light particles (such as axions or sterile neutrinos) or gravitational waves. The gravitational effects of any such light relics can be observed through their unique imprint in the cosmic microwave background (CMB), the large-scale structure, and the primordial light element abundances, and are important in determining the initial conditions of the universe. We argue that future cosmological observations, in particular improved maps of the CMB on small angular scales, can be orders of magnitude more sensitive for probing the thermal history of the early universe than current experiments. These observations offer a unique and broad discovery space for new physics in the dark sector and beyond, even when its effects would not be visible in terrestrial experiments or in astrophysical environments. A detection of an excess light relic abundance would be a clear indication of new physics and would provide the first direct information about the universe between the times of reheating and neutrino decoupling one second later.

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

We compute the tree-level bispectrum of dark matter in the presence of massive neutrinos in the mildly non-linear regime in the context of the effective field theory of large-scale structure (EFTofLSS). For neutrinos, whose typical free streaming wavenumber ($k_{\rm fs}$) is longer than the non-linear scale ($k_{\mathrm{NL}}$), we solve a Boltzmann equation coupled to the effective fluid equation for dark matter. We solve perturbatively the coupled system by expanding in powers of the neutrino density fraction ($f_{\nu}$) and the ratio of the wavenumber of interest over the non-linear scale ($k/k_{\mathrm{NL}}$) and add suitable counterterms to remove the dependence from short distance physics. For equilateral configurations, we find that the total-matter tree-level bispectrum is approximately $16f_{\nu}$ times the dark matter one on short scales ($k > k_{\rm fs}$). The largest contribution stems from the back-reaction of massive neutrinos on the dark matter growth factor. On large scales ($k < k_{\rm fs}$) the contribution of neutrinos to the bispectrum is smaller by up to two orders of magnitude.

Effective field theory methods suggest that some rather-general extensions of General Relativity include, or are mimicked by, certain higher-order curvature corrections, with coupling constants expected to be small but otherwise arbitrary. Thus, the tantalizing prospect to test the fundamental nature of gravity with gravitational-wave observations, in a systematic way, emerges naturally. Here, we build black hole solutions in such a framework and study their main properties. Once rotation is included, we find the first purely gravitational example of geometries without $\mathbb{Z}_2$-symmetry. Despite the higher-order operators of the theory, we show that linearized fluctuations of such geometries obey second-order differential equations. We find nonzero tidal Love numbers. We study and compute the quasinormal modes of such geometries. These results are of interest to gravitational-wave science but also potentially relevant for electromagnetic observations of the galactic center or $X$-ray binaries.

We develop an analytic method for implementing the IR-resummation of arXiv:1404.5954, which allows one to correctly and consistently describe the imprint of baryon acoustic oscillations (BAO) on statistical observables in large-scale structure. We show that the final IR-resummed correlation function can be computed analytically without relying on numerical integration, thus allowing for an efficient and accurate use of these predictions on real data in cosmological parameter fitting. In this work we focus on the one-loop correlation function and the BAO peak. We show that, compared with the standard numerical integration method of IR-resummation, the new method is accurate to better than 0.2 %, and is quite easily improvable. We also give an approximate resummation scheme which is based on using the linear displacements of a fixed fiducial cosmology, which when combined with the method described above, is about six times faster than the standard numerical integration. Finally, we show that this analytic method is generalizable to higher loop computations.

The existence of light particles with spin during inflation is prohibited by the Higuchi bound. This conclusion can be evaded if one considers states with a sizeable coupling with the inflaton foliation, since this breaks the de Sitter isometries. The action for these states can be constructed within the Effective Field Theory of Inflation, or using a CCWZ procedure. Light particles with spin have prescribed couplings with soft inflaton perturbations, which are encoded in consistency relations. We study the phenomenology of light states with spin 2. These mix with the graviton changing the tensor power spectrum and can lead to sizeable tensor non-Gaussianities. They also give rise to a scalar bispectrum and trispectrum with a characteristic angle-dependent non-Gaussianity.

We further develop the description of redshift space distortions within the Effective Field Theory of Large Scale Structures. First, we generalize the counterterms to include the effect of baryonic physics and primordial non-Gaussianity. Second, we evaluate the IR-resummation of the dark matter power spectrum in redshift space. This requires us to identify a controlled approximation that makes the numerical evaluation straightforward and efficient. Third, we compare the predictions of the theory at one loop with the power spectrum from numerical simulations up to $\ell=6$. We find that the IR-resummation allows us to correctly reproduce the BAO peak. The $k$-reach, or equivalently the precision for a given $k$, depends on additional counterterms that need to be matched to simulations. Since the non-linear scale for the velocity is expected to be longer than the one for the overdensity, we consider a minimal and a non-minimal set of counterterms. The quality of our numerical data makes it hard to firmly establish the performance of the theory at high wavenumbers. Within this limitation, we find that the theory at redshift $z=0.56$ and up to $\ell=2$ matches the data to percent level approximately up to $k \sim 0.13 \, h { \rm Mpc^{-1}}$ or $k \sim 0.18 \, h { \rm Mpc^{-1}}$, depending on the number of counterterms used, with potentially large improvement over former analytical techniques.

The Effective Field Theory of Large-Scale Structure (EFTofLSS) provides a consistent perturbative framework for describing the statistical distribution of cosmological large-scale structure. In a previous EFTofLSS calculation that involved the one-loop power spectra and tree-level bispectra, it was shown that the $k$-reach of the prediction for biased tracers is comparable for all investigated masses if suitable higher-derivative biases, which are less suppressed for more massive tracers, are added. However, it is possible that the non-linear biases grow faster with tracer mass than the linear bias, implying that loop contributions could be the leading correction to the bispectra. To check this, we include the one-loop contributions in a fit to numerical data in the limit of strongly enhanced higher-order biases. We show that the resulting one-loop power spectra and higher-derivative plus leading one-loop bispectra fit the two- and three-point functions respectively up to $k\simeq 0.19\ h\ \rm{Mpc}^{-1}$ and $k\simeq 0.14\ h\ \rm{Mpc}^{-1}$ at the percent level. We find that the higher-order bias coefficients are not strongly enhanced, and we argue that the gain in perturbative reach due to the leading one-loop contributions to the bispectra is relatively small. Thus, we conclude that higher-derivative biases provide the leading correction to the bispectra for tracers of a very wide range of masses.

We propose a simplification for the IR-resummation scheme of Senatore and Zaldarriaga, and also include its next-to-leading order corrections coming from the tree-level three-point function of the long displacement field. First we show that the new simplified formula shares the same properties of the resummation of Senatore and Zaldarriaga. In Fourier space, the IR-resummed power spectrum has no residual wiggles and the two-loop calculation matches the non-linear power spectrum of the Dark Sky simulation at $z=0$ up to $k\simeq0.34\,h\,\text{Mpc}^{-1}$ within cosmic variance. Then, we find that the additional subleading terms (although parametrically infrared-enhanced) modify the leading-order IR-resummed correlation function only in a marginal way, implying that the IR-resummation scheme can robustly predict the shape of the BAO peak.

The recent direct observation of gravitational waves (GW) from merging black holes opens up the possibility of exploring the theory of gravity in the strong regime at an unprecedented level. It is therefore interesting to explore which extensions to General Relativity (GR) could be detected. We construct an Effective Field Theory (EFT) satisfying the following requirements. It is testable with GW observations; it is consistent with other experiments, including short distance tests of GR; it agrees with widely accepted principles of physics, such as locality, causality and unitarity; and it does not involve new light degrees of freedom. The most general theory satisfying these requirements corresponds to adding to the GR Lagrangian operators constructed out of powers of the Riemann tensor, suppressed by a scale comparable to the curvature of the observed merging binaries. The presence of these operators modifies the gravitational potential between the compact objects, as well as their effective mass and current quadrupoles, ultimately correcting the waveform of the emitted GW.

In the next few years, we are going to probe the low-redshift universe with unprecedented accuracy. Among the various fruits that this will bear, it will greatly improve our knowledge of the dynamics of dark energy, though for this there is a strong theoretical preference for a cosmological constant. We assume that dark energy is described by the so-called Effective Field Theory of Dark Energy, which assumes that dark energy is the Goldstone boson of time translations. Such a formalism makes it easy to ensure that our signatures are consistent with well-established principles of physics. Since most of the information resides at high wavenumbers, it is important to be able to make predictions at the highest wavenumber that is possible. The Effective Field Theory of Large-Scale Structure (EFTofLSS) is a theoretical framework that has allowed us to make accurate predictions in the mildly non-linear regime. In this paper, we derive the non-linear equations that extend the EFTofLSS to include the effect of dark energy both on the matter fields and on the biased tracers. For the specific case of clustering quintessence, we then perturbatively solve to cubic order the resulting non-linear equations and construct the one-loop power spectrum of the total density contrast.

Because large-scale structure surveys may very well be the next leading sources of cosmological information, it is important to have a precise understanding of the cosmological observables; for this reason, the Effective Field Theory of Large-Scale Structure (EFTofLSS) was developed. So far, most results in the EFTofLSS have used the so-called Einstein-de Sitter approximation, an approximation of the time dependence which is known to be accurate to better than one percent. However, in order to reach even higher accuracy, the full time dependence must be used. The computation with exact time dependence is sensitive to both infrared (IR) and ultraviolet (UV) effects in the loop integrands, and while these effects must cancel because of diffeomorphism invariance, they make numerical computation much less efficient. We provide a formulation of the one-loop, equal-time exact-time-dependence power spectrum of density perturbations which is manifestly free of these spurious IR and UV divergences at the level of the integrand. We extend our results to the total matter mode with clustering quintessence, show that IR and UV divergences cancel, and provide the associated IR- and UV-safe integrand. This also establishes that the consistency conditions are satisfied in this system. We then use our one-loop result to do an improved precision comparison of the two-loop dark-matter power spectrum with the Dark Sky N-body simulation.

We develop a formalism to analytically describe the clustering of matter in the mildly non-linear regime in the presence of massive neutrinos. Neutrinos, whose free streaming wavenumber ($k_{\rm fs}$) is typically longer than the non-linear scale ($k_{\rm NL}$) are described by a Boltzmann equation coupled to the effective fluid-like equations that describe dark matter. We solve the equations expanding in the neutrino density fraction $(f_\nu)$ and in $k/ k_{\rm NL}$, and add suitable counterterms to renormalize the theory. This allows us to describe the contribution of short distances to long-distance observables. Equivalently, we construct an effective Boltzmann equation where we add additional terms whose coefficients renormalize the contribution from short-distance physics. We argue that neutrinos with $k_{\rm fs}\gtrsim k_{\rm NL}$ require an additional counterterm similar to the speed of sound ($c_s$) for dark matter. We compute the one-loop total-matter power spectrum, and find that it is roughly equal to $16f_\nu$ times the dark matter one for $k$'s larger that the typical $k_{\rm fs}$. It is about half of that for smaller $k$'s. The leading contribution results from the back-reaction of the neutrinos on the dynamics of the dark matter. The counterterms contribute in a hierarchical way: the leading ones can either be computed in terms of $c_s$, or can be accounted for by shifting $c_s$ by an amount proportional to $f_\nu$.

The most effective use of data from current and upcoming large scale structure~(LSS) and CMB observations requires the ability to predict the clustering of LSS with very high precision. The Effective Field Theory of Large Scale Structure (EFTofLSS) provides an instrument for performing analytical computations of LSS observables with the required precision in the mildly nonlinear regime. In this paper, we develop efficient implementations of these computations that allow for an exploration of their dependence on cosmological parameters. They are based on two ideas. First, once an observable has been computed with high precision for a reference cosmology, for a new cosmology the same can be easily obtained with comparable precision just by adding the difference in that observable, evaluated with much less precision. Second, most cosmologies of interest are sufficiently close to the Planck best-fit cosmology that observables can be obtained from a Taylor expansion around the reference cosmology. These ideas are implemented for the matter power spectrum at two loops and are released as public codes. When applied to cosmologies that are within 3$\sigma$ of the Planck best-fit model, the first method evaluates the power spectrum in a few minutes on a laptop, with results that have 1\% or better precision, while with the Taylor expansion the same quantity is instantly generated with similar precision. The ideas and codes we present may easily be extended for other applications or higher-precision results.

If time-translations are spontaneously broken, so are boosts. This symmetry breaking pattern can be non-linearly realized by either just the Goldstone boson of time translations, or by four Goldstone bosons associated with time translations and boosts. In this paper we extend the Effective Field Theory of Multifield Inflation to consider the case in which the additional Goldstone bosons associated with boosts are light and coupled to the Goldstone boson of time translations. The symmetry breaking pattern forces a coupling to curvature so that the mass of the additional Goldstone bosons is predicted to be equal to $\sqrt{2}H$ in the vast majority of the parameter space where they are light. This pattern therefore offers a natural way of generating self-interacting particles with Hubble mass during inflation. After constructing the general effective Lagrangian, we study how these particles mix and interact with the curvature fluctuations, generating potentially detectable non-Gaussian signals.

In homogeneous and isotropic Friedmann-Robertson-Walker cosmology, the topology of the universe determines its ultimate fate. If the Weak Energy Condition is satisfied, open and flat universes must expand forever, while closed cosmologies can recollapse to a Big Crunch. A similar statement holds for homogeneous but anisotropic (Bianchi) universes. Here, we prove that $arbitrarily$ inhomogeneous and anisotropic cosmologies with "flat" (including toroidal) and "open" (including compact hyperbolic) spatial topology that are initially expanding must continue to expand forever at least in some region at a rate bounded from below by a positive number, despite the presence of arbitrarily large density fluctuations and/or the formation of black holes. Because the set of 3-manifold topologies is countable, a single integer determines the ultimate fate of the universe, and, in a specific sense, most 3-manifolds are "flat" or "open". Our result has important implications for inflation: if there is a positive cosmological constant (or suitable inflationary potential) and initial conditions for the inflaton, cosmologies with "flat" or "open" topology must expand forever in some region at least as fast as de Sitter space, and are therefore very likely to begin inflationary expansion eventually, regardless of the scale of the inflationary energy or the spectrum and amplitude of initial inhomogeneities and gravitational waves. Our result is also significant for numerical general relativity, which often makes use of periodic (toroidal) boundary conditions.

The Effective Field Theory of Large-Scale Structure (EFTofLSS) provides a novel formalism that is able to accurately predict the clustering of large-scale structure (LSS) in the mildly non-linear regime. Here we provide the first computation of the power spectrum of biased tracers in redshift space at one loop order, and we make the associated code publicly available. We compare the multipoles $\ell=0,2$ of the redshift-space halo power spectrum, together with the real-space matter and halo power spectra, with data from numerical simulations at $z=0.67$. For the samples we compare to, which have a number density of $\bar n=3.8 \cdot 10^{-2}(h \ {\rm Mpc}^{-1})^3$ and $\bar n=3.9 \cdot 10^{-4}(h \ {\rm Mpc}^{-1})^3$, we find that the calculation at one-loop order matches numerical measurements to within a few percent up to $k\simeq 0.43 \ h \ {\rm Mpc}^{-1}$, a significant improvement with respect to former techniques. By performing the so-called IR-resummation, we find that the Baryon Acoustic Oscillation peak is accurately reproduced. Based on the results presented here, long-wavelength statistics that are routinely observed in LSS surveys can be finally computed in the EFTofLSS. This formalism thus is ready to start to be compared directly to observational data.

We construct the Supersymmetric Effective Field Theory of Inflation, that is the most general theory of inflationary fluctuations when time-translations and supersymmetry are spontaneously broken. The non-linear realization of these invariances allows us to define a complete SUGRA multiplet containing the graviton, the gravitino, the Goldstone of time translations and the Goldstino, with no auxiliary fields. Going to a unitary gauge where only the graviton and the gravitino are present, we write the most general Lagrangian built out of the fluctuations of these fields, invariant under time-dependent spatial diffeomorphisms, but softly-breaking time diffeomorphisms and gauged SUSY. With a suitable St\"uckelberg transformation, we introduce the Goldstone boson of time translation and the Goldstino of SUSY. No additional dynamical light field is needed. In the high energy limit, larger than the inflationary Hubble scale for the Goldstino, these fields decouple from the graviton and the gravitino, greatly simplifying the analysis in this regime. We study the phenomenology of this Lagrangian. The Goldstino can have a non-relativistic dispersion relation. Gravitino and Goldstino affect the primordial curvature perturbations at loop level. The UV modes running in the loops generate three-point functions which are degenerate with the ones coming from operators already present in the absence of supersymmetry. Their size is potentially as large as corresponding to $f_{\rm NL}^{\rm equil.,\,orthog.}\sim 1$ or, for particular operators, even $\gg 1$. The non-degenerate contribution from modes of order $H$ is estimated to be very small.

Planning to explore the beginning of the Universe? A lightweight introductory guide to the theory of Inflation.

Most of the upcoming cosmological information will come from analyzing the clustering of the Large Scale Structures (LSS) of the universe through LSS or CMB observations. It is therefore essential to be able to understand their behavior with exquisite precision. The Effective Field Theory of Large Scale Structures (EFTofLSS) provides a consistent framework to make predictions for LSS observables in the mildly non-linear regime. In this paper we focus on biased tracers. We argue that in calculations at a given order in the dark matter perturbations, highly biased tracers will underperform because of their larger higher derivative biases. A natural prediction of the EFTofLSS is therefore that by simply adding higher derivative biases, all tracers should perform comparably well. We implement this prediction for the halo-halo and the halo-matter power spectra at one loop, and the halo-halo-halo, halo-halo-matter, and halo-matter-matter bispectra at tree-level, and compare with simulations. We find good agreement with the prediction: for all tracers, we are able to match the bispectra up to $k\simeq0.17\,h/$Mpc at $z=0$ and the power spectra to a higher wavenumber.

Using numerical solutions of the full Einstein field equations coupled to a scalar inflaton field in 3+1 dimensions, we study the conditions under which a universe that is initially expanding, highly inhomogeneous and dominated by gradient energy can transition to an inflationary period. If the initial scalar field variations are contained within a sufficiently flat region of the inflaton potential, and the universe is spatially flat or open on average, inflation will occur following the dilution of the gradient and kinetic energy due to expansion. This is the case even when the scale of the inhomogeneities is comparable to the initial Hubble length, and overdense regions collapse and form black holes, because underdense regions continue expanding, allowing inflation to eventually begin. This establishes that inflation can arise from highly inhomogeneous initial conditions and solve the horizon and flatness problems, at least as long as the variations in the scalar field do not include values that exceed the inflationary plateau.

We study the prediction of the dark matter power spectrum at two-loop order in the Effective Field Theory of Large Scale Structures (EFTofLSS) using high precision numerical simulations. In our universe, short distance non-linear perturbations, not under perturbative control, affect long distance fluctuations through an effective stress tensor that needs to be parametrized in terms of counterterms that are functions of the long distance fluctuating fields. We find that at two-loop order it is necessary to include three counterterms: a linear term in the over density, $\delta$, a quadratic term, $\delta^2$, and a higher derivative term, $\partial^2\delta$. After the inclusion of these three terms, the EFTofLSS at two-loop order matches simulation data up to $k\simeq 0.34 \,h\, {\rm Mpc}^{-1}$ at redshift $z=0$, up to $k\simeq 0.55\,h\, {\rm Mpc}^{-1}$ at $z=1$, and up to $k\simeq 1.1\,h\, {\rm Mpc}^{-1}$ at $z=2$. At these wavenumbers, the cosmic variance of the simulation is at least as small as $10^{-3}$, providing a high precision comparison between theory and data. The actual reach of the theory is affected by theoretical uncertainties associated to not having included higher order terms in perturbation theory, for which we provide an estimate, and by potentially overfitting the data, which we also try to address. Since in the EFTofLSS the coupling constants associated with the counterterms are unknown functions of time, we show how a simple parametrization gives a sensible description of their time-dependence. Overall, the $k$-reach of the EFTofLSS is much larger than previous analytical techniques, showing that the amount of cosmological information amenable to high-precision analytical control might be much larger than previously believed.

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_{\rm NL}$, where $k$ is the wavenumber of interest and $k_{\rm 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\simeq 0.3 h\,$Mpc$^{-1}$ and $k\simeq 0.6 h\,$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_{\rm NL}$ is smaller than in the previously considered cases ($\sigma_8=0.9$), we find that the prediction from the EFTofLSS agrees very well with $N$-body simulations up to $k\simeq 0.25 h\,$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\simeq 0.25 h\,$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.

We analyze the shape and amplitude of oscillatory features in the primordial power spectrum and non-Gaussianity induced by periodic production of heavy degrees of freedom coupled to the inflaton $\phi$. We find that non-adiabatic production of particles can contribute effects which are detectable or constrainable using cosmological data even if their time-dependent masses are always heavier than the scale $\dot \phi^{1/2}$, much larger than the Hubble scale. This provides a new role for UV completion, consistent with the criteria from effective field theory for when heavy fields cannot be integrated out. This analysis is motivated in part by the structure of axion monodromy, and leads to an additional oscillatory signature in a subset of its parameter space. At the level of a quantum field theory model that we analyze in detail, the effect arises consistently with radiative stability for an interesting window of couplings up to of order $\lesssim 1$. The amplitude of the bispectrum and higher-point functions can be larger than that for Resonant Non-Gaussianity, and its signal/noise may be comparable to that of the corresponding oscillations in the power spectrum (and even somewhat larger within a controlled regime of parameters). Its shape is distinct from previously analyzed templates, but was partly motivated by the oscillatory equilateral searches performed recently by the {\it Planck} collaboration. We also make some general comments about the challenges involved in making a systematic study of primordial non-Gaussianity.

We study the prediction of the Effective Field Theory of Large Scale Structures (EFTofLSS) for the matter power spectrum at different redshifts. In previous work, we found that the two-loop prediction can match the nonlinear power spectrum measured from $N$-body simulations at redshift zero within approximately 2% up to $k\sim 0.6\,h\, {\rm Mpc}^{-1}$ after fixing a single free parameter, the so-called "speed of sound". We determine the time evolution of this parameter by matching the EFTofLSS prediction to simulation output at different redshifts, and find that it is well-described by a fitting function that only includes one additional parameter. After the two free parameters are fixed, the prediction agrees with nonlinear data within approximately 2% up to at least $k\sim 1\,h\, {\rm Mpc}^{-1}$ at $z\geq 1$, and also within approximately 5% up to $k\sim 1.2\,h\, {\rm Mpc}^{-1}$ at $z=1$ and $k\sim 2.3\,h\, {\rm Mpc}^{-1}$ at $z=3$, a major improvement with respect to other perturbative techniques. We also develop an accurate way to estimate where the EFTofLSS predictions at different loop orders should fail, based on the sizes of the next-order terms that are neglected, and find agreement with the actual comparisons to data. Finally, we use our matter power spectrum results to perform analytical calculations of lensing potential power spectra corresponding to both CMB and galaxy lensing. This opens the door to future direct applications of the EFTofLSS to observations of gravitational clustering on cosmic scales.

With the completion of the Planck mission, in order to continue to gather cosmological information it has become crucial to understand the Large Scale Structures (LSS) of the universe to percent accuracy. The Effective Field Theory of LSS (EFTofLSS) is a novel theoretical framework that aims to develop an analytic understanding of LSS at long distances, where inhomogeneities are small. We further develop the description of biased tracers in the EFTofLSS to account for the effect of baryonic physics and primordial non-Gaussianities, finding that new bias coefficients are required. Then, restricting to dark matter with Gaussian initial conditions, we describe the prediction of the EFTofLSS for the one-loop halo-halo and halo-matter two-point functions, and for the tree-level halo-halo-halo, matter-halo-halo and matter-matter-halo three-point functions. Several new bias coefficients are needed in the EFTofLSS, even though their contribution at a given order can be degenerate and the same parameters contribute to multiple observables. We develop a method to reduce the number of biases to an irreducible basis, and find that, at the order at which we work, seven bias parameters are enough to describe this extremely rich set of statistics. We then compare with the output of $N$-body simulations. For the lowest mass bin, we find percent level agreement up to $k\simeq 0.3\,h\,{\rm Mpc}^{-1}$ for the one-loop two-point functions, and up to $k\simeq 0.15\,h\,{\rm Mpc}^{-1}$ for the tree-level three-point functions, with the $k$-reach decreasing with higher mass bins. This is consistent with the theoretical estimates, and suggests that the cosmological information in LSS amenable to analytical control is much more than previously believed.

We study the local response to long wavelength fluctuations in cosmological $N$-body simulations, focusing on the matter and halo power spectra, halo abundance and non-linear transformations of the density field. The long wavelength mode is implemented using an effective curved cosmology and a mapping of time and distances. The method provides an alternative, most probably more precise, way to measure the isotropic halo biases. Limiting ourselves to the linear case, we find generally good agreement between the biases obtained from the curvature method and the traditional power spectrum method at the level of a few percent. We also study the response of halo counts to changes in the variance of the field and find that the slope of the relation between the responses to density and variance differs from the naive derivation assuming a universal mass function by 18%. This has implications for measurements of the amplitude of local non-Gaussianity using scale dependent bias. We also analyze the halo power spectrum and halo-dark matter cross-spectrum response to long wavelength fluctuations and derive second order halo bias from it, as well as the super-sample variance contribution to the galaxy power spectrum covariance matrix.

The Standard Model Higgs potential becomes unstable at large field values. After clarifying the issue of gauge dependence of the effective potential, we study the cosmological evolution of the Higgs field in presence of this instability throughout inflation, reheating and the present epoch. We conclude that anti-de Sitter patches in which the Higgs field lies at its true vacuum are lethal for our universe. From this result, we derive upper bounds on the Hubble constant during inflation, which depend on the reheating temperature and on the Higgs coupling to the scalar curvature or to the inflaton. Finally we study how a speculative link between Higgs meta-stability and consistence of quantum gravity leads to a sharp prediction for the Higgs and top masses, which is consistent with measured values.

We revisit alternative mechanisms of gravitational wave production during inflation and argue that they generically emit a non-negligible amount of scalar fluctuations. We find the scalar power is larger than the tensor power by a factor of order $1/\epsilon^2$. For an appreciable tensor contribution the associated scalar emission completely dominates the zero-point fluctuations of inflaton, resulting in a tensor-to-scalar ratio $r\sim \epsilon^2$. A more quantitative result can be obtained if one further assumes that gravitational waves are emitted by localized sub-horizon processes, giving $r_{\rm max} \simeq 0.3 \epsilon^2$. However, $\epsilon$ is generally time dependent, and this result for $r$ depends on its instantaneous value during the production of the sources, rather than just its average value, somewhat relaxing constraints from the tilt $n_s$. We calculate the scalar 3-point correlation function in the same class of models and show that non-Gaussianity cannot be made arbitrarily small, i.e. $f_{NL} \geq 1$, independently of the value of $r$. Possible exceptions in multi-field scenarios are discussed.

We present a new method to resum the effect of large scale motions in the Effective Field Theory of Large Scale Structures. Because the linear power spectrum in $\Lambda$CDM is not scale free the effects of the large scale flows are enhanced. Although previous EFT calculations of the equal-time density power spectrum at one and two loops showed a remarkable agreement with numerical results, they also showed a 2% residual which appeared related to the BAO oscillations. We show that this was indeed the case, explain the physical origin and show how a Lagrangian based calculation removes this differences. We propose a simple method to upgrade existing Eulerian calculations to effectively make them Lagrangian and compare the new results with existing fits to numerical simulations. Our new two-loop results agrees with numerical results up to $k\sim 0.6 h/$Mpc to within 1% with no oscillatory residuals. We also compute power spectra involving momentum which is significantly more affected by the large scale flows. We show how keeping track of these velocities significantly enhances the UV reach of the momentum power spectrum in addition to removing the BAO related residuals. We compute predictions for the real space correlation function around the BAO scale and investigate its sensitivity to the EFT parameters and the details of the resummation technique.

Last year we argued that if slow-roll inflation followed the decay of a false vacuum in a large landscape, the steepening of the scalar potential between the inflationary plateau and the barrier generically leads to a potentially observable suppression of the scalar power spectrum at large distances. Here we revisit this analysis in light of the recent BICEP2 results. Assuming that both the BICEP2 B-mode signal and the Planck analysis of temperature fluctuations hold up, we find that the data now discriminate more sharply between our scenario and $\Lambda$CDM. Nonzero tensor modes exclude standard $\Lambda$CDM with notable but not yet conclusive confidence: at $\sim 3.8\,\sigma$ if $r\approx0.2$, or at $\sim 3.5\,\sigma$ if $r=0.15$. Of the two steepening models of our previous work, one is now ruled out by existing bounds on spatial curvature. The other entirely reconciles the tension between BICEP2 and Planck. Upcoming $EE$ polarization measurements have the potential to rule out unmodified $\Lambda$CDM decisively. Next generation Large Scale Structure surveys can further increase the significance. More precise measurements of $BB$ at low $\ell$ will help distinguish our scenario from other explanations. If steepening is confirmed, the prospects for detecting open curvature increase but need not be large.

We develop a general framework for data analysis and phenomenology of the CMB four-point function or trispectrum. To lowest order in the derivative expansion, the inflationary action admits three quartic operators consistent with symmetry: $\dot\sigma^4$, $\dot\sigma^2 (\partial\sigma^2)$, and $(\partial\sigma)^4$. In single field inflation, only the first of these operators can be the leading non-Gaussian signal. A Fisher matrix analysis shows that there is one near-degeneracy among the three CMB trispectra, so we parameterize the trispectrum with two coefficients $g_{NL}^{\dot\sigma^4}$ and $g_{NL}^{(\partial\sigma)^4}$, in addition to the coefficient $g_{NL}^{\rm loc}$ of $\zeta^3$-type local non-Gaussianity. This three-parameter space is analogous to the parameter space $(f_{NL}^{\rm loc}, f_{NL}^{\rm equil}, f_{NL}^{\rm orth})$ commonly used to parameterize the CMB three-point function. We next turn to data analysis and show how to represent these trispectra in a factorizable form which leads to computationally fast operations such as evaluating a CMB estimator or simulating a non-Gaussian CMB. We discuss practical issues in CMB analysis pipelines, and perform an optimal analysis of WMAP data. Our minimum-variance estimates are $g_{NL}^{\rm loc} = (-3.80 \pm 2.19) \times 10^5$, $g_{NL}^{\dot\sigma^4} = (-3.20 \pm 3.09) \times 10^6$, and $g_{NL}^{(\partial\sigma)^4} = (-10.8 \pm 6.33) \times 10^5$ after correcting for the effects of CMB lensing. No evidence of a nonzero inflationary four-point function is seen.

The large scale structures of the universe will likely be the next leading source of cosmological information. It is therefore crucial to understand their behavior. The Effective Field Theory of Large Scale Structures provides a consistent way to perturbatively predict the clustering of dark matter at large distances. The fact that baryons move distances comparable to dark matter allows us to infer that baryons at large distances can be described in a similar formalism: the backreaction of short-distance non-linearities and of star-formation physics at long distances can be encapsulated in an effective stress tensor, characterized by a few parameters. The functional form of baryonic effects can therefore be predicted. In the power spectrum the leading contribution goes as $\propto k^2 P(k)$, with $P(k)$ being the linear power spectrum and with the numerical prefactor depending on the details of the star-formation physics. We also perform the resummation of the contribution of the long-wavelength displacements, allowing us to consistently predict the effect of the relative motion of baryons and dark matter. We compare our predictions with simulations that contain several implementations of baryonic physics, finding percent agreement up to relatively high wavenumbers such as $k\simeq 0.3\,h\, Mpc^{-1}$ or $k\simeq 0.6\, h\, Mpc^{-1}$, depending on the order of the calculation. Our results open a novel way to understand baryonic effects analytically, as well as to interface with simulations.

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.

We study how to describe collapsed objects, such as galaxies, in the context of the Effective Field Theory of Large Scale Structures. The overdensity of galaxies at a given location and time is determined by the initial tidal tensor, velocity gradients and spatial derivatives of the regions of dark matter that, during the evolution of the universe, ended up at that given location. Similarly to what recently done for dark matter, we show how this Lagrangian space description can be recovered by upgrading simpler Eulerian calculations. We describe the Eulerian theory. We show that it is perturbatively local in space, but non-local in time, and we explain the observational consequences of this fact. We give an argument for why to a certain degree of accuracy the theory can be considered as quasi time-local and explain what the operator structure is in this case. We describe renormalization of the bias coefficients so that, after this and after upgrading the Eulerian calculation to a Lagrangian one, the perturbative series for galaxies correlation functions results in a manifestly convergent expansion in powers of $k/k_{\rm NL}$ and $k/k_{\rm M}$, where $k$ is the wavenumber of interest, $k_{\rm NL}$ is the wavenumber associated to the non-linear scale, and $k_{\rm M}$ is the comoving wavenumber enclosing the mass of a galaxy.

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, 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, 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.

We point out that detectable inflationary tensor modes can be generated by particle or string sources produced during inflation, consistently with the requirements for inflation and constraints from scalar fluctuations. We show via examples that this effect can dominate over the contribution from quantum fluctuations of the metric, occurring even when the inflationary potential energy is too low to produce a comparable signal. Thus a detection of tensor modes from inflation does not automatically constitute a determination of the inflationary Hubble scale.

We introduce a formalism, valid both for dark matter and collapsed objects, that allows us to describe redshift space distortions in the context of the Effective Field Theory of Large Scale Structures (EFTofLSS). Expressing density perturbations in redshift space corresponds to performing a change of coordinates and the resulting expressions contain products of density perturbations and velocity fields evaluated at the same location. These terms are sensitive to non-perturbative short-distance physics and in order to correctly treat them they need to be renormalized by adding suitable counterterms. Therefore more counterterms are required in redshift space expressions compared to their real space analogs. In particular in the expression for the one-loop matter power spectrum there are two new counterterms. Just as in real space, long wavelength displacements affect correlation functions in redshift space and need to be resummed. We generalize the real space formulas for IR resummation to this case: the final expressions are conceptually similar but are more challenging to compute numerically due to their reduced symmetry.

Given the fantastic experimental effort, it is important to thoroughly explore the signature space of inflationary models. The fact that higher derivative operators do not renormalize lower derivative ones allows us to find a large class of technically natural single-clock inflationary models where, in the context of the Effective Field Theory of Inflation, the leading interactions have many derivatives. We systematically explore the 3-point function induced by these models and their overlap with the standard equilateral and orthogonal templates. We find that in order to satisfactorily cover the signature space of these models, two new additional templates need to be included. We then perform the optimal analysis of the WMAP 9-year data for the resulting four templates, finding that the overall significance of a non-zero signal is between 2--2.5$\sigma$, depending on the choice of parameter space, partially driven by the preference for nonzero $f_{NL}^{\rm orth}$ in WMAP9.

Fluctuations in the intensity and polarization of the cosmic microwave background (CMB) and the large-scale distribution of matter in the universe each contain clues about the nature of the earliest moments of time. The next generation of CMB and large-scale structure (LSS) experiments are poised to test the leading paradigm for these earliest moments---the theory of cosmic inflation---and to detect the imprints of the inflationary epoch, thereby dramatically increasing our understanding of fundamental physics and the early universe. A future CMB experiment with sufficient angular resolution and frequency coverage that surveys at least 1% of the sky to a depth of 1 uK-arcmin can deliver a constraint on the tensor-to-scalar ratio that will either result in a 5-sigma measurement of the energy scale of inflation or rule out all large-field inflation models, even in the presence of foregrounds and the gravitational lensing B-mode signal. LSS experiments, particularly spectroscopic surveys such as the Dark Energy Spectroscopic Instrument, will complement the CMB effort by improving current constraints on running of the spectral index by up to a factor of four, improving constraints on curvature by a factor of ten, and providing non-Gaussianity constraints that are competitive with the current CMB bounds.

We introduce a Lagrangian-space Effective Field Theory (LEFT) formalism for the study of cosmological large scale structures. Unlike the previous Eulerian-space construction, it is naturally formulated as an effective field theory of extended objects in Lagrangian space. In LEFT the resulting finite size effects are described using a multipole expansion parameterized by a set of time dependent coefficients and organized in powers of the ratio of the wavenumber of interest $k$ over the non-linear scale $k_{\rm NL}$. The multipoles encode the effects of the short distance modes on the long-wavelength universe and absorb UV divergences when present. There are no IR divergences in LEFT. Some of the parameters that control the perturbative approach are not assumed to be small and can be automatically resummed. We present an illustrative one-loop calculation for a power law universe. We describe the dynamics both at the level of the equations of motion and through an action formalism.

In single-field models of inflation the effect of a long mode with momentum q reduces to a diffeomorphism at zeroth and first order in q. This gives the well-known consistency relations for the n-point functions. At order q^2 the long mode has a physical effect on the short ones, since it induces curvature, and we expect that this effect is the same as being in a curved FRW universe. In this paper we verify this intuition in various examples of the three-point function, whose behaviour at order q^2 can be written in terms of the power spectrum in a curved universe. This gives a simple alternative understanding of the level of non-Gaussianity in single-field models. Non-Gaussianity is always parametrically enhanced when modes freeze at a physical scale k_{ph, f} shorter than H: f_{NL} \sim (k_{ph, f}/H)^2.

We assess two potential signals of the formation of our universe by the decay of a false vacuum. Negative spatial curvature is one possibility, but the window for its detection is now small. However, another possible signal is a suppression of the CMB power spectrum at large angles. This arises from the steepening of the effective potential as it interpolates between a flat inflationary plateau and the high barrier separating us from our parent vacuum. We demonstrate that these two effects can be parametrically separated in angular scale. Observationally, the steepening effect appears to be excluded at large l; but it remains consistent with the slight lack of power below l about 30 found by the WMAP and Planck collaborations. We give two simple models which improve the fit to the Planck data; one with observable curvature and one without. Despite cosmic variance, we argue that future CMB polarization and most importantly large-scale structure observations should be able to corroborate the Planck anomaly if it is real. If we further assume the specific theoretical setting of a landscape of metastable vacua, as suggested by string theory, we can estimate the probability of seeing a low-l suppression in the CMB. There are significant theoretical uncertainties in such calculations, but we argue the probability for a detectable suppression may be as large as O(1), and in general is significantly larger than the probability of seeing curvature.

We present algorithms for searching for azimuthally symmetric features in CMB data. Our algorithms are fully optimal for masked all-sky data with inhomogeneous noise, computationally fast, simple to implement, and make no approximations. We show how to implement the optimal analysis in both Bayesian and frequentist cases. In the Bayesian case, our algorithm for evaluating the posterior likelihood is so fast that we can do a brute-force search over parameter space, rather than using a Monte Carlo Markov chain. Our motivating example is searching for bubble collisions, a pre-inflationary signal which can be generated if multiple tunneling events occur in an eternally inflating spacetime, but our algorithms are general and should be useful in other contexts.

An appealing theory is that our current patch of universe was born as a nucleation bubble from a phase of false vacuum eternal inflation. We search for evidence for this theory by looking for the signal imprinted on the CMB that is generated when another bubble "universe" collides with our own. We create an efficient and optimal estimator for the signal in the WMAP 7-year data. We find no detectable signal, and constrain the amplitude, a, of the initial curvature perturbation that would be generated by a collision: -4.66 \times 10^{-8} < a (\sin{\thetabubble})^{4/3} < 4.73 \times 10^{-8} [Mpc^{-1}] at 95% confidence where \thetabubble is the angular radius of the bubble signal.

We apply the Effective Field Theory of Inflation to study the case where the continuous shift symmetry of the Goldstone boson \pi is softly broken to a discrete subgroup. This case includes and generalizes recently proposed String Theory inspired models of Inflation based on Axion Monodromy. The models we study have the property that the 2-point function oscillates as a function of the wavenumber, leading to oscillations in the CMB power spectrum. The non-linear realization of time diffeomorphisms induces some self-interactions for the Goldstone boson that lead to a peculiar non-Gaussianity whose shape oscillates as a function of the wavenumber. We find that in the regime of validity of the effective theory, the oscillatory signal contained in the n-point correlation functions, with n>2, is smaller than the one contained in the 2-point function, implying that the signature of oscillations, if ever detected, will be easier to find first in the 2-point function, and only then in the higher order correlation functions. Still the signal contained in higher-order correlation functions, that we study here in generality, could be detected at a subleading level, providing a very compelling consistency check for an approximate discrete shift symmetry being realized during inflation.

In single clock models of inflation the coupling between modes of very different scales does not have any significant dynamical effect during inflation. It leads to interesting projection effects. Larger and smaller modes change the relation between the scale a mode of interest will appear in the post-inflationary universe and will also change the time of horizon crossing of that mode. We argue that there are no infrared projection effects in physical questions, that there are no effects from modes of longer wavelength than the one of interest. These potential effects cancel when computing fluctuations as a function of physically measurable scales. Modes on scales smaller than the one of interest change the mapping between horizon crossing time and scale. The correction to the mapping computed in the absence of fluctuations is enhanced by a factor N_e, the number of e-folds of inflation between horizon crossing and reheating. The new mapping is stochastic in nature but its variance is not enhanced by N_e.

We analyze the signatures of inflationary models that are coupled to strongly interacting field theories, a basic class of multifield models also motivated by their role in providing dynamically small scales. Near the squeezed limit of the bispectrum, we find a simple scaling behavior determined by operator dimensions, which are constrained by the appropriate unitarity bounds. Specifically, we analyze two simple and calculable classes of examples: conformal field theories (CFTs), and large-N CFTs deformed by relevant time-dependent double-trace operators. Together these two classes of examples exhibit a wide range of scalings and shapes of the bispectrum, including nearly equilateral, orthogonal and local non-Gaussianity in different regimes. Along the way, we compare and contrast the shape and amplitude with previous results on weakly coupled fields coupled to inflation. This signature provides a precision test for strongly coupled sectors coupled to inflation via irrelevant operators suppressed by a high mass scale up to 1000 times the inflationary Hubble scale.

Studying loop corrections to inflationary perturbations, with particular emphasis on infrared factors, is important to understand the consistency of the inflationary theory, its predictivity and to establish the existence of the slow-roll eternal inflation phenomena and its recently found volume bound. In this paper we show that \zeta-correlators are time-independent at large distances at all-loop level in single clock inflation. We write the n-th order correlators of \dot\zeta\ as the time-integral of Green's functions times the correlators of local sources that are function of the lower order fluctuations. The Green's functions are such that only non-vanishing correlators of the sources at late times can lead to non-vanishing correlators for \dot\zeta\ at long distances. When the sources are connected by high wavenumber modes, the correlator is peaked at short distances, and these diagrams cannot lead to a time-dependence by simple diff. invariance arguments. When the sources are connected by long wavenumber modes one can use similar arguments once the constancy of \zeta\ at lower orders was established. Therefore the conservation of \zeta\ at a given order follows from the conservation of \zeta\ at the lower orders. Since at tree-level \zeta\ is constant, this implies constancy at all-loops by induction.

Large scale structure surveys will likely become the next leading cosmological probe. In our universe, matter perturbations are large on short distances and small at long scales, i.e. strongly coupled in the UV and weakly coupled in the IR. To make precise analytical predictions on large scales, we develop an effective field theory formulated in terms of an IR effective fluid characterized by several parameters, such as speed of sound and viscosity. These parameters, determined by the UV physics described by the Boltzmann equation, are measured from N-body simulations. We find that the speed of sound of the effective fluid is c_s^2 10^(-6) and that the viscosity contributions are of the same order. The fluid describes all the relevant physics at long scales k and permits a manifestly convergent perturbative expansion in the size of the matter perturbations \delta(k) for all the observables. As an example, we calculate the correction to the power spectrum at order \delta(k)^4. The predictions of the effective field theory are found to be in much better agreement with observation than standard cosmological perturbation theory, already reaching percent precision at this order up to a relatively short scale k \sim 0.24 h/Mpc.

We generalize the Effective Field Theory of Inflation to include additional light scalar degrees of freedom that are in their vacuum at the time the modes of interest are crossing the horizon. In order to make the scalars light in a natural way we consider the case where they are the Goldstone bosons of a global symmetry group or are partially protected by an approximate supersymmetry. We write the most general Lagrangian that couples the scalar mode associated to the breaking of time translation during inflation to the additional light scalar fields. This Lagrangian is constrained by diffeomorphism invariance and the additional symmetries that keep the new scalars light. This Lagrangian describes the fluctuations around the time of horizon crossing and it is supplemented with a general parameterization describing how the additional fluctuating fields can affect cosmological perturbations. We find that multifield inflation can reproduce the non-Gaussianities that can be generated in single field inflation but can also give rise to new kinds of non-Gaussianities. We find several new three-point function shapes. We show that in multifield inflation it is possible to naturally suppress the three-point function making the four-point function the leading source of detectable non-Gaussianities. We find that under certain circumstances, i.e. if specific shapes of non-Gaussianities are detected in the data, one could distinguish between single and multifield inflation and sometimes even among the various mechanisms that kept the additional fields light.

Studying loop corrections to inflationary perturbations, with particular emphasis on infrared factors, is important to understand the consistency of the inflationary theory, its predictivity and to establish the existence of the slow-roll eternal inflation phenomena and its recently found volume bound. In this paper we prove that the zeta correlation function is time-independent at one-loop level in single clock inflation. While many of the one-loop diagrams lead to a time-dependence when considered individually, the time-dependence beautifully cancels out in the overall sum. We identify two subsets of diagrams that cancel separately due to different physical reasons. The first cancellation is related to the change of the background cosmology due to the renormalization of the stress tensor. It results in a cancellation between the non-1PI diagrams and some of the diagrams made with quartic vertices. The second subset of diagrams that cancel is made up of cubic operators, plus the remaining quartic ones. We are able to write the sum of these diagrams as the integral over a specific three-point function between two very short wavelengths and one very long one. We then apply the consistency condition for this three-point function in the squeezed limit to show that the sum of these diagrams cannot give rise to a time dependence. This second cancellation is thus a consequence of the fact that in single clock inflation the attractor nature of the solution implies that a long wavelength zeta perturbation is indistinguishable from a trivial rescaling of the background, and so results in no physical effect on short wavelength modes.

We show that the squeezed limit of (N+1)-point functions of primordial correlation functions in which one of the modes has a very small wavenumber can be inferred from the spatial variation of locally measured N-point function. We then show how in single clock inflation a long wavelength perturbation can be re-absorbed in the background cosmology and how in computing correlation functions the integrals of the interaction Hamiltonian are dominated by conformal times of order of the short wavelength modes, when the long mode is already outside of the horizon. This allows us to generalize the consistency condition for N-point functions to the case in which the short wavelength fluctuations are inside the horizon and derivatives acts on them. We further discuss the consistency condition in the soft internal squeezed limit in which in an (N+M)-point function with (N+M) short modes the sum of the first N modes is a very soft momentum. These results are very useful to study infrared effects in Inflation.

We generalize the effective field theory of single clock inflation to include dissipative effects. Working in unitary gauge we couple a set of composite operators in the effective action which is constrained solely by invariance under time-dependent spatial diffeomorphisms. We restrict ourselves to situations where the degrees of freedom responsible for dissipation do no contribute to the density perturbations at late time. The dynamics of the perturbations is then modified by the appearance of `friction' and noise terms, and assuming certain locality properties for the Green's functions of these composite operators, we show that there is a regime characterized by a large friction term \gamma >> H in which the \zeta-correlators are dominated by the noise and the power spectrum can be significantly enhanced. We also compute the three point function <\zeta\zeta\zeta> for a wide class of models and discuss under which circumstances large friction leads to an increased level of non-Gaussianities. In particular, under our assumptions, we show that strong dissipation together with the required non-linear realization of the symmetries implies |f_NL| ~ \gamma/(c_s^2H) >> 1. As a paradigmatic example we work out a variation of the `trapped inflation' scenario with local response functions and perform the matching with our effective theory. A detection of the generic type of signatures that result from incorporating dissipative effects during inflation, as we describe here, would teach us about the dynamics of the early universe and also extend the parameter space of inflationary models.

It has recently been shown that in single field slow-roll inflation the total volume cannot grow by a factor larger than e^(S_dS/2) without becoming infinite. The bound is saturated exactly at the phase transition to eternal inflation where the probability to produce infinite volume becomes non zero. We show that the bound holds sharply also in any space-time dimensions, when arbitrary higher-dimensional operators are included and in the multi-field inflationary case. The relation with the entropy of de Sitter and the universality of the bound strengthen the case for a deeper holographic interpretation. As a spin-off we provide the formalism to compute the probability distribution of the volume after inflation for generic multi-field models, which might help to address questions about the population of vacua of the landscape during slow-roll inflation.

Length scales probed by large scale structure surveys are becoming closer to the horizon scale. Further, it has been recently understood that non-Gaussianity in the initial conditions could show up in a scale dependence of the bias of galaxies at the largest distances. It is therefore important to include General Relativistic effects. Here we provide a General Relativistic generalization of the bias, valid both for Gaussian and non-Gaussian initial conditions. The collapse of objects happens on very small scales, while long-wavelength modes are always in the quasi linear regime. Around every collapsing region, it is therefore possible to find a reference frame that is valid for all times and where the space time is almost flat: the Fermi frame. Here the Newtonian approximation is applicable and the equations of motion are the ones of the N-body codes. The effects of long-wavelength modes are encoded in the mapping from the cosmological frame to the local frame. For the linear bias, the effect of the long-wavelength modes on the dynamics is encoded in the local curvature of the Universe, which allows us to define a General Relativistic generalization of the bias in the standard Newtonian setting. We show that the bias due to this effect goes to zero as the squared ratio of the physical wavenumber with the Hubble scale for modes longer than the horizon, as modes longer than the horizon have no dynamical effects. However, the bias due to non-Gaussianities does not need to vanish for modes longer than the Hubble scale, and for non-Gaussianities of the local kind it goes to a constant. As a further application, we show that it is not necessary to perform large N-body simulations to extract information on long-wavelength modes: N-body simulations can be done on small scales and long-wavelength modes are encoded simply by adding curvature to the simulation and rescaling the coordinates.

The bispectrum vanishes for linear Gaussian fields and is thus a sensitive probe of non-linearities and non-Gaussianities in the cosmic density field. Hence, a detection of the bispectrum in the halo density field would enable tight constraints on non-Gaussian processes in the early Universe and allow inference of the dynamics driving inflation. We present a tree level derivation of the halo bispectrum arising from non-linear clustering, non-linear biasing and primordial non-Gaussianity. A diagrammatic description is developed to provide an intuitive understanding of the contributing terms and their dependence on scale, shape and the non-Gaussianity parameter fNL. We compute the terms based on a multivariate bias expansion and the peak-background split method and show that non-Gaussian modifications to the bias parameters lead to amplifications of the tree level bispectrum that were ignored in previous studies. Our results are in a good agreement with published simulation measurements of the halo bispectrum. Finally, we estimate the expected signal to noise on fNL and show that the constraint obtainable from the bispectrum analysis significantly exceeds the one obtainable from the power spectrum analysis.

The effective theory of adiabatic fluctuations around arbitrary Friedmann-Robertson-Walker backgrounds - both expanding and contracting - allows for more than one way to obtain scale-invariant two-point correlations. However, as we show in this paper, it is challenging to produce scale-invariant fluctuations that are weakly coupled over the range of wavelengths accessible to cosmological observations. In particular, requiring the background to be a dynamical attractor, the curvature fluctuations are scale-invariant and weakly coupled for at least 10 e-folds only if the background is close to de Sitter space. In this case, the time-translation invariance of the background guarantees time-independent n-point functions. For non-attractor solutions, any predictions depend on assumptions about the evolution of the background even when the perturbations are outside of the horizon. For the simplest such scenario we identify the regions of the parameter space that avoid both classical and quantum mechanical strong coupling problems. Finally, we present extensions of our results to backgrounds in which higher-derivative terms play a significant role.

Using the recently developed effective field theory of inflation, we argue that the size and the shape of the non-Gaussianities generated by single-field inflation are generically well described by two parameters: f_NL^equil, which characterizes the size of the signal that is peaked on equilateral configurations, and f_NL^orthog, which instead characterizes the size of the signal which is peaked both on equilateral configurations and flat-triangle configurations (with opposite signs). The shape of non-Gaussianities associated with f_NL^orthog is orthogonal to the one associated to f_NL^equil, and former analysis have been mostly blind to it. We perform the optimal analysis of the WMAP 5-year data for both of these parameters. We find no evidence of non-Gaussianity, and we have the following constraints: -125 < f_NL^equil < 435, -369 < f_NL^orthog < 71 at 95% CL. We show that both of these constraints can be translated into limits on parameters of the Lagrangian of single-field inflation. For one of them, the speed of sound of the inflaton fluctuations, we find that it is either bounded to be c_s > 0.011 at 95% CL. or alternatively to be so small that the higher-derivative kinetic term dominate at horizon crossing. We are able to put similar constraints on the other operators of the inflaton Lagrangian.

The universe is smooth on large scales but very inhomogeneous on small scales. Why is the spacetime on large scales modeled to a good approximation by the Friedmann equations? Are we sure that small-scale non-linearities do not induce a large backreaction? Related to this, what is the effective theory that describes the universe on large scales? In this paper we make progress in addressing these questions. We show that the effective theory for the long-wavelength universe behaves as a viscous fluid coupled to gravity: integrating out short-wavelength perturbations renormalizes the homogeneous background and introduces dissipative dynamics into the evolution of long-wavelength perturbations. The effective fluid has small perturbations and is characterized by a few parameters like an equation of state, a sound speed and a viscosity parameter. These parameters can be matched to numerical simulations or fitted from observations. We find that the backreaction of small-scale non-linearities is very small, being suppressed by the large hierarchy between the scale of non-linearities and the horizon scale. The effective pressure of the fluid is always positive and much too small to significantly affect the background evolution. Moreover, we prove that virialized scales decouple completely from the large-scale dynamics, at all orders in the post-Newtonian expansion. We propose that our effective theory be used to formulate a well-defined and controlled alternative to conventional perturbation theory, and we discuss possible observational applications. Finally, our way of reformulating results in second-order perturbation theory in terms of a long-wavelength effective fluid provides the opportunity to understand non-linear effects in a simple and physically intuitive way.

Non-Gaussianities of the primordial density perturbations have emerged as a very powerful possible signal to test the dynamics that drove the period of inflation. While in general the most sensitive observable is the three-point function in this paper we show that there are technically natural inflationary models where the leading source of non-Gaussianity is the four-point function. Using the recently developed Effective Field Theory of Inflation, we are able to show that it is possible to impose an approximate parity symmetry and an approximate continuos shift symmetry on the inflaton fluctuations that allow, when the dispersion relation is of the form $\omega\sim c_s k$, for a unique quartic operator, while approximately forbidding all the cubic ones. The resulting shape for the four-point function is unique. In the models where the dispersion relation is of the form $\omega\sim k^2/M$ a similar construction can be carried out and additional shapes are possible.

We study the spherical collapse model in the presence of quintessence with negligible speed of sound. This case is particularly motivated for w<-1 as it is required by stability. As pressure gradients are negligible, quintessence follows dark matter during the collapse. The spherical overdensity behaves as a separate closed FLRW universe, so that its evolution can be studied exactly. We derive the critical overdensity for collapse and we use the extended Press-Schechter theory to study how the clustering of quintessence affects the dark matter mass function. The effect is dominated by the modification of the linear dark matter growth function. A larger effect occurs on the total mass function, which includes the quintessence overdensities. Indeed, here quintessence constitutes a third component of virialized objects, together with baryons and dark matter, and contributes to the total halo mass by a fraction ~ (1+w) Omega_Q / Omega_m. This gives a distinctive modification of the total mass function at low redshift.

We study loop corrections to correlation functions of inflationary perturbations. Previous calculations have found that the two-point function can have a logarithmic running of the form log(k/mu), where k is the wavenumber of the perturbation, and mu is the renormalization scale. We highlight that this result would have profound consequences for both eternal inflation and the predictivity of standard inflation. We find a different result. We consider two sets of theories: one where the inflaton has a large cubic self-interaction and one where the inflaton interacts gravitationally with N massless spectator scalar fields. We find that there is a logarithmic running but of the form log(H/mu), where H is the Hubble constant during inflation. We find this result in three independent ways: by performing the calculation with a sharp cutoff in frequency-momentum space, in dimensional regularization and by the simple procedure of making the loop integral dimensionless. For the simplest of our theories we explicitly renormalize the correlation function proving that the divergencies can be reabsorbed and that the correlation function for super-horizon modes does not depend on time (once the tadpole terms have been properly taken into account). We prove the time-independence of the super-horizon correlation function in several additional ways: by doing the calculation of the correlation function at finite time using both the regularizations and by developing a formalism which expresses loop corrections directly in terms of renormalized quantities at each time. We find this last formalism particularly helpful to develop intuition which we then use to generalize our results to higher loops and different interactions.

We have applied the optimal estimator for f_{NL}^{local} to the 5 year WMAP data. Marginalizing over the amplitude of foreground templates we get -4 < f_{NL}^{local} < 80 at 95% CL. Error bars of previous (sub-optimal) analyses are roughly 40% larger than these. The probability that a Gaussian simulation, analyzed using our estimator, gives a result larger in magnitude than the one we find is 7%. Our pipeline gives consistent results when applied to the three and five year WMAP data releases and agrees well with the results from our own sub-optimal pipeline. We find no evidence of any residual foreground contamination.

We approximately compute the bispectrum induced on the CMB temperature by fluctuations in the standard recombination epoch. Of all the second order sources that can induce non-Gaussianity during recombination, we concentrate on those proportional to the perturbation in the free electron density, which is about a factor of 5 larger than the other first order perturbations. This term induces some non-Gaussianity by delaying the time of recombination and by changing the photon diffusion scale. We find that the signal is not scale invariant, peaked on squeezed triangles with the smaller multipole around the scale of the first acoustic peak, and that its size corresponds to an effective f_NL ~ -3.5, which could be marginally detected by Planck if both temperature and polarization are measured.

A new and powerful probe of the origin and evolution of structures in the Universe has emerged and been actively developed over the last decade. In the coming decade, non-Gaussianity, i.e., the study of non-Gaussian contributions to the correlations of cosmological fluctuations, will become an important probe of both the early and the late Universe. Specifically, it will play a leading role in furthering our understanding of two fundamental aspects of cosmology and astrophysics: (i) the physics of the very early universe that created the primordial seeds for large-scale structures, and (ii) the subsequent growth of structures via gravitational instability and gas physics at later times. To date, observations of fluctuations in the Cosmic Microwave Background (CMB) and the Large-Scale Structure of the Universe (LSS) have focused largely on the Gaussian contribution as measured by the two-point correlations (or the power spectrum) of density fluctuations. However, an even greater amount of information is contained in non-Gaussianity and a large discovery space therefore still remains to be explored. Many observational probes can be used to measure non-Gaussianity, including CMB, LSS, gravitational lensing, Lyman-alpha forest, 21-cm fluctuations, and the abundance of rare objects such as clusters of galaxies and high-redshift galaxies. Not only does the study of non-Gaussianity maximize the science return from a plethora of present and future cosmological experiments and observations, but it also carries great potential for important discoveries in the coming decade.

We summarize the utility of precise cosmic microwave background (CMB) polarization measurements as probes of the physics of inflation. We focus on the prospects for using CMB measurements to differentiate various inflationary mechanisms. In particular, a detection of primordial B-mode polarization would demonstrate that inflation occurred at a very high energy scale, and that the inflaton traversed a super-Planckian distance in field space. We explain how such a detection or constraint would illuminate aspects of physics at the Planck scale. Moreover, CMB measurements can constrain the scale-dependence and non-Gaussianity of the primordial fluctuations and limit the possibility of a significant isocurvature contribution. Each such limit provides crucial information on the underlying inflationary dynamics. Finally, we quantify these considerations by presenting forecasts for the sensitivities of a future satellite experiment to the inflationary parameters.

How did the universe evolve? The fine angular scale (l>1000) temperature and polarization anisotropies in the CMB are a Rosetta stone for understanding the evolution of the universe. Through detailed measurements one may address everything from the physics of the birth of the universe to the history of star formation and the process by which galaxies formed. One may in addition track the evolution of the dark energy and discover the net neutrino mass. We are at the dawn of a new era in which hundreds of square degrees of sky can be mapped with arcminute resolution and sensitivities measured in microKelvin. Acquiring these data requires the use of special purpose telescopes such as the Atacama Cosmology Telescope (ACT), located in Chile, and the South Pole Telescope (SPT). These new telescopes are outfitted with a new generation of custom mm-wave kilo-pixel arrays. Additional instruments are in the planning stages.

We analyze a distinctive mechanism for inflation in which particle production slows down a scalar field on a steep potential, and show how it descends from angular moduli in string compactifications. The analysis of density perturbations -- taking into account the integrated effect of the produced particles and their quantum fluctuations -- requires somewhat new techniques that we develop. We then determine the conditions for this effect to produce sixty e-foldings of inflation with the correct amplitude of density perturbations at the Gaussian level, and show that these requirements can be straightforwardly satisfied. Finally, we estimate the amplitude of the non-Gaussianity in the power spectrum and find a significant equilateral contribution.

We consider the contribution to the three-point function of matter density fluctuations from nonlinear growth after modes re-enter the horizon, and discuss effects that must be included in order to predict the three-point function with an accuracy comparable to primordial nongaussianities with f_NL ~ few. In particular, we note that the shortest wavelength modes measured in galaxy surveys entered the horizon during the radiation era, and, as a result, the radiation era modifies their three-point function by a magnitude equivalent to f_NL ~ O(4). On longer wavelengths, where the radiation era is negligible, we find that the corrections to the nonlinear growth from relativistic effects become important at the level f_NL ~ few. We implement a simple method for numerically calculating the three-point function, by solving the second-order equations of motion for the perturbations with the first order perturbations providing a source.

We derive the full set of second-order equations governing the evolution of cosmological perturbations, including the effects of the first-order electron number density perturbations, \delta_e. We provide a detailed analysis of the perturbations to the recombination history of the universe and show that a perturbed version of the Peebles effective 3-level atom is sufficient for obtaining the evolution of \delta_e for comoving wavenumbers smaller than 1Mpc^{-1}. We calculate rigorously the perturbations to the Ly\alpha escape probability and show that to a good approximation it is governed by the local baryon velocity divergence. For modes shorter than the photon diffusion scale, we find that \delta_e is enhanced during recombination by a factor of roughly 5 relative to other first-order quantities sourcing the CMB anisotropies at second order. Using these results, in a companion paper we calculate the CMB bispectrum generated during recombination.

We calculate the probability distribution for the volume of the Universe after slow-roll inflation both in the eternal and in the non-eternal regime. Far from the eternal regime the probability distribution for the number of e-foldings, defined as one third of the logarithm of the volume, is sharply peaked around the number of e-foldings of the classical inflaton trajectory. At the transition to the eternal regime this probability is still peaked (with the width of order one e-folding) around the average, which gets twice larger at the transition point. As one enters the eternal regime the probability for the volume to be finite rapidly becomes exponentially small. In addition to developing techniques to study eternal inflation, our results allow us to establish the quantum generalization of a recently proposed bound on the number of e-foldings in the non-eternal regime: the probability for slow-roll inflation to produce a finite volume larger than e^(S_dS/2), where S_dS is the de Sitter entropy at the end of the inflationary stage, is smaller than the uncertainty due to non-perturbative quantum gravity effects. The existence of such a bound provides a consistency check for the idea of de Sitter complementarity.

For slow-roll inflation we study the phase transition to the eternal regime. Starting from a finite inflationary volume, we consider the volume of the universe at reheating as order parameter. We show that there exists a critical value for the classical inflaton speed, \dot\phi^2/H^4 = 3/(2 \pi^2), where the probability distribution for the reheating volume undergoes a sharp transition. In particular, for sub-critical inflaton speeds all distribution moments become infinite. We show that at the same transition point the system develops a non-vanishing probability of having a strictly infinite reheating volume, while retaining a finite probability for finite values. Our analysis represents the exact quantum treatment of the system at lowest order in the slow-roll parameters and H^2/M_Pl^2.

Even if nothing but a light Higgs is observed at the LHC, suggesting that the Standard Model is unmodified up to scales far above the weak scale, Higgs physics can yield surprises of fundamental significance for cosmology. As has long been known, the Standard Model vacuum may be metastable for low enough Higgs mass, but a specific value of the decay rate holds special significance: for a very narrow window of parameters, our Universe has not yet decayed but the current inflationary period can not be future eternal. Determining whether we are in this window requires exquisite but achievable experimental precision, with a measurement of the Higgs mass to 0.1 GeV at the LHC, the top mass to 60 MeV at a linear collider, as well as an improved determination of alpha_s by an order of magnitude on the lattice. If the parameters are observed to lie in this special range, particle physics will establish that the future of our Universe is a global big crunch, without harboring pockets of eternal inflation, strongly suggesting that eternal inflation is censored by the fundamental theory. This conclusion could be drawn even more sharply if metastability with the appropriate decay rate is found in the MSSM, where the physics governing the instability can be directly probed at the TeV scale.

We study the effective field theory of inflation, i.e. the most general theory describing the fluctuations around a quasi de Sitter background, in the case of single field models. The scalar mode can be eaten by the metric by going to unitary gauge. In this gauge, the most general theory is built with the lowest dimension operators invariant under spatial diffeomorphisms, like g^{00} and K_{mu nu}, the extrinsic curvature of constant time surfaces. This approach allows us to characterize all the possible high energy corrections to simple slow-roll inflation, whose sizes are constrained by experiments. Also, it describes in a common language all single field models, including those with a small speed of sound and Ghost Inflation, and it makes explicit the implications of having a quasi de Sitter background. The non-linear realization of time diffeomorphisms forces correlation among different observables, like a reduced speed of sound and an enhanced level of non-Gaussianity.

The consistency relation for the 3-point function of the CMB is a very powerful observational signature which is believed to be true for every inflationary model in which there is only one dynamical degree of freedom. Its importance relies on the fact that deviations from it might be detected in next generation experiments, allowing us to rule out all single field inflationary models. After making more precise the already existing proof of the consistency relation, we use a recently developed effective field theory for inflationary perturbations to provide an alternative and very explicit proof valid at leading non trivial order in slow roll parameters.

We present a bouncing cosmology which evolves from the contracting to the expanding phase in a smooth way, without developing instabilities or pathologies and remaining in the regime of validity of 4d effective field theory. A nearly scale invariant spectrum of perturbations is generated during the contracting phase by an isocurvature scalar with a negative exponential potential and then converted to adiabatic. The model predicts a slightly blue spectrum, n_S >~ 1, no observable gravitational waves and a high (but model dependent) level of non-Gaussianities with local shape. The model represents an explicit and predictive alternative to inflation, although, at present, it is clearly less compelling.

We study the Likelihood function of data given f_NL for the so-called local type of non-Gaussianity. In this case the curvature perturbation is a non-linear function, local in real space, of a Gaussian random field. We compute the Cramer-Rao bound for f_NL and show that for small values of f_NL the 3-point function estimator saturates the bound and is equivalent to calculating the full Likelihood of the data. However, for sufficiently large f_NL, the naive 3-point function estimator has a much larger variance than previously thought. In the limit in which the departure from Gaussianity is detected with high confidence, error bars on f_NL only decrease as 1/ln Npix rather than Npix^-1/2 as the size of the data set increases. We identify the physical origin of this behavior and explain why it only affects the local type of non-Gaussianity, where the contribution of the first multipoles is always relevant. We find a simple improvement to the 3-point function estimator that makes the square root of its variance decrease as Npix^-1/2 even for large f_NL, asymptotically approaching the Cramer-Rao bound. We show that using the modified estimator is practically equivalent to computing the full Likelihood of f_NL given the data. Thus other statistics of the data, such as the 4-point function and Minkowski functionals, contain no additional information on f_NL. In particular, we explicitly show that the recent claims about the relevance of the 4-point function are not correct. By direct inspection of the Likelihood, we show that the data do not contain enough information for any statistic to be able to constrain higher order terms in the relation between the Gaussian field and the curvature perturbation, unless these are orders of magnitude larger than the size suggested by the current limits on f_NL.

We analyze the 3-year WMAP data and look for a deviation from Gaussianity in the form of a 3-point function that has either of the two theoretically motivated shapes: local and equilateral. There is no evidence of departure from Gaussianity and the analysis gives the presently tightest bounds on the parameters f_ NL^local and f_NL^equil., which define the amplitude of respectively the local and the equilateral non-Gaussianity: -36 < f_NL^local < 100, -256 < f_NL^equil. < 332 at 95% C.L.

We present a consistent effective theory that violates the null energy condition (NEC) without developing any instabilities or other pathological features. The model is the ghost condensate with the global shift symmetry softly broken by a potential. We show that this system can drive a cosmological expansion with dH/dt > 0. Demanding the absence of instabilities in this model requires dH/dt <~ H^2. We then construct a general low-energy effective theory that describes scalar fluctuations about an arbitrary FRW background, and argue that the qualitative features found in our model are very general for stable systems that violate the NEC. Violating the NEC allows dramatically non-standard cosmological histories. To illustrate this, we construct an explicit model in which the expansion of our universe originates from an asymptotically flat state in the past, smoothing out the big-bang singularity within control of a low-energy effective theory. This gives an interesting alternative to standard inflation for solving the horizon problem. We also construct models in which the present acceleration has w < -1; a periodic ever-expanding universe and a model with a smooth ``bounce'' connecting a contracting and expanding phase.

We study a recently proposed mechanism to solve the hierarchy problem in the context of the landscape, where the solution of the hierarchy problem is connected to the requirement of having baryons in our universe via Electroweak Baryogenesis. The phase transition is triggered by the fermion condensation of a new gauge sector which becomes strong at a scale Lambda determined by dimensional transmutation, and it is mediated to the standard model by a new singlet field. In a ``friendly'' neighborhood of the landscape, where only the relevant operators are ``scanned'' among the vacua, baryogenesis is effective only if the higgs mass m_h is comparable to this low scale Lambda, forcing m_h to be of order Lambda, and solving the hierarchy problem. A new CP violating phase is needed coupling the new singlet and the higgs field to new matter fields. We study the constraints on this model given by baryogenesis and by the electron electric dipole moment (EDM), and we briefly comment on gauge coupling unification and on dark matter relic abundance. We find that next generation experiments on the EDM will be sensitive to essentially the entire viable region of the parameter space, so that absence of a signal would effectively rule out the model.

Gauge coupling unification and the success of TeV-scale weakly interacting dark matter are usually taken as evidence of low energy supersymmetry (SUSY). However, if we assume that the tuning of the higgs can be explained in some unnatural way, from environmental considerations for example, SUSY is no longer a necessary component of any Beyond the Standard Model theory. In this paper we study the minimal model with a dark matter candidate and gauge coupling unification. This consists of the SM plus fermions with the quantum numbers of SUSY higgsinos, and a singlet. It predicts thermal dark matter with a mass that can range from 100 GeV to around 2 TeV and generically gives rise to an electric dipole moment that is just beyond current experimental limits, with a large portion of its allowed parameter space accessible to next generation EDM and direct detection experiments. We study precision unification in this model by embedding it in a 5-D orbifold GUT where certain large threshold corrections are calculable, achieving gauge coupling and b-tau unification, and predicting a rate of proton decay just beyond current limits.

We develop a method to constrain the level of non-Gaussianity of density perturbations when the 3-point function is of the "equilateral" type. Departures from Gaussianity of this form are produced by single field models such as ghost or DBI inflation and in general by the presence of higher order derivative operators in the effective Lagrangian of the inflaton. We show that the induced shape of the 3-point function can be very well approximated by a factorizable form, making the analysis practical. We also show that, unless one has a full sky map with uniform noise, in order to saturate the Cramer-Rao bound for the error on the amplitude of the 3-point function, the estimator must contain a piece that is linear in the data. We apply our technique to the WMAP data obtaining a constraint on the amplitude f_NL^equil of "equilateral" non-Gaussianity: -366 < f_NL^equil < 238 at 95% C.L. We also apply our technique to constrain the so-called "local" shape, which is predicted for example by the curvaton and variable decay width models. We show that the inclusion of the linear piece in the estimator improves the constraint over those obtained by the WMAP team, to -27 < f_NL^local < 121 at 95% C.L.

In a ghost inflationary scenario, we study the observational consequences of a tilt in the potential of the ghost condensate. We show how the presence of a tilt tends to make contact between the natural predictions of ghost inflation and the ones of slow roll inflation. In the case of positive tilt, we are able to build an inflationary model in which the Hubble constant H is growing with time. We compute the amplitude and the tilt of the 2-point function, as well as the 3-point function, for both the cases of positive and negative tilt. We find that a good fraction of the parameter space of the model is within experimental reach.

In recently introduced Split Susy theories, in which the scale of Susy breaking is very high, the requirement that the relic abundance of the Lightest SuperPartner (LSP) provides the Dark Matter of the Universe leads to the prediction of fermionic superpartners around the weak scale. This is no longer obviously the case if the LSP is a hidden sector field, such as a Gravitino or an other hidden sector fermion, so, it is interesting to study this scenario. We consider the case in which the Next-Lightest SuperPartner (NLSP) freezes out with its thermal relic abundance, and then it decays to the LSP. We use the constraints from BBN and CMB, together with the requirement of attaining Gauge Coupling Unification and that the LSP abundance provides the Dark Matter of the Universe, to infer the allowed superpartner spectrum. As very good news for a possible detaction of Split Susy at LHC, we find that if the Gravitino is the LSP, than the only allowed NLSP has to be very purely photino like. In this case, a photino from 700 GeV to 5 TeV is allowed, which is difficult to test at LHC. We also study the case where the LSP is given by a light fermion in the hidden sector which is naturally present in Susy breaking in Extra Dimensions. We find that, in this case, a generic NLSP is allowed to be in the range 1-20 TeV, while a Bino NLSP can be as light as tens of GeV.

We show that the isocurvature perturbations imply that the gravitationally produced superheavy dark matter must have masses larger than few times the Hubble expansion rate at the end of inflation. This together with the bound on tensor to scalar contribution to the CMB induces a lower bound on the reheating temperature for superheavy dark matter to be about 10^7 GeV. Hence, if the superheavy dark matter scenario is embedded in supergravity models with gravity mediated SUSY breaking, the gravitino bound will squeeze this scenario. Furthermore, the CMB constraint strengthens the statement that gravitationally produced superheavy dark matter scenario prefers a relatively large tensor mode amplitude if the reheating temperature must be less than 10^9 GeV.