We review current theoretical cosmology, including fundamental and mathematical cosmology and physical cosmology (as well as cosmology in the quantum realm), with an emphasis on open questions.

It is necessary to make assumptions in order to derive models to be used for cosmological predictions and comparison with observational data. In particular, in standard cosmology the spatial curvature is assumed to be constant and zero (or at least very small). But there is, as yet, no fully independent constraint with an appropriate accuracy that gaurentees a value for the magnitude of the effective normalized spatial curvature $\Omega_{k}$ of less than approximately $0.01$. Moreover, a small non-zero measurement of $\Omega_{k}$ at such a level perhaps indicates that the assumptions in the standard model are not satisfied. It has also been increasingly emphasised that spatial curvature is, in general, evolving in relativistic cosmological models. We review the current situation, and conclude that the possibility of such a non-zero value of $\Omega_k$ should be taken seriously.

In the era of precision cosmology, even percentage level effects are significant on cosmological observables. The recent tension between the local and global values of $H_0$ is much more significant than this, and any possible solution might rely on us going beyond the standard $\Lambda$CDM cosmological model. For much smaller, yet potentially significant effects, spatial curvature from averaging and cosmological backreaction on observational predictions could play a role. This is especially true with the higher precision of new observational data and improved statistical techniques. In this paper, we discuss the observational viability of a class of physically motivated cosmologies which can be parametrized by a phenomenological two-scale backreaction model with decoupled spatial curvature parameters and two Hubble scales. Using the latest JLA Supernovae data together with some of the latest BAO data, we perform a Bayesian model selection analysis and find that the phenomenological models are not favoured over the standard $\Lambda$CDM cosmological model. Although there is still a preference for non-zero and unequal dynamic and geometric spatial curvatures, there is little evidence for differing Hubble scales within these phenomenological template models.

In this paper we explore the idea that black holes can persist in a universe that collapses to a big crunch and then bounces into a new phase of expansion. We use a scalar field to model the matter content of such a universe {near the time} of the bounce, and look for solutions that represent a network of black holes within a dynamical cosmology. We find exact solutions to Einstein's constraint equations that provide the geometry of space at the minimum of expansion and that can be used as initial data for the evolution of hyperspherical cosmologies. These solutions illustrate that there exist models in which multiple distinct black holes can persist through a bounce, and allow for concrete computations of quantities such as the black hole filling factor. We then consider solutions in flat cosmologies, as well as in higher-dimensional spaces (with up to nine spatial dimensions). We derive conditions for the black holes to remain distinct (i.e. avoid merging) and hence persist into the new expansion phase. Some potentially interesting consequences of these models are also discussed.

General relativity is a set of physical and geometric principles, which lead to a set of (Einstein) field equations that determine the gravitational field, and to the geodesic equations that describe light propagation and the motion of particles on the background. But open questions remain, including: What is the scale on which matter and geometry are dynamically coupled in the Einstein equations? Are the field equations valid on small and large scales? What is the largest scale on which matter can be coarse grained while following a geodesic of a solution to Einstein's equations? We address these questions. If the field equations are causal evolution equations, whose average on cosmological scales is not an exact solution of the Einstein equations, then some simplifying physical principle is required to explain the statistical homogeneity of the late epoch Universe. Such a principle may have its origin in the dynamical coupling between matter and geometry at the quantum level in the early Universe. This possibility is hinted at by diverse approaches to quantum gravity which find a dynamical reduction to two effective dimensions at high energies on one hand, and by cosmological observations which are beginning to strongly restrict the class of viable inflationary phenomenologies on the other. We suggest that the foundational principles of general relativity will play a central role in reformulating the theory of spacetime structure to meet the challenges of cosmology in the 21st century.

An important problem in precision cosmology is the determination of the effects of averaging and backreaction on observational predictions, particularly in view of the wealth of new observational data and improved statistical techniques. In this paper, we discuss the observational viability of a class of averaged cosmologies which consist of a simple parametrized phenomenological two-scale backreaction model with decoupled spatial curvature parameters. We perform a Bayesian model selection analysis and find that this class of averaged phenomenological cosmological models is favored with respect to the standard $\Lambda$CDM cosmological scenario when a joint analysis of current SNe Ia and BAO data is performed. In particular, the analysis provides observational evidence for non-trivial spatial curvature.

We summarise some of the main observational challenges for the standard Friedmann-Lemaitre-Robertson-Walker cosmological model and describe how results recently presented in the parallel session `Large--scale Structure and Statistics' (DE3) at the `Fourteenth Marcel Grossman Meeting on General Relativity' are related to these challenges.

We investigate spherically symmetric cosmological models in Einstein-aether theory with a tilted (non-comoving) perfect fluid source. We use a 1+3 frame formalism and adopt the comoving aether gauge to derive the evolution equations, which form a well-posed system of first order partial differential equations in two variables. We then introduce normalized variables. The formalism is particularly well-suited for numerical computations and the study of the qualitative properties of the models, which are also solutions of Horava gravity. We study the local stability of the equilibrium points of the resulting dynamical system corresponding to physically realistic inhomogeneous cosmological models and astrophysical objects with values for the parameters which are consistent with current constraints. In particular, we consider dust models in ($\beta-$) normalized variables and derive a reduced (closed) evolution system and we obtain the general evolution equations for the spatially homogeneous Kantowski-Sachs models using appropriate bounded normalized variables. We then analyse these models, with special emphasis on the future asymptotic behaviour for different values of the parameters. Finally, we investigate static models for a mixture of a (necessarily non-tilted) perfect fluid with a barotropic equations of state and a scalar field.

No. In a number of papers Green and Wald argue that the standard FLRW model approximates our Universe extremely well on all scales, except close to strong field astrophysical objects. In particular, they argue that the effect of inhomogeneities on average properties of the Universe (backreaction) is irrelevant. We show that this latter claim is not valid. Specifically, we demonstrate, referring to their recent review paper, that (i) their two-dimensional example used to illustrate the fitting problem differs from the actual problem in important respects, and it assumes what is to be proven; (ii) the proof of the trace-free property of backreaction is unphysical and the theorem about it fails to be a mathematically general statement; (iii) the scheme that underlies the trace-free theorem does not involve averaging and therefore does not capture crucial non-local effects; (iv) their arguments are to a large extent coordinate-dependent, and (v) many of their criticisms of backreaction frameworks do not apply to the published definitions of these frameworks. It is therefore incorrect to infer that Green and Wald have proven a general result that addresses the essential physical questions of backreaction in cosmology.

Averaging in general relativity is a complicated operation, due to the general covariance of the theory and the non-linearity of Einstein's equations. The latter of these ensures that smoothing spacetime over cosmological scales does not yield the same result as solving Einstein's equations with a smooth matter distribution, and that the smooth models we fit to observations need not be simply related to the actual geometry of spacetime. One specific consequence of this is a decoupling of the geometrical spatial curvature term in the metric from the dynamical spatial curvature in the Friedmann equation. Here we investigate the consequences of this decoupling by fitting to a combination of HST, CMB, SNIa and BAO data sets. We find that only the geometrical spatial curvature is tightly constrained, and that our ability to constrain dark energy dynamics will be severely impaired until we gain a thorough understanding of the averaging problem in cosmology.

We use cosmological perturbation theory to study the backreaction effects of a self-consistent and well-defined cosmological averaging on the dynamics and the evolution of the Universe. Working with a perturbed Friedman-Lemaitre-Robertson-Walker Einstein-de Sitter cosmological solution in a comoving volume-preserving gauge, we compute the expressions for the expansion scalar and deceleration parameter to second order, which we use to characterize the backreaction. We find that the fractional shift in the Hubble parameter with respect to the input background cosmological model is Delta~10^{-5}, which leads to an effective energy density of the order of a few times 10^{-5}. In addition, we find that an appropriate measure of the fractional shift in the deceleration parameter Q is very large.

The cosmological backreaction from perturbations is clearly gauge-dependent, and obviously depends on the choice of averaged Hubble rate. We consider two common choices of Hubble rate and advocate the use of comoving volume-preserving gauges. We highlight two examples valid to an appropriate order in perturbation theory, uniform curvature gauge, which is as close to volume-preserving as possible, and a spatially-traceless uniform cold dark matter gauge which preserves the volume to linear order. We demonstrate the strong gauge- and frame-dependences in averaging. In traceless uniform CDM gauge the backreaction exhibits a strong ultra-violet divergence and can be tuned to an arbitrary magnitude with an appropriate choice of smoothing scale. In uniform curvature gauge we find that for a choice of Hubble rate locked to the spatial surface the backreaction vanishes identically, while for a Hubble rate defined from a fluid's expansion scalar the effective energy density at the current epoch in an Einstein-de Sitter universe is Omega_eff~5e-4, slightly bigger than but in broad agreement with previous results in conformal Newtonian gauge.

We use a dynamical systems analysis to investigate the future behaviour of Einstein-Aether cosmological models with a scalar field coupling to the expansion of the aether and a non-interacting perfect fluid. The stability of the equilibrium solutions are analysed and the results are compared with the standard inflationary cosmological solutions and previously studied cosmological Einstein-Aether models.

We discuss the construction of cosmological models within the framework of Macroscopic Gravity (MG), which is a theory that models the effects of averaging the geometry of space-time on large scales. We find new exact spatially homogeneous and isotropic FLRW solutions to the MG field equations, and investigate large-scale perturbations around them. We find that any inhomogeneous perturbations to the averaged geometry are severely restricted, but that possible anisotropies in the correlation tensor can have dramatic consequences for the measurement of distances. These calculations are a first step within the MG approach toward developing averaged cosmological models to a point where they can be used to interpret real cosmological data, and hence to provide a working alternative to the "concordance" LCDM model.

We discuss whether black holes could persist in a universe which recollapses and then bounces into a new expansion phase. Whether the bounce is of classical or quantum gravitational origin, such cosmological models are of great current interest. In particular, we investigate the mass range in which black holes might survive a bounce and ways of differentiating observationally between black holes formed just after and just before the last bounce. We also discuss the consequences of the universe going through a sequence of dimensional changes as it passes through a bounce.

The averaging problem in cosmology is of considerable importance for the correct interpretation of cosmological data. We review cosmological observations and discuss some of the issues regarding averaging. We present a precise definition of a cosmological model and a rigorous mathematical definition of averaging, based entirely in terms of scalar invariants.

The averaging problem in cosmology is of considerable importance for the correct interpretation of cosmological data. A rigorous mathematical definition of averaging in a cosmological model is necessary. In general, a spacetime is completely characterized by its scalar curvature invariants, and this suggests a particular spacetime averaging scheme based entirely on scalars. We clearly identify the problems of averaging in a cosmological model. We then present a precise definition of a cosmological model, and based upon this definition, we propose an averaging scheme in terms of scalar curvature invariants. This scheme is illustrated in a simple static spherically symmetric perfect fluid cosmological spacetime, where the averaging scales are clearly identified.

The universe is not isotropic or spatially homogeneous on local scales. The averaging of local inhomogeneities in general relativity can lead to significant dynamical effects on the evolution of the universe and on the interpretation of cosmological data. In particular, all deductions about cosmology are based on light paths; averaging can have an important effect on photon propagation and hence cosmological observations. It would be desirable to describe the physical effects of averaging in terms of observational quantities and focussing on the behaviour of light. Data (e.g., matter terms, such as the density or galaxy number counts, which are already expressed as averaged quantities) is given on the null cone. Therefore, it is observationally meaningful to consider light-cone averages of quantities. In principle, we wish to describe the cosmological equations on the null cone, and hence we need to construct the averaged geometry on the null cone. However, we argue that it is still necessary to average the full Einstein field equations to obtain suitably averaged equations on the null cone. Since it is not the geometry per se that appears in the observational relations, we discuss whether it is possible to covariantly `average' just a subset of the evolution equations on the null cone, focussing on relevant observational quantities. We present an averaged version of the scalar null Raychaudhuri equation, which may be a useful first step in this regard.

The Universe is not isotropic or spatially homogeneous on local scales. The averaging of local inhomogeneities in general relativity can lead to significant dynamical effects on the evolution of the Universe, and even if the effects are at the 1% level they must be taken into account in a proper interpretation of cosmological observations. We discuss the effects that averaging (and inhomogeneities in general) can have on the dynamical evolution of the Universe and the interpretation of cosmological data. All deductions about cosmology are based on the paths of photons. We discuss some qualitative aspects of the motion of photons in an averaged geometry, particularly within the context of the luminosity distance-redshift relation in the simple case of spherical symmetry.

We study the late-time behaviour of tilted perfect fluid Bianchi type III models using a dynamical systems approach. We consider models with dust, and perfect fluids stiffer than dust, and eludicate the late-time behaviour by studying the centre manifold which dominates the behaviour of the model at late times. In the dust case, this centre manifold is 3-dimensional and can be considered as a double bifurcation as the 2 parameters ($h$ and $\gamma$) of the type VI$_h$ model are varied. We calculate the decay rates and show that for dust or stiffer the models approach a vacuum spacetime, however, it does so rather slowly: $\rho/H^2\sim 1/\ln t$.

We study tilted perfect fluid cosmological models with a constant equation of state parameter in spatially homogeneous models of Bianchi type VI$_{-1/9}$ using dynamical systems methods and numerical simulations. We study models with and without vorticity, with an emphasis on their future asymptotic evolution. We show that for models with vorticity there exists, in a small region of parameter space, a closed curve acting as the attractor.

We study tilted perfect fluid cosmological models with a constant equation of state parameter in spatially homogeneous models of Bianchi type VI_h using dynamical systems methods and numerical experimentation, with an emphasis on their future asymptotic evolution. We determine all of the equilibrium points of the type VI_h state space (which correspond to exact self-similar solutions of the Einstein equations, some of which are new), and their stability is investigated. We find that there are vacuum plane-wave solutions that act as future attractors. In the parameter space, a `loophole' is shown to exist in which there are no stable equilibrium points. We then show that a Hopf-bifurcation can occur resulting in a stable closed orbit (which we refer to as the Mussel attractor) corresponding to points both inside the loophole and points just outside the loophole; in the former case the closed curves act as late-time attractors while in the latter case these attracting curves will co-exist with attracting equilibrium points. In the special Bianchi type III case, centre manifold theory is required to determine the future attractors. Comprehensive numerical experiments are carried out to complement and confirm the analytical results presented. We note that the Bianchi type VI_h case is of particular interest in that it contains many different subcases which exhibit many of the different possible future asymptotic behaviours of Bianchi cosmological models.

In the macroscopic gravity approach to the averaging problem in cosmology, the Einstein field equations on cosmological scales are modified by appropriate gravitational correlation terms. We study the averaging problem within the class of spherically symmetric cosmological models. That is, we shall take the microscopic equations and effect the averaging procedure to determine the precise form of the correlation tensor in this case. In particular, by working in volume preserving coordinates, we calculate the form of the correlation tensor under some reasonable assumptions on the form for the inhomogeneous gravitational field and matter distribution. We find that the correlation tensor in a FLRW background must be of the form of a spatial curvature. Inhomogeneities and spatial averaging, through this spatial curvature correction term, can have a very significant dynamical effect on the dynamics of the Universe and cosmological observations; in particular, we discuss whether spatial averaging might lead to a more conservative explanation of the observed acceleration of the Universe (without the introduction of exotic dark matter fields). We also find that the correlation tensor for a non-FLRW background can be interpreted as the sum of a spatial curvature and an anisotropic fluid. This may lead to interesting effects of averaging on astrophysical scales. We also discuss the results of averaging an inhomogeneous Lemaitre-Tolman-Bondi solution as well as calculations of linear perturbations (that is, the back-reaction) in an FLRW background, which support the main conclusions of the analysis.

We investigate the effect of dark energy on the limits on the shear anisotropy in spatially homogeneous Bianchi cosmological models obtained from measurements of the temperature anisotropies in the cosmic microwave background. We shall primarily assume that the dark energy is modelled by a cosmological constant. In general, we find that there are tighter bounds on the shear than in models with no cosmological constant, although the limits are (Bianchi) model dependent. In addition, there are special spatially homogeneous cosmological models whose rate of expansion is highly anisotropic, but whose cosmic microwave background temperature is measured to be exactly isotropic at one instant of time.

We study the properties of a future singularity encountered by a perfect fluid observer in tilting spatially homogeneous Bianchi cosmologies. We derive the boost formulae for the Weyl tensor to establish that, for two observers that are asymptotically null with respect to each other, their respective Weyl parameters generally both tend to zero, constant, or infinity together. We examine three classes of typical examples and one exceptional class. Given the behaviour of the Weyl parameter, we can predict that the singularity encountered is a Weyl singularity or a kinematic singularity. The analysis suggests that the kinematic variables are also useful in indicating a singularity in these models.

We show that in tilting perfect fluid cosmological models with an ultra-radiative equation of state, generically the tilt becomes extreme at late times and, as the tilt instability sets in, observers moving with the tilting fluid will experience singular behaviour in which infinite expansion is reached within a finite proper time, similar to that of phantom cosmology (but without the need for exotic forms of matter).

We discuss the averaging problem in general relativity, using the form of the macroscopic gravity equations in the case of spherical symmetry in volume preserving coordinates. In particular, we calculate the form of the correlation tensor under some reasonable assumptions on the form for the inhomogeneous gravitational field and matter distribution. On cosmological scales, the correlation tensor in a Friedmann-Lema\^{\i}tre-Robertson-Walker (FLRW) background is found to be of the form of a spatial curvature. On astrophysical scales the correlation tensor can be interpreted as the sum of a spatial curvature and an anisotropic fluid. We briefly discuss the physical implications of these results.

We use expansion-normalised variables to investigate the Bianchi type VII$_0$ model with a tilted $\gamma$-law perfect fluid. We emphasize the late-time asymptotic dynamical behaviour of the models and determine their asymptotic states. Unlike the other Bianchi models of solvable type, the type VII$_0$ state space is unbounded. Consequently we show that, for a general non-inflationary perfect fluid, one of the curvature variables diverges at late times, which implies that the type VII$_0$ model is not asymptotically self-similar to the future. Regarding the tilt velocity, we show that for fluids with $\gamma<4/3$ (which includes the important case of dust, $\gamma=1$) the tilt velocity tends to zero at late times, while for a radiation fluid, $\gamma=4/3$, the fluid is tilted and its vorticity is dynamically significant at late times. For fluids stiffer than radiation ($\gamma>4/3$), the future asymptotic state is an extremely tilted spacetime with vorticity.

Self-similar models are important in general relativity and other fundamental theories. In this paper we shall discuss the ``similarity hypothesis'', which asserts that under a variety of physical circumstances solutions of these theories will naturally evolve to a self-similar form. We will find there is good evidence for this in the context of both spatially homogenous and inhomogeneous cosmological models, although in some cases the self-similar model is only an intermediate attractor. There are also a wide variety of situations, including critical pheneomena, in which spherically symmetric models tend towards self-similarity. However, this does not happen in all cases and it is it is important to understand the prerequisites for the conjecture.

In the macroscopic gravity approach to the averaging problem in cosmology, the Einstein field equations on cosmological scales are modified by appropriate gravitational correlation terms. We present exact cosmological solutions to the equations of macroscopic gravity for a spatially homogeneous and isotropic macroscopic space-time and find that the correlation tensor is of the form of a spatial curvature term. We briefly discuss the physical consequences of these results.

Recently it was shown that the exact cosmological solutions known as the vacuum plane-wave solutions are late-time attractors for an open set of the spatially homogeneous Bianchi universes containing a non-inflationary $\gamma$-law perfect fluid. In this paper we study inhomogeneous perturbations of these plane-wave spacetimes. By using expansion-normalised scale-invariant variables we show that these solutions are unstable to generic inhomogeneous perturbations. The crucial observation for establishing this result is a divergence of the expansion-normalised frame variables which ultimately leads to unstable modes.

Brane-world cosmology is motivated by recent developments in string/M-theory and offers a new perspective on the hierarchy problem. In the brane-world scenario, our Universe is a four-dimensional subspace or {\em brane} embedded in a higher-dimensional {\em bulk} spacetime. Ordinary matter fields are confined to the brane while the gravitational field can also propagate in the bulk, leading to modifications of Einstein's theory of general relativity at high energies. In particular, the Randall-Sundrum-type models are self-consistent and simple and allow for an investigation of the essential non-linear gravitational dynamics. The governing field equations induced on the brane differ from the general relativistic equations in that there are nonlocal effects from the free gravitational field in the bulk, transmitted via the projection of the bulk Weyl tensor, and the local quadratic energy-momentum corrections, which are significant in the high-energy regime close to the initial singularity. In this review we discuss the asymptotic dynamical evolution of spatially homogeneous brane-world cosmological models containing both a perfect fluid and a scalar field close to the initial singularity. Using dynamical systems techniques it is found that, for models with a physically relevant equation of state, an isotropic singularity is a past-attractor in all orthogonal spatially homogeneous models (including Bianchi type IX models). In addition, we describe the dynamics in a class of inhomogeneous brane-world models, and show that these models also have an isotropic initial singularity. These results provide support for the conjecture that typically the initial cosmological singularity is isotropic in brane-world cosmology.

It has recently been suggested that an isotropic singularity may be a generic feature of brane cosmologies, even in the inhomogeneous case. Using the covariant and gauge-invariant approach we present a detailed analysis of linear perturbations of the isotropic model ${\cal F}_b$ which is a past attractor in the phase space of homogeneous Bianchi models on the brane. We find that for matter with an equation of state parameter $\gamma > 1$, the dimensionless variables representing generic anisotropic and inhomogeneous perturbations decay as $t\to 0$, showing that the model ${\cal F}_b$ is asymptotically stable in the past. We conclude that brane universes are born with isotropy naturally built-in, contrary to standard cosmology. The observed large-scale homogeneity and isotropy of the universe can therefore be explained as a consequence of the initial conditions if the brane-world paradigm represents a description of the very early universe.

The Observed Perlmutter-Riess Acceleration (OPRA) implies that the expansion of the Universe is currently increasing and is motivation for the so-called phantom energy models. We consider the dynamics of phantom scalar field models. An important physical time constraint, which can be used to rule out many cosmological models, is obtained from the condition that all forms of energy density, including the field causing OPRA (e.g., the phantom field), must be non-negligible for an extended period, which is conservatively estimated to be of the order of a few $Gyr$. We find that this physical time constraint cannot be satisfied in conventional phantom cosmological models.

In brane-world cosmology gravitational waves can propagate in the higher dimensions (i.e., in the `bulk'). In some appropriate regimes the bulk gravitational waves may be approximated by plane waves. We systematically study five-dimensional gravitational waves that are algebraically special and of type N. In the most physically relevant case the projected non-local stress tensor on the brane is formally equivalent to the energy-momentum tensor of a null fluid. Some exact solutions are studied to illustrate the features of these branes; in particular, we show explicity that any plane wave brane can be embedded into a 5-dimensional Siklos spacetime. More importantly, it is possible that in some appropriate regime the bulk can be approximated by gravitational plane waves and thus may act as initial conditions for the gravitational field in the bulk (thereby enabling the field equations to be integrated on the brane).

Using dynamical systems theory and a detailed numerical analysis, the late-time behaviour of tilting perfect fluid Bianchi models of types IV and VII$_h$ are investigated. In particular, vacuum plane-wave spacetimes are studied and the important result that the only future attracting equilibrium points for non-inflationary fluids are the plane-wave solutions in Bianchi type VII$_h$ models is discussed. A tiny region of parameter space (the loophole) in the Bianchi type IV model is shown to contain a closed orbit which is found to act as an attractor (the Mussel attractor). From an extensive numerical analysis it is found that at late times the normalised energy-density tends to zero and the normalised variables 'freeze' into their asymptotic values. A detailed numerical analysis of the type VII$_h$ models then shows that there is an open set of parameter space in which solution curves approach a compact surface that is topologically a torus.

We use a dynamical systems approach to analyse the tilting spatially homogeneous Bianchi models of solvable type (e.g., types VI$_h$ and VII$_h$) with a perfect fluid and a linear barotropic $\gamma$-law equation of state. In particular, we study the late-time behaviour of tilted Bianchi models, with an emphasis on the existence of equilibrium points and their stability properties. We briefly discuss the tilting Bianchi type V models and the late-time asymptotic behaviour of irrotational Bianchi VII$_0$ models. We prove the important result that for non-inflationary Bianchi type VII$_h$ models vacuum plane-wave solutions are the only future attracting equilibrium points in the Bianchi type VII$_h$ invariant set. We then investigate the dynamics close to the plane-wave solutions in more detail, and discover some new features that arise in the dynamical behaviour of Bianchi cosmologies with the inclusion of tilt. We point out that in a tiny open set of parameter space in the type IV model (the loophole) there exists closed curves which act as attracting limit cycles. More interestingly, in the Bianchi type VII$_h$ models there is a bifurcation in which a set of equilibrium points turn into closed orbits. There is a region in which both sets of closed curves coexist, and it appears that for the type VII$_h$ models in this region the solution curves approach a compact surface which is topologically a torus.

We present the complete set of propagation and constraint equations for the kinematic and non-local first order quantities which describe general linear inhomogeneous and anisotropic perturbations of a flat FRW braneworld with vanishing cosmological constant and decompose them in the standard way into their scalar, vector and tensor contributions. A detailed analysis of the perturbation dynamics is performed using dimensionless variables that are specially tailored for the different regimes of interest; namely, the low energy GR regime, the high energy regime and the dark energy regime. Tables are presented for the evolution of all the physical quantities, making it easy to do a detailed comparison of the past asymptotic behaviour of the perturbations of these models. We find results that exactly match those obtained in the analysis of the spatially inhomogeneous $G_{2}$ braneworld cosmologies presented recently; i.e., that isotropization towards the ${\cal F}_b$ model occurs for $\gamma > 4/3$.

We propose a new, self-consistent and dynamical scenario which gives rise to well-defined initial conditions for five-dimensional brane-world cosmologies with radion stabilization. At high energies, the five-dimensional effective theory is assumed to have a scale invariance so that it admits an expanding scaling solution as a future attractor. The system automatically approaches the scaling solution and, hence, the initial condition for the subsequent low-energy brane cosmology is set by the scaling solution. At low energies, the scale invariance is broken and a radion stabilization mechanism drives the dynamics of the brane-world system. We present an exact, analytic scaling solution for a class of scale-invariant effective theories of five-dimensional brane-world models which includes the five-dimensional reduction of the Horava-Witten theory, and provide convincing evidence that the scaling solution is a future attractor.

We use the cosmic microwave background temperature anisotropy to place limits on large-scale magnetic fields in an inhomogeneous (perturbed Friedmann) universe. If no assumptions are made about the spacetime geometry, only a weak limit can be deduced directly from the CMB. In the special case where spatial inhomogeneity is neglected to first order, the upper limit is much stronger, i.e. a few nano-G

We study the behavior of spatially homogeneous brane-worlds close to the initial singularity in the presence of both local and nonlocal stresses. It is found that the singularity in these brane-worlds can be locally either isotropic or anisotropic. We then investigate the Weyl curvature conjecture, according to which some measure of the Weyl curvature is related to a gravitational entropy. In particular, we study the Weyl curvature conjecture on the brane with respect to the dimensionless ratio of the Weyl invariant and the Ricci square and the measure proposed by Gr{\o}n and Hervik. We also argue that the Weyl curvature conjecture should be formulated on brane (i.e., in the four-dimensional context).

We examine, from both a qualitative and a numerical point of view, the evolution of Kantowski--Sachs cosmological models whose source is a mixture of a gas of weakly interacting massive particles (WIMP's) and a radiative gas made up of a ``tightly coupled'' mixture of electrons, baryons and photons from the end of nucleosynthesis up to the duration of radiative interactions. In this cosmic era annihilation processes are negligible, while the WIMP's only interact gravitationally with the radiative gas and the latter behaves as a single dissipative fluid that can be studied within a hydrodynamical framework. Applying the full transport equations of Extended Irreversible Thermodynamics, coupled with the field and balance equations, we obtain a set of governing equations that becomes an autonomous system of ordinary differential equations once the shear viscosity relaxation time is specified; assuming the latter to be proportional to the Hubble time, the qualitative analysis indicates that models begin in the radiation dominated epoch close to an isotropic equilibrium point. We show how the relaxation time governs the relaxation of the models towards an equilibrium photon entropy, leading to ``near-Eckart'' and transient regimes associated with `abrupt' and `smooth' relaxation processes, respectively. Assuming the WIMP particle to be a super-symmetric neutralino of a 100$ GeV, the numerical analysis reveals that a physically plausible evolution, compatible with a stable equilibrium state and with observed bounds on CMB anisotropies and neutralino abundance, is only possible for models characterized by initial conditions associated with nearly zero spatial curvature and total initial energy density very close to unity.

We study the asymptotic behaviour of scaling solutions with a dissipative fluid and we show that, contrary to recent claims, the existence of stable accelerating attractor solution which solves the `energy' coincidence problem depends crucially on the chosen equations of state for the thermodynamical variables. We discuss two types of equations of state, one which contradicts this claim, and one which supports it.

We study the qualitative properties of the class of spatially homogeneous Bianchi VI_o cosmological models containing a perfect fluid with a linear equation of state, a scalar field with an exponential potential and a uniform cosmic magnetic field, using dynamical systems techniques. We find that all models evolve away from an expanding massless scalar field model in which the matter and the magnetic field are negligible dynamically. We also find that for a particular range of parameter values the models evolve towards the usual power-law inflationary model (with no magnetic field) and, furthermore, we conclude that inflation is not fundamentally affected by the presence of a uniform primordial magnetic field. We investigate the physical properties of the Bianchi I magnetic field models in some detail.

Observations of the high degree of isotropy of the cosmic microwave background are commonly believed to indicate that the Universe is `almost' Friedmann-Lemaitre-Robertson-Walker (at least since the time of last scattering). Theoretical support for this belief comes from the so-called Ehlers-Geren-Sachs Theorem. We show that a generalization of this theorem rules out any strong magnetic fields in the Universe. Our theoretical result is model-independent and includes the case of inhomogeneous magnetic fields, complementing previous results. We thus prove that cosmic microwave background observations severely constrain all types of primordial and protogalactic magnetic fields in the universe.

We investigate the dynamical properties of a class of spatially homogeneous and isotropic cosmological models containing a barotropic perfect fluid and multiple scalar fields with independent exponential potentials. We show that the assisted inflationary scaling solution is the global late-time attractor for the parameter values for which the model is inflationary, even when curvature and barotropic matter are included. For all other parameter values the multi-field curvature scaling solution is the global late-time attractor (in these solutions asymptotically the curvature is not dynamically negligible). Consequently, we find that in general all of the scalar fields in multi-field models with exponential potentials are non-negligible in late-time behaviour, contrary to what is commonly believed. The early-time and intermediate behaviour of the models is also studied. In particular, n-scalar field models are investigated and the structure of the saddle equilibrium points corresponding to inflationary m-field scaling solutions and non-inflationary m-field matter scaling solutions are also studied (where m<n), leading to interesting transient dynamical behaviour with new physical scenarios of potential importance.

The qualitative properties of spatially homogeneous stiff perfect fluid and minimally coupled massless scalar field models within general relativity are discussed. Consequently, by exploiting the formal equivalence under conformal transformations and field redefinitions of certain classes of theories of gravity, the asymptotic properties of spatially homogeneous models in a class of scalar-tensor theories of gravity that includes the Brans-Dicke theory can be determined. For example, exact solutions are presented, which are analogues of the general relativistic Jacobs stiff perfect fluid solutions and vacuum plane wave solutions, which act as past and future attractors in the class of spatially homogeneous models in Brans-Dicke theory.

Dynamical systems theory is especially well-suited for determining the possible asymptotic states (at both early and late times) of cosmological models, particularly when the governing equations are a finite system of autonomous ordinary differential equations. We begin with a brief review of dynamical systems theory. We then discuss cosmological models as dynamical systems and point out the important role of self-similar models. We review the asymptotic properties of spatially homogeneous perfect fluid models in general relativity. We then discuss some results concerning scalar field models with an exponential potential (both with and without barotropic matter). Finally, we discuss some isotropic cosmological models derived from the string effective action.

We investigate spatially flat isotropic cosmological models which contain a scalar field with an exponential potential and a perfect fluid with a linear equation of state. We include an interaction term, through which the energy of the scalar field is transferred to the matter fields, consistent with a term that arises from scalar--tensor theory under a conformal transformation and field redefinition. The governing ordinary differential equations reduce to a dynamical system when appropriate normalized variables are defined. We analyse the dynamical system and find that the interaction term can significantly affect the qualitative behaviour of the models. The late-time behaviour of these models may be of cosmological interest. In particular, for a specific range of values for the model parameters there are late-time attracting solutions, corresponding to a novel attracting equilib rium point, which are inflationary and in which the scalar field's energy-density remains a fixed fraction of the matter field's energy density. These scalar field models may be of interest as late-time cosmologies, part icularly in view of the recent observations of the current accelerated cosmic expansion. For appropriate values of the interaction coupling parameter, this equilibrium point is an attracting focus, and hence as inflating solutions approach this late-time attractor the scalar field oscillates. Hence these models may also be of importance in the study of inflation in the early universe.

A complete global analysis of spatially-flat, four-dimensional cosmologies derived from the type IIA string and M-theory effective actions is presented. A non--trivial Ramond-Ramond sector is included. The governing equations are written as a dynamical system. Asymptotically, the form fields are dynamically negligible, but play a crucial role in determining the possible intermediate behaviour of the solutions (i.e. the nature of the equilibrium points). The only past-attracting solution (source in the system) may be interpreted in the eleven-dimensional setting in terms of flat space. This source is unstable to the introduction of spatial curvature.

A qualitative analysis is presented for spatially flat, isotropic and homogeneous cosmologies derived from the string effective action when the combined effects of a dilaton, modulus, two-form potential and central charge deficit are included. The latter has significant effects on the qualitative dynamics. The analysis is also directly applicable to the anisotropic Bianchi type I cosmology.

We classify all spherically symmetric perfect fluid solutions of Einstein's equations with equation of state p/mu=a which are self-similar in the sense that all dimensionless variables depend only upon z=r/t. For a given value of a, such solutions are described by two parameters and they can be classified in terms of their behaviour at large and small distances from the origin; this usually corresponds to large and small values of z but (due to a coordinate anomaly) it may also correspond to finite z. We base our analysis on the demonstration that all similarity solutions must be asymptotic to solutions which depend on either powers of z or powers of lnz. We show that there are only three similarity solutions which have an exact power-law dependence on z: the flat Friedmann solution, a static solution and a Kantowski-Sachs solution (although the latter is probably only physical for a<0). For a>1/5, there are also two families of solutions which are asymptotically (but not exactly) Minkowski: the first is asymptotically Minkowski as z tends to infinity and is described by one parameter; the second is asymptotically Minkowski at a finite value of z and is described by two parameters. A complete analysis of the dust solutions is given, since these can be written down explicitly and elucidate the link between the z>0 and z<0 solutions. Solutions with pressure are then discussed in detail; these share many of the characteristics of the dust solutions but they also exhibit new features.

We review recent lattice QCD results for the surface tension at the finite temperature quark-hadron phase transition and discuss their implications on the possible scale of inhomogeneities. In the quenched approximation the average distance between nucleating centers is smaller than the diffusion length of a protron, so that inhomogeneities are washed out by the time nucleosynthesis sets in. Consequently the baryon density fluctuations formed by a QCD phase transition in the early universe cannot significantly affect standard big-bang nucleosynthesis calculations and certainly cannot allow baryons to close the universe. At present lattice results are inconclusive when dynamical fermions are included.