This work is devoted to the thermodynamics of gravitational clustering, a collective phenomenon with a great relevance in the $N$-body cosmological problem. We study a classical self-gravitating gas of identical non-relativistic particles defined on the sphere $\mathbb{S}^{3}\subset \mathbb{R}^{4}$ by considering gravitational interaction that corresponds to this geometric space. The analysis is performed within microcanonical description of an isolated Hamiltonian system by combining continuum approximation and steepest descend method. According to numerical solution of resulting equations, the gravitational clustering can be associated with two microcanonical phase transitions. A first phase transition with a continuous character is associated with breakdown of $SO(4)$ symmetry of this model. The second one is the gravitational collapse, whose continuous or discontinuous character crucially depends on the regularization of short-range divergence of gravitation potential. We also derive the thermodynamic limit of this model system, the astrophysical counterpart of Gibbs-Duhem relation, the order parameters that characterize its phase transitions and the equation of state. Other interesting behavior is the existence of states with negative heat capacities, which appear when the effects of gravitation turn dominant for energies sufficiently low. Finally, we comment the relevance of some of these results in the study of astrophysical and cosmological situations. Special interest deserves the gravitational modification of the equation of state due to the local inhomogeneities of matter distribution. Although this feature is systematically neglected in studies about Universe expansion, the same one is able to mimic an effect that is attributed to the dark energy: a negative pressure.

General aspects about the thermodynamics of astrophysical systems are discussed, overall, those concerning to astrophysical systems in mutual interaction (or the called \emph{open astrophysical systems}). A special interest is devoted along the paper to clarify several misconceptions that are still common in the recent literature, such as the direct application to the astrophysical scenario of notions and theoretical frameworks that were originally conceived to deal with extensive systems of the everyday practice (large systems with short-range interactions).

We introduce a parametric family of models to characterize the properties of astrophysical systems in a quasi-stationary evolution under the incidence evaporation. We start from an one-particle distribution $f_{\gamma}\left(\mathbf{q},\mathbf{p}|\beta,\varepsilon_{s}\right)$ that considers an appropriate deformation of Maxwell-Boltzmann form with inverse temperature $\beta$, in particular, a power-law truncation at the scape energy $\varepsilon_{s}$ with exponent $\gamma>0$. This deformation is implemented using a generalized $\gamma$-exponential function obtained from the \emph{fractional integration} of ordinary exponential. As shown in this work, this proposal generalizes models of tidal stellar systems that predict particles distributions with \emph{isothermal cores and polytropic haloes}, e.g.: Michie-King models. We perform the analysis of thermodynamic features of these models and their associated distribution profiles. A nontrivial consequence of this study is that profiles with isothermal cores and polytropic haloes are only obtained for low energies whenever deformation parameter $\gamma<\gamma_{c}\simeq 2.13$.

The present effort addresses the question about the existence of a well-defined thermodynamic limit for the astrophysical systems with the following power law form: to tend the number of particles, N, the total energy, E, and the characteristic linear dimension of the system, L, to infinity, keeping constant E/N^{\Lambda_{E}} and L/N^{\Lambda_{L}}, being \Lambda_{E} and \Lambda_{L} certain scaling exponent constant. This study is carried out for a system constituted by a non-rotating fluid under the influence of its own Newtonian gravitational interaction. The analysis yields that a thermodynamic limit of the above form will only appear when the local pressure depends on the energy density of fluid as $p=\gamma\epsilon$, being $\gamma$ certain constant. Therefore, a thermodynamic limit with a power law form can be only satisfied by a reduced set of models, such as the selfgravitating gas of fermions and the Antonov isothermal model.

Astrophysical systems will never be in a real Thermodynamic equilibrium: they undergo an evaporation process due to the fact that the gravity is not able to confine the particles. Ordinarily, this difficulty is overcome by enclosing the system in a rigid container which avoids the evaporation. We proposed an energetic prescription which is able to confine the particles, leading in this way to an alternative version of the Antonov isothermal model which unifies the well-known isothermal and polytropic profiles. Besides of the main features of the isothermal sphere model: the existence of the gravitational collapse and the energetic region with a negative specific heat, this alternative model has the advantage that the system size naturally appears as a consequence of the particles evaporation.

We reconsider some general aspects about the mean field thermodynamical description of the astrophysical systems based on the microcanonical ensemble. Starting from these basis, we devote a special attention to the analysis of the scaling laws of the thermodynamical variables and potentials in the thermodynamic limit. Geometrical considerations motivate a way by means of which could be carried out a well-defined generalized canonical-like description for this kind of systems, even being nonextensive. This interesting possibility allows us to extend the applicability of the Standard Thermodynamic methods, even in the cases in which the system exhibits a negative specific heat. As example of application, we reconsider the classical Antonov problem of the isothermal spheres.

According to self-similarity hypothesis, the thermodynamic limit could be defined from the scaling laws for the system self-similarity by using the microcanonical ensemble. This analysis for selfgravitating systems yields the following thermodynamic limit: send N to infinity, keeping constant E/N^{(7/3)} and LN^{(1/3)}, in which is ensured the extensivity of the Boltzmann entropy S_{B}=lnW(E,N). It is shown how the consideration of this thermodynamic limit allows us to explain the origin of dynamical anomalies in numerical simulations of selfgravitating systems.