Cosmological reconstruction of Little Rip model in $f(R, T)$ gravity is investigated, where $R$ is the curvature scalar and $T$ the trace of the energy momentum tensor. The model perfectly reproduces the present stage of the universe, characterized by the $\Lambda CDM$ model, without singularity at future finite-time (without the Big Rip). The input parameters are determined according to Supernovae Cosmology data and perfectly fit with the WMAP around the Little Rip. Moreover, the thermodynamics is considered in this Little Rip cosmology and it is illustrated that the second law of thermodynamics is always satisfied around the Little Rip universe for the temperature inside the horizon being the same as that of the apparent horizon. Moreover, we show the existence of a stable fixed point in the Little Rip universe which confirms that this is actually a late-time attractor in the phantom-dominated universe. The linear perturbation analysis is performed around the critical points, showing that the Little Rip model obtained is stable.

We obtain a new solution of the Einstein-anti-Maxwell theory with cosmological constant, called anti-Reissner-Nordstrom-(A)de Sitter (anti-RN-(A)dS) solution. The basic properties of this solution is reviewed. Its thermodynamics is consistently established, with the extreme cases and phase transitions, where the analysis is performed through two methods, the usual one and that of Geometrothermodynamics . The Geometrothermodynamics analysis does not provide a result in agreement with the usual method or by the specific heat. We establish local and global thermodynamic stabilities of anti-RN-AdS solution through the specific heat and the canonical and grand-canonical ensembles.

The dark energy component of the cosmic budget is represented by a self-interacting scalar field. The violation of the null energy condition is allowed. Hence, such component can also represent a phantom fluid. The model is tested using supernova type Ia and matter power spectrum data. The supernova test leads to preferred values for configurations representing the phantom fluid. The matter power spectrum constraints for the dark energy equation of state parameter are highly degenerated. In both cases, values for the equation of state parameter corresponding to the phantom fluid are highly admitted if no particular prior is used.

Scalar perturbations can grow during a phantomic cosmological phase as the big rip is approached, in spite of the high accelerated expansion regime, if the equation of state is such that $\frac{p}{\rho} = \alpha < - {5/3}$. It is shown that such result is independent of the spatial curvature. The perturbed equations are exactly solved for any value of the curvature parameter $k$ and of the equation of state parameter $\alpha$. Growing modes are found asymptotically under the condition $\alpha < - {5/3}$. Since the Hubble radius decreases in a phantom universe, such result indicates that a phantom scenario may not survive longtime due to gravitational instability.