---
title: Einstein Cluster Formalism
url: https://www.emergentmind.com/topics/einstein-cluster-formalism
type: topic
---

# Einstein Cluster Formalism

The Einstein cluster formalism is a relativistic description of a collisionless ensemble of identical particles moving on circular timelike geodesics in a static, spherically symmetric spacetime. In Maeda et al., “Einstein Cluster as Central Spiky Distribution of Galactic Dark Matter” [2410.04175], the formalism is used to construct a fully relativistic, spherically symmetric, spiky structure of matter distribution near a supermassive black hole, including three simple toy models and a more realistic Hernquist-type model. In this framework, the matter distribution is not represented by an isotropic fluid: the radial pressure vanishes identically, while the angular directions carry pressure. A central result is that the self-gravity of the cluster modifies strong-field geodesic stability and shifts the innermost stable circular orbit (ISCO) of the combined black-hole-plus-cluster system inward from the Schwarzschild value.

## 1. Geometric setup and definition of the cluster

The spacetime is assumed to be static and spherically symmetric, with line element
$$
ds^2=-e^{2\Phi(r)}dt^2+\frac{dr^2}{1-2m(r)/r}+r^2(d\theta^2+\sin^2\theta\,d\phi^2).
$$
Here $t\in(-\infty,+\infty)$ is a timelike coordinate, $(r,\theta,\phi)$ are standard spherical coordinates, $m(r)$ is half the Misner–Sharp mass inside radius $r$, and the metric coefficient satisfies
$$
g_{rr}=\left[1-\frac{2m(r)}{r}\right]^{-1}.
$$
The function $\Phi(r)$ is a redshift potential, with asymptotic flatness imposed through $\Phi(r\to\infty)\to 0$ [2410.04175].

The matter source is an Einstein cluster: a collisionless ensemble of identical particles of rest mass $\mu_0$, each moving on a circular timelike geodesic in the equatorial plane. The configuration is therefore intrinsically anisotropic. Because all particles are on circular orbits, the radial pressure vanishes identically, $P_r=0$, while the tangential directions support nonzero pressure. This gives the formalism its characteristic structure: spherical symmetry is retained at the level of the coarse-grained stress tensor, but the microscopic matter content is orbital rather than hydrostatic.

A common misconception is to identify the Einstein cluster with an ordinary perfect fluid. The formalism explicitly excludes that identification: the effective matter source has $p_r(r)=0$ by construction and equal tangential pressures in the $\theta$ and $\phi$ directions, so its stress tensor is anisotropic rather than isotropic.

## 2. Stress–energy tensor and effective matter variables

For $N$ point particles, the stress–energy tensor is written as
$$
T^\mu{}_\nu=\mu_0\sum_{I=1}^N\int d\tau_I\,u_I^\mu u_{I\nu}\,
\frac{\delta^4[x-z_I(\tau_I)]}{\sqrt{-g}},
$$
where $z_I(\tau)$ is the worldline of particle $I$ and $u_I^\mu\equiv dz_I^\mu/d\tau_I$ is its four-velocity. After smoothing over many particles, one obtains the anisotropic perfect-fluid form
$$
T^\mu{}_\nu=\mathrm{diag}[-\rho(r),\,0,\,p_t(r),\,p_t(r)]
$$
in the $(t,r,\theta,\phi)$ frame [2410.04175].

The energy density $\rho(r)$ is measured by a static observer. The radial pressure vanishes, and the tangential pressure $p_t(r)$ is equal in the $\theta$ and $\phi$ directions. If $n(r)$ denotes the number density of particles in coordinate volume, and if each particle on a circular orbit at radius $r$ has conserved energy per unit mass $E$ and angular momentum magnitude $L$, then
$$
\rho(r)=u^0u_0\,n(r)=\frac{E^2}{e^{2\Phi}}\,n(r)=\left(1+\frac{L^2}{r^2}\right)n(r),
$$
while the average tangential pressure is
$$
p_t(r)=\langle u^\phi u_\phi\rangle\,n(r)=\frac{L^2}{2r^2}n(r).
$$
The factor $1/2$ arises from averaging over the two independent angular directions $\theta$ and $\phi$.

These relations make clear that the effective stress tensor is a coarse-grained encoding of orbital motion. The pressure is not thermodynamic in origin; it is the macroscopic manifestation of angular momentum support in a collisionless distribution.

## 3. Einstein equations, anisotropy condition, and angular-momentum representation

For the diagonal stress tensor above, Einstein’s equations reduce to two ordinary differential equations in $r$ together with an algebraic anisotropy condition. The first is the “Hamiltonian” equation,
$$
m'(r)=4\pi r^2\rho(r).
$$
The second is a Tolman–Oppenheimer–Volkoff-type equation,
$$
\Phi'(r)=\frac{m(r)+4\pi r^3 p_t(r)}{r[r-2m(r)]}.
$$
Conservation, $\nabla_\mu T^\mu{}_r=0$, or equivalently the $\theta$–$\theta$ Einstein equation, yields the anisotropy condition enforcing $P_r=0$:
$$
p_t(r)=\frac{1}{2}\frac{r\,m'(r)}{r-2m(r)}
      =\frac{m(r)}{2[r-2m(r)]}\rho(r).
$$
Once $\rho(r)$ is specified, $p_t(r)$ is therefore fixed algebraically by $m(r)$, and $\Phi(r)$ follows by quadrature [2410.04175].

The same formalism can be recast in terms of the angular-momentum distribution. Circular timelike geodesics in the metric above satisfy
$$
E^2=e^{2\Phi}\left(1+\frac{L^2}{r^2}\right), \qquad
L^2=\frac{m(r)\,r^2}{r-3m(r)}.
$$
One may therefore choose $n(r)$ or directly choose $L(r)$ and solve for the remaining variables through $m'(r)=4\pi r^2\rho$. Then $p_t(r)=L^2 n(r)/(2r^2)$ follows automatically.

This dual description is one of the formalism’s defining features. One may regard the cluster either as a density profile that determines the geometry or as an orbital distribution whose angular-momentum content generates the effective anisotropic matter sector.

## 4. Analytic model families

Maeda et al. present three elementary analytic models and a more realistic Hernquist-type construction in which $m(r)$, and hence $\rho$, $p_t$, and $\Phi$, can be written in closed form [2410.04175].

| Model | Mass function | Characteristic feature |
|---|---|---|
| Model I | $m(r)=M_{\rm BH}+\alpha(r-r_I)$ | “isothermal-like” |
| Model II | $m(r)=\dfrac{r}{3[1+(r/r_c)^2]}$ | “everywhere marginally stable” |
| Model III | $m(r)=m_p r^p$ | power-law cluster |
| Hernquist-type | $m(r)=\dfrac{m_0+m_1 r+M r^2}{(r+r_*)^2}$ | characteristic scale $r_*\gg M_{\rm BH}$ |

In Model I, defined for $r_I\le r\le r_O$, one takes
$$
m(r)=M_{\rm BH}+\alpha(r-r_I), \qquad 0<\alpha<\frac13.
$$
Inside $r_I$ the geometry is vacuum Schwarzschild with $m(r)=M_{\rm BH}$, and outside $r_O$ one glues to Schwarzschild of total mass $M\equiv m(r_O)$. The energy density is
$$
\rho(r)=\frac{\alpha}{4\pi r^2},
$$
and the metric function $f\equiv e^{2\Phi}$ obeys
$$
\frac{d\ln f}{dr}=\frac{2m}{r(r-2m)},
$$
with solution
$$
f(r)=c_0\,r^{-1}[r-r_*]^{1/(1-2\alpha)},
$$
where
$$
r_*=\frac{2(1-\alpha r_I)}{1-2\alpha},
$$
and $c_0$ is fixed by continuity at $r_I$ and $r_O$.

In Model II, the mass profile is chosen so that the ISCO condition
$$
r^2m'(r)+r\,m(r)-6m(r)^2=0
$$
holds identically throughout the cluster. This gives
$$
m(r)=\frac{r}{3[1+(r/r_c)^2]},
$$
with
$$
\rho(r)=\frac{1-(r/r_c)^2}{12\pi r^2[1+(r/r_c)^2]^2},
\qquad
f(r)=c_0\,\frac{r^2}{1+3(r/r_c)^2}.
$$
Vacuum is glued at $r_I$ where $m(r_I)=M_{\rm BH}$, and at $r_O\le r_c$ to obtain a finite-mass configuration.

In Model III, one chooses
$$
m(r)=m_p r^p, \qquad 0<p<1,
$$
for $r_I\le r\le r_O$, and imposes $m(r_I)=M_{\rm BH}$ together with the ISCO condition at $r_I$. This yields
$$
p=\frac{6-r_I}{r_I}, \qquad m_p=r_I^{-p}.
$$
The density becomes
$$
\rho(r)=\frac{p\,m_p\,r^{p-3}}{4\pi},
$$
and
$$
f(r)=c_0\left[1-\frac{2m_p}{r^{1-p}}\right]^{1/(1-p)},
$$
with $c_0$ again fixed by continuity at $r_I$ and $r_O$.

The Hernquist-type model introduces a rational mass function with typical galaxy scale $r_*\gg M_{\rm BH}$,
$$
m(r)=\frac{m_0+m_1r+Mr^2}{(r+r_*)^2}.
$$
Imposing
$$
m(r_I)=M_{\rm BH}, \qquad
r_I^2m'(r_I)+r_I m(r_I)-6m(r_I)^2=0
$$
determines $m_0$ and $m_1$ in terms of $r_I$, $r_*$, and total galaxy mass $M$. Three broad classes arise. Type A is Hernquist-like, with density $\rho\propto r^{-1}$ for $r\ll r_*$ and $\rho\propto r^{-4}$ for $r\gg r_*$. Type $B_0$ has $\rho\propto r^{-2}\to r^{-5}$. The variants $B_+$ and $B_-$ differ according to whether $\rho$ vanishes at a finite $r_O$ or extends to infinity. In all cases,
$$
p_t(r)=\frac{m(r)\rho(r)}{2[r-2m(r)]},
$$
and $\Phi(r)$ can be obtained in closed form in terms of logarithms and arctangents once the three roots of the cubic denominator are known.

## 5. ISCO condition and strong-field behavior

In pure Schwarzschild vacuum of mass $M_{\rm BH}$, timelike circular geodesics exist for $r>3M_{\rm BH}$, and the ISCO is located at
$$
r_{\rm ISCO}^{\rm vac}=6M_{\rm BH}.
$$
With a self-gravitating Einstein cluster present, the ISCO of the combined system is determined by
$$
r^2m'(r)+r\,m(r)-6m(r)^2=0,
$$
subject to
$$
r>3m(r).
$$
The cluster therefore changes the strong-field stability criterion through its mass profile rather than through an externally imposed perturbation [2410.04175].

For Model I, substituting
$$
m(r)=M_{\rm BH}+\alpha(r-r_I)
$$
into the ISCO condition gives
$$
\alpha r_I^2+M_{\rm BH}r_I-6M_{\rm BH}^2=0,
$$
so that
$$
r_I=\frac{-M_{\rm BH}+\sqrt{M_{\rm BH}^2+24\alpha M_{\rm BH}^2}}{2\alpha},
$$
with $0<\alpha<1/3$. As $\alpha\to 0$, one recovers $r_I\to 6M_{\rm BH}$. As $\alpha\to 1/3$, one finds $r_I\to 3M_{\rm BH}$, the photon radius.

Across the model families discussed, the resulting ISCO satisfies
$$
3M_{\rm BH}<r_{\rm ISCO}<6M_{\rm BH}.
$$
The formalism therefore shows that a spiky, self-gravitating halo of particles on circular orbits can move the ISCO inward from the vacuum Schwarzschild value toward the photon sphere. A common misconception is that the Schwarzschild ISCO at $6M_{\rm BH}$ is unchanged so long as spherical symmetry is maintained. The Einstein-cluster construction demonstrates that spherical symmetry alone does not preserve the vacuum ISCO when self-gravitating matter occupies the strong-field region.

## 6. Interpretation, scope, and modeling significance

The formalism treats dark-matter or stellar halos as collisionless particles on circular geodesics. Operationally, one trades the unknown distribution function $n(r)$, or equivalently the angular-momentum profile $L(r)$, for the mass profile $m(r)$; from $m(r)$ one obtains $\rho(r)$ and $p_t(r)$ algebraically, and then $\Phi(r)$ by a single quadrature [2410.04175]. This structure is what makes the toy models analytically tractable and allows the more realistic Hernquist-type construction to remain explicit.

The emphasis on circular orbits is essential. Because $P_r=0$ is built into the Einstein-cluster formalism, the resulting matter model is adapted specifically to orbit-supported configurations. It is therefore not a generic prescription for arbitrary collisionless halos, nor a general anisotropic-fluid ansatz with independently specifiable radial and tangential pressures. A plausible implication is that the formalism is best viewed as a controlled relativistic idealization of centrally concentrated, orbit-supported matter distributions.

Within that idealization, the framework isolates a sharp physical effect: self-gravity of a central spiky distribution modifies geodesic stability in the strong-field region. The analytic solvability of Models I–III and the Hernquist-type example shows that this effect can be traced directly to the mass function $m(r)$ and its derivative, rather than to additional phenomenological assumptions. In that sense, the Einstein cluster formalism provides a compact relativistic mechanism for studying how a galactic dark-matter spike or analogous collisionless structure alters the near-hole orbital structure of a supermassive black hole environment.

Source: https://www.emergentmind.com/topics/einstein-cluster-formalism