We compute the one-loop density power spectrum including Newtonian and relativistic contributions, as well as the primordial non-Gaussianity contributions from $f_{\rm NL}$ and $g_{\rm NL}$ in the local configuration. To this end we take solutions to the Einstein equations in the long-wavelength approximation and provide expressions for the matter density perturbation at second and third order. These solutions have shown to be complementary to the usual Newtonian cosmological perturbations. We confirm a sub-dominant effect from pure relativistic terms, manifested at scales dominated by cosmic variance, but find that a sizable effect of order one comes from $g_{\rm NL}$ values allowed by Planck-2018 constraints, manifested at scales probed by forthcoming galaxy surveys like DESI and Euclid. As a complement, we present the matter bispectrum at the tree-level including the mentioned contributions.

A variety of gauges are used in cosmological perturbation theory. These are often chosen in order to attribute physical properties to a particular choice of coordinates, or otherwise to simplify the form of the resultant equations. Calculations are then performed with the understanding that they could have been done in any gauge, and that transformations between different gauges can be made at will. We show that this logic can be extended to the domain of large density contrasts, where different types of perturbative expansion are required, but that the way in which gauges can be chosen in the presence of such structures is severely constrained. In particular, most gauges that are commonly considered in the cosmology literature are found to be unviable in the presence of non-linear structures. This includes spatially flat gauge, synchronous gauge, comoving orthogonal gauge, total matter gauge, N-body gauge, and the uniform density gauge. In contrast, we find that the longitudinal gauge and the Newtonian motion gauge are both viable choices in both standard cosmological perturbation theory, and in the post-Newtonian perturbative expansions that are required in order to model non-linear structures.

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.

Next generation surveys will be capable of determining cosmological parameters beyond percent level. To match this precision, theoretical descriptions should look beyond the linear perturbations to approximate the observables in large scale structure. A quantity of interest is the Number density of galaxies detected by our instruments. This has been focus of interest recently, and several efforts have been made to explain relativistic effects theoretically, thereby testing the full theory. However, the results at nonlinear level from previous works are in disagreement. We present a new and independent approach to computing the relativistic galaxy number counts to second order in cosmological perturbation theory. We derive analytical expressions for the full second order relativistic observed redshift, for the angular diameter distance and for the volume spanned by a survey. Finally, we compare our results with previous works which compute the general distance-redshift relation, finding that our result is in agreement at linear order.

The generation of magnetic fields is a natural consequence of the existence of vortical currents in the pre-recombination era. This has been confirmed in detail for the case of adiabatic initial conditions, using second-order Boltzmann solvers, but has not been fully explored in the presence of isocurvatures. In this work, we use a modified version of the second-order Boltzmann code SONG to compute the magnetic field generated by vortical currents for general initial conditions. A mild enhancement of the generated magnetic field is found in the presence of general isocurvature modes, when compared to the adiabatic case. A particularly interesting case is that of the compensated isocurvature mode, for which the enhancement increases by several orders of magnitude due to the observationally allowed large amplitude of those modes. We show in this particular case how these compensated modes can influence observables at second order, such as the magnetic fields, and produce interesting effects which may be used to constrain these modes in the future.

We study a class of almost scale-invariant modified gravity theories, using a particular form of $f(R, G) = \alpha R^2 + \beta G \log G$ where $R$ and $G$ are the Ricci and Gauss-Bonnet scalars, respectively and $\alpha$, $\beta$ are arbitrary constants. We derive the Einstein-like field equations to first order in cosmological perturbation theory in longitudinal gauge.

We study how to set the initial evolution of general cosmological fluctuations at second order, after neutrino decoupling. We compute approximate initial solutions for the transfer functions of all the relevant cosmological variables sourced by quadratic combinations of adiabatic and isocurvature modes. We perform these calculations in synchronous gauge, assuming a Universe described by the $\Lambda$CDM model and composed of neutrinos, photons, baryons and dark matter. We highlight the importance of mixed modes, which are sourced by two different isocurvature or adiabatic modes and do not exist at the linear level. In particular, we investigate the so-called compensated isocurvature mode and find non-trivial initial evolution when it is mixed with the adiabatic mode, in contrast to the result at linear order and even at second order for the unmixed mode. Non-trivial evolution also arises when this compensated isocurvature is mixed with the neutrino density isocurvature mode. Regarding the neutrino velocity isocurvature mode, we show it unavoidably generates non-regular (decaying) modes at second order. Our results can be applied to second order Boltzmann solvers to calculate the effects of isocurvatures on non-linear observables.

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 introduce a double power series method for finding approximate analytical solutions for systems of differential equations commonly found in cosmological perturbation theory. The method was set out, in a non-cosmological context, by Feshchenko, Shkil' and Nikolenko (FSN) in 1966, and is applicable to cases where perturbations are on sub-horizon scales. The FSN method is essentially an extension of the well known Wentzel-Kramers-Brillouin (WKB) method for finding approximate analytical solutions for ordinary differential equations. The FSN method we use is applicable well beyond perturbation theory to solve systems of ordinary differential equations, linear in the derivatives, that also depend on a small parameter, which here we take to be related to the inverse wave-number. We use the FSN method to find new approximate oscillating solutions in linear order cosmological perturbation theory for a flat radiation-matter universe. Together with this model's well known growing and decaying M\'esz\'aros solutions, these oscillating modes provide a complete set of sub-horizon approximations for the metric potential, radiation and matter perturbations. Comparison with numerical solutions of the perturbation equations shows that our approximations can be made accurate to within a typical error of 1%, or better. We also set out a heuristic method for error estimation. A Mathematica notebook which implements the double power series method is made available online.

We propose and construct a two-parameter perturbative expansion around a Friedmann-Lema\^{i}tre-Robertson-Walker geometry that can be used to model high-order gravitational effects in the presence of non-linear structure. This framework reduces to the weak-field and slow-motion post-Newtonian treatment of gravity in the appropriate limits, but also includes the low-amplitude large-scale fluctuations that are important for cosmological modelling. We derive a set of field equations that can be applied to the late Universe, where non-linear structure exists on supercluster scales, and perform a detailed investigation of the associated gauge problem. This allows us to identify a consistent set of perturbed quantities in both the gravitational and matter sectors, and to construct a set of gauge-invariant quantities that correspond to each of them. The field equations, written in terms of these quantities, take on a relatively simple form, and allow the effects of small-scale structure on the large-scale properties of the Universe to be clearly identified. We find that inhomogeneous structures source the global expansion, that there exist new field equations at new orders, and that there is vector gravitational potential that is a hundred times larger than one might naively expect from cosmological perturbation theory. Finally, we expect our formalism to be of use for calculating relativistic effects in upcoming ultra-large-scale surveys, as the form of the gravitational coupling between small and large scales depends on the non-linearity of Einstein's equations, and occurs at what is normally thought of as first order in cosmological perturbations.

We derive the evolution equation for the second order curvature perturbation using standard techniques of cosmological perturbation theory. We do this for different definitions of the gauge invariant curvature perturbation, arising from different splits of the spatial metric, and compare the expressions. The results are valid at all scales and include all contributions from scalar, vector and tensor perturbations, as well as anisotropic stress, with all our results written purely in terms of gauge invariant quantities. Taking the large-scale approximation, we find that a conserved quantity exists only if, in addition to the non-adiabatic pressure, the transverse traceless part of the anisotropic stress tensor is also negligible. We also find that the version of the gauge invariant curvature perturbation which is exactly conserved is the one defined with the determinant of the spatial part of the inverse metric.

We use gauge-invariant cosmological perturbation theory to calculate the displacement field that sets the initial conditions for $N$-body simulations. Using first and second-order fully relativistic perturbation theory in the synchronous-comoving gauge, allows us to go beyond the Newtonian predictions and to calculate relativistic corrections to it. We use an Einstein--de Sitter model, including both growing and decaying modes in our solutions. The impact of our results should be assessed through the implementation of the featured displacement in cosmological $N$-body simulations.

Working with perturbations about an FLRW spacetime, we compute the gauge-invariant curvature perturbation to second order solely in terms of scalar field fluctuations. Using the curvature perturbation on uniform density hypersurfaces as our starting point, we give our results in terms of field fluctuations in the flat gauge, incorporating both large and small scale behaviour. For ease of future numerical implementation we give our result in terms of the scalar field fluctuations and their time derivatives.

We study linear perturbations to a Lema{\^\i}tre-Tolman-Bondi (LTB) background spacetime. Studying the transformation behaviour of the perturbations under gauge transformations, we construct gauge invariant quantities. We show, using the perturbed energy conservation equation, that there are conserved quantities in LTB, in particular a spatial metric trace perturbation, \zeta_{SMTP}, which is conserved on all scales. We then briefly extend our discussion to the Lema{\^\i}tre spacetime, and construct gauge-invariant perturbations in this extension of LTB spacetime.

Magnetic fields are present on all scales in the Universe. While we understand the processes which amplify the fields fairly well, we do not have a "natural" mechanism to generate the small initial seed fields. By using fully relativistic cosmological perturbation theory and going beyond the usual confines of linear theory we show analytically how magnetic fields are generated. This is the first analytical calculation of the magnetic field at second order, using gauge-invariant cosmological perturbation theory, and including all the source terms. To this end, we have rederived the full set of governing equations independently. Our results suggest that magnetic fields of the order of $10^{-30}$ G can be generated (although this depends on the small scale cut-off of the integral), which is largely in agreement with previous results that relied upon numerical calculations. These fields are likely too small to act as the primordial seed fields for dynamo mechanisms.

We study the production of primordial black holes (PBHs) during the preheating stage that follows a chaotic inflationary phase. The scalar fields present in the process are evolved numerically using a modified version of the HLATTICE code. From the output of the numerical simulation, we compute the probability distribution of curvature fluctuations, paying particular attention to sub-horizon scales. We find that in some specific models these modes grow to large amplitudes developing highly non-Gaussian probability distributions. We then calculate PBH abundances using the standard Press-Schechter criterion and find that overproduction of PBHs is likely in some regions of the chaotic preheating parameter space.

We revisit the gauge issue in cosmological perturbation theory, and highlight its relation to the notion of covariance in general relativity. We also discuss the similarities and differences of the covariant approach in perturbation theory to the Bardeen or metric approach in a non-technical fashion.

We calculate the relativistic constraint equation which relates the curvature perturbation to the matter density contrast at second order in cosmological perturbation theory. This relativistic "second order Poisson equation" is presented in a gauge where the hydrodynamical inhomogeneities coincide with their Newtonian counterparts exactly for a perfect fluid with constant equation of state. We use this constraint to introduce primordial non-Gaussianity in the density contrast in the framework of General Relativity. We then derive expressions that can be used as initial conditions of N-body codes for structure formation which probe the observable signature of primordial non-Gaussianity in the statistics of the evolved matter density field.

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.

Currently, most of the numerical simulations of structure formation use Newtonian gravity. When modelling pressureless dark matter, or `dust', this approach gives the correct results for scales much smaller than the cosmological horizon, but for scenarios in which the fluid has pressure this is no longer the case. In this article, we present the correspondence of perturbations in Newtonian and cosmological perturbation theory, showing exact mathematical equivalence for pressureless matter, and giving the relativistic corrections for matter with pressure. As an example, we study the case of scalar field dark matter which features non-zero pressure perturbations. We discuss some problems which may arise when evolving the perturbations in this model with Newtonian numerical simulations and with CMB Boltzmann codes.

Vorticity is ubiquitous in nature however, to date, studies of vorticity in cosmology and the early universe have been quite rare. In this paper, based on a talk in session CM1 of the 13th Marcel Grossmann Meeting, we consider vorticity generation from scalar cosmological perturbations of a perfect fluid system. We show that, at second order in perturbation theory, vorticity is sourced by a coupling between energy density and entropy gradients, thus extending a well-known feature of classical fluid dynamics to a relativistic cosmological framework. This induced vorticity, sourced by isocurvature perturbations, may prove useful in the future as an additional discriminator between inflationary models.

We study the effect of curvaton decay on the primordial power spectrum. Using analytical approximations and also numerical calculations, we find that the power spectrum is enhanced during the radiation dominated era after the curvaton decay. The amplitude of the Bardeen potential is controlled by the fraction of the energy density in the curvaton at the time of curvaton decay. We show that the enhancement in the amplitude of the primordial curvature perturbation is, however, not large enough to lead to primordial black hole overproduction on scales which re-enter the horizon after the time of curvaton decay.

Understanding the non-adiabatic pressure, or relative entropy, perturbation is crucial for studies of early-universe vorticity and Cosmic Microwave Background observations. We calculate the evolution of the linear non-adiabatic pressure perturbation from radiation domination to late times, numerically solving the linear governing equations for a wide range of wavenumbers. Using adiabatic initial conditions consistent with WMAP seven year data, we find nevertheless that the non-adiabatic pressure perturbation is non-zero and grows at early times, peaking around the epoch of matter/radiation equality and decaying in matter domination. At early times or large redshifts (z=10,000) its power spectrum peaks at a comoving wavenumber k~0.2h/Mpc, while at late times (z=500) it peaks at k~0.02 h/Mpc.

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.

We compare and contrast two different metric based formulations of non- linear cosmological perturbation theory: the MW2009 approach in [K. A. Malik and D. Wands, Phys. Rept. 475 (2009), 1.] following Bardeen and the recent approach of the paper KN2010 [K. Nakamura, Advances in Astronomy 2010 (2010), 576273]. We present each formulation separately. In the MW2009 approach, one considers the gauge transformations of perturbative quantities, choosing a gauge by requiring that certain quantities vanish, rendering all other variables gauge invariant. In the KN2010 formalism, one decomposes the metric tensor into a gauge variant and gauge invariant part from the outset. We compare the two approaches in both the longitudinal and uniform curvature gauges. In the longitudinal gauge, we find that Nakamura's gauge invariant variables correspond exactly to those in the longitudinal gauge (i.e., for scalar perturbations, to the Bardeen potentials), and in the uniform curvature gauge we obtain the usual relationship between gauge invariant variables in the flat and longitudinal gauge. Thus, we show that these two approaches are equivalent.

In this paper we analyse three models of the early universe, for which the respective mechanisms for generating the curvature perturbation are considered disparate. We find that in fact the mechanisms are very similar, and hence explain why they give rise to a large non-gaussianity. We show that the mechanism for generating the primordial curvature perturbation, and hence the observable non-gaussianity, is similar in both the Curvaton and Modulated Reheating models. In both cases the model can be written in terms of an energy transfer between the constituting fluids. We then show that this is also true for the mechanism of generating the curvature perturbation by symmetry breaking the end of inflation. We then relate this to the non-gaussian contribution to the curvature perturbation and find that it is inversely proportional to the efficiency with which the curvature perturbation is transferred between the fluids. For the first time, we generalise models of modulated reheating to allow for a non-linear energy transfer rate.

In this concise, albeit subjective review of structure formation, I shall introduce the cosmological standard model and its theoretical and observational underpinnings. I will focus on recent results and current issues in theoretical cosmology, in particular in cosmological perturbation theory and its applications.

In this special issue article, based on the talk with the same title in session B5 (Theoretical and Mathematical Cosmology) at GR19, we review the case of vorticity generation in cosmology using cosmological perturbation theory. We show that, while at linear order the vorticity evolution equation has no source term in the absence of anisotropic stress, at second order vorticity is sourced by gradients in entropy and energy density perturbations. We then present some estimates for the magnitude and scale dependence of the vorticity power spectrum using simple input power spectra for the energy density and entropy perturbations. Finally, we close with possible directions for future work followed by some hints toward the observational importance of the vorticity so generated, and the possibility of primordial magnetic field generation.

We estimate the amount of vorticity generated at second order in cosmological perturbation theory from the coupling between first order energy density and non-adiabatic pressure, or entropy, perturbations. Assuming power law input spectra for the source terms, and working in a radiation background, we calculate the wave number dependence of the vorticity power spectrum and its amplitude. We show that the vorticity generated by this mechanism is non-negligible on small scales, and hence should be taken into consideration in current and future CMB experiments.

We present a formalism for spatial averaging in cosmology applicable to general spacetimes and coordinates, and allowing the easy incorporation of a wide variety of matter sources. We apply this formalism to a Friedmann-LeMaitre-Robertson-Walker universe perturbed to second-order and present the corrections to the background in an unfixed gauge. We then present the corrections that arise in uniform curvature and conformal Newtonian gauges.

We discuss cosmological perturbation theory at third order, deriving the gauge transformation rules for metric and matter perturbations, and constructing third order gauge invariant quantities. We present the Einstein tensor components, the evolution equations for a perfect fluid, and the Klein-Gordon equation at third order, including scalar, vector and tensor perturbations. In doing so, we also give all second order tensor components and evolution equations in full, exhilarating generality.

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.

We show that at second order in cosmological perturbation theory vorticity generation is sourced by entropy gradients. This is an extension of Crocco's theorem to a cosmological setting.

We review the study of inhomogeneous perturbations about a homogeneous and isotropic background cosmology. We adopt a coordinate based approach, but give geometrical interpretations of metric perturbations in terms of the expansion, shear and curvature of constant-time hypersurfaces and the orthogonal timelike vector field. We give the gauge transformation rules for metric and matter variables at first and second order. We show how gauge invariant variables are constructed by identifying geometric or matter variables in physically-defined coordinate systems, and give the relations between many commonly used gauge-invariant variables. In particular we show how the Einstein equations or energy-momentum conservation can be used to obtain simple evolution equations at linear order, and discuss extensions to non-linear order. We present evolution equations for systems with multiple interacting fluids and scalar fields, identifying adiabatic and entropy perturbations. As an application we consider the origin of primordial curvature and isocurvature perturbations from field perturbations during inflation in the very early universe.

Primordial black holes (PBHs) can form in the early Universe via the collapse of large density perturbations. There are tight constraints on the abundance of PBHs formed due to their gravitational effects and the consequences of their evaporation. These abundance constraints can be used to constrain the primordial power spectrum, and hence models of inflation, on scales far smaller than those probed by cosmological observations. We compile, and where relevant update, the constraints on the abundance of PBHs before calculating the constraints on the curvature perturbation, taking into account the growth of density perturbations prior to horizon entry. We consider two simple parameterizations of the curvature perturbation spectrum on the scale of interest: constant and power-law. The constraints from PBHs on the amplitude of the power spectrum are typically in the range 10^{-2}-10^{-1} with some scale dependence.

We discuss the non-adiabatic or entropy perturbation, which controls the evolution of the curvature perturbation in the uniform density gauge, for a scalar field system minimally coupled to gravity with non-canonical action. We highlight the differences between the sound and the phase speed in these systems, and show that the non-adiabatic pressure perturbation vanishes in the single field case, resulting in the conservation of the curvature perturbation on large scales.

We give a concise, self-contained introduction to perturbation theory in cosmology at linear and second order, striking a balance between mathematical rigour and usability. In particular we discuss gauge issues and the active and passive approach to calculating gauge transformations. We also construct gauge-invariant variables, including the second order tensor perturbation on uniform curvature hypersurfaces.

We give the governing equations for multiple scalar fields in a flat Friedmann-Robertson-Walker (FRW) background spacetime on all scales, allowing for metric and field perturbations up to second order. We then derive the Klein-Gordon equation at second order in closed form in terms of gauge-invariant perturbations of the fields in the uniform curvature gauge. We also give a simplified form of the Klein-Gordon equation using the slow-roll approximation.

We calculate the tree-level bispectrum of the inflaton field perturbation directly from the field equations, and construct the corresponding f_NL parameter. Our results agree with previous ones derived from the Lagrangian. We argue that quantum theory should only be used to calculate the correlators when they first become classical a few Hubble times after horizon exit, the classical evolution taking over thereafter.

We derive the Klein--Gordon equation for a single scalar field coupled to gravity at second order in perturbation theory and leading order in slow-roll. This is done in two ways: we derive the Klein--Gordon equation first using the Einstein field equations, and then directly from the action after integrating out the constraint equations. We also point out an unexpected result regarding the treatment of the field equations.

Oscillating moduli fields can support a cosmological scaling solution in the presence of a perfect fluid when the scalar field potential satisfies appropriate conditions. We examine when such conditions arise in higher-dimensional, non-linear sigma-models that are reduced to four dimensions under a generalized Scherk-Schwarz compactification. We show explicitly that scaling behaviour is possible when the higher-dimensional action exhibits a global SL(n,R) or O(2,2) symmetry. These underlying symmetries can be exploited to generate non-trivial scaling solutions when the moduli fields have non-canonical kinetic energy. We also consider the compactification of eleven-dimensional vacuum Einstein gravity on an elliptic twisted torus.

We calculate the constraints on the primordial curvature perturbation at the end of inflation from the present day abundance of Primordial Black Holes (PBHs), as a function of the reheat temperature T_{\rm RH}. We first extend recent work on the formation of PBHs on scales which remain within the horizon during inflation and calculate the resulting constraints on the curvature perturbation. We then evaluate the constraint from PBHs that form, more conventionally, from super-horizon perturbations. The constraints apply for T_{\rm RH} < 10^{8} GeV and the inclusion of sub-horizon PBHs leads to a limit which is roughly three times tighter than the bound from super-horizon PBHs.

We study the curvaton scenario using gauge-invariant second order perturbation theory and solving the governing equations numerically. Focusing on large scales we calculate the non-linearity parameter f_nl in the two-fluid curvaton model and compare our results with previous analytical studies employing the sudden decay approximation. We find good agreement of the two approaches for large curvaton energy densities at curvaton decay, Omega_dec, but significant differences of up to 10 percent for small Omega_dec.

We derive the governing equations for multiple scalar fields minimally coupled to gravity in a flat Friedmann-Robertson-Walker (FRW) background spacetime on large scales. We include scalar perturbations up to second order and write the equations in terms of physically transparent gauge-invariant variables at first and second order. This allows us to write the perturbed Klein-Gordon equation at second order solely in terms of the field fluctuations on flat slices at first and second order.

Previous authors have calculated the mass function of primordial black holes only on scales which are well outside the horizon at the end of inflation. Here we extend the calculation to sub-horizon scales, on which the density perturbation never becomes classical. Regarding the formation of black holes as a `measurement' of the (high peaks) of the density perturbation, we estimate a mass function by assuming that black holes form as soon as inflation ends, in those rare regions where the Bardeen potential exceeds a threshold value of $Psi_c\simeq 0.5$.

Without invoking a perturbative expansion, we define the cosmological curvature perturbation, and consider its behaviour assuming that the universe is smooth over a sufficiently large comoving scale. The equations are simple, resembling closely the first-order equations, and they lead to results which generalise those already proven in linear perturbation theory and (in part) in second-order perturbation theory. In particular, the curvature perturbation is conserved provided that the pressure is a unique function of the energy density.

We develop a gauge-invariant formalism for the study of density perturbations in a Friedmann-Robertson-Walker universe with multiple interacting fluids and/or scalar fields. We show how N scalar fields may be described by N kinetic fluids (with maximally stiff equation of state) interacting with a non-dynamical potential (with vacuum equation of state). We split generic perturbations into adiabatic and entropic parts, and give the coupled first-order evolution equations on all scales, including energy and momentum exchange. We identify the non-adiabatic effects on large scales, and define adiabatic initial conditions in the presence of multiple fluids and fields.

We revisit the calculation of the abundance of primordial black holes (PBHs) formed from primordial density perturbations, using a formation criterion derived by Shibata and Sasaki which refers to a metric perturbation variable rather than the usual density contrast. We implement a derivation of the PBH abundance which uses peaks theory, and compare it to the standard calculation based on a Press--Schechter-like approach. We find that the two are in reasonable agreement if the Press--Schechter threshold is in the range $\Delta_{{\rm th}} \simeq 0.3$ to 0.5, but advocate use of the peaks theory expression which is based on a sounder theoretical footing.

We present a method for constructing gauge-invariant cosmological perturbations which are gauge-invariant up to second order. As an example we give the gauge-invariant definition of the second-order curvature perturbation on uniform density hypersurfaces. Using only the energy conservation equation we show that this curvature perturbation is conserved at second order on large scales for adiabatic perturbations.

In the curvaton model of inflation, where a second scalar field, the "curvaton", is responsible for the observed inhomogeneity, a non-zero neutrino degeneracy may lead to a characteristic pattern of isocurvature perturbations in the neutrino, cold dark matter and baryon components. We find the current data can only place upper limits on the level of isocurvature perturbations. These can be translated into upper limits on the neutrino degeneracy parameter. In the case that lepton number is created before curvaton decay, we find that the limit on the neutrino degeneracy parameter is comparable with that obtained from Big-bang nucleosynthesis. For the case that lepton number is created by curvaton decay we find that the absolute value of the non-Gaussianity parameter, |f_nl|, must be less than 10 (95% confidence interval).

We study the evolution of the cosmological perturbations after inflation in curvaton models where the non-relativistic curvaton decays into both radiation and a cold dark matter component. We calculate the primordial curvature and correlated isocurvature perturbations inherited by the radiation and cold dark matter after the curvaton has decayed. We give the transfer coefficient in terms of the initial curvaton density relative to the curvaton decay rate.

We extend the Kelvin-Helmholtz instability to an expanding background. We study the evolution of a non-viscous irrotational fluid and find that for wavelengths much smaller than the Hubble scale small perturbations of the fluid are unstable for wavenumbers larger than a critical value. We then apply this result in the early universe, treating cold dark matter as a classical fluid with vanishing background pressure.

We study the definition of perturbations in the presence of a submanifold, like e.g. a brane. In the standard theory of cosmological perturbations, one compares quantities at the same coordinate points in the non-perturbed and the perturbed manifolds, identified via a (non-unique) mapping between the two manifolds. In the presence of a physical submanifold one needs to modify this definition in order to evaluate perturbations of quantities at the submanifold location. As an application, we compute the perturbed metric and the extrinsic curvature tensors at the brane position in a general gauge.

We present a gauge-invariant formalism to study the evolution of curvature perturbations in a Friedmann-Robertson-Walker universe filled by multiple interacting fluids. We resolve arbitrary perturbations into adiabatic and entropy components and derive their coupled evolution equations. We demonstrate that perturbations obeying a generalised adiabatic condition remain adiabatic in the large-scale limit, even when one includes energy transfer between fluids. As a specific application we study the recently proposed curvaton model, in which the curvaton decays into radiation. We use the coupled evolution equations to show how an initial isocurvature perturbation in the curvaton gives rise to an adiabatic curvature perturbation after the curvaton decays.

We study cosmological perturbations in a brane-world scenario where the matter fields live on a four-dimensional brane and gravity propagates in the five-dimensional bulk. We present the equations of motion in an arbitrary gauge for metric perturbations in the bulk and matter perturbations on the brane. Gauge-invariant perturbations are then constructed corresponding to perturbations in longitudinal and Gaussian normal gauges. Longitudinal gauge metric perturbations may be directly derived from three master variables (separately describing scalar, vector and tensor metric perturbations) which obey five-dimensional wave-equations. Gaussian normal gauge perturbations are directly related to the induced metric perturbations on the brane with the additional bulk degrees of freedom interpreted as an effective Weyl energy-momentum tensor on the brane. We construct gauge-invariant perturbations describing the effective density, momentum and pressures of this Weyl fluid at the brane and throughout the bulk. We show that there exist gauge-invariant curvature perturbations on the brane and in the bulk that are conserved on large-scales when three-dimensional spatial gradients are negligible.

During the preheating process at the end of inflation the amplification of field fluctuations can lead to the amplification of curvature perturbations. If the curvature perturbations on small scales are sufficiently large, primordial black holes (PBHs) will be overproduced. In this paper we study PBH production in the two-field preheating model with quadratic inflaton potential. We show that for many values of the inflaton mass m, and coupling g, small scale perturbations will be amplified sufficiently, before backreaction can shut off preheating, so that PBHs will be overproduced during the subsequent radiation dominated era.

After introducing gauge-invariant cosmological perturbation theory we give an improved set of governing equations for multiple fluids including energy transfer. Having defined adiabatic and entropic perturbations we derive the ``conservation law'' for the curvature perturbation on large scales using only the energy conservation equation. We then investigate the dynamics of assisted inflation. By choosing an appropriate rotation in field space we can write down explicitly the potential for the weighted mean field along the scaling solution and for fields orthogonal to it. This allows us to present analytic solutions describing homogeneous and inhomogeneous perturbations about the attractor solution without resorting to slow-roll approximations. Finally we analyze the simplest model of preheating analytically, and show that in linear perturbation theory the effect of preheating on the amplitude of the curvature perturbation on large scales is negligible. We end with some concluding remarks, possible extensions and an outlook to future work.

We discuss the evolution of linear perturbations about a Friedmann-Robertson-Walker background metric, using only the local conservation of energy-momentum. We show that on sufficiently large scales the curvature perturbation on spatial hypersurfaces of uniform-density is constant when the non-adiabatic pressure perturbation is negligible. We clarify the conditions under which super-horizon curvature perturbations may vary, using preheating as an example.

We discuss the evolution of linear perturbations about a Friedmann-Robertson-Walker background metric, using only the local conservation of energy-momentum. We show that on sufficiently large scales the curvature perturbation on spatial hypersurfaces of uniform-density is conserved when the non-adiabatic pressure perturbation is negligible. This is the first time that this result has been demonstrated independently of the gravitational field equations. A physical picture of long-wavelength perturbations as being composed of separate Robertson-Walker universes gives a simple understanding of the possible evolution of the curvature perturbation, in particular clarifying the conditions under which super-horizon curvature perturbations may vary.

It has recently been claimed by Bassett et al that preheating after inflation may affect the amplitude of curvature perturbations on large scales, undermining the usual inflationary prediction. We analyze the simplest model, and confirm the results of Jedamzik and Sigl and of Ivanov that in linear perturbation theory the effect is negligible. However the dominant effect is second-order in the field perturbation and we show that this too is negligible, and hence conclude that preheating has no significant influence on large-scale perturbations in this model. We briefly discuss the likelihood of an effect in other models.

We investigate the dynamics of the recently proposed model of assisted inflation. In this model an arbitrary number of scalar fields with exponential potentials evolve towards an inflationary scaling solution, even if each of the individual potentials is too steep to support inflation on its own. By choosing an appropriate rotation in field space we can write down explicitly the potential for the weighted mean field along the scaling solution and for fields orthogonal to it. This demonstrates that the potential has a global minimum along the scaling solution. We show that the potential close to this attractor in the rotated field space is analogous to a hybrid inflation model, but with the vacuum energy having an exponential dependence upon a dilaton field. We present analytic solutions describing homogeneous and inhomogeneous perturbations about the attractor solution without resorting to slow-roll approximations. We discuss the curvature and isocurvature perturbation spectra produced from vacuum fluctuations during assisted inflation.

We show how gauge-invariant cosmological perturbations may be constructed by an unambiguous choice of hypersurface-orthogonal time-like vector field (i.e., time-slicing). This may be defined either in terms of the metric quantities such as curvature or shear, or using some matter field. As an example, we show how linear perturbations in the covariant fluid-flow approach can then be presented in an explicitly gauge-invariant form in the coordinate based formalism.