[Paper Review] On the dynamical evolution of globular clusters

M. Hénon|arXiv (Cornell University)|Mar 17, 2011

Mathematics and Applications12 citations

TL;DR

This paper presents a theoretical model of globular cluster evolution using a self-consistent, time-dependent, spherically symmetric system of equations under the assumption of velocity isotropy and equal stellar masses. The key result is the discovery of a stable, unique 'homology model' where mass loss occurs via escape of stars (30% of negative energy), while the remaining 70% accumulates in the core as multiple stars, leading to core collapse and finite central density growth over time, in agreement with observational data on M3 and 47 Tuc.

ABSTRACT

This paper is an English translation of Michel Hénon's thesis, "Sur l'évolution dynamique des amas globulaires" originally published in French in the Annales d'Astrophysique, Vol. 24, p.369 (1961).

Motivation & Objective

• To develop a unified theoretical framework for the dynamical evolution of globular clusters, treating structure and evolution simultaneously.
• To model the long-term evolution of stellar systems under self-gravity, assuming spherical symmetry, velocity isotropy, and equal masses.
• To determine whether a stable, self-similar solution (homology model) exists that describes the final state of globular clusters.
• To compare theoretical predictions with observational data, including projected densities, mass functions, and escape rates.

Proposed method

• Derives a system of fundamental equations (2.25) governing the evolution of the distribution function in phase space under spherical symmetry and isotropic velocity dispersion.
• Applies a time-dependent homology transformation to reduce the model to a canonical form, seeking a self-similar solution invariant under scaling.
• Imposes physical constraints: mass conservation, energy flux toward the center, and finite mass and radius, leading to a unique solution—the homology model.
• Performs numerical integration using an electronic calculator to solve the system of equations, yielding the homology model with finite central density and mass loss.
• Introduces mass-segregated stars into the model to study their spatial distribution and escape rates, showing mass-dependent concentration and reduced escape for heavier stars.
• Compares theoretical predictions with observational data: projected density profiles (M3, 47 Tuc), mass-luminosity relations, and cluster mass functions to validate the model.

Experimental results

Research questions

• RQ1 Can a self-consistent, time-evolving model of globular clusters be constructed that satisfies both structural and dynamical constraints?
• RQ2 Does a unique, stable, self-similar solution (homology model) exist that describes the final state of cluster evolution?
• RQ3 What is the fate of energy in the cluster—how is it distributed between escaping stars and the core?
• RQ4 How do stars of different masses behave in the evolving cluster, and what is their escape rate?
• RQ5 To what extent do the theoretical predictions match observational data on cluster structure and mass loss?

Key findings

• The homology model is the unique solution satisfying all physical constraints: finite mass, finite radius, and self-similarity under time-dependent scaling.
• The cluster's mass decreases linearly over time at a rate of $2.3 \times 10^{-6}\, \mathrm{M}_{\odot}/\mathrm{yr}$, consistent with observed mass functions.
• One-third of the cluster's negative energy is carried off by escaping stars, while two-thirds accumulate in the core, leading to core collapse.
• Stars with mass greater than 1.5 times the average mass are nearly all concentrated in the center, with zero escape rate.
• Theoretical projected density profiles for M3 and 47 Tuc agree well with observations, validating the model's predictive power.
• The central density grows to infinity in finite time, and the homology model is stable against small perturbations, suggesting it is the final evolutionary state.