Recent theoretical, numerical developments supported by observational evidence strongly suggest that many Globular Clusters host a black hole (BH) population in their centers. This is in contrast to previous long-standing belief that a BH subcluster would evaporate after undergoing core collapse and decoupling from the cluster. We propose that one mechanism which may add a stabilizing pressure to a BH population is the inhomogeneous Brownian motion generated by fluctuations of the stellar gravitational field. We argue that the diffusion equation for Brownian motion in inhomogeneous medium with spatially varying diffusion coefficient and temperature, discovered firstly by Van Kampen, applies to self-gravitating systems. Applying the stationary phase space probability distribution to a single BH immersed in a Plummer Globular Cluster we infer it may wanderer as far as $\sim 0.05,\,0.1,\,0.5{\rm pc}$ for mass $m_{\rm b} \sim 10^3, \,10^2,\,10{\rm M}_\odot$, respectively. We further find that the fluctuations of a fixed stellar mean gravitational field are sufficient to stabilize a BH population above the Spitzer instability threshold. We, nevertheless, identify an instability, whose onset depends on the Spitzer parameter $S = (M_{\rm b}/M_\star) (m_{\rm b}/m_\star)^{3/2} $ and parameter $B = \rho_{\rm b}(0) (4\pi r_c^3/M_b)(m_\star/m_{\rm b})^{3/2} $, where $\rho_{\rm b}(0)$ is the Brownian population central density. For a Plummer sphere the instability occurs at $(B,S) = (140,0.25)$. For $B > 140$ we get very cuspy BH subcluster profiles, unstable under the support of fluctuations alone. For $S > 0.25$ there do not exist stationary states of the BH population inhomogeneous diffusion equation.

Supernova theory suggests that black holes of a stellar origin cannot attain masses in the range of 50-135 solar masses in isolation. We argue here that this mass gap is filled in by black holes that grow by gas accretion in dense stellar clusters, such as protoglobular clusters. The accretion proceeds rapidly, during the first 10 megayears of the cluster life, before the remnant gas is depleted. We predict that binaries of black holes within the mass gap can be observed by LIGO.

I consider a self-gravitating, N-body system assuming that the N constituents follow regular orbits about the center of mass of the cluster, where a central massive object may be present. I calculate the average over a characteristic timescale of the full, N-body Hamiltonian including all kinetic and potential energy terms. The resulting effective system allows for the identification of the orbital planes with N rigid, disk-shaped tops, that can rotate about their fixed common centre and are subject to mutual gravitational torques. The time-averaging imposes boundaries on the canonical generalized momenta of the resulting canonical phase space. I investigate the statistical mechanics induced by the effective Hamiltonian on this bounded phase space and calculate the thermal equilibrium states. These are a result of the relaxation of spins' directions, identified with orbital planes' orientations, which is called vector resonant relaxation. I calculate the dependence of spins' angular velocity dispersion on temperature and calculate the velocity distribution functions. I argue that the range of validity of the gravitational phase transitions, identified in the special case of zero kinetic term by Roupas, Kocsis & Tremaine, is expanded to non-zero values of the ratio of masses between the cluster of N-bodies and the central massive object. The relevance with astrophysics is discussed focusing on stellar clusters. The same analysis performed on an unbounded phase space accounts for continuous rigid tops.

Thermal energy points toward a disordered, completely uniform state acting counter to gravity's tendency to generate order and structure through gravitational collapse. It is, therefore, expected to contribute to the stabilization of a self-gravitating, classical ideal gas over collapse. However, I identified an instability that always occurs at sufficiently high energies: the high-energy or relativistic gravothermal instability. I argue here that this instability presents an analogous core--halo structure as its Newtonian counterpart, the Antonov instability. The main difference is that in the former case the core is dominated by the gravitation of thermal energy and not rest mass energy. A relativistic generalization of Antonov's instability---the low-energy gravothermal instability---also occurs. The two turning points, which make themselves evident as a double spiral of the caloric curve, approach each other as relativistic effects become more intense and eventually merge in a single point. Thus, the high and low-energy cases may be realized as two aspects of a single phenomenon---the gravothermal instability---which involves a core--halo separation and an intrinsic heat flow. Finally, I argue that the core formed during a core-collapse supernova is subject to the relativistic gravothermal instability if it becomes sufficiently hot and compactified at the time of the bounce. In this case, it will continue to collapse towards the formation of a black hole.

We describe microcanonical phase transitions and instabilities of the ideal Fermi gas in general relativity at nonzero temperature confined in the interior of a spherical shell. The thermodynamic behaviour is governed by the compactness of rest mass, namely of the total rest mass over radius of the system. For a fixed value of rest mass compactness, we study the caloric curves as a function of the size of the spherical box. At low compactness values, low energies and for sufficiently big systems the system is subject to a gravothermal catastrophe, which cannot be halted by quantum degeneracy pressure, and the system collapses towards the formation of a black hole. For small systems, there appears no instability at low energies. For intermediate sizes, between two marginal values, gravothermal catastrophe is halted and a microcanonical phase transition occurs from a gaseous phase to a condensed phase with a nearly degenerate core. The system is subject to a relativistic instability at low energy, when the core gets sufficiently condensed above the Oppenheimer-Volkoff limit. For sufficiently high values of rest mass compactness the microcanonical phase transitions are suppressed. They are replaced either by an Antonov type gravothermal catastrophe for sufficiently big systems or by stable equilibria for small systems. At high energies the system is subject to the `relativistic gravothermal instability', identified by Roupas in [1], for all values of compactness and any size.

We show that binaries of stellar-mass black holes formed inside a young protoglobular cluster, can grow rapidly inside the cluster's core by accretion of the intracluster gas, before the gas may be depleted from the core. A black hole with mass of the order of eight solar masses can grow to values of the order of thirty five solar masses in accordance with recent gravitational waves signals observed by LIGO. Due to the black hole mass increase, a binary may also harden. The growth of binary black holes in a dense protoglobular cluster through mass accretion indicates a potentially important formation and hardening channel.

The gravitational instability, responsible for the formation of the structure of the Universe, occurs below energy thresholds and above spatial scales of a self-gravitating expanding region, when thermal energy can no longer counterbalance self-gravity. I argue that at sufficiently-large scales, dark energy may restore thermal stability. This stability re-entrance of an isothermal sphere defines a turnaround radius, which dictates the maximum allowed size of any structure generated by gravitational instability. On the opposite limit of high energies and small scales, I will show that an ideal, quantum or classical, self-gravitating gas is subject to a high-energy relativistic gravothermal instability. It occurs at sufficiently-high energy and small radii, when thermal energy cannot support its own gravitational attraction. Applications of the phenomenon include neutron stars and core-collapse supernovae. I also extend the original Oppenheimer--Volkov calculation of the maximum mass limit of ideal neutron cores to the non-zero temperature regime, relevant to the whole cooling stage from a hot proto-neutron star down to the final cold state.

We examine dense self-gravitating stellar systems dominated by a central potential, such as nuclear star clusters hosting a central supermassive black hole. Different dynamical properties of these systems evolve on vastly different timescales. In particular, the orbital-plane orientations are typically driven into internal thermodynamic equilibrium by vector resonant relaxation before the orbital eccentricities or semimajor axes relax. We show that the statistical mechanics of such systems exhibit a striking resemblance to liquid crystals, with analogous ordered-nematic and disordered-isotropic phases. The ordered phase consists of bodies orbiting in a disk in both directions, with the disk thickness depending on temperature, while the disordered phase corresponds to a nearly isotropic distribution of the orbit normals. We show that below a critical value of the total angular momentum, the system undergoes a first-order phase transition between the ordered and disordered phases. At the critical point the phase transition becomes second-order while for higher angular momenta there is a smooth crossover. We also find metastable equilibria containing two identical disks with mutual inclinations between $90^{\circ}$ and $180^\circ$.

The thermodynamic instabilities of the self-gravitating, classical ideal gas are studied in the case of static, spherically symmetric configurations in General Relativity taking into account the Tolman-Ehrenfest effect. One type of instabilities is found at low energies, where thermal energy becomes too weak to halt gravity and another at high energies, where gravitational attraction of thermal pressure overcomes its stabilizing effect. These turning points of stability are found to depend on the total rest mass $\mathcal{M}$ over the radius $R$. The low energy instability is the relativistic generalization of Antonov instability, which is recovered in the limit $G\mathcal{M} \ll R c^2$ and low temperatures, while in the same limit and high temperatures, the high energy instability recovers the instability of the radiation equation of state. In the temperature versus energy diagram of series of equilibria, the two types of gravothermal instabilities make themselves evident as a double spiral! The two energy limits correspond also to radius limits. So that, stable static configurations exist only in between two marginal radii for any fixed energy with negative thermal plus gravitational energy. Ultimate limits of rest mass, as well as total mass-energy, are reported. Applications to neutron cores are discussed.

In [1], the thermal equilibrium of static, spherically symmetric perfect fluids in General Relativity was studied. I would like to elaborate three points relevant to the results of [1]. The first point is only a clarification, summarized in theorem 1 below, of results that appear in [1]. The following two points correct the error in [1], stating that the condition for thermodynamic stability, found in [1], is referring to the microcanonical ensemble, while it was referring to the canonical one. In theorems 2 and 3, specific cases for which equivalence of dynamical and thermodynamic stability holds are specified.

Static thermal equilibrium of a quantum self-gravitating ideal gas in general relativity is studied at any temperature, taking into account the Tolman-Ehrenfest effect. Thermal contribution to the gravitational stability of static neutron cores is quantified. The curve of maximum mass with respect to temperature is reported. At low temperatures the Oppenheimer-Volkoff calculation is recovered, while at high temperatures the recently reported classical gas calculation is recovered. An ultimate upper mass limit $M = 2.43M_\odot$ of all maximum values is found to occur at Tolman temperature $ T = 1.27mc^2$ with radius $R = 15.2km$.

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.

Thermodynamical stability of fluid spheres is studied in the presence of a cosmological constant, both in the Newtonian limit, as well as in General Relativity. In all cases, an increase of the cosmological constant tends to stabilize the system, making asymptotically de Sitter space more thermodynamically stable than anti-de Sitter at the purely classical level. In addition, in the Newtonian case reentrant phase transitions are observed for a positive cosmological constant, due to its repelling property in this case. In General Relativity is studied the case of radiation, for which is found that the critical radius, at which an instability sets in, is always bigger than the black hole radius of the system and furthermore, at some value of the cosmological constant this critical radius hits at the cosmological horizon.

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.