We derive quantum Boltzmann equations for preheating by means of the density matrix formalism, which account for both the non-adiabatic particle production and the leading collisional processes between the produced particles. In so doing, we illustrate the pivotal role played by pair correlations in mediating the particle production. In addition, by numerically solving the relevant system of Boltzmann equations, we demonstrate that collisional processes lead to a suppression of the growth of the number density even at the very early stages of preheating.

Light degrees of freedom that modify gravity on cosmological scales must be "screened" on solar system scales in order to be compatible with data. The Vainshtein mechanism achieves this through a breakdown of classical perturbation theory, as large interactions involving new degrees of freedom become important below the so-called Vainshtein radius. We begin to develop an extension of the Parameterized Post-Newtonian (PPN) formalism that is able to handle Vainshteinian corrections. We argue that theories with a unique Vainshtein scale must be expanded using two small parameters. In this Parameterized Post-Newtonian-Vainshteinian (PPNV) expansion, the primary expansion parameter that controls the PPN order is, as usual, the velocity $v$. The secondary expansion parameter, $\alpha$, controls the strength of the Vainshteinian correction and is a theory-specific combination of the Schwarzschild radius and the Vainshtein radius of the source that is independent of its mass. We present the general framework and apply it to Cubic Galileon theory both inside and outside the Vainshtein radius. The PPNV framework can be used to determine the compatibility of such theories with solar system and other strong-field data.

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.

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.

Current data are in good agreement with the predictions of single field inflation. However, the hemispherical asymmetry seen in the cosmic microwave background data, may hint at a potential problem. Generalizing to multi-field models may provide one possible explanation. A useful way of modeling perturbations in multi-field inflation is to investigate the projection of the perturbation along and perpendicular to the background fields' trajectory. These correspond to the adiabatic and isocurvature perturbations. However, it is important to note that in general there are no corresponding adiabatic and isocurvature fields. The purpose of this article is to highlight the distinction between a field redefinition and a perturbation projection. We provide a detailed derivation of the evolution of the isocurvature perturbation to show that no assumption of an adiabatic or isocurvature field is needed. We also show how this evolution equation is consistent with the field covariant evolution equations for the isocurvature perturbation in the flat field space limit.

We reconsider cosmological perturbation theory for multi-component scalars, enforcing covariance in field-space, and ensuring that phyical observations are independent of field re-definitions. We use the formalism to clarify some issues in the literature, and use pseudo-supersymmetry to derive exact expressions for terms of interest.

We continue to introduce bi-galileon theory, the generalisation of the single galileon model introduced by Nicolis et al. The theory contains two coupled scalar fields and is described by a Lagrangian that is invariant under Galilean shifts in those fields. This paper is the second of two, and focuses on the phenomenology of the theory. We are particularly interesting in models that admit solutions that are asymptotically self accelerating or asymptotically self tuning. In contrast to the single galileon theories, we find examples of self accelerating models that are simultaneously free from ghosts, tachyons and tadpoles, able to pass solar system constraints through Vainshtein screening, and do not suffer from problems with superluminality, Cerenkov emission or strong coupling. We also find self tuning models and discuss how Weinberg's no go theorem is evaded by breaking Poincar\'e invariance in the scalar sector. Whereas the galileon description is valid all the way down to solar system scales for the self-accelerating models, unfortunately the same cannot be said for self tuning models owing to the scalars backreacting strongly on to the geometry.

We show how to generalize the classical duals found by Gabadadze {\it et al} to a very large class of self-interacting theories. This enables one to adopt a perturbative description beyond the scale at which classical perturbation theory breaks down in the original theory. This is particularly relevant if we want to test modified gravity scenarios that exhibit Vainshtein screening on solar system scales. We recognise the duals as being related to the Legendre transform of the original Lagrangian, and present a practical method for finding the dual in general; our methods can also be applied to self-interacting theories with a hierarchy of strong coupling scales, and with multiple fields. We find the classical dual of the full quintic galileon theory as an example.

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.

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.

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.

The evolution of multiple scalar fields in cosmology has been much studied, particularly when the potential is formed from a series of exponentials. For a certain subclass of such systems it is possible to get `assisted` behaviour, where the presence of multiple terms in the potential effectively makes it shallower than the individual terms indicate. It is also known that when compactifying on coset spaces one can achieve a consistent truncation to an effective theory which contains many exponential terms, however, if there are too many exponentials then exact scaling solutions do not exist. In this paper we study the potentials arising from such compactifications of eleven dimensional supergravity and analyse the regions of parameter space which could lead to scaling behaviour.

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

Within the context of finding the initial conditions of the universe we consider gravitational instantons falling into the Bianchi IX classification. That is, a Euclidean four-manifold with a metric that satisfies Einstein's equations with an induced metric on S^3 submanifolds that is homogeneous but anisotropic. As well as finding regular solutions to the field equations with a tunnelling scalar field, we also look at the case of singular instantons with a view to applying the results to generic potentials. The study is in agreement with the prejudice that instantons with higher symmetry have a lower Euclidean action, even when we consider the singular class of solutions. It is also found that the Euclidean action can diverge for simple potentials, showing that the Hawking Turok instanton had finite action owing to its symmetry.

We report some new exact instantons in general relativity. These solutions are K\"ahler and fall into the symmetry classes of Bianchi types VI0 and VII0, with matter content of a stiff fluid. The qualitative behaviour of the solutions is presented, and we compare it to the known results of the corresponding self-dual Bianchi solutions. We also give axisymmetric Bianchi VII0 solutions with an electromagnetic field.

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.