The swampland conjectures seek to distinguish effective field theories which can be consistently embedded in a theory of quantum gravity from those which can not (and are hence referred to as being in the swampland). We consider two such conjectures, known as the Swampland Distance and de Sitter Conjectures, showing that taken together they place bounds on the amplitude of primordial gravitational waves generated during single field slow-roll inflation. The bounds depend on two parameters which for reasonable estimates restrict the tensor-to-scalar ratio to be within reach of future surveys.

A second order pole in the scalar kinetic term can lead to a class of inflation models with universal predictions referred to as pole inflation or $\alpha$-attractors. While this kinetic structure is ubiquitous in supergravity effective field theories, realising a consistent UV complete model in e.g. string theory is a non-trivial task. For one, one expects quantum corrections arising in the vicinity of the pole which may spoil the typical attractor dynamics. As a conservative estimate of the range of validity of supergravity models of pole inflation we employ the weak gravity conjecture (WGC). We find that this constrains the accessible part of the inflationary plateau by limiting the decay constant of the axion partner. For the original single complex field models, the WGC does not even allow the inflaton to reach the inflationary plateau region. We analyze if evoking the assistance of $N$ scalar fields from the open string moduli helps addressing these problems. Pole $N$-flation could improve radiative control by reducing the required range of each individual field. However, the WGC bound prohibiting pole inflation for a single such field persists even for a collective motion of $N$ such scalars if we impose the sublattice WGC. Finally, we outline steps towards an embedding of pole N-flation in type IIB string theory on fibred Calabi-Yau manifolds.

In this work we propose a statistical approach to handling sources of theoretical uncertainty in string theory models of inflation. By viewing a model of inflation as a probabilistic graph, we show that there is an inevitable information bottleneck between the ultraviolet input of the theory and observables, as a simple consequence of the data processing theorem. This information bottleneck can result in strong hierarchies in the sensitivity of observables to the parameters of the underlying model and hence universal predictions with respect to at least some microphysical considerations. We also find other intriguing behaviour, such as sharp transitions in the predictions when certain hyperparameters cross a critical value. We develop a robust numerical approach to studying these behaviours by adapting methods often seen in the context of machine learning. We first test our approach by applying it to well known examples of universality, sharp transitions, and concentration phenomena in random matrix theory. We then apply the method to inflation with axion monodromy. We find universality with respect to a number of model parameters and that consistency with observational constraints implies that with very high probability certain perturbative corrections are non-negligible.

We study inflation in models with many interacting fields subject to randomly generated scalar potentials. We use methods from non-equilibrium random matrix theory to construct the potentials and an adaption of the 'transport method' to evolve the two-point correlators during inflation. This construction allows, for the first time, for an explicit study of models with up to 100 interacting fields supporting a period of 'approximately saddle-point' inflation. We determine the statistical predictions for observables by generating over 30,000 models with 2-100 fields supporting at least 60 efolds of inflation. These studies lead us to seven lessons: i) Manyfield inflation is not single-field inflation, ii) The larger the number of fields, the simpler and sharper the predictions, iii) Planck compatibility is not rare, but future experiments may rule out this class of models, iv) The smoother the potentials, the sharper the predictions, v) Hyperparameters can transition from stiff to sloppy, vi) Despite tachyons, isocurvature can decay, vii) Eigenvalue repulsion drives the predictions. We conclude that many of the 'generic predictions' of single-field inflation can be emergent features of complex inflation models.

We study the sensitivity of cosmological observables to the reheating phase following inflation driven by many scalar fields. We describe a method which allows semi-analytic treatment of the impact of perturbative reheating on cosmological perturbations using the sudden decay approximation. Focusing on $\mathcal{N}$-quadratic inflation, we show how the scalar spectral index and tensor-to-scalar ratio are affected by the rates at which the scalar fields decay into radiation. We find that for certain choices of decay rates, reheating following multiple-field inflation can have a significant impact on the prediction of cosmological observables.

We present a complete framework for numerical calculation of the power spectrum and bispectrum in canonical inflation with an arbitrary number of light or heavy fields. Our method includes all relevant effects at tree-level in the loop expansion, including (i) interference between growing and decaying modes near horizon exit; (ii) correlation and coupling between species near horizon exit and on superhorizon scales; (iii) contributions from mass terms; and (iv) all contributions from coupling to gravity. We track the evolution of each correlation function from the vacuum state through horizon exit and the superhorizon regime, with no need to match quantum and classical parts of the calculation; when integrated, our approach corresponds exactly with the tree-level Schwinger or 'in-in' formulation of quantum field theory. In this paper we give the equations necessary to evolve all two- and three-point correlation functions together with suitable initial conditions. The final formalism is suitable to compute the amplitude, shape, and scale dependence of the bispectrum in models with |fNL| of order unity or less, which are a target for future galaxy surveys such as Euclid, DESI and LSST. As an illustration we apply our framework to a number of examples, obtaining quantitatively accurate predictions for their bispectra for the first time. Two accompanying reports describe publicly-available software packages that implement the method.

We construct ensembles of random scalar potentials for $N_f$ interacting scalar fields using non-equilibrium random matrix theory, and use these to study the generation of observables during small-field inflation. For $N_f={\cal O}({\rm few})$, these heavily featured scalar potentials give rise to power spectra that are highly non-linear, at odds with observations. For $N_f\gg 1$, the superhorizon evolution of the perturbations is generically substantial, yet the power spectra simplify considerably and become more predictive, with most realisations being well approximated by a linear power spectrum. This provides proof of principle that complex inflationary physics can give rise to simple emergent power spectra. We explain how these results can be understood in terms of large $N_f$ universality of random matrix theory.

We give an explicit relation, up to second-order terms, between scalar-field fluctuations defined on spatially-flat slices and the curvature perturbation on uniform-density slices. This expression is a necessary ingredient for calculating observable quantities at second-order and beyond in multiple-field inflation. We show that traditional cosmological perturbation theory and the `separate universe' approach yield equivalent expressions for superhorizon wavenumbers, and in particular that all nonlocal terms can be eliminated from the perturbation-theory expressions.

Even simple inflationary scenarios have many free parameters. Beyond the variables appearing in the inflationary action, these include dynamical initial conditions, the number of fields, and couplings to other sectors. These quantities are often ignored but cosmological observables can depend on the unknown parameters. We use Bayesian networks to account for a large set of inflationary parameters, deriving generative models for the primordial spectra that are conditioned on a hierarchical set of prior probabilities describing the initial conditions, reheating physics, and other free parameters. We use $N_f$--quadratic inflation as an illustrative example, finding that the number of $e$-folds $N_*$ between horizon exit for the pivot scale and the end of inflation is typically the most important parameter, even when the number of fields, their masses and initial conditions are unknown, along with possible conditional dependencies between these parameters.

We describe how to apply the transport method to compute inflationary observables in a broad range of multiple-field models. The method is efficient and encompasses scenarios with curved field-space metrics, violations of slow-roll conditions and turns of the trajectory in field space. It can be used for an arbitrary mass spectrum, including massive modes and models with quasi-single-field dynamics. In this note we focus on practical issues. It is accompanied by a Mathematica code which can be used to explore suitable models, or as a basis for further development.

Power suppression of the cosmic microwave background on the largest observable scales could provide valuable clues about the particle physics underlying inflation. Here we consider the prospect of power suppression in the context of the multifield landscape. Based on the assumption that our observable universe emerges from a tunnelling event and that the relevant features originate purely from inflationary dynamics, we find that the power spectrum not only contains information on single-field dynamics, but also places strong con- straints on all scalar fields present in the theory. We find that the simplest single-field models giving rise to power suppression do not generalise to multifield models in a straightforward way, as the resulting superhorizon evolution of the curvature perturbation tends to erase any power suppression present at horizon crossing. On the other hand, multifield effects do present a means of generating power suppression which to our knowledge has so far not been considered. We propose a mechanism to illustrate this, which we dub flume inflation.

We study the tensor spectral index $n_t$ and the tensor-to-scalar ratio $r$ in the simplest multifield extension to single-field, slow-roll inflation models. We show that multifield models with potentials $V \sim \sum_i \lambda_i |\phi_i|^p$ have different predictions for $n_t/r$ than single-field models, even when all the couplings are equal $\lambda_i=\lambda_j$, due to the probabilistic nature of the fields' initial values. We analyze well-motivated prior probabilities for the $\lambda_i$ and initial conditions to make detailed predictions for the marginalized probability distribution of $n_t/r$. With $\mathcal O(100)$ fields and $p>3/4$, we find that $n_t/r$ differs from the single-field result of $n_t/r=-1/8$ at the 5$\sigma$ level. This gives a novel and testable prediction for the simplest multifield inflation models.

This paper presents a method for obtaining an analytic expression for the density function of observables in multifield models of inflation with sum-separable potentials. The most striking result is that the density function in general possesses a sharp peak and the location of this peak is only mildly sensitive to the distribution of initial conditions. A simple argument is given for why this result holds for a more general class of models than just those with sum-separable potentials and why for such models, it is possible to obtain robust predictions for observable quantities. As an example, the joint density function of the spectral index and running in double quadratic inflation is computed. For scales leaving the horizon 55 e-folds before the end of inflation, the density function peaks at n_{s}=0.967 and \alpha=0.0006 for the spectral index and running respectively.

We present MultiModeCode, a Fortran 95/2000 package for the numerical exploration of multifield inflation models. This program facilitates efficient Monte Carlo sampling of prior probabilities for inflationary model parameters and initial conditions and is the first publicly available code that can efficiently generate large sample-sets for inflation models with $\mathcal O(100)$ fields. The code numerically solves the equations of motion for the background and first-order perturbations of multi-field inflation models with canonical kinetic terms and arbitrary potentials, providing the adiabatic, isocurvature, and tensor power spectra at the end of inflation. For models with sum-separable potentials MultiModeCode also computes the slow-roll prediction via the $\delta N$ formalism for easy model exploration and validation. We pay particular attention to the isocurvature perturbations as the system approaches the adiabatic limit, showing how to avoid numerical instabilities that affect some other approaches to this problem. We demonstrate the use of MultiModeCode by exploring a few toy models. Finally, we give a concise review of multifield perturbation theory and a user's manual for the program.

N-flation is a radiatively stable scenario for chaotic inflation in which the displacements of $N \gg 1$ axions with decay constants $f_1 \le \ldots \le f_N < M_P$ lead to a super-Planckian effective displacement equal to the Pythagorean sum $f_{Py}$ of the $f_i$. We show that mixing in the axion kinetic term generically leads to the phenomenon of kinetic alignment, allowing for effective displacements as large as $\sqrt{N} f_{N} \ge f_{Py}$, even if $f_1, \ldots, f_{N-1}$ are arbitrarily small. At the level of kinematics, the necessary alignment occurs with very high probability, because of eigenvector delocalization. We present conditions under which inflation can take place along an aligned direction. Our construction sharply reduces the challenge of realizing N-flation in string theory.

We explore whether multifield inflationary models make unambiguous predictions for fundamental cosmological observables. Focusing on $N$-quadratic inflation, we numerically evaluate the full perturbation equations for models with 2, 3, and $\mathcal{O}(100)$ fields, using several distinct methods for specifying the initial values of the background fields. All scenarios are highly predictive, with the probability distribution functions of the cosmological observables becoming more sharply peaked as $N$ increases. For $N=100$ fields, 95% of our Monte Carlo samples fall in the ranges $n_s \in (0.9455,0.9534)$; $\alpha \in (-9.741,-7.047)\times 10^{-4}$; $r\in(0.1445,0.1449)$; and $r_\mathrm{iso} \in (0.02137,3.510)\times 10^{-3}$ for the spectral index, running, tensor-to-scalar ratio, and isocurvature-to-adiabatic ratio, respectively. The expected amplitude of isocurvature perturbations grows with $N$, raising the possibility that many-field models may be sensitive to post-inflationary physics and suggesting new avenues for testing these scenarios.

We analyse the multifield behaviour in D-brane inflation when contributions from the bulk are taken into account. For this purpose, we study a large number of realisations of the potential; we find the nature of the inflationary trajectory to be very consistent despite the complex construction. Inflation is always canonical and occurs in the vicinity of an inflection point. Extending the transport method to non-slow-roll and to calculate the running, we obtain distributions for observables. The spectral index is typically blue and the running positive, putting the model under moderate pressure from WMAP7 constraints. The local f_NL and tensor-to-scalar ratio are typically unobservably small, though we find approximately 0.5% of realisations to give observably large local f_NL. Approximating the potential as sum-separable, we are able to give fully analytic explanations for the trends in observed behaviour. Finally we find the model suffers from the persistence of isocurvature perturbations, which can be expected to cause further evolution of adiabatic perturbations after inflation. We argue this is a typical problem for models of multifield inflation involving inflection points and renders models of this type technically unpredictive without a description of reheating.

A pressing problem in comparing inflationary models with observation is the accurate calculation of correlation functions. One approach is to evolve them using ordinary differential equations ("transport equations"), analogous to the Schwinger-Dyson hierarchy of in-out quantum field theory. We extend this approach to the complete set of momentum space correlation functions. A formal solution can be obtained using raytracing techniques adapted from geometrical optics. We reformulate inflationary perturbation theory in this language, and show that raytracing reproduces the familiar "delta N" Taylor expansion. Our method produces ordinary differential equations which allow the Taylor coefficients to be computed efficiently. We use raytracing methods to express the gauge transformation between field fluctuations and the curvature perturbation, zeta, in geometrical terms. Using these results we give a compact expression for the nonlinear gauge-transform part of fNL in terms of the principal curvatures of uniform energy-density hypersurfaces in field space.

We explore the super-horizon evolution of the two-point and three-point correlation functions of the primordial density perturbation in randomly-generated multi-field potentials. We use the Transport method to evolve perturbations and give full evolutionary histories for observables. Identifying the separate universe assumption as being analogous to a geometrical description of light rays, we give an expression for the width of the bundle, thereby allowing us to monitor evolution towards the adiabatic limit, as well as providing a useful means of understanding the behaviour in $f_NL$. Finally, viewing our random potential as a toy model of inflation in the string landscape, we build distributions for observables by evolving trajectories for a large number of realisations of the potential and comment on the prospects for testing such models. We find the distributions for observables to be insensitive to the number of fields over the range 2 to 6, but that these distributions are highly sensitive to the scale of features in the potential. Most sensitive to the scale of features is the spectral index, with more than an order of magnitude increase in the dispersion of predictions over the range of feature scales investigated. Least sensitive was the non-Gaussianity parameter $f_NL$, which was consistently small; we found no examples of realisations whose non-Gaussianity is capable of being observed by any planned experiment.

We explore inflationary trajectories within randomly-generated two-dimensional potentials, considered as a toy model of the string landscape. Both the background and perturbation equations are solved numerically, the latter using the two-field formalism of Peterson and Tegmark which fully incorporates the effect of isocurvature perturbations. Sufficient inflation is a rare event, occurring for only roughly one in $10^5$ potentials. For models generating sufficient inflation, we find that the majority of runs satisfy current constraints from WMAP. The scalar spectral index is less than 1 in all runs. The tensor-to-scalar ratio is below the current limit, while typically large enough to be detected by next-generation CMB experiments and perhaps also by Planck. In many cases the inflationary consistency equation is broken by the effect of isocurvature modes.

We explore the stability properties of multi-field solutions in the presence of a perfect fluid, as appropriate to assisted quintessence scenarios. We show that the stability condition for multiple fields $\phi_i$ in identical potentials $V_i$ is simply $d^2V_i/d \phi_i^2 > 0$, exactly as in the absence of a fluid. A possible new instability associated with the fluid is shown not to arise in situations of cosmological interest.