The theory of the dynamical systems is a very complex subject which has brought several surprises in the recent past in connection with the theory of chaos and fractals. The application of the tools of the dynamical systems in cosmological settings is less known in spite of the amount of published scientific papers on this subject. In this paper a -- mostly pedagogical -- introduction to the application in cosmology of the basic tools of the dynamical systems theory is presented. It is shown that, in spite of their amazing simplicity, these allow to extract essential information on the asymptotic dynamics of a wide variety of cosmological models. The power of these tools is illustrated within the context of the so called $Lambda$CDM and scalar field models of dark energy. This paper is suitable for teachers, undergraduate and postgraduate students from physics and mathematics disciplines.
Here we investigate the cosmic dynamics of Friedmann-Robertson-Walker universes -- flat spatial sections -- which are driven by nonlinear electrodynamics (NLED) Lagrangians. We pay special attention to the check of the sign of the square sound speed since, whenever the latter quantity is negative, the corresponding cosmological model is classically unstable against small perturbations of the background energy density. Besides, based on causality arguments, one has to require that the mentioned small perturbations of the background should propagate at most at the local speed of light. We also look for the occurrence of curvature singularities. Our results indicate that several cosmological models which are based in known NLED Lagrangians, either are plagued by curvature singularities of the sudden and/or big rip type, or are violently unstable against small perturbations of the cosmological background -- due to negative sign of the square sound speed -- or both. In addition, causality issues associated with superluminal propagation of the background perturbations may also arise.
