The self-gravitating gas in the Newtonian limit is studied in the presence of dark energy with a linear and constant equation of state. Entropy extremization associates to the isothermal Boltzmann distribution an effective density that includes `dark energy particles', which either strengthen or weaken mutual gravitational attraction, in case of quintessence or phantom dark energy, respectively, that satisfy a linear equation of state. Stability is studied for microcanonical (fixed energy) and canonical (fixed temperature) ensembles. Compared to the previously studied cosmological constant case, in the present work it is found that quintessence increases, while phantom dark energy decreases the instability domain under gravitational collapse. Thus, structures are more easily formed in a quintessence rather than in a phantom dominated Universe. Assuming that galaxy clusters are spherical, nearly isothermal and in hydrostatic equilibrium we find that dark energy with a linear and constant equation of state, for fixed radius, mass and temperature, steepens their total density profile. In case of a cosmological constant, this effect accounts for a 1.5% increase in the density contrast, that is the center to edge density ratio of the cluster. We also propose a method to constrain phantom dark energy.

We present here how the gravothermal or Antonov's instability, which was originally formulated in the microcanonical ensemble, is modified in the presence of a cosmological constant and in the canonical ensemble. In contrast to the microcanonical ensemble, there is a minimum, and not maximum, radius for which metastable states exist. In addition this critical radius is decreasing, and not increasing, with increasing cosmological constant. The minimum temperature for which metastable states exist is decreasing with increasing cosmological constant, while above some positive value of the cosmological constant, there appears a second critical temperature. For lower temperatures than the second critical temperature value, metastable states reappear, indicating a typical reentrant phase transition. The two critical temperatures merge when the cosmological density equals one half the mean density of the system.

Gravitational instabilities of isothermal spheres are studied in the presence of a positive or negative cosmological constant, in the Newtonian limit. In gravity, the statistical ensembles are not equivalent. We perform the analysis both in the microcanonical and the canonical ensembles, for which the corresponding instabilities are known as `gravothermal catastrophe' and `isothermal collapse', respectively. In the microcanonical ensemble, no equilibria can be found for radii larger than a critical value, which is increasing with increasing cosmological constant. In contrast, in the canonical ensemble, no equilibria can be found for radii smaller than a critical value, which is decreasing with increasing cosmological constant. For a positive cosmological constant, characteristic reentrant behaviour is observed.

We investigate the effect of a cosmological constant on the gravothermal catastrophe in the Newtonian limit. A negative cosmological constant acts as a thermodynamic `destabilizer'. The Antonov radius gets smaller and the instability occurs, not only for negative but also for positive energy values. A positive cosmological constant acts as a `stabilizer' of the system, which, in this case, exhibits a novel `reentrant behaviour'. In addition to the Antonov radius we find a second critical radius, where an `inverse Antonov transition' occurs; a series of local entropy maxima is restored.