No Arabic abstract
The supermassive black holes found at the centers of galaxies are often surrounded by dense star clusters. The ages of these clusters are generally longer than the resonant-relaxation time and shorter than the two-body relaxation time over a wide range of radii. We explore the thermodynamic equilibria of such clusters using a simple self-similar model. We find that the cluster exhibits a phase transition between a high-temperature spherical equilibrium and a low-temperature equilibrium in which the stars are on high-eccentricity orbits with nearly the same orientation. In the absence of relativistic precession, the spherical equilibrium is metastable below the critical temperature and the phase transition is first-order. When relativistic effects are important, the spherical equilibrium is linearly unstable below the critical temperature and the phase transition is continuous. A similar phase transition has recently been found in a model cluster composed of stars with a single semimajor axis. The presence of the same phenomenon in two quite different cluster models suggests that lopsided equilibria may form naturally in a wide variety of black-hole star clusters.
Supermassive black holes at the centres of galaxies are often surrounded by dense star clusters. For a wide range of cluster properties and orbital radii the resonant relaxation times in these clusters are much shorter than the Hubble time. Since resonant relaxation conserves semimajor axes, these clusters should be in the maximum-entropy state consistent with the given semimajor axis distribution. We determine these maximum-entropy equilibria in a simplified model in which all of the stars have the same semimajor axes. We find that the cluster exhibits a phase transition from a disordered, spherical, high-temperature equilibrium to an ordered low-temperature equilibrium in which the stellar orbits have a preferred orientation or line of apsides. Here `temperature is a measure of the non-Keplerian or self-gravitational energy of the cluster; in the spherical state, temperature is a function of the rms eccentricity of the stars. We explore a simple two-parameter model of black-hole star clusters -- the two parameters are semimajor axis and black-hole mass --- and find that clusters are susceptible to the lopsided phase transition over a range of ~100 in semimajor axis, mostly for black-hole masses less than $10^{7.5}$ solar masses.
The centers of most galaxies contain massive black holes surrounded by dense star clusters. The structure of these clusters determines the rate and properties of observable transient events, such as flares from tidally disrupted stars and gravitational-wave signals from stars spiraling into the black hole. Most estimates of these rates enforce spherical symmetry on the cluster. Here we show that, in the course of generic evolutionary processes, a star cluster surrounding a black hole can undergo a robust phase transition from a spherical thermal equilibrium to a lopsided equilibrium, in which most stars are on high-eccentricity orbits with aligned orientations. The rate of transient events is expected to be much higher in the ordered phase. Better models of cluster formation and evolution are needed to determine whether clusters should be found in the ordered or disordered phase.
Recent research has been constraining the retention fraction of black holes (BHs) in globular clusters by comparing the degree of mass segregation with $N$-body simulations. They are consistent with an upper limit of the retention fraction being $50,%$ or less. In this work, we focus on direct simulations of the dynamics of BHs in star clusters. We aim to constrain the effective distribution of natal kicks that BHs receive during supernova (SN) explosions and to estimate the BH retention fraction. We used the collisional $N$-body code nbody6 to measure the retention fraction of BHs for a given set of parameters, which are: the initial mass of a star cluster, the initial half-mass radius, and $sigma_mathrm{BH}$, which sets the effective Maxwellian BH velocity kick distribution. We compare these direct $N$-body models with our analytic estimates and newest observational constraints. The numerical simulations show that for the one-dimensional (1D) velocity kick dispersion $sigma_mathrm{BH} < 50,mathrm{km,s^{-1}}$, clusters with radii of 2 pc and that are initially more massive than $5 times 10^3,M_odot$ retain more than $20,%$ of BHs within their half-mass radii. Our simple analytic model yields a number of retained BHs that is in good agreement with the $N$-body models. Furthermore, the analytic estimates show that ultra-compact dwarf galaxies (UCDs) should have retained more than $80,%$ of their BHs for $sigma_mathrm{BH} leq 190,mathrm{km,s^{-1}}$. Although our models do not contain primordial binaries, in the most compact clusters with $10^3$ stars, we have found evidence of delayed SN explosions producing a surplus of BHs compared to the IMF due to dynamically formed binary stars. These cases do not occur in the more populous or expanded clusters.
Hierarchical triples are expected to be produced by the frequent binary-mediated interactions in the cores of globular clusters. In some of these triples, the tertiary companion can drive the inner binary to merger following large eccentricity oscillations, as a result of the eccentric Kozai-Lidov mechanism. In this paper, we study the dynamics and merger rates of black hole (BH) hierarchical triples, formed via binary--binary encounters in the CMC Cluster Catalog, a suite of cluster simulations with present-day properties representative of the Milky Ways globular clusters. We compare the properties of the mergers from triples to the other merger channels in dense star clusters, and show that triple systems do not produce significant differences in terms of mass and effective spin distribution. However, they represent an important pathway for forming eccentric mergers, which could be detected by LIGO--Virgo/KAGRA (LVK), and future missions such as LISA and DECIGO. We derive a conservative lower limit for the merger rate from this channel of $0.35$ Gpc$^{-3}$yr$^{-1}$ in the local Universe and up to $sim9%$ of these events may have a detectable eccentricity at LVK design sensitivity. Additionally, we find that triple systems could play an important role in retaining second-generation BHs, which can later merge again in the core of the host cluster.
The detection of gravitational waves emitted during a neutron star - black hole merger and the associated electromagnetic counterpart will provide a wealth of information about stellar evolution nuclear matter, and General Relativity. While the theoretical framework about neutron star - black hole binaries formed in isolation is well established, the picture is loosely constrained for those forming via dynamical interactions. Here, we use N-body simulations to show that mergers forming in globular and nuclear clusters could display distinctive marks compared to isolated mergers, namely larger masses, heavier black holes, and the tendency to have no associated electromagnetic counterpart. These features could represent a useful tool to interpreting forthcoming observations. In the Local Universe, gravitational waves emitted from dynamical mergers could be unravelled by detectors sensitive in the decihertz frequency band, while those occurring at the distance range of Andromeda and the Virgo Cluster could be accessible to lower-frequency detectors like LISA.