Given the recent advances in gravitational-wave detection technologies, the detection and characterisation of gravitational-wave backgrounds (GWBs) with the Laser Interferometer Space Antenna (LISA) is a real possibility. To assess the abilities of the LISA satellite network to reconstruct anisotropies of different angular scales and in different directions on the sky, we develop a map-maker based on an optimal quadratic estimator. The resulting maps are maximum likelihood representations of the GWB intensity on the sky integrated over a broad range of frequencies. We test the algorithm by reconstructing known input maps with different input distributions and over different frequency ranges. We find that, in an optimal scenario of well understood noise and high frequency, high SNR signals, the maximum scales LISA may probe are $\ell_{\rm max} \lesssim 15$. The map-maker also allows to test the directional dependence of LISA noise, providing insight on the directional sky sensitivity we may expect.

The detection of gravitational waves (GWs) propagating through cosmic structures can provide invaluable information on the geometry and content of our Universe, as well as on the fundamental theory of gravity. In order to test possible departures from General Relativity, it is essential to analyse, in a modified gravity setting, how GWs propagate through a perturbed cosmological space-time. Working within the framework of geometrical optics, we develop tools to address this topic for a broad class of scalar-tensor theories, including scenarios with non-minimal, derivative couplings between scalar and tensor modes. We determine the corresponding evolution equations for the GW amplitude and polarization tensor. The former satisfies a generalised evolution equation that includes possible effects due to a variation of the effective Planck scale; the latter can fail to be parallely transported along GW geodesics unless certain conditions are satisfied. We apply our general formulas to specific scalar-tensor theories with unit tensor speed, and then focus on GW propagation on a perturbed space-time. We determine corrections to standard formulas for the GW luminosity distance and for the evolution of the polarization tensor, which depend both on modified gravity and on the effects of cosmological perturbations. Our results can constitute a starting point to disentangle among degeneracies from different sectors that can influence GW propagation through cosmological space-times.

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.

The possibility that primordial black holes constitute a fraction of dark matter motivates a detailed study of possible mechanisms for their production. Black holes can form by the collapse of primordial curvature fluctuations, if the amplitude of their small scale spectrum gets amplified by several orders of magnitude with respect to CMB scales. Such enhancement can for example occur in single-field inflation that exhibit a transient non-attractor phase: in this work, we make a detailed investigation of the shape of the curvature spectrum in this scenario. We make use of an analytical approach based on a gradient expansion of curvature perturbations, which allows us to follow the changes in slope of the spectrum during its way from large to small scales. After encountering a dip in its amplitude, the spectrum can acquire steep slopes with a spectral index up to $n_s-1\,=\,8$, to then relax to a more gentle growth with $n_s-1\,\lesssim\,3$ towards its peak, in agreement with the results found in previous literature. For scales following the peak associated with the presence of the non-attractor phase, the spectrum amplitude then mildly decays, during a transitional stage from non-attractor back to attractor evolution. Our analysis indicates that this gradient approach offers a transparent understanding of the contributions controlling the slope of the curvature spectrum. As an application of our findings, we characterise the slope in frequency of a stochastic gravitational wave background generated at second order from curvature fluctuations, using the more accurate information we gained on the shape of curvature power spectrum.

An observable stochastic background of gravitational waves is generated whenever primordial black holes are created in the early universe thanks to a small-scale enhancement of the curvature perturbation. We calculate the anisotropies and non-Gaussianity of such stochastic gravitational waves background which receive two contributions, the first at formation time and the second due to propagation effects. The former contribution can be generated if the distribution of the curvature perturbation is characterized by a local and scale-invariant shape of non-Gaussianity. Under such an assumption, we conclude that a sizeable magnitude of anisotropy and non-Gaussianity in the gravitational waves would suggest that primordial black holes may not comply the totality of the dark matter.

The stochastic gravitational wave background (SGWB) contains a wealth of information on astrophysical and cosmological processes. A major challenge of upcoming years will be to extract the information contained in this background and to disentangle the contributions of different sources. In this paper we provide the formalism to extract, from the correlation of three signals in the Laser Interferometer Space Antenna (LISA), information about the tensor three-point function, which characterizes the non-Gaussian properties of the SGWB. This observable can be crucial to discriminate whether a SGWB has a primordial or astrophysical origin. Compared to the two-point function, the SGWB three-point function has a richer dependence on the gravitational wave momenta and chiralities. It can be used therefore as a powerful discriminator between different models. For the first time we provide the response functions of LISA to a general SGWB three-point function. As examples, we study in full detail the cases of an equilateral and squeezed SGWB bispectra, and provide the explicit form of the response functions, ready to be convoluted with any theoretical prediction of the bispectrum to obtain the observable signal. We further derive the optimal estimator to compute the signal-to-noise ratio. Our formalism covers general shapes of non-Gaussianity, and can be extended straightaway to other detector geometries. Finally, we provide a short overview of models of the early universe that can give rise to a non-Gaussian SGWB.

A future detection of the Stochastic Gravitational Wave Background (SGWB) with GW experiments is expected to open a new window on early universe cosmology and on the astrophysics of compact objects. In this paper we study SGWB anisotropies, that can offer new tools to discriminate between different sources of GWs. In particular, the cosmological SGWB inherits its anisotropies both (i) at its production and (ii) during its propagation through our perturbed universe. Concerning (i), we show that it typically leads to anisotropies with order one dependence on frequency. We then compute the effect of (ii) through a Boltzmann approach, including contributions of both large-scale scalar and tensor linearized perturbations. We also compute for the first time the three-point function of the SGWB energy density, which can allow one to extract information on GW non-Gaussianity with interferometers. Finally, we include non-linear effects associated with long wavelength scalar fluctuations, and compute the squeezed limit of the 3-point function for the SGWB density contrast. Such limit satisfies a consistency relation, conceptually similar to what found in the literature for the case of CMB perturbations.

Parity violating interactions in the early Universe can source a stochastic gravitational wave background (SGWB) with a net circular polarization. In this paper, we study possible ways to search for circular polarization of the SGWB with interferometers. Planar detectors are unable to measure the net circular polarization of an isotropic SGWB. We discuss the possibility of using the dipolar anisotropy kinematically induced by the motion of the solar system with respect to the cosmic reference frame to measure the net circular polarization of the SGWB with planar detectors. We apply this approach to LISA, re-assessing previous analyses by means of a more detailed computation and using the most recent instrument specifications, and to the Einstein Telescope (ET), estimating for the first time its sensitivity to circular polarization. We find that both LISA and ET, despite operating at different frequencies, could detect net circular polarization with a signal-to-noise ratio of order one in a SGWB with amplitude $h^2 \Omega_\text{GW} \simeq 10^{-11}$. We also investigate the case of a network of ground based detectors. We present fully analytical, covariant formulas for the detector overlap functions in the presence of circular polarization. Our formulas do not rely on particular choices of reference frame, and can be applied to interferometers with arbitrary angles among their arms.

Dimensional flow, the scale dependence of the dimensionality of spacetime, is a feature shared by many theories of quantum gravity (QG). We present the first study of the consequences of QG dimensional flow for the luminosity distance scaling of gravitational waves in the frequency ranges of LIGO and LISA. We find generic modifications with respect to the standard general-relativistic scaling, largely independent of specific QG proposals. We constrain these effects using two examples of multimessenger standard sirens, the binary neutron-star merger GW170817 and a simulated supermassive black-hole merger event detectable with LISA. We apply these constraints to various QG candidates, finding that the quantum geometries of group field theory, spin foams and loop quantum gravity can give rise to observable signals in the gravitational-wave spin-2 sector. Our results complement and improve GW propagation-speed bounds on modified dispersion relations. Under more model-dependent assumptions, we also show that bounds on quantum geometry can be strengthened by solar-system tests.

We investigate possible signatures of quantum gravity which could be tested with current and future gravitational-wave (GW) observations. In particular, we analyze how quantum gravity can influence the GW luminosity distance, the time dependence of the effective Planck mass and the instrumental strain noise of interferometers. Using both model-dependent and model-independent formulae, we show that these quantities can encode a non-perturbative effect typical of all quantum-gravity theories, namely the way the dimension of spacetime changes with the probed scale. Effects associated with such dimensional flow might be tested with GW observations and constrained significantly in theories with a microscopically discrete spacetime geometry, more strongly than from propagation-speed constraints. Making use of public LIGO data as well as of a simulated higher-redshift LISA source, we impose the first, respectively, actual and mock constraints on quantum-gravity parameters affecting the GW luminosity distance and discuss specific theoretical examples. If also the Newtonian potential is modified but light geodesics are not, then solar-system bounds may be stronger than GW ones. Yet, for some theories GW astronomy can give unique information not available from solar-system tests.

Since the very beginning of astronomy the location of objects on the sky has been a fundamental observational quantity that has been taken for granted. While precise two dimensional positional information is easy to obtain for observations in the electromagnetic spectrum, the positional accuracy of current and near future gravitational wave detectors is limited to between tens and hundreds of square degrees, which makes it extremely challenging to identify the host galaxies of gravitational wave events or to confidently detect any electromagnetic counterparts. Gravitational wave observations provide information on source properties and distances that is complementary to the information in any associated electromagnetic emission and that is very hard to obtain in any other way. Observing systems with multiple messengers thus has scientific potential much greater than the sum of its parts. A gravitational wave detector with higher angular resolution would significantly increase the prospects for finding the hosts of gravitational wave sources and triggering a multi-messenger follow-up campaign. An observatory with arcminute precision or better could be realised within the Voyage 2050 programme by creating a large baseline interferometer array in space and would have transformative scientific potential. Precise positional information of standard sirens would enable precision measurements of cosmological parameters and offer new insights on structure formation; a high angular resolution gravitational wave observatory would allow the detection of a stochastic background and resolution of the anisotropies within it; it would also allow the study of accretion processes around black holes; and it would have tremendous potential for tests of modified gravity and the discovery of physics beyond the Standard Model.

The Stochastic Gravitational Wave Background (SGWB) is expected to be a key observable for Gravitational Wave (GW) interferometry. Its detection will open a new window on early universe cosmology and on the astrophysics of compact objects. Using a Boltzmann approach, we study the angular anisotropies of the GW energy density, which is an important tool to disentangle the different cosmological and astrophysical contributions to the SGWB. Anisotropies in the cosmological background are imprinted both at its production, and by GW propagation through the large-scale scalar and tensor perturbations of the universe. The first contribution is not present in the Cosmic Microwave Background (CMB) radiation (as the universe is not transparent to photons before recombination), causing an order one dependence of the anisotropies on frequency. Moreover, we provide a new method to characterize the cosmological SGWB through its possible deviation from a Gaussian statistics. In particular, the SGWB will become a new probe of the primordial non-Gaussianity of the large-scale cosmological perturbations.

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]

We present a set of tools to assess the capabilities of LISA to detect and reconstruct the spectral shape and amplitude of a stochastic gravitational wave background (SGWB). We first provide the LISA power-law sensitivity curve and binned power-law sensitivity curves, based on the latest updates on the LISA design. These curves are useful to make a qualitative assessment of the detection and reconstruction prospects of a SGWB. For a quantitative reconstruction of a SGWB with arbitrary power spectrum shape, we propose a novel data analysis technique: by means of an automatized adaptive procedure, we conveniently split the LISA sensitivity band into frequency bins, and fit the data inside each bin with a power law signal plus a model of the instrumental noise. We apply the procedure to SGWB signals with a variety of representative frequency profiles, and prove that LISA can reconstruct their spectral shape. Our procedure, implemented in the code SGWBinner, is suitable for homogeneous and isotropic SGWBs detectable at LISA, and it is also expected to work for other gravitational wave observatories.

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.

The LIGO discoveries have rekindled suggestions that primordial black holes (BHs) may constitute part to all of the dark matter (DM) in the Universe. Such suggestions came from 1) the observed merger rate of the BHs, 2) their unusual masses, 3) their low/zero spins, and 4) also from the independently uncovered cosmic infrared background (CIB) fluctuations signal of high amplitude and coherence with unresolved cosmic X-ray background (CXB). Here we summarize the prospects to resolve this important issue with electromagnetic observations using the instruments and tools expected in the 2020's. These prospects appear promising to make significant, and potentially critical, advances. We demonstrate that in the next decade, new space- and ground-borne electromagnetic instruments, combined with concurrent theoretical efforts, should shed critical light on the long-considered link between primordial BHs and DM. Specifically the new data and methodologies under this program will involve: I) Probing with high precision the spatial spectrum of source-subtracted CIB with Euclid and WFIRST, and its coherence with unresolved cosmic X-ray background using eROSITA and Athena, II) Advanced searches for microlensing of Galactic stars by the intervening Galactic Halo BHs with OGLE, Gaia, LSST and WFIRST, III) Supernovae (SNe) lensing in the upcoming surveys with WFIRST, LSST and also potentially with Euclid and JWST, IV) Advanced theoretical work to understand the details of PBH accretion and evolution and their influence on cosmic microwave background (CMB) anisotropies in light of the next generation CMB experiments, V) Better new samples and theoretical understanding involving stability and properties of ultra faint dwarf galaxies, pulsar timing, and cosmological quasar lensing.

We discuss a new mechanism to enhance the spectrum of primordial tensor fluctuations in single field inflationary scenarios. The enhancement relies on a transitory non-attractor inflationary phase, which amplifies the would-be decaying tensor mode, and gives rise to a growth of tensor fluctuations at superhorizon scales. We show that the enhancement produced during this phase can be neatly treated via a tensor duality between an attractor and non-attractor phase, which we introduce. We illustrate the mechanism and duality in a kinetically driven scenario of inflation, with non-minimal couplings between the scalar and the metric.

We investigate primordial tensor non-Gaussianity in single field inflation, during a phase of non-attractor evolution when the spectrum of primordial tensor modes can be enhanced to a level detectable at interferometer scales. Making use of a tensor duality we introduced in arXiv:1808.10475, we analytically compute the full bispectrum of primordial tensor fluctuations during the non-attractor era. During this epoch the shape of the tensor bispectrum is enhanced in the squeezed limit, its amplitude can be amplified with respect to slow-roll models, and tensor non-Gaussianity can exhibit a scale dependence distinctive of our set-up. We prove that our results do not depend on the frame used for the calculations. Squeezed tensor non-Gaussianity induces a characteristic quadrupolar anisotropy on the power spectrum of the stochastic background of primordial tensor perturbations. As a step to make contact with gravitational wave experiments, we discuss the response function of a ground based Michelson interferometer to a gravitational wave background with such a feature.

We propose a simple, well-motivated and robust alternative to a metastable de Sitter vacuum in string theory, consistent with current observations of dark energy and naturally satisfying conjectured swampland constraints. Inflation ends in a supersymmetric Minkowski minimum, with a flat direction that is protected by non-renormalisation theorems. At some scale non-perturbative effects kick in, inducing a runaway scalar potential. The tail of this runaway potential cannot sustain a late-time dominating, slow-roll quintessence. However, the potential always contains a dS maximum. If the modulus starts close to the hilltop, it remains frozen there by Hubble friction for much of the cosmological history, at first sourcing a tiny classical vacuum energy and then constituting a rolling quintessence, with observable consequences. So long as the modulus is localised away from the Standard Model in the extra dimensions, there are no observable fifth forces nor is there time-variation of fundamental constants, and the modulus mass is protected from radiative corrections. We revisit concrete string models based on heterotic orbifold compactifications, and show that their de Sitter extrema satisfy the recently refined dS Swampland Conjecture.

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.

We forecast the prospective of detection for a stochastic gravitational wave background sourced by cosmological first-order phase transitions. We focus on first-order phase transitions with negligible plasma effects, and consider the experimental infrastructures built by the end of the LISA mission. We make manifest the synergy among LISA, pulsar time array experiments, and ground-based interferometers. For phase transitions above the TeV scale or below the electroweak scale, LISA can detect the corresponding gravitational wave signal together with Einstein Telescope, SKA or even aLIGO-aVIRGO-KAGRA. For phase transitions at the electroweak scale, instead, LISA can be the only experiment observing the gravitational wave signal. In case of detection, by using a parameter reconstruction method that we anticipate in this work, we show that LISA on its own has the potential to determine when the phase transition occurs and, consequently, the energy scale above which the standard model of particle physics needs to be modified. The result may likely guide the collider community in the post-LHC era.

We consider mechanisms for producing a significant population of primordial black holes (PBHs) within string inspired single field models of inflation. The production of PBHs requires a large amplification in the power spectrum of curvature perturbations between scales associated with CMB and PBH formation. In principle, this can be achieved by temporarily breaking the slow-roll conditions during inflation. In this work, we identify two string setups that can realise this process. In string axion models of inflation, subleading non-perturbative effects can superimpose steep cliffs and gentle plateaus onto the leading axion potential. The cliffs can momentarily violate the slow-roll conditions, and the plateaus can lead to phases of ultra slow-roll inflation. We thus achieve a string motivated model which both matches the Planck observations at CMB scales and produces a population of light PBHs, which can account for an order one fraction of dark matter. In DBI models of inflation, a sharp increase in the speed of sound sourced by a steep downward step in the warp factor can drive the amplification. In this scenario, discovery of PBHs could indicate non-trivial dynamics in the bulk, such as flux-antibrane annihilation at the tip of a warped throat.

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 develop a scenario of inflation with spontaneously broken time and space diffeomorphisms, with distinctive features for the primordial tensor modes. Inflationary tensor fluctuations are non adiabatic, and can acquire a mass during the inflationary epoch. They can evade the Higuchi bound around de Sitter space, thanks to interactions with the fields driving expansion. Correspondingly, the primordial stochastic gravitational wave background (SGWB) is characterised by a tuneable scale dependence, and can be detectable at interferometer scales. In this set-up, tensor non-Gaussianity can be parametrically enhanced in the squeezed limit. This induces a coupling between long and short tensor modes, leading to a specific quadrupolar anisotropy in the primordial SGWB spectrum, which can be used to build estimators for tensor non-Gaussianity. We analyse how our inflationary system can be tested with interferometers, also discussing how an interferometer can be sensitive to a primordial anisotropic SGWB.

Future observations of cosmic microwave background (CMB) polarisation have the potential to answer some of the most fundamental questions of modern physics and cosmology. In this paper, we list the requirements for a future CMB polarisation survey addressing these scientific objectives, and discuss the design drivers of the CORE space mission proposed to ESA in answer to the "M5" call for a medium-sized mission. The rationale and options, and the methodologies used to assess the mission's performance, are of interest to other future CMB mission design studies. CORE is designed as a near-ultimate CMB polarisation mission which, for optimal complementarity with ground-based observations, will perform the observations that are known to be essential to CMB polarisation scienceand cannot be obtained by any other means than a dedicated space mission.

Supersolid inflation is a class of inflationary theories that simultaneously breaks time and space reparameterization invariance during inflation, with distinctive features for the dynamics of cosmological fluctuations. We investigate concrete realizations of such a scenario, including non-minimal couplings between gravity and the fields driving inflation. We focus in particular on the dynamics of primordial gravitational waves and discuss how their properties depend on the pattern of symmetry breaking that we consider. Tensor modes can have a blue spectrum, and for the first time we build models in which the squeezed limit of primordial tensor bispectra can be parametrically enhanced with respect to standard single-field scenarios. At leading order in a perturbative expansion, the tensor-to-scalar ratio depends only on the parameter controlling the breaking of space-reparameterization. It is independent from the quantities controlling the breaking of time-reparameterization, and this represents a difference with respect to standard single-field inflationary models.

We forecast the scientific capabilities to improve our understanding of cosmic inflation of CORE, a proposed CMB space satellite submitted in response to the ESA fifth call for a medium-size mission opportunity. The CORE satellite will map the CMB anisotropies in temperature and polarization in 19 frequency channels spanning the range 60-600 GHz. CORE will have an aggregate noise sensitivity of $1.7 \mu$K$\cdot \,$arcmin and an angular resolution of 5' at 200 GHz. We explore the impact of telescope size and noise sensitivity on the inflation science return by making forecasts for several instrumental configurations. This study assumes that the lower and higher frequency channels suffice to remove foreground contaminations and complements other related studies of component separation and systematic effects, which will be reported in other papers of the series "Exploring Cosmic Origins with CORE." We forecast the capability to determine key inflationary parameters, to lower the detection limit for the tensor-to-scalar ratio down to the $10^{-3}$ level, to chart the landscape of single field slow-roll inflationary models, to constrain the epoch of reheating, thus connecting inflation to the standard radiation-matter dominated Big Bang era, to reconstruct the primordial power spectrum, to constrain the contribution from isocurvature perturbations to the $10^{-3}$ level, to improve constraints on the cosmic string tension to a level below the presumptive GUT scale, and to improve the current measurements of primordial non-Gaussianities down to the $f_{NL}^{\rm local} < 1$ level. For all the models explored, CORE alone will improve significantly on the present constraints on the physics of inflation. Its capabilities will be further enhanced by combining with complementary future cosmological observations.

The direct detection of gravitational waves opens new perspectives for measuring properties of gravitationally bound compact objects. It is then important to investigate black holes and neutron stars in alternative theories of gravity, since they can have features that make them observationally distinguishable from their General Relativity (GR) counterparts. In this work, we examine a special case of vector Galileons, a vector-tensor theory of gravity with interesting cosmological properties, which consists of a one parameter modification of the Einstein-Maxwell action. Within this theory, we study configurations describing asymptotically flat, spherically symmetric black holes and neutron stars. The set of black hole solutions in this theory is surprisingly rich, generalising results found in GR or in related scalar-tensor theories. We investigate the properties and conserved charges of black holes, using both analytical and numerical techniques, highlighting configurations that are more compact than in GR. We then study properties of neutron stars, showing how the vector profile can influence the star internal structure. Depending on properties of matter and fields inside the star, neutron stars can be more massive than in GR, and they can be even more compact than Schwarzschild black holes, making these objects observationally interesting. We also comment on possible extensions of our configurations to magnetically charged or rotating configurations.

We investigate the potential for the LISA space-based interferometer to detect the stochastic gravitational wave background produced from different mechanisms during inflation. Focusing on well-motivated scenarios, we study the resulting contributions from particle production during inflation, inflationary spectator fields with varying speed of sound, effective field theories of inflation with specific patterns of symmetry breaking and models leading to the formation of primordial black holes. The projected sensitivities of LISA are used in a model-independent way for various detector designs and configurations. We demonstrate that LISA is able to probe these well-motivated inflationary scenarios beyond the irreducible vacuum tensor modes expected from any inflationary background.

We present all scalar-tensor Lagrangians that are cubic in second derivatives of a scalar field, and that are degenerate, hence avoiding Ostrogradsky instabilities. Thanks to the existence of constraints, they propagate no more than three degrees of freedom, despite having higher order equations of motion. We also determine the viable combinations of previously identified quadratic degenerate Lagrangians and the newly established cubic ones. Finally, we study whether the new theories are connected to known scalar-tensor theories such as Horndeski and beyond Horndeski, through conformal or disformal transformations.

Measurements of the non-Gaussianity of the primordial density field have the power to considerably improve our understanding of the physics of inflation. Indeed, if we can increase the precision of current measurements by an order of magnitude, a null-detection would rule out many classes of scenarios for generating primordial fluctuations. Large-scale galaxy redshift surveys represent experiments that hold the promise to realise this goal. Thus, we model the galaxy bispectrum and forecast the accuracy with which it will probe the parameter $f_{\rm NL}$, which represents the degree of primordial local-type non Gaussianity. Specifically, we address the problem of modelling redshift space distortions (RSD) in the tree-level galaxy bispectrum including $f_{\rm NL}$. We find novel contributions associated with RSD, with the characteristic large scale amplification induced by local-type non-Gaussianity. These RSD effects must be properly accounted for in order to obtain un-biased measurements of $f_{\rm NL}$ from the galaxy bispectrum. We propose an analytic template for the monopole which can be used to fit against data on large scales, extending models used in the recent measurements. Finally, we perform idealised forecasts on $\sigma_{f_{\rm NL}}$ -- the accuracy of the determination of local non-linear parameter $f_{\rm NL}$ -- from measurements of the galaxy bispectrum. Our findings suggest that current surveys can in principle provide $f_{\rm NL}$ constraints competitive with Planck, and future surveys could improve them further.

We study new consistent scalar-tensor theories of gravity recently introduced by Langlois and Noui with potentially interesting cosmological applications. We derive the conditions for the existence of a primary constraint that prevents the propagation of an additional dangerous mode associated with higher order equations of motion. We then classify the most general, consistent scalar-tensor theories that are at most quadratic in the second derivatives of the scalar field. In addition, we investigate the possible connection between these theories and (beyond) Horndeski through conformal and disformal transformations. Finally, we point out that these theories can be associated with new operators in the effective field theory of dark energy, which might open up new possibilities to test dark energy models in future surveys.

Determining the most general, consistent scalar tensor theory of gravity is important for building models of inflation and dark energy. In this work we investigate the number of degrees of freedom present in the theory of beyond Horndeski. We discuss how to construct the theory from the extrinsic curvature of the constant scalar field hypersurface, and find a simple expression for the action which guarantees the existence of the primary constraint necessary to avoid the Ostrogradsky instability. Our analysis is completely gauge-invariant. However we confirm that, mixing together beyond Horndeski with a different order of Horndeski, obstructs the construction of this primary constraint. Instead, when the mixing is between actions of the same order, the theory can be mapped to Horndeski through a generalised disformal transformation. This mapping however is impossible with beyond Horndeski alone, since we find that the theory is invariant under such a transformation. The picture that emerges is that beyond Horndeski is a healthy but isolated theory: combined with Horndeski, it either becomes Horndeski, or likely propagates a ghost.

An attractive candidate for the inflaton is an axion slowly rolling down a flat potential protected by a perturbative shift symmetry. Realisations of this idea within large field, natural and monomial inflation have been disfavoured by observations and are difficult to embed in string theory. We show that subleading, but significant non-perturbative corrections can superimpose sharp cliffs and gentle plateaus into the potential, whose overall effect is to enhance the number of e-folds of inflation. Sufficient e-folds are therefore achieved for smaller field ranges compared to the potential without such corrections. Thus, both single-field natural and monomial inflation in UV complete theories like string theory, can be restored into the favour of current observations, with distinctive signatures. Tensor modes result un-observably small, but there is a large negative running of the spectral index. Remarkably, natural inflation can be achieved with a single field whose axion decay constant is sub-Planckian.

Though simple inflationary models describe the CMB well, their corrections are often plagued by infrared effects that obstruct a reliable calculation of late-time behaviour. We adapt to cosmology tools designed to address similar issues in other physical systems with the goal of making reliable late-time inflationary predictions. The main such tool is Open EFTs which reduce in the inflationary case to Stochastic Inflation plus calculable corrections. We apply this to a simple inflationary model that is complicated enough to have dangerous IR behaviour yet simple enough to allow the inference of late-time behaviour. We find corrections to standard Stochastic Inflationary predictions for the noise and drift, and we find these corrections ensure the IR finiteness of both these quantities. The late-time probability distribution, ${\cal P}(\phi)$, for super-Hubble field fluctuations are obtained as functions of the noise and drift and so these too are IR finite. We compare our results to other methods (such as large-$N$ models) and find they agree when these models are reliable. In all cases we can explore in detail we find IR secular effects describe the slow accumulation of small perturbations to give a big effect: a significant distortion of the late-time probability distribution for the field. But the energy density associated with this is only of order $H^4$ at late times and so does {\em not} generate a dramatic gravitational back-reaction.

The effective field theory of inflation is a powerful tool for obtaining model independent predictions common to large classes of inflationary models. It requires only information about the symmetries broken during the inflationary era, and on the number and nature of fields that drive inflation. In this paper, we consider the case for scenarios that simultaneously break time reparameterization and spatial diffeomorphisms during inflation. We examine how to analyse such systems using an effective field theory approach, and we discuss several observational consequences for the statistics of scalar and tensor modes. For example, examining the three point functions, we show that this symmetry breaking pattern can lead to an enhanced amplitude for the squeezed bispectra, and to a distinctive angular dependence between their three wavevectors. We also discuss how our results indicate prospects for constraining the level of spatial diffeomorphism breaking during inflation.

Vector theories with non-linear derivative self-interactions that break gauge symmetries have been shown to have interesting cosmological applications. In this paper we introduce a way to spontaneously break the gauge symmetry and construct these theories via a Higgs mechanism. In addition to the purely gauge field interactions, our method generates new ghost-free scalar-vector interactions between the Higgs field and the gauge boson. We show how these additional terms are found to reduce, in a suitable decoupling limit, to scalar bi-Galileon interactions between the Higgs field and Goldstone bosons. Our formalism is first developed in the context of abelian symmetry, which allows us to connect with earlier work on the extension of the Proca action. We then show how this formalism is straightforwardly generalised to generate theories with non-abelian symmetry.

Primordial non-Gaussianity can lead to a scale-dependent bias in the density of collapsed halos relative to the underlying matter density. The galaxy power spectrum already provides constraints on local-type primordial non-Gaussianity complementary those from the cosmic microwave background (CMB), while the bispectrum contains additional shape information and has the potential to outperform CMB constraints in future. We develop the bias model for the halo density contrast in the presence of local-type primordial non-Gaussianity, deriving a bivariate expansion up to second order in terms of the local linear matter density contrast and the local gravitational potential in Lagrangian coordinates. Non-linear evolution of the matter density introduces a non-local tidal term in the halo model. Furthermore, the presence of local-type non-Gaussianity in the Lagrangian frame leads to a novel non-local convective term in the Eulerian frame, that is proportional to the displacement field when going beyond the spherical collapse approximation. We use an extended Press-Schechter approach to evaluate the halo mass function and thus the halo bispectrum. We show that including these non-local terms in the halo bispectra can lead to corrections of up to $25\%$ for some configurations, on large scales or at high redshift.

We study the phenomenon of discrete symmetry breaking during the inflationary epoch, using a model-independent approach based on the effective field theory of inflation. We work in a context where both time reparameterization symmetry and spatial diffeomorphism invariance can be broken during inflation. We determine the leading derivative operators in the quadratic action for fluctuations that break parity and time-reversal. Within suitable approximations, we study their consequences for the dynamics of linearized fluctuations. Both in the scalar and tensor sectors, we show that such operators can lead to new direction-dependent phases for the modes involved. They do not affect the power spectra, but can have consequences for higher correlation functions. Moreover, a small quadrupole contribution to the sound speed can be generated.

The basic building block for Lorentz invariant and ghost free massive gravity is the square root of the combination $g^{-1}\eta\,$, where $g^{-1}$ is the inverse of the physical metric and $\eta$ is a reference metric. Since the square root of a matrix is not uniquely defined, it is possible to have physically inequivalent potentials corresponding to different branches. We show that around Minkowski background the only perturbatively well defined branch is the potential proposed by de Rham, Gabadadze and Tolley. On the other hand, if Lorentz symmetry is broken spontaneously, other potentials exist with a standard perturbative expansion. We show this explicitly building new Lorentz invariant, ghost-free massive gravity potentials for theories that in the background preserve rotational invariance, but break Lorentz boosts.

We study cosmology in the bigravity formulation of the dRGT model where matter couples to both metrics. At linear order in perturbation theory two mass scales emerge: an hard one from the dRGT potential, and an environmental dependent one from the coupling of bigravity with matter. At early time, the dynamics is dictated by the second mass scale which is of order of the Hubble scale. The set of gauge invariant perturbations that couples to matter follow closely the same behaviour as in GR. The remaining perturbations show no issue in the scalar sector, while problems arise in the tensor and vector sectors. During radiation domination, a tensor mode grows power-like at super-horizon scales. More dangerously, the only propagating vector mode features an exponential instability on sub-horizon scales. We discuss the consequences of such instabilities and speculate on possible ways to deal with them.

Cosmological magnetic fields pervade the entire universe, from small to large scales. Since they apparently extend into the intergalactic medium, it is tantalizing to believe that they have a primordial origin, possibly being produced during inflation. However, finding consistent scenarios for inflationary magnetogenesis is a challenging theoretical problem. The requirements to avoid an excessive production of electromagnetic energy, and to avoid entering a strong coupling regime characterized by large values for the electromagnetic coupling constant, typically allow one to generate only a tiny amplitude of magnetic field during inflation. We propose a scenario for building gauge-invariant models of inflationary magnetogenesis potentially free from these issues. The idea is to derivatively couple a dynamical scalar, not necessarily the inflaton, to fermionic and electromagnetic fields during the inflationary era. Such couplings give additional freedom to control the time-dependence of the electromagnetic coupling constant during inflation. This fact allows us to find conditions to avoid the strong coupling problems that affect many of the existing models of magnetogenesis. We do not need to rely on a particular inflationary set-up for developing our scenario, that might be applied to different realizations of inflation. On the other hand, specific requirements have to be imposed on the dynamics of the scalar derivatively coupled to fermions and electromagnetism, that we are able to satisfy in an explicit realization of our proposal.

Using an effective field theory approach to inflation, we examine novel properties of the spectrum of inflationary tensor fluctuations, that arise when breaking some of the symmetries or requirements usually imposed on the dynamics of perturbations. During single-clock inflation, time-reparameterization invariance is broken by a time-dependent cosmological background. In order to explore more general scenarios, we consider the possibility that spatial diffeomorphism invariance is also broken by effective mass terms or by derivative operators for the metric fluctuations in the Lagrangian. We investigate the cosmological consequences of the breaking of spatial diffeomorphisms, focussing on operators that affect the power spectrum of fluctuations. We identify the operators for tensor fluctuations that can provide a blue spectrum without violating the null energy condition, and operators for scalar fluctuations that lead to non-conservation of the comoving curvature perturbation on superhorizon scales even in single-clock inflation. In the last part of our work, we also examine the consequences of operators containing more than two spatial derivatives, discussing how they affect the sound speed of tensor fluctuations, and showing that they can mimic some of the interesting effects of symmetry breaking operators, even in scenarios that preserve spatial diffeomorphism invariance.

We investigate some cosmological consequences of a vector-tensor theory where an Abelian symmetry in the vector sector is slightly broken by a mass term and by ghost-free derivative self-interactions. When studying cosmological expansion in the presence of large bare cosmological constant $\Lambda_{cc}$, we find that the theory admits branches of de Sitter solutions in which the scale of the Hubble parameter is inversely proportional to a power of $\Lambda_{cc}$. Hence, a large value of $\Lambda_{cc}$ leads to a small size for the Hubble scale. In an appropriate limit, in which the symmetry breaking parameters are small, the theory recovers the Abelian symmetry plus an additional Galileon symmetry acting on the longitudinal vector polarization. The approximate Galileon symmetry can make the structure of this theory stable at the energy scales we are interested in. We also analyze the dynamics of linearized cosmological fluctuations around the de Sitter solutions, showing that no manifest instabilities arise, and that the transverse vector polarizations become massless around these configurations.

In the context of a cubic Galileon model in which the Vainshtein mechanism suppresses the scalar field interactions with matter, we study low-density stars with slow rotation and static relativistic stars. We develop an expansion scheme to find approximated solutions inside the Vainshtein radius, and show that deviations from General Relativity (GR), while considering rotation, are also suppressed by the Vainshtein mechanism. In a quadratic coupling model, in which the scalarisation effect can significantly enhance deviations from GR in normal scalar tensor gravity, the Galileon term successfully suppresses the large deviations away from GR. Moreover, using a realistic equation of state, we construct solutions for a relativistic star, and show that deviations from GR are more suppressed for higher density objects. However, we found that the scalar field solution ceases to exist above a critical density, which roughly corresponds to the maximum mass of a neutron star. This indicates that, for a compact object described by a polytropic equation of state, the configuration that would collapse into a black hole cannot support a non-trivial scalar field.

We identify the effective theory describing inflationary super-Hubble scales and show it to be a special case of effective field theories appropriate to open systems. Open systems allow information to be exchanged between the degrees of freedom of interest and those that are integrated out, such as for particles moving through a fluid. Strictly speaking they cannot in general be described by an effective lagrangian; rather the appropriate `low-energy' limit is instead a Lindblad equation describing the evolution of the density matrix of the slow degrees of freedom. We derive the equation relevant to super-Hubble modes of quantum fields in near-de Sitter spacetimes and derive two implications. We show the evolution of the diagonal density-matrix elements quickly approaches the Fokker-Planck equation of Starobinsky's stochastic inflationary picture. This provides an alternative first-principles derivation of this picture's stochastic noise and drift, as well as its leading corrections. (An application computes the noise for systems with a sub-luminal sound speed.) We argue that the presence of interactions drives the off-diagonal density-matrix elements to zero in the field basis. This shows why the field basis is the `pointer basis' for the decoherence of primordial quantum fluctuations while they are outside the horizon, thus allowing them to re-enter as classical fluctuations, as assumed when analyzing CMB data. The decoherence process is efficient, occurring after several Hubble times even for interactions as weak as gravitational-strength. Crucially, the details of the interactions largely control only the decoherence time and not the nature of the final late-time stochastic state, much as interactions can control the equilibration time for thermal systems but are largely irrelevant to the properties of the resulting equilibrium state.

We study static, spherically symmetric solutions in a recently proposed ghost-free model of non-linear massive gravity. We focus on a branch of solutions where the helicity-0 mode can be strongly coupled within certain radial regions, giving rise to the Vainshtein effect. We truncate the analysis to scales below the gravitational Compton wavelength, and consider the weak field limit for the gravitational potentials, while keeping all non-linearities of the helicity-0 mode. We determine analytically the number and properties of local solutions which exist asymptotically on large scales, and of local (inner) solutions which exist on small scales. We find two kinds of asymptotic solutions, one of which is asymptotically flat, while the other one is not, and also two types of inner solutions, one of which displays the Vainshtein mechanism, while the other exhibits a self-shielding behaviour of the gravitational field. We analyse in detail in which cases the solutions match in an intermediate region. The asymptotically flat solutions connect only to inner configurations displaying the Vainshtein mechanism, while the non asymptotically flat solutions can connect with both kinds of inner solutions. We show furthermore that there are some regions in the parameter space where global solutions do not exist, and characterise precisely in which regions of the phase space the Vainshtein mechanism takes place.

The statistics of large-scale structure in the Universe can be used to probe non-Gaussianity of the primordial density field, complementary to existing constraints from the cosmic microwave background. In particular, the scale dependence of halo bias, which affects the halo distribution at large scales, represents a promising tool for analyzing primordial non-Gaussianity of local form. Future observations, for example, may be able to constrain the trispectrum parameter gNL that is difficult to study and constrain using the CMB alone. We investigate how galaxy and matter bispectra can distinguish between the two non-Gaussian parameters fNL and gNL, whose effects give nearly degenerate contributions to the power spectra. We use a generalization of the univariate bias approach, making the hypothesis that the number density of halos forming at a given position is a function of the local matter density contrast and of its local higher-order statistics. Using this approach, we calculate the halo-matter bispectra and analyze their properties. We determine a connection between the sign of the halo bispectrum on large scales and the parameter gNL. We also construct a combination of halo and matter bispectra that is sensitive to fNL, with little contamination from gNL. We study both the case of single and multiple sources to the primordial gravitational potential, discussing how to extend the concept of stochastic halo bias to the case of bispectra. We use a specific halo mass-function to calculate numerically the bispectra in appropriate squeezed limits, confirming our theoretical findings.

We discuss a consistent theory for a self-interacting vector field, breaking an Abelian symmetry in such a way to obtain an interesting behavior for its longitudinal polarization. In an appropriate decoupling limit, the dynamics of the longitudinal mode is controlled by Galileon interactions. The full theory away from the decoupling limit does not propagate ghost modes, and can be investigated in regimes where non-linearities become important. When coupled to gravity, this theory provides a candidate for dark energy, since it admits de Sitter cosmological solutions characterized by a technically natural value for the Hubble parameter. We also consider the homogeneous evolution when, besides the vector, additional matter in the form of perfect fluids is included. We find that the vector can have an important role in characterizing the universe expansion.

In slow roll inflation, the running of the spectral tilt is generically proportional to the square of the deviation from scale invariance, $\alpha_s \propto (n_s-1)^2$, and it is therefore currently undetectable. We present a mechanism able to generate a much larger running within slow roll. The mechanism is based on a curvaton field with a large mass term, and a time evolving normalization. This may happen for instance to the angular direction of a complex field in presence of an evolving radial direction. At the price of a single tuning between the mass term and the rate of change of the normalization, the curvaton can be made effectively light at the CMB scales, giving a spectral tilt in agreement with observations. The lightness is not preserved at later times, resulting in a detectable running of the spectral tilt. This mechanism shows that fields with a large mass term do not necessarily decouple from the inflationary physics, and provides a new tool for model building in inflation.

Gravitational vector degrees of freedom typically arise in many examples of modified gravity models. We start to systematically explore their role in these scenarios, studying the effects of coupling gravitational vector and scalar degrees of freedom. We focus on set-ups that enjoy a Galilean symmetry in the scalar sector and an Abelian gauge symmetry in the vector sector. These symmetries, together with the requirement that the equations of motion contain at most two space-time derivatives, only allow for a small number of operators in the Lagrangian for the gravitational fields. We investigate the role of gravitational vector fields for two broad classes of phenomena that characterize modified gravity scenarios. The first is self-acceleration: we analyze in general terms the behavior of vector fluctuations around self-accelerating solutions, and show that vanishing kinetic terms of vector fluctuations lead to instabilities on cosmological backgrounds. The second phenomenon is the screening of long range fifth forces by means of Vainshtein mechanism. We show that if gravitational vector fields are appropriately coupled to a spherically symmetric source, they can play an important role for defining the features of the background solution and the scale of the Vainshtein radius. Our general results can be applied to any concrete model of modified gravity, whose low-energy vector and scalar degrees of freedom satisfy the symmetry requirements that we impose.

We study linear perturbations around time dependent spherically symmetric solutions in the Lambda_3 massive gravity theory, which self-accelerate in the vacuum. We find that the dynamics of the scalar perturbations depend on the coordinate choice for the background solutions. For particular choices of coordinates there is a symmetry enhancement, leaving no propagating scalar degrees of freedom at linear order in perturbations. In contrast, any other coordinate choice propagates a single scalar mode. We find that the Hamiltonian of this scalar mode is unbounded from below for all self-accelerating solutions, signalling an instability.

Starting from the general Horndeski action, we derive the most general effective theory for scalar perturbations around flat space that allows us to screen fifth forces via the Vainshtein mechanism. The effective theory is described by a generalization of the Galileon Lagrangian, which we use to study the stability of spherically symmetric configurations exhibiting the Vainshtein effect. In particular, we discuss the phenomenological consequences of a scalar-tensor coupling that is absent in the standard Galileon Lagrangian. This coupling controls the superluminality and stability of fluctuations inside the Vainshtein radius in a way that depends on the density profile of a matter source. Particularly we find that the vacuum solution is unstable due to this coupling.

The new Planck constraints on the local bispectrum parameter fnl are about 10^5 times tighter than the current constraints on the trispectrum parameter gnl, which means that the allowed numerical values of the second and third order terms in the perturbative expansion of the curvature perturbation are comparable. We show that a consequence of this is that if gnl is large enough to be detectable, then it will induce a large variation between the observable value of fnl and its value in a larger inflated volume. Even if there were only a few extra efoldings between the beginning of inflation and horizon crossing of our Hubble horizon, an observably large gnl means that fnl is unlikely to be as small as its current constraint, regardless of its true background value. This result is very general, it holds regardless of how many fields contributed to the curvature perturbation. We also generalise this result to other shapes of non-Gaussianity, beyond the local model. We show that the variance of the 3-point function in the squeezed limit is bounded from below by the square of the squeezed limit of the 4-point function.

We study a model of multiple-field DBI inflation leading to mixed form of primordial non-Gaussianity, including equilateral and local bispectrum shapes. We present a general formalism based on the Hamilton-Jacobi approach, allowing us to go beyond slow-roll, combining the three-point function for the fields at Hubble-exit with the non-linear evolution of super-Hubble scales. We are able to obtain analytic results by taking a separable Ansatz for the Hubble rate. We find general expressions for both the equilateral and local type non-Gaussianity parameter fNL. The equilateral non-Gaussianity includes the usual enhancement for small sound speeds, but multiplied by an analytic factor which can lead to a suppression. We illustrate our results with two scenarios. In the first model, previously found to have detectable local non-Gaussianity, we find that the equilateral signal is not sufficiently suppressed to evade current observational bounds. In our second scenario we construct a model which exhibits both a detectable equilateral fNL and a negative local fNL.

Massive gravity is a good theoretical laboratory to study modifications of General Relativity. The theory offers a concrete set-up to study models of dark energy, since it admits cosmological self-accelerating solutions in the vacuum, in which the size of the acceleration depends on the graviton mass. Moreover, non-linear gravitational self-interactions, in the proximity of a matter source, manage to mimic the predictions of linearised General Relativity, hence agreeing with solar-system precision measurements. In this article, we review our work in the subject, classifying, on one hand, static solutions, and on the other hand, self-accelerating backgrounds. For what respects static solutions we exhibit black hole configurations, together with other solutions that recover General Relativity near a source via the Vainshtein mechanism. For the self-accelerating solutions we describe a wide class of cosmological backgrounds, including an analysis of their stability.

Inflation is the leading paradigm for explaining the origin of primordial density perturbations and the observed temperature fluctuations of the cosmic microwave background. However many open questions remain, in particular whether one or more scalar fields were present during inflation and how they contributed to the primordial density perturbation. We propose a new observational test of whether multiple fields, or only one (not necessarily the inflaton) generated the perturbations. We show that our test, relating the bispectrum and trispectrum, is protected against loop corrections at all orders, unlike previous relations.

Primordial perturbations with wavelengths greater than the observable universe shift the effective background fields in our observable patch from their global averages over the inflating space. This leads to a landscape picture where the properties of our observable patch depend on its location and may significantly differ from the expectation values predicted by the underlying fundamental inflationary model. We show that if multiple fields are present during inflation, this may happen even if our horizon exit would be preceded by only a few e-foldings of inflation. Non-Gaussian statistics are especially affected: for example models of local non-Gaussianity predicting |f_NL|>> 10 over the entire inflating volume can have a probability up to a few tens of percent to generate a non-detectable bispectrum in our observable patch |fNL^{obs.}|<10. In this work we establish systematic connections between the observable local properties of primordial perturbations and the global properties of the inflating space which reflect the underlying high energy physics. We study in detail the implications of both a detection and non-detection of primordial non-Gaussianity by Planck, and discover novel ways of characterising the naturalness of different observational configurations.

We propose a new model of slow-roll inflation in string cosmology, based on warped throat supergravity solutions displaying `walking' dynamics, i.e. the coupling constant of the dual gauge theory slowly varies over a range of energy scales. The features of the throat geometry are sourced by a rich field content, given by the dilaton and RR and NS fluxes. By considering the motion of a D3-brane probe in this geometry, we are able to analytically calculate the brane potential in a physically interesting regime. This potential has an inflection point: in its proximity we realize a model of inflation lasting sixty e-foldings, and whose robust predictions are in agreement with current observations. We are also able to interpret some of the most interesting aspects of this scenario in terms of the properties of the QFT dual theory.

We investigate vector contributions to the Lagrangian of $\Lambda_3-$massive gravity in the decoupling limit, the less explored sector of this theory. The main purpose is to understand the stability of maximally symmetric %self-accelerating vacuum solutions. Around self-accelerating configurations, vector degrees of freedom become strongly coupled since their kinetic terms vanish, so their dynamics is controlled by higher order interactions. Even in the decoupling limit, the vector Lagrangian contains an infinite number of terms. We develop a systematic method to covariantly determine the vector Lagrangian at each order in perturbations, fully manifesting the symmetries of the system. We show that, around self-accelerating solutions, the structure of higher order $p$-form Galileons arise, avoiding the emergence of a sixth BD ghost mode. However, a careful analysis shows that there are directions along which the Hamiltonian is unbounded from below. This instability can be interpreted as one of the available fifth physical modes behaving as a ghost. Therefore, we conclude that self-accelerating configurations, in the decoupling limit of $\Lambda_3$-massive gravity, are generically unstable.

We propose a method to probe higher-order correlators of the primordial density field through the inhomogeneity of local non-Gaussian parameters, such as f_NL, measured within smaller patches of the sky. Correlators between n-point functions measured in one patch of the sky and k-point functions measured in another patch depend upon the (n+k)-point functions over the entire sky. The inhomogeneity of non-Gaussian parameters may be a feasible way to detect or constrain higher-order correlators in local models of non-Gaussianity, as well as to distinguish between single and multiple-source scenarios for generating the primordial density perturbation, and more generally to probe the details of inflationary physics.

We study the effect of non-trivial sound speeds on local-type non-Gaussianity during multiple-field inflation. To this end, we consider a model of multiple-field DBI and use the deltaN formalism to track the super-horizon evolution of perturbations. By adopting a sum separable Hubble parameter we derive analytic expressions for the relevant quantities in the two-field case, valid beyond slow variation. We find that non-trivial sound speeds can, in principle, curve the trajectory in such a way that significant local-type non-Gaussianity is produced. Deviations from slow variation, such as rapidly varying sound speeds, enhance this effect. To illustrate our results we consider two-field inflation in the tip regions of two warped throats and find large local-type non-Gaussianity produced towards the end of the inflationary process.

We expose a new symmetry for linear perturbations around a solution of non-linear Fierz-Pauli massive gravity plus a bare cosmological constant. The cosmological constant is chosen such that the background metric is flat while the Stuckelberg fields have a non-trivial profile. Around this background, at linear order the new symmetry reduces the propagating degrees of freedom to those of General Relativity, namely the massless helicity 2 modes only. We discuss the physical consequences and possible applications of these findings.

We introduce a natural model of quintessence in string theory where the light rolling scalar is radiatively stable and couples to Standard Model matter with weaker-than- Planckian strength. The model is embedded in an anisotropic type IIB compactification with two exponentially large extra dimensions and TeV-scale gravity. The bulk turns out to be nearly supersymmetric since the scale of the gravitino mass is of the order of the observed value of the cosmological constant. The quintessence field is a modulus parameterising the size of an internal four-cycle which naturally develops a potential of the order (gravitino mass)^4, leading to a small dark energy scale without tunings. The mass of the quintessence field is also radiatively stable since it is protected by supersymmetry in the bulk. Moreover, this light scalar couples to ordinary matter via its mixing with the volume mode. Due to the fact that the quintessence field is a flat direction at leading order, this mixing is very small, resulting in a suppressed coupling to Standard Model particles which avoids stringent fifth-force constraints. On the other hand, if dark matter is realised in terms of Kaluza-Klein states, unsuppressed couplings between dark energy and dark matter can emerge, leading to a scenario of coupled quintessence within string theory. We study the dynamics of quintessence in our set-up, showing that its main features make it compatible with observations.

A generic feature of the known string inflationary models is that the same physics that makes the inflaton lighter than the Hubble scale during inflation often also makes other scalars this light. These scalars can acquire isocurvature fluctuations during inflation, and given that their VEVs determine the mass spectrum and the coupling constants of the effective low-energy field theory, these fluctuations give rise to couplings and masses that are modulated from one Hubble patch to another. These seem just what is required to obtain primordial adiabatic fluctuations through conversion into density perturbations through the `modulation mechanism', wherein reheating takes place with different efficiency in different regions of our Universe. Fluctuations generated in this way can generically produce non-gaussianity larger than obtained in single-field slow-roll inflation; potentially observable in the near future. We provide here the first explicit example of the modulation mechanism at work in string cosmology, within the framework of LARGE Volume Type-IIB string flux compactifications. The inflationary dynamics involves two light Kaehler moduli: a fibre divisor plays the role of the inflaton whose decay rate to visible sector degrees of freedom is modulated by the primordial fluctuations of a blow-up mode (which is made light by the use of poly-instanton corrections). We find the challenges of embedding the mechanism into a concrete UV completion constrains the properties of the non-gaussianity that is found, since for generic values of the underlying parameters, the model predicts a local bi-spectrum with fNL of order `a few'. However, a moderate tuning of the parameters gives also rise to explicit examples with fNL O(20) potentially observable by the Planck satellite.

We investigate strong coupling effects in a covariant massive gravity model, which is a candidate for a ghost free non-linear completion of Fierz-Pauli. We analyse the conditions to recover general relativity via Vainshtein mechanism in the weak field limit, and find three main cases depending on the choice of parameters. In the first case, the potential is such that all non-linearities disappear and the vDVZ discontinuity cannot be avoided. In the second case, the Vainshtein mechanism allows to recover general relativity within a macroscopic radius from a source. In the last case, the strong coupling of the scalar graviton completely shields the massless graviton, and weakens gravity when approaching the source. In the second part of the paper, we explore new exact vacuum solutions, that asymptote de Sitter or anti de Sitter space depending on the parameter choice. The curvature of the space is proportional to the mass of the graviton, thus providing a cosmological background which may explain the present day acceleration in terms of the graviton mass. Moreover, by expressing the potential for non-linear massive gravity in a convenient form, we also suggest possible connections with a higher dimensional framework.

We study the non-Gaussianity generated during multiple-field inflation. We provide an exact expression for the bispectrum parameter f_NL which is valid beyond the slow-roll regime, valid for certain classes of inflationary models. We then study a new, exact multi-field inflationary model considering a case where the bispectrum grows to observable values at the end of inflation. We show that in this case the trispectrum is also large and may even provide the dominant signal of non-Gaussianity.

We propose a new inflationary scenario in type IIB Calabi-Yau compactifications, where the inflaton is a K\"ahler modulus parameterising the volume of an internal four-cycle. The inflaton potential is generated via poly-instanton corrections to the superpotential which give rise to a naturally flat direction due to their double exponential suppression. Given that the volume mode is kept stable during inflation, all the inflaton-dependent higher dimensional operators are suppressed. Moreover, string loop effects can be shown to be negligible throughout all the inflationary dynamics for natural values of the underlying parameters. The model is characterised by a reheating temperature of the order $T_{\rm rh}\simeq 10^6$ GeV which requires $N_e \simeq 54$ e-foldings of inflation. All the inflationary observables are compatible with current observations since the spectral index is $n_s \simeq 0.96$, while the tensor-to-scalar ratio is $r\simeq 10^{-5}$. The volume of the Calabi-Yau is of order $10^3$ in string units, corresponding to an inflationary scale around $10^{15}$ GeV

We explore the possibility of realising self-accelerated expansion of the Universe taking into account the vector components of a massive graviton. The effective action in the decoupling limit contains an infinite number of terms, once the vector degrees of freedom are included. These can be re-summed in physically interesting situations, which result in non-polynomial couplings between the scalar and vector modes. We show there are self-accelerating background solutions for this effective action, with the possibility of having a non-trivial profile for the vector fields. We then study fluctuations around these solutions and show that there is always a ghost, if a background vector field is present. When the background vector field is switched off, the ghost can be avoided, at the price of entering into a strong coupling regime, in which the vector fluctuations have vanishing kinetic terms. Finally we show that the inclusion of a bare cosmological constant does not change the previous conclusions and it does not lead to a ghost mode in the absence of a background vector field.

We study spherically symmetric solutions in a covariant massive gravity model, which is a candidate for a ghost-free non-linear completion of the Fierz-Pauli theory. There is a branch of solutions that exhibits the Vainshtein mechanism, recovering General Relativity below a Vainshtein radius given by $(r_g m^2)^{1/3}$, where $m$ is the graviton mass and $r_g$ is the Schwarzschild radius of a matter source. Another branch of exact solutions exists, corresponding to Schwarzschild-de Sitter spacetimes where the curvature scale of de Sitter space is proportional to the mass squared of the graviton.

Inflationary correlation functions are potentially affected by infrared divergences. For example, the two-point correlator of curvature perturbation at momentum k receives corrections ~ln(kL), where L is the size of the region in which the measurement is performed. We define infrared-safe correlation functions which have no sensitivity to the size L of the box used for the observation. The conventional correlators with their familiar log-enhanced corrections (both from scalar and tensor long-wavelength modes) are easily recovered from our IR-safe correlation functions. Among other examples, we illustrate this by calculating the corrections to the non-Gaussianity parameter f_NL coming from long-wavelength tensor modes. In our approach, the IR corrections automatically emerge in a resummed, all-orders form. For the scalar corrections, the resulting all-orders expression can be evaluated explicitly.

We study the inflationary dynamics in a model of slow-roll inflation in warped throat. Inflation is realized by the motion of a D-brane along the radial direction of the throat, and at later stages instabilities develop in the angular directions. We closely investigate both the single field potential relevant for the slow-roll phase, and the full multi-field one including the angular modes which becomes important at later stages. We study the main features of the instability process, discussing its possible consequences and identifying the vacua towards which the angular modes are driven.

We calculate the scale dependence of the bispectrum and trispectrum in (quasi) local models of non-Gaussian primordial density perturbations, and characterize this scale dependence in terms of new observable parameters. They can help to discriminate between models of inflation, since they are sensitive to properties of the inflationary physics that are not probed by the standard observables. We find consistency relations between these parameters in certain classes of models. We apply our results to a scenario of modulated reheating, showing that the scale dependence of non-Gaussianity can be significant. We also discuss the scale dependence of the bispectrum and trispectrum, in cases where one varies the shape as well as the overall scale of the figure under consideration. We conclude providing a formulation of the curvature perturbation in real space, which generalises the standard local form by dropping the assumption that f_NL and g_NL are constants.

We describe a simple way of incorporating fluctuations of the Hubble scale during the horizon exit of scalar perturbations into the delta N formalism. The dominant effect comes from the dependence of the Hubble scale on low-frequency modes of the inflaton. This modifies the coefficient of the log-enhanced term appearing in the curvature spectrum at second order in field fluctuations. With this modification, the relevant coefficient turns out to be proportional to the second derivative of the tree-level spectrum with respect to the inflaton phi at horizon exit. A logarithm with precisely the same coefficient appears in a calculation of the log-enhancement of the curvature spectrum based purely on the geometry of the reheating surface. We take this agreement as strong support for the proposed implementation of the delta N formalism. Moreover, our analysis makes it apparent that the log-enhancement of the inflationary power-spectrum is indeed physical if this quantity is defined using a global coordinate system on the reheating surface (or any other post-inflationary surface of constant energy density). However, it can be avoided by defining the spectrum using invariant distances on this surface.

Inflationary scenarios in string theory often involve a large number of light scalar fields, whose presence can enrich the post-inflationary evolution of primordial fluctuations generated during the inflationary epoch. We provide a simple example of such post-inflationary processing within an explicit string-inflationary construction, using a Kahler modulus as the inflaton within the framework of LARGE Volume Type-IIB string flux compactifications. We argue that inflationary models within this broad category often have a selection of scalars that are light enough to be cosmologically relevant, whose contributions to the primordial fluctuation spectrum can compete with those generated in the standard way by the inflaton. These models consequently often predict nongaussianity at a level, f_NL ~ O(10), potentially observable by the Planck satellite, with a bi-spectrum maximized by triangles with squeezed shape in a string realisation of the curvaton scenario. We argue that the observation of such a signal would robustly prefer string cosmologies such as these that predict a multi-field dynamics during the very early universe.

We consider possible scale-dependence of the non-linearity parameter f_NL in local and quasi-local models of non-Gaussian primordial density perturbations. In the simplest model where the primordial perturbations are a quadratic local function of a single Gaussian field then f_NL is scale-independent by construction. However scale-dependence can arise due to either a local function of more than one Gaussian field, or due to non-linear evolution of modes after horizon-exit during inflation. We show that the scale dependence of f_NL is typically first order in slow-roll. For some models this may be observable with experiments such as Planck provided that f_NL is close to the current observational bounds.

We adopt a covariant formalism to derive exact evolution equations for nonlinear perturbations, in a universe dominated by two scalar fields. These scalar fields are characterized by non-canonical kinetic terms and an arbitrary field space metric, a situation typically encountered in inflationary models inspired by string theory. We decompose the nonlinear scalar perturbations into adiabatic and entropy modes, generalizing the definition adopted in the linear theory, and we derive the corresponding exact evolution equations. We also obtain a nonlinear generalization of the curvature perturbation on uniform density hypersurfaces, showing that on large scales it is sourced only by the nonlinear version of the entropy perturbation. We then expand these equations to second order in the perturbations, using a coordinate based formalism. Our results are relatively compact and elegant and enable one to identify the new effects coming from the non-canonical structure of the scalar fields Lagrangian. We also explain how to analyze, in our formalism, the interesting scenario of multifield Dirac-Born-Infeld inflation.

We explore models of intersecting brane-worlds with induced gravity terms on codimension one branes and on their intersection. Maximally symmetric solutions for the branes and the intersection are found. We find new self-accelerating solutions. In a 6d spacetime, the solutions realize the see-saw modification of gravity where the UV scale of the modification to 4d gravity is determined by 6d Planck scale given by $M_6 \sim 10^{-3}$eV and the IR scale of the modification is determined by $M_6^2/M_4 \sim H_0 \sim 10^{-42}$ GeV where $H_0$ is present-day Hubble scale. We find that it is increasingly difficult to construct phenomenologically viable models in higher-dimensional spacetime due to the necessity to have the lower value for the fundamental Planck scale to realize the late time acceleration. It is found that the system also admits self-tuning solutions where the tension at the intersection does not change the geometry of the intersection. The induced gravity terms can avoid the necessity to compactify the extra dimensions. Finally, we discuss the possibility to have ordinary matter at the intersection, without introducing any regularisation, using the induced gravity terms.

We explore cosmology of intersecting braneworlds with induced gravity on the branes. We find the cosmological equations that control the evolution of a moving codimension one brane and a codimension two brane that sits at the intersection. We study the Friedmann equation at the intersection, finding new contributions from the six dimensional bulk. These higher dimensional contributions allow us to find new examples of self-accelerating configurations for the codimension two brane at the intersection and we discuss their features.

We study the cosmological implications of including angular motion in the DBI brane inflation scenario. The non-canonical kinetic terms of the Dirac-Born-Infeld action give an interesting alternative to slow roll inflation, and cycling branes can drive periods of accelerated expansion in the Universe. We present explicit numerical solutions demonstrating brane inflation in the Klebanov-Strassler throat. We find that demanding sufficient inflation takes place in the throat is in conflict with keeping the brane's total energy low enough so that local gravitational backreaction on the Calabi-Yau manifold can be safely ignored. We deduce that spinflation (brane inflation with angular momentum) can ease this tension by providing extra e-foldings at the start of inflation. Cosmological expansion rapidly damps the angular momentum causing an exit to a more conventional brane inflation scenario. Finally, we set up a general framework for cosmological perturbation theory in this scenario, where we have multi-field non-standard kinetic term inflation.

We analyse the dynamics of a probe D3-(anti-)brane propagating in a warped string compactification, making use of the Dirac-Born-Infeld action approximation. We also examine the time dependent expansion of such moving branes from the ``mirage cosmology'' perspective, where cosmology is induced by the brane motion in the background spacetime. A range of physically interesting backgrounds are considered: AdS5, Klebanov-Tseytlin and Klebanov-Strassler. Our focus is on exploring what new phenomenology is obtained from giving the brane angular momentum in the extra dimensions. We find that in general, angular momentum creates a centrifugal barrier, causing bouncing cosmologies. More unexpected, and more interesting, is the existence of bound orbits, corresponding to cyclic universes.

We study the cosmological properties of a codimension two brane world that sits at the intersection between two four branes, in the framework of six dimensional Einstein-Gauss-Bonnet gravity. Due to contributions of the Gauss-Bonnet terms, the junction conditions require the presence of localized energy density on the codimension two defect. The induced metric on this surface assumes a FRW form, with a scale factor associated to the position of the brane in the background; we can embed on the codimension two defect the preferred form of energy density. We present the cosmological evolution equations for the three brane, showing that, for the case of pure AdS$_6$ backgrounds, they acquire the same form of the ones for the Randall-Sundrum II model. When the background is different from pure AdS$_6$, the cosmological behavior is potentially modified in respect to the typical one of codimension one brane worlds. We discuss, in a particular model embedded in an AdS$_6$ black hole, the conditions one should satisfy in order to obtain standard cosmology at late epochs.

We re-consider the self tuning idea in brane world models of finite volume. We notice that in a large class of self tuning models, the four dimensional physics is sensitive to the vacuum energy on the brane. In particular the compactification volume changes each time the tension of the brane is modified: consequently, observable constants, as the effective Planck mass or masses of matter fields, change as well. We notice that the self tuning mechanism and the stabilization mechanism of the size of the extra dimensions are generically in apparent conflict. We focus on a self tuning model in six spacetime dimensions to analyze how the above considerations are explicitely realized.

We present several higher-dimensional spacetimes for which observers living on 3-branes experience an induced metric which bounces. The classes of examples include boundary branes on generalised S-brane backgrounds and probe branes in D-brane/anti D-brane systems. The bounces we consider normally would be expected to require an energy density which violates the weak energy condition, and for our co-dimension one examples this is attributable to bulk curvature terms in the effective Friedmann equation. We examine the features of the acceleration which provides the bounce, including in some cases the existence of positive acceleration without event horizons, and we give a geometrical interpretation for it. We discuss the stability of the solutions from the point of view of both the brane and the bulk. Some of our examples appear to be stable from the bulk point of view, suggesting the possible existence of stable bouncing cosmologies within the brane-world framework.

In supersymmetric theories, the occurrence of universal singlets is a delicate issue, because they usually induce tadpoles that destabilize the hierarchy. We study the effects of these tadpoles in supersymmetric hybrid inflation models. The resulting scenario is generically modified, but it is still possible to achieve inflation in a natural way. It is argued that singlets, despite the problems associated with their presence, can lead to interesting cosmological consequences.