Given current cluster mass, orbital period, tidal stripping rate, and two-body relaxation time, compute OC's future mass curve and dissolution epoch.
Two processes deplete the cluster mass simultaneously. Tidal stripping removes a fraction f_strip of the remaining mass at each pericentre passage. After i passages, the tidal contribution gives M_tidal(i) = M_cl × (1 − f_strip)^i. Two-body evaporation is a continuous exponential: M_evap(t) = exp(−t / t_rh), where t_rh is the half-mass relaxation time. Combined:
M(t) = M_cl × (1 − f_strip)^(t/T_orb) × exp(−t / t_rh)
This treats the number of pericentre passages as a continuous variable (t/T_orb), which is an approximation; the discrete version produces a staircase that averages to the same curve.
T_diss(1%) is the smallest t where M(t)/M_cl < 0.01. T₁/₂ is where M(t)/M_cl < 0.5. These are solved numerically by scanning t from 0 to 50 Gyr in 10 Myr steps.
While the cluster exists, the IMBH undergoes energy-equipartition Brownian motion (Merritt 2001): v_BH = σ × √(⟨m_★⟩ / M_BH), where σ ~ 20 km/s and ⟨m_★⟩ ~ 0.5 M☉. After dissolution, the stellar restoring force vanishes and the IMBH wanders freely through the Galactic halo at a velocity comparable to its orbital speed in the Milky Way (~100–200 km/s).
Baumgardt & Makino 2003 (ApJ 600:204) — N-body dissolution timescales for tidally limited clusters. King 1962 (AJ 67:471) — tidal radius and King profile for globular clusters. Baumgardt & Hilker 2018 (MNRAS 478:1520) — photometric mass for Omega Centauri: 3.55 × 10⁶ M☉. Gieles et al. 2021 (MNRAS 507:4788) — OC survival timescale and tidal stripping. Merritt 2001 (AJ 121:2385) — Brownian motion of IMBH in stellar cusp.