The equilibrium properties of classical self-gravitating systems in the grand canonical ensemble are studied by using the correspondence with an euclidean field theory with infrared and ultraviolet cutoffs. It is shown that the system developes a first order phase transition between a low and a high density regime. In addition, due to the long range of the gravitational potential, the system is close to criticality within each phase, with the exponents of mean field theory. The coexistence of a sharp first order transition and critical behavior can explain both the presence of voids in large regions of the universe as well as the self-similar density correlations in terms of self-gravity alone.

It is shown that the argument presented in astro-ph/0307318 against the proof that the diluted limit of a self-gravitating system is ill behaved [Nucl. Phys. B668, 409 (2003) (astro-ph/0303301)] is misleading.

It is shown that the diluted thermodynamical limit of a self-gravitating system proposed by de Vega and Sanchez suffers from the same problems as the usual thermodynamical limit and leads to divergent thermodynamical functions. This question is also discussed from the point of view of mean field theory.

The microcanonical statistical mechanics of a set of self-gravitating particles is analyzed in mean-field approach. In order to deal with an upper bounded entropy functional, a softened gravitational potential is used. The softening is achieved by truncating to N terms an expansion of the Newtonian potential in spherical Bessel functions. The order N is related to the softening at short distances. This regularization has the remarkable property that it allows for an exact solution of the mean field equation. It is found that for N not too large the absolute maximum of the entropy coincides to high accuracy with the solution of the Lane-Emden equation, which determines the mean field mass distribution for the Newtonian potential for energies larger than $E_c\approx -0.335 G M^2/R$. Below this energy a collapsing phase transition, with negative specific heat, takes place. The dependence of this result on the regularizing parameter N is discussed.

The effect of angular momentum conservation in microcanonical thermodynamics is considered. This is relevant in gravitating systems, where angular momentum is conserved and the collapsing nature of the forces makes the microcanonical ensemble the proper statistical description of the physical processes. The microcanonical distribution function with non-vanishing angular momentum is obtained as a function of the coordinates of the particles. As an example, a simple model of gravitating particles, introduced by Thirring long ago, is worked out. The phase diagram contains three phases: for low values of the angular momentum $L$ the system behaves as the original model, showing a complete collapse at low energies and an entropy with a convex intruder. For intermediate values of $L$ the collapse at low energies is not complete and the entropy still has a convex intruder. For large $L$ there is neither collapse nor anomalies in the thermodynamical quantities. A short discussion of the extension of these results to more realistic situations is exposed.