No Arabic abstract
The classical theory of cluster relaxation is unsatisfactory because it involves the Coulomb logarithm. The Balescu-Lenard (BL) equation provides a rigorous alternative that has no ill-defined parameter. Moreover, the BL equation, unlike classical theory, includes the clusters self-gravity. A heuristic argument is given that indicates that relaxation does not occur predominantly through two-particle scattering and is enhanced by self-gravity. The BL equation is adapted to a spherical system and used to estimate the flux through the action space of isochrone clusters with different velocity anisotropies. A range of fairly different secular behaviours is found depending on the fraction of radial orbits. Classical theory is also used to compute the corresponding classical fluxes. The BL and classical fluxes are very different because (a) the classical theory materially under-estimates the impact of large-scale collectively amplified fluctuations and (b) only the leading terms in an infinite sum for the BL flux are computed. A complete theory of cluster relaxation likely requires that the sum in the BL equation be decomposed into a sum over a finite number of small wavenumbers complemented by an integral over large wavenumbers analogous to classical theory.
Globular clusters contain a finite number of stars. As a result, they inevitably undergo secular evolution (`relaxation) causing their mean distribution function (DF) to evolve on long timescales. On one hand, this long-term evolution may be interpreted as driven by the accumulation of local deflections along each stars mean field trajectory -- so-called `non-resonant relaxation. On the other hand, it can be thought of as driven by non-local, collectively dressed and resonant couplings between stellar orbits, a process termed `resonant relaxation. In this paper we consider a model globular cluster represented by a spherical, isotropic isochrone DF, and compare in detail the predictions of both resonant and non-resonant relaxation theories against tailored direct $N$-body simulations. In the space of orbital actions (namely the radial action and total angular momentum), we find that both resonant and non-resonant theories predict the correct morphology for the secular evolution of the clusters DF, although non-resonant theory over-estimates the amplitude of the relaxation rate by a factor ${sim 2}$. We conclude that the secular relaxation of hot isotropic spherical clusters is not dominated by collectively amplified large-scale potential fluctuations, despite the existence of a strong ${ell = 1}$ damped mode. Instead, collective amplification affects relaxation only marginally even on the largest scales. The predicted contributions to relaxation from smaller scale fluctuations are essentially the same from resonant and non-resonant theories.
We perform a systematic combined photometric and kinematic analysis of a sample of globular clusters under different relaxation conditions, based on their core relaxation time (as listed in available catalogs), by means of two well-known families of spherical stellar dynamical models. Systems characterized by shorter relaxation time scales are expected to be better described by isotropic King models, while less relaxed systems might be interpreted by means of non-truncated, radially-biased anisotropic f^( u) models, originally designed to represent stellar systems produced by a violent relaxation formation process and applied here for the first time to the study of globular clusters. The comparison between dynamical models and observations is performed by fitting simultaneously surface brightness and velocity dispersion profiles. For each globular cluster, the best-fit model in each family is identified, along with a full error analysis on the relevant parameters. Detailed structural properties and mass-to-light ratios are also explicitly derived. We find that King models usually offer a good representation of the observed photometric profiles, but often lead to less satisfactory fits to the kinematic profiles, independently of the relaxation condition of the systems. For some less relaxed clusters, f^( u) models provide a good description of both observed profiles. Some derived structural characteristics, such as the total mass or the half-mass radius, turn out to be significantly model-dependent. The analysis confirms that, to answer some important dynamical questions that bear on the formation and evolution of globular clusters, it would be highly desirable to acquire larger numbers of accurate kinematic data-points, well distributed over the cluster field.
Internal rotation is considered to play a major role in the dynamics of some globular clusters. However, in only few cases it has been studied by quantitative application of realistic and physically justified global models. Here we present a dynamical analysis of the photometry and three-dimensional kinematics of omega Cen, 47 Tuc, and M15, by means of a recently introduced family of self-consistent axisymmetric rotating models. The three clusters, characterized by different relaxation conditions, show evidence of differential rotation and deviations from sphericity. The combination of line-of-sight velocities and proper motions allows us to determine their internal dynamics, predict their morphology, and estimate their dynamical distance. The well-relaxed cluster 47 Tuc is very well interpreted by our model; internal rotation is found to explain the observed morphology. For M15, we provide a global model in good agreement with the data, including the central behavior of the rotation profile and the shape of the ellipticity profile. For the partially relaxed cluster omega Cen, the selected model reproduces the complex three-dimensional kinematics; in particular the observed anisotropy profile, characterized by a transition from isotropy, to weakly-radial anisotropy, and then to tangential anisotropy in the outer parts. The discrepancy found for the steep central gradient in the observed line-of-sight velocity dispersion profile and for the ellipticity profile is ascribed to the condition of only partial relaxation of this cluster and the interplay between rotation and radial anisotropy.
Blue hook (BHk) stars are a rare class of horizontal branch stars that so far have been found in only very few Galactic globular clusters (GCs). The dominant mechanism for producing these objects is currently still unclear. In order to test if the presence of BHk populations in a given GC is linked to specific physical or structural cluster properties, we have constructed a parent sample of GCs for which existing data is sufficient to establish the presence or absence of BHk populations with confidence. We then compare the properties of those clusters in our parent sample that do contain a BHk population to those that do not. We find that there is only one compelling difference between BHk and non-BHk clusters: all known BHk clusters are unusually massive. However, we also find that the BHk clusters are consistent with being uniformly distributed within the cumulative mass distribution of the parent sample. Thus, while it is attractive to suggest there is is a lower mass cut-off for clusters capable of forming BHk stars, the data do not require this. Instead, the apparent preference for massive clusters could still be a purely statistical effect: intrinsically rare objects can only be found by searching a sufficiently large number of stars.
It is widely believed that globular clusters evolve over many two-body relaxation times toward a state of energy equipartition, so that velocity dispersion scales with stellar mass as sigma ~ m^{-eta} with eta = 0.5. We show that this is incorrect, using direct N-body simulations with a variety of realistic IMFs and initial conditions. No simulated system ever reaches a state close to equipartition. Near the center, the luminous main-sequence stars reach a maximum eta_{max} ~ 0.15 pm 0.03. At large times, all radial bins convergence on an asymptotic value eta_{infty} ~ 0.08 pm 0.02. The development of this partial equipartition is strikingly similar across our simulations, despite the range of initial conditions employed. Compact remnants tend to have higher eta than main-sequence stars (but still eta < 0.5), due to their steeper (evolved) mass function. The presence of an intermediate-mass black hole (IMBH) decreases eta, consistent with our previous findings of a quenching of mass segregation under these conditions. All these results can be understood as a consequence of the Spitzer instability for two-component systems, extended by Vishniac to a continuous mass spectrum. Mass segregation (the tendency of heavier stars to sink toward the core) has often been studied observationally, but energy equipartition has not. Due to the advent of high-quality proper motion datasets from the Hubble Space Telescope, it is now possible to measure eta. Detailed data-model comparisons open up a new observational window on globular cluster dynamics and evolution. Comparison of our simulations to Omega Cen observations yields good agreement, confirming that globular clusters are not generally in energy equipartition. Modeling techniques that assume equipartition by construction (e.g., multi-mass Michie-King models) are approximate at best.