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.

In a recent paper [Phys. Rev. D{\bf 85}, 023528 (2012)] the authors proposed a generalized Born-Infeld electrodynamics coupled to general relativity which produces a nonsingular bouncing universe. For a magnetic universe the resulting cosmic evolution inevitably interpolates between asymptotic de Sitter states. Here we shall show that (i) the above theory does not have the standard weak field Maxwell limit, (ii) a sudden curvature singularity -- not better than the big bang -- arises, (iii) the speed of sound squared is a negative quantity signaling instability against small perturbations of the background energy density, and that (iv) the conclusion about the inevitability of the asymptotic vacuum regime in a magnetic universe is wrong.

By means of a combined use of the type Ia supernovae and H(z) data tests, together with the study of the asymptotic properties in the equivalent phase space - through the use of the dynamical systems tools - we demonstrate that the bulk viscous matter-dominated scenario is not a good model to explain the accepted cosmological paradigm, at least, under the parametrization of bulk viscosity considered in this paper. The main objection against such scenarios is the absence of conventional radiation and matter-dominated critical points in the phase space of the model. This entails that radiation and matter dominance are not generic solutions of the cosmological equations, so that these stages can be implemented only by means of unique and very specific initial conditions, i. e., of very unstable particular solutions. Such a behavior is in marked contradiction with the accepted cosmological paradigm which requires of an earlier stage dominated by relativistic species, followed by a period of conventional non-relativistic matter domination, during which the cosmic structure we see was formed. Also, we found that the bulk viscosity is positive just until very late times in the cosmic evolution, around z < 1. For earlier epochs it is negative, been in tension with the local second law of thermodynamics.

We study the phase space asymptotics of the so called Veneziano ghost dark energy models. Models where the ghost field's energy density: i) $\rho_{ghost}\propto H$, and ii) $\rho_{ghost}\propto H+H^2$, are investigated. Both, cases with and without additional non-gravitational interaction between cold dark matter and ghost dark energy, are subject to scrutiny. We pay special attention to the choice of phase space variables leading to bounded and compact phase space so that no critical point of physical interest is missing. A rich asymptotic structure is revealed: depending on the kind of non-minimal coupling critical points associated with radiation dominance, matter dominance, cold dark matter/ghost dark energy scaling, and ghost dark energy dominance, are found. Past and future attractors, as well as saddle equilibrium points, are identified in the corresponding phase spaces.

We investigate in detail the asymptotic properties of tachyon cosmology for a broad class of self-interaction potentials. The present approach relies in an appropriate re-definition of the tachyon field, which, in conjunction with a method formerly applied in the bibliography in a different context, allows to generalize the dynamical systems study of tachyon cosmology to a wider class of self-interaction potentials beyond the (inverse) square-law one. It is revealed that independent of the functional form of the potential, the matter-dominated solution and the ultra-relativistic (also matter-dominated) solution, are always associated with equilibrium points in the phase space of the tachyon models. The latter is always the past attractor, while the former is a saddle critical point. For inverse power-law potentials $V\propto\phi^{-2\lambda}$ the late-time attractor is always the de Sitter solution, while for sinh-like potentials $V\propto\sinh^{-\alpha}(\lambda\phi)$, depending on the region of parameter space, the late-time attractor can be either the inflationary tachyon-dominated solution or the matter-scaling (also inflationary) phase. In general, for most part of known quintessential potentials, the late-time dynamics will be associated either with de Sitter inflation, or with matter-scaling, or with scalar field-dominated solutions.