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.

In this thesis we present research into linear perturbations in Lemaitre-Tolman-Bondi (LTB) and Assisted Coupled Quintessence (ACQ) Cosmologies. First we give a brief overview of the standard model of cosmology. We then introduce Cosmological Perturbation Theory (CPT) at linear order for a flat Friedmann-Robertson-Walker (FRW) cosmology. Next we study linear perturbations to a Lemaitre-Tolman-Bondi (LTB) background spacetime. Studying the transformation behaviour of the perturbations under gauge transformations, we construct gauge invariant quantities in LTB. We show, using the perturbed energy conservation equation, that there is a conserved quantity in LTB which is conserved on all scales. We then briefly extend our discussion to the Lemaitre spacetime, and construct gauge-invariant perturbations in this extension of LTB spacetime. We also study the behaviour of linear perturbations in assisted coupled quintessence models in a FRW background. We provide the full set of governing equations for this class of models, and solve the system numerically. The code written for this purpose is then used to evolve growth functions for various models and parameter values, and we compare these both to the standard $\Lambda$CDM model and to current and future observational bounds. We also 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 results of the full equations deviates from the approximation by more than the experimental uncertainty for these future surveys. The construction of the numerical code, PYESSENCE, written in Python to solve the system of background and perturbed evolution equations for assisted coupled quintessence, is also discussed.

This paper is a guide to the installation and use of the Python package PYESSENCE. PYESSENCE is designed to evolve linear perturbations to Coupled Quintessence models with a arbitrary number of cold dark matter (CDM) fluids and dark energy (DE) scalar fields as dictated by a given model. The equations are sufficiently general to allow for more exotic dark matter with a non-zero equation of state. Several example uses are included in order to demonstrate typical functionality to the potential user. PYESSENCE is released under an open source modified BSD license and is available on Bitbucket.

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.