In this paper we study the observational constraints that can be imposed on the coupling parameter, $\hat \alpha$, of the recently proposed 4-dimensional Einstein-Gauss-Bonnet theory of gravity. We consider the field equations that can be derived from regularization or dimensional reduction, and use them to perform a thorough investigation of the slow-motion and weak-field limit of this theory. We apply our results to observations of a wide array of physical systems that admit such a description, and find that the LAGEOS satellites are the most constraining, requiring $| \hat \alpha | \lesssim 10^{10} \,{\rm m}^2$. We then consider constraints that can be imposed on this theory from cosmology, black hole systems, and table-top experiments. It is found that early universe inflation prohibits all but the smallest negative values of $\hat \alpha$, while observations of binary black hole systems are likely to offer the tightest constraints on positive values, leading to overall bounds $0 \lesssim \hat \alpha \lesssim 10^8 \, {\rm m}^2$.

We consider the steepest rate at which the power spectrum from single field inflation can grow, with the aim of providing a simple explanation for the $k^4$ growth found recently. With this explanation in hand we show that a slightly steeper $k^5 (\log k )^2$ growth is in fact possible. Moreover, we argue that the power spectrum after a steep growth cannot immediately decay, but must remain large for the $k$ modes which exit during a $\sim2$ e-fold period. We also briefly consider how a strong growth can affect the spectral index of longer wavelengths preceding the growth, and show that even the conversion of isocurvature modes likely cannot lead to a stronger growth. These results have implications for the formation of primordial black holes, and other phenomena which require a large amplitude of power spectrum at short scales.

We study the contribution to the primordial curvature perturbation on observational scales generated by the reheating field in massless preheating. To do so we use lattice simulations and a recent extension to the $\delta N$ formalism. The work demonstrates the functionality of these techniques for calculating the observational signatures of models in which non-perturbative reheating involves a light scalar field.

We revisit the question of how to calculate correlations of the curvature perturbation, $\zeta$, using the $\delta N$ formalism when one cannot employ a truncated Taylor expansion of $N$. This problem arises when one uses lattice simulations to probe the effects of isocurvature modes on models of reheating. Working in real space, we use an expansion in the cross-correlation between fields at different positions, and present simple expressions for observables such as the power spectrum and the reduced bispectrum, $f_{\rm NL}$. These take the same form as those of the usual $\delta N$ expressions, but with the derivatives of $N$ replaced by non-perturbative $\delta N$ coefficients. We test the validity of this expansion and argue that our expressions are particularly well suited for use with simulations.

We study multifield inflation in scenarios where the fields are coupled non-minimally to gravity via $\xi_I(\phi^I)^n g^{\mu\nu}R_{\mu\nu}$, where $\xi_I$ are coupling constants, $\phi^I$ the fields driving inflation, $g_{\mu\nu}$ the space-time metric, $R_{\mu\nu}$ the Ricci tensor, and $n>0$. We consider the so-called $\alpha$-attractor models in two formulations of gravity: in the usual metric case where $R_{\mu\nu}=R_{\mu\nu}(g_{\mu\nu})$, and in the Palatini formulation where $R_{\mu\nu}$ is an independent variable. As the main result, we show that, regardless of the underlying theory of gravity, the field-space curvature in the Einstein frame has no influence on the inflationary dynamics at the limit of large $\xi_I$, and one effectively retains the single-field case. However, the gravity formulation does play an important role: in the metric case the result means that multifield models approach the single-field $\alpha$-attractor limit, whereas in the Palatini case the attractor behaviour is lost also in the case of multifield inflation. We discuss what this means for distinguishing between different models of inflation.

We extend the transport framework for numerically evaluating the power spectrum and bispectrum in multi-field inflation to the case of a curved field-space metric. This method naturally accounts for all sub- and super-horizon tree level effects, including those induced by the curvature of the field-space. We present an open source implementation of our equations in an extension of the publicly available PyTransport code. Finally we illustrate how our technique is applied to examples of inflationary models with a non-trivial field-space metric.

PyTransport constitutes a straightforward code written in C++ together with Python scripts which automatically edit, compile and run the C++ code as a Python module. It has been written for Unix-like systems (OS X and Linux). Primarily the module employs the transport approach to inflationary cosmology to calculate the tree-level power-spectrum and bispectrum of user specified models of multi-field inflation, accounting for all sub and super-horizon effects. The transport method we utilise means only coupled differential equations need to be solved, and the implementation presented here combines the speed of C++ with the functionality and convenience of Python. This document details the code and illustrates how to use it with a worked example. It has been updated to be a companion to the second version of the code, PyTransport 2.0, which includes functionality to deal with models of inflation with a curved field space metric.

We study the behaviour of linear perturbations in multifield coupled quintessence models. Using gauge invariant linear cosmological perturbation theory we provide the full set of governing equations for this class of models, and solve the system numerically. We apply the numerical code to generate growth functions for various examples, and compare these both to the standard $\Lambda$CDM model and to current and future observational bounds. Finally, we examine the applicability of the "small scale approximation", often used to calculate growth functions in quintessence models, in light of upcoming experiments such as SKA and Euclid. We find the deviation of the full equation results for large k modes from the approximation exceeds the experimental uncertainty for these future surveys. The numerical code, PYESSENCE, written in Python will be publicly available.

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 develop a formalism for calculating soft limits of $n$-point inflationary correlation functions using separate universe techniques. Our method naturally allows for multiple fields and leads to an elegant diagrammatic approach. As an application we focus on the trispectrum produced by inflation with multiple light fields, giving explicit formulae for all possible single- and double-soft limits. We also investigate consistency relations and present an infinite tower of inequalities between soft correlation functions which generalise the Suyama-Yamaguchi inequality.

In this work, we present a method for implementing the $\delta N$ formalism to study the primordial non-Gaussianity produced in multiple three-form field inflation. Using a dual description relating three-form fields to noncanonical scalar fields, and employing existing results, we produce expressions for the bispectrum of the curvature perturbation in terms of three-form quantities. We study the bispectrum generated in a two three-form field inflationary scenario for a particular potential that for suitable values of the parameters was found in earlier work to give values of the spectral index and ratio of tensor to scalar perturbations compatible with current bounds. We calculate the reduced bispectrum for this model, finding an amplitude in equilateral and orthogonal configurations of ${\cal O}(1)$ and in the squeezed limit of ${\cal O}(10^{-3})$. We confirm, therefore, that this three-form inflationary scenario is compatible with present observational constraints.

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.

We calculate the squeezed limit of the bispectrum produced by inflation with multiple light fields. To achieve this we allow for different horizon exit times for each mode and calculate the intrinsic field-space three-point function in the squeezed limit using soft-limit techniques. We then use the $\delta N$ formalism from the time the last mode exits the horizon to calculate the bispectrum of the primordial curvature perturbation. We apply our results to calculate the spectral index of the halo bias, $n_{\delta b}$, an important observational probe of the squeezed limit of the primordial bispectrum and compare our results with previous formulae. We give an example of a curvaton model with $n_{\delta b} \sim {\cal O}(n_s-1)$ for which we find a 20% correction to observable parameters for squeezings relevant to future experiments. For completeness, we also calculate the squeezed limit of three-point correlation functions involving gravitons for multiple field models.

We consider a higher order term in the $\delta N$ expansion for the CMB power asymmetry generated by a superhorizon isocurvature field fluctuation. The term can generate the asymmetry without requiring a large value of $f_{NL}$. Instead it produces a non-zero value of $g_{NL}$. A combination of constraints lead to an allowed region in $f_{NL}-g_{NL}$ space. To produce the asymmetry with this term without a large value of $f_{NL}$ we find that the isocurvature field needs to contribute less than the inflaton towards the power spectrum of the curvature perturbation.

We discuss correlations among spectral observables as a new tool for differentiating between models for the primordial perturbation. We show that if generated in the isocurvature sector, a running of the scalar spectral index is correlated with the statistical properties of non-Gaussianities. In particular, we find a large running will inevitably be accompanied by a large running of $f_{\rm NL}$ and enhanced $g_{\rm NL}$, with $g_{\rm NL}\gg f_{\rm NL}^2$. If the tensor to scalar ratio is large, a large negative running must turn positive on smaller scales. Interestingly, the characteristic scale of the transition could potentially distinguish between the inflaton and isocurvature fields.

The observational signatures of multi-field inflation will generally evolve as the Universe reheats. We introduce a general analytic formalism for tracking this evolution through perturbative reheating, applicable to two field models with arbitrary separable potentials. The various transitions, including the onset of scalar field oscillations and the reheating of each field, can happen in different orders and on arbitrary hypersurfaces. The effective equations of state of the oscillating fields are also arbitrary. Nevertheless, our results are surprisingly simple. Our formalism encapsulates and generalises a huge range of previous calculations including two-field inflation, spectator models, the inhomogeneous end of inflation scenario and numerous generalised curvaton scenarios.

We explore inflation via the effective potential of the minimal Wess-Zumino model, considering both the real and imaginary components of the complex field. Using transport techniques, we calculate the full allowed range of $n_s$, $r$ and $f_{\rm NL}$ for different choices of the single free parameter, $v$, and present the probability distribution of these signatures given a simple choice for the prior distribution of initial conditions. Our work provides a case study of multi-field inflation in a simple but realistic setting, with important lessons that are likely to apply more generally. For example, we find that there are initial conditions consistent with observations of $n_s$ and $r$ for values of $v$ that would be excluded if only evolutions in the real field direction were to be considered, and that these may yield enhanced values of $f_{\rm NL}$. Moreover, we find that initial conditions fixed at high energy density, where the potential is close to quartic in form, can still lead to evolutions in a concave region of the potential during the observable number of e-folds, as preferred by present data. The Wess-Zumino model therefore provides an illustration that multi-field dynamics must be taken into account when seeking to understand fully the phenomenology of such models of inflation.

Planck data has not found the 'smoking gun' of non-Gaussianity that would have necessitated consideration of inflationary models beyond the simplest canonical single field scenarios. This raises the important question of what these results do imply for more general models, and in particular, multi-field inflation. In this paper we revisit four ways in which two-field scenarios can behave differently from single field models; two-field slow-roll dynamics, curvaton-type behaviour, inflation ending on an inhomogeneous hypersurface and modulated reheating. We study the constraints that Planck data puts on these classes of behaviour, focusing on the latter two which have been least studied in the recent literature. We show that these latter classes are almost equivalent, and extend their previous analyses by accounting for arbitrary evolution of the isocurvature mode which, in particular, places important limits on the Gaussian curvature of the reheating hypersurface. In general, however, we find that Planck bispectrum results only constrain certain regions of parameter space, leading us to conclude that inflation sourced by more than one scalar field remains an important possibility.

This proceedings contribution provides a brief update on the transport approach to calculating the statistics of perturbations produced during inflation. It is based in particular on work in collaboration with Gemma Anderson and David Seery.

We extend the `moment transport method' for calculating the statistics of inflationary perturbations to the quantum phase of evolution on sub-horizon scales. The quantum transport equations form a set of coupled ordinary differential equations for the evolution of quantum correlation functions during inflation, which are valid on sub- and super-horizon scales, and reduce to the known classical transport equations after horizon crossing. The classical and quantum equations follow directly from the field equations of cosmological perturbation theory. In this paper, we focus on how the evolution equations arise, and explore how transport methods relate to other approaches, and in particular how formal integral solutions to the transport equations connect to those of the In-In formalism.

We calculate the perturbed action, at second and third order, for a massive three-form field minimally coupled to gravity, and use it to explore the observational predictions of three-form inflation. One intriguing result is that the value of the spectral index is nearly independent of the three-form potential, being fixed solely by the number of e-folds of inflation, with n_s=0.97 for the canonical number of 60. Considering the bispectrum, we employ standard techniques to give explicit results for two models, one of which produces a large non-Gaussianity. Finally, we confirm our results by employing a duality relating the three-form theory to a non-canonical scalar field theory and explicitly re-computing results in this dual picture.

We use transport techniques to calculate the trispectrum produced in multiple-field inflationary models with canonical kinetic terms. Our method allows the time evolution of the local trispectrum parameters, tauNL and gNL, to be tracked throughout the inflationary phase. We illustrate our approach using examples. We give a simplified method to calculate the superhorizon part of the relation between field fluctuations on spatially flat hypersurfaces and the curvature perturbation on uniform density slices, and obtain its third-order part for the first time. We clarify how the 'backwards' formalism of Yokoyama et al. relates to our analysis and other recent work. We supply explicit formulae which enable each inflationary observable to be computed in any canonical model of interest, using a suitable first-order ODE solver.

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 study inflationary perturbations in multiple-field models, for which zeta typically evolves until all isocurvature modes decay--the "adiabatic limit". We use numerical methods to explore the sensitivity of the nonlinear parameter fNL to the process by which this limit is achieved, finding an appreciable dependence on model-specific data such as the time at which slow-roll breaks down or the timescale of reheating. In models with a sum-separable potential where the isocurvature modes decay before the end of the slow-roll phase we give an analytic criterion for the asymptotic value of fNL to be large. Other examples can be constructed using a waterfall field to terminate inflation while fNL is transiently large, caused by descent from a ridge or convergence into a valley. We show that these two types of evolution are distinguished by the sign of the bispectrum, and give approximate expressions for the peak fNL.

We calculate the conditions required to produce a large local trispectrum during two-field slow-roll inflation. This is done by extending and simplifying the 'heatmap' approach developed by Byrnes et al. The conditions required to generate a large trispectrum are broadly the same as those that can produce a large bispectrum. We derive a simple relation between tauNL and fNL for models with separable potentials, and furthermore show that gNL and tauNL can be related in specific circumstances. Additionally, we interpret the heatmaps dynamically, showing how they can be used as qualitative tools to understand the evolution of non-Gaussianity during inflation. We also show how fNL, tauNL and gNL are sourced by generic shapes in the inflationary potential, namely ridges, valleys and inflection points.

We study the hybrid inflationary potential in a regime where the defect field is light, and more than 60 e-folds of accelerated expansion occur after the symmetry breaking transition. Using analytic and numerical techniques, we then identify parameter values within this regime for which the statistics of the primordial curvature perturbation are significantly non-Gaussian. Focusing on this range of parameters, we provide a specific example which leads to an observationally consistent power spectrum, and a level of non-Gaussianity within current WMAP bounds and in reach of the Planck satellite. An interesting feature of this example is that the initial conditions at horizon crossing appear quite natural.

We review the current observational status of string cosmology when confronted with experimental datasets. We begin by defining common observational parameters and discuss how they are determined for a given model. Then we review the observable footprints of several string theoretic models, discussing the significance of various potential signals. Throughout we comment on present and future prospects of finding evidence for string theory in cosmology, and on significant issues for the future.

We study the evolution of non-Gaussianity in multiple-field inflationary models, focusing on three fundamental questions: (a) How is the sign and peak magnitude of the non-linearity parameter fNL related to generic features in the inflationary potential? (b) How sensitive is fNL to the process by which an adiabatic limit is reached, where the curvature perturbation becomes conserved? (c) For a given model, what is the appropriate tool -- analytic or numerical -- to calculate fNL at the adiabatic limit? We summarise recent results obtained by the authors and further elucidate them by considering an inflection point model.

In a recent publication, we proposed that inflationary perturbation theory can be reformulated in terms of a probability transport equation, whose moments determine the correlation properties of the primordial curvature perturbation. In this paper we generalize this formulation to an arbitrary number of fields. We deduce ordinary differential equations for the evolution of the moments of zeta on superhorizon scales, which can be used to obtain an evolution equation for the dimensionless bispectrum, fNL. Our equations are covariant in field space and allow identification of the source terms responsible for evolution of fNL. In a model with M scalar fields, the number of numerical integrations required to obtain solutions of these equations scales like O(M^3). The performance of the moment transport algorithm means that numerical calculations with M >> 1 fields are straightforward. We illustrate this performance with a numerical calculation of fNL in Nflation models containing M ~ 10^2 fields, finding agreement with existing analytic calculations. We comment briefly on extensions of the method beyond the slow-roll approximation, or to calculate higher order parameters such as gNL.

We study the dynamics of states perturbatively expanded about a harmonic system of loop quantum cosmology, exhibiting a bounce. In particular, the evolution equations for the first and second order moments of the system are analyzed. These moments back-react on the trajectories of the expectation values of the state and hence alter the energy density at the bounce. This analysis is performed for isotropic loop quantum cosmology coupled to a scalar field with a small but non-zero constant potential, hence in a regime in which the kinetic energy of matter dominates. Analytic restrictions on the existence of dynamical coherent states and the meaning of semi-classicality within these systems are discussed. A numerical investigation of the trajectories of states that remain semi-classical across the bounce demonstrates that, at least for such states, the bounce persists and that its properties are similar to the standard case, in which the moments of the states are entirely neglected. However the bounce density does change, implying that a quantum bounce may not be guaranteed to happen when the potential is no longer negligible.

In hybrid inflationary models, inflation ends by a sudden instability associated with a steep ridge in the potential. Here we argue that this feature can generate a large contribution to the curvature perturbation on observable scales. This contribution is almost scale-invariant but highly non-Gaussian. The degree of non-Gaussianity can exceed current observational bounds, unless the inflationary scale is extremely low or the hybrid potential contains very large coupling constants. Non-linear effects on small scales may quench the non-Gaussian signal, and while we find no compelling evidence that this occurs, full lattice simulations are required to definitively address this issue. Note added: We now believe that nonlinear effects will invalidate the original computation in this paper essentially instantaneously after the short-wavelength modes reach the minimum of their potential. This means that the mechanism described in this paper will not lead to appreciatable curvature perturbations on long wavelengths, and no useful constraints on hybrid inflation will result. We have inserted a brief calculation on p2 of this manuscript to explain this fact, but have otherwise left the manuscript unchanged.

We present a novel method for calculating the primordial non-Gaussianity produced by super-horizon evolution during inflation. Our method evolves the distribution of coarse-grained inflationary field values using a transport equation. We present simple evolution equations for the moments of this distribution, such as the variance and skewness. This method possesses some advantages over existing techniques. Among them, it cleanly separates multiple sources of primordial non-Gaussianity, and is computationally efficient when compared with popular alternatives, such as the "delta N" framework. We adduce numerical calculations demonstrating that our new method offers good agreement with those already in the literature. We focus on two fields and the fNL parameter, but we expect our method will generalize to multiple scalar fields and to moments of arbitrarily high order. We present our expressions in a field-space covariant form which we postulate to be valid for any number of fields.

Non-local equations of motion contain an infinite number of derivatives and commonly appear in a number of string theory models. We review how these equations can be rewritten in the form of a diffusion-like equation with non-linear boundary conditions. Moreover, we show that this equation can be solved as an initial value problem once a set of non-trivial initial conditions that satisfy the boundary conditions is found. We find these initial conditions by looking at the linear approximation to the boundary conditions. We then numerically solve the diffusion-like equation, and hence the non-local equations, as an initial value problem for the full non-linear potential and subsequently identify the cases when inflation is attained.

We investigate the behaviour of tensor fluctuations in Loop Quantum Cosmology, focusing on a class of scaling solutions which admit a near scale-invariant scalar field power spectrum. We obtain the spectral index of the gravitational field perturbations, and find a strong blue tilt in the power spectrum with $n_t \approx 2$. The amplitude of tensor modes are, therefore, suppressed by many orders of magnitude on large scales compared to those predicted by the standard inflationary scenario where $n_t \approx 0$.

There has been considerable recent interest in solving non-local equations of motion which contain an infinite number of derivatives. Here, focusing on inflation, we review how the problem can be reformulated as the question of finding solutions to a diffusion-like partial differential equation with non-linear boundary conditions. Moreover, we show that this diffusion-like equation, and hence the non-local equations, can be solved as an initial value problem once non-trivial initial data consistent with the boundary conditions is found. This is done by considering linearised equations about any field value, for which we show that obtaining solutions using the diffusion-like equation is equivalent to solving a local but infinite field cosmology. These local fields are shown to consist of at most two canonically normalized or phantom fields together with an infinite number of quintoms. We then numerically solve the diffusion-like equation for the full non-linear case for two string field theory motivated models.

A number of scalar field models proposed as alternatives to the standard inflationary scenario involve contracting phases which precede the universe's present phase of expansion. An important question concerning such models is whether there are effects which could potentially distinguish them from purely expanding cosmologies. Vector perturbations have recently been considered in this context. At first order such perturbations are not supported by a scalar field. In this paper, therefore, we consider second order vector perturbations. We show that such perturbations are generated by first order scalar mode-mode couplings, and give an explicit expression for them. We compare the magnitude of vector perturbations produced in collapsing models with the corresponding amplitudes produced during inflation, using a number of suitable power-law solutions to model the inflationary and collapsing scenarios. We conclude that the ratios of the magnitudes of these perturbations depend on the details of the collapsing scenario as well as on how the hot big bang is recovered, but for certain cases could be large, growing with the duration of the collapse.

The computation of the spectrum of primordial perturbations, generated by a scalar field during the super-inflationary phase of Loop Quantum Cosmology, is revisited. The calculation is performed for two different cases. The first considers the dynamics of a massless field and it is found that scale invariance can only be achieved under a severe fine tuning. The second assumes that the field evolves with a constant ratio between kinetic and potential energy, i.e. in a scaling solution. In this case, near scale invariance is a generic feature of the theory if the field rolls in a steep self interaction potential.

Positively-curved, oscillatory universes have recently been shown to have important consequences for the pre-inflationary dynamics of the early universe. In particular, they may allow a self-interacting scalar field to climb up its potential during a very large number of these cycles. The cycles are naturally broken when the potential reaches a critical value and the universe begins to inflate, thereby providing a `graceful entrance' to early universe inflation. We study the dynamics of this behaviour within the context of braneworld scenarios which exhibit a bounce from a collapsing phase to an expanding one. The dynamics can be understood by studying a general class of braneworld models that are sourced by a scalar field with a constant potential. Within this context, we determine the conditions a given model must satisfy for a graceful entrance to be possible in principle. We consider the bouncing braneworld model proposed by Shtanov and Sahni and show that it exhibits the features needed to realise a graceful entrance to inflation for a wide region of parameter space.

Closed, singularity-free, inflationary cosmological models have recently been studied in the context of general relativity. Despite their appeal, these so called emergent models suffer from a number of limitations. These include the fact that they rely on an initial Einstein static state to describe the past eternal phase of the universe. Given the instability of such a state within the context of general relativity, this amounts to a very severe fine tuning. Also in order to be able to study the dynamics of the universe within the context of general relativity, they set the initial conditions for the universe in the classical phase. Here we study the existence and stability of such models in the context of Loop Quantum Cosmology and show that both these limitations can be partially remedied, once semi-classical effects are taken into account. An important consequence of these effects is to give rise to a static solution (not present in GR), which dynamically is a centre equilibrium point and located in the more natural semi-classical regime. This allows the construction of emergent models in which the universe oscillates indefinitely about such an initial static state. We construct an explicit emergent model of this type, in which a non-singular past eternal oscillating universe enters a phase where the symmetry of the oscillations is broken, leading to an emergent inflationary epoch, while satisfying all observational and semi-classical constraints. We also discuss emergent models in which the universe possesses both early- and late-time accelerating phases.

Positively-curved, oscillatory universes are studied within the context of Loop Quantum Cosmology subject to a consistent semi-classical treatment. The semi-classical effects are reformulated in terms of an effective phantom fluid with a variable equation of state. In cosmologies sourced by a massless scalar field, these effects lead to a universe that undergoes ever-repeating cycles of expansion and contraction. The presence of a self-interaction potential for the field breaks the symmetry of the cycles and can enable the oscillations to establish the initial conditions for successful slow-roll inflation, even when the field is initially at the minimum of its potential with a small kinetic energy. The displacement of the field from its minimum is enhanced for lower and more natural values of the parameter that sets the effective quantum gravity scale. For sufficiently small values of this parameter, the universe can enter a stage of eternal self-reproduction.

In loop quantum gravity, modifications to the geometrical density cause a self-interacting scalar field to accelerate away from a minimum of its potential. In principle, this mechanism can generate the conditions that subsequently lead to slow-roll inflation. The consequences for this mechanism of various quantization ambiguities arising within loop quantum cosmology are considered. For the case of a quadratic potential, it is found that some quantization procedures are more likely to generate a phase of slow--roll inflation. In general, however, loop quantum cosmology is robust to ambiguities in the quantization and extends the range of initial conditions for inflation.