Searches for primordial gravitational waves have resulted in constraints in a large frequency range from a variety of sources. The standard Cosmic Microwave Background (CMB) technique is to parameterise the tensor power spectrum in terms of the tensor-to-scalar ratio, $r$, and spectral index, $n_{\rm t}$, and constrain these using measurements of the temperature and polarization power spectra. Another method, applicable to modes well inside the cosmological horizon at recombination, uses the shortwave approximation, under which gravitational waves behave as an effective neutrino species. In this paper we give model-independent CMB constraints on the energy density of gravitational waves, $\Omega_\textrm{gw} h^2$, for the entire range of observable frequencies. On large scales, $f \lesssim 10^{-16}\, \text{Hz}$, we reconstruct the initial tensor power spectrum in logarithmic frequency bins, finding maximal sensitivity for scales close to the horizon size at recombination. On small scales, $f \gtrsim10^{-15}\,\mbox{Hz}$, we use the shortwave approximation, finding $\Omega_\textrm{gw} h^2 < 1.7 \times10^{-6}$ for adiabatic initial conditions and $\Omega_\textrm{gw} h^2 < 2.9 \times10^{-7}$ for homogeneous initial conditions (both $2\sigma$ upper limits). For scales close to the horizon size at recombination, we use second-order perturbation theory to calculate the back-reaction from gravitational waves, finding $\Omega_\textrm{gw} h^2 < 8.4 \times10^{-7}$, in the absence of neutrino anisotropic stress and $\Omega_\textrm{gw} h^2 < 8.6 \times10^{-7}$ when including neutrino anisotropic stress. These constraints are valid for $ 10^{-15}\, \text{Hz} \gtrsim f \gtrsim 3 \times 10^{-16}\, \text{Hz}$.

A small fraction, $f$, of cosmic string loops can collapse to form Primordial Black Holes (PBHs). Constraints on the abundance of PBHs can therefore be used to constrain $f$. We update these calculations, taking into account the PBH extended mass function, and find $f < 10^{-31} (G \mu/ c^2)^{-3/2}$. This is roughly two orders of magnitude tighter than previous constraints. The improvement from the tighter constraints on the abundance of PBHs is partly off-set by refinements to the theoretical calculation of the cosmic string loop formation rate.

We revisit the constraints on scalar tensor theories of modified gravity following the purge of GW170817. We pay particular attention to dynamical loopholes where the anomalous speed of propagation of the gravitational wave can vanish on-shell, when we impose the dynamical field equations. By working in the effective field theory formalism we are able to improve on previous analyses, scanning a much wider class of theories, including Beyond Horndeski and DHOST. Furthermore, the formalism is well adapted to consider the effect of inhomogeneous perturbations, where we explicitly take into account the fact that the galactic overdensities are pressureless to leading order.

Models in which scalar field dark energy interacts with dark matter via a pure momentum coupling have previously been found to potentially ease the structure formation tension between early- and late-universe observations. In this article we explore the physical mechanism underlying this feature. We argue analytically that the perturbation growth equations imply the suppression of structure growth, illustrating our discussion with numerical calculations. Then we generalise the previously studied quadratic coupling between the dark energy and dark matter to a more general power law case, also allowing for the slope of the dark energy exponential potential to vary. We find that the structure growth suppression is a generic feature of power law couplings and it can, for a range of parameter values, be larger than previously found.

We investigate cosmological models in which dynamical dark energy consists of a scalar field whose present-day value is controlled by a coupling to the neutrino sector. The behaviour of the scalar field depends on three functions: a kinetic function, the scalar field potential, and the scalar field-neutrino coupling function. We present an analytic treatment of the background evolution during radiation- and matter-domination for exponential and inverse power law potentials, and find a relaxation of constraints compared to previous work on the amount of early dark energy in the exponential case. We then carry out a numerical analysis of the background cosmology for both types of potential and various illustrative choices of the kinetic and coupling functions. By applying bounds from Planck on the amount of early dark energy, we are able to constrain the magnitude of the kinetic function at early times.

Scalar fields coupled to gravity through the Ricci scalar have been considered both as dark matter candidates and as a possible modified gravity explanation for galactic dynamics. It has recently been demonstrated that the dynamics of baryonic matter in disk galaxies may be explained, in the absence of particle dark matter, by a symmetron scalar field that mediates a fifth force. The symmetron provides a concrete and archetypal field theory within which to explore how large a role modifications of gravity can play on galactic scales. In this article, we extend these previous works by asking whether the same symmetron field can explain the difference between the baryonic and lens masses of galaxies through a modification of gravity. We consider the possibilities for minimal modifications of the model and find that this difference cannot be explained entirely by the symmetron fifth force without extending the field content of the model. Instead, we are pushed towards a regime of parameter space where one scalar field both mediates a fifth force and stores enough energy density that it also contributes to the galaxy's gravitational potential as a dark matter component, a regime which remains to be fully explored.

We present a class of two-field inflationary models, known as`shift-symmetric orbital inflation', whose behaviour is strongly multi-field but whose predictions are remarkably close to those of single-field inflation. In these models, the field space metric and potential are such that the inflaton trajectory is along an `angular' isometry direction whose `radius' is constant but arbitrary. As a result, the radial (isocurvature) perturbations away from the trajectory are exactly massless and they freeze on superhorizon scales. These models are the first exact realization of the `ultra-light isocurvature' scenario, previously described in the literature, where a combined shift symmetry emerges between the curvature and isocurvature perturbations and results in primordial perturbation spectra that are entirely consistent with current observations. Due to the turning trajectory, the radial perturbation sources the tangential (curvature) perturbation and makes it grow linearly in time. As a result, only one degree of freedom (i.e. the one from isocurvature modes) is responsible for the primordial observables at the end of inflation, which yields the same phenomenology as in single-field inflation. In particular, isocurvature perturbations and local non-Gaussianity are highly suppressed here, even if the inflationary dynamics is truly multi-field. We comment on the generalization to models with more than two fields.

The accelerated expansion of the universe motivates a wide class of scalar field theories that modify gravity on large scales. In regions where the weak field limit of General Relativity has been confirmed by experiment, such theories need a screening mechanism to suppress the new force. We have measured the acceleration of an atom toward a macroscopic test mass inside a high vacuum chamber, where the new force is unscreened in some theories. Our measurement, made using atom interferometry, shows that the attraction between atoms and the test mass does not differ appreciably from Newtonian gravity. This result places stringent limits on the free parameters in chameleon and symmetron theories of modified gravity.

The Nobel Prize winning confirmation in 1998 of the accelerated expansion of our Universe put into sharp focus the need of a consistent theoretical model to explain the origin of this acceleration. As a result over the past two decades there has been a huge theoretical and observational effort into improving our understanding of the Universe. The cosmological equations describing the dynamics of a homogeneous and isotropic Universe are systems of ordinary differential equations, and one of the most elegant ways these can be investigated is by casting them into the form of dynamical systems. This allows the use of powerful analytical and numerical methods to gain a quantitative understanding of the cosmological dynamics derived by the models under study. In this review we apply these techniques to cosmology. We begin with a brief introduction to dynamical systems, fixed points, linear stability theory, Lyapunov stability, centre manifold theory and more advanced topics relating to the global structure of the solutions. Using this machinery we then analyse a large number of cosmological models and show how the stability conditions allow them to be tightly constrained and even ruled out on purely theoretical grounds. We are also able to identify those models which deserve further in depth investigation through comparison with observational data. This review is a comprehensive and detailed study of dynamical systems applications to cosmological models focusing on the late-time behaviour of our Universe, and in particular on its accelerated expansion. In self contained sections we present a large number of models ranging from canonical and non-canonical scalar fields, interacting models and non-scalar field models through to modified gravity scenarios. Selected models are discussed in detail and interpreted in the context of late-time cosmology.

We study the relationship between the strength of fifth forces and the origin of scale breaking in the Standard Model (SM) of particle physics. We start with a light scalar field that is conformally coupled to a toy SM matter sector through a Weyl rescaling of the metric. After appropriately normalizing the fields, the conformally coupled scalar only interacts directly with the would-be Higgs field through kinetic-mixing and Higgs-portal terms. Thus, for the first time, we describe the equivalence of conformally coupled scalar-tensor modifications of gravity and Higgs-portal theories, and we find that the usual tree-level fifth forces only emerge if there is mass mixing between the conformally coupled scalar and the Higgs field. The strength of the fifth force, mediated by the light scalar, then depends on whether the mass of the Higgs arises from an explicit symmetry-breaking term or a spontaneous mechanism of scale breaking. Solar System tests of gravity and the non-observation of fifth forces therefore have the potential to provide information about the structure of the Higgs sector and the origin of its symmetry breaking, setting an upper bound on the magnitude of any explicit scale-breaking terms. These results demonstrate the phenomenological importance (both for cosmology and high-energy physics) of considering how scalar-tensor modifications of gravity are embedded within extensions of the SM.

It has been recently suggested that oscillons produced in the early universe from certain asymmetric potentials continue to emit gravitational waves for a number of $e$-folds of expansion after their formation, leading to potentially detectable gravitational wave signals. We revisit this claim by conducting a convergence study using graphics processing unit (GPU)-accelerated lattice simulations and show that numerical errors accumulated with time are significant in low-resolution scenarios, or in scenarios where the run-time causes the resolution to drop below the relevant scales in the problem. Our study determines that the dominant, growing high frequency peak of the gravitational wave signals in the fiducial "hill-top model" in [arXiv:1607.01314] is a numerical artifact. This finding prompts the need for a more careful analysis of the numerical validity of other similar results related to gravitational waves from oscillon dynamics.

Chameleon scalar fields can screen their associated fifth forces from detection by changing their mass with the local density. These models are an archetypal example of a screening mechanism, and have become an important target for both cosmological surveys and terrestrial experiments. In particular there has been much recent interest in searching for chameleon fifth forces in the laboratory. It is known that the chameleon force is less screened around non-spherical sources, but only the field profiles around a few simple shapes are known analytically. In this work we introduce a numerical code that solves for the chameleon field around arbitrary shapes with azimuthal symmetry placed in a spherical vacuum chamber. We find that deviations from spherical symmetry can increase the chameleon acceleration experienced by a test particle by up to a factor of $\sim 3$, and that the least screened objects are those which minimize some internal dimension.

We show that the radial acceleration relation for rotationally-supported galaxies may be explained, in the absence of cold dark matter, by a non-minimally coupled scalar field, whose fifth forces are partially screened on galactic scales by the symmetron mechanism. In addition, we show that sufficient energy is stored in the symmetron field to explain the dynamic stability of galactic disks.

We describe a symmetron model in which the screening of fifth forces arises at the one-loop level through the Coleman-Weinberg mechanism of spontaneous symmetry breaking. We show that such a theory can avoid current constraints on the existence of fifth forces but still has the potential to give rise to observable deviations from general relativity, which could be seen in cold atom experiments.

We present the first complete Markov chain Monte Carlo analysis of cosmological models with evolving cosmic (super)string networks, using the unconnected segment model in the unequal-time correlator formalism. For ordinary cosmic string networks, we derive joint constraints on Lambda cold dark matter (CDM) and string network parameters, namely the string tension Gmu, the loop-chopping efficiency c_r and the string wiggliness \alpha. For cosmic superstrings, we obtain joint constraints on the fundamental string tension Gmu_F, the string coupling g_s, the self-interaction coefficient c_s, and the volume of compact extra dimensions w. This constitutes the most comprehensive CMB analysis of LambdaCDM cosmology + strings to date. For ordinary cosmic string networks our updated constraint on the string tension is, in relativistic units, Gmu<1.1x10^-7, while for cosmic superstrings our constraint on the fundamental string tension is Gmu_F<2.8x10^-8, both obtained using Planck2015 temperature and polarisation data.

Light scalar fields coupled to matter are a common consequence of theories of dark energy and attempts to solve the cosmological constant problem. The chameleon screening mechanism is commonly invoked in order to suppress the fifth forces mediated by these scalars, suficiently to avoid current experimental constraints, without fine tuning. The force is suppressed dynamically by allowing the mass of the scalar to vary with the local density. Recently it has been shown that near future cold atoms experiments using atom-interferometry have the ability to access a large proportion of the chameleon parameter space. In this work we demonstrate how experiments utilising asymmetric parallel plates can push deeper into the remaining parameter space available to the chameleon.

We review the tantalising prospect that the first evidence for the dark energy driving the observed acceleration of the Universe on giga-parsec scales may be found through metre scale laboratory based atom interferometry experiments. To do that, we first introduce the idea that scalar fields could be responsible for dark energy and show that in order to be compatible with fifth force constraints these fields must have a screening mechanism which hides their effects from us within the solar system. Particular emphasis is placed on one such screening mechanism known as the chameleon effect where the field's mass becomes dependent on the environment. The way the field behaves in the presence of a spherical source is determined and we then go on to show how in the presence of the kind of high vacuum associated with atom interferometry experiments, and when the test particle is an atom, it is possible to use the associated interference pattern to place constraints on the acceleration due to the fifth force of the chameleon field - this has already been used to rule out large regions of the chameleon parameter space and maybe one day will be able to detect the force due to the dark energy field in the laboratory.

We revisit Starobinsky inflation in a quantum gravitational context, by means of the exact Renormalisation Group (RG). We calculate the non-perturbative beta functions for Newton's `constant' G and the dimensionless R^2 coupling, and show that there exists an attractive UV fixed point where the latter one vanishes but not the former one, and we provide the corresponding beta functions. The smallness of the R^2 coupling, required for agreement with inflationary observables, is naturally ensured by its vanishing at the UV fixed point, ensuring the smallness of the primordial fluctuations, as well as providing a theoretical motivation for the initial conditions needed for successful inflation in this context. We discuss the corresponding RG dynamics, showing both how inflationary and classical observations define the renormalisation conditions for the couplings, and also how the UV regime is connected with lower energies along the RG flow. Finally, we discuss the consistency of our results when higher order curvature corrections are included, and show that they are robust to the inclusion of R^3 corrections.

We present the most general parametrisation of models of dark energy in the form of a scalar field which is explicitly coupled to dark matter. We follow and extend the Parameterized Post-Friedmannian approach, previously applied to modified gravity theories, in order to include interacting dark energy. We demonstrate its use through a number of worked examples and show how the initially large parameter space of free functions can be significantly reduced and constrained to include only a few non-zero coefficients. This paves the way for a model-independent approach to classify and test interacting dark energy theories.

Theories of dark energy require a screening mechanism to explain why the associated scalar fields do not mediate observable long range fifth forces. The archetype of this is the chameleon field. Here we show that individual atoms are too small to screen the chameleon field inside a large high-vacuum chamber, and therefore can detect the field with high sensitivity. We derive new limits on the chameleon parameters from existing experiments, and show that most of the remaining chameleon parameter space is readily accessible using atom interferometry.

The chameleon mechanism enables a long range fifth force to be screened in dense environments when non-trivial self interactions of the field cause its mass to increase with the local density. To date, chameleon fifth forces have mainly been studied for spherically symmetric sources, however the non-linear self interactions mean that the chameleon responds to changes in the shape of the source differently to gravity. In this work we focus on ellipsoidal departures from spherical symmetry and compute the full form of the chameleon force, comparing it's shape dependence to that of gravity. Enhancement of the chameleon force by up to 40% is possible when deforming a sphere to an ellipsoid of the same mass, with an ellipticity $\simeq 0.99$.

Q-balls are non-topological solitonic solutions to a wide class of field theories that possess global symmetries. Here we show that in these same theories there also exists a tower of novel composite Q-ball solutions where, within one composite Q-ball, positive and negative charges co-exist and swap at a frequency lower than the natural frequency of an individual Q-ball. These charge-swapping Q-balls are constructed by assembling Q-balls and anti-Q-balls tightly such that their nonlinear cores overlap. We explain why charge-swapping Q-balls can form and why they swap charges.

We have recently proposed a novel self tuning mechanism to alleviate the famous cosmological constant problem, based on the general scalar tensor theory proposed by Horndeski. The self-tuning model ends up consisting of four geometric terms in the action, with each term containing a free potential function of the scalar field; the four together being labeled as the Fab-Four. In this paper we begin the important task of deriving the cosmology associated with the Fab-Four Lagrangian. Performing a phase plane analysis of the system we are able to obtain a number of fixed points for the system, with some remarkable new solutions emerging from the trade-off between the various potentials. As well as obtaining inflationary solutions we also find conventional radiation/matter-like solutions, but in regimes where the energy density is dominated by a cosmological constant, and where we do not have any explicit forms of radiation or matter. Stability conditions for matter solutions are obtained and we show how it is possible for there to exist an extended period of `matter domination' opening up the possibility that we can generate cosmological structures, and recover a consistent cosmology even in the presence of a large cosmological constant.

We propose a completely new parametrisation of the dark energy equation of state, which uses the dark energy density, $\Omega_e$ as a cosmic clock. We expand the equation of state in a series of orthogonal polynomials, with $\Omega_e$ as the expansion parameter and determine the expansion coefficients by fitting to SNIa and $H(z)$ data. Assuming that $\Omega_e$ is a monotonic function of time, we show that our parametrisation performs better than the popular Chevallier--Polarski--Linder (CPL) and Gerke and Efstathiou (GE) parametrisations, and we demonstrate that it is robust to the choice of prior. Expanding in orthogonal polynomials allows us to relate models of dark energy directly to our parametrisation, which we illustrate by placing constraints on the expansion coefficients extracted from two popular quintessence models. Finally, we comment on how this parametrisation could be modified to accommodate high redshift data, where any non--monotonicity of $\Omega_e$ would need to be accounted for.

We present three distinct types of models of dark energy in the form of a scalar field which is explicitly coupled to dark matter. Our construction draws from the pull-back formalism for fluids and generalises the fluid action to involve couplings to the scalar field. We investigate the cosmology of each class of model both at the background and linearly perturbed level. We choose a potential for the scalar field and a specific coupling function for each class of models and we compute the Cosmic Microwave Background and matter power spectra.

Oscillons are long-lived, localized excitations of nonlinear scalar fields which may be copiously produced during preheating after inflation, leading to a possible oscillon-dominated phase in the early Universe. For example, this can happen after axion monodromy inflation, on which we run our simulations. We investigate the stochastic gravitational wave background associated with an oscillon-dominated phase. An isolated oscillon is spherically symmetric and does not radiate gravitational waves, and we show that the flux of gravitational radiation generated between oscillons is also small. However, a significant stochastic gravitational wave background may be generated during preheating itself (i.e, when oscillons are forming), and in this case the characteristic size of the oscillons is imprinted on the gravitational wave power spectrum, which has multiple, distinct peaks.

We study the evolution of the non-linear curvature perturbation during perturbative reheating, and hence how observables evolve to their final values which we may compare against observations. Our study includes the evolution of the two trispectrum parameters, $\gnl$ and $\taunl$, as well as the scale dependence of both $\fnl$ and $\taunl$. In general the evolution is significant and must be taken into account, which means that models of multifield inflation cannot be compared to observations without specifying how the subsequent reheating takes place. If the trispectrum is large at the end of inflation, it normally remains large at the end of reheating. In the classes of models we study, it is very hard to generate $\taunl\gg\fnl^2$, regardless of the decay rates of the fields. Similarly, for the classes of models in which $\gnl\simeq\taunl$ during slow--roll inflation, we find the relation typically remains valid during reheating. Therefore it is possible to observationally test such classes of models without specifying the parameters of reheating, even though the individual observables are sensitive to the details of reheating. It is hard to generate an observably large $\gnl$ however. The runnings, $\nfnl$ and $\ntaunl$, tend to satisfy a consistency relation $\ntaunl=(3/2)\nfnl$, but are in general too small to be observed for the class of models considered regardless of reheating timescale.

Primordial Black Holes (PBH's) can form in the early Universe from the collapse of large density fluctuations. Tight observational limits on their abundance constrain the amplitude of the primordial fluctuations on very small scales which can not otherwise be constrained, with PBH's only forming from the extremely rare large fluctuations. The number of PBH's formed is therefore sensitive to small changes in the shape of the tail of the fluctuation distribution, which itself depends on the amount of non-Gaussianity present. We study, for the first time, how quadratic and cubic local non-Gaussianity of arbitrary size (parameterised by f_nl and g_nl respectively) affects the PBH abundance and the resulting constraints on the amplitude of the fluctuations on very small scales. Intriguingly we find that even non-linearity parameters of order unity have a significant impact on the PBH abundance. The sign of the non-Gaussianity is particularly important, with the constraint on the allowed fluctuation amplitude tightening by an order of magnitude as f_nl changes from just -0.5 to 0.5. We find that if PBH's are observed in the future, then regardless of the amplitude of the fluctuations, non-negligible negative f_nl would be ruled out. Finally we show that g_nl can have an even larger effect on the number of PBH's formed than f_nl.

We present analytic expressions for the cosmic string unequal time correlator (UETC) in the context of the Unconnected Segment Model (USM). This eliminates the need to simulate the many thousands of network realisations needed to estimate the UETC numerically. With our approach we can compute the UETC very accurately, over all scales of interest, in the order of ~20-30 seconds on a single CPU. Our formalism facilitates an efficient approach to computing Cosmic Microwave Background (CMB) anisotropies for cosmic strings. Discretising the UETC and performing an eigen-decomposition to act as sources in the CAMB CMB code, the power spectrum can be calculated by summing over a finite number of eigenmodes. A much smaller number of eigenmodes are required compared to the conventional approach of averaging power spectra over a finite number of realisations of the string network. With the additional efficiency and performance improvements offered by the OpenMP CAMB code, the time required to compute string power spectra is significantly reduced compared to the standard serial CMBACT code. The latter takes ~30 hours on a modern single threaded CPU for 2000 network realisations. Similar percent level accuracy can be achieved with our approach on a moderately threaded CPU (8 threads) in only ~15 minutes. If accuracy is only required at the 10 percent level and the CPU is more highly threaded, cosmic string power spectra are now possible in ~2-3 minutes. This makes exploration of the string parameter space now possible for Markov-Chain-Monte-Carlo (MCMC) analysis.

This paper is concerned with string cosmology and the dynamics of multiple scalar fields in potentials that can become negative, and their features as (Early) Dark Energy models. Our point of departure is the "String Axiverse", a scenario that motivates the existence of cosmologically light axion fields as a generic consequence of string theory. We couple such an axion to its corresponding modulus. We give a detailed presentation of the rich cosmology of such a model, ranging from the setting of initial conditions on the fields during inflation, to the asymptotic future. We present some simplifying assumptions based on the fixing of the axion decay constant $f_a$, and on the effective field theory when the modulus trajectory is adiabatic, and find the conditions under which these assumptions break down. As a by-product of our analysis, we find that relaxing the assumption of fixed $f_a$ leads to the appearance of a new meta-stable de-Sitter region for the modulus without the need for uplifting by an additional constant. A dynamical systems analysis reveals the existence of many fixed point attractors, repellers and saddle points, which we analyse in detail. We also provide geometric interpretations of the phase space. The fixed points can be used to bound the couplings in the model. A systematic scan of certain regions of parameter space reveals that the future evolution of the universe in this model can be rich, containing multiple epochs of accelerated expansion.

We study the effects of perturbative reheating on the evolution of the curvature perturbation \zeta, in two-field inflation models. We use numerical methods to explore the sensitivity of f_NL, n_s and r to the reheating process, and present simple qualitative arguments to explain our results. In general, if a large non-Gaussian signal exists at the start of reheating, it will remain non zero at the end of reheating. Unless all isocurvature modes have completely decayed before the start of reheating, we find that the non-linearity parameter, f_NL, can be sensitive to the reheating timescale, and that this dependence is most appreciable for `runaway' inflationary potentials that only have a minimum in one direction. For potentials with a minimum in both directions, f_NL can also be sensitive to reheating if a mild hierarchy exists between the decay rates of each field. Within the class of models studied, we find that the spectral index n_s, is fairly insensitive to large changes in the field decay rates, indicating that n_s is a more robust inflationary observable, unlike the non-linearity parameter f_NL. Our results imply that the statistics of \zeta, especially f_NL, can only be reliably used to discriminate between models of two-field inflation if the physics of reheating are properly accounted for.

Galileon models arise in certain braneworld scenarios as modifications to General Relativity, and are also interesting as field theories in their own right. We show how the galileon model can be naturally generalized to include local gauge symmetries, by allowing for couplings to Yang-Mills fields. The resulting theories have at most second order spacetime derivatives in any representation of the gauge group, thereby avoiding Ostrogradski ghosts. We also extend the models to include curved space, and show how in that case we need to include non-minimal couplings between the galileons and the curvature tensors for the theory to retain its second order nature.

We have recently proposed a special class of scalar tensor theories known as the Fab Four. These arose from attempts to analyse the cosmological constant problem within the context of Horndeski's most general scalar tensor theory. The Fab Four together give rise to a model of self-tuning, with the relevant solutions evading Weinberg's no-go theorem by relaxing the condition of Poincare invariance in the scalar sector. The Fab Four are made up of four geometric terms in the action with each term containing a free potential function of the scalar field. In this paper we rigorously derive this model from the general model of Horndeski, proving that the Fab Four represents the only classical scalar tensor theory of this type that has any hope of tackling the cosmological constant problem. We present the full equations of motion for this theory, and give an heuristic argument to suggest that one might be able to keep radiative corrections under control. We also give the Fab Four in terms of the potentials presented in Deffayet et al's version of Horndeski.

A sufficiently light scalar field slowly evolving in a potential can account for the dark energy that presently dominates the universe. This quintessence field is expected to couple directly to matter components, unless some symmetry of a more fundamental theory protects or suppresses it. Such a coupling would leave distinctive signatures in the background expansion history of the universe and on cosmic structure formation, particularly at galaxy cluster scales. Using semi--analytic expressions for the CDM halo mass function, we make predictions for halo abundance in models where the quintessence scalar field is coupled to cold dark matter, for a variety of quintessence potentials. We evaluate the linearly extrapolated density contrast at the redshift of collapse using the spherical collapse model and we compare this result to the corresponding prediction obtained from the non--linear perturbation equations in the Newtonian limit. For all the models considered in this work, if there is a continuous flow of energy from the quintessence scalar field to the CDM component, then the predicted number of CDM haloes can only lie below that of $\Lambda$CDM, when each model shares the same cosmological parameters today. In the last stage of our analysis we perform a global MCMC fit to data to find the best fit values for the cosmological model parameters. We find that for some forms of the quintessence potential, coupled dark energy models can offer a viable alternative to $\Lambda$CDM in light of the recent detections of massive high--$z$ galaxy clusters, while other models of coupled quintessence predict a smaller number of massive clusters at high redshift compared to $\Lambda$CDM.

We study signatures of cosmic superstring networks containing strings of multiple tensions and Y-junctions, on the cosmic microwave background (CMB) temperature and polarisation spectra. Focusing on the crucial role of the string coupling constant $g_s$, we show that the number density and energy density of the scaling network are dominated by different types of string in the $g_s \sim 1$ and $g_s \ll 1$ limits. This can lead to an observable shift in the position of the B-mode peak --- a distinct signal leading to a direct constraint on $g_s$. We forecast the joint bounds on $g_s$ and the fundamental string tension $\mu_F$ from upcoming and future CMB polarisation experiments, as well as the signal to noise in detecting the difference between B-mode signals in the limiting cases of large and small $g_s$. We show that such a detectable shift is within reach of planned experiments.

We review the existence, formation and properties of cosmic strings in string theory, the wide variety of observational techniques that are being employed to detect them, and the constraints that current observations impose on string theory models.

Starting from the most general scalar-tensor theory with second order field equations in four dimensions, we establish the unique action that will allow for the existence of a consistent self-tuning mechanism on FLRW backgrounds, and show how it can be understood as a combination of just four base Lagrangians with an intriguing geometric structure dependent on the Ricci scalar, the Einstein tensor, the double dual of the Riemann tensor and the Gauss-Bonnet combination. Spacetime curvature can be screened from the net cosmological constant at any given moment because we allow the scalar field to break Poincar\'e invariance on the self-tuning vacua, thereby evading the Weinberg no-go theorem. We show how the four arbitrary functions of the scalar field combine in an elegant way opening up the possibility of obtaining non-trivial cosmological solutions.

We study the cosmic microwave background temperature and polarisation spectra sourced by multi-tension cosmic superstring networks. First we obtain solutions for the characteristic length scales and velocities associated with the evolution of a network of F-D strings, allowing for the formation of junctions between strings of different tensions. We find two distinct regimes describing the resulting scaling distributions for the relative densities of the different types of strings, depending on the magnitude of the fundamental string coupling g_s. In one of them, corresponding to the value of the coupling being of order unity, the network's stress-energy power spectrum is dominated by populous light F and D strings, while the other regime, at smaller values of g_s, has the spectrum dominated by rare heavy D strings. These regimes are seen in the CMB anisotropies associated with the network. We focus on the dependence of the shape of the B-mode polarisation spectrum on g_s and show that measuring the peak position of the B-mode spectrum can point to a particular value of the string coupling. Finally, we assess how this result, along with pulsar bounds on the production of gravitational waves from strings, can be used to constrain a combination of g_s and the fundamental string tension mu_F. Since CMB and pulsar bounds constrain different combinations of the string tensions and densities, they result in distinct shapes of bounding contours in the (mu_F, g_s) parameter plane, thus providing complementary constraints on the properties of cosmic superstrings.

Modified gravity theories with the Gauss-Bonnet term $G=R^2-4R^{\mu\nu}R_{\mu\nu}+R^{\mu\nu\rho\sigma}R_{\mu\nu\rho\sigma}$ have recently gained a lot of attention as a possible explanation of dark energy. We perform a thorough phase space analysis on the so-called $f(G)$ models, where $f(G)$ is some general function of the Gauss-Bonnet term, and derive conditions for the cosmological viability of $f(G)$ dark energy models. Following the $f(R)$ case, we show that these conditions can be nicely presented as geometrical constraints on the derivatives of $f(G)$. We find that for general $f(G)$ models there are two kinds of stable accelerated solutions, a de Sitter solution and a phantom-like solution. They co-exist with each other and which solution the universe evolves to depends on the initial conditions. Finally, several toy models of $f(G)$ dark energy are explored. Cosmologically viable trajectories that mimic the $\Lambda$CDM model in the radiation and matter dominated periods, but have distinctive signatures at late times, are obtained.

We analyze the dynamics of a Dirac-Born-Infeld (DBI) field in a cosmological set-up which includes a perfect fluid. Introducing convenient dynamical variables, we show the evolution equations form an autonomous system when the potential and the brane tension of the DBI field are arbitrary power-law or exponential functions of the DBI field. In particular we find scaling solutions can exist when powers of the field in the potential and warp-factor satisfy specific relations. A new class of fixed-point solutions are obtained corresponding to points which initially appear singular in the evolution equations, but on closer inspection are actually well defined. In all cases, we perform a phase-space analysis and obtain the late-time attractor structure of the system. Of particular note when considering cosmological perturbations in DBI inflation is a fixed-point solution where the Lorentz factor is a finite large constant and the equation of state parameter of the DBI field is $w=-1$. Since in this case the speed of sound $c_s$ becomes constant, the solution can be thought to serve as a good background to perturb about.

We consider the evolution of a network of strings in an expanding universe, allowing for the formation of junctions between strings of different tensions. By explicitly including, in the velocity-dependent evolution equations for the network, kinematic constraints associated with the formation of Y-shaped string junctions, we show how they lead to scaling solutions in regimes where they would not otherwise be found, thereby extending the range of parameters which lead to scaling. By incorporating these constraints we are able to study their general behaviour for networks with cosmic superstring interaction rules, and predict the scaling densities expected by these networks.

Cosmic strings are predicted by many field-theory models, and may have been formed at a symmetry-breaking transition early in the history of the universe, such as that associated with grand unification. They could have important cosmological effects. Scenarios suggested by fundamental string theory or M-theory, in particular the popular idea of brane inflation, also strongly suggest the appearance of similar structures. Here we review the reasons for postulating the existence of cosmic strings or superstrings, the various possible ways in which they might be detected observationally, and the special features that might discriminate between ordinary cosmic strings and superstrings.

We discuss some hitherto puzzling features of the small-scale structure of cosmic strings. We argue that kinks play a key role, and that an important quantity to study is their sharpness distribution. In particular we suggest that for very small scales the two-point correlation function of the string tangent vector varies linearly with the separation and not as a fractional power, as proposed by Polchinski and Rocha [Phys. Rev. D 74, 083504 (2006)]. However, our results are consistent with theirs, because the range of scales to which this linearity applies shrinks as evolution proceeds.

We study the classical and absolute stability of Q-balls in scalar field theories with flat potentials arising in both gravity-mediated and gauge-mediated models. We show that the associated Q-matter formed in gravity-mediated potentials can be stable against decay into their own free particles as long as the coupling constant of the nonrenormalisable term is small, and that all of the possible three-dimensional Q-ball configurations are classically stable against linear fluctuations. Three-dimensional gauge-mediated Q-balls can be absolutely stable in the "thin-wall-like" limit, but are completely unstable in the "thick-wall" limit.

We analyze the dynamics of a single scalar field in Friedmann-Robertson-Walker universes with spatial curvature. We obtain the fixed point solutions which are shown to be late time attractors. In particular, we determine the corresponding scalar field potentials which correspond to these stable solutions. The analysis is quite general and incorporates expanding and contracting universes with both positive and negative scalar potentials. We demonstrate that the known power law, exponential, and de-Sitter solutions are certain limits of our general set of solutions.

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

Superstring theory predicts the potential formation of string networks with bound states ending in junctions. Kinematic constraints for junction formation have been derived within the Nambu-Goto thin string approximation. Here we test these constraints numerically in the framework of the Abelian-Higgs model in the Type-I regime and report on good agreement with the analytical predictions. We also demonstrate that strings can effectively pass through each other when they meet at speeds slightly above the critical velocity permitting bound state formation. This is due to reconnection effects that are beyond the scope of the Nambu-Goto approximation.

We study the formation of three-string junctions between (p,q)-cosmic superstrings, and collisions between such strings and show that kinematic constraints analogous to those found previously for collisions of Nambu-Goto strings apply here too, with suitable modifications to take account of the additional requirements of flux conservation. We examine in detail several examples involving collisions between strings with low values of p and q, and also examine the rates of growth or shrinkage of strings at a junction. Finally, we briefly discuss the formation of junctions for strings in a warped space, specifically with a Klebanov-Strassler throat, and show that similar constraints still apply with changes to the parameters taking account of the warping and the background flux.

We derive a new class of exact time dependent solutions in a warped six dimensional supergravity model. Under the assumptions we make for the form of the underlying moduli fields, we show that the only consistent time dependent solutions lead to all six dimensions evolving in time, implying the eventual decompactification or collapse of the extra dimensions. We also show how the dynamics affects the quantization of the deficit angle.

We consider the constraints on string networks with junctions in which the strings may all be different, as may be found for example in a network of $(p,q)$ cosmic superstrings. We concentrate on three aspects of junction dynamics. First we consider the propagation of small amplitude waves across a static three-string junction. Then, generalizing our earlier work, we determine the kinematic constraints on two colliding strings with different tensions. As before, the important conclusion is that strings do not always reconnect with a third string; they can pass straight through one another (or in the case of non-abelian strings become stuck in an X configuration), the constraint depending on the angle at which the strings meet, on their relative velocity, and on the ratios of the string tensions. For example, if the two colliding strings have equal tensions, then for ultra-relativistic initial velocities they pass through one another. However, if their tensions are sufficiently different they can reconnect. Finally, we consider the global properties of junctions and strings in a network. Assuming that, in a network, the incoming waves at a junction are independently randomly distributed, we determine the r.m.s. velocities of strings and calculate the average speed at which a junction moves along each of the three strings from which it is formed. Our findings suggest that junction dynamics may be such as to preferentially remove the heavy strings from the network leaving a network of predominantly light strings. Furthermore the r.m.s. velocity of strings in a network with junctions is smaller than 1/\sqrt{2}, the result for conventional Nambu-Goto strings without junctions in Minkowski spacetime.

We investigate the evolution of scalar metric perturbations across a sudden cosmological transition, allowing for an inhomogeneous surface stress at the transition leading to a discontinuity in the local expansion rate, such as might be expected in a big crunch/big bang event. We assume that the transition occurs when some function of local matter variables reaches a critical value, and that the surface stress is also a function of local matter variables. In particular we consider the case of a single scalar field and show that a necessary condition for the surface stress tensor to be perturbed at the transition is the presence of a non-zero intrinsic entropy perturbation of the scalar field. We present the matching conditions in terms of gauge-invariant variables assuming a sudden transition to a fluid-dominated universe with barotropic equation of state. For adiabatic perturbations the comoving curvature perturbation is continuous at the transition, while the Newtonian potential may be discontinuous if there is a discontinuity in the background Hubble expansion.

We study the dynamics of Nambu--Goto strings with junctions at which three strings meet. In particular, we exhibit one simple exact solution and examine the process of intercommuting of two straight strings, in which they exchange partners but become joined by a third string. We show that there are important kinematical constraints on this process. The exchange cannot occur if the strings meet with very large relative velocity. This may have important implications for the evolution of cosmic superstring networks and non-abelian string networks.

In this paper we review in detail a number of approaches that have been adopted to try and explain the remarkable observation of our accelerating Universe. In particular we discuss the arguments for and recent progress made towards understanding the nature of dark energy. We review the observational evidence for the current accelerated expansion of the universe and present a number of dark energy models in addition to the conventional cosmological constant, paying particular attention to scalar field models such as quintessence, K-essence, tachyon, phantom and dilatonic models. The importance of cosmological scaling solutions is emphasized when studying the dynamical system of scalar fields including coupled dark energy. We study the evolution of cosmological perturbations allowing us to confront them with the observation of the Cosmic Microwave Background and Large Scale Structure and demonstrate how it is possible in principle to reconstruct the equation of state of dark energy by also using Supernovae Ia observational data. We also discuss in detail the nature of tracking solutions in cosmology, particle physics and braneworld models of dark energy, the nature of possible future singularities, the effect of higher order curvature terms to avoid a Big Rip singularity, and approaches to modifying gravity which leads to a late-time accelerated expansion without recourse to a new form of dark energy.

Braneworld scenarios are motivated by string/M-theory and can be characterized by the way in which they modify the conventional Friedmann equations of Einstein gravity. An alternative approach to quantum gravity, however, is the loop quantum cosmology program. In the semi-classical limit, the cosmic dynamics in this scenario can also be described by a set of modified Friedmann equations. We demonstrate that a dynamical correspondence can be established between these two paradigms at the level of the effective field equations. This allows qualitatively similar features between the two approaches to be compared and contrasted as well as providing a framework for viewing braneworld scenarios in terms of constrained Hamiltonian systems. As concrete examples of this correspondence, we illustrate the relationships between different cosmological backgrounds representing scaling solutions.

We consider the constraints that can be imposed on a wide class of Inflation models in modified gravity scenarios in which the Friedmann equation is modified by the inclusion of $\rho^2$ terms, where $\rho$ is the total energy density. In particular we obtain the reheating temperature and gravitino abundance associated with the end of inflation. Whereas models of chaotic inflation and natural inflation can easily avoid the conventional gravitino overproduction problem, we show that supersymmetric hybrid inflation models (driven by both F and D-terms) do not work in the $\rho^2$ dominated era. We also study inflation driven by exponetial potentials in this modified background, and show that the gravitino production is suppressed enough to avoid there being a problem, although other conditions severely constrain these models.

We model the behaviour of a network of interacting (p,q) strings from IIB string theory by considering a field theory containing multiple species of string, allowing us to study the effect of non-intercommuting events due to two different species crossing each other. This then has the potential for a string dominated Universe with the network becoming so tangled that it freezes. We give numerical evidence, explained by a one-scale model, that such freezing does not take place, with the network reaching a scaling limit where its density relative to the background increases with N, the number of string types.

We investigate the end of the inflationary period in the recently proposed scenario of locked inflation, and consider various constraints arising from density perturbations, loop corrections, parametric resonance and defect formation. We show that in a scenario where there is one long period of locked inflation, it is not possible to satisfy all of these constraints without having a period of saddle inflation afterwards, which would wipe out all observable signatures. On the other hand, if one does not insist on satisfying the loop correction constraint, saddle inflation can be avoided, but then inflation must have ended through parametric resonance.

We study a dark energy scenario in the presence of a tachyon field $\phi$ with potential $V(\phi)$ and a barotropic perfect fluid. The cosmological dynamics crucially depends on the asymptotic behavior of the quantity $\lambda=-M_pV_\phi/V^{3/2}$. If $\lambda$ is a constant, which corresponds to an inverse square potential $V(\phi) \propto \phi^{-2}$, there exists one stable critical point that gives an acceleration of the universe at late times. When $\lambda \to 0$ asymptotically, we can have a viable dark energy scenario in which the system approaches an ``instantaneous'' critical point that dynamically changes with $\lambda$. If $|\lambda|$ approaches infinity asymptotically, the universe does not exhibit an acceleration at late times. In this case, however, we find an interesting possibility that a transient acceleration occurs in a regime where $|\lambda|$ is smaller than of order unity.

Motivated by the recent interest in cosmologies arising from energy density modifications to the Friedmann equation, we analyse the scaling behaviour for a broad class of these cosmologies comprised of scalar fields and background barotropic fluid sources. In particular, we determine the corresponding scalar field potentials which lead to attractor scaling solutions in a wide class of braneworld and dark energy scenarios. We show how a number of recent proposals for modifying the Friedmann equation can be thought of as being dual to one another, and determine the conditions under which such dualities arise.

Detecting dark energy dynamics is the main quest of current dark energy research. Addressing the issue demands a fully consistent analysis of CMB, large scale structure and SN-Ia data with multi-parameter freedom valid for all redshifts. Here we undertake a ten parameter analysis of general dark energy confronted with the first year WMAP, 2dF galaxy survey and latest SN-Ia data. Despite the huge freedom in dark energy dynamics there are no new degeneracies with standard cosmic parameters apart from a mild degeneracy between reionisation and the redshift of acceleration, both of which effectively suppress small scale power. Breaking this degeneracy will help significantly in detecting dynamics, if it exists. Our best-fit model to the data has significant late-time evolution at z<1.5. Phantom models are also considered and we find that the best-fit crosses w=-1 which, if confirmed, would be a clear signal for radically new physics. Treatment of such rapidly varying models requires careful integration of the dark energy density usually not implemented in standard codes, leading to crucial errors of up to 5%. Nevertheless cosmic variance means that standard LCDM models are still a very good fit to the data and evidence for dynamics is currently very weak. Independent tests of reionisation or the epoch of acceleration (e.g. ISW-LSS correlations) or reduction of cosmic variance at large scales (e.g. cluster polarisation at high redshift) may prove key in the hunt for dynamics.

By combining the recent WMAP measurements of the cosmic microwave background anisotropies and the results of the recent luminosity distance measurements to type-Ia supernovae, we find that the normalization of the matter power spectrum on cluster scales, sigma_8, can be used to discriminate between dynamical models of dark energy (quintessence models) and a conventional cosmological constant model (LCDM).

Dark energy affects the CMB through its perturbations and affects both CMB and Sn-Ia through its background evolution. Using recent CMB and Sn-Ia data sets, together with the most general parameterization of the dark energy equation of state available, we find that today w < -0.8 (2 sigma). We also find that the value of the normalization of the power spectrum on cluster scales, sigma8, can be used to discriminate between dynamical models of dark energy (Quintessence models) and a cosmological constant model (Lambda CDM).

We study the cosmology of a general scalar field and barotropic fluid during the early stage of a brane-world where the Friedmann constraint is dominated by the square of the energy density. Assuming both the scalar field and fluid are confined to the brane, we find a range of behaviour depending on the form of the potential. Generalising an approach developed for a standard Friedmann cosmology, in \cite{delaMacorra:1999ff}, we show that the potential dependence $V(\phi)$ can be described through a parameter $\lambda \equiv -\sqrt{2} m_5^{3/2} V'/(\sqrt{H}V)$, where $m_5$ is the 5-dimensional Planck mass, $H$ is the Hubble parameter and $V' \equiv \frac{dV}{d\phi}$. For the case where $\lambda$ asymptotes to zero, we show that the solution exhibits stable inflationary behaviour. On the other hand if it approaches a finite constant, then $V(\phi) \propto \frac{1}{\phi^2}$. For $\lambda \to \infty$ asymptotically, we find examples where it does so both with and without oscillating. In the latter case, the barotropic fluid dominates the scalar filed asymptotically. Finally we point out an interesting duality which leads to identical evolution equations in the high energy $\rho^2$ dominated regime and the low energy $\rho$ dominated regime.

Macroscopic fundamental and Dirichlet strings have several potential instabilities: breakage, tachyon decays, and confinement by axion domain walls. We investigate the conditions under which metastable strings can exist, and we find that such strings are present in many models. There are various possibilities, the most notable being a network of (p,q) strings. Cosmic strings give a potentially large window into string physics.

We explore the types of slow-roll inflationary potentials that result in scalar perturbations with a constant spectral index, i.e., perturbations that may be described by a single power-law spectrum over all observable scales. We devote particular attention to the type of potentials that result in the Harrison--Zel'dovich spectrum.

We study the change of the effective fine structure constant in the cosmological models of a scalar field with a non-vanishing coupling to the electromagnetic field. Combining cosmological data and terrestrial observations we place empirical constraints on the size of the possible coupling and explore a large class of models that exhibit tracking behavior. The change of the fine structure constant implied by the quasar absorption spectra together with the requirement of tracking behavior impose a lower bound of the size of this coupling. Furthermore, the transition to the quintessence regime implies a narrow window for this coupling around $10^{-5}$ in units of the inverse Planck mass. We also propose a non-minimal coupling between electromagnetism and quintessence which has the effect of leading only to changes of alpha determined from atomic physics phenomena, but leaving no observable consequences through nuclear physics effects. In doing so we are able to reconcile the claimed cosmological evidence for a changing fine structure constant with the tight constraints emerging from the Oklo natural nuclear reactor.

We discuss the onset of the non-linear regime in the context of unified dark matter models involving a generalised Chaplygin gas. We show that the transition from dark matter-like to dark energy-like behaviour will never be smooth. In some regions of space the transition will never take place while in others it may happen sooner or later than naively expected. As a result the linear theory used in previous studies may break down late in the matter dominated era even on large cosmological scales. We study the importance of this effect showing that its magnitude depends on the exact form of the equation of state in the low density regime. We expect that our results will be relevant for other unified dark matter scenarios particularly those where the quartessence candidate is a perfect fluid.

K-essence has been proposed as a possible means of explaining the coincidence problem of the Universe beginning to accelerate only at the present epoch. We carry out a comprehensive dynamical systems analysis of the k-essence models given so far in the literature. We numerically study the basin of attraction of the tracker solutions and we highlight the behaviour of the field close to sound speed divergences. We find that, when written in terms of parameters with a simple dynamical interpretation, the basins of attraction represent only a small region of the phase space.

We analyse the mechanism of quintessential inflation in brane world scenarios for a number of particle physics inspired scalar potentials. We constrain the parameter space of those scalar potentials and comment on the likelihood that we could discriminate these models from standard inflation, based on upcoming large scale structure and cosmic microwave background observations.

K-essence models, relying on scalar fields with non-canonical kinetic terms, have been proposed as an alternative to quintessence in explaining the observed acceleration of the Universe. We consider the use of field redefinitions to cast k-essence in a more familiar form. While k-essence models cannot in general be rewritten in the form of quintessence models, we show that in certain dynamical regimes an equivalence can be made, which in particular can shed light on the tracking behaviour of k-essence. In several cases, k-essence cannot be observationally distinguished from quintessence using the homogeneous evolution, though there may be small effects on the perturbation spectrum. We make a detailed analysis of two k-essence models from the literature and comment on the nature of the fine tuning arising in the models.

We study the integrated Sachs-Wolfe (ISW) effect using a model-independent parameterization of the dark energy equation of state. Cosmic variance severely restricts the class of models distinguishable from Lambda-CDM. In particular if the present value of the equation of state satisfies w_Q^o <-0.8, the only distinguishable models are those with a rapidly varying equation of state, characterized by a late transition occurring at redshifts z<3. Due to the degeneracy with other cosmological parameters, models with a slowly varying behavior cannot be differentiated from each other or from a cosmological constant using even cosmic variance-limited CMB measurements. This may place a fundamental limit on our understanding of the origin of the presently observed acceleration.

The consensus of opinion in cosmology is that the Universe is currently undergoing a period of accelerated expansion. With current and proposed high precision experiments it offers the hope of being able to discriminate between the two competing models that are being suggested to explain the observations, namely a cosmological constant or a time dependent `Quintessence' model. The latter suffers from a plethora of scalar field potentials all leading to similar late time behaviour of the universe, hence to a lack of predictability. In this paper, we develop a model independent approach which simply involves parameterizing the dark energy equation of state in terms of known observables. This allows to analyse the impact dark energy has had on cosmology without the need to refer to particular scalar field models and opens up the possibility that future experiments will be able to constrain the dark energy equation of state in a model independent manner.

We study the instability of a scalar field at the end of hybrid inflation, using both analytical techniques and numerical simulations. We improve previous studies by taking the inflaton field fully into account, and show that the range of unstable modes depends sensitively on the velocity of the inflaton field, and thereby on the Hubble rate, at the end of inflation. If topological defects are formed, their number density is determined by the shortest unstable wavelength. Finally, we show that the oscillations of the inflaton field amplify the inhomogeneities in the energy density, leading to local symmetry restoration and faster thermalization. We believe this explains why tachyonic preheating is so effective in transferring energy away from the inflaton zero mode.

Quintessence has been introduced as an alternative to the cosmological constant scenario to account for the current acceleration of the universe. This new dark energy component allows values of the equation of state parameter $w_{Q}^0\geq-1$, and in principle measurements of cosmological distances to Type Ia supernovae can be used to distinguish between these two types of models. Assuming a flat universe, we use the supernovae data and measurements of the position of the acoustic peaks in the Cosmic Microwave Background (CMB) spectra to constrain a rather general class of Quintessence potentials, including inverse power law models and recently proposed Supergravity inspired potentials. In particular we use a likelihood analysis, marginalizing over the dark energy density $\Omega_{Q}$, the physical baryon density $\Omega_{b}h^2$ and the scalar spectral index $n$, to constrain the slopes of our Quintessence potential. Considering only the first Doppler peak the best fit in our range of models gives $w_{Q}^0\sim-0.8$. However, including the SnIa data and the three peaks, we find an upper limit on the present value of the equation of state parameter, $-1\leq w_{Q}^0\leq-0.93$ at $2\sigma$, a result that appears to rule out a class of recently proposed potentials.

We consider the construction of inverted hybrid inflation models in which the vacuum energy during inflation is at the TeV scale, and the inflaton couples to the Higgs field. Such models are of interest in the context of some recently proposed models of electroweak baryogenesis. We demonstrate how constraints on these models arise from quantum corrections, and how self-consistent examples may be constructed, albeit at the expense of fine-tuning. We discuss two possible ways in which the baryon asymmetry of the universe may be produced in these models. One of them is based on preheating and a consequent non-thermal electroweak symmetry restoration and the other on the formation of Higgs winding configurations by the Kibble mechanism at the end of inflation.

In a class of non-singular cosmologies derived from higher-order corrections to the low-energy bosonic string action, we derive evolution equations for the most general cosmological scalar, vector and tensor perturbations. In the large scale limit, the evolutions of both scalar and tensor perturbations are characterised by conserved quantities, the usual curvature perturbation in the uniform-field gauge and the tensor-type perturbed metric. The vector perturbation is not affected, being described by the conservation of the angular momentum of the fluid component in the absence of any additional dissipative process. For the scalar- and tensor-type perturbations, we show how, given a background evolution during kinetic driven inflation of the dilaton field, we can obtain the final power spectra generated from the vacuum quantum fluctuations of the metric and the dilaton field during the inflation.

We demonstrate how exponential potentials that could arise in the early Universe as a result of Kaluza-Klein type compactifications of string theory, can lead to cosmological solutions which correspond to the currently observed accelerating Universe. The idea is simple, relying solely on the known scaling properties associated with exponential potentials. In particular we show that the existence of stable attractor solutions implies that the results hold for a wide range of coupling constants and initial conditions.

We study the evolution of tensor metric fluctuations in a class of non-singular models based on the string effective action, by including in the perturbation equation the higher-derivative and loop corrections needed to regularise the background solutions. We discuss the effects of such higher-order corrections on the final graviton spectrum, and we compare the results of analytical and numerical computations.

We propose a scenario for inflation based upon the braneworld picture, in which high-energy corrections to the Friedmann equation permit inflation to take place with potentials ordinarily too steep to sustain it. Inflation ends when the braneworld corrections begin to lose their dominance. Reheating may naturally be brought about via gravitational particle production, rather than the usual inflaton decay mechanism; the reheat temperature may be low enough to satisfy the gravitino bound and the Universe becomes radiation dominated early enough for nucleosynthesis. We illustrate the idea by considering steep exponential potentials, and show they can give satisfactory density perturbations (both amplitude and slope) and reheat successfully. The scalar field may survive to the present epoch without violating observational bounds, and could be invoked in the quintessential inflation scenario of Peebles and Vilenkin.

Scalar field models of quintessence typically require that the expectation value of the field today is of order the Planck mass, if we want them to explain the observed acceleration of the Universe. This suggests that we should be considering models in the context of supergravity. We discuss a particular class of supergravity models and analyze their behavior under different choices of the Kahler metric.

We show by studying the Abelian Higgs model with numerical lattice simulations that non-thermal phase transitions arising out of preheating after inflation are possible in gauge-Higgs models under rather general circumstances. This may lead to the formation of gauged topological defects and, if the scale at which inflation ends is low enough, to electroweak baryogenesis after preheating.

We Study instanton solutions in general relativity with a scalar field. The metric ansatz we use is composed of a particular warp product of general Einstein metrics, such as those found in a number of cosmological settings, including string cosmology, supergravity compactifications and general Kaluza Klein reductions. Using the Hartle-Hawking prescription the instantons we obtain determine whether metrics involving extra compact dimensions of this type are favoured as initial conditions for the universe. Specifically, we find that these product metric instantons, viewed as constrained instantons, do have a local minima in the action. These minima are then compared with the higher dimensional version of the Hawking-Turok instantons, and we argue that the latter always have lower action than those associated with these product metrics.

Aspects of superstring cosmology are reviewed with an emphasis on the cosmological implications of duality symmetries in the theory. The string effective actions are summarized and toroidal compactification to four dimensions reviewed. Global symmetries that arise in the compactification are discussed and the duality relationships between the string effective actions are then highlighted. Higher-dimensional Kasner cosmologies are presented and interpreted in both string and Einstein frames, and then given in dimensionally reduced forms. String cosmologies containing both non-trivial Neveu-Schwarz/Neveu-Schwarz and Ramond-Ramond fields are derived by employing the global symmetries of the effective actions. Anisotropic and inhomogeneous cosmologies in four-dimensions are also developed. The review concludes with a detailed analysis of the pre-big bang inflationary scenario. The generation of primordial spectra of cosmological perturbations in such a scenario is discussed. Possible future directions offered in the Horava-Witten theory are outlined.

We perform a new and accurate study of the dependence on cosmological parameters of structure formation with local cosmic strings. The crucial new ingredients are the inclusion of the effects of gravitational backreaction on the evolution of the network, and the accurate evolution of the network through the radiation to matter transition. Our work re-iterates the fact that expanding Universe numerical simulations only probe a transient regime, and we incorporate our results into the unequal time correlators recently measured. We then compute the CMB and CDM fluctuations' power spectra for various values of the Hubble constant $H_0$ and baryon fraction $\Omega_b$. We find that, whereas the dependence on $\Omega_b$ is negligible, there is still a strong dependence on $H_0$.

The stability of binary orbits can significantly shape the gravity wave signal which future Earth-based interferometers hope to detect. The inner most stable circular orbit has been of interest as it marks the transition from the late inspiral to final plunge. We consider purely relativistic orbits beyond the circular assumption. Homoclinic orbits are of particular importance to the question of stability as they lie on the boundary between dynamical stability and instability. We identify these, estimate their rate of energy loss to gravity waves, and compute their gravitational waveforms.

We obtain a new class of exact cosmological solutions for multi-scalar fields with exponential potentials. We generalize the assisted inflation solutions previously obtained, and demonstrate how they are modified when there exist cross-couplings between the fields, such as occur in supergravity inspired cosmological models.

We present a phase-plane analysis of cosmologies containing a barotropic fluid with equation of state $p_\gamma = (\gamma-1) \rho_\gamma$, plus a scalar field $\phi$ with an exponential potential $V \propto \exp(-\lambda \kappa \phi)$ where $\kappa^2 = 8\pi G$. In addition to the well-known inflationary solutions for $\lambda^2 < 2$, there exist scaling solutions when $\lambda^2 > 3\gamma$ in which the scalar field energy density tracks that of the barotropic fluid (which for example might be radiation or dust). We show that the scaling solutions are the unique late-time attractors whenever they exist. The fluid-dominated solutions, where $V(\phi)/\rho_\gamma \to 0$ at late times, are always unstable (except for the cosmological constant case $\gamma = 0$). The relative energy density of the fluid and scalar field depends on the steepness of the exponential potential, which is constrained by nucleosynthesis to $\lambda^2 > 20$. We show that standard inflation models are unable to solve this `relic density' problem.

Accurate estimation of cosmological parameters from microwave background anisotropies requires high-accuracy understanding of the cosmological model. Normally, a power-law spectrum of density perturbations is assumed, in which case the spectral index $n$ can be measured to around $\pm 0.004$ using microwave anisotropy satellites such as MAP and Planck. However, inflationary models generically predict that the spectral index $n$ of the density perturbation spectrum will be scale-dependent. We carry out a detailed investigation of the measurability of this scale dependence by Planck, including the influence of polarization on the parameter estimation. We also estimate the increase in the uncertainty in all other parameters if the scale dependence has to be included. This increase applies even if the scale dependence is too small to be measured unless it is assumed absent, but is shown to be a small effect. We study the implications for inflation models, beginning with a brief examination of the generic slow-roll inflation situation, and then move to a detailed examination of a recently-devised hybrid inflation model for which the scale dependence of $n$ may be observable.

Pre--big bang models of inflation based on string cosmology produce a stochastic gravitational wave background whose spectrum grows with decreasing wavelength, and which may be detectable using interferometers such as LIGO. We point out that the gravitational wave spectrum is closely tied to the density perturbation spectrum, and that the condition for producing observable gravitational waves is very similar to that for producing an observable density of primordial black holes. Detection of both would provide strong support to the string cosmology scenario.

If primordial scalar and tensor perturbation spectra can be inferred from observations of the cosmic background radiation and large-scale structure, then one might hope to reconstruct a unique single-field inflaton potential capable of generating the observed spectra. In this paper we examine conditions under which such a potential can be reliably reconstructed. For it to be possible at all, the spectra must be well fit by a Taylor series expansion. A complete reconstruction requires a statistically-significant tensor mode to be measured in the microwave background. We find that the observational uncertainties dominate the theoretical error from use of the slow-roll approximation, and conclude that the reconstruction procedure will never insidiously lead to an irrelevant potential.

We review the relation between the inflationary potential and the spectra of density (scalar) perturbations and gravitational waves (tensor perturbations) produced, with particular emphasis on the possibility of reconstructing the inflaton potential from observations. The spectra provide a potentially powerful test of the inflationary hypothesis; they are not independent but instead are linked by consistency relations reflecting their origin from a single inflationary potential. To lowest-order in a perturbation expansion there is a single, now familiar, relation between the tensor spectral index and the relative amplitude of the spectra. We demonstrate that there is an infinite hierarchy of such consistency equations, though observational difficulties suggest only the first is ever likely to be useful. We also note that since observations are expected to yield much better information on the scalars than on the tensors, it is likely to be the next-order version of this consistency equation which will be appropriate, not the lowest-order one. If inflation passes the consistency test, one can then confidently use the remaining observational information to constrain the inflationary potential, and we survey the general perturbative scheme for carrying out this procedure. Explicit expressions valid to next-lowest order in the expansion are presented. We then briefly assess the prospects for future observations reaching the quality required, and consider a simulated data set that is motivated by this outlook.

In an Abelian gauge symmetry, spontaneously broken at a first order phase transition, we investigate the evolution of two and three bubbles of the broken symmetry phase. The full field equations are evolved and we concentrate in particular on gauge invariant quantities, such as the magnetic field and the integral around a closed loop of the phase gradient. An intriguing feature emerges, namely, the geodesic rule, commonly used in numerical simulations to determine the density of defects formed is shown not to hold in a number of circumstances. It appears to be a function of the initial separation of the bubbles, and the coupling strength of the gauge field. The reason for the breakdown is that in the collision region the radial mode {\it can} be excited and {\it often} oscillates about its symmetry restoring value rather than settling to its broken symmetry value. This can lead to extra windings being induced in these regions, hence extra defects (anti-defects) being formed.

Currently favoured extensions of the Standard Model typically contain `flaton fields' defined as fields with large vacuum expectation values (vevs) and almost flat potentials. If a flaton field is trapped at the origin in the early universe, one expects `thermal inflation' to take place before it rolls away to the true vacuum, because the finite-temperature correction to the potential will hold it at the origin until the temperature falls below $1\TeV$ or so. In the first part of the paper, that expectation is confirmed by an estimate of the finite temperature corrections and of the tunneling rate to the true vacuum, paying careful attention to the validity of the approximations that are used. The second part of the paper considers topological defects which may be produced at the end of an era of thermal inflation. If the flaton fields associated with the era are GUT higgs fields, then its end corresponds to the GUT phase transition. In that case monopoles (as well as GUT higgs particles) will have to be diluted by a second era of thermal inflation. Such an era will not affect the cosmology of GUT strings, for which the crucial parameter is the string mass per unit length. Because of the flat Higgs potential, the GUT symmetry breaking scale required for the strings to be a candidate for the origin of large scale structure and the cmb anisotropy is about three times bigger than usual, but given the uncertainties it is still compatible with the one required by the unification of the Standard Model gauge couplings. The cosmology of textures and of global monopoles is unaffected by the flatness of the potential.

Oscillons are localized, non-singular, time-dependent, spherically-symmetric solutions of nonlinear scalar field theories which, although unstable, are extremely long-lived. We show that they naturally appear during the collapse of subcritical bubbles in models with symmetric and asymmetric double-well potentials. By a combination of analytical and numerical work we explain several of their properties, including the conditions for their existence, their longevity, and their final demise. We discuss several contexts in which we expect oscillons to be relevant. In particular, their nucleation during cosmological phase transitions may have wide-ranging consequences.

Using the formalism developed in a previous paper we analyze the cosmological implications of our conclusions concerning the scaling behaviour of a network of cosmic strings, in particular the idea that once gravitational back-reaction becomes important the transient scaling regime so far explored by numerical simulations will be replaced by a new `full scaling' regime, with slightly different parameters. This has consequences for the normalization of the string tension on all scales. We show how the new scaling solutions affect the gravitational wave bound from nucleosynthesis and the milli-second pulsar, and also describe the consequences from large scale structure up to COBE scales.

A review is presented of recent work by the authors concerning the use of large scale structure and microwave background anisotropy data to determine the potential of the inflaton field. The importance of a detection of the stochastic gravitational wave background is emphasised, and some preliminary new results of tests of the method on simulated data sets with uncertainties are described. (Proceedings of ``Unified Symmetry in the Small and in the Large'', Coral Gables, 1994)

We use numerical simulations to calculate the cosmic microwave background anisotropy induced by the evolution of a global texture field, with special emphasis on individual textures. Both spherically symmetric and general configurations are analysed, and in the latter case we consider field configurations which exhibit unwinding events and also ones which do not. We compare the results given by evolving the field numerically under both the expanded core (XCORE) and non-linear sigma model (NLSM) approximations with the analytic predictions of the NLSM exact solution for a spherically symmetric self-similar (SSSS) unwinding. We find that the random unwinding configuration spots' typical peak height is 60--75\% and angular size typically only 10\% of those of the SSSS unwinding, and that random configurations without an unwinding event nonetheless may generate indistinguishable hot and cold spots. The influence of these results on analytic estimates of texture induced microwave anisotropies is examined, and comparison made with other numerical work.

We investigate chaotic inflation models with two scalar fields, such that one field (the inflaton) rolls while the other is trapped in a false vacuum state. The false vacuum becomes unstable when the inflaton field falls below some critical value, and a first or second order transition to the true vacuum ensues. Particular attention is paid to Linde's second-order `Hybrid Inflation'; with the false vacuum dominating, inflation differs from the usual true vacuum case both in its cosmology and in its relation to particle physics. The spectral index of the adiabatic density perturbation can be very close to 1, or it can be around ten percent higher. The energy scale at the end of inflation can be anywhere between $10^{16}$\,GeV and $10^{11}$\,GeV, though reheating is prompt so the reheat temperature can't be far below $10^{11}\,$GeV. Topological defects are almost inevitably produced at the end of inflation, and if the inflationary energy scale is near its upper limit they can have significant effects. Because false vacuum inflation occurs with the inflaton field far below the Planck scale, it is easier to implement in the context of supergravity than standard chaotic inflation. That the inflaton mass is small compared with the inflationary Hubble parameter is still a problem for generic supergravity theories, but remarkably this can be avoided in a natural way for a class of supergravity models which follow from orbifold compactification of superstrings. This opens up the prospect of a truly realistic, superstring

One method to reconstruct the scalar field potential of inflation is a perturbative approach, where the values of the potential and its derivatives are calculated as an expansion in departures from the slow-roll approximation. They can then be expressed in terms of observable quantities, such as the square of the ratio of the gravitational wave amplitude to the density perturbation amplitude, the deviation of the spectral index from the Harrison--Zel'dovich value, etc. Here, we calculate complete expressions for the second-order contributions to the coefficients of the expansion by including for the first time corrections to the standard expressions for the perturbation spectra. As well as offering an improved result, these corrections indicate the expected accuracy of the reconstruction. Typically the corrections are only a few percent.

We extend and develop our previous work on the evolution of a network of cosmic strings. The new treatment is based on an analysis of the probability distribution of the end-to-end distance of a randomly chosen segment of left-moving string of given length. The description involves three distinct length scales: $\xi$, related to the overall string density, $\bar\xi$, the persistence length along the string, and $\zeta$, describing the small-scale structure, which is an important feature of the numerical simulations that have been done of this problem. An evolution equation is derived describing how the distribution develops in time due to the combined effects of the universal expansion, of intercommuting and loop formation, and of gravitational radiation. With plausible assumptions about the unknown parameters in the model, we confirm the conclusions of our previous study, that if gravitational radiation and small-scale structure effects are neglected, the two dominant length scales both scale in proportion to the horizon size. When the extra effects are included, we find that while $\xi$ and $\bar\xi$ grow, $\zeta$ initially does not. Eventually, however, it does appear to scale, at a much lower level, due to the effects of gravitational back-reaction.

We show how observations of the density perturbation (scalar) spectrum and the gravitational wave (tensor) spectrum allow a reconstruction of the potential responsible for cosmological inflation. A complete functional reconstruction or a perturbative approximation about a single scale are possible; the suitability of each approach depends on the data available. Consistency equations between the scalar and tensor spectra are derived, which provide a powerful signal of inflation.

Generalizing the original work by Hodges and Blumenthal, we outline a formalism which allows one, in principle, to reconstruct the potential of the inflaton field from knowledge of the tensor gravitational wave spectrum or the scalar density fluctuation spectrum, with special emphasis on the importance of the tensor spectrum. We provide some illustrative examples of such reconstruction. We then discuss in some detail the question of whether one can use real observations to carry out this procedure. We conclude that in practice, a full reconstruction of the functional form of the potential will not be possible within the foreseeable future. However, with a knowledge of the dark matter components, it should soon be possible to combine intermediate-scale data with measurements of large-scale cosmic microwave background anisotropies to yield useful information regarding the potential.