Why is the Universe so homogeneous and isotropic? We summarize a general study of a $\gamma$-law perfect fluid alongside an inhomogeneous, massless scalar gauge field (with homogeneous gradient) in anisotropic spaces with General Relativity. The anisotropic matter sector is implemented as a $j$-form (field-strength level), where $j\,\in\,\{1,3\}$, and the spaces studied are Bianchi space-times of solvable type. Wald's no-hair theorem is extended to include the $j$-form case. We highlight three new self-similar space-times: the Edge, the Rope and Wonderland. The latter solution is so far found to exist in the physical state space of types I,II, IV, VI$_0$, VI$_h$, VII$_0$ and VII$_h$, and is a global attractor in I and V. The stability analysis of the other types has not yet been performed. This paper is a summary of ~[1], with some remarks towards new results which will be further laid out in upcoming work.

In this paper the dynamics of free gauge fields in Bianchi type I-VII$_{h}$ space-times is investigated. The general equations for a matter sector consisting of a $p$-form field strength ($p\,\in\,\{1,3\}$), a cosmological constant ($4$-form) and perfect fluid in Bianchi type I-VII$_{h}$ space-times are computed using the orthonormal frame method. The number of independent components of a $p$-form in all Bianchi types I-IX are derived and, by means of the dynamical systems approach, the behaviour of such fields in Bianchi type I and V are studied. Both a local and a global analysis are performed and strong global results regarding the general behaviour are obtained. New self-similar cosmological solutions appear both in Bianchi type I and Bianchi type V, in particular, a one-parameter family of self-similar solutions,"Wonderland ($\lambda$)" appears generally in type V and in type I for $\lambda=0$. Depending on the value of the equation of state parameter other new stable solutions are also found ("The Rope" and "The Edge") containing a purely spatial field strength that rotates relative to the co-moving inertial tetrad. Using monotone functions, global results are given and the conditions under which exact solutions are (global) attractors are found.

We apply the causal Israel-Stewart theory of irreversible thermodynamics to model the matter content of the universe as a dissipative fluid with bulk and shear viscosity. Along with the full transport equations we consider their widely used truncated version. By implementing a dynamical systems approach to Bianchi type IV and V cosmological models with and without cosmological constant, we determine the future asymptotic states of such universes and show that the truncated Israel-Stewart theory leads to solutions essentially different from the full theory. The solutions of the truncated theory may also manifest unphysical properties. Finally, we find that in the full theory shear viscosity can give a substantial rise to dissipative fluxes, driving the fluid extremely far from equilibrium, where the linear Israel-Stewart theory ceases to be valid.

We investigate the three types of class B Bianchi cosmologies filled with a tilted perfect fluid undergoing velocity diffusion in a scalar field background. We consider the two most importantcases: dust and radiation. A complete numerical integration of the Einstein field equations coupled with the diffusion equations is done to demonstrate how the presence of diffusion can affect the dynamics of cosmological evolution, where the most attention is paid to changes to the late-time behaviour. We show that aside from quantitative effects, diffusion can result in significant qualitative differences. For example, the cosmologies may recollapse if diffusion is sufficiently strong, or evolve towards the de Sitter state otherwise. In constrast to the diffusionless case, radiation isotropizes in presence of diffusion, and the tilt decreases exponentially at later times: $V\sim e^{-0.25\tau}$; also, we determine the decay rates of energy density, which become slower when the diffusion term is non-zero.

We apply the dynamical systems approach to ever-expanding Bianchi type VIII cosmologies filled with a tilted $\gamma$-fluid undergoing velocity diffusion on a scalar field. We determine the future attractors and investigate the late-time behaviour of the models. We find that at late times the normalized energy density $\Omega$ tends to zero, while the scalar potential $\Phi$ approaches 1 and dominates the evolution. Moreover, we demonstrate that in presence of diffusion fluids with $\gamma<3/2$, which includes physically important cases of dust $(\gamma=1)$ and radiation $(\gamma=4/3)$, are asymptotically non-tilted; the velocity of the fluid with $\gamma=3/2$ tends to a constant value $0<\bar{V}<1$; and stiffer fluids evolve towards a state of extreme tilt. Finally, we show that diffusion significantly reduces the decay rates of energy density for dust and fluids stiffer than dust $(\gamma \geq 1)$; for example, at $\gamma=4/3$ (radiation) we obtain $\rho/H^2 \propto e^{-3H_0 t}$ at late times, while $\rho/H^2 \propto e^{-4H_0 t}$ when diffusion is absent.

We investigate a simple inhomogeneous anisotropic cosmology (plane symmetric $G_2$ model) filled with a tilted perfect fluid undergoing velocity diffusion on a scalar field. Considered are two types of fluid: dust and radiation. We solve the system of Einstein field equations and diffusion equations numerically and demonstrate how the universe evolves towards its future asymptotic state. Also, typical time scales of characteristic processes are determined. The obtained results for dust- and radiation-filled cosmologies are compared to each other and to those in the diffusionless case, giving a hint on which effects can be the result of including diffusion in more complicated models. For example, diffusion causes the accelerated expansion stage to arrive at later times.

Motivated by the couplings of the dilaton in four-dimensional effective actions, we investigate the cosmological consequences of a scalar field coupled both to matter and a Maxwell-type vector field. The vector field has a background isotropy-violating component. New anisotropic scaling solutions which can be responsible for the matter and dark energy dominated epochs are identified and explored. For a large parameter region the universe expands almost isotropically. Using that the CMB quadrupole is extremely sensitive to shear, we constrain the ratio of the matter coupling to the vector coupling to be less than 10^(-5). Moreover, we identify a large parameter region, corresponding to a strong vector coupling regime, yielding exciting and viable cosmologies close to the LCDM limit.

Recently an inflationary model with a vector field coupled to the inflaton was proposed and the phenomenology studied for the Bianchi type I spacetime. It was found that the model demonstrates a counter-example to the cosmic no-hair theorem since there exists a stable anisotropically inflationary fix-point. One of the great triumphs of inflation, however, is that it explains the observed flatness and isotropy of the universe today without requiring special initial conditions. Any acceptable model for inflation should thus explain these observations in a satisfactory way. To check whether the model meets this requirement, we introduce curvature to the background geometry and consider axisymmetric spacetimes of Bianchi type II,III and the Kantowski-Sachs metric. We show that the anisotropic Bianchi type I fix-point is an attractor for the entire family of such spacetimes. The model is predictive in the sense that the universe gets close to this fix-point after a few e-folds for a wide range of initial conditions. If inflation lasts for N e-folds, the curvature at the end of inflation is typically of order exp(-2N). The anisotropy in the expansion rate at the end of inflation, on the other hand, while being small on the one-percent level, is highly significant. We show that after the end of inflation there will be a period of isotropization lasting for about 2N/3 e-folds. After that the shear scales as the curvature and becomes dominant around N e-folds after the end of inflation. For plausible bounds on the reheat temperature the minimum number of e-folds during inflation, required for consistency with the isotropy of the supernova Ia data, lays in the interval (21,48). Thus the results obtained for our restricted class of spacetimes indicates that inflation with anisotropic hair is cosmologically viable.

We display some simple cosmological solutions of gravity theories with quadratic Ricci curvature terms added to the Einstein-Hilbert lagrangian which exhibit anisotropic inflation. The Hubble expansion rates are constant and unequal in three orthogonal directions. We describe the evolution of the simplest of these homogeneous and anisotropic cosmological models from its natural initial state and evaluate the deviations they will create from statistical isotropy in the fluctuations produced during a period of anisotropic inflation. The anisotropic inflation is not a late-time attractor in these models but the rate of approach to a final isotropic de Sitter state is slow and is conducive to the creation of observable anisotropic statistical effects in the microwave background. The statistical anisotropy would not be scale invariant and the level of statistical anisotropy will grow with scale.

We consider the perihelion precession and bending of light in a class of Kaluza-Klein models and show that the "electric redshift" model, proposed in Zhang (2006) to explain the redshift of Quasars, does not agree with observations. As Zhang's model only considers the Jordan frame, we also compute the perihelion precession as seen in the Einstein frame and show that, to lowest order, the perihelion precession can only be consistent with observation in the Einstein frame. However, when we consider the corresponding "electric redshift" in the Einstein frame, we find that the redshift is significantly lower than for the Jordan frame and is therefore very unlikely to explain the redshift of quasars.

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.

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 use a dynamical systems approach to investigate Bianchi type I and II universes in quadratic theories of gravity. Due to the complicated nature of the equations of motion we focus on the stability of exact solutions and find that there exists an isotropic FRW universe acting as a past attractor. This may indicate that there is an isotropisation mechanism at early times for these kind of theories. We also discuss the Kasner universes, elucidate the associated centre manifold structure, and show that there exists a set of non-zero measure which has the Kasner solutions as a past attractor. Regarding the late-time behaviour, the stability shows a dependence of the parameters of the theory. We give the conditions under which the de Sitter solution is stable and also show that for certain values of the parameters there is a possible late-time behaviour with phantom-like behaviour. New types of anisotropic inflationary behaviour are found which do not have counterparts in general relativity.

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 use the dynamical systems approach to investigate the Bianchi type VIII models with a tilted $\gamma$-law perfect fluid. We introduce expansion-normalised variables and investigate the late-time asymptotic behaviour of the models and determine the late-time asymptotic states. For the Bianchi type VIII models the state space is unbounded and consequently, for all non-inflationary perfect fluids, one of the curvature variables grows without bound. Moreover, we show that for fluids stiffer than dust ($1<\gamma<2$), the fluid will in general tend towards a state of extreme tilt. For dust ($\gamma=1$), or for fluids less stiff than dust ($0<\gamma< 1$), we show that the fluid will in the future be asymptotically non-tilted. Furthermore, we show that for all $\gamma\geq 1$ the universe evolves towards a vacuum state but does so rather slowly, $\rho/H^2\propto 1/\ln t$.

We show that in theories of gravity that add quadratic curvature invariants to the Einstein-Hilbert action there exist expanding vacuum cosmologies with positive cosmological constant which do not approach the de Sitter universe. Exact solutions are found which inflate anisotropically. This behaviour is driven by the Ricci curvature invariant and has no counterpart in the general relativistic limit. These examples show that the cosmic no-hair theorem does not hold in these higher-order extensions of general relativity and raises new questions about the ubiquity of inflation in the very early universe and the thermodynamics of gravitational fields.

We generalize the predictions for the CMB anisotropy patterns arising in Bianchi type VIIh universes to include a dark energy component. We consider these models in light of the result of Jaffe et al. (2005a,b) in which a correlation was found on large angular scales between the WMAP data and the anisotropy structure in a low density Bianchi universe. We find that by including a term Omega_L > 0, the same best-fit anisotropy pattern is reproduced by several combinations of cosmological parameters. This sub-set of models can then be further constrained by current observations that limit the values of various cosmological parameters. In particular, we consider the so-called geometric degeneracy in these parameters imposed by the peak structure of the WMAP data itself. Apparently, despite the additional freedom allowed by the dark energy component, the modified Bianchi models are ruled out at high significance.

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.

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.

We show that a recent comment where Palle criticises works regarding CMBR fluctuations for a particular type of Bianchi models, is incorrect. In particular, in contrast to the claims of Palle, we point out that Bianchi models do allow for vorticity.

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 study the asymptotic behaviour of the Bianchi type VI$_0$ universes with a tilted $\gamma$-law perfect fluid. The late-time attractors are found for the full 7-dimensional state space and for several interesting invariant subspaces. In particular, it is found that for the particular value of the equation of state parameter, $\gamma=6/5$, there exists a bifurcation line which signals a transition of stability between a non-tilted equilibrium point to an extremely tilted equilibrium point. The initial singular regime is also discussed and we argue that the initial behaviour is chaotic for $\gamma<2$.

We consider the stability of spatially homogeneous plane-wave spacetimes. We carry out a full analysis for plane-wave spacetimes in (4+1) dimensions, and find there are two cases to consider; what we call non-exceptional and exceptional. In the non-exceptional case the plane waves are stable to (spatially homogeneous) vacuum perturbations as well as a restricted set of matter perturbations. In the exceptional case we always find an instability. Also we consider the Milne universe in arbitrary dimensions and find it is also stable provided the strong energy condition is satisfied. This implies that there exists an open set of stable plane-wave solutions in arbitrary dimensions.

In this paper we derive homogeneous vacuum plane-wave solutions to Einstein's field equations in 4+1 dimensions. The solutions come in five different types of which three generalise the vacuum plane-wave solutions in 3+1 dimensions to the 4+1 dimensional case. By doing a Kaluza-Klein reduction we obtain solutions to the Einstein-Maxwell equations in 3+1 dimensions. The solutions generalise the vacuum plane-wave spacetimes of Bianchi class B to the non-vacuum case and describe spatially homogeneous spacetimes containing an extremely tilted fluid. Also, using a similar reduction we obtain 3+1 dimensional solutions to the Einstein equations with a scalar field.

An asymptotic stability analysis of spatially homogeneous models of Bianchi type containing tilted perfect fluids is performed. Using the known attractors for the non-tilted Bianchi type universes, we check whether they are stable against perturbations with respect to tilted perfect fluids. We perform the analysis for all Bianchi class B models and the Bianchi type VI_0 model. In particular, we find that none of the non-tilted equilibrium points are stable against tilted perfect fluids stiffer than radiation. We also indicate parts of the phase space where new tilted exact solutions might be found.

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 show that generic anisotropic universes arbitrarily close to the open Friedmann universe allow information processing to continue into the infinite future if there is no cosmological constant or stable gravitationally repulsive stress, and the spatial topology is non-compact. An infinite amount of information can be processed by ``civilisations'' who harness the temperature gradients created by gravitational tidal energy. These gradients are driven by the gravitational waves that sustain the expansion shear and three-curvature anisotropy.

The early time behaviour of brane-world models is analysed in the presence of anisotropic stresses. It is shown that that the initial singularity cannot be isotropic, unless there is also an isotropic fluid stiffer than radiation present. Also, a magnetic Bianchi type I brane-world is analysed in detail. It is known that the Einstein equations for the magnetic Bianchi type I models are in general oscillatory and are believed to be chaotic, but in the brane-world model this chaotic behaviour does not seem to be possible.

In this paper we classify all 4+1 cosmological models where the spatial hypersurfaces are connected and simply connected homogeneous Riemannian manifolds. These models come in two categories, multiply transitive and simply transitive models. There are in all five different multiply transitive models which cannot be considered as a special case of a simply transitive model. The classification of simply transitive models, relies heavily upon the classification of the four dimensional (real) Lie algebras. For the orthogonal case, we derive all the equations of motion and give some examples of exact solutions. Also the problem of how these models can be compactified in context with the Kaluza-Klein mechanism, is addressed.

We study the evolution of the Weyl curvature invariant in all spatially homogeneous universe models containing a non-tilted gamma-law perfect fluid. We investigate all the Bianchi and Thurston type universe models and calculate the asymptotic evolution of Weyl curvature invariant for generic solutions to the Einstein field equations. The influence of compact topology on Bianchi types with hyperbolic space sections is also considered. Special emphasis is placed on the late-time behaviour where several interesting properties of the Weyl curvature invariant occur. The late-time behaviour is classified into five distinctive categories. It is found that for a large class of models, the generic late-time behaviour the Weyl curvature invariant is to dominate the Ricci invariant at late times. This behaviour occurs in universe models which have future attractors that are plane-wave spacetimes, for which all scalar curvature invariants vanish. The overall behaviour of the Weyl curvature invariant is discussed in relation to the proposal that some function of the Weyl tensor or its invariants should play the role of a gravitational 'entropy' for cosmological evolution. In particular, it is found that for all ever-expanding models the measure of gravitational entropy proposed by Gron and Hervik increases at late times.

In this paper we review Penrose's Weyl curvature conjecture which states that the concept of gravitational entropy and the Weyl tensor is somehow linked, at least in a cosmological setting. We give a description of a certain entity constructed from the Weyl tensor, from the very early history of our universe until the present day. Inflation is an important mechanism in our early universe for homogenisation and isotropisation, and thus it must cause large effects upon the evolution of the gravitational entropy. Therefore the effects from inflationary fluids and a cosmological constant are studied in detail.

We investigate brane-worlds with a pure magnetic field and a perfect fluid. We extend earlier work to brane-worlds, and find new properties of the Bianchi type I brane-world. We find new asymptotic behaviours on approach to the singularity and classify the critical points of the dynamical phase space. It is known that the Einstein equations for the magnetic Bianchi type I models are in general oscillatory and are believed to be chaotic, but in the brane-world model this chaotic behaviour does not seem to be possible.