Singularity-free cosmology from interactions in the dark sector
We study the dynamics of Friedmann-Lemaître-Robertson-Walker models where a dark energy component with a quadratic equation of state (EoS) nonlinearly interacts with cold dark matter. Thus, two energy scales naturally come into play: $ρ_*$ is the scale at which the nonlinearity of the EoS becomes relevant; $ρ_i$ is the energy scale around which the interaction starts to play an important role in the dynamics. Our focus is to understand whether there are parameter ranges for this system that can produce non-singular bouncing and emergent cosmologies for any initial condition. We complete a dynamical systems analysis, and find the parameter range such that trajectories always expand from a high energy non-singular de Sitter state. For flat and negative curvature models this de Sitter state is represented by a fixed point, the asymptotic past from which the universe emerges from. We find a subset of positive curvature models that during contraction get arbitrarily close to the de Sitter state, thus having a quasi-de Sitter bounce, then emerge from the bounce and expand, evolving toward spatial flatness. We find that the dimensionless parameter $q\equiv ρ_*/ρ_i$, which measures the relative strength of the nonlinear terms in the system, plays a crucial role in the topology of the phase space. When $q < 3$, some trajectories expand toward a singularity, while others evolve toward a low energy cosmological constant at late-times, with a subset going through a decelerated matter dominate era before the final acceleration. When $q > 3$, all trajectories are non-singular, and evolve toward a late-time cosmological constant. We find a subclass in this case in which all trajectories have at least one decelerated matter dominated phase, and accelerate at late-times.