The quantum field theoretic prediction for the vacuum energy density leads to a value for the effective cosmological constant that is incorrect by between 60 to 120 orders of magnitude. We review an old proposal of replacing Einstein's Field Equations by their trace-free part (the Trace-Free Einstein Equations), together with an independent assumption of energy--momentum conservation by matter fields. While this does not solve the fundamental issue of why the cosmological constant has the value that is observed cosmologically, it is indeed a viable theory that resolves the problem of the discrepancy between the vacuum energy density and the observed value of the cosmological constant. However, one has to check that, as well as preserving the standard cosmological equations, this does not destroy other predictions, such as the junction conditions that underlie the use of standard stellar models. We confirm that no problems arise here: hence, the Trace-Free Einstein Equations are indeed viable for cosmological and astrophysical applications.

The thoughts expressed in this article are based on remarks made by J\"urgen Ehlers at the Albert-Einstein-Institut, Golm, Germany in July 2007. The main objective of this article is to demonstrate, in terms of plausible order-of-magnitude estimates for geometrical scalars, the relevance of spatial curvature in realistic models of the Universe that describe the dynamics of structure formation since the epoch of matter-radiation decoupling. We introduce these estimates with a commentary on the use of a quasi-Newtonian metric form in this context.

In this paper we investigate asymptotic isotropization. We derive the asymptotic dynamics of spatially inhomogeneous cosmological models with a perfect fluid matter source and a positive cosmological constant near the de Sitter equilibrium state at late times, and near the flat FL equilibrium state at early times. Our results show that there exists an open set of solutions approaching the de Sitter state at late times, consistent with the cosmic no-hair conjecture. On the other hand, solutions that approach the flat FL state at early times are special and admit a so-called isotropic initial singularity. For both classes of models the asymptotic expansion of the line element contains an arbitrary spatial metric at leading order, indicating asymptotic spatial inhomogeneity. We show, however, that in the asymptotic regimes this spatial inhomogeneity is significant only at super-horizon scales.

We present a general framework for analyzing spatially inhomogeneous cosmological dynamics. It employs Hubble-normalized scale-invariant variables which are defined within the orthonormal frame formalism, and leads to the formulation of Einstein's field equations with a perfect fluid matter source as an autonomous system of evolution equations and constraints. This framework incorporates spatially homogeneous dynamics in a natural way as a special case, thereby placing earlier work on spatially homogeneous cosmology in a broader context, and allows us to draw on experience gained in that field using dynamical systems methods. One of our goals is to provide a precise formulation of the approach to the spacelike initial singularity in cosmological models, described heuristically by Belinski\v{\i}, Khalatnikov and Lifshitz. Specifically, we construct an invariant set which we conjecture forms the local past attractor for the evolution equations. We anticipate that this new formulation will provide the basis for proving rigorous theorems concerning the asymptotic behavior of spatially inhomogeneous cosmological models.

We give a careful general relativistic and (1+3)-covariant analysis of cosmological peculiar velocities induced by matter density perturbations in the presence of a cosmological constant. In our quasi-Newtonian approach, constraint equations arise to maintain zero shear of the non-comoving fundamental worldlines which define a Newtonian-like frame, and these lead to the (1+3)-covariant dynamical equations, including a generalized Poisson-type equation. We investigate the relation between peculiar velocity and peculiar acceleration, finding the conditions under which they are aligned. In this case we find (1+3)-covariant relativistic generalizations of well-known Newtonian results.

We employ the extended 1+3 orthonormal frame formalism for fluid spacetime geometries $({\cal M}, {\bf g}, {\bf u})$, which contains the Bianchi field equations for the Weyl curvature, to derive a 44-D evolution system of first-order symmetric hyperbolic form for a set of geometrically defined dynamical field variables. Describing the matter source fields phenomenologically in terms of a barotropic perfect fluid, the propagation velocities $v$ (with respect to matter-comoving observers that Fermi-propagate their spatial reference frames) of disturbances in the matter and the gravitational field, represented as wavefronts by the characteristic 3-surfaces of the system, are obtained. In particular, the Weyl curvature is found to account for two (non-Lorentz-invariant) Coulomb-like characteristic eigenfields propagating with $v = 0$ and four transverse characteristic eigenfields propagating with $|v| = 1$, which are well known, and four (non-Lorentz-invariant) longitudinal characteristic eigenfields propagating with $|v| = \sfrac{1}{2}$. The implications of this result are discussed in some detail and a parallel is drawn to the propagation of irregularities in the matter distribution. In a worked example, we specialise the equations to cosmological models in locally rotationally symmetric class II and include the constraints into the set of causally propagating dynamical variables.

Exact dynamical equations for a generic dust matter source field in a cosmological context are formulated with respect to a non-comoving Newtonian-like timelike reference congruence and investigated for internal consistency. On the basis of a lapse function $N$ (the relativistic acceleration scalar potential) which evolves along the reference congruence according to $\dot{N} = \alpha \Theta N$ ($\alpha = {const}$), we find that consistency of the quasi-Newtonian dynamical equations is not attained at the first derivative level. We then proceed to show that a self-consistent set can be obtained by linearising the dynamical equations about a (non-comoving) FLRW background. In this case, on properly accounting for the first-order momentum density relating to the non-relativistic peculiar motion of the matter, additional source terms arise in the evolution and constraint equations describing small-amplitude energy density fluctuations that do not appear in similar gravitational instability scenarios in the standard literature.

We revisit the issue of integrability conditions for the irrotational silent cosmological models. We formulate the problem both in 1+3 covariant and 1+3 orthonormal frame notation, and show there exists a series of constraint equations that need to be satisfied. These conditions hold identically for FLRW-linearised silent models, but not in the general exact non-linear case. Thus there is a linearisation instability, and it is highly unlikely that there is a large class of silent models. We conjecture that there are no spatially inhomogeneous solutions with Weyl curvature of Petrov type I, and indicate further issues that await clarification.

After a brief overview of the so-called silent models and their present status, we consider the subclass of Bianchi Type--I models with a magnetic field source. Due to the presence of the magnetic field, the initial singularity shows ``oscillatory'' features reminiscent of the Bianchi Type--IX case. The Bianchi Type--I models with a magnetic field are therefore a counterexample to the folklore that matter fields can be neglected in the vicinity of the singularity.

The dynamics of perfect fluid spacetime geometries which exhibit {\em Local Rotational Symmetry} (LRS) are reformulated in the language of a $1+\,3$ "threading" decomposition of the spacetime manifold, where covariant fluid and curvature variables are used. This approach presents a neat alternative to the orthonormal frame formalism. The dynamical equations reduce to a set of differential relations between purely scalar quantities. The consistency conditions are worked out in a transparent way. We discuss their various subcases in detail and focus in particular on models with higher symmetries within the class of expanding spatially inhomogeneous LRS models, via a consideration of functional dependencies between the dynamical variables.