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.

We present the first application to real data from the BOSS DR12 dataset of the Extractor procedure to determine the acoustic scale imprinted on Baryonic Acoustic Oscillations (BAO). We show that, being largely insensitive to the broadband shape of the Power Spectrum, this procedure requires a lower number of nuisance parameters than those used by the BOSS collaboration, For non-reconstructed data our analysis improves the accuracy on the acoustic scale by about 20 %, while for reconstructed ones we get essentially the same level of accuracy as the BOSS analysis.

The Einstein Telescope (ET), a proposed European ground-based gravitational-wave detector of third-generation, is an evolution of second-generation detectors such as Advanced LIGO, Advanced Virgo, and KAGRA which could be operating in the mid 2030s. ET will explore the universe with gravitational waves up to cosmological distances. We discuss its main scientific objectives and its potential for discoveries in astrophysics, cosmology and fundamental physics.

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.

Motivated by some of the recent swampland conjectures, we study the implementation for the late time acceleration of the Universe of a mechanism developed by Anber and Sorbo in the context of primordial inflation, in which an axion field can slowly roll in a steep potential due to additional friction provided by its coupling to some U(1) gauge field. We first study the realization of this mechanism in N = 2 supergravity models resulting from string compactifications on Calabi--Yau manifolds. We then study the transition between matter domination and the axion domination, and show that indeed the backreaction of the produced gauge field can sufficiently slow the motion of the axion, so to produce the present accelerated era. We finally study the transition from a pre-inflationary matter or radiation domination to primordial inflation. In the regime that we could explore numerically, the evolution is characterized by stages of faster axion roll (and consequent bursts of gauge field amplification) intermitted by stages of slower roll, with a pattern that "oscillates'' about the steady state Anber and Sorbo solution, but that does not appear to relax to it.

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.

We study the production of fermions through a derivative coupling with a pseudoscalar inflaton and the effects of the produced fermions on the scalar primordial perturbations. We present analytic results for the modification of the scalar power spectrum due to the produced fermions, and we estimate the amplitude of the non-Gaussianities in the equilateral regime. Remarkably, we find a regime where the effect of the fermions gives the dominant contribution to the scalar spectrum while the amplitude of the bispectrum is small and in agreement with observation. We also note the existence of a regime in which the backreaction of the fermions on the evolution of the zero-mode of the inflaton can lead to inflation even if the potential of the inflaton is steep and does not satisfy the slow-roll conditions.

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.

We study the formation of primordial black holes when they are generated by the collapse of large overdensities in the early universe. Since the density contrast is related to the comoving curvature perturbation by a nonlinear relation, the overdensity statistics is unavoidably non-Gaussian. We show that the abundance of primordial black holes at formation may not be captured by a perturbative approach which retains the first few cumulants of the non-Gaussian probability distribution. We provide two techniques to calculate the non-Gaussian abundance of primordial black holes at formation, one based on peak theory and the other on threshold statistics. Our results show that the unavoidable non-Gaussian nature of the inhomogeneities in the energy density makes it harder to generate PBHs. We provide simple (semi-)analytical expressions to calculate the non-Gaussian abundances of the primordial black holes and show that for both narrow and broad power spectra the gaussian case from threshold statistics is reproduced by increasing the amplitude of the power spectrum by a factor ${\cal O}(2\div 3)$.

We study a variant of the Chromo-Natural Inflation (CNI) mechanism in which the inflaton interacts only gravitationally with the CNI fields. Integrating out all the non-dynamical scalar fields of the model results in a coupling between the perturbations of the inflaton and of the CNI pseudo-scalar which is significantly greater than the one obtained in the absence of the gauge CNI dynamics. We compute how this greater coupling impacts the power spectrum of the inflaton perturbations that are sourced nonlinearly by the unstable (tensor) gauge CNI modes, and we require that the amplitude of these modes is well below that of the linear perturbations. Combining this result with various constraints, including backreaction effects, the requirement of having observable and dominant sourced gravitational waves (GW), and the current upper bound on the tensor-to-scalar ratio, significantly constrains the range of parameter space where this model can produce an interesting GW signal.

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.

The idea that primordial black holes (PBHs) can comprise most of the dark matter of the universe has recently reacquired a lot of momentum. Observational constraints, however, rule out this possibility for most of the PBH masses, with a notable exception around $10^{-12} M_\odot$. These light PBHs may be originated when a sizeable comoving curvature perturbation generated during inflation re-enters the horizon during the radiation phase. During such a stage, it is unavoidable that gravitational waves (GWs) are generated. Since their source is quadratic in the curvature perturbations, these GWs are generated fully non-Gaussian. Their frequency today is about the mHz, which is exactly the range where the LISA mission has the maximum of its sensitivity. This is certainly an impressive coincidence. We show that this scenario of PBHs as dark matter can be tested by LISA by measuring the GW two-point correlator. On the other hand, we show that the short observation time (as compared to the age of the universe) and propagation effects of the GWs across the perturbed universe from the production point to the LISA detector suppress the bispectrum to an unobservable level. This suppression is completely general and not specific to our model.

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.

We present analytic results for the gravitational wave power spectrum induced in models where the inflaton is coupled to a fermionic pseudocurrent. We show that although such a coupling creates helically polarized fermions, the polarized component of the resulting gravitational waves is parametrically suppressed with respect to the non-polarized one. We also show that the amplitude of the gravitational wave signal associated to this production cannot exceed that generated by the standard mechanism of amplification of vacuum fluctuations. We previously found that this model allows for a regime in which the backreaction of the produced fermions allows for slow-roll inflation even for a steep inflaton potential, and still leads to Gaussian primordial scalar perturbations. The present analysis shows that this regime also results in a gravitational wave signal compatible with the current bounds.

There has recently been renewed interest in the possibility that the dark matter in the universe consists of primordial black holes (PBHs). Current observational constraints leave only a few PBH mass ranges for this possibility. One of them is around $10^{-12} M_\odot$. If PBHs with this mass are formed due to an enhanced scalar-perturbation amplitude, their formation is inevitably accompanied by the generation of gravitational waves (GWs) with frequency peaked in the mHz range, precisely around the maximum sensitivity of the LISA mission. We show that, if these primordial black holes are the dark matter, LISA will be able to detect the associated GW power spectrum. Although the GW source signal is intrinsically non-Gaussian, the signal measured by LISA is a sum of the signal from a large number of independent sources suppressing the non-Gaussianity at detection to an unobservable level. We also discuss the effect of the GW propagation in the perturbed universe. PBH dark matter generically leads to a detectable, purely isotropic, Gaussian and unpolarised GW signal, a prediction that is testable with LISA.

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

Several models of inflation employing a triplet of SU(2) vectors with spatially orthogonal vacuum expectation values (VEVs) have been recently proposed. One (tensor) combination $t$ of the vector modes is amplified in some momentum range during inflation. Due to the vector VEVs, this combination mixes with gravitational waves (GW) at the linear level, resulting in a GW amplification that has been well studied in the literature. Scalar perturbations in this class of models have been so far studied only at the linear level. We perform a first step toward the nonlinear computation using as an example the original model of Chromo-Natural Inflation. We compute the contribution to the scalar power spectrum arising from the coupling of the combination $t$ to the inflaton. This contribution is mostly controlled by a single parameter of the model (namely, the ratio between the mass of the fluctuations of the vector field and the Hubble rate), and, for a wide range of this parameter, it can significantly affect the phenomenology obtained from the linear theory. This nonlinear contribution is significantly blue, improving the comparison between the two-point function and the Cosmic Microwave Background (CMB) data. This growth can be also relevant for smaller scale phenomenology, such as large scale structure, CMB distortions, and primordial black holes.

We generalize models invoking a spin-2 particle as a mediator between the dark sector and the Standard Model. We show that a massive spin-2 messenger can efficiently play the role of a portal between the two sectors. The dark matter is then produced via a freeze-in mechanism during the reheating epoch. In a large part of the parameter space, production through the exchange of a massive spin-2 mediator dominates over processes involving a graviton with Planck suppressed couplings. We perform a systematic analysis of such models for different values of the spin-2 mass relative to the maximum and the final temperature attained at reheating.

We study a new procedure to measure the sound horizon scale via Baryonic Acoustic Oscillations (BAO). Instead of fitting the measured power spectrum (PS) to a theoretical model containing the cosmological informations and all the nonlinear effects, we define a procedure to project out (or to "extract") the oscillating component from a given nonlinear PS. We show that the BAO scale extracted in this way is extremely robust and, moreover, can be reproduced by simple theoretical models at any redshift. By using N-body simulations, we discuss the effect of the nonlinear evolution of the matter field, of redshift space distortions and of scale-dependent halo bias, showing that all these effects can be reproduced with sub-percent accuracy. We give a one-parameter theoretical model based on a simple (IR) modification of 1-loop perturbation theory, which reproduces the BAO scale from measurements of halo clustering in redshift space at better than $0.1\%$ level and does not need any external UV input, such as coefficients measured from N-body simulations.

We compute the gravitational leptogenesis generated from the parity-violating gravitational waves sourced by an abelian gauge field coupled to a pseudo-scalar inflation. We show that, once the CMB bound on the tensor-to-scalar ratio is enforced, the lepton asymmetry produced by this mechanism during inflation is too small to account for the observed baryon asymmetry of the universe, irrespectively of the inflaton potential, the strength of its coupling to the gauge field, and the details of reheating.

In the first stages of inflationary reheating, the temperature of the radiation produced by inflaton decays is typically higher than the commonly defined reheating temperature $T_{RH} \sim (\Gamma_\phi M_P)^{1/2}$ where $\Gamma_\phi$ is the inflaton decay rate. We consider the effect of particle production at temperatures at or near the maximum temperature attained during reheating. We show that the impact of this early production on the final particle abundance depends strongly on the temperature dependence of the production cross section. For $\langle \sigma v \rangle \sim T^n/M^{n+2}$, and for $n < 6$, any particle produced at $T_{\rm max}$ is diluted by the later generation of entropy near $T_{RH}$. This applies to cases such as gravitino production in low scale supersymmetric models ($n=0$) or NETDM models of dark matter ($n=2$). However, for $n\ge6$ the net abundance of particles produced during reheating is enhanced by over an order of magnitude, dominating over the dilution effect. This applies, for instance to gravitino production in high scale supersymmetry models where $n=6$.

We define a procedure to extract the oscillating part of a given nonlinear Power Spectrum, and derive an equation describing its evolution including the leading effects at all scales. The intermediate scales are taken into account by standard perturbation theory, the long range (IR) displacements are included by using consistency relations, and the effect of small (UV) scales is included via effective coefficients computed in simulations. We show that the UV effects are irrelevant in the evolution of the oscillating part, while they play a crucial role in reproducing the smooth component. Our "extractor" operator can be applied to simulations and real data in order to extract the Baryonic Acoustic Oscillations (BAO) without any fitting function and nuisance parameter. We conclude that the nonlinear evolution of BAO can be accurately reproduced at all scales down to $z=0$ by our fast analytical method, without any need of extra parameters fitted from simulations.

Primordial Black Holes (PBH) could be the cold dark matter of the universe. They could have arisen from large (order one) curvature fluctuations produced during inflation that reentered the horizon in the radiation era. At reentry, these fluctuations source gravitational waves (GW) via second order anisotropic stresses. These GW, together with those (possibly) sourced during inflation by the same mechanism responsible for the large curvature fluctuations, constitute a primordial stochastic GW background (SGWB) that unavoidably accompanies the PBH formation. We study how the amplitude and the range of frequencies of this signal depend on the statistics (Gaussian versus $\chi^2$) of the primordial curvature fluctuations, and on the evolution of the PBH mass function due to accretion and merging. We then compare this signal with the sensitivity of present and future detectors, at PTA and LISA scales. We find that this SGWB will help to probe, or strongly constrain, the early universe mechanism of PBH production. The comparison between the peak mass of the PBH distribution and the peak frequency of this SGWB will provide important information on the merging and accretion evolution of the PBH mass distribution from their formation to the present era. Different assumptions on the statistics and on the PBH evolution also result in different amounts of CMB $\mu$-distortions. Therefore the above results can be complemented by the detection (or the absence) of $\mu$-distortions with an experiment such as PIXIE.

We revisit the study of the phenomenology associated to a burst of particle production of a field whose mass is controlled by the inflaton field and vanishes at one given instance during inflation. This generates a bump in the correlators of the primordial scalar curvature. We provide a unified formalism to compute various effects that have been obtained in the literature and confirm that the dominant effects are due to the rescattering of the produced particles on the inflaton condensate. We improve over existing results (based on numerical fits) by providing exact analytic expressions for the shape and height of the bump, both in the power spectrum and the equilateral bispectrum. We then study the regime of validity of the perturbative computations of this signature. Finally, we extend these computations to the case of a burst of particle production in a sector coupled only gravitationally to the inflaton.

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 study the prospects of detection at terrestrial and space interferometers, as well as at pulsar timing array experiments, of a stochastic gravitational wave background which can be produced in models of axion inflation. This potential signal, and the development of these experiments, open a new window on inflation on scales much smaller than those currently probed with Cosmic Microwave Background and Large Scale Structure measurements. The sourced signal generated in axion inflation is an ideal candidate for such searches, since it naturally grows at small scales, and it has specific properties (chirality and non-gaussianity) that can distinguish it from an astrophysical background. We study under which conditions such a signal can be produced at an observable level, without the simultaneous overproduction of scalar perturbations in excess of what is allowed by the primordial black hole limits. We also explore the possibility that scalar perturbations generated in a modified version of this model may provide a distribution of primordial black holes compatible with the current bounds, that can act as a seeds of the present black holes in the universe.

Trapped inflation is a mechanism in which particle production from the moving inflaton is the main source of friction in the inflaton equation of motion. The produced fields source inflaton perturbations, which dominate over the vacuum ones. We employ the set of equations for the inflaton zero mode and its perturbations which was developed in the original work on trapped inflation, and which we extend to second order in the perturbations. We build on this study by updating the experimental constraints, and by replacing the existing approximate solutions with more accurate ones. We obtain a different numerical value for the amplitude of the scalar power spectrum, and a parametrically different result for the bispectrum. This has implications for the allowed region of parameter space in models of trapped inflation, and for some of the phenomenological results obtained in this region. The main results in the allowed region are the following: monomial inflaton potentials, such as $V \propto \varphi ,\, \varphi^2$ can be compatible with the data, and (in a portion of the allowed region) the inflaton can be sub-Planckian over all the "observable" stage of inflation; gravitational waves from this mechanism are too small to be observed in the foreseeable future; a non-gaussianity parameter $-60 < f_{\rm NL} < -20 $ is obtained in the allowed region for exactly equilateral configurations, with a mild dependence on the parameters of the model and of the inflaton potential.

This book lays out the scientific goals to be addressed by the next-generation ground-based cosmic microwave background experiment, CMB-S4, envisioned to consist of dedicated telescopes at the South Pole, the high Chilean Atacama plateau and possibly a northern hemisphere site, all equipped with new superconducting cameras. CMB-S4 will dramatically advance cosmological studies by crossing critical thresholds in the search for the B-mode polarization signature of primordial gravitational waves, in the determination of the number and masses of the neutrinos, in the search for evidence of new light relics, in constraining the nature of dark energy, and in testing general relativity on large scales.

Many of the methods proposed so far to go beyond Standard Perturbation Theory break invariance under time-dependent boosts (denoted here as extended Galilean Invariance, or GI). This gives rise to spurious large scale effects which spoil the small scale predictions of these approximation schemes. By using consistency relations we derive fully non-perturbative constraints that GI imposes on correlation functions. We then introduce a method to quantify the amount of GI breaking of a given scheme, and to correct it by properly tailored counterterms. Finally, we formulate resummation schemes which are manifestly GI, discuss their general features, and implement them in the so called Time-Flow, or TRG, equations.

The motion of a pseudo-scalar field $X$ during inflation naturally induces a significant amplification of the gauge fields to which it is coupled. The amplified gauge fields can source characteristic scalar and tensor primordial perturbations. Several phenomenological implications have been discussed in the cases in which (i) $X$ is the inflation, and (ii) $X$ is a field different from the inflation, that experiences a temporary speed up during inflation. In this second case, visible sourced gravitational waves (GW) can be produced at the CMB scales without affecting the scalar perturbations, even if the scale of inflation is several orders of magnitude below what is required to produce a visible vacuum GW signal. Perturbativity considerations can be used to limit the regime in which these results are under perturbative control. We revised limits recently claimed for the case (i), and we extend these considerations to the case (ii). We show that, in both cases, these limits are satisfied by the applications that generate signals at CMB scales. Applications that generate gravitational waves and primordial black holes at much smaller scales are at the limit of the validity of this perturbativity analysis, so we expect those results to be valid up to possibly order one corrections.

We revisit gravitino production following inflation. As a first step, we review the standard calculation of gravitino production in the thermal plasma formed at the end of post-inflationary reheating when the inflaton has completely decayed. Next we consider gravitino production prior to the completion of reheating, assuming that the inflaton decay products thermalize instantaneously while they are still dilute. We then argue that instantaneous thermalization is in general a good approximation, and also show that the contribution of non-thermal gravitino production via the collisions of inflaton decay products prior to thermalization is relatively small. Our final estimate of the gravitino-to-entropy ratio is approximated well by a standard calculation of gravitino production in the post-inflationary thermal plasma assuming total instantaneous decay and thermalization at a time $t \simeq 1.2/\Gamma_\phi$. Finally, in light of our calculations, we consider potential implications of upper limits on the gravitino abundance for models of inflation, with particular attention to scenarios for inflaton decays in supersymmetric Starobinsky-like models.

We consider a model in which a pseudo-scalar field $\sigma$ rolls for some e-folds during inflation, sourcing one helicity of a gauge field. These fields are only gravitationally coupled to the inflaton, and therefore produce scalar and tensor primordial perturbations only through gravitational interactions. These sourced signals are localized on modes that exit the horizon while the roll of $\sigma$ is significant. We focus our study on cases in which the model can simultaneously produce (i) a large gravitational wave signal, resulting in observable B-modes of the CMB polarizations, and (ii) sufficiently small scalar perturbations, so to be in agreement with the current limits from temperature anisotropies. Different choice of parameters can instead lead to a localized and visible departure from gaussianity in the scalar sector, either at CMB or LSS scales.

We examine the cosmological correlators induced by the simultaneous breaking of parity and of statistical isotropy, e.g., in presence of the coupling ${\cal L} = f(\phi) ( - \frac{1}{4} F^2 + \frac{\gamma}{4} F \tilde{F} )$ between the inflaton $\phi$ and a vector field with vacuum expectation value ${\bf A}$. For a suitably chosen function $f$, the energy in the vector field $\rho_{\rm A}$ does not decay during inflation. This results in nearly scale-invariant signatures of broken statistical isotropy and parity. Specifically, we find that the scalar-scalar correlator of primordial curvature perturbations includes a quadrupolar anisotropy, $P_\zeta ( {\bf k}) = P(k)[ 1 + g_* ( \hat{\bf k} \cdot \hat{\bf A})^2]$, and a (angle-averaged) scalar bispectrum that is a linear combination of the first $3$ Legendre polynomials, $B_\zeta(k_1, k_2, k_3) = \sum_L c_L P_L (\hat{\bf k}_1 \cdot \hat{\bf k}_2) P(k_1) P(k_2) + 2~{\rm perms} $, with $c_0 : c_1 : c_2 = 2 : -3 : 1$ ($c_1 \neq 0$ is a consequence of parity violation, corresponding to the constant $\gamma \neq 0$). The latter is one of the main results of this paper, which provides for the first time a clear example of an inflationary model where a non-negligible $c_1$ contribution to the bispectrum is generated. The scalar-tensor and tensor-tensor correlators induce characteristic signatures in the Cosmic Microwave Background temperature anisotropies (T) and polarization (E/B modes); namely, non-diagonal contributions to $\langle a_{\ell_1 m_1} a_{\ell_2 m_2}^* \rangle$, with $|\ell_1 - \ell_2| = 1$ in TT, TE, EE and BB, and $|\ell_1 - \ell_2| = 2$ in TB and EB. The latest CMB bounds on the scalar observables ($g_*$, $c_0$, $c_1$ and $c_2$), translate into the upper limit $\rho_{\rm A} / \rho_\phi \lesssim 10^{-9}$ at $\gamma=0$. We find that the upper limit on the vector energy density becomes much more stringent as $\gamma$ grows.

During inflation, scalar fields, including the Higgs boson, may acquire a nonzero vacuum expectation value, which must later relax to the equilibrium value during reheating. In the presence of the time-dependent condensate, the vacuum state can evolve into a state with a nonzero particle number. We show that, in the presence of lepton number violation in the neutrino sector, the particle production can explain the observed matter-antimatter asymmetry of the universe. We find that this form of leptogenesis is particularly effective when the Higgs condensate decays rapidly and at low reheat scale. As part of the calculation, we present some exact results for the Bogoliubov transformations for Majorana fermions with a nonzero time-dependent chemical potential, in addition to a time-dependent mass.

We study the impact of neutrino masses on the shape and height of the BAO peak of the matter correlation function, both in real and redshift space. In order to describe the nonlinear evolution of the BAO peak we run N-body simulations and compare them with simple analytic formulae. We show that the evolution with redshift of the correlation function and its dependence on the neutrino masses is well reproduced in a simplified version of the Zel'dovich approximation, in which the mode-coupling contribution to the power spectrum is neglected. While in linear theory the BAO peak decreases for increasing neutrino masses, the effect of nonlinear structure formation goes in the opposite direction, since the peak broadening by large scale flows is less effective. As a result of this combined effect, the peak decreases by $\sim 0.6 \%$ for $ \sum m_\nu = 0.15$ eV and increases by $\sim1.2 \%$ for $ \sum m_\nu = 0.3$ eV, with respect to a massless neutrino cosmology with equal value of the other cosmological parameters. We extend our analysis to redshift space and to halos, and confirm the agreement between simulations and the analytic formulae. We argue that all analytical approaches having the Zel'dovich propagator in their lowest order approximation should give comparable performances, irrespectively to their formulation in Lagrangian or in Eulerian space.

In Aligned Natural Inflation, an alignment between different potential terms produces an inflaton excursion greater than the axion scales in the potential. We show that, starting from a general potential of two axions with two aligned potential terms, the effective theory for the resulting light direction is characterized by four parameters: an effective potential scale, an effective axion constant, and two extra parameters (related to ratios of the axion scales and the potential scales in the $2-$field theory). For all choices of these extra parameters, the model can support inflation along valleys (in the $2-$field space) that end in minima of the potential. This leads to a phenomenology similar to that of single field Natural Inflation. For a significant range of the extra two parameters, the model possesses also higher altitude inflationary trajectories passing through saddle points of the $2-$field potential, and disconnected from any minimum. These plateaus end when the heavier direction becomes unstable, and therefore all of inflation takes place close to the saddle point, where - due to the higher altitude - the potential is flatter (smaller $\epsilon$ parameter). As a consequence, a tensor-to-scalar ratio $r = {\rm O } \left( 10^{-4} - 10^{-2} \right)$ can be easily achieved in the allowed $n_s$ region, well within the latest $1 \sigma$ CMB contours.

Motivated by the coincidence between the Hubble scale during inflation and the typical see-saw neutrino mass scale, we present a supergravity model where the inflaton is identified with a linear combination of right-handed sneutrino fields. The model accommodates an inflaton potential that is flatter than quadratic chaotic inflation, resulting in a measurable but not yet ruled out tensor-to-scalar ratio. Small CP-violation in the neutrino mass matrix and supersymmetry breaking yield an evolution in the complex plane for the sneutrino fields. This induces a net lepton charge that, via the Affleck-Dine mechanism, can be the origin of the observed baryon asymmetry of the universe.

A pseudo-scalar inflaton field can have interesting phenomenological signatures associated with parity violation. The existing analyses of these signatures typically assume statistical isotropy. In the present work we instead investigate the possibility that a pseudo-scalar inflaton is coupled to a vector field carrying a small but non-negligible vacuum expectation value (vev) coherent over our Hubble patch. We show that, in such case, correlators involving the primordial curvature perturbations and gravitational waves violate both statistical isotropy and parity symmetry. We compute the Cosmic Microwave Background (CMB) temperature anisotropies (T) and polarization (E/B) generated by these primordial modes. The CMB two-point correlation functions present distinct signals of broken rotational and parity invariance. Specifically, we find non-vanishing TT, TE, EE and BB correlators between $\ell_1$ and $\ell_2 = \ell_1 \pm 1$ multipoles, and non-vanishing TB and EB correlators between $\ell_1$ and $\ell_2 = \ell_1 \pm 2$ multipoles. Such signatures are specific of the models under consideration and they cannot be generated if one of parity and isotropy is preserved. As a specific example we consider the simple case in which the vector field has just an "electric" background component decaying in the standard way as $a^{-2}$. In this case a strong scale-dependent quadrupolar modulation of the primordial power spectra is generated and we find that almost noiseless data of the large-scale temperature and E-mode polarization anisotropies (like, e.g., the ones provided by WMAP or $Planck$) should be able to constrain the quadrupolar amplitude coefficients $g_{2M}$ of the primordial scalar power spectrum (normalized at the pivot scale comparable to the present horizon size $k_0^{-1} = 14~{\rm Gpc}$) down to $g_{2M} = 30$ (68%CL).

Standard cosmological perturbation theory (SPT) for the Large Scale Structure (LSS) of the Universe fails at small scales (UV) due to strong nonlinearities and to multistreaming effects. In Pietroni et al. 2011 a new framework was proposed in which the large scales (IR) are treated perturbatively while the information on the UV, mainly small scale velocity dispersion, is obtained by nonlinear methods like N-body simulations. Here we develop this approach, showing that it is possible to reproduce the fully nonlinear power spectrum (PS) by combining a simple (and fast) 1-loop computation for the IR scales and the measurement of a single, dominant, correlator from N-body simulations for the UV ones. We measure this correlator for a suite of seven different cosmologies, and we show that its inclusion in our perturbation scheme reproduces the fully non-linear PS with percent level accuracy, for wave numbers up to $k\sim 0.4\, h~{\rm Mpc^{-1}}$ down to $z=0$. We then show that, once this correlator has been measured in a given cosmology, there is no need to run a new simulation for a different cosmology in the suite. Indeed, by rescaling this correlator by a proper function computable in SPT, the reconstruction procedure works also for the other cosmologies and for all redshifts, with comparable accuracy. Finally, we clarify the relation of this approach to the Effective Field Theory methods recently proposed in the LSS context.

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

Solid inflation is an effective field theory of inflation in which isotropy and homogeneity are accomplished via a specific combination of anisotropic sources (three scalar fields that individually break isotropy). This results in specific observational signatures that are not found in standard models of inflation: a non-trivial angular dependence for the squeezed bispectrum, and a possibly long period of anisotropic inflation (to drive inflation, the "solid" must be very insensitive to any deformation, and thus background anisotropies are very slowly erased). In this paper we compute the expected level of statistical anisotropy in the power spectrum of the curvature perturbations of this model. To do so, we account for the classical background values of the three scalar fields that are generated on large (superhorizon) scales during inflation via a random walk sum, as the perturbation modes leave the horizon. Such an anisotropy is unavoidably generated, even starting from perfectly isotropic classical initial conditions. The expected level of anisotropy is related to the duration of inflation and to the amplitude of the squeezed bispectrum. If this amplitude is close to its current observational limit (so that one of the most interesting predictions of the model can be observed in the near future), we find that a level of statistical anisotropy $\gtrsim 3\%$ in the power spectrum is to be expected, if inflation lasted $\gtrsim 20-30$ e-folds more than the final $50-60$ efolds required to generare the CMB modes. We also comment and point out various similarities between solid inflation and models of inflation where a suitable coupling of the inflaton to a vector kinetic term $F^{2}$ gives frozen and scale invariant vector perturbations on superhorizon scales.

We discuss a mechanism of particle production during inflation that can result in a blue gravitational wave (GW) spectrum, compatible with the BICEP2 result and with the r < 0.11 limit on the tensor-to-scalar ratio at the Planck pivot scale. The mechanism is based on the production of vector quanta from a rolling pseudo-scalar field. Both the vector and the pseudo-scalar are only gravitationally coupled to the inflaton, to keep the production of inflaton quanta at an unobservable level (the overproduction of non-gaussian scalar perturbations is a generic difficulty for mechanisms that aim to generate a visible GW signal from particle production during inflation). This mechanism can produce a detectable amount of GWs for any inflationary energy scale. The produced GWs are chiral and non-gaussian; both these aspects can be tested with large-scale polarization data (starting from Planck). We study how to reconstruct the pseudo-scalar potential from the GW spectrum.

The BICEP2 collaboration has reported a strong B mode signal in the CMB polarization, which is well fit by a tensor-to-scalar ratio of r ~ 0.2. This is greater than the upper limit r < 0.11 obtained from the temperature anisotropies under the assumption of a constant scalar spectral index n_s. This discrepancy can be reduced once the statistical error and the contamination from polarized dust are accounted for. If however a large value for r will be confirmed, it will need to be reconciled with the temperature anisotropies data. The most advocated explanation involves a variation of n_s with scales (denoted as running) that has a magnitude significantly greater than the generic slow roll predictions. We instead study the possibility that the large scale temperature anisotropies are not enhanced because of a suppression of the scalar power at large scales. Such a situation can be achieved for instance by a sudden change of the speed of the inflaton (by about 14 %), and we show that it fits the temperature anisotropies and polarization data considerably better than a constant running (its chi^2 improves by ~ 7.5 over that of the constant running, at the cost of one more parameter). We also consider the possibility that the large scale temperature fluctuations are suppressed by an anti-correlation between tensor and scalar modes. Unfortunately, while such effect does affect the temperature fluctuations at large scales, it does not affect the temperature power spectrum and cannot, therefore, help in reconciling a large value of r with the limits from temperature fluctuations.

We present fully nonlinear consistency relations for the squeezed bispectrum of Large Scale Structure. These relations hold when the matter component of the Universe is composed of one or more species, and generalize those obtained in Peloso, Pietroni '13 and Kehagias, Riotto '13 in the single species case. The multi-species relations apply to the standard dark matter + baryons scenario, as well as to the case in which some of the fields are auxiliary quantities describing a particular population, such as dark matter halos or a specific galaxy class. If a large scale velocity bias exists between the different populations new terms appear in the consistency relations with respect to the single species case. As an illustration, we discuss two physical cases in which such a velocity bias can exist: (1) a new long range scalar force in the dark matter sector (resulting in a violation of the equivalence principle in the dark matter-baryon system), and (2) the distribution of dark matter halos relative to that of the underlying dark matter field.

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.

Soft limits of $N$-point correlation functions, in which one wavenumber is much smaller than the others, play a special role in constraining the physics of inflation. Anisotropic sources such as a vector field during inflation generate distinct angular dependence in all these correlators. In this paper we focus on the four-point correlator (the trispectrum $T$). We adopt a parametrization motivated by models in which the inflaton $\phi$ is coupled to a vector field through a $I^2 \left( \phi \right) F^2$ interaction, namely $T_{\zeta}({\bf k}_1, {\bf k}_2, {\bf k}_3, {\bf k}_4) \equiv \sum_n d_n [ P_n(\hat{\bf k}_1 \cdot \hat{\bf k}_3) + P_n(\hat{\bf k}_1 \cdot \hat{\bf k}_{12}) + P_n(\hat{\bf k}_3 \cdot \hat{\bf k}_{12}) ] P_{\zeta}(k_1) P_{\zeta}(k_3) P_\zeta(k_{12}) + (23~{\rm perm})$, where $P_n$ denotes the Legendre polynomials. This shape is enhanced when the wavenumbers of the diagonals of the quadrilateral are much smaller than the sides, ${\bf k}_i$. The coefficient of the isotropic part, $d_0$, is equal to $\tau_{\rm NL}/6$ discussed in the literature. A $I^2 \left( \phi \right) F^2$ interaction generates $d_2 = 2 d_0$ which is, in turn, related to the quadrupole modulation parameter of the power spectrum, $g_*$, as $d_2 \approx 14 |g_*| N^2$ with $N \approx 60$. We show that $d_0$ and $d_2$ can be equally well-constrained: the expected $68 \%$ CL error bars on these coefficients from a cosmic-variance-limited experiment measuring temperature anisotropy of the cosmic microwave background up to $\ell_{\rm max}=2000$ are $\delta d_2 \approx 4 \delta d_0 = 105$. Therefore, we can reach $|g_*|=10^{-3}$ by measuring the angle-dependent trispectrum. The current upper limit on $\tau_{\rm NL}$ from the ${\it Planck}$ temperature maps yields $|g_*|<0.02$ ($95 \%$ CL).

Gauge-flation is a recently proposed model in which inflation is driven solely by a non-Abelian gauge field thanks to a specific higher order derivative operator. The nature of the operator is such that it does not introduce ghosts. We compute the cosmological scalar and tensor perturbations for this model, improving over an existing computation. We then confront these results with the Planck data. The model is characterized by the quantity \gamma = (g^2 Q^2)/H^2 (where g is the gauge coupling constant, Q the vector vev, and H the Hubble rate). For \gamma < 2, the scalar perturbations show a strong tachyonic instability. In the stable region, the scalar power spectrum n_s is too low at small \gamma, while the tensor-to-scalar ratio r is too high at large \gamma. No value of \gamma leads to acceptable values for n_s and r, and so the model is ruled out by the CMB data. The same behavior with \gamma was obtained in Chromo-natural inflation, a model in which inflation is driven by a pseudo-scalar coupled to a non-Abelian gauge field. When the pseudo-scalar can be integrated out, one recovers the model of Gauge-flation plus corrections. It was shown that this identification is very accurate at the background level, but differences emerged in the literature concerning the perturbations of the two models. On the contrary, our results show that the analogy between the two models continues to be accurate also at the perturbative level.

In the model of solid / elastic inflation, inflation is driven by a source that has the field theoretical description of a solid. To allow for prolonged slow roll inflation, the solid needs to be extremely insensitive to the spatial expansion. We point out that, because of this property, the solid is also rather inefficient in erasing anisotropic deformations of the geometry. This allows for a prolonged inflationary anisotropic solution, providing the first example with standard gravity and scalar fields only which evades the conditions of the so called cosmic no-hair conjecture. We compute the curvature perturbations on the anisotropic solution, and the corresponding phenomenological bound on the anisotropy. Finally, we discuss the analogy between this model and the f (phi) F^2 model, which also allows for anisotropic inflation thanks to a suitable coupling between the inflaton phi and a vector field. We remark that the bispectrum of the curvature perturbations in solid inflation is enhanced in the squeezed limit and presents a nontrivial angular dependence, as had previously been found for the f (phi) F^2 model.

Because the inflationary mechanism is extremely sensitive to UV-physics, the construction of theoretically robust models of inflation provides a unique window on Planck-scale physics. We review efforts to use an axion with a shift symmetry to ensure a prolonged slow-roll background evolution. The symmetry dictates which operators are allowed, and these in turn determine the observational predictions of this class of models, which include observable gravitational waves (potentially chiral), oscillations in all primordial correlators, specific deviations from scale invariance and Gaussianity and primordial black holes. We discuss the constraints on this class of models in light of the recent Planck results and comment on future perspectives. The shift symmetry is very useful in models of large-field inflation, which typically have monomial potentials, but it cannot explain why two or more terms in the potential are fine-tuned against each other, as needed for typical models of small-field inflation. Therefore some additional symmetries or fine-tuning will be needed if forthcoming experiments will constrain the tensor-to-scalar ratio to be r < 0.01.

The bispectrum of primordial curvature perturbations in the squeezed configuration, in which one wavenumber, $k_3$, is much smaller than the other two, $k_3\ll k_1\approx k_2$, plays a special role in constraining the physics of inflation. In this paper we study a new phenomenological signature in the squeezed-limit bispectrum: namely, the amplitude of the squeezed-limit bispectrum depends on an angle between ${\bf k}_1$ and ${\bf k}_3$ such that $B_\zeta(k_1, k_2, k_3) \to 2 \sum_L c_L P_L(\hat{\bf k}_1 \cdot \hat{\bf k}_3) P_\zeta(k_1)P_\zeta(k_3)$, where $P_L$ are the Legendre polynomials. While $c_0$ is related to the usual local-form $f_{\rm NL}$ parameter as $c_0=6f_{\rm NL}/5$, the higher-multipole coefficients, $c_1$, $c_2$, etc., have not been constrained by the data. Primordial curvature perturbations sourced by large-scale magnetic fields generate non-vanishing $c_0$, $c_1$, and $c_2$. Inflation models whose action contains a term like $I(\phi)^2 F^2$ generate $c_2=c_0/2$. A recently proposed "solid inflation" model generates $c_2\gg c_0$. A cosmic-variance-limited experiment measuring temperature anisotropy of the cosmic microwave background up to $\ell_{\rm max}=2000$ is able to measure these coefficients down to $\delta c_0=4.4$, $\delta c_1=61$, and $\delta c_2=13$ (68% CL). We also find that $c_0$ and $c_1$, and $c_0$ and $c_2$, are nearly uncorrelated. Measurements of these coefficients will open up a new window into the physics of inflation such as the existence of vector fields during inflation or non-trivial symmetry structure of inflaton fields. Finally, we show that the original form of the Suyama-Yamaguchi inequality does not apply to the case involving higher-spin fields, but a generalized form does.

We discuss the constraints imposed on the nonlinear evolution of the Large Scale Structure (LSS) of the universe by galilean invariance, the symmetry relevant on subhorizon scales. Using Ward identities associated to the invariance, we derive fully nonlinear consistency relations between statistical correlators of the density and velocity perturbations, such as the power spectrum and the bispectrum. These relations are valid up to O (f_{NL}^2) corrections. We then show that most of the semi-analytic methods proposed so far to resum the perturbative expansion of the LSS dynamics fail to fulfill the constraints imposed by galilean invariance, and are therefore susceptible to non-physical infrared effects. Finally, we identify and discuss a nonperturbative semi-analytical scheme which is manifestly galilean invariant at any order of its expansion.

We perform the complete stability study of the model of chromo-natural inflation (Adshead and Wyman '12), where, due to its coupling to a SU(2) vector, a pseudo-scalar inflaton chi slowly rolls on a steep potential. As a typical example, one can consider an axion with a sub-Planckian decay constant f. The phenomenology of the model was recently studied (Dimastrogiovanni, Fasiello, and Tolley '12) in the m_g >> H limit, where m_g is the mass of the fluctuations of the vector field, and H the Hubble rate. We show that the inflationary solution is stable for m_g > 2 H, while it otherwise experiences a strong instability due to scalar perturbations in the sub-horizon regime. The tensor perturbations are instead standard, and the vector ones remain perturbatively small. Depending on the parameters, this model can give a gravity wave signal that can be detected in ongoing or forthcoming CMB experiments. This detection can occur even if, during inflation, the inflaton spans an interval of size Delta chi = O (f) which is some orders of magnitude below the Planck scale, evading a well known bound that holds for a free inflaton (Lyth '97).

A stochastic gravitational-wave background (SGWB) is expected to arise from the superposition of many independent and unresolved gravitational-wave signals, of either cosmological or astrophysical origin. Some cosmological models (characterized, for instance, by a pseudo-scalar inflaton, or by some modification of gravity) break parity, leading to a polarized SGWB. We present a new technique to measure this parity violation, which we then apply to the recent results from LIGO to produce the first upper limit on parity violation in the SGWB, assuming a generic power-law SGWB spectrum across the LIGO sensitive frequency region. We also estimate sensitivity to parity violation of the future generations of gravitational-wave detectors, both for a power-law spectrum and for a model of axion inflation. This technique offers a new way of differentiating between the cosmological and astrophysical sources of the isotropic SGWB, as astrophysical sources are not expected to produce a polarized SGWB.

A suitable coupling of the inflaton phi to a vector kinetic term F^2 gives frozen and scale invariant vector perturbations. We compute the cosmological perturbations zeta that result from such coupling by taking into account the classical vector field that unavoidably gets generated at large scales during inflation. This generically results in a too anisotropic power spectrum of zeta. Specifically, the anisotropy exceeds the 1% level (10% level) if inflation lasted ~5 e-folds (~50 e-folds) more than the minimal amount required to produce the CMB modes. This conclusion applies, among others, to the application of this mechanism for magnetogenesis, for anisotropic inflation, and for the generation of anisotropic perturbations at the end of inflation through a waterfall field coupled to the vector (in this case, the unavoidable contribution that we obtain is effective all throughout inflation, and it is independent of the waterfall field). For a tuned duration of inflation, a 1% (10%) anisotropy in the power spectrum corresponds to an anisotropic bispectrum which is enhanced like the local one in the squeezed limit, and with an effective local f_{NL} ~3 (~30). More in general, a significant anisotropy of the perturbations may be a natural outcome of all models that sustain higher than 0 spin fields during inflation.

We study the possibility that particle production during inflation could source observable gravity waves on scales relevant for Cosmic Microwave Background experiments. A crucial constraint on such scenarios arises because particle production can also source inflaton perturbations, and might ruin the usual predictions for a nearly scale invariant spectrum of nearly Gaussian curvature fluctuations. To minimize this effect, we consider two models of particle production in a sector that is only gravitationally coupled to the inflaton. For a single instantaneous burst of massive particle production, we find that localized features in the scalar spectrum and bispectrum might be observable, but gravitational wave signatures are unlikely to be detectable (due to the suppressed quadrupole moment of non-relativistic quanta) without invoking some additional effects. We also consider a model with a rolling pseudoscalar that leads to a continuous production of relativistic gauge field fluctuations during inflation. Here we find that gravitational waves from particle production can actually exceed the usual inflationary vacuum fluctuations in a regime where non-Gaussianity is consistent with observational limits. In this model observable B-mode polarization can be obtained for any choice of inflaton potential, and the amplitude of the signal is not necessarily correlated with the scale of inflation.

In any realistic particle physics model of inflation, the inflaton can be expected to couple to other fields. We consider a model with a dilaton-like coupling between a U(1) gauge field and a scalar inflaton. We show that this coupling can result in observable non-gaussianity, even in the conventional regime where inflation is supported by a single scalar slowly rolling on a smooth potential: the time dependent inflaton condensate leads to amplification of the large-scale gauge field fluctuations, which can feed-back into the scalar/tensor cosmological perturbations. In the squeezed limit, the resulting bispectrum is close to the local one, but it shows a sizable and characteristic quadrupolar dependence on the angle between the shorter and the larger modes in the correlation. Observable non-gaussianity is obtained in a regime where perturbation theory is under control. If the gauge field is identified with the electromagnetic field, the model that we study is a realization of the magnetogenesis idea originally proposed by Ratra, and widely studied. This identification (which is not necessary for the non-gaussianity production) is however problematic in light of a strong coupling problem already noted in the literature.

The reported spatial variation in the fine-structure constant at high redshift, if physical, could be due to the presence of dilatonic domains, and one or more domain walls inside our horizon. An absorption spectrum of an object in a different domain from our own would be characterized by a different value of alpha. We show that while a single wall solution is statically comparable to a dipole fit, and is a big improvement over a weighted mean (despite adding 3 parameters), a two-wall solution is a far better fit (despite adding 3 parameters over the single wall solution). We derive a simple model accounting for the two-domain wall solution. The goodness of these fits is however dependent on the extra random error which was argued to account for the large scatter in most of the data. When this error is omitted, all the above solutions are poor fits to the data. When included, the solutions that exhibit a spatial dependence agree with the data much more significantly than the Standard Model; however, the Standard Model itself is not a terrible fit to the data, having a p-value of ~ 20 %.

Models of inflation based on axions, which owe their popularity to the robustness against UV corrections, have also a very distinct class of signatures. The relevant interactions of the axion are a non-perturbative oscillating contribution to the potential and a shift-symmetric coupling to gauge fields. We review how these couplings affect the cosmological perturbations via a unified study based on the in-in formalism. We then note that, when the inflaton coupling to gauge fields is high enough to lead to interesting observational results, the backreaction of the produced gauge quanta on the inflaton dynamics becomes relevant during the final stage of inflation, and prolongs its duration by about 10 e-foldings. We extend existing results on gravity wave production in these models to account for this late inflationary phase. The strong backreaction phase results in an enhancement of the gravity wave signal at the interferometer scales. As a consequence, the signal is potentially observable at Advanced LIGO/VIRGO for the most natural duration of inflation in such models. Finally, we explicitly compute the axion couplings to gauge fields in string theory construction of axion monodromy inflation and identify cases where they can trigger interesting phenomenological effects.

Many controlled realizations of chaotic inflation employ pseudo-scalar axions. Pseudo-scalars \phi are naturally coupled to gauge fields through c \phi F \tilde{F}. In the presence of this coupling, gauge field quanta are copiously produced by the rolling inflaton. The produced gauge quanta, in turn, source inflaton fluctuations via inverse decay. These new cosmological perturbations add incoherently with the "vacuum" perturbations, and are highly nongaussian. This provides a natural mechanism to generate large nongaussianity in single or multi field slow-roll inflation. The resulting phenomenological signatures are highly distinctive: large nongaussianity of (nearly) equilateral shape, in addition to detectably large values of both the scalar spectral tilt and tensor-to-scalar ratio (both being typical of large field inflation). The WMAP bound on nongaussianity implies that the coupling, c, of the pseudo-scalar inflaton to any gauge field must be smaller than about 10^{2} M_p^{-1}.

The inflationary paradigm has enjoyed phenomenological success, however, a compelling particle physics realization is still lacking. The key obstruction is that the requirement of a suitably flat scalar potential is sensitive to Ultra-Violet (UV) physics. Axions are among the best-motivated inflaton candidates, since the flatness of their potential is naturally protected by a shift symmetry. We re-consider the cosmological perturbations in axion inflation, consistently accounting for the coupling to gauge fields \phi F \tilde{F}, which is generically present in these models. This coupling leads to production of gauge quanta, which provide a new source of inflaton fluctuations, \delta\phi. For an axion decay constant < 10^{-2} M_p, this effect typically dominates over the standard fluctuations from the vacuum, and saturates the current observational bounds on nongaussianity of the CMB anisotropies. Since sub-Planckian values of the decay constant are typical for concrete realizations that admit a UV completion (such as N-flation and axion monodromy), we conclude that large nongaussianity is easily obtained in very minimal and natural realizations of inflation.

We consider the signatures of a domain wall produced in the spontaneous symmetry breaking involving a dilaton-like scalar field coupled to electromagnetism. Domains on either side of the wall exhibit slight differences in their respective values of the fine-structure constant, alpha. If such a wall is present within our Hubble volume, absorption spectra at large redshifts may or may not provide a variation in alpha relative to the terrestrial value, depending on our relative position with respect to the wall. This wall could resolve the ``contradiction'' between claims of a variation of alpha based on Keck/Hires data and of the constancy of alpha based on VLT data. We derive the properties of the wall and the parameters of the underlying microscopic model required to reproduce the possible spatial variation of alpha. We discuss the constraints on the existence of the low-energy domain wall and describe its observational implications concerning the variation of the fundamental constants.

We consider gravitino couplings in theories with broken gauge symmetries. In particular, we compute the single gravitino production cross section in W+ W- fusion processes. Despite recent claims to the contrary, we show that this process is always subdominant to gluon fusion processes in the high energy limit. The full calculation is performed numerically; however, we give analytic expressions for the cross section in the supersymmetric and electroweak limits. We also confirm these results with the use of the effective theory of goldstino interactions.

We compute the phenomenological signatures of a model (Watanabe et al' 09) of anisotropic inflation driven by a scalar and a vector field. The action for the vector is U(1) invariant, and the model is free of ghost instabilities. A suitable coupling of the scalar to the kinetic term of the vector allows for a slow roll evolution of the vector vev, and hence for a prolonged anisotropic expansion; this provides a counter example to the cosmic no hair conjecture. We compute the nonvanishing two point correlation functions between physical modes of the system, and express them in terms of power spectra with angular dependence. The anisotropy parameter g_* for the scalar-scalar spectrum (defined as in the Ackerman et al '07 parametrization) turns out to be negative in the simplest realization of the model, which, therefore, cannot account for the angular dependence emerged in some analyses of the WMAP data. A g_* of order -0.1 is achieved when the energy of the vector is about 6-7 orders of magnitude smaller than that of the scalar during inflation. For such values of the parameters, the scalar-tensor correlation (which is in principle a distinctive signature of anisotropic spaces) is smaller than the tensor-tensor correlation.

We prove that many cosmological models characterized by vectors nonminimally coupled to the curvature (such as the Turner-Widrow mechanism for the production of magnetic fields during inflation, and models of vector inflation or vector curvaton) contain ghosts. The ghosts are associated with the longitudinal vector polarization present in these models, and are found from studying the sign of the eigenvalues of the kinetic matrix for the physical perturbations. Ghosts introduce two main problems: (1) they make the theories ill-defined at the quantum level in the high energy/sub horizon regime (and create serious problems for finding a well behaved UV completion); (2) they create an instability already at the linearized level. This happens because the eigenvalue corresponding to the ghost crosses zero during the cosmological evolution. At this point the linearized equations for the perturbations become singular (we show that this happens for all the models mentioned above). We explicitly solve the equations in the simplest cases of a vector without vev in a FRW geometry, and of a vector with vev plus a cosmological constant, and we show that indeed the solutions of the linearized equations diverge when these equations become singular.

Because the source term for the equations of motion of a conformally coupled scalar field, such as the dilaton, is given by the trace of the matter energy momentum tensor, it is commonly assumed to vanish during the radiation dominated epoch in the early universe. As a consequence, such fields are generally frozen in the early universe. Here we compute the finite temperature radiative correction to the source term and discuss its consequences on the evolution of such fields in the early universe. We discuss in particular, the case of scalar tensor theories of gravity which have general relativity as an attractor solution. We show that, in some cases, the universe can experience an early phase of contraction, followed by a non-singular bounce, and standard expansion. This can have interesting consequences for the abundance of thermal relics; for instance, it can provide a solution to the gravitino problem. We conclude by discussing the possible consequences of the quantum corrections to the evolution of the dilaton.

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

We prove that the anisotropic inflationary background of the Ackerman-Carroll-Wise model, characterized by a fixed-norm vector field, is unstable. We found the instability by explicitly solving the linearized equations for the most general set of perturbations around this background, and by noticing that the solutions diverge close to horizon crossing. This happens because one perturbation becomes a ghost at that moment. A simplified computation, with only the perturbations of the vector field included, shows the same instability, clarifying the origin of the problem. We then discuss several other models, with a particular emphasis on the case of a nonminimal coupling to the curvature, in which vector fields are used either to support an anisotropic expansion, or to generate cosmological perturbations on an isotropic background. In many cases, the mass term of the vector needs to have the ``wrong'' sign; we show that, as a consequence, the longitudinal vector mode is a ghost (a field with negative kinetic term, and negative energy; not simply a tachyon). We comment on problems that arise at the quantum level. In particular, the presence of a ghost can be a serious difficulty for the UV completion that such models require in the sub-horizon regime.

Models with vector fields acquiring a non-vanishing vacuum expectation value along one spatial direction have been proposed to sustain a prolonged stage of anisotropic accelerated expansion. Such models have been used for realizations of early time inflation, with a possible relation to the large scale CMB anomalies, or of the late time dark energy. We show that, quite generally, the concrete realizations proposed so far are plagued by instabilities (either ghosts, or unstable growth of the linearized perturbations) which can be ultimately related to the longitudinal vector polarization present in them. Phenomenological results based on these models are therefore unreliable.

We show that expanding or contracting Kasner universes are unstable due to the amplification of gravitational waves (GW). As an application of this general relativity effect, we consider a pre-inflationary anisotropic geometry characterized by a Kasner-like expansion, which is driven dynamically towards inflation by a scalar field. We investigate the evolution of linear metric fluctuations around this background, and calculate the amplification of the long-wavelength GW of a certain polarization during the anisotropic expansion (this effect is absent for another GW polarization, and for scalar fluctuations). These GW are superimposed to the usual tensor modes of quantum origin from inflation, and are potentially observable if the total number of inflationary e-folds exceeds the minimum required to homogenize the observable universe only by a small margin. Their contribution to the temperature anisotropy angular power spectrum decreases with the multipole l as l^(-p), where p depends on the slope of the initial GW power-spectrum. Constraints on the long-wavelength GW can be translated into limits on the total duration of inflation and the initial GW amplitude. The instability of classical GW (and zero-vacuum fluctuations of gravitons) during Kasner-like expansion (or contraction) may have other interesting applications. In particular, if GW become non-linear, they can significantly alter the geometry before the onset of inflation.

The presence of Kaluza-Klein particles in the universe is a potential manifestation of string theory cosmology. In general, they can be present in the high temperature bath of the early universe. In particular examples, string theory inflation often ends with brane-antibrane annihilation followed by the energy cascading through massive closed string loops to KK modes which then decay into lighter standard model particles. However, massive KK modes in the early universe may become dangerous cosmological relics if the inner manifold contains warped throat(s) with approximate isometries. In the complimentary picture, in the AdS/CFT dual gauge theory with extra symmetries, massive glueballs of various spins become the dangerous cosmological relics. The decay of these angular KK modes/glueballs, located around the tip of the throat, is caused by isometry breaking which results from gluing the throat to the compact CY manifold. We address the problem of these angular KK particles/glueballs, studying their interactions and decay channels, from the theory side, and the resulting cosmological constraints on the warped compactification parameters, from the phenomenology side. The abundance and decay time of the long-lived non-relativistic angular KK modes depend strongly on the parameters of the warped geometry, so that observational constraints rule out a significant fraction of the parameter space. In particular, the coupling of the angular KK particles can be weaker than gravitational.

We compute the nonperturbative decay of supersymmetric flat directions due to their D-term potential. Flat directions can develop large vacuum expectation values (vevs) during inflation, and, if they are long-lived, this can strongly affect the reheating and thermalization stages after the inflation. We study a generic system of two U(1) or SU(2) flat directions which are cosmologically evolving after inflation. After proper gauge fixing, we show that the excitations of the fields around this background can undergo exponential amplification, at the expense of the energy density of the flat directions. We compute this effect for several values of the masses and the initial vevs of the two flat directions, through a combination of analytical methods and extensive numerical simulations. For a wide range of parameters the flat directions decay within their first few rotations.

We extend the standard theory of cosmological perturbations to homogeneous but anisotropic universes. We present an exhaustive computation for the case of a Bianchi I model, with a residual isotropy between two spatial dimensions, which is undergoing complete isotropization at the onset of inflation; we also show how the computation can be further extended to more general backgrounds. In presence of a single inflaton field, there are three physical perturbations (precisely as in the isotropic case), which are obtained (i) by removing gauge and nondynamical degrees of freedom, and (ii) by finding the combinations of the remaining modes in terms of which the quadratic action of the perturbations is canonical. The three perturbations, which later in the isotropic regime become a scalar mode and two tensor polarizations (gravitational wave), are coupled to each other already at the linearized level during the anisotropic phase. This generates nonvanishing correlations between different modes of the CMB anisotropies, which can be particularly relevant at large scales (and, potentially, be related to the large scale anomalies in the WMAP data). As an example, we compute the spectrum of the perturbations in this Bianchi I geometry, assuming that the inflaton is in a slow roll regime also in the anisotropic phase. For this simple set-up, fixing the initial conditions for the perturbations appears more difficult than in the standard case, and additional assumptions seem to be needed to provide predictions for the CMB anisotropies.

Most of the analysis of the Cosmic Microwave Background relies on the assumption of statistical isotropy. However, given some recent evidence pointing against isotropy, as for instance the observed alignment of different multipoles on large scales, it is worth testing this assumption against the increasing amount of available data. As a pivot model, we assume that the spectrum of the primordial perturbations depends also on their directionality (rather than just on the magnitude of their momentum, as in the standard case). We explicitly compute the correlation matrix for the temperature anisotropies in the simpler case in which there is a residual isotropy between two spatial directions. As a concrete example, we consider a different initial expansion rate along one direction, and the following isotropization which takes place during inflation. Depending on the amount of inflation, this can lead to broken statistical isotropy on the largest observable scales.

We consider the role of supersymmetric flat directions in reheating the Universe after inflation. One or more flat directions can develop large vevs during inflation, which can potentially affect reheating by slowing down scattering processes among inflaton decay products or by coming to dominate the energy density of the Universe. Both effects occur only if flat directions are sufficiently long-lived. The computation of their perturbative decay rate, and a simple estimate of their nonperturbative decay have led to the conclusion that this is indeed the case. In contrast, we show that flat directions can decay quickly through nonperturbative channels in realistic models. The mass matrix for MSSM excitations around flat directions has nondiagonal entries, which vary with the phase of the (complex) flat directions. The quasi-periodic motion of the flat directions results in a strong parametric resonance, leading to the rapid depletion of the flat direction within its first few rotations. This may preclude any significant role for the flat directions in reheating the Universe after inflation in models in which the inflaton decays perturbatively.

We investigate the effects of bosonic trilinear interactions in preheating after chaotic inflation. A trilinear interaction term allows for the complete decay of the massive inflaton particles, which is necessary for the transition to radiation domination. We found that typically the trilinear term is subdominant during early stages of preheating, but it actually amplifies parametric resonance driven by the four-legs interaction. In cases where the trilinear term does dominate during preheating, the process occurs through periodic tachyonic amplifications with resonance effects, which is so effective that preheating completes within a few inflaton oscillations. We develop an analytic theory of this process, which we call tachyonic resonance. We also study numerically the influence of trilinear interactions on the dynamics after preheating. The trilinear term eventually comes to dominate after preheating, leading to faster rescattering and thermalization than could occur without it. Finally, we investigate the role of non-renormalizable interaction terms during preheating. We find that if they are present they generally dominate (while still in a controllable regime) in chaotic inflation models. Preheating due to these terms proceeds through a modified form of tachyonic resonance.

We study the out-of-equilibrium nonlinear dynamics of fields after post-inflationary preheating. During preheating, the energy in the homogeneous inflaton is exponentially rapidly transfered into highly occupied out-of-equilibrium inhomogeneous modes, which subsequently evolve towards equilibrium. The infrared modes excited during preheating evolve towards a saturated distribution long before thermalization completes. We compute the equation of state during and immediately after preheating. It rapidly evolves towards radiation domination long before the actual thermal equilibrium is established. The exact time of this transition is a non-monotonic function of the coupling between the inflaton and the decay products, and it varies only very weakly (around 10^(-35) s) as this coupling changes over several orders of magnitude. This result is applied to refine the relation between the number of efoldings N and the physical wavelength of perturbations generated during inflation. We also discuss the implications for the theory of modulated perturbations from preheating. We finally argue that many questions of the thermal history of the universe should be addressed in terms of pre-thermalization, illustrating this point with a calculation of perturbative production of gravitinos immediately after chaotic inflation. We also highlight the effects of three-legs inflaton interactions on the dynamics of preheating and thermalization in an expanding universe.

Any attempt to regularize a negative tension brane through a bulk scalar requires that this field is a ghost. One can try to improve in this aspect in a number of ways. For instance, it has been suggested to employ a field whose kinetic term is not sign definite, in the hope that the background may be overall stable. We show that this is not the case; the physical perturbations (gravity included) of the system do not extend across the zeros of the kinetic term; hence, all the modes are entirely localized either where the kinetic term is positive, or where it is negative; this second type of modes are ghosts. We show that this conclusion does not depend on the specific choice for the kinetic and potential functions for the bulk scalar.

We present a covariant nonlinear completion of the Fierz-Pauli (FP) mass term for the graviton. The starting observation is that the FP mass is immediately obtained by expanding the cosmological constant term, i.e. the determinant of the vielbein, around Minkowski space to second order in the vielbein perturbations. Since this is an unstable expansion in the standard case, we consider an extended theory of gravity which describes two vielbeins that give rise to chiral spin--connections (consequently, fermions of a definite chirality only couple to one of the gravitational sectors). As for Einstein gravity with a cosmological constant, a single fine-tuning is needed to recover a Minkowski background; the two sectors then differ only by a constant conformal factor. The spectrum of this theory consists of a massless and a massive graviton, with FP mass term. The theory possesses interesting limits in which only the massive graviton is coupled to matter at the linearized level.

If the inflaton is a pseudo-scalar axion, the axion shift symmetry can protect the flatness of its potential from too large radiative corrections. This possibility, known as natural inflation, requires an axion scale which is greater than the (reduced) Planck scale. It is unclear whether such a high value is compatible with an effective field theoretical description, and if the global axionic symmetry survives quantum gravity effects. We propose a mechanism which provides an effective large axion scale, although the original one is sub-Planckian. The mechanism is based on the presence of two axions, with a potential provided by two anomalous gauge groups. The effective large axion scale is due to an almost exact symmetry between the couplings of the axions to the anomalous groups. We also comment on a possible implementation in heterotic string theory.

We describe an algorithm and a C++ implementation that we have written and made available for calculating the fully nonlinear evolution of 5D braneworld models with scalar fields. Bulk fields allow for the stabilization of the extra space. However, they complicate the dynamics of the system, so that analytic calculations (performed within an effective 4D theory) are typically only reliable close to stabilized configurations or when the evolution of the extra space is negligible. In the general case, a numerical study of the 5D equations is necessary, and the algorithm and code we describe are the first ones designed for this task. The program and its full documentation are available on the Web at http://www.cita.utoronto.ca/~jmartin/BRANECODE/. In this paper we provide a brief overview of what the program does and how to use it.

Ghost condensation has been recently proposed as a mechanism inducing the spontaneous breaking of Lorentz symmetry. Corrections to the Newton potential generated by a static source have been computed: they yield a limit M < 10 MeV on the symmetry breaking scale, and - if the limit is saturated - they are maximal at a distance L ~ 1000 km from the source. However, these corrections propagate at a tiny velocity, v_s ~ 10^{-12} m/s, many orders of magnitude smaller than the velocity of any plausible source. We compute the gravitational potential taking the motion of the source into account: the standard Newton law is recovered in this case, with negligible corrections for any distance from the source up to astrophysical scales. Still, the vacuum of the theory is unstable, and requiring stability over the lifetime of the Universe gives a limit for M which is not too far from the one given above. In the absence of a direct coupling of the ghost to matter, signatures of this model will have to be searched in the form of exotic astrophysical events.

We consider warped compactifications in (4+d)-dimensional theories, with four dimensional de Sitter dS_4 vacua (with Hubble parameter H) and with a compact internal space. After introducing a gauge-invariant formalism for the generic metric perturbations of these backgrounds, we focus on modes which are scalar with respect to dS_4. The physical eigenmasses of these modes acquire a large universal tachyonic contribution -12d/(d+2) H^2, independently of the stabilization mechanism for the compact space, in addition to the usual KK masses, which instead encode the effects of the stabilization. General arguments, as well as specific examples, lead us to conjecture that, for sufficiently large dS curvature, the compactified geometry becomes gravitationally unstable due to the tachyonic growth of the scalar perturbations. This mean that for any stabilization mechanism the curvature of the dS geometry cannot exceed some critical value. We relate this effect to the anisotropy of the bulk geometry and suggest the end points of the instability. Of relevance for inflationary cosmology, the perturbations of the bulk metric inevitably induce a new modulus field, which describes the conformal fluctuations of the 4 dimensional metric. If this mode is light during inflation, the induced conformal fluctuations will be amplified with a scale free spectrum and with an amplitude which is disentangled from the standard result of slow-roll inflation. The conformal 4d metric fluctuations give rise to a very generic realization of the mechanism of modulated cosmological fluctuations, related to spatial variation of couplings during (p)reheating after inflation.

Braneworld models in extra dimensions can be tested in laboratory by the coupling of the radion to the Standard Model fields. The identification of the radion as a canonically normalized field involves a careful General Relativity treatment: if a bulk scalar is responsible for the stabilization of the system, its fluctuations are entangled with the perturbations of the metric and they also have to be taken into account (similarly to the well-developed theory of scalar metric perturbations in 4D cosmology with a scalar field). Extracting a proper dynamical variable in a warped geometry/scalar setting is a nontrivial task, performed so far only in the limit of negligible backreaction of the scalar field on the background geometry. We perform the general calculation, diagonalizing the action up to second order in the perturbations and identifying the physical eigenmodes of the system for any amplitude of the bulk scalar. This computation allows us to derive a very simple expression for the exact coupling of the eigenmodes to the Standard Model fields on the brane, valid for an arbitrary background configuration. As an application, we discuss the Goldberger-Wise mechanism for the stabilization of the radion in the Randall-Sundrum type models. The existing studies, limited to small amplitude of the bulk scalar field, are characterized by a radion mass which is significantly below the physical scale at the observable brane. We extend them beyond the small backreaction regime. For intermediate amplitudes, the radion mass approaches the electroweak scale, while its coupling to the observable brane remains nearly constant. At very high amplitudes, the radion mass instead decreases, while the coupling sharply increases. Severe experimental constraints are expected in this regime.

We give a full nonlinear numerical treatment of time-dependent 5d braneworld geometry, which is determined self-consistently by potentials for the scalar field in the bulk and at two orbifold branes, supplemented by boundary conditions at the branes. We describe the BraneCode, an algorithm which we designed to solve the dynamical equations numerically. We applied the BraneCode to braneworld models and found several novel phenomena of the brane dynamics. Starting with static warped geometry with de Sitter branes, we found numerically that this configuration is often unstable due to a tachyonic mass of the radion during inflation. If the model admits other static configurations with lower values of de Sitter curvature, this effect causes a violent re-structuring towards them, flattening the branes, which appears as a lowering of the 4d effective cosmological constant. Braneworld dynamics can often lead to brane collisions. We found that in the presence of the bulk scalar field, the 5d geometry between colliding branes approaches a universal, homogeneous, anisotropic strong gravity Kasner-like asymptotic, irrespective of the bulk/brane potentials. The Kasner indices of the brane directions are equal to each other but different from that of the extra dimension.

We show that shape moduli in sub-millimeter extra dimensional scenarios, addressing the gauge hierarchy problem, can dominate the energy density of the universe today. In our scenario, the volume of the extra dimensions is stabilized at a sufficiently high scale to avoid conflicts with nucleosynthesis and solar-system precision gravity experiments, while the shape moduli remain light but couple extremely weakly to brane-localized matter and easily avoid these bounds. Nonlocal effects in the bulk of the extra dimension generate a potential for the shape moduli. The potential has the right form and order of magnitude to account for the present day cosmic acceleration, in a way analogous to models of quintessence as a pseudo Nambu-Goldstone boson.

The Cosmic Microwave Background (CMB) anisotropy power on the largest angular scales observed both by WMAP and COBE DMR appears to be lower than the one predicted by the standard model of cosmology with almost scale free primordial perturbations arising from a period of inflation \cite{cobe,Bennett:2003bz,Spergel,Peiris}. One can either interpret this as a manifestation of cosmic variance or as a physical effect that requires an explanation. We discuss various mechanisms that could be responsible for the suppression of such low $\ell$ multipoles. Features in the late time evolution of metric fluctuations may do this via the integral Sachs-Wolfe effect. Another possibility is a suppression of power at large scales in the primordial spectrum induced by a fast rolling stage in the evolution of the inflaton field at the beginning of the last 65 e-folds of inflation. We illustrate this effect in a simple model of inflation and fit the resulting CMB spectrum to the observed temperature-temperature (TT) power spectrum. We find that the WMAP observations suggest a cutoff at $k_c=4.9^{+1.3}_{-1.6}\times 10^{-4}$Mpc$^{-1}$ at 68% confidence, while only an upper limit of $k_c < 7.4\times 10^{-4}$Mpc$^{-1}$ at 95%. Thus, although it improves the fit of the data, the presence of a cutoff in power spectrum is only required at a level close to $2\sigma$. This is obtained with a prior which corresponds to equal distribution wrt $k_c$. We discuss how other choices (such as an equal distribution wrt $\ln k_c$ which is natural in the context of inflation) can affect the statistical interpretation.

We show that photons with energies above the GZK cutoff can reach us from very distant sources if they mix with light axions in extragalactic magnetic fields. The effect which enables this is the conversion of photons into axions, which are sufficiently weakly coupled to travel large distances unimpeded. These axions then convert back into high energy photons close to the Earth. We show that photon-axion mixing facilitates the survival of super-GZK photons most efficiently with a photon-axion coupling scale of order 10^11 GeV, which is in the same range as the scale for the photon-axion mixing explanation for the dimming of supernovae without cosmic acceleration. We discuss possible observational consequences of this effect.

We consider a D3-brane as boundary of a five dimensional charged anti de Sitter black hole. We show that the charge of the black hole induces a regular cosmological evolution for the scale factor of the brane, with a smooth transition between a contracting and an eventual expanding phase. Simple analytical solutions can be obtained in the case of a vanishing effective cosmological constant on the brane. A nonvanishing cosmological constant, or the inclusion of radiation on the brane, does not spoil the regularity of these solutions at small radii, and observational constraints such as the ones from primordial nucleosynthesis can be easily met. Fluctuations of brane fields remain in the linear regime provided the minimal size of the scale factor is sufficiently large. We conclude with an analysis of the Cardy-Verlinde formula in this set up.

We investigate the possible effects on the evolution of perturbations in the inflationary epoch due to short distance physics. We introduce a suitable non local action for the inflaton field, suggested by Noncommutative Geometry, and obtained by adopting a generalized star product on a Friedmann-Robertson-Walker background. In particular, we study how the presence of a length scale where spacetime becomes noncommutative affects the gaussianity and isotropy properties of fluctuations, and the corresponding effects on the Cosmic Microwave Background spectrum.

The neutrino energy density of the Universe can be conveniently parametrized in terms of the so-called effective number of neutrinos, N_nu^eff. This parameter enters in several cosmological observables. In particular it is an important input in those numerical codes, like CMBFAST, which are used to study the Cosmic Microwave Background anisotropy spectrum. By studying the neutrino decoupling with Boltzmann equations, one can show that this quantity differs from the number of massless neutrino species for an additional contribution due to a partial heating of neutrinos during the electron-positron annihilations, leading to non thermal features in their final distributions. In this paper we review the different results obtained in the literature and perform a new analysis which takes into account, in a fully consistent way, the QED corrections at finite temperature to the photon and e^+- plasma equation of state. The value found for three massless active neutrinos is N_nu^eff=3.0395, in perfect agreement with the recommended value used in CMBFAST, N_nu^eff=3.04. We also discuss the case of additional relativistic relics and massive active neutrinos.

We study gravitational creation of light fermions in the presence of classical scalar metric perturbations about a flat Friedmann-Lemaitre- Robertson-Walker (FLRW) background. These perturbations can be large during preheating, breaking the conformal flatness of the background spacetime. We compute numerically the total number of particles generated by the modes of the metric perturbations which have grown sufficiently to become classical. In the absence of inhomogeneities massless fermions are not gravitationally produced, and then this effect may be relevant for abundance estimates of light gravitational relics.

We discuss the possibility that a perturbative reheating stage after inflation produces a scalar particle gas in a Bose condensate state, emphasizing the possible cosmological role of this phenomenon for symmetry restoration.

We consider the potential problems due to the production of inflatinos and gravitinos after inflation. Inflationary models with a single scale set by the microwave background anisotropies have a low enough reheat temperature to avoid problems with the thermal production of gravitinos. Moreover, the nonthermal production of gravitinos has been shown to be sufficiently small if the sector ultimately responsible for supersymmetry breaking is coupled only gravitationally to the inflationary sector. Still, in some models, inflatinos can be created during preheating with a substantial abundance. The main contribution to the gravitino abundance may thus come from their decay into the inflaton, or into its scalar partner, as well as from the inverse processes. We show that this production needs to be strongly suppressed. This suppression can be realized in the simplest scenarios which typically have a sufficiently high inflationary scale.

We provide the formalism for the quantization of systems of coupled bosonic and fermionic fields in a time dependent classical background. The occupation numbers of the particle eigenstates can be clearly defined and computed, through a generalization of the standard procedure valid for a single field in which Bogolyubov coefficients are employed. We apply our formalism to the problem of nonthermal gravitino production in a two-fields model where supersymmetry is broken gravitationally in the vacuum. Our explicit calculations show that this production is strongly suppressed in the model considered, due to the weak coupling between the sector which drives inflation and the one responsible for supersymmetry breakdown.

We explicitly calculate nonthermal gravitino production during the preheating period in the inflationary Universe. Contrary to earlier investigations, we consider a two--field model to separate the mechanisms of supersymmetry breaking and inflation. We show that the superpartner of the inflaton is significantly generated, while the gravitino production is considerably smaller. Nonthermal production of gravitinos seems thus less worrisome than recently claimed.

We calculate the exact kinetic evolution of cosmic neutrinos until complete decoupling, in the case when a large neutrino asymmetry exists. While not excluded by present observations, this large asymmetry can have relevant cosmological consequences and in particular may be helpful in reconciling Primordial Nucleosynthesis with a high baryon density as suggested by the most recent observations of the Cosmic Microwave Background Radiation. By solving numerically the Boltzmann kinetic equations for the neutrino distribution functions, we find the momentum-dependent corrections to the equilibrium spectra and briefly discuss their phenomenological implications.

The very recent Boomerang results give an estimate of unprecedented precision of the Cosmic Microwave Background anisotropies on sub--degree scales. A puzzling feature for theoretical cosmology is the low amplitude of the second acoustic peak. Through a qualitative discussion, we argue that a scarcely considered category of flat models, with a leptonic asymmetry, a high baryon density and a low cosmological constant seems to be in very good agreement with the data, while still being compatible with big bang nucleosynthesis and some other observational constraints. Although this is certainly not the only way to explain the data, we believe that these models deserve to be included in forthcoming likelihood analyses.

Non-perturbative production of fermions after chaotic inflation has been the object of several studies in the very recent past. However, the results in the most interesting case of production of massive fermions in an expanding Universe were so far known only numerically. We provide very simple and readable analytical formulae, both for the spectra of the created fermions and for their total energy density. Their derivation is closely related to the one adopted for bosons and exploits the fact that the production occurs during very short intervals of nonadiabatical change of the fermionic frequency. Our formulae show the presence of resonance bands if the expansion of the Universe is neglected, and their disappearance when the latter is included. As in the bosonic case, this last effect is due to the stochastic character that the expansion gives to the whole process. Backreaction is considered in the last part of the work. All our analytical results are in excellent agreement with the previous numerical ones in the regime of validity of the latter. However, a more accurate scaling for the energy density of the produced fermions is here found.

We provide the first complete analysis of cosmological evolution in the Randall-Sundrum model with stabilized dilaton. We give the exact expansion law for matter densities on the two branes with arbitrary equations of state. The effective four-dimensional theory leads to standard cosmology at low energy. The limit of validity of the low energy theory and possible deviations from the ordinary expansion law in the high energy regime are finally discussed.

We propose a new scenario of baryogenesis in the context of theories with large extra dimensions. The baryon number is almost conserved at zero temperature by means of a localization mechanism recently analyzed by Arkani-Hamed and Schmaltz: leptons and quarks are located at two slightly displaced positions in the extra space, and this naturally suppresses the interactions which ``convert'' the latter in the former. We show that this is expected to be no longer true when finite temperature effects are taken into account. The whole scenario is first presented in its generality, without referring to the bulk geometry or to the specific mechanism which may generate the baryon asymmetry. As an example, we then focus on a baryogenesis model reminiscent of GUT baryogenesis. The Sakharov out of equilibrium condition is satisfied by assuming nonthermal production of the bosons that induce baryon number violation.

Attention has been recently drawn towards models in which inflation and quintessence schemes are unified. In such `quintessential inflation' models, a unique scalar field is required to play both the role of the inflaton and of the late-time dynamical cosmological constant. We address the issue of the initial conditions for quintessence in this context and find that, in the two explicit examples provided, inflation can uniquely fix them to be in the allowed range for a present day tracking.

The standard CDM model fails to describe the power spectrum of fluctuations since it gives too much power at small scales. Among other possible improvements, it has been suggested that an agreement with observations can be achieved with the addition of a late decaying particle, through the injection of non-thermal radiation and the consequent increase of the horizon length at the equivalence time. We analyze the possibility of implementing this idea in some extensions of the electroweak standard model, discussing the cosmological and astrophysical bounds to which these schemes are subject.

Leptogenesis constitues a very simple scenario to achieve the baryon asymmetry that we observe today. It requires only the presence of right handed neutrinos (which arise very naturally in many extensions of the Standard Model) and depends crucially on the mechanism responsible for their production. In particular, when their mass exceeds the inflaton mass (10^13 GeV in chaotic inflation) only non perturbative production can occur. It is shown here that non perturbative production of fermions in an expanding Universe is a very efficient mechanism up to masses of order of (10^17 - 10^18) GeV, and that it can be easily applied to solve the baryon asymmetry problem.

We present a complete computation of the inflaton decay into very massive fermions during preheating. We show that heavy fermions are produced very efficiently up to masses of order 10^{17}-10^{18} GeV; the accessible mass range is thus even broader than the one for heavy bosons. We apply our findings to the leptogenesis scenario, proposing a new version of it, in which the massive right-handed neutrinos, responsible for the generation of the baryon asymmetry, are produced during preheating. We also discuss other production mechanisms of right-handed neutrinos in the early Universe, identifying the neutrino mass parameters for which the observed baryon asymmetry is reproduced.