We examine numerically and qualitatively the Lemaître–Tolman–Bondi (LTB) inhomogeneous dust solutions as a three-dimensional dynamical system characterized by six critical points. One of the coordinates of the phase space is an average density parameter, ⟨Ω⟩, which behaves as the ordinary Ω in Friedman–Lemaître–Robertson–Walker (FLRW) dust spacetimes. The other two coordinates, a shear parameter and a density contrast function, convey the effects of inhomogeneity. As long as shell crossing singularities are absent, this phase space is bounded or it can be trivially compactified. This space contains several invariant subspaces which define relevant particular cases, such as ‘parabolic’ evolution, FLRW dust and the Schwarzschild–Kruskal vacuum limit. We examine in detail the phase-space evolution of several dust configurations: a low-density void formation scenario, high-density re-collapsing universes with open, closed and wormhole topologies, a structure formation scenario with a black hole surrounded by an expanding background and the Schwarzschild–Kruskal vacuum case. Solution curves (except regular centers) start expanding from a past attractor (source) in the plane ⟨Ω⟩ = 1, associated with self-similar regime at an initial singularity. Depending on the initial conditions and specific configurations, the curves approach several saddle points as they evolve between this past attractor and two other possible future attractors: perpetually expanding curves terminate at a line of sinks at ⟨Ω⟩ = 0, while collapsing curves reach maximal expansion as ⟨Ω⟩ diverges and end up in a sink that coincides with the past attractor and is also associated with self-similar behavior.
Read full abstract