BKL have conjectured that one should expect a generic singularity in spatially inhomogeneous cosmologies to be locally of the Mixmaster type (local Mixmaster dynamics (LMD)) . For a review of homogeneous cosmologies, inhomogeneous cosmologies, and the relation between them, see . The main difficulty with the acceptance of this conjecture has been the controversy over whether the required time slicing can be constructed globally [13, 122]. Montani , Belinskii , and Kirillov and Kochnev [175, 174] have pointed out that if the BKL conjecture is correct, the spatial structure of the singularity could become extremely complicated as bounces occur at different locations at different times. BKL seem to imply  that LMD should only be expected to occur in completely general spacetimes with no spatial symmetries. However, LMD is actually possible in any spatially inhomogeneous cosmology with a local MSS with a “closed” potential (in the sense of the standard triangular potentials of Bianchi VIII and IX). This closure may be provided by spatial curvature, matter fields, or rotation. A class of cosmological models which appear to have local MD are vacuum universes on with a symmetry . Even simpler plane symmetric Gowdy cosmologies [121, 23] have “open” local MSS potentials. However, these models are interesting in their own right since they have been conjectured to possess an AVTD singularity . One way to obtain these Gowdy models is to allow spatial dependence in one direction in Bianchi I homogeneous cosmologies . It is well-known that addition of matter terms to homogeneous Bianchi I, Bianchi VI, and other AVTD models can yield Mixmaster behavior [168, 178, 177]. Allowing spatial dependence in one direction in such models might then yield a spacetime with LMD. Application of this procedure to magnetic Bianchi VI models yields magnetic Gowdy models [250, 247]. Of course, Gowdy cosmologies are not the most general symmetric vacuum spacetimes [121, 82, 34]. Restoring the “twists” introduces a centrifugal wall to close the MSS. Magnetic Gowdy and general symmetric models appear to admit LMD [247, 248, 41].
The past few years have seen the development of strong numerical evidence in support of the BKL claims . The Method of Consistent Potentials (MCP)  has been used to organize the data obtained in simulations of spatially inhomogeneous cosmologies [42, 36, 250, 43, 37, 32, 41]. The main idea is to obtain a Kasner-like velocity term dominated (VTD) solution at every spatial point by solving Einstein’s equations truncated by removing all terms containing spatial derivatives. If the spacetime really is AVTD, all the neglected terms will be subdominant (exponentially small in variables where the VTD solution is linear in the time ) when the VTD solution is substituted back into them. For the MCP to successfully predict whether or not the spacetime is AVTD, the dynamics of the full solution must be dominated at (almost) every spatial point by the VTD solution behavior. Surprisingly, MCP predictions have proved valid in numerical simulations of cosmological spacetimes with one  and two [42, 36, 41] spatial symmetries. In the case of symmetric models, a comparison between the observed behavior  and that in a vacuum, diagonal Bianchi IX model written in terms of variables provides strong support for LMD  since the phenomenology of the inhomogeneous cosmologies can be reproduced by this rewriting of the standard Bianchi IX MD.
Polarized plane symmetric cosmologies have been evolved numerically using standard techniques by Anninos, Centrella, and Matzner [5, 6]. The long-term project involving Berger, Garfinkle, and Moncrief and their students to study the generic cosmological singularity numerically has been reviewed in detail elsewhere [38, 37, 31] but will be discussed briefly here.
The Gowdy model on is described by gravitational wave amplitudes and which propagate in a spatially inhomogeneous background universe described by . (We note that the physical behavior of a Gowdy spacetime can be computed from the effect of the metric evolution on a test cylinder .) We impose and periodic boundary conditions. The time variable measures the area in the symmetry plane with being a curvature singularity.
Einstein’s equations split into two groups. The first is nonlinearly coupled wave equations for dynamical variables and (where ) obtained from the variation of 
For the special case of the polarized Gowdy model (), satisfies a linear wave equation whose exact solution is well-known . For this case, it has been proven that the singularity is AVTD . This has also been conjectured to be true for generic Gowdy models .
AVTD behavior is defined in  as follows: Solve the Gowdy wave equations neglecting all terms containing spatial derivatives. This yields the VTD solution . If the approach to the singularity is AVTD, the full solution comes arbitrarily close to a VTD solution at each spatial point as . As , the VTD solution becomes. We have shown that, except at isolated spatial points, the nonlinear terms in the wave equation for drive into this range [36, 38]. The exceptional points occur when coefficients of the nonlinear terms vanish and are responsible for the growth of spiky features seen in the wave forms [42, 36]. We conclude that generic Gowdy cosmologies have an AVTD singularity except at isolated spatial points [36, 38]. This has been confirmed by Hern and Stuart  and by van Putten . After the nature of the solutions became clear through numerical experiments, it became possible to use Fuchsian asymptotic methods to prove that Gowdy solutions with and AVTD behavior almost everywhere are generic . These methods have recently been applied to Gowdy spacetimes with and topologies with similar conclusions .
One striking property of the Gowdy models are the development of “spiky features” at isolated spatial points where the coefficient of a local “potential term” vanishes [42, 36]. Recently, Rendall and Weaver have shown analytically how to generate such spikes from a Gowdy solution without spikes .
Addition of a magnetic field to the vacuum Gowdy models (plus a topology change) which yields the inhomogeneous generalization of magnetic Bianchi VI models provides an additional potential which grows exponentially if . Local Mixmaster behavior has recently been observed in these magnetic Gowdy models [250, 247].
Garfinkle has used a vacuum Gowdy model with spatial topology to test an algorithm for axis regularity . Along the way, he has shown that these models are also AVTD with behavior at generic spatial points that is eventually identical to that in the case. Comparison of the two models illustrates that topology or other global or boundary conditions are important early in the simulation but become irrelevant as the singularity is approached.
Gowdy spacetimes are not the most general symmetric vacuum cosmologies. Certain off-diagonal metric components (the twists which are , in the notation of (12)) have been set to zero . Restoring these terms (see [83, 34]) yields spacetimes that are not AVTD but rather appear to exhibit a novel type of LMD [41, 249]. The LMD in these models is an inhomogeneous generalization of non-diagonal Bianchi models with “centrifugal” MSS potential walls [227, 169] in addition to the usual curvature walls. In , remarkable quantitative agreement is found between predictions of the MCP and numerical simulation of the full Einstein equations. A version of the code with AMR has been developed . (Asymptotic methods have been used to prove that the polarized version of these spacetimes have AVTD solutions .)
Moncrief has shown  that cosmological models on with a spatial symmetry can be described by five degrees of freedom and their respective conjugate momenta . All variables are functions of spatial variables , and time . Einstein’s equations can be obtained by variation of is used. A spectral evaluation of derivatives  which has been shown to work in Gowdy simulations  does not appear to be helpful in the case.) A particular solution to the initial value problem is used since the general solution is not available .
Current limitations of the code do not affect simulations for the polarized case since problematic spiky features do not develop. Polarized models have . Grubisić and Moncrief  have conjectured that these polarized models are AVTD. The numerical simulations provide strong support for this conjecture [38, 44]. Asymptotic methods have been used to prove that an open set of AVTD solutions exist for this case .
The MCP indicates that the term containing gradients of in (14) acts as a Mixmaster-like potential to drive the system away from AVTD behavior in generic models . Numerical simulations provide support for this suggestion [38, 43]. Whether this potential term grows or decays depends on a function of the field momenta. This in turn is restricted by the Hamiltonian constraint. However, failure to enforce the constraints can cause an erroneous relationship among the momenta to yield qualitatively wrong behavior. There is numerical evidence that this error tends to suppress Mixmaster-like behavior leading to apparent AVTD behavior in extended spatial regions of symmetric cosmologies [27, 28]. In fact, it has been found , that when the Hamiltonian constraint is enforced at every time step, the predicted local oscillatory behavior of the approach to the singularity is observed. (The momentum constraint is not enforced.) (Note that in a numerical study of vacuum Bianchi IX homogeneous cosmologies, Zardecki obtained a spurious enhancement of Mixmaster oscillations due to constraint violation [253, 147]. In this case, the constraint violation introduced negative energy.)
Mixmaster simulations with the new algorithm  can easily evolve more than 250 bounces reaching . This compares to earlier simulations yielding 30 or so bounces with . The larger number of bounces quickly reveals that it is necessary to enforce the Hamiltonian constraint. An explicitly constraint enforcing code was developed some years ago by Ove (see  and references therein).
It is well known  that a scalar field can suppress Mixmaster oscillations in homogeneous cosmologies. BKL argued that the suppression would also occur in spatially inhomogeneous models. This was demonstrated numerically for magnetic Gowdy and symmetric spacetimes . Andersson and Rendall proved that completely general cosmological (spatially ) spacetimes (no symmetries) with sufficiently strong scalar fields have generic AVTD solutions . Garfinkle  has constructed a 3D harmonic code which, so far, has found AVTD solutions with a scalar field present. Work on generic vacuum models is in progress.
Cosmological models inspired by string theory contain higher derivative curvature terms and exotic matter fields. Damour and Henneaux have applied the BKL approach to such models and conclude that their approach to the singularity exhibits LMD .
Finally, there has been a study of the relationship between the “long wavelength approximation” and the BKL analyses by Deruelle and Langlois .
© Max Planck Society and the author(s)