No Arabic abstract
We study the dynamical evolution of globular clusters using our Henon-type Monte Carlo code for stellar dynamics including all relevant physics such as two-body relaxation, single and binary stellar evolution, Galactic tidal stripping, and strong interactions such as physical collisions and binary mediated scattering. We compute a large database of several hundred models starting from broad ranges of initial conditions guided by observations of young and massive star clusters. We show that these initial conditions very naturally lead to present day clusters with properties including the central density, core radius, half-light radius, half-mass relaxation time, and cluster mass, that match well with those of the old Galactic globular clusters. In particular, we can naturally reproduce the bimodal distribution in observed core radii separating the core-collapsed vs the non core-collapsed clusters. We see that the core-collapsed clusters are those that have reached or are about to reach the equilibrium binary burning phase. The non core-collapsed clusters are still undergoing gravo-thermal contraction.
Dynamical evolution drives globular clusters toward core collapse, which strongly shapes their internal properties. Diagnostics of core collapse have so far been based on photometry only, namely on the study of the concentration of the density profiles. Here we present a new method to robustly identify core-collapsed clusters based on the study of their stellar kinematics. We introduce the textit{kinematic concentration} parameter, $c_k$, the ratio between the global and local degree of energy equipartition reached by a cluster, and show through extensive direct $N$-body simulations that clusters approaching core collapse and in the post-core collapse phase are strictly characterized by $c_k>1$. The kinematic concentration provides a suitable diagnostic to identify core-collapsed clusters, independent from any other previous methods based on photometry. We also explore the effects of incomplete radial and stellar mass coverage on the calculation of $c_k$ and find that our method can be applied to state-of-art kinematic datasets.
Numerical and observational evidence suggests that massive white dwarfs dominate the innermost regions of core-collapsed globular clusters by both number and total mass. Using NGC 6397 as a test case, we constrain the features of white dwarf populations in core-collapsed clusters, both at present day and throughout their lifetimes. The dynamics of these white dwarf subsystems have a number of astrophysical implications. We demonstrate that the collapse of globular cluster cores is ultimately halted by the dynamical burning of white dwarf binaries. We predict core-collapsed clusters in the local universe yield a white dwarf merger rate of $mathcal{O}(10rm{),Gpc}^{-3},rm{yr}^{-1}$, roughly $0.1-1%$ of the observed Type Ia supernova rate. We show that prior to merger, inspiraling white dwarf binaries will be observable as gravitational wave sources at milli- and decihertz frequencies. Over $90%$ of these mergers have a total mass greater than the Chandrasekhar limit. If the merger/collision remnants are not destroyed completely in an explosive transient, we argue the remnants may be observed in core-collapsed clusters as either young neutron stars/pulsars/magnetars (in the event of accretion-induced collapse) or as young massive white dwarfs offset from the standard white dwarf cooling sequence. Finally, we show collisions between white dwarfs and main sequence stars, which may be detectable as bright transients, occur at a rate of $mathcal{O}(100rm{),Gpc}^{-3},rm{yr}^{-1}$ in the local universe. We find that these collisions lead to depletion of blue straggler stars and main sequence star binaries in the centers of core-collapsed clusters.
A self-interacting dark matter halo can experience gravothermal collapse, resulting in a central core with an ultrahigh density. It can further contract and collapse into a black hole, a mechanism proposed to explain the origin of supermassive black holes. We study dynamical instability of the core in general relativity. We use a truncated Maxwell-Boltzmann distribution to model the dark matter distribution and solve the Tolman-Oppenheimer-Volkoff equation. For given model parameters, we obtain a series of equilibrium configurations and examine their dynamical instability based on considerations of total energy, binding energy, fractional binding energy, and adiabatic index. The numerical results from our semi-analytical method are in good agreement with those from fully relativistic N-body simulations. We further show for the instability to occur in the classical regime, the boundary temperature of the core should be at least $10%$ of the mass of dark matter particles; for a $10^9~{rm M_odot}$ seed black hole, the particle mass needs to be larger than a few keV. These results can be used to constrain different collapse models, in particular, those with dissipative dark matter interactions.
We performed deep photometry of the central region of Galactic globular cluster M15 from archival Hubble Space Telescope data taken on the High Resolution Channel and Solar Blind Channel of the Advanced Camera for Surveys. Our data set consists of images in far-UV (FUV$_{140}$; F140LP), near-UV (NUV$_{220}$; F220W), and blue (B$_{435}$; F435W) filters. The addition of an optical filter complements previous UV work on M15 by providing an additional constraint on the UV-bright stellar populations. Using color-magnitude diagrams (CMDs) we identified several populations that arise from non-canonical evolution including candidate blue stragglers, extreme horizontal branch stars, blue hook stars (BHks), cataclysmic variables (CVs), and helium-core white dwarfs (He WDs). Due to preliminary identification of several He WD and BHk candidates, we add M15 as a cluster containing a He WD sequence and suggest it be included among clusters with a BHk population. We also investigated a subset of CV candidates that appear in the gap between the main sequence (MS) and WDs in FUV$_{140}-$NUV$_{220}$ but lie securely on the MS in NUV$_{220}-$B$_{435}$. These stars may represent a magnetic CV or detached WD-MS binary population. Additionally, we analyze our candidate He WDs using model cooling sequences to estimate their masses and ages and investigate the plausibility of thin vs. thick hydrogen envelopes. Finally, we identify a class of UV-bright stars that lie between the horizontal branch and WD cooling sequences, a location not usually populated on cluster CMDs. We conclude these stars may be young, low-mass He WDs.
Globular clusters (GCs) in the Milky Way exhibit a well-observed bimodal distribution in core radii separating the so-called core-collapsed and non-core-collapsed clusters. Here, we use our Henon-type Monte Carlo code, CMC, to explore initial cluster parameters that map into this bimodality. Remarkably, we find that by varying the initial size of clusters (specified in our initial conditions in terms of the initial virial radius, $r_v$) within a relatively narrow range consistent with the measured radii of young star clusters in the local universe ($r_v approx 0.5-5$ pc), our models reproduce the variety of present-day cluster properties. Furthermore, we show that stellar-mass black holes (BHs) play an intimate role in this mapping from initial conditions to the present-day structural features of GCs. We identify best-fit models for three GCs with known observed BH candidates, NGC 3201, M22, and M10, and show that these clusters harbor populations of $sim 50-100$ stellar-mass BHs at present. As an alternative case, we also compare our models to the core-collapsed cluster NGC 6752 and show that this cluster likely contains few BHs at present. Additionally, we explore the formation of BH binaries in GCs and demonstrate that these systems form naturally in our models in both detached and mass-transferring configurations with a variety of companion stellar types, including low-mass main sequence stars, white dwarfs, and sub-subgiants.