No Arabic abstract
We present the results of realistic N-body modelling of massive star clusters in the Magellanic Clouds, aimed at investigating a dynamical origin for the radius-age trend observed in these systems. We find that stellar-mass black holes, formed in the supernova explosions of the most massive cluster stars, can constitute a dynamically important population. If a significant number of black holes are retained (here we assume complete retention), these objects rapidly form a dense core where interactions are common, resulting in the scattering of black holes into the cluster halo, and the ejection of black holes from the cluster. These two processes heat the stellar component, resulting in prolonged core expansion of a magnitude matching the observations. Significant core evolution is also observed in Magellanic Cloud clusters at early times. We find that this does not result from the action of black holes, but can be reproduced by the effects of mass-loss due to rapid stellar evolution in a primordially mass segregated cluster.
Our current understanding of the stellar initial mass function and massive star evolution suggests that young globular clusters may have formed hundreds to thousands of stellar-mass black holes, the remnants of stars with initial masses from $sim 20 - 100, M_odot$. Birth kicks from supernova explosions may eject some black holes from their birth clusters, but most should be retained. Using a Monte Carlo method we investigate the long-term dynamical evolution of globular clusters containing large numbers of stellar black holes. We describe numerical results for 42 models, covering a range of realistic initial conditions, including up to $1.6times10^6$ stars. In almost all models we find that significant numbers of black holes (up to $sim10^3$) are retained all the way to the present. This is in contrast to previous theoretical expectations that most black holes should be ejected dynamically within a few Gyr. The main reason for this difference is that core collapse driven by black holes (through the Spitzer mass segregation instability) is easily reverted through three-body processes, and involves only a small number of the most massive black holes, while lower-mass black holes remain well-mixed with ordinary stars far from the central cusp. Thus the rapid segregation of stellar black holes does not lead to a long-term physical separation of most black holes into a dynamically decoupled inner core, as often assumed previously. Combined with the recent detections of several black hole X-ray binary candidates in Galactic globular clusters, our results suggest that stellar black holes could still be present in large numbers in many globular clusters today, and that they may play a significant role in shaping the long-term dynamical evolution and the present-day dynamical structure of many clusters.
Globular clusters should be born with significant numbers of stellar-mass black holes (BHs). It has been thought for two decades that very few of these BHs could be retained through the cluster lifetime. With masses ~10 MSun, BHs are ~20 times more massive than an average cluster star. They segregate into the cluster core, where they may eventually decouple from the remainder of the cluster. The small-N core then evaporates on a short timescale. This is the so-called Spitzer instability. Here we present the results of a full dynamical simulation of a globular cluster containing many stellar-mass BHs with a realistic mass spectrum. Our Monte Carlo simulation code includes detailed treatments of all relevant stellar evolution and dynamical processes. Our main finding is that old globular clusters could still contain many BHs at present. In our simulation, we find no evidence for the Spitzer instability. Instead, most of the BHs remain well-mixed with the rest of the cluster, with only the innermost few tens of BHs segregating significantly. Over the 12 Gyr evolution, fewer than half of the BHs are dynamically ejected through strong binary interactions in the cluster core. The presence of BHs leads to long-term heating of the cluster, ultimately producing a core radius on the high end of the distribution for Milky Way globular clusters (and those of other galaxies). A crude extrapolation from our model suggests that the BH--BH merger rate from globular clusters could be comparable to the rate in the field.
Black holes formed in dense star clusters, where dynamical interactions are frequent, may have fundamentally different properties than those formed through isolated stellar evolution. Theoretical models for single star evolution predict a gap in the black hole mass spectrum from roughly $40-120,M_{odot}$ caused by (pulsational) pair-instability supernovae. Motivated by the recent LIGO/Virgo event GW190521, we investigate whether black holes with masses within or in excess of this upper-mass gap can be formed dynamically in young star clusters through strong interactions of massive stars in binaries. We perform a set of $N$-body simulations using the CMC cluster-dynamics code to study the effects of the high-mass binary fraction on the formation and collision histories of the most massive stars and their remnants. We find that typical young star clusters with low metallicities and high binary fractions in massive stars can form several black holes in the upper-mass gap and often form at least one intermediate-mass black hole. These results provide strong evidence that dynamical interactions in young star clusters naturally lead to the formation of more massive black hole remnants.
Massive black holes (MBHs) are nowadays recognized as integral parts of galaxy evolution. Both the approximate proportionality between MBH and galaxy mass, and the expected importance of feedback from active MBHs in regulating star formation in their host galaxies point to a strong interplay between MBHs and galaxies. MBHs must form in the first galaxies and be fed by gas in these galaxies, with continuous or intermittent inflows that, at times, can be larger than the Eddington rate. Feedback from supernovae and from the MBHs themselves modulates the growth of the first MBHs. While current observational data only probe the most massive and luminous MBHs, the tip of the iceberg, we will soon be able to test theoretical models of MBH evolution on more normal MBHs: the MBHs that are indeed relevant in building the population that we observe in local galaxies, including our own Milky Way.
Intermediate-mass black holes (IMBHs) could form via runaway merging of massive stars in a young massive star cluster (YMC). We combine a suite of numerical simulations of YMC formation with a semi-analytic model for dynamical friction and merging of massive stars and evolution of a central quasi-star, to predict how final quasi-star and relic IMBH masses scale with cluster properties (and compare with observations). The simulations argue that inner YMC density profiles at formation are steep (approaching isothermal), producing some efficient merging even in clusters with relatively low effective densities, unlike models which assume flat central profiles resembling those of globular clusters (GCs) {em after} central relaxation. Our results can be approximated by simple analytic scalings, with $M_{rm IMBH} propto v_{rm cl}^{3/2}$ where $v_{rm cl}^{2} = G,M_{rm cl}/r_{rm h}$ is the circular velocity in terms of initial cluster mass $M_{rm cl}$ and half-mass radius $r_{rm h}$. While this suggests IMBH formation is {em possible} even in typical clusters, we show that predicted IMBH masses for these systems are small, $sim 100-1000,M_{odot}$ or $sim 0.0003,M_{rm cl}$, below even the most conservative observational upper limits in all known cases. The IMBH mass could reach $gtrsim 10^{4},M_{odot}$ in the centers nuclear star clusters, ultra-compact dwarfs, or compact ellipticals, but in all these cases the prediction remains far below the present observed supermassive BH masses in these systems.