Motivated by tension between the predictions of ordinary cold dark matter (CDM) and observations at galactic scales, ultralight axionlike particles (ULALPs) with mass of the order $10^{-22}~{\rm eV}$ have been proposed as an alternative CDM candidate. We consider cold and collisionless ULALPs produced in the early Universe by the vacuum realignment mechanism and constituting most of CDM. The ULALP fluid is commonly described by classical field equations. However, we show that, like QCD axions, the ULALPs thermalize by gravitational self-interactions and form a Bose-Einstein condensate, a quantum phenomenon. ULALPs, like QCD axions, explain the observational evidence for caustic rings of dark matter because they thermalize and go to the lowest energy state available to them. This is one of rigid rotation on the turnaround sphere. By studying the heating effect of infalling ULALPs on galactic disk stars and the thickness of the nearby caustic ring as observed from a triangular feature in the IRAS map of our galactic disk, we obtain lower-mass bounds on the ULALP mass of order $10^{-23}$ and $10^{-20}~{\rm eV}$, respectively.

Modified theories of gravity have been invoked recently as an alternative to dark energy, in an attempt to explain the apparent accelerated expansion of the universe at the present time. In order to describe inhomogeneities in cosmological models, cosmological perturbation theory is used, of which two formalisms exist: the metric approach and the covariant approach. In this paper I present the relationship between the metric and covariant approaches for modeling $f(R)$ theories of gravity. This provides a useful resource that researchers primarily working with one formalism can use to compare or translate their results to the other formalism.

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.

We introduce a new mechanism for generating magnetic fields in the recombination era. This Harrison-like mechanism utilizes vorticity in baryons that is sourced through the Bose-Einstein condensate of axions via gravitational interactions. The magnetic fields generated are on the galactic scales $\sim 10\,{\rm kpc}$ and have a magnitude of the order of $B\sim10^{-23}\,{\rm G}$ today. The field has a greater magnitude than those generated from other mechanisms relying on second order perturbation theory, and is sufficient to provide a seed for battery mechanisms.

The axions produced during the QCD phase transition by vacuum realignment, string decay and domain wall decay thermalize as a result of their gravitational self-interactions when the photon temperature is approximately 500 eV. They then form a Bose-Einstein condensate (BEC). Because the axion BEC rethermalizes on time scales shorter than the age of the universe, it has properties that distinguish it from other forms of cold dark matter. The observational evidence for caustic rings of dark matter in galactic halos is explained if the dark matter is axions, at least in part, but not if the dark matter is entirely WIMPs or sterile neutrinos.

We consider a model of the early universe which consists of two scalar fields: the inflaton, and a second field which drives the stabilisation of the Planck mass (or gravitational constant). We show that the non-minimal coupling of this second field to the Ricci scalar sources a non-adiabatic pressure perturbation. By performing a fully numerical calculation we find, in turn, that this boosts the amplitude of the primordial power spectrum after inflation.

We obtain solutions to the coupled Schr\"odinger-Poisson equations. The solutions describe the evolution of cold dark matter density perturbations in an otherwise homogeneous expanding Friedmann universe. We discuss the relationships between descriptions of cold dark matter in terms of a pressureless fluid, in terms of a wavefunction, of a classical scalar field, and a quantum scalar field. We identify the regimes where the various descriptions coincide and where they differ.

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.

In this paper I review some recent, interlinked, work undertaken using cosmological perturbation theory -- a powerful technique for modelling inhomogeneities in the Universe. The common theme which underpins these pieces of work is the presence of non-adiabatic pressure, or entropy, perturbations. After a brief introduction covering the standard techniques of describing inhomogeneities in both Newtonian and relativistic cosmology, I discuss the generation of vorticity. As in classical fluid mechanics, vorticity is not present in linearized perturbation theory (unless included as an initial condition). Allowing for entropy perturbations, and working to second order in perturbation theory, I show that vorticity is generated, even in the absence of vector perturbations, by purely scalar perturbations, the source term being quadratic in the gradients of first order energy density and isocurvature, or non-adiabatic pressure perturbations. This generalizes Crocco's theorem to a cosmological setting. I then introduce isocurvature perturbations in different models, focusing on the entropy perturbation in standard, concordance cosmology, and in inflationary models involving two scalar fields. As the final topic, I investigate magnetic fields, which are a potential observational consequence of vorticity in the early universe. I briefly review some recent work on including magnetic fields in perturbation theory in a consistent way. I show, using solely analytical techniques, that magnetic fields can be generated by higher order perturbations, albeit too small to provide the entire primordial seed field, in agreement with some numerical studies. I close with a summary and some potential extensions of this work.

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.

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.

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.

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.

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.

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.

Cosmological perturbation theory is a key tool to study the universe. The linear or first order theory is well understood, however, developing and applying the theory beyond linear order is at the cutting edge of current research in theoretical cosmology. In this article, I will describe some signatures of non-linear perturbation theory that do not exist at linear order, focusing on vorticity generation at second order. In doing so, we discuss why this, among other features such as induced gravitational waves and non-Gaussianities, shows that cosmological perturbation theory is crucial for testing models of the universe.

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.

Cosmological perturbation theory is crucial for our understanding of the universe. The linear theory has been well understood for some time, however developing and applying the theory beyond linear order is currently at the forefront of research in theoretical cosmology. This thesis studies the applications of perturbation theory to cosmology and, specifically, to the early universe. Starting with some background material introducing the well-tested 'standard model' of cosmology, we move on to develop the formalism for perturbation theory up to second order giving evolution equations for all types of scalar, vector and tensor perturbations, both in gauge dependent and gauge invariant form. We then move on to the main result of the thesis, showing that, at second order in perturbation theory, vorticity is sourced by a coupling term quadratic in energy density and entropy perturbations. This source term implies a qualitative difference to linear order. Thus, while at linear order vorticity decays with the expansion of the universe, the same is not true at higher orders. This will have important implications on future measurements of the polarisation of the Cosmic Microwave Background, and could give rise to the generation of a primordial seed magnetic field. Having derived this qualitative result, we then estimate the scale dependence and magnitude of the vorticity power spectrum, finding, for simple power law inputs a small, blue spectrum. The final part of this thesis concerns higher order perturbation theory, deriving, for the first time, the metric tensor, gauge transformation rules and governing equations for fully general third order perturbations. We close with a discussion of natural extensions to this work and other possible ideas for off-shooting projects in this continually growing field.

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 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.

When analysing a system consisting of both dark matter and dark energy, an often used practice in the literature is to neglect the perturbations in the dark energy component. However, it has recently been argued, through the use of numerical simulations, that one cannot do so. In this work we show that by neglecting such perturbations one is implicitly making a choice of gauge. As such, one no longer has the freedom to choose, for example, a gauge comoving with the dark matter -- in fact doing so will give erroneous, gauge dependent results. We obtain results consistent with the numerical simulations by using the formalism of cosmological perturbation theory, and thus without resorting to involved numerical calculations.

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 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 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.