Inflationary models involving more than one scalar field naturally produce isocurvature perturbations. However, while these are fairly well studied, less is known about their evolution through the reheating epoch, when the inflationary fields decay into the standard constituents of the present universe. In this paper, by modelling reheating perturbatively, we calculate the power spectrum of the non-adiabatic pressure perturbation in three different inflationary models. We show that the isocurvature can grow large initially, but decays faster than the pressure perturbations. When reheating ends, the isocurvature is negligible for the double quadratic and double quartic inflationary models. For the product exponential potential, which features large isocurvature at the end of inflation, the isocurvature decays during reheating and is around five orders of magnitudes smaller than the pressure perturbation at the end of reheating.

How much does the curvature perturbation change after it leaves the horizon, and when should one evaluate the power spectrum? To answer these questions we study single field inflation models numerically, and compare the evolution of different curvature perturbations from horizon crossing to the end of inflation. In particular we calculate the number of efolds it takes for the curvature perturbation at a given wavenumber to settle down to within a given fraction of their value at the end of inflation. We find that e.g. in chaotic inflation, the amplitude of the comoving and the curvature perturbation on uniform density hypersurfaces differ by up to 180 % at horizon crossing assuming the same amplitude at the end of inflation, and that it takes approximately 3 efolds for the curvature perturbation to be within 1 % of its value at the end of inflation.

How much does the curvature perturbation change after it leaves the horizon, and when should one evaluate the power spectrum? To answer these questions we study single field inflation models numerically, and compare the evolution of different curvature perturbations from horizon crossing to the end of inflation. We find that e.g. in chaotic inflation, the amplitude of the comoving and the curvature perturbation on uniform density hypersurfaces differ by up to 180 \% at horizon crossing assuming the same amplitude at the end of inflation, and that it takes approximately 3 efolds for the curvature perturbation to be within 1 \% of its value at the end of inflation.

Non-adiabatic pressure perturbations naturally occur in models of inflation consisting of more than one scalar field. The amount of non-adiabatic pressure present at the end of inflation can have observational consequences through changes in the curvature perturbation, the generation of vorticity and subsequently the sourcing of B-mode polarisation. In this work, based on a presentation at the 13th Marcel Grossmann Meeting, we give a very brief overview of non-adiabatic pressure perturbations in multi-field inflationary models and describe our recent calculation of the spectrum of isocurvature perturbations generated at the end of inflation for different models which have two scalar fields.

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.

Isocurvature perturbations naturally occur in models of inflation consisting of more than one scalar field. In this paper we calculate the spectrum of isocurvature perturbations generated at the end of inflation for three different inflationary models consisting of two canonical scalar fields. The amount of non-adiabatic pressure present at the end of inflation can have observational consequences through the generation of vorticity and subsequently the sourcing of B-mode polarisation. We compare two different definitions of isocurvature perturbations and show how these quantities evolve in different ways during inflation. Our results are calculated using the open source Pyflation numerical package which is available to download.

We numerically calculate the evolution of second order cosmological perturbations for an inflationary scalar field without resorting to the slow-roll approximation or assuming large scales. In contrast to previous approaches we therefore use the full non-slow-roll source term for the second order Klein-Gordon equation which is valid on all scales. The numerical results are consistent with the ones obtained previously where slow-roll is a good approximation. We investigate the effect of localised features in the scalar field potential which break slow-roll for some portion of the evolution. The numerical package solving the second order Klein-Gordon equation has been released under an open source license and is available for download.

In this paper we summarise the status of single field models of inflation in light of the WMAP 7 data release. We find little has changed since the 5 year release, and results are consistent with previous findings. The increase in the upper bound on the running of the spectral index impacts on the status of the production of Primordial Black Holes from single field models. The lower bound on the equilateral configuration of the non-gaussianity parameter is reduced and thus the bounds on the theoretical parameters of (UV) DBI single brane models are weakened. In the case of multiple coincident branes the bounds are also weakened and the two, three or four brane cases will produce a tensor-signal that could possibly be observed in the future.

Inflationary cosmology is the leading explanation of the very early universe. Many different models of inflation have been constructed which fit current observational data. In this work theoretical and numerical methods for constraining the parameter space of a wide class of such models are described. First, string-theoretic models with large non-Gaussian signatures are investigated. An upper bound is placed on the amplitude of primordial gravitational waves produced by ultra-violet Dirac-Born-Infeld inflation. In all but the most finely tuned cases, this bound is incompatible with a lower bound derived for inflationary models which exhibit a red spectrum and detectable non-Gaussianity. By analysing general non-canonical actions, a class of models is found which can evade the upper bound when the phase speed of perturbations is small. The multi-coincident brane scenario with a finite number of branes is one such model. For models with a potentially observable gravitational wave spectrum the number of coincident branes is shown to take only small values. The second method of constraining inflationary models is the numerical calculation of second order perturbations for a general class of single field models. The Klein-Gordon equation at second order, written in terms of scalar field variations only, is numerically solved. The slow roll version of the second order source term is used and the method is shown to be extendable to the full equation. This procedure allows the evolution of second order perturbations in general and the calculation of the non-Gaussianity parameter in cases where there is no analytical solution available.

We numerically solve the Klein-Gordon equation at second order in cosmological perturbation theory in closed form for a single scalar field, describing the method employed in detail. We use the slow-roll version of the second order source term and argue that our method is extendable to the full equation. We consider two standard single field models and find that the results agree with previous calculations using analytic methods, where comparison is possible. Our procedure allows the evolution of second order perturbations in general and the calculation of the non-linearity parameter f_NL to be examined in cases where there is no analytical solution available.

A class of non-canonical inflationary models is identified, where the leading-order contribution to the non-Gaussianity of the curvature perturbation is determined by the sound speed of the fluctuations in the inflaton field. Included in this class of models is the effective action for multiple coincident branes in the finite n limit. The action for this configuration is determined using a powerful iterative technique, based upon the fundamental representation of SU(2). In principle the upper bounds on the tensor-scalar ratio that arise in the standard, single-brane DBI inflationary scenario can be relaxed in such multi-brane configurations if a large and detectable non-Gaussianity is generated. Moreover models with a small number of coincident branes could generate a gravitational wave background that will be observable to future experiments.

An upper bound on the amplitude of the primordial gravitational wave spectrum generated during ultra-violet DBI inflation is derived. The bound is insensitive to the form of the inflaton potential and the warp factor of the compactified dimensions and can be expressed entirely in terms of observational parameters once the volume of the five-dimensional sub-manifold of the throat has been specified. For standard type IIB compactification schemes, the bound predicts undetectably small tensor perturbations with a tensor-scalar ratio $r < 10^{-7}$. This is incompatible with a corresponding lower limit of $r > 0.1 (1-n_s)$, which applies to any model that generates a red spectral index $n_s <1$ and a potentially detectable non-Gaussianity in the curvature perturbation. Possible ways of evading these bounds in more general DBI-type scenarios are discussed and a multiple-brane model is investigated as a specific example.