Abstract. Quasi-local scalar variables approach is applied to a spherically symmetric inhomogeneous Lemaître-Tolman-Bondi metric containing a mixture of non-relativistic cold dark matter and coupled dark energy with constant equation of state. The quasi-local coupling term considered is proportional to the quasi-local cold dark matter energy density and a quasi-local Hubble factor-like scalar via a coupling constant α. The autonomous numerical system obtained from the evolution equations is classified for different choices of the free parameters: the adiabatic constant of the dark energy w and α. The presence of a past attractor in a non-physical region of the energy densities phase-space of the system makes the coupling term non physical when the energy flows from the matter to the dark energy in order to avoid negative values of the dark energy density in the past. On the other hand, if the energy flux goes from dark energy to dark matter, the past attractor lays in a physical region. The system is also numerically solved for some interesting initial profiles leading to different configurations: an ever expanding mixture, a scenario where the dark energy is completely consumed by the nonrelativistic matter by means of the coupling term, a scenario where the dark energy disappears in the inner layers while the outer layers expand as a mixture of both sources, and, finally, a structure formation toy model scenario, where the inner shells containing the mixture collapse while the outer shells expand.PACS numbers: 98.80.-k, 04.20.-q, 95.36.+x, 95.35.+d

We examine the evolution of an inhomogeneous mixture of non-relativistic pressureless cold dark matter (CDM), coupled to dark energy (DE) characterised by the equation of state parameter w < −1/3, with the interaction term proportional to the DE density. This coupled mixture is the source of a spherically symmetric Lemaître-Tolman-Bondi (LTB) metric admitting an asymptotic Friedman-Lemaître-Robertson-Walker (FLRW) background. Einstein's equations reduce to a 5-dimensional autonomous dynamical system involving quasi-local variables related to suitable averages of covariant scalars and their fluctuations. The phase space evolution around the critical points (past/future attractors and five saddles) is examined in detail. For all parameter values and both directions of energy flow (CDM to DE and DE to CDM) the phase space trajectories are compatible with a physically plausible early cosmic times behaviour near the past attractor. This result compares favourably with mixtures with interaction driven by the CDM density, whose past evolution is unphysical for DE to CDM energy flow. Numerical examples are provided describing the evolution of an initial profile that can be associated with idealised structure formation scenarios.

In this paper, we consider inhomogeneous spherically symmetric models based on the Lemaître–Tolman–Bondi (LTB) metric, assuming as its source an interactive mixture of ordinary baryonic matter, cold dark matter and dark energy with a coupling term proportional to the addition of energy densities of both dark fluids. We reduce Einstein’s field equations to a first-order seven-dimensional autonomous dynamical system of evolution equations and algebraic constraints. We study in detail the evolution of the energy density and spatial curvature profiles along the phase–space by means of two subspace projections: a three-dimensional projection associated with the solutions of the Friedman–Lemaître–Robertson–Walker metric (invariant subspace) and a four-dimensional projection describing the evolution of the inhomogeneous fluctuations. We also classify and study the critical points of the system in comparison with previous work on similar sources, as well as solving numerically the equations for initial energy density and curvature profiles that lead to a spherical bounce whose collapsing time we estimate appropriately.