In a recent article, Kleidis and Spyrou (2015) proposed that both dark matter (DM) and dark energy (DE) can be treated as a single component, if accommodated in the context of a polytropic DM fluid with thermodynamical content. Depending only on the polytropic exponent, $-0.103 < \Gamma \leq 0$, this unified DM model reproduces to high accuracy the distance measurements performed with the aid of the supernovae Type Ia (SNe Ia) standard candles, without suffering either from the age or from the coincidence problem. To demonstrate also its compatibility with current observational data concerning structure formation, in the present article we discuss the evolution of cosmological perturbations in the $\Lambda$ CDM-like (i.e., $\Gamma = 0$) limit of the polytropic DM model. The corresponding results are quite encouraging, since, such a model reproduces every major effect already known from conventional (i.e., pressureless cold dark matter - CDM) structure formation theory, such as the constancy of metric perturbations in the vicinity of recombination and the (late-time) Meszaros effect on their rest-mass density counterparts (Meszaros 1974). The non-zero (polytropic) pressure, on the other hand, drives the evolution of small-scale velocity perturbations along the lines of the root-mean-square velocity law of conventional Statistical Physics. As a consequence, in this model "peculiar velocities" slightly increase, instead of being redshifted away by cosmic expansion. What is more important is that, upon consideration of scale-invariant metric perturbations, the spectrum of their rest-mass density counterparts exhibits an effective power-law dependence on the (physical) wavenumber, with the associated scalar spectral index being equal to $n_s = 0.970$; a theoretical value that actually reproduces the corresponding observational (Planck) result (Ade et al. 2016).

In this article, we review a series of recent theoretical results regarding a conventional approach to the dark energy (DE) concept. This approach is distinguished among others for its simplicity and its physical relevance. By compromising General Relativity (GR) and Thermodynamics at cosmological scale, we end up with a model without DE. Instead, the Universe we are proposing is filled with a perfect fluid of self-interacting dark matter (DM), the volume elements of which perform hydrodynamic flows. To the best of our knowledge, it is the first time in a cosmological framework that the energy of the cosmic fluid internal motions is also taken into account as a source of the universal gravitational field. As we demonstrate, this form of energy may compensate for the DE needed to compromise spatial flatness, while, depending on the particular type of thermodynamic processes occurring in the interior of the DM fluid (isothermal or polytropic), the Universe depicts itself as either decelerating or accelerating (respectively). In both cases, there is no disagreement between observations and the theoretical prediction of the distant supernovae (SNe) Type Ia distribution. In fact, the cosmological model with matter content in the form of a thermodynamically-involved DM fluid not only interprets the observational data associated with the recent history of Universe expansion, but also confronts successfully with every major cosmological issue (such as the age and the coincidence problems). In this way, depending on the type of thermodynamic processes in it, such a model may serve either for a conventional DE cosmology or for a viable alternative one.

Currently, a large amount of data implies that the matter constituents of the cosmological dark sector might be collisional. An attractive feature of such a possibility is that, it can reconcile dark matter (DM) and dark energy (DE) in terms of a single component, accommodated in the context of a polytropic-DM fluid. Accordingly, we explore the time evolution and the dynamical characteristics of a spatially-flat cosmological model, in which, in principle, there is no DE at all. Instead, in this model, the DM itself possesses some sort of fluid-like properties, i.e., the fundamental units of the Universe matter-energy content are the volume elements of a DM fluid, performing polytropic flows. In this case, the energy of this fluid's internal motions is also taken into account as a source of the universal gravitational field. This form of energy can compensate for the extra energy needed to compromise spatial flatness, namely, to justify that, today, the total-energy density parameter is exactly unity. The polytropic cosmological model, depends on only one free parameter, the corresponding exponent, \Gamma. What makes this model particularly interesting, is that, for \Gamma < 0.541, the (conventional) pressure becomes negative enough, so that the Universe accelerates its expansion at cosmological redshifts below a transition value. Several physical reasons impose further constraints on the value of \Gamma, which, eventually, is settled down to the range -0.089 < \Gamma < 0. Such a cosmological model does not suffer either from the age problem or from the coincidence problem. At the same time, this model reproduces to high accuracy the distance measurements performed with the aid of the supernovae Type Ia standard candles, and most naturally interprets, not only when, but also, why the Universe transits from deceleration to acceleration, thus arising as a mighty contestant for a DE model.

Motivated by results implying that the constituents of dark matter (DM) might be collisional, we consider a cosmological (toy-) model, in which the DM itself possesses some sort of thermodynamic properties. In this case, not only can the matter content of the Universe be treated as a classical gravitating fluid of positive pressure, but, together with all its other physical characteristics, the energy of this fluid's internal motions should be taken into account as a source of the universal gravitational field. This form of energy can compensate for the extra (dark) energy, needed to compromise spatial flatness, while the post-recombination Universe remains ever-decelerating. At the same time (i.e., in the context of the collisional-DM approach), the theoretical curve representing the distance modulus as a function of the cosmological redshift, {\mu}(z), fits the Hubble diagram of a multi-used sample of supernova Ia events quite accurately. A cosmological model filled with collisional DM could accommodate the majority of the currently-available observational data (including, also, those from baryon acoustic oscillations), without the need for either any dark energy (DE) or the cosmological constant. However, as we demonstrate, this is not the case for someone who, although living in a Universe filled with self-interacting DM, insists on adopting the traditional, collisionless-DM approach. From the point of view of this observer, the cosmologically-distant light-emitting sources seem to lie farther (i.e., they appear to be dimmer) than expected, while the Universe appears to be either accelerating or decelerating, depending on the value of the cosmological redshift. This picture, which, nowadays, represents the common perception in observational cosmology, acquires a more conventional interpretation within the context of the collisional-DM approach.

We use covariant methods to analyse the nonlinear evolution of self-gravitating, non-relativistic media. The formalism is first applied to imperfect fluids, aiming at the kinematic effects of viscosity, before extended to inhomogeneous magnetised environments. The nonlinear electrodynamic formulae are derived and successively applied to electrically resistive and to highly conductive fluids. By nature, the covariant equations isolate the magnetic effects on the kinematics and the dynamics of the medium, combining mathematical transparency and physical clarity. Employing the Newtonian analogue of the relativistic 1+3 covariant treatment, also facilitates the direct comparison with the earlier relativistic studies and helps to identify the differences in an unambiguous way. The purpose of this work is to set the framework and take a first step towards the detailed analytical study of complex nonlinear systems, like non-relativistic astrophysical plasmas and collapsing protogalactic clouds.

We use covariant techniques to examine the implications of the dynamical equivalence between geodesic motions and adiabatic hydrodynamic flows. Assuming that the metrics of a geodesically and a non-geodesically moving fluid are conformally related, we calculate and compare their mass densities. The density difference is then expressed in terms of the fundamental physical quantities of the fluid, such as its energy and isotropic pressure. Both the relativistic and the non-relativistic case are examined and their differences identified. Our analysis suggests that observational determinations of astrophysical masses based on purely Keplerian motions could underestimate the available amount of matter.

Using a highly accurate numerical code, we study the spin down of rotating relativistic stars, undergoing a quark deconfinement phase transition. Such phase transitions have been suggested to yield an observable signal in the braking index of spinning-down pulsars, which is based on a ``backbending'' behaviour of the moment of inertia. We focus on a particular equation of state that has been used before to study this behaviour, and find that for the population of normal pulsars the moment of inertia does not exhibit a backbending behaviour. In contrast, for supramassive millisecond pulsars a very strong backbending behaviour is found. Essentially, once a quark core appears in a spinning-down supramassive millisecond pulsar, the star spins up and continues to do so until it reaches the instability to collapse. This strong spin-up behaviour makes it easier to distinguish a phase transition in such pulsars: a negative first time-derivative of the rotational period suffices and one does not have to measure the braking index. In the spin-up era, the usually adopted spin-down power law fails to describe the evolution of the angular velocity. We adopt a general-relativistic spin-down power law and derive the equations that describe the angular velocity and braking index evolution in rapidly rotating pulsars.

The spin-down of a millisecond pulsar in a compact binary system leads to self-energy losses, that cause the binary's orbital period to increase (the effect being of first post-Newtonian (1PN) order). If the pulsar's period-derivative is not exceedingly small, this effect can become measurable in long-term high-precision timing measurements of the binary's orbital motion. We use rotating compressible spheroids to obtain an approximate, explicit expression for the orbital period derivative, due to spin-down, valid to 1PN order. This expression can be used to observationally constrain the pulsar's moment of inertia, and, combined with other observations, the high-density equation of state of compact stars. We apply our expression to representative models of millisecond pulsars in binary systems and demonstrate the importance of including higher-order rotational effects as well as the importance of the choice of equation of state, which both have been neglected, so far, in the literature. The computed increase in orbital period is larger, by as much as 32%, when the higher-order rotational effects are included, while it can change by more than a factor of two between different plausible equations of state.