We study the stochastic formalism of inflation beyond the usual slow-roll approximation. We verify that the assumptions on which the stochastic formalism relies still hold even far from the slow-roll attractor. This includes demonstrating the validity of the separate universe approach to evolving long-wavelength scalar field perturbations beyond slow roll. We also explain that, in general, there is a gauge correction to the amplitude of the stochastic noise. This is because the amplitude is usually calculated in the spatially-flat gauge, while the number of e-folds is used as the time variable (hence one works in the uniform-$N$ gauge) in the Langevin equations. We show that these corrections vanish in the slow-roll limit, but we also explain how to calculate them in general. We compute them in difference cases, including ultra-slow roll and the Starobinsky model that interpolates between slow roll and ultra-slow roll, and find the corrections to be negligible in practice. This confirms the validity of the stochastic formalism for studying quantum backreaction effects in the very early universe beyond slow roll.

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

The relative probability to decay towards different vacua during inflation is studied. The calculation is performed in single-field slow-roll potentials using the stochastic inflation formalism. Various situations are investigated, including falling from a local maximum of the potential and escaping from a local minimum. In the latter case, our result is consistent with that of Hawking and Moss, but is applicable to any potential. The decay rates are also computed, and the case of a generic potential with multiple minima and maxima is discussed.

It is often claimed that the ultra-slow-roll regime of inflation, where the dynamics of the inflaton field are friction dominated, is a non-attractor and/or transient. In this work we carry out a phase-space analysis of ultra-slow roll in an arbitrary potential, $V(\phi)$. We show that while standard slow roll is always a dynamical attractor whenever it is a self-consistent approximation, ultra-slow roll is stable for an inflaton field rolling down a convex potential with $M_{\scriptscriptstyle{\mathrm{Pl}}} V''>|V'|$ (or for a field rolling up a concave potential with $M_{\scriptscriptstyle{\mathrm{Pl}}} V''<-|V'|$). In particular, when approaching a flat inflection point, ultra-slow roll is always stable and a large number of $e$-folds may be realised in this regime. However, in ultra-slow roll, $\dot{\phi}$ is not a unique function of $\phi$ as it is in slow roll and dependence on initial conditions is retained. We confirm our analytical results with numerical examples.

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

We calculate the full probability density function (PDF) of inflationary curvature perturbations, even in the presence of large quantum backreaction. Making use of the stochastic-$\delta N$ formalism, two complementary methods are developed, one based on solving an ordinary differential equation for the characteristic function of the PDF, and the other based on solving a heat equation for the PDF directly. In the classical limit where quantum diffusion is small, we develop an expansion scheme that not only recovers the standard Gaussian PDF at leading order, but also allows us to calculate the first non-Gaussian corrections to the usual result. In the opposite limit where quantum diffusion is large, we find that the PDF is given by an elliptic theta function, which is fully characterised by the ratio between the squared width and height (in Planck mass units) of the region where stochastic effects dominate. We then apply these results to the calculation of the mass fraction of primordial black holes from inflation, and show that no more than $\sim 1$ $e$-fold can be spent in regions of the potential dominated by quantum diffusion. We explain how this requirement constrains inflationary potentials with two examples.

Stochastic effects in generic scenarios of inflation with multiple fields are investigated. First passage time techniques are employed to calculate the statistical moments of the number of inflationary $e$-folds, which give rise to all correlation functions of primordial curvature perturbations through the stochastic $\delta N$ formalism. The number of fields is a critical parameter. The probability of exploring arbitrarily large-field regions of the potential becomes non-vanishing when more than two fields are driving inflation. The mean number of $e$-folds can be infinite, depending on the number of fields; for plateau potentials, this occurs even with one field. In such cases, correlation functions of curvature perturbations are infinite. They can, however, be regularised if a reflecting (or absorbing) wall is added at large energy or field value. The results are found to be independent of the exact location of the wall and this procedure is, therefore, well-defined for a wide range of cutoffs, above or below the Planck scale. Finally, we show that, contrary to single-field setups, multi-field models can yield large stochastic corrections even at sub-Planckian energy, opening interesting prospects for probing quantum effects on cosmological fluctuations.

Stochastic effects in multi-field inflationary scenarios are investigated. A hierarchy of diffusion equations is derived, the solutions of which yield moments of the numbers of inflationary $e$-folds. Solving the resulting partial differential equations in multi-dimensional field space is more challenging than the single-field case. A few tractable examples are discussed, which show that the number of fields is, in general, a critical parameter. When more than two fields are present for instance, the probability to explore arbitrarily large-field regions of the potential, otherwise inaccessible to single-field dynamics, becomes non-zero. In some configurations, this gives rise to an infinite mean number of $e$-folds, regardless of the initial conditions. Another difference with respect to single-field scenarios is that multi-field stochastic effects can be large even at sub-Planckian energy. This opens interesting new possibilities for probing quantum effects in inflationary dynamics, since the moments of the numbers of $e$-folds can be used to calculate the distribution of primordial density perturbations in the stochastic-$\delta N$ formalism.

We explore the hemispherical asymmetry predicted in cosmic microwave background polarization when there is an asymmetry in temperature anisotropies due to primordial perturbations. We consider the cases of asymmetries due to adiabatic and isocurvature modes, and tensor perturbations. We show that the asymmetry in the TE, EE and/or BB correlations can be substantially larger than those in the TT power spectrum in certain cases. The relative asymmetry in the different cross-correlations, as well as the angular scale dependence, can in principle distinguish between different origins for the asymmetry.

We investigate whether non-adiabatic perturbations from inflation could produce an asymmetric distribution of temperature anisotropies on large angular scales in the cosmic microwave background (CMB). We use a generalised non-linear $\delta N$ formalism to calculate the non-Gaussianity of the primordial density and isocurvature perturbations due to the presence of non-adiabatic, but approximately scale-invariant field fluctuations during multi-field inflation. This local-type non-Gaussianity leads to a correlation between very long wavelength inhomogeneities, larger than our observable horizon, and smaller scale fluctuations in the radiation and matter density. Matter isocurvature perturbations contribute primarily to low CMB multipoles and hence can lead to a hemispherical asymmetry on large angular scales, with negligible asymmetry on smaller scales. In curvaton models, where the matter isocurvature perturbation is partly correlated with the primordial density perturbation, we are unable to obtain a significant asymmetry on large angular scales while respecting current observational constraints on the observed quadrupole. However in the axion model, where the matter isocurvature and primordial density perturbations are uncorrelated, we find it may be possible to obtain a significant asymmetry due to isocurvature modes on large angular scales. Such an isocurvature origin for the hemispherical asymmetry would naturally give rise to a distinctive asymmetry in the CMB polarisation on large scales.

We study primordial density perturbations generated by the late decay of a curvaton field whose decay rate may be modulated by the local value of another isocurvature field, analogous to models of modulated reheating at the end of inflation. We calculate the primordial density perturbation and its local-type non-Gaussianity using the sudden-decay approximation for the curvaton field, recovering standard curvaton and modulated reheating results as limiting cases. We verify the Suyama-Yamaguchi inequality between bispectrum and trispectrum parameters for the primordial density field generated by multiple field fluctuations, and find conditions for the bound to be saturated.

We study the primordial density perturbations and non-Gaussianities generated from the combined effects of an inhomogeneous end of inflation and curvaton decay in hybrid inflation. This dual role is played by a single isocurvature field which is massless during inflation but acquire a mass at the end of inflation via the waterfall phase transition. We calculate the resulting primordial non-Gaussianity characterized by the non-linearity parameter, $f_{NL}$, recovering the usual end-of-inflation result when the field decays promptly and the usual curvaton result if the field decays sufficiently late.

We consider the stochastic background of gravitational waves produced during the radiation-dominated hot big bang as a constraint on the primordial density perturbation on comoving length scales much smaller than those directly probed by the cosmic microwave background or large-scale structure. We place weak upper bounds on the primordial density perturbation from current data. Future detectors such as BBO and DECIGO will place much stronger constraints on the primordial density perturbation on small scales.

We compute the second-order matching conditions for tensor metric perturbations at an abrupt change in the equation of state. For adiabatic perturbations on large scales the matching hypersurface coincides with a uniform-density hypersurface. We show that in the uniform-density gauge both the tensor perturbation and its time-derivative are continuous in this case. For non-adiabatic perturbations, the matching hypersurface need not coincide with a uniform-density hypersurface and the tensor perturbation in the uniform-density gauge may be discontinuous. However, we show that in the Poisson gauge both the tensor perturbation and its time-derivative are continuous for adiabatic or non-adiabatic perturbations. As an application we solve the evolution equation for second-order tensor perturbations on large scales for a constant equation of state and we use the matching conditions to evolve the solutions through the transition from an inflationary era to a radiation era. We show that in the radiation era the resulting free part of the large-scale tensor perturbation (constant mode) is slow-roll suppressed in both the uniform-density and Poisson gauges. Thus, we conclude that second-order gravitational waves from slow-roll inflation are suppressed.

We investigate the generation of gravitational waves due to the gravitational instability of primordial density perturbations in an early matter-dominated era which could be detectable by experiments such as LIGO and LISA. We use relativistic perturbation theory to give analytic estimates of the tensor perturbations generated at second order by linear density perturbations. We find that large enhancement factors with respect to the naive second-order estimate are possible due to the growth of density perturbations on sub-Hubble scales. However very large enhancement factors coincide with a breakdown of linear theory for density perturbations on small scales. To produce a primordial gravitational wave background that would be detectable with LIGO or LISA from density perturbations in the linear regime requires primordial comoving curvature perturbations on small scales of order 0.02 for Advanced LIGO or 0.005 for LISA, otherwise numerical calculations of the non-linear evolution on sub-Hubble scales are required.

We study a model where two scalar fields, that are subdominant during inflation, decay into radiation some time after inflation has ended but before primordial nucleosynthesis. Perturbations of these two curvaton fields can be responsible for the primordial curvature perturbation. We write down the full non-linear equations that relate the primordial perturbation to the curvaton perturbations on large scales, and solve them in a sudden-decay approximation. We calculate the power spectrum of the primordial perturbation, and finally go to second order to find the non-linearity parameter, fNL. Not surprisingly, we find large positive values of fNL if the energy densities of the curvatons are sub-dominant when they decay, as in the single curvaton case. But we also find a novel effect, which can be present only in multi-curvaton models: fNL becomes large even if the curvatons dominate the total energy density in the case when the inhomogeneous radiation produced by the first curvaton decay is diluted by the decay of a second nearly homogeneous curvaton. The minimum value min(fNL)=-5/4 which we find is the same as in the single-curvaton case. Using (non-)Gaussianity observations, Planck can be able to distinguish between single-field inflation and curvaton model. Hence it is important to derive theoretical predictions for curvaton model. From particle physics point of view it is more natural to assume multiple scalar fields (rather than just one ``curvaton'' in addition to inflaton). Our work updates the theoretical predictions of curvaton model to this case.

We study a model where two scalar fields, that are subdominant during inflation, decay into radiation some time after inflation has ended but before primordial nucleosynthesis. Perturbations of these two curvaton fields can be responsible for the primordial curvature perturbation. We write down the full non-linear equations that relate the primordial perturbation to the curvaton perturbations on large scales, calculate the power spectrum of the primordial perturbation, and finally go to second order to find the non-linearity parameter, fNL. We find large positive values of fNL if the energy densities of the curvatons are sub-dominant when they decay, as in the single curvaton case. But we also find a large fNL even if the curvatons dominate the total energy density in the case when the inhomogeneous radiation produced by the first curvaton decay is diluted by the decay of a second nearly homogeneous curvaton. The minimum value min(fNL)=-5/4 which we find is the same as in the single-curvaton case.