Characterizing the physical properties of the stochastic gravitational waves background (SGWB) is a key step towards identifying the nature of its possible origin. We focus our analysis on SGWB anisotropies. The existence of a non-trivial primordial scalar-tensor-tensor (STT) correlation in the squeezed configuration may be inferred from the effect that a long wavelength scalar mode has on the gravitational wave power spectrum: an anisotropic contribution. Crucially, such contribution is correlated with temperature anisotropies in the cosmic microwave background (CMB). We show that, for inflationary models that generate suitably large STT non-Gaussianities, cross-correlating the CMB with the stochastic background of gravitational waves is a very effective probe of early universe physics. The resulting signal can be a smoking-gun for primordial SGWB anisotropies.

With an energy scale that can be as high as $10^{14}\,{\rm GeV}$, inflation may provide a unique probe of high-energy physics. Both scalar and tensor fluctuations generated during this early accelerated expansion contain crucial information about the particle content of the primordial universe. The advent of ground- and space-based interferometers enables us to probe primordial physics at length-scales much smaller than those corresponding to current CMB constraints. One key prediction of single-field slow-roll inflation is a red-tilted gravitational wave spectrum, currently inaccessible at interferometer scales. Interferometers probe directly inflationary physics that deviates from the minimal scenario and in particular additional particle content with sizeable couplings to the inflaton field. We adopt here an effective description for such fields and focus on the case of extra spin-2 fields. We find that a time-dependent sound speed for the helicity-2 modes can generate primordial gravitational waves with a blue-tilted spectrum, potentially detectable at interferometer scales.

Modifications of General Relativity leave their imprint both on the cosmic expansion history through a non-trivial dark energy equation of state, and on the evolution of cosmological perturbations in the scalar and in the tensor sectors. In particular, the modification in the tensor sector gives rise to a notion of gravitational-wave (GW) luminosity distance, different from the standard electromagnetic luminosity distance, that can be studied with standard sirens at GW detectors such as LISA or third-generation ground based experiments. We discuss the predictions for modified GW propagation from some of the best studied theories of modified gravity, such as Horndeski or the more general degenerate higher order scalar-tensor (DHOST) theories, non-local infrared modifications of gravity, bigravity theories and the corresponding phenomenon of GW oscillation, as well as theories with extra or varying dimensions. We show that modified GW propagation is a completely generic phenomenon in modified gravity. We then use a simple parametrization of the effect in terms of two parameters $(\Xi_0,n)$, that is shown to fit well the results from a large class of models, to study the prospects of observing modified GW propagation using supermassive black hole binaries as standard sirens with LISA. We construct mock source catalogs and perform detailed Markov Chain Monte Carlo studies of the likelihood obtained from LISA standard sirens alone, as well as by combining them with CMB, BAO and SNe data to reduce the degeneracies between cosmological parameters. We find that the combination of LISA with the other cosmological datasets allows one to measure the parameter $\Xi_0$ that characterizes modified GW propagation to the percent level accuracy, sufficient to test several modified gravity theories. [Abridged]

Evidence for the presence of extra fields during inflation may be found in the anisotropies of the scalar and tensor spectra across a vast range of scales. Indeed, beyond the single-field slow-roll paradigm, a long tensor mode can modulate the power spectrum inducing a sizable quadrupolar anisotropy. We investigate how this dynamics plays out for the tensor two-point correlator. The resulting quadrupole stores information on squeezed tensor non-Gaussianities, specifically those sourced by the extra field content and responsible for the breaking of so-called consistency relations. We underscore the potential of anisotropies as a probe of new physics: testable at CMB scales through the detection of B-modes, they are accessible at smaller scales via interferometers and pulsar timing arrays.

Following the pioneering observations with COBE in the early 1990s, studies of the cosmic microwave background (CMB) have focused on temperature and polarization anisotropies. CMB spectral distortions - tiny departures of the CMB energy spectrum from that of a perfect blackbody - provide a second, independent probe of fundamental physics, with a reach deep into the primordial Universe. The theoretical foundation of spectral distortions has seen major advances in recent years, which highlight the immense potential of this emerging field. Spectral distortions probe a fundamental property of the Universe - its thermal history - thereby providing additional insight into processes within the cosmological standard model (CSM) as well as new physics beyond. Spectral distortions are an important tool for understanding inflation and the nature of dark matter. They shed new light on the physics of recombination and reionization, both prominent stages in the evolution of our Universe, and furnish critical information on baryonic feedback processes, in addition to probing primordial correlation functions at scales inaccessible to other tracers. In principle the range of signals is vast: many orders of magnitude of discovery space could be explored by detailed observations of the CMB energy spectrum. Several CSM signals are predicted and provide clear experimental targets, some of which are already observable with present-day technology. Confirmation of these signals would extend the reach of the CSM by orders of magnitude in physical scale as the Universe evolves from the initial stages to its present form. The absence of these signals would pose a huge theoretical challenge, immediately pointing to new physics.

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 next generation of instruments designed to measure the polarization of the cosmic microwave background (CMB) will provide a historic opportunity to open the gravitational wave window to the primordial Universe. Through high sensitivity searches for primordial gravitational waves, and tighter limits on the energy released in processes like phase transitions, the CMB polarization data of the next decade has the potential to transform our understanding of the laws of physics underlying the formation of the Universe.

We consider a system where dark matter dynamics is enriched by the presence of clustering quintessence in the approximation where the system is effectively reduced to one degree of freedom. We study the corresponding observables up to one-loop order and then point out similarities between the power spectrum of the reduced system and the behaviour of non-equal time pure dark matter correlators. We then focus on the one-loop total density power spectrum in the IR limit as a diagnostic tool for consistency relations breaking. Unlike the non-equal time case, the reduced system does still obey consistency relations; we illustrate this by explicitly verifying the 1-loop IR cancellation. A more general setup, obtained by relaxing the assumption of a vanishing sound speed, is also analyzed. In this and similar scenarios the presence of additional dynamics, typical of dark energy and modified gravity models, implies that one may no longer gauge away the squeezed contribution of observables such as the dark matter bispectrum. We show how these effects propagate all the way to biased tracers.

The perturbative treatment of dark matter in structure formation relies on the existence of a well-defined expansion parameter, $k/k_{\rm NL}$, with $k_{\rm NL}$ signalling the onset and ultimately the leading role of non-linearities in the system. Cosmologies beyond the {\Lambda}CDM model often come with additional degree(s) of freedom. The scale $k_{\rm V}$ at which non-linearities become important in the additional sector(s) can be rather different from $k_{\rm NL}$. For theories endowed with a Vainshtein-type screening mechanism, $k_{\rm V}$ sets the scale where screening becomes efficient and restores continuity with the predictions of general relativity. This is precisely the dynamics that allows such theories to pass existing observational tests at scales where general relativity has been tested with exquisite precision (e.g. solar system scales). We consider here the mildly-non-linear scales of a dark matter component coupled to a galileon-type field and focus in particular on the case of a $k_{\rm V}$ < $k_{\rm NL}$ hierarchy. We put forward a phenomenological framework that describes the effects of screening dynamics on large scale structure observables.

Tensor non-Gaussianity represents an important future probe of the physics of inflation. Inspired by recent works, we elaborate further on the possibility of significant primordial tensor non-Gaussianities sourced by extra fields during inflation. Unitarity constraints limit the impact of extra (spinning) particle content by means of a lower bound on the corresponding mass spectrum. For spin-2 particles, this takes the form of the well-known Higuchi bound. Massive ($m\gtrsim H$) particles will typically decay during inflation unless they are non-minimally coupled to the inflaton sector: the inflating field "lifts" the dynamics of the extra field(s), effectively getting around the limits imposed by unitarity. There exist several models that realize such a mechanism, but we focus here on the set-up of [1] where, through an EFT approach, one is able to capture the essential features common to an entire class of theories. In the presence of an extra massive spin-2 particle, the interactions in the tensor sector mimic very closely those in the scalar sector of quasi-single-field inflationary models. We calculate the tensor bispectrum in different configurations and extract its dependence on the extra tensor sound speed. We show in detail how one may obtain significant tensor non-Gaussianities whose shape-function interpolates between local and equilateral, depending on the mass of the extra field. We also estimate the LISA response functions to a tensor bispectrum supporting the intermediate-type shapes we find.

The detection and characterization of primordial gravitational waves through their impact on the polarization anisotropies of the cosmic microwave background (CMB) is a primary science goal of current and future observations of the CMB. An ancillary dataset that will become accessible with the great leaps in sensitivity of CMB experiments is the polarized Sunyaev Zel'dovich (pSZ) effect, small-scale CMB polarization anisotropies induced by scattering from free electrons in the post-reionization Universe. The cross correlation of the pSZ effect with galaxy surveys, a technique known as pSZ tomography, can be used to reconstruct the remote quadrupole field: the CMB quadrupole observed from different locations in the Universe. Primordial gravitational waves leave a distinct imprint on the remote quadrupole field, making pSZ tomography a potential new method to characterize their properties. Building on previous work, we explore the utility of the full set of correlations between the primary CMB and the reconstructed remote quadrupole field to both provide exclusion limits on the amplitude of primordial gravitational waves, as well as to provide constraints on several phenomenological models of the tensor sector: axion gauge field inflation, general models with chiral tensors, and models with modified late-time decay of tensors. We find that relatively futuristic experimental requirements are necessary to provide competitive exclusion limits compared with the primary CMB. However, pSZ tomography can be a powerful probe of the late-time evolution of tensors and, through cross-correlations with the primary CMB, can provide mild improvements on parameter constraints in various models with chiral primordial gravitational waves.

We study the scalar-tensor-tensor non-Gaussian signal in an inflationary model comprising also an axion coupled with SU(2) gauge fields. In this set-up, metric fluctuations are sourced by the gauge fields already at the linear level providing an enhanced chiral gravitational waves spectrum. The same mechanism is at work in generating an amplitude for the three-point function that is parametrically larger than in standard single-field inflation.

We study how inflationary observables associated with primordial tensor modes are affected by coupling the minimal field content with an extra spin-2 particle during inflation. We work with a model that is ghost-free at the fully non-linear level and show how the new degrees of freedom modify standard consistency relations for the tensor bispectrum. The extra interacting spin-2 field is necessarily massive and unitarity dictates its mass be in the $m \gtrsim H$ range. Despite the fact that this bound selects a decaying solution for the corresponding tensor mode, cosmological correlators still carry the imprints of such "fossil" fields. Remarkably, fossil(s) of spin $\geq 1$ generate distinctive anisotropies in observables such as the tensor power spectrum. We show how this plays out in our set-up.

We study the imprints of a massive spin-2 field on inflationary observables, and in particular on the breaking of consistency relations. In this setup, the minimal inflationary field content interacts with the massive spin-2 field through dRGT interactions, thus guaranteeing the absence of Boulware-Deser ghostly degrees of freedom. The unitarity requirement on spinning particles, known as Higuchi bound, plays a crucial role for the size of the observable signal.

A specific value for the cosmological constant, \Lambda, can account for late-time cosmic acceleration. However, motivated by the so-called cosmological constant problem(s), several alternative mechanisms have been explored. To date, a host of well-studied dynamical dark energy and modified gravity models exists. Going beyond \Lambda CDM often comes with additional degrees of freedom (dofs). For these to pass existing observational tests, an efficient screening mechanism must be in place. The linear and quasi-linear regimes of structure formation are ideal probes of such dofs and can capture the onset of screening. We propose here a semi-phenomenological treatment to account for screening dynamics on LSS observables, with special emphasis on Vainshtein-type screening.

The perturbative approach to structure formation has recently received a lot of attention in the literature. In such setups the final predictions for observables like the power spectrum is often derived under additional approximations such as a simplified time dependence. Here we provide all-order perturbative integral solutions for density and velocity fields in generalized cosmologies, with a direct application to clustering quintessence. We go beyond the standard results based on extending the EdS-like approximations. As an illustrative example, we apply our findings to the calculation of the one-loop power spectrum of density and momentum fields. We find corrections close to $1\%$ in the mildly non-linear regime of $\Lambda$CDM cosmologies for the density power spectrum, while in the case of the density-momentum power spectrum effects can reach up to $1.5\%$ for $k\sim 0.2h/$Mpc.

Inspired by the chromo-natural inflation model of Adshead&Wyman, we reshape its scalar content to relax the tension with current observational bounds. Besides an inflaton, the setup includes a spectator sector in which an axion and SU(2) gauge fields are coupled via a Chern-Simons-type term. The result is a viable theory endowed with an alternative production mechanism for gravitational waves during inflation. The gravitational wave signal sourced by the spectator fields can be much larger than the contribution from standard vacuum fluctuations, it is distinguishable from the latter on the basis of its chirality and, depending on the theory parameters values, also its tilt. This production process breaks the well-known relation between the tensor-to-scalar ratio and the energy scale of inflation. As a result, even if the Hubble rate is itself too small for the vacuum to generate a tensor amplitude detectable by upcoming experiments, this model still supports observable gravitational waves.

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.

Attention has focussed recently on models of inflation that involve a second or more fields with a mass near the inflationary Hubble parameter $H$, as may occur in supersymmetric theories if the supersymmetry-breaking scale is not far from $H$. Quasi-single-field (QsF) inflation is a relatively simple family of phenomenological models that serve as a proxy for theories with additional fields with masses $m\sim H$. Since QsF inflation involves fields in addition to the inflaton, the consistency conditions (ccs) between correlations that arise in single-clock inflation are not necessarily satisfied. As a result, correlation functions in the squeezed limit may be larger than in single-field inflation. Scalar non-Gaussianities mediated by the massive isocurvature field in QsF have been shown to be potentially observable. These are especially interesting since they would convey information about the mass of the isocurvature field. Here we consider non-Gaussian correlators involving tensor modes and their observational signatures. A physical correlation between a (long-wavelength) tensor mode and two scalar modes (tss), for instance, may give rise to local departures from statistical isotropy or, in other words, a non-trivial four-point function. The presence of the tensor mode may moreover be inferred geometrically from the shape dependence of the four-point function. We compute tss and stt (one soft curvature mode and two hard tensors) bispectra in QsF inflation, identifying the conditions necessary for these to "violate" the ccs. We find that while ccs are violated by stt correlations, they are preserved by the tss in the minimal QsF model. Our study of primordial correlators which include gravitons in seeking imprints of additional fields with masses $m\sim H$ during inflation can be seen as complementary to the recent "cosmological collider physics" proposal.

Inflation models make specific predictions for a tensor-scalar-scalar three-point correlation, or bispectrum, between one gravitational-wave (tensor) mode and two density-perturbation (scalar) modes. This tensor-scalar-scalar correlation leads to a local power quadrupole, an apparent departure from statistical isotropy in our Universe, as well as characteristic four-point correlations in the current mass distribution in the Universe. So far, the predictions for these observables have been worked out only for single-clock models in which certain consistency conditions between the tensor-scalar-scalar correlation and tensor and scalar power spectra are satisfied. Here we review the requirements on inflation models for these consistency conditions to be satisfied. We then consider several examples of inflation models, such as non-attractor and solid inflation models, in which these conditions are put to the test. In solid inflation the simplest consistency conditions are already violated whilst in the non-attractor model we find that, contrary to the standard scenario, the tensor-scalar-scalar correlator probes directly relevant model-dependent information. We work out the predictions for observables in these models. For non-attractor inflation we find an apparent local quadrupolar departure from statistical isotropy in large-scale structure but that this power quadrupole decreases very rapidly at smaller scales. The consistency of the CMB quadrupole with statistical isotropy then constrains the distance scale that corresponds to the transition from the non-attractor to attractor phase of inflation to be larger than the currently observable horizon. Solid inflation predicts clustering fossils signatures in the current galaxy distribution that may be large enough to be detectable with forthcoming, and possibly even current, galaxy surveys.

If the amplitude of primordial gravitational waves is measured in the near-future, what could it tell us about bigravity? To address this question, we study massive bigravity theories by focusing on a region in parameter space which is safe from known instabilities. Similarly to investigations on late time constraints, we implicitly assume there is a successful implementation of the Vainshtein mechanism which guarantees that standard cosmological evolution is largely unaffected. We find that viable bigravity models are subject to far less stringent constraints than massive gravity, where there is only one set of (massive) tensor modes. In principle sensitive to the effective graviton mass at the time of recombination, we find that in our setup the primordial tensor spectrum is more responsive to the dynamics of the massless tensor sector rather than its massive counterpart. We further show there are intriguing windows in the parameter space of the theory which could potentially induce distinctive signatures in the B-modes spectrum.

The most conventional mechanism for gravitational waves (gw) production during inflation is the amplification of vacuum metric fluctuations. In this case the gw production can be uniquely related to the inflationary expansion rate $H$. For example, a gw detection close to the present experimental limit (tensor-to-scalar ratio $r \sim 0.1$) would indicate an inflationary expansion rate close to $10^{14} \, {\rm GeV}$. This conclusion, however, would be invalid if the observed gw originated from a different source. We construct and study one of the possible covariant formulations of the mechanism suggested in [43], where a spectator field $\sigma$ with a sound speed $c_{s} \ll 1$ acts as a source for gw during inflation. In our formulation $\sigma$ is described by a so-called $P(X)$ Lagrangian and a non-minimal coupling to gravity. This field interacts only gravitationally with the inflaton, which has a standard action. We compute the amount of scalar and tensor density fluctuations produced by $\sigma$ and find that, in our realization, $r$ is not enhanced with respect to the standard result but it is strongly sensitive to $c_s$, thus breaking the direct $r \leftrightarrow H$ connection.

We present exact FLRW solutions in generalized massive gravity where the mass parameters are naturally promoted to Lorentz-invariant functions of the Stuckelberg fields. This new dependence relaxes the constraint that would otherwise prevent massive gravity from possessing exact FLRW solutions. It does so without the need to introduce additional degrees of freedom. We find self-accelerating cosmological solutions and show that, with a mild restriction on the region of phase space, these cosmological solutions exhibit full stability, i.e. absence of ghosts and gradient instabilities for all the tensor, vector and scalar modes, for all cosmic time. We perform the full decoupling limit analysis, including vector degrees of freedom, which can be used to confirm the existence of an active Vainshtein mechanism about these solutions.

We consider the possible signatures of a recently introduced class of healthy theories beyond Horndeski models on higher-order correlators of the inflationary curvature fluctuation. Despite the apparent large number and complexity of the cubic interactions, we show that the leading-order bispectrum generated by the Generalized Horndeski (also called $G^3$) interactions can be reduced to a linear combination of two well known $k$-inflationary shapes. We conjecture that said behavior is not an accident of the cubic order but a consequence dictated by the requirements on the absence of Ostrogradski instability, the general covariance and the linear dispersion relation in these theories.

We show that every Galileon theory admits a dual formulation as a Galileon theory with new operator coefficients. In n dimensions a free scalar field in Minkowski spacetime is dual to a (n+1)-th order Galileon theory which exhibits the Vainshtein mechanism when coupled to sources and superluminal propagation even on-shell. This demonstrates that superluminal propagation is compatible with an analytic S-matrix and causality. For point sources, the duality interchanges the strongly coupled Vainshtein regime with the weakly coupled asymptotic regime. The duality is made manifest in the context of the decoupling limit of bigravity, but is independent of this.

We show that a large contribution to tensor modes during inflation can be generated by a spectator scalar field with speed of sound lower than unity.

We give a simple derivation of a cosmological bound on the graviton mass for spatially flat FRW solutions in massive gravity with an FRW reference metric and for bigravity theories. This bound comes from the requirement that the kinetic term of the helicity zero mode of the graviton is positive definite. The bound is dependent only on the parameters in the massive gravity potential and the Hubble expansion rate for the two metrics. We derive the decoupling limit of bigravity and FRW massive gravity, and use this to give an independent derivation of the cosmological bound. We recover our previous results that the tension between satisfying the Friedmann equation and the cosmological bound is sufficient to rule out all observationally relevant FRW solutions for massive gravity with an FRW reference metric. In contrast, in bigravity this tension is resolved due to different nature of the Vainshtein mechanism. We find that in bigravity theories there exists an FRW solution with late time self-acceleration for which the kinetic terms for the helicity-2, helicity-1 and helicity-0 are generically nonzero and positive making this a compelling candidate for a model of cosmic acceleration. We confirm that the generalized bound is saturated for the candidate partially massless (bi)gravity theories but the existence of helicity-1/helicity-0 interactions implies the absence of the conjectured partially massless symmetry for both massive gravity and bigravity.

A generalized theory of multi-field galileons has been recently put forward. This model stems from the ongoing effort to embed generic galileon theories within brane constructions. Such an approach has proved very useful in connecting interesting and essential features of these theories with geometric properties of the branes embedding. We investigate the cosmological implications of a very restrictive multi-field galileon theory whose leading interaction is solely quartic in the scalar field and lends itself nicely to an interesting cosmology. The trispectrum of curvature fluctuations has features which are quite distinctive with respect to their P(X,\phi) counterpart. We also show that, despite an absent cubic Lagrangian in the full theory, non-Gaussianities in this model cannot produce the combination a of small bispectrum alongside with a large trispectrum. We further expand on this point to draw a lesson on what having a symmetry in the full background independent theory entails at the level of fluctuations and vice-versa.

We present a detailed study of the trispectrum of the curvature perturbation generated within a stable, well defined and predictive theory which comprises an inflationary phase. In this model the usual shift symmetry is enhanced up to the so-called Galileon symmetry. The appeal of this type of theories rests on being unitary and stable under quantum corrections. Furthermore, in the specific model under consideration here, these properties have been shown to approximately hold in realistic scenarios which account for curved spacetime and the coupling with gravity. In the literature, the analysis of the bispectrum of the curvature perturbation for this theory revealed non-Gaussian features which are shared by a number of inflationary models, including stable ones. It is therefore both timely and useful to investigate further and turn to observables such as the trispectrum. We find that, in a number of specific momenta configurations, the trispectrum shape-functions present strikingly different features as compared to, for example, the entire class of the so-called $P(X,\phi)$ inflationary models.

We prove that the linear instability in a non-degenerate higher derivative theory, the Ostrogradski instability, can only be removed by the addition of constraints if the original theory's phase space is reduced.

We calculate the trispectrum of curvature perturbations for a model of inflation endowed with Galilean symmetry at the level of the fluctuations around an FRW background. Such a model has been shown to posses desirable properties such as unitarity (up to a certain scale) and non-renormalization of the leading operators, all of which point towards the reasonable assumption that a full theory whose fluctuations reproduce the one here might exist as well as be stable and predictive. The cubic curvature fluctuations of this model produce quite distinct signatures at the level of the bispectrum. Our analysis shows how this holds true at higher order in perturbations. We provide a detailed study of the trispectrum shape-functions in different configurations and a comparison with existent literature. Most notably, predictions markedly differ from their P(X,\phi) counterpart in the so called equilateral trispectrum configuration. The zoo of inflationary models characterized by somewhat distinctive predictions for higher order correlators is already quite populated; what makes this model more compelling resides in the above mentioned stability properties.

Chromo-natural inflation is a novel model of inflation which relies on the existence of non-abelian gauge fields interacting with an axion. In its simplest realization, an SU(2) gauge field is assumed to begin inflation in a rotationally invariant VEV. The dynamics of the gauge fields significantly modifies the equations of motion for the axion, providing an additional damping term that supports slow-roll inflation, without the need to fine tune the axion decay constant. We demonstrate that in an appropriate slow-roll limit it is possible to integrate out the massive gauge field fluctuations whilst still maintaining the nontrivial modifications of the gauge field to the axion. In this slow-roll limit, chromo-natural inflation is exactly equivalent to a single scalar field effective theory with a non-minimal kinetic term, i.e. a P(X,\chi) model. This occurs through a precise analogue of the gelaton mechanism, whereby heavy fields can have unsuppressed effects on the light field dynamics without contradicting decoupling. The additional damping effect of the gauge fields can be completely captured by the non-minimal kinetic term of the single scalar field effective theory. We utilize the single scalar field effective theory to infer the power spectrum and non-gaussianities in chromo-natural inflation and confirm that the mass squared of all the gauge field fluctuations is sufficiently large and positive that they completely decouple during inflation. These results confirm that chromo-natural inflation is a viable, stable and compelling model for the generation of inflationary perturbations.

In de Sitter spacetime there exists an absolute minimum for the mass of a spin-2 field set by the Higuchi bound m^2 \geq 2H^2. We generalize this bound to arbitrary spatially flat FRW geometries in the context of the recently proposed ghost-free models of Massive Gravity with an FRW reference metric, by performing a Hamiltonian analysis for cosmological perturbations. We find that the bound generically indicates that spatially flat FRW solutions in FRW massive gravity, which exhibit a Vainshtein mechanism in the background as required by consistency with observations, imply that the helicity zero mode is a ghost. In contradistinction to previous works, the tension between the Higuchi bound and the Vainshtein mechanism is equally strong regardless of the equation of state for matter.

This is a version of the author's Ph.D. thesis. The methods of effective field theory are used to study generic theories of inflation with a single inflaton field and to perform a general analysis of the associated non-Gaussianities. We investigate the amplitudes and shapes of the various three and four-point correlators which are generated by different classes of single-field inflationary models. Besides the well-known results for the so called P(X,\phi) model and the ghost inflationary theories, which we recover, we point out that extrinsic curvature-generated interactions may give rise to large non-Gaussianities with distinctive features in the form of specific shape-functions (e.g. flat, orthogonal etc..) for the correlators.

Employing an effective field theory approach to inflationary perturbations, we analyze in detail the effect of curvature-generated Lagrangian operators on various observables, focusing on their running with scales. At quadratic order, we solve the equation of motion at next-to-leading leading order in a generalized slow-roll approximation for a very general theory of single-field inflation. We derive the resulting power spectrum, its tilt and running. We then focus on the contribution to the primordial non-Gaussianity amplitude f_{NL} sourced by a specific interaction term. We show that the running of f_{NL} can be substantially larger than what dictated by the slow-roll parameters.

The methods of effective field theory are used to study generic theories of inflation with a single inflaton field and to perform a general analysis of the associated non-Gaussianities. We investigate the amplitudes and shapes of the various generic three-point correlators, the bispectra, which may be generated by different classes of single-field inflationary models. Besides the well-known results for the DBI-like models and the ghost inflationary theories, we point out that curvature-related interactions may give rise to large non-Gaussianities in the form of bispectra characterized by a flat shape which, quite interestingly, is independently produced by several interaction terms. In a subsequent work, we will perform a similar general analysis for the non-Gaussianities generated by the generic four-point correlator, the trispectrum.

We perform the analysis of the trispectrum of curvature perturbations generated by the interactions characterizing a general theory of single-field inflation obtained by effective field theory methods. We find that curvature-generated interaction terms, which can in general give an important contribution to the amplitude of the four-point function, show some new distinctive features in the form of their trispectrum shape-function. These interesting interactions are invariant under some recently proposed symmetries of the general theory and, as shown explicitly, do allow for a large value of the trispectrum.